A class of efficient high-order time-stepping methods for the anisotropic phase-field dendritic crystal growth model

Abstract

In this paper, we develop and analyze a novel class of arbitrarily high-order and unconditionally energy stable schemes for the anisotropic phase-field dendritic crystal growth model, which is a highly nonlinear system that combines the anisotropic Allen–Cahn equation with the thermal equation. The proposed schemes are based on an extrapolated and linearized Runge–Kutta method for an auxiliary variable reformulation of the crystal growth model. A delicate implementation demonstrates that the proposed method can be realized in a very efficient way, requiring only the solution of a coupled linear elliptic system at each time step. We prove that the constructed schemes satisfy the energy dissipation property and demonstrate the convergence order through a consistency error analysis. To the best of our knowledge, this is the first unconditionally stable scheme of arbitrarily high-order for the anisotropic phase-field dendritic crystal growth model. Numerical experiments for two and three dimensional dendritic crystal growth are simulated to verify our theoretical accuracy as well as the efficiency of the proposed method. Furthermore, to emphasize the superiority of the high-order schemes, we present a numerical comparison that directly contrasts its performance with the lower-order schemes.