Abstract
The modified phase field crystal (MPFC) model is a sixth-order evolutive nonlinear partial differential equation which can describe many crystal phenomena. In this paper, we propose a highly efficient and accurate numerical method to construct linear and unconditionally energy stable schemes for the MPFC model. In recent years, the scalar auxiliary variable (SAV) and SAV-based methods have attracted much attention in numerical solution for dissipative systems due to their inherent advantage of preserving certain discrete analogues of the energy dissipation law. The considered numerical schemes are based on a new SAV-type approach which named new scalar auxiliary variable (nSAV) approach. We first give a first-order energy stable numerical scheme by introducing a new SAV R(t). Then, the high-order nSAV schemes based on the k-step backward differentiation formula (BDFk) are constructed. The considered nSAV schemes allow us to construct high-order schemes for both the phase variable ϕ and ψ while only using a first-order approximation of the energy balance equation. To our knowledge, there is no careful research to give high-order accurate but energy stable schemes for the MPFC model. Meanwhile, the proposed approach only needs to solve linear equation with constant coefficients in one time step which is easy to use fast Fourier transform (FFT) to save more CPU time in calculation. Some numerical simulations are demonstrated to verify the accuracy and efficiency of our proposed schemes.
Similar content being viewed by others
References
Baskaran, A., Hu, Z., Lowengrub, J.S., Wang, C., Wise, S.M., Zhou, P.: Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation. J. Comput. Phys. 250, 270–292 (2013)
Baskaran, A., Lowengrub, J.S., Wang, C., Wise, S.M.: Convergence analysis of a second order convex splitting scheme for the modified phase field crystal equation. SIAM J. Numer. Anal. 51, 2851–2873 (2013)
Berry, J., Grant, M., Elder, K.R.: Diffusive atomistic dynamics of edge dislocations in two dimensions. Physical Review E Statistical Nonlinear & Soft Matter Physics 73, 031609 (2006)
Cheng, Q., Liu, C., Shen, J.: A new lagrange multiplier approach for gradient flows. Comput. Methods Appl. Mech. Eng. 367, 113070 (2020)
Elder, K.R., Katakowski, M., Haataja, M., Grant, M.: Modeling elasticity in crystal growth. Phys. Rev. Lett. 88, 245701 (2002)
Huang, F., Shen, J., Yang, Z.: A highly efficient and accurate new scalar auxiliary variable approach for gradient flows. SIAM J. Sci. Comput. 42, A2514–A2536 (2020)
Lee, H.G., Shin, J., Lee, J.-Y.: First-and second-order energy stable methods for the modified phase field crystal equation. Comput. Methods Appl. Mech. Eng. 321, 1–17 (2017)
Li, Q., Mei, L., Yang, X., Li, Y.: Efficient numerical schemes with unconditional energy stabilities for the modified phase field crystal equation. Adv. Comput. Math. 45, 1551–1580 (2019)
Li, X., Shen, J.: Efficient linear and unconditionally energy stable schemes for the modified phase field crystal equation. arXiv:2004.04319 (2020)
Li, X., Shen, J., Liu, Z.: New SAV-pressure correction methods for the Navier-Stokes equations: stability and error analysis. arXiv:2002.09090 (2020)
Li, Y., Kim, J.: An efficient and stable compact fourth-order finite difference scheme for the phase field crystal equation. Comput. Methods Appl. Mech. Eng. 319, 194–216 (2017)
Liu, Z., Li, X.: Efficient modified stabilized invariant energy quadratization approaches for phase-field crystal equation. Numerical Algorithms 85, 107–132 (2020)
Liu, Z., Li, X.: The exponential scalar auxiliary variable (e-SAV) approach for phase field models and its explicit computing. SIAM J. Sci. Comput. 42, B630–B655 (2020)
Liu, Z., Li, X.: Two fast and efficient linear semi-implicit approaches with unconditional energy stability for nonlocal phase field crystal equation. Appl. Numer. Math. 150, 491–506 (2020)
Liu, Z., Li, X.: A highly efficient and accurate exponential semi-implicit scalar auxiliary variable (esi-sav) approach for dissipative system. Journal of Computational Physics 447, 110703 (2021)
Liu, Z., Li, X.: Step-by-step solving schemes based on scalar auxiliary variable and invariant energy quadratization approaches for gradient flows. Numerical Algorithms 89, 65–86 (2022)
Shin, J., Lee, H.G., Lee, J.-Y.: First and second order numerical methods based on a new convex splitting for phase-field crystal equation. J. Comput. Phys. 327, 519–542 (2016)
Stefanovic, P., Haataja, M., Provatas, N.: Phase-field crystals with elastic interactions. Phys. Rev. Lett. 96, 225504 (2006)
Stefanovic, P.N.P., Haataja, M.: Phase field crystal study of deformation and plasticity in nanocrystalline materials. Phys. Rev. E 80, 046107 (2009)
Wang, C., Wise, S. M.: An energy stable and convergent finite-difference scheme for the modified phase field crystal equation. SIAM J. Numer. Anal. 49, 945–969 (2011)
Wang, M., Huang, Q., Wang, C.: A second order accurate scalar auxiliary variable (sav) numerical method for the square phase field crystal equation. J. Sci. Comput. 88, 1–36 (2021)
Wu, K.-A., Adland, A., Karma, A.: Phase-field-crystal model for fcc ordering. Physical Review E Statistical Nonlinear & Soft Matter Physics 81, 061601 (2010)
Xia, B., Mei, C., Yu, Q., Li, Y.: A second order unconditionally stable scheme for the modified phase field crystal model with elastic interaction and stochastic noise effect. Comput. Methods Appl. Mech. Eng. 363, 112795 (2020)
Yang, J., Kim, J.: Energy dissipation–preserving time-dependent auxiliary variable method for the phase-field crystal and the swift–hohenberg models. Numerical Algorithms, pp. 1–30 (2021)
Yang, X.: Efficient linear, stabilized, second-order time marching schemes for an anisotropic phase field dendritic crystal growth model. Comput. Methods Appl. Mech. Eng. 347, 316–339 (2019)
Yang, X., Han, D.: Linearly first-and second-order, unconditionally energy stable schemes for the phase field crystal model. J. Comput. Phys. 330, 1116–1134 (2017)
Yang, X., Ju, L.: Efficient linear schemes with unconditional energy stability for the phase field elastic bending energy model. Comput. Methods Appl. Mech. Eng. 315, 691–712 (2017)
Zhang, J., Yang, X.: Efficient second order unconditionally stable time marching numerical scheme for a modified phase-field crystal model with a strong nonlinear vacancy potential. Comput. Phys. Commun. 245, 106860 (2019)
Zhang, J., Yang, X.: Numerical approximations for a new l2-gradient flow based phase field crystal model with precise nonlocal mass conservation. Comput. Phys. Commun. 243, 51–67 (2019)
Acknowledgements
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Funding
This work is supported by the National Natural Science Foundation of China (Grant Nos. 12001336, 11971276), by the Postdoctoral Science Foundation of China under grant number 2020M672111, and by Shandong Province Natural Science Foundation (Grant No. ZR2020QA030).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Liu, Z., Zheng, N. & Zhou, Z. A highly efficient and accurate new SAV approach for the modified phase field crystal model. Numer Algor 93, 543–562 (2023). https://doi.org/10.1007/s11075-022-01426-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-022-01426-4