An efficient and accurate method for the conservative Swift-Hohenberg equation and its numerical implementation

Abstract

The conservative Swift-Hohenberg equation was introduced to reformulate the phase-field crystal model. A challenge in solving the conservative Swift-Hohenberg equation numerically is how to treat the nonlinear term to preserve mass conservation without compromising efficiency and accuracy. To resolve this problem, we present a linear, high-order, and mass conservative method by placing the linear and nonlinear terms in the implicit and explicit parts, respectively, and employing the implicit-explicit Runge-Kutta method. We show analytically that the method inherits the mass conservation. Numerical experiments are presented demonstrating the efficiency and accuracy of the proposed method. In particular, long time simulation for pattern formation in 2D is carried out, where the phase diagram can be observed clearly. The MATLAB code for numerical implementation of the proposed method is provided in Appendix .