A fully-decoupled arbitrarily high-order time-stepping scheme based on matrix diagonalization for the anisotropic phase-field dendritic crystal growth model
For researchers simulating dendritic crystal growth, this work provides a more efficient high-order scheme that decouples the linear system, enabling parallel computation and faster simulations.
The paper proposes a fully-decoupled, arbitrarily high-order time-stepping scheme for the anisotropic phase-field dendritic crystal growth model, achieving discrete energy dissipation and improved computational efficiency via matrix diagonalization. Numerical experiments verify convergence, stability, and efficiency, with comparisons showing effectiveness over coupled formulations.
We propose a fully-decoupled arbitrarily high-order time-stepping scheme for the anisotropic phase-field dendritic crystal growth model. The scheme combines an auxiliary-variable formulation with algebraically stable Runge-Kutta methods and satisfies a discrete energy dissipation law. To address the computational bottleneck arising from the coupled linear system in existing high-order schemes, a matrix diagonalization technique is introduced to transform the coupled linear elliptic system into a set of independent constant-coefficient elliptic equations. The resulting equations can be solved separately and in parallel, thereby improving computational efficiency. Numerical experiments in both two and three dimensions are presented to verify the convergence, energy stability, and efficiency of the proposed scheme. Comparisons with the original coupled formulation demonstrate the effectiveness of the matrix diagonalization strategy, while additional tests illustrate the advantages of high-order temporal discretizations. Simulations under different anisotropy coefficients, latent heat parameters, rotation angles, and initial nucleus shapes are also presented to investigate their effects on dendritic morphology.