High-order time-accurate, efficient, and structure-preserving numerical methods for the conservative Swift-Hohenberg model

In this study, we develop high-order time-accurate, efficient, and energy stable schemes for solving the conservative Swift-Hohenberg equation that can be used to describe the L-2-gradient flow based phase-field crystal dynamics. By adopting a modified exponential scalar auxiliary variable approach, we first transform the original equations into an expanded system. Based on the expanded system, the first-, second-, and third-order time-accurate schemes are constructed using the backward Euler formula, second-order backward difference formula (BDF2), and third-order backward difference formula (BDF3), respectively. The energy dissipation law can be easily proved with respect to a modified energy. In each time step, the local variable is updated by solving one elliptic type equation and the non-local variables are explicitly computed. The whole algorithm is totally decoupled and easy to implement. Extensive numerical experiments in two- and three-dimensional spaces are performed to show the accuracy, energy stability, and practicability of the proposed schemes.

Automated summary

In this study, high-order time-accurate, efficient, and energy stable schemes have been developed for solving the conservative Swift-Hohenberg equation, which describes phase-field crystal dynamics. The original equations are transformed into an expanded system, from which first-, second-, and third-order time-accurate schemes are constructed. The algorithm is decoupled and easy to implement, showing accuracy, energy stability, and practicability in numerical experiments in two- and three-dimensional spaces.