Linear Second Order Energy Stable Schemes of Phase Field Model with Nonlocal Constraints for Crystal Growth
This work provides efficient numerical methods for a phase field model with nonlocal constraints, relevant to researchers in computational materials science and crystal growth simulation.
The authors develop linear, second-order, unconditionally energy stable numerical schemes for the nonlocal Allen-Cahn model for crystal growth, which conserves mass. The nonlocal model exhibits dynamics slower than the classic Allen-Cahn but faster than Cahn-Hilliard, offering an alternative for crystal growth simulations.
We present a set of linear, second order, unconditionally energy stable schemes for the Allen-Cahn model with a nonlocal constraint for crystal growth that conserves the mass of each phase. Solvability conditions are established for the linear systems resulting from the linear schemes. Convergence rates are verified numerically. Dynamics obtained using the nonlocal Allen-Cahn model are compared with the one obtained using the classic Allen-Cahn model as well as the Cahn-Hilliard model, demonstrating slower dynamics than that of the Allen-Cahn model but faster dynamics than that of the Cahn-Hillard model. Thus, the nonlocal Allen-Cahn model can be an alternative to the Cahn-Hilliard model in simulating crystal growth. Two Benchmark examples are presented to illustrate the prediction made with the nonlocal Allen-Cahn model in comparison to those made with the Allen-Cahn model and the Cahn- Hillard model.