In this paper, we propose high order and unconditionally energy stable methods for a modified phase field crystal equation by applying the strategy of the energy quadratization Runge-Kutta methods. We transform the original model into an equivalent system with auxiliary variables and quadratic free energy. The modified system preserves the laws of mass conservation and energy dissipation with the associated energy functional. We present rigorous proofs of the mass conservation and energy dissipation properties of the proposed numerical methods and present numerical experiments conducted to demonstrate their accuracy and energy stability. Finally, we compare long-term simulations using an indicator function to characterize the pattern formation.
|Journal||Modelling and Simulation in Materials Science and Engineering|
|State||Published - Mar 2022|
- energy quadratization Runge-Kutta method
- mass conservation
- modified phase field crystal equation
- unconditional energy stability