L-stable spectral deferred correction methods and applications to phase field models

This paper presents the L-stable spectral deferred correction (SDC) methods with low stages. These schemes are initiated by the Crank-Nicolson method. We adopt the linear stabilization approach for the phase field models to obtain the linear implicit SDC scheme. This is done by adding and subtracting the linear stabilization operators that are provided for the different phase field problems. Without loss of the low-stage property, the extrapolation technique is also used in the prediction step of the semi-implicit SDC method. Numerical experiments are given to validate the high-order accuracy and the energy decay property of the proposed semi-implicit SDC methods for the Allen-Cahn, Cahn-Hilliard, and molecular beam epitaxy equations.