ENERGY-DECAYING EXTRAPOLATED RK-SAV METHODS FOR THE ALLEN-CAHN AND CAHN-HILLIARD EQUATIONS

We construct and analyze a class of extrapolated and linearized Runge-Kutta (RK) methods, which can be of arbitrarily high order, for the time discretization of the Allen-Cahn and Cahn-Hilliard phase field equations, based on the scalar auxiliary variable (SAV) formulation. We prove that the proposed q-stage RK-SAV methods have qth-order convergence in time and satisfy a discrete version of the energy decay property. Numerical examples are provided to illustrate the discrete energy decay property and accuracy of the proposed methods.