Analysis of the Energy Stability for Stabilized Semi-implicit Schemes of the Functionalized Cahn-Hilliard Mass-conserving Gradient Flow Equation

A stabilized semi-implicit scheme was designed in [ 3 ] to solve the Functionalized Cahn-Hilliard (CCH) equation, but there is a lack of theoretical analysis of the energy stability. In this paper, we generalize this scheme to solve the general FCH mass-conserving gradient flow (FCH-MCGF) equation and show the theoretical analysis results about the unique solvability and energy stability. We successfully prove that this scheme is uniquely solvable and energy stable in theory by rewriting the double-well potential function to satisfy the Lipschitz-type condition. The range of stabilization parameters is theoretically given as well. In addition, another similar energy stable scheme is proposed, which slightly widens the range of stabilization parameters in theory and has almost the same precision as the previous one. Both the detailed numerical procedure and the selection of stabilization parameters are presented. Finally, several numerical experiments are performed for the FCH-MCGF equation based on these schemes. Specially, the adaptive time step size is considered in the scheme for the simulations of the phase separation in 2D and 3D, since any time step size can be used according to our theoretical results. Numerical results show that these schemes are energy stable and the large time step size indeed can be used in computations. Moreover, by comprehensive comparisons of stability and accuracy among the stabilized semi-implicit scheme, the convex splitting scheme, and the fully implicit scheme, we conclude that the performance of the stabilized semi-implicit scheme is the best, and the convex splitting scheme performs better than the fully implicit scheme.

[1]  Keith Promislow,et al.  Pearling Bifurcations in the Strong Functionalized Cahn-Hilliard Free Energy , 2017, SIAM J. Math. Anal..

[2]  Cheng Wang,et al.  Preconditioned steepest descent methods for some nonlinear elliptic equations involving p-Laplacian terms , 2016, J. Comput. Phys..

[3]  Axel Voigt,et al.  A new phase-field model for strongly anisotropic systems , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[4]  Andrew M. Stuart,et al.  Model Problems in Numerical Stability Theory for Initial Value Problems , 1994, SIAM Rev..

[5]  Jie Shen,et al.  Efficient energy stable schemes for isotropic and strongly anisotropic Cahn-Hilliard systems with the Willmore regularization , 2018, J. Comput. Phys..

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

[7]  Keith Promislow,et al.  Meander and Pearling of Single-Curvature Bilayer Interfaces in the Functionalized Cahn-Hilliard Equation , 2014, SIAM J. Math. Anal..

[8]  Xiaofeng Yang,et al.  Convergence Analysis for the Invariant Energy Quadratization (IEQ) Schemes for Solving the Cahn–Hilliard and Allen–Cahn Equations with General Nonlinear Potential , 2020, J. Sci. Comput..

[9]  Shibin Dai,et al.  Geometric evolution of bilayers under the functionalized Cahn–Hilliard equation , 2013, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[10]  Cheng Wang,et al.  A Second-Order, Weakly Energy-Stable Pseudo-spectral Scheme for the Cahn–Hilliard Equation and Its Solution by the Homogeneous Linear Iteration Method , 2016, J. Sci. Comput..

[11]  Lin Wang,et al.  On Efficient Second Order Stabilized Semi-implicit Schemes for the Cahn–Hilliard Phase-Field Equation , 2017, J. Sci. Comput..

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

[13]  Zhonghua Qiao,et al.  An Adaptive Time-Stepping Strategy for the Cahn-Hilliard Equation , 2012 .

[14]  Xiao Li,et al.  Stabilized linear semi-implicit schemes for the nonlocal Cahn-Hilliard equation , 2018, J. Comput. Phys..

[15]  Frank S Bates,et al.  On the Origins of Morphological Complexity in Block Copolymer Surfactants , 2003, Science.

[16]  Katarzyna P. Adamala,et al.  Photochemically driven redox chemistry induces protocell membrane pearling and division , 2012, Proceedings of the National Academy of Sciences.

[17]  Frank S. Bates,et al.  Consequences of Nonergodicity in Aqueous Binary PEO-PB Micellar Dispersions , 2004 .

[18]  Cheng Wang,et al.  Second order convex splitting schemes for periodic nonlocal Cahn-Hilliard and Allen-Cahn equations , 2014, J. Comput. Phys..

[19]  Keith Promislow,et al.  PEM Fuel Cells: A Mathematical Overview , 2009, SIAM J. Appl. Math..

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

[21]  Ying Chen,et al.  A Uniquely Solvable, Energy Stable Numerical Scheme for the Functionalized Cahn–Hilliard Equation and Its Convergence Analysis , 2018, J. Sci. Comput..

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

[23]  Xiaofeng Yang,et al.  Energy stable schemes for Cahn-Hilliard phase-field model of two-phase incompressible flows , 2010 .

[24]  Steven M. Wise,et al.  Solving the regularized, strongly anisotropic Cahn-Hilliard equation by an adaptive nonlinear multigrid method , 2007, J. Comput. Phys..

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

[26]  Jie Shen,et al.  A Phase-Field Model and Its Numerical Approximation for Two-Phase Incompressible Flows with Different Densities and Viscosities , 2010, SIAM J. Sci. Comput..

[27]  Héctor D. Ceniceros,et al.  Computation of multiphase systems with phase field models , 2002 .

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

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

[30]  Keith Promislow,et al.  Curvature driven flow of bi-layer interfaces , 2011 .

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

[32]  Jie Shen,et al.  Efficient spectral-Galerkin methods for systems of coupled second-order equations and their applications , 2012, J. Comput. Phys..

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

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

[35]  Yan Xu,et al.  Local Discontinuous Galerkin Methods for the Functionalized Cahn–Hilliard Equation , 2015, J. Sci. Comput..

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

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

[38]  Jian Zhang,et al.  Fast Explicit Integration Factor Methods for Semilinear Parabolic Equations , 2015, J. Sci. Comput..

[39]  Cheng Wang,et al.  Global-in-time Gevrey regularity solutions for the functionalized Cahn-Hilliard equation , 2020, Discrete & Continuous Dynamical Systems - S.

[40]  Keith Promislow,et al.  Existence of pearled patterns in the planar functionalized Cahn–Hilliard equation , 2014, 1410.0447.

[41]  Ying Chen,et al.  Efficient, adaptive energy stable schemes for the incompressible Cahn-Hilliard Navier-Stokes phase-field models , 2016, J. Comput. Phys..

[42]  Jie Shen,et al.  Applications of semi-implicit Fourier-spectral method to phase field equations , 1998 .

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

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