Numerical approximations of the Cahn-Hilliard and Allen-Cahn Equations with general nonlinear potential using the Invariant Energy Quadratization approach

In this paper, we carry out stability and error analyses for two first-order, semi-discrete time stepping schemes, which are based on the newly developed Invariant Energy Quadratization approach, for solving the well-known Cahn-Hilliard and Allen-Cahn equations with general nonlinear bulk potentials. Some reasonable sufficient conditions about boundedness and continuity of the nonlinear functional are given in order to obtain optimal error estimates. These conditions are naturally satisfied by two commonly used nonlinear potentials including the double-well potential and regularized logarithmic Flory-Huggins potential. The well-posedness, unconditional energy stabilities and optimal error estimates of the numerical schemes are proved rigorously.

[1]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[2]  J. Cahn,et al.  A Microscopic Theory for Domain Wall Motion and Its Experimental Verification in Fe‐Al Alloy Domain Growth Kinetics , 1977 .

[3]  Charles M. Elliott,et al.  Numerical analysis of the Cahn-Hilliard equation with a logarithmic free energy , 1992 .

[4]  C. M. Elliott,et al.  On the Cahn-Hilliard equation with degenerate mobility , 1996 .

[5]  Yunzhi Wang,et al.  The continuum field approach to modeling microstructural evolution , 1996 .

[6]  J. Lowengrub,et al.  Quasi–incompressible Cahn–Hilliard fluids and topological transitions , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

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

[8]  Andreas Prohl,et al.  Numerical analysis of the Allen-Cahn equation and approximation for mean curvature flows , 2003, Numerische Mathematik.

[9]  Jie Shen,et al.  A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method , 2003 .

[10]  Q. Du,et al.  A phase field approach in the numerical study of the elastic bending energy for vesicle membranes , 2004 .

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

[12]  Xiaofeng Yang,et al.  Dynamics of Defect Motion in Nematic Liquid Crystal Flow: Modeling and Numerical Simulation , 2007 .

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

[14]  Axel Voigt,et al.  Phase-field modeling of the dynamics of multicomponent vesicles: Spinodal decomposition, coarsening, budding, and fission. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  Cheng Wang,et al.  Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation , 2009, J. Comput. Phys..

[16]  Christian Miehe,et al.  A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits , 2010 .

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

[18]  Xiaofeng Yang,et al.  Journal of Non-newtonian Fluid Mechanics Shear Cell Rupture of Nematic Liquid Crystal Droplets in Viscous Fluids , 2011 .

[19]  Steven M. Wise,et al.  An Energy Stable and Convergent Finite-Difference Scheme for the Modified Phase Field Crystal Equation , 2011, SIAM J. Numer. Anal..

[20]  Franck Boyer,et al.  Numerical schemes for a three component Cahn-Hilliard model , 2011 .

[21]  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..

[22]  Xesús Nogueira,et al.  An unconditionally energy-stable method for the phase field crystal equation , 2012 .

[23]  Junseok Kim.  Phase-Field Models for Multi-Component Fluid Flows , 2012 .

[24]  Sébastian Minjeaud.  An unconditionally stable uncoupled scheme for a triphasic Cahn–Hilliard/Navier–Stokes model , 2013 .

[25]  Francisco Guillén-González,et al.  On linear schemes for a Cahn-Hilliard diffuse interface model , 2013, J. Comput. Phys..

[26]  Francisco Guillén-González,et al.  Second order schemes and time-step adaptivity for Allen-Cahn and Cahn-Hilliard models , 2014, Comput. Math. Appl..

[27]  van der Kg Kristoffer Zee,et al.  Stabilized second‐order convex splitting schemes for Cahn–Hilliard models with application to diffuse‐interface tumor‐growth models , 2014, International journal for numerical methods in biomedical engineering.

[28]  Xiaofeng Yang,et al.  Decoupled energy stable schemes for phase-field vesicle membrane model , 2015, J. Comput. Phys..

[29]  Jia Zhao,et al.  Modeling the Excess Cell Surface Stored in a Complex Morphology of Bleb-Like Protrusions , 2016, PLoS Comput. Biol..

[30]  Xiaofeng Yang,et al.  Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends , 2016, J. Comput. Phys..

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

[32]  Jun Li,et al.  Energy Stable Numerical Schemes for a Hydrodynamic Model of Nematic Liquid Crystals , 2016, SIAM J. Sci. Comput..

[33]  Xiaofeng Yang,et al.  Numerical approximations for a three-component Cahn–Hilliard phase-field model based on the invariant energy quadratization method , 2017 .

[34]  Xiaofeng Yang,et al.  Numerical approximations for a phase field dendritic crystal growth model based on the invariant energy quadratization approach , 2017 .

[35]  Lili Ju,et al.  Efficient linear schemes with unconditional energy stability for the phase field elastic bending energy model , 2017 .

[36]  Daozhi Han,et al.  Numerical Analysis of Second Order, Fully Discrete Energy Stable Schemes for Phase Field Models of Two-Phase Incompressible Flows , 2017, J. Sci. Comput..

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

[38]  Xiaofeng Yang,et al.  Numerical approximations for a phase-field moving contact line model with variable densities and viscosities , 2017, J. Comput. Phys..

[39]  Xiaofeng Yang,et al.  Linear, Second order and Unconditionally Energy Stable schemes for a phase-field moving contact line Model , 2017 .

[40]  Lili Ju,et al.  Linear and unconditionally energy stable schemes for the binary fluid–surfactant phase field model , 2017, 1701.07446.

[41]  Jia Zhao,et al.  Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method , 2017, J. Comput. Phys..

[42]  Xiaofeng Yang,et al.  Second Order, Linear, and Unconditionally Energy Stable Schemes for a Hydrodynamic Model of Smectic-A Liquid Crystals , 2017, SIAM J. Sci. Comput..

[43]  Daozhi Han,et al.  Linearly first- and second-order, unconditionally energy stable schemes for the phase field crystal model , 2017, J. Comput. Phys..

[44]  Jie Shen,et al.  Efficient and accurate numerical schemes for a hydro-dynamically coupled phase field diblock copolymer model , 2017, J. Comput. Phys..

[45]  Xiaofeng Yang,et al.  A novel linear second order unconditionally energy stable scheme for a hydrodynamic Q-tensor model of liquid crystals , 2017 .

[46]  Xiaofeng Yang,et al.  Numerical Approximations for the Cahn–Hilliard Phase Field Model of the Binary Fluid-Surfactant System , 2017, Journal of Scientific Computing.

[47]  Xiaofeng Yang,et al.  Efficient Second Order Unconditionally Stable Schemes for a Phase Field Moving Contact Line Model Using an Invariant Energy Quadratization Approach , 2018, SIAM J. Sci. Comput..

[48]  Jiang Yang,et al.  The scalar auxiliary variable (SAV) approach for gradient flows , 2018, J. Comput. Phys..