Abstract
It is well-known that the Allen–Cahn equation not only satisfies the energy dissipation law but also possesses the maximum bound principle (MBP) in the sense that the absolute value of its solution is pointwise bounded for all time by some specific constant under appropriate initial/boundary conditions. In recent years, the scalar auxiliary variable (SAV) method and many of its variants have attracted much attention in numerical solutions for gradient flow problems due to their inherent advantage of preserving certain discrete analogues of the energy dissipation law. However, existing SAV schemes usually fail to preserve the MBP when applied to the Allen–Cahn equation. In this paper, we develop and analyze new first- and second-order stabilized exponential-SAV schemes for a class of Allen–Cahn type equations, which are shown to simultaneously preserve the energy dissipation law and MBP in discrete settings. In addition, optimal error estimates for the numerical solutions are rigorously obtained for both schemes. Extensive numerical tests and comparisons are also conducted to demonstrate the performance of the proposed schemes.
Similar content being viewed by others
Data Availability
Enquiries about data availability should be directed to the authors.
References
Allen, S.M., Cahn, J.W.: A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening. Acta Metall. 27, 1085–1095 (1979)
Akrivis, G., Li, B.Y., Li, D.F.: Energy-decaying extrapolated RK-SAV methods for the Allen-Cahn and Cahn-Hilliard equations. SIAM J. Sci. Comput. 41, A3703–A3727 (2019)
Bates, P.W.: On some nonlocal evolution equations arising in materials science. Fields Inst. Commun. 48, 13–52 (2006)
Chen, C.J., Yang, X.F.: Fast, provably unconditionally energy stable, and second-order accurate algorithms for the anisotropic Cahn-Hilliard model. Comput. Methods Appl. Mech. Engrg. 351, 35–59 (2019)
Chen, W.B., Wang, X.M., Yan, Y., Zhang, Z.Y.: A second order BDF numerical scheme with variable steps for the Cahn-Hilliard equation. SIAM J. Numer. Anal. 57, 495–525 (2019)
Cheng, Q., Liu, C., Shen, J.: A new Lagrange multiplier approach for gradient flows. Comput. Methods Appl. Mech. Engrg. 367, 113070 (2020)
Cheng, Q., Liu, C., Shen, J.: Generalized SAV approaches for gradient systems. J. Comput. Appl. Math. 394, 113532 (2021)
Du, Q., Gunzburger, M., Lehoucq, R.B., Zhou, K.: Analysis and approximation of nonlocal diffusion problems with volume constraints. SIAM Rev. 54, 667–696 (2012)
Du, Q., Ju, L., Li, X., Qiao, Z.H.: Maximum principle preserving exponential time differencing schemes for the nonlocal Allen-Cahn equation. SIAM J. Numer. Anal. 57, 875–898 (2019)
Du, Q., Ju, L., Li, X., Qiao, Z.H.: Maximum bound principles for a class of semilinear parabolic equations and exponential time-differencing schemes. SIAM Rev. 63, 317–359 (2021)
Du, Q., Yang, J., Zhou, Z.: Time-fractional Allen-Cahn equations: analysis and numerical methods. J. Sci. Comput. 85, 42 (2020)
Feng, X.L., Tang, T., Yang, J.: Stabilized Crank-Nicolson/Adams-Bashforth schemes for phase field models. East Asian J. Appl. Math. 3, 59–80 (2013)
Guan, Z., Wang, C., Wise, S.M.: A convergent convex splitting scheme for the periodic nonlocal Cahn-Hilliard equation. Numer. Math. 128, 377–406 (2014)
Gui, C.F., Zhao, M.F.: Traveling wave solutions of Allen-Cahn equation with a fractional Laplacian. Ann. Inst. H. Poincaré-An. 32, 785–812 (2015)
Hou, D.M., Azaiez, M., Xu, C.J.: A variant of scalar auxiliary variable approaches for gradient flows. J. Comput. Phys. 395, 307–332 (2019)
Hou, T.L., Leng, H.T.: Numerical analysis of a stabilized Crank-Nicolson/Adams-Bashforth finite difference scheme for Allen-Cahn equations. Appl. Math. Lett. 102, 106150 (2020)
Hou, T.L., Tang, T., Yang, J.: Numerical analysis of fully discretized Crank-Nicolson scheme for fractional-in-space Allen-Cahn equations. J. Sci. Comput. 72, 1214–1231 (2017)
Huang, F.K., Shen, J., Yang, Z.G.: A highly efficient and accurate new scalar auxiliary variable approach for gradient flows. SIAM J. Sci. Comput. 42, A2514–A2536 (2020)
Ju, L., Li, X., Qiao, Z.H., Yang, J.: Maximum bound principle preserving integrating factor Runge-Kutta methods for semilinear parabolic equations. J. Comput. Phys. 439, 110405 (2021)
Ju, L., Li, X., Qiao, Z.H., Zhang, H.: Energy stability and error estimates of exponential time differencing schemes for the epitaxial growth model without slope selection. Math. Comp. 87, 1859–1885 (2018)
Ju, L., Zhang, J., Du, Q.: Fast and accurate algorithms for simulating coarsening dynamics of Cahn-Hilliard equations. Comput. Mater. Sci. 108, 272–282 (2015)
Li, J.W., Ju, L., Cai, Y.Y., Feng, X.L.: Unconditionally maximum bound principle preserving linear schemes for the conservative Allen-Cahn equation with nonlocal constraint. J. Sci. Comput. 87, 98 (2021)
Li, J.W., Li, X., Ju, L., Feng, X.L.: Stabilized integrating factor Runge-Kutta method and unconditional preservation of maximum bound principle. SIAM J. Sci. Comput. 43, A1780–A1802 (2021)
Li, X., Qiao, Z.H., Wang, C.: Convergence analysis for a stabilized linear semi-implicit numerical scheme for the nonlocal Cahn-Hilliard equation. Math. Comp. 90, 171–188 (2021)
Liao, H.L., Tang, T., Zhou, T.: On energy stable, maximum-principle preserving, second-order BDF scheme with variable steps for the Allen-Cahn equation. SIAM J. Numer. Anal. 58, 2294–2314 (2020)
Liu, Z.G., Li, X.L.: The exponential scalar auxiliary variable (E-SAV) approach for phase field models and its explicit computing. SIAM J. Sci. Comput. 42, B630–B655 (2020)
Qiao, Z.H., Zhang, Z.R., Tang, T.: An adaptive time-stepping strategy for the molecular beam epitaxy models. SIAM J. Sci. Comput. 33, 1395–1414 (2011)
Samko, S.G., Kilbas, A.A., Marichev, O.I.: Fractional Integrals and Derivatives. Gordon and Breach, Yverdon (1993)
Shen, J., Tang, T., Yang, J.: On the maximum principle preserving schemes for the generalized Allen-Cahn equation. Commun. Math. Sci. 14, 1517–1534 (2016)
Shen, J., Wang, C., Wang, X.M., Wise, S.M.: Second-order convex splitting schemes for gradient flows with Ehrlich-Schwoebel type energy: application to thin film epitaxy. SIAM J. Numer. Anal. 50, 105–125 (2012)
Shen, J., Xu, J.: Convergence and error analysis for the scalar auxiliary variable (SAV) schemes to gradient flows. SIAM J. Numer. Anal. 56, 2895–2912 (2018)
Shen, J., Xu, J., Yang, J.: The scalar auxiliary variable (SAV) approach for gradient flows. J. Comput. Phys. 353, 407–416 (2018)
Shen, J., Xu, J., Yang, J.: A new class of efficient and robust energy stable schemes for gradient flows. SIAM Rev. 61, 474–506 (2019)
Shen, J., Yang, X.F.: Numerical approximations of Allen-Cahn and Cahn-Hilliard equations. Discrete Contin. Dyn. Syst. 28, 1669–1691 (2010)
Tang, T., Yang, J.: Implicit-explicit scheme for the Allen-Cahn equation preserves the maximum principle. J. Comput. Math. 34, 471–481 (2016)
Wise, S.M., Wang, C., Lowengrub, J.S.: An energy stable and convergent finite difference scheme for the phase field crystal equation. SIAM J. Numer. Anal. 47, 2269–2288 (2009)
Xiao, X.F., Feng, X.L., Yuan, J.Y.: The stabilized semi-implicit finite element method for the surface Allen-Cahn equation. Discrete Contin. Dyn. Syst. Ser. B 22, 2857–2877 (2017)
Xu, C.J., Tang, T.: Stability analysis of large time-stepping methods for epitaxial growth models. SIAM J. Numer. Anal. 44, 1759–1779 (2006)
Xu, Z., Yang, X.F., Zhang, H., Xie, Z.Q.: Efficient and linear schemes for anisotropic Cahn-Hilliard model using the stabilized-invariant energy quadratization (S-IEQ) approach. Comput. Phys. Commun. 238, 36–49 (2019)
Yan, Y., Chen, W.B., Wang, C., Wise, S.M.: A second-order energy stable BDF numerical scheme for the Cahn-Hilliard equation. Commun. Comput. Phys. 23, 572–602 (2018)
Yang, X.F.: Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends. J. Comput. Phys. 327, 294–316 (2016)
Yang, X.F., Zhang, G.D.: Convergence analysis for the invariant energy quadratization (IEQ) schemes for solving the Cahn-Hilliard and Allen-Cahn equations with general nonlinear potential. J. Sci. Comput. 82, 55 (2020)
Funding
The authors have not disclosed any funding.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Competing Interests
The authors have not disclosed any competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This work is supported by the CAS AMSS-PolyU Joint Laboratory of Applied Mathematics. L. Ju’s work is partially supported by US National Science Foundation grant DMS-2109633 and US Department of Energy grant DE-SC0020270. X. Li’s work is partially supported by the Hong Kong Research Council GRF grant 15300821 and the Hong Kong Polytechnic University internal grants 4-ZZMK and 1-BD8N. Z. Qiao’s work is partially supported by the Hong Kong Research Council RFS grant RFS2021-5S03 and GRF grants 15300417 and 15302919 and the Hong Kong Polytechnic University internal grant 4-ZZKK.
Rights and permissions
About this article
Cite this article
Ju, L., Li, X. & Qiao, Z. Stabilized Exponential-SAV Schemes Preserving Energy Dissipation Law and Maximum Bound Principle for The Allen–Cahn Type Equations. J Sci Comput 92, 66 (2022). https://doi.org/10.1007/s10915-022-01921-9
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-022-01921-9