Why large time-stepping methods for the Cahn-Hilliard equation is stable

Li, Dong

doi:10.1090/mcom/3768

Why large time-stepping methods for the Cahn-Hilliard equation is stable
HTML articles powered by AMS MathViewer

by Dong Li HTML | PDF

Math. Comp. 91 (2022), 2501-2515 Request permission

Abstract:

We consider the Cahn-Hilliard equation with standard double-well potential. We employ a prototypical class of first order in time semi-implicit methods with implicit treatment of the linear dissipation term and explicit extrapolation of the nonlinear term. When the dissipation coefficient is held small, a conventional wisdom is to add a judiciously chosen stabilization term in order to afford relatively large time stepping and speed up the simulation. In practical numerical implementations it has been long observed that the resulting system exhibits remarkable stability properties in the regime where the stabilization parameter is $\mathcal O(1)$, the dissipation coefficient is vanishingly small and the size of the time step is moderately large. In this work we develop a new stability theory to address this perplexing phenomenon.

References

J. W. Cahn and J. E. Hilliard, Free energy of a nonuniform system. I, Interfacial energy free energy, J. Chem. Phys. 28 (1958) 258–267.
Buyang Li, Jiang Yang, and Zhi Zhou, Arbitrarily high-order exponential cut-off methods for preserving maximum principle of parabolic equations, SIAM J. Sci. Comput. 42 (2020), no. 6, A3957–A3978. MR 4186541, DOI 10.1137/20M1333456
Yinnian He, Yunxian Liu, and Tao Tang, On large time-stepping methods for the Cahn-Hilliard equation, Appl. Numer. Math. 57 (2007), no. 5-7, 616–628. MR 2322435, DOI 10.1016/j.apnum.2006.07.026
Dong Li, Effective maximum principles for spectral methods, Ann. Appl. Math. 37 (2021), no. 2, 131–290. MR 4294331, DOI 10.4208/aam.oa-2021-0003
Dong Li, Chaoyu Quan, and Tao Tang, Stability and convergence analysis for the implicit-explicit method to the Cahn-Hilliard equation, Math. Comp. 91 (2022), no. 334, 785–809. MR 4379976, DOI 10.1090/mcom/3704
Dong Li and Tao Tang, Stability of the semi-implicit method for the Cahn-Hilliard equation with logarithmic potentials, Ann. Appl. Math. 37 (2021), no. 1, 31–60. MR 4284064, DOI 10.4208/aam.OA-2020-0003
Huailing Song and Chi-Wang Shu, Unconditional energy stability analysis of a second order implicit-explicit local discontinuous Galerkin method for the Cahn-Hilliard equation, J. Sci. Comput. 73 (2017), no. 2-3, 1178–1203. MR 3719623, DOI 10.1007/s10915-017-0497-5
Dong Li, Xinwei Yu, and Zhichun Zhai, On the Euler-Poincaré equation with non-zero dispersion, Arch. Ration. Mech. Anal. 210 (2013), no. 3, 955–974. MR 3116009, DOI 10.1007/s00205-013-0662-4
Dong Li, On a frequency localized Bernstein inequality and some generalized Poincaré-type inequalities, Math. Res. Lett. 20 (2013), no. 5, 933–945. MR 3207362, DOI 10.4310/MRL.2013.v20.n5.a9
Wenbin Chen, Sidafa Conde, Cheng Wang, Xiaoming Wang, and Steven M. Wise, A linear energy stable scheme for a thin film model without slope selection, J. Sci. Comput. 52 (2012), no. 3, 546–562. MR 2948706, DOI 10.1007/s10915-011-9559-2
David J. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation, Computational and mathematical models of microstructural evolution (San Francisco, CA, 1998) Mater. Res. Soc. Sympos. Proc., vol. 529, MRS, Warrendale, PA, 1998, pp. 39–46. MR 1676409, DOI 10.1557/PROC-529-39
Cheng Wang, Xiaoming Wang, and Steven M. Wise, Unconditionally stable schemes for equations of thin film epitaxy, Discrete Contin. Dyn. Syst. 28 (2010), no. 1, 405–423. MR 2629487, DOI 10.3934/dcds.2010.28.405
Jie Shen and Xiaofeng Yang, Numerical approximations of Allen-Cahn and Cahn-Hilliard equations, Discrete Contin. Dyn. Syst. 28 (2010), no. 4, 1669–1691. MR 2679727, DOI 10.3934/dcds.2010.28.1669
Chuanju Xu and Tao Tang, Stability analysis of large time-stepping methods for epitaxial growth models, SIAM J. Numer. Anal. 44 (2006), no. 4, 1759–1779. MR 2257126, DOI 10.1137/050628143
Jie Shen, Jie Xu, and Jiang Yang, The scalar auxiliary variable (SAV) approach for gradient flows, J. Comput. Phys. 353 (2018), 407–416. MR 3723659, DOI 10.1016/j.jcp.2017.10.021