Stability analysis for the Implicit-Explicit discretization of the Cahn-Hilliard equation

Implicit-Explicit methods have been widely used for the efficient numerical simulation of phase field problems such as the Cahn-Hilliard equation or thin film type equations. Due to the lack of maximum principle and stiffness caused by the effect of small dissipation coefficient, most existing theoretical analysis relies on adding additional stabilization terms, mollifying the nonlinearity or introducing auxiliary variables which implicitly either changes the structure of the problem or trades accuracy for stability in a subtle way. In this work we introduce a robust theoretical framework to analyze directly the stability of the standard implicit-explicit approach without stabilization or any other modification. We take the Cahn-Hilliard equation as a model case and prove energy stability under natural time step constraints which are optimal with respect to energy scaling. These settle several questions which have been open since the work of Chen and Shen \cite{CS98}.

[1]  R. Nicolaides,et al.  Numerical analysis of a continuum model of phase transition , 1991 .

[2]  Bo Li,et al.  Center for Scientific Computation And Mathematical Modeling , 2003 .

[3]  Charles M. Elliott,et al.  The global dynamics of discrete semilinear parabolic equations , 1993 .

[4]  Jie Shen,et al.  Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[5]  Fei Liu,et al.  Stabilized semi‐implicit spectral deferred correction methods for Allen–Cahn and Cahn–Hilliard equations , 2015 .

[6]  Feng Chen,et al.  Efcient Energy Stable Schemes with Spectral Discretization in Space for Anisotropic Cahn-Hilliard Systems , 2020 .

[7]  Zhonghua Qiao,et al.  Characterizing the Stabilization Size for Semi-Implicit Fourier-Spectral Method to Phase Field Equations , 2014, SIAM J. Numer. Anal..

[8]  Andreas Prohl,et al.  Error analysis of a mixed finite element method for the Cahn-Hilliard equation , 2004, Numerische Mathematik.

[9]  Tao Tang,et al.  Stability Analysis of Large Time-Stepping Methods for Epitaxial Growth Models , 2006, SIAM J. Numer. Anal..

[10]  Xiaofeng Yang,et al.  Numerical approximations of Allen-Cahn and Cahn-Hilliard equations , 2010 .

[11]  Jie Shen,et al.  A New Class of Efficient and Robust Energy Stable Schemes for Gradient Flows , 2017, SIAM Rev..

[12]  Cheng Wang,et al.  An $H^2$ convergence of a second-order convex-splitting, finite difference scheme for the three-dimensional Cahn–Hilliard equation , 2016 .

[13]  Yinnian He,et al.  On large time-stepping methods for the Cahn--Hilliard equation , 2007 .

[14]  Keith Promislow,et al.  Variational Models of Network Formation and Ion Transport: Applications to Perfluorosulfonate Ionomer Membranes , 2012 .

[15]  D. J. Eyre,et al.  An Unconditionally Stable One-Step Scheme for Gradient Systems , 1997 .

[16]  Thomas J. R. Hughes,et al.  Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models , 2011, J. Comput. Phys..

[17]  Dong Li,et al.  On Second Order Semi-implicit Fourier Spectral Methods for 2D Cahn–Hilliard Equations , 2017, J. Sci. Comput..

[18]  Keith Promislow,et al.  High accuracy solutions to energy gradient flows from material science models , 2014, J. Comput. Phys..

[19]  A. Bertozzi,et al.  Unconditionally stable schemes for higher order inpainting , 2011 .

[20]  Dong Li,et al.  On the stabilization size of semi-implicit fourier-spectral methods for 3D Cahn-Hilliard equations , 2017 .

[21]  Robert L. Pego,et al.  Front migration in the nonlinear Cahn-Hilliard equation , 1989, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[22]  T. Tang.  ON EFFECTIVE NUMERICAL METHODS FOR PHASE-FIELD MODELS , 2019, Proceedings of the International Congress of Mathematicians (ICM 2018).

[23]  Zhonghua Qiao,et al.  Gradient bounds for a thin film epitaxy equation , 2014 .

[24]  Qiang Du,et al.  A Fourier Spectral Moving Mesh Method for the Cahn-Hilliard Equation with Elasticity , 2009 .

[25]  Jie Shen,et al.  Second-order Convex Splitting Schemes for Gradient Flows with Ehrlich-Schwoebel Type Energy: Application to Thin Film Epitaxy , 2012, SIAM J. Numer. Anal..

[26]  Steven M. Wise,et al.  Unconditionally stable schemes for equations of thin film epitaxy , 2010 .

[27]  D. J. Eyre.  Unconditionally Gradient Stable Time Marching the Cahn-Hilliard Equation , 1998 .

[28]  Shuonan Wu,et al.  On the stability and accuracy of partially and fully implicit schemes for phase field modeling , 2016, Computer Methods in Applied Mechanics and Engineering.

[29]  Jie Shen,et al.  Decoupled Energy Stable Schemes for Phase-Field Models of Two-Phase Complex Fluids , 2014, SIAM J. Sci. Comput..

[30]  Xinlong Feng,et al.  Stabilized Crank-Nicolson/Adams-Bashforth Schemes for Phase Field Models , 2013 .

[31]  Cheng Wang,et al.  A convergent convex splitting scheme for the periodic nonlocal Cahn-Hilliard equation , 2014, Numerische Mathematik.

[32]  Tao Tang,et al.  Long Time Numerical Simulations for Phase-Field Problems Using p-Adaptive Spectral Deferred Correction Methods , 2015, SIAM J. Sci. Comput..

[33]  Chert,et al.  Applications of semi-implicit Fourier-spectral method to phase field equations , 2004 .

[34]  Chi-Wang Shu,et al.  Unconditional Energy Stability Analysis of a Second Order Implicit–Explicit Local Discontinuous Galerkin Method for the Cahn–Hilliard Equation , 2017, Journal of Scientific Computing.

[35]  Héctor D. Ceniceros,et al.  Three-dimensional, fully adaptive simulations of phase-field fluid models , 2010, J. Comput. Phys..