On the Accuracy of Partially Implicit Schemes for Phase Field Modeling

Jinchao Xu,Yukun Li,Shuonan Wu

2016-01-01

Abstract:Partially implicit (or partially explicit) schemes, such as the convex splitting schemes (CSS in short), are among the most popular numerical schemes used in phase-field modeling. Using both theoretical analysis and numerical experiments, we demonstrate in this paper that all existing partially implicit schemes for phase-filed simulations may lack convergence accuracy. As the main example, we first show that some well-known CSS can be interpreted as some fully implicit schemes (FIS in short) in disguise. For the Allen-Cahn model, we prove that the standard CSS is exactly the same as the standard FIS but with a (much) smaller time step size and as a result, it would provide an approximation to the original solution of the Allen-Cahn model at a delayed time (although the magnitude of the delay is reduced when the time step size is reduced). Such time delay is also observed for other partially implicit schemes when time step size is not sufficient small. For the Cahn-Hilliard model, we prove that the standard CSS is exactly the same as the standard FIS for a different model that is a (nontrivial) perturbation of the original Cahn-Hilliard model. We then show that each numerical scheme on each time step can be interpreted as the discrete energy minimization problem, the convexity of which is obtained by certain restriction pertaining to the time step size. In this sense, we argue that CSS or other partially implicit scheme can be viewed as adding an artificial convexity into the discretized system. Motivated by the equivalence between CSS and FIS, we propose a modification of a typical FIS for the Allen-Cahn model so that the maximum principle will be valid on the discrete level and we further rigorously prove that, the linearization of such a modified FIS can be uniformly preconditioned by a Poisson-like operator.