修正相场晶体模型的带松弛技术的ESAV方法
ESAV Method with Relaxation Technique for the Modified Phase-Field Crystal Model
摘要: 本文首先提出了修正相场晶体(MPFC)模型的一个有效的松弛指数标量辅助变量方法(R-ESAV)。该方法通过引入松弛机制,使离散能量更接近原始自由能,同时保持线性化与无条件能量稳定性的优势。本文构造了一阶向后欧拉格式,严格证明了其无条件能量稳定性,并给出了一阶误差估计的理论分析。数值算例验证了方法的准确性与有效性。
Abstract: This paper first proposes an efficient relaxed exponential scalar auxiliary variable (R-ESAV) method for the modified phase field crystal (MPFC) model. By introducing a relaxation mechanism, the proposed method enables the discrete energy to more closely approximate the original free energy, while preserving the advantages of linearization and unconditional energy stability. A first-order backward Euler scheme is constructed, and its unconditional energy stability is rigorously proved. Furthermore, a theoretical analysis of the first-order error estimate is provided. Numerical examples validate the accuracy and effectiveness of the proposed method.
文章引用:刘苗, 贾宏恩. 修正相场晶体模型的带松弛技术的ESAV方法[J]. 应用数学进展, 2026, 15(4): 538-555. https://doi.org/10.12677/aam.2026.154181

参考文献

[1] Elder, K. and Grant, M. (2004) Modeling Elastic and Plastic Deformations in Nonequilibrium Processing Using Phase Field Crystals. Physical Review E, 70, Article 051605.
[2] Elder, K.R., Katakowski, M., Haataja, M. and Grant, M. (2002) Modeling Elasticity in Crystal Growth. Physical Review Letters, 88, Article 245701. [Google Scholar] [CrossRef] [PubMed]
[3] Stefanovic, P., Haataja, M. and Provatas, N. (2006) Phase-field Crystals with Elastic Interactions. Physical Review Letters, 96, Article 225504. [Google Scholar] [CrossRef] [PubMed]
[4] Stefanovic, P., Haataja, M. and Provatas, N. (2009) Phase Field Crystal Study of Deformation and Plasticity in Nanocrystalline Materials. Physical Review E, 80, Article 046107. [Google Scholar] [CrossRef] [PubMed]
[5] Baskaran, A., Hu, Z., Lowengrub, J.S., Wang, C., Wise, S.M. and Zhou, P. (2013) Energy Stable and Efficient Finite-Difference Nonlinear Multigrid Schemes for the Modified Phase Field Crystal Equation. Journal of Computational Physics, 250, 270-292. [Google Scholar] [CrossRef
[6] Wu, K.A., Adland, A. and Karma, A. (2010) Phase-Field-Crystal Model for FCC Ordering. Physical Review E, 81, Article 061601. [Google Scholar] [CrossRef] [PubMed]
[7] Elder, K.R., Provatas, N., Berry, J., Stefanovic, P. and Grant, M. (2007) Phase-Field Crystal Modeling and Classical Density Functional Theory of Freezing. Physical Review B, 75, Article 064107. [Google Scholar] [CrossRef
[8] Lee, H.G., Shin, J. and Lee, J. (2017) First-and Second-Order Energy Stable Methods for the Modified Phase Field Crystal Equation. Computer Methods in Applied Mechanics and Engineering, 321, 1-17. [Google Scholar] [CrossRef
[9] Shin, J., Lee, H.G. and Lee, J. (2016) First and Second Order Numerical Methods Based on a New Convex Splitting for Phase-Field Crystal Equation. Journal of Computational Physics, 327, 519-542. [Google Scholar] [CrossRef
[10] Li, Q., Mei, L., Yang, X., et al. (2019) Efficient Numerical Schemes with Unconditional Energy Stabilities for the Modified Phase Field Crystal Equation. Advances in Computational Mathematics, 45, 1551-1580. [Google Scholar] [CrossRef
[11] Liu, Z. and Li, X. (2020) Efficient Modified Stabilized Invariant Energy Quadratization Approaches for Phase-Field Crystal Equation. Numerical Algorithms, 85, 107-132. [Google Scholar] [CrossRef
[12] Yu, H. and Yang, X. (2017) Numerical Approximations for a Phase-Field Moving Contact Line Model with Variable Densities and Viscosities. Journal of Computational Physics, 334, 665-686. [Google Scholar] [CrossRef
[13] Fu, Z. and Yang, J. (2022) Energy-Decreasing Exponential Time Differencing Runge-Kutta Methods for Phase-Field Models. Journal of Computational Physics, 454, Article 110943. [Google Scholar] [CrossRef
[14] Li, X. and Shen, J. (2021) Efficient Linear and Unconditionally Energy Stable Schemes for the Modified Phase Field Crystal Equation. Science China Mathematics, 64, 1-18.
[15] Shen, J., Xu, J. and Yang, J. (2018) The Scalar Auxiliary Variable (SAV) Approach for Gradient Flows. Journal of Computational Physics, 353, 407-416. [Google Scholar] [CrossRef
[16] Wang, M., Huang, Q. and Wang, C. (2021) A Second Order Accurate Scalar Auxiliary Variable (SAV) Numerical Method for the Square Phase Field Crystal Equation. Journal of Scientific Computing, 88, Article No. 33. [Google Scholar] [CrossRef
[17] Yang, J. and Kim, J. (2022) Energy Dissipation-Preserving Time-Dependent Auxiliary Variable Method for the Phase-Field Crystal and the Swift-Hohenberg Models. Numerical Algorithms, 89, 1389-1418.
[18] Zhang, J. and Yang, X. (2019) Efficient Second Order Unconditionally Stable Time Marching Numerical Scheme for a Modified Phase-Field Crystal Model with a Strong Nonlinear Vacancy Potential. Computer Physics Communications, 245, Article 106860. [Google Scholar] [CrossRef
[19] Liu, Z. and Li, X. (2020) The Exponential Scalar Auxiliary Variable (E-SAV) Approach for Phase Field Models and Its Explicit Computing. SIAM Journal on Scientific Computing, 42, B630-B655. [Google Scholar] [CrossRef
[20] Zhang, Y. and Shen, J. (2022) A Generalized SAV Approach with Relaxation for Dissipative Systems. Journal of Computational Physics, 464, Article 111311. [Google Scholar] [CrossRef
[21] Liu, Z. and Li, X. (2021) A Highly Efficient and Accurate Exponential Semi-Implicit Scalar Auxiliary Variable (ESI-SAV) Approach for Dissipative System. Journal of Computational Physics, 447, Article 110703. [Google Scholar] [CrossRef
[22] Liu, Z., Zheng, N. and Zhou, Z. (2023) A Highly Efficient and Accurate New SAV Approach for the Modified Phase Field Crystal Model. Numerical Algorithms, 93, 543-562. [Google Scholar] [CrossRef
[23] Wang, C. and Wise, S.M. (2010) Global Smooth Solutions of the Three-Dimensional Modified Phase Field Crystal Equation. Methods and Applications of Analysis, 17, 191-212. [Google Scholar] [CrossRef
[24] Yang, X., Zhao, J. and Wang, Q. (2017) Numerical Approximations for the Molecular Beam Epitaxial Growth Model Based on the Invariant Energy Quadratization Method. Journal of Computational Physics, 333, 104-127.