In this article, we develop a linear, unconditionally energy stable computational scheme for solving the dendritic crystal growth model with the orientational field. We apply the phase field model to describe the evolution of crystal with rotation. The model, which couples the heat equation and anisotropic Allen–Cahn type equation, is a complicated nonlinear system. The time integration is based on the second-order Crank–Nicolson method. The anisotropic coefficient is treated by using the invariant energy quadratization. We mathematically prove that the proposed method is unconditionally energy stable. The second-order spatial and temporal accuracy will be preserved for the numerical approximation. Various computational tests are performed to show the accuracy, stability, and efficiency of the proposed scheme.
Bibliographical noteFunding Information:
Y.B. Li is supported by the Fundamental Research Funds for the Central Universities (No. XTR042019005 ). The corresponding author (J.S. Kim) was supported by Korea University Grant. The authors are grateful to the reviewers whose valuable suggestions and comments significantly improved the quality of this paper.
© 2022 IMACS
- Crystal growth model
- Orientational field model
- Phase field model
- Unconditionally energy-stable
ASJC Scopus subject areas
- Numerical Analysis
- Computational Mathematics
- Applied Mathematics