A high-order and energy stable scheme is developed to simulate phase-field models by combining the semi-implicit spectral deferred correction (SDC) method and the energy stable convex splitting technique. The convex splitting scheme we use here is a linear unconditionally stable method but is only of first-order accuracy, so the SDC method can be used to iteratively improve the rate of convergence. However, it is found that the accuracy improvement may affect the overall energy stability which is intrinsic to the phase-field models. To compromise the accuracy and stability, a local p-adaptive strategy is proposed to enhance the accuracy by sacrificing some local energy stability in an acceptable level. The proposed strategy is found very useful for producing accurate numerical solutions at small time (dynamics) as well as long time (steady state) with reasonably large time stepsizes. Numerical experiments are carried out to demonstrate the high effectiveness of the proposed numerical strategy.
In this paper, stabilized Crank-Nicolson/Adams-Bashforth schemes are presented for the Allen-Cahn and Cahn-Hilliard equations. It is shown that the proposed time discretization schemes are either unconditionally energy stable, or conditionally energy stable under some reasonable stability conditions. Optimal error estimates for the semi-discrete schemes and fully-discrete schemes will be derived. Numerical experiments are carried out to demonstrate the theoretical results.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.