A weakly nonlinear, energy stable scheme for the strongly anisotropic Cahn-Hilliard equation and its convergence analysis

Abstract In this paper we propose and analyze a weakly nonlinear, energy stable numerical scheme for the strongly anisotropic Cahn-Hilliard model. In particular, a highly nonlinear and singular anisotropic surface energy makes the PDE system very challenging at both the analytical and numerical levels. To overcome this well-known difficulty, we perform a convexity analysis on the anisotropic interfacial energy, and a careful estimate reveals that all its second order functional derivatives stay uniformly bounded by a global constant. This subtle fact enables one to derive an energy stable numerical scheme. Moreover, a linear approximation becomes available for the surface energy part, and a detailed estimate demonstrates the corresponding energy stability. Its combination with an appropriate treatment for the nonlinear double well potential terms leads to a weakly nonlinear, energy stable scheme for the whole system. In particular, such an energy stability is in terms of the interfacial energy with respect to the original phase variable, and no auxiliary variable needs to be introduced. This has important implications, for example, in the case that the method needs to satisfy a maximum principle. More importantly, with a careful application of the global bound for the second order functional derivatives, an optimal rate convergence analysis becomes available for the proposed numerical scheme, which is the first such result in this area. Meanwhile, for a Cahn-Hilliard system with a sufficiently large degree of anisotropy, a Willmore or biharmonic regularization has to be introduced to make the equation well-posed. For such a physical model, all the presented analyses are still available; the unique solvability, energy stability and convergence estimate can be derived in an appropriate manner. In addition, the Fourier pseudo-spectral spatial approximation is applied, and all the theoretical results could be extended for the fully discrete scheme. Finally, a few numerical results are presented, which confirm the robustness and accuracy of the proposed scheme.

[1]  Martin Burger,et al.  A level set approach to anisotropic flows with curvature regularization , 2007, J. Comput. Phys..

[2]  Cheng Wang,et al.  Convergence analysis and error estimates for a second order accurate finite element method for the Cahn–Hilliard–Navier–Stokes system , 2016, Numerische Mathematik.

[3]  D. W. Hoffman,et al.  A Vector Thermodynamics for Anisotropic Surfaces—II. Curved and Faceted Surfaces , 1974 .

[4]  Jun Zhang,et al.  Efficient, non-iterative, and second-order accurate numerical algorithms for the anisotropic Allen-Cahn Equation with precise nonlocal mass conservation , 2020, J. Comput. Appl. Math..

[5]  Cheng Wang,et al.  Error analysis of a mixed finite element method for a Cahn–Hilliard–Hele–Shaw system , 2017, Numerische Mathematik.

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

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

[8]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[9]  Cheng Wang,et al.  Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation , 2012, J. Comput. Phys..

[10]  Cheng Wang,et al.  A Linear Iteration Algorithm for a Second-Order Energy Stable Scheme for a Thin Film Model Without Slope Selection , 2014, J. Sci. Comput..

[11]  Peter W Voorhees,et al.  A phase-field model for highly anisotropic interfacial energy , 2001 .

[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]  Steven M. Wise,et al.  Analysis of a Mixed Finite Element Method for a Cahn-Hilliard-Darcy-Stokes System , 2013, SIAM J. Numer. Anal..

[14]  Jan S. Hesthaven,et al.  Spectral Methods for Time-Dependent Problems: Contents , 2007 .

[15]  John W. Cahn,et al.  Diffuse interfaces with sharp corners and facets: phase field models with strongly anisotropic surfaces , 1998 .

[16]  Wenqiang Feng,et al.  A Fourier pseudospectral method for the “good” Boussinesq equation with second‐order temporal accuracy , 2014, 1401.6327.

[17]  Cheng Wang,et al.  A second order energy stable scheme for the Cahn-Hilliard-Hele-Shaw equations , 2016, Discrete & Continuous Dynamical Systems - B.

[18]  Steven M. Wise,et al.  Analysis of a Darcy-Cahn-Hilliard Diffuse Interface Model for the Hele-Shaw Flow and Its Fully Discrete Finite Element Approximation , 2011, SIAM J. Numer. Anal..

[19]  S. M. Wise,et al.  Unconditionally Stable Finite Difference, Nonlinear Multigrid Simulation of the Cahn-Hilliard-Hele-Shaw System of Equations , 2010, J. Sci. Comput..

[20]  Morton E. Gurtin,et al.  A regularized equation for anisotropic motion-by-curvature , 1992 .

[21]  Robert F. Sekerka,et al.  Analytical criteria for missing orientations on three-dimensional equilibrium shapes , 2005 .

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

[23]  John W. Cahn,et al.  On Spinodal Decomposition , 1961 .

[24]  Wenqiang Feng,et al.  Convergence analysis and numerical implementation of a second order numerical scheme for the three-dimensional phase field crystal equation , 2016, Comput. Math. Appl..

[25]  Jie Shen,et al.  Efficient energy stable schemes with spectral discretization in space for anisotropic , 2013 .

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

[27]  Alain Miranville,et al.  Existence of solutions for anisotropic Cahn-Hilliard and Allen-Cahn systems in higher space dimensions , 2016 .

[28]  Axel Voigt,et al.  Faceting of Equilibrium and Metastable Nanostructures: A Phase-Field Model of Surface Diffusion Tackling Realistic Shapes , 2015 .

[29]  Cheng Wang,et al.  Long Time Stability of High Order MultiStep Numerical Schemes for Two-Dimensional Incompressible Navier-Stokes Equations , 2016, SIAM J. Numer. Anal..

[30]  Cheng Wang,et al.  Stability and Convergence Analysis of Fully Discrete Fourier Collocation Spectral Method for 3-D Viscous Burgers’ Equation , 2012, J. Sci. Comput..

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

[32]  Zhonghua Qiao,et al.  A Third Order Exponential Time Differencing Numerical Scheme for No-Slope-Selection Epitaxial Thin Film Model with Energy Stability , 2019, Journal of Scientific Computing.

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

[34]  Cheng Wang,et al.  An Energy-Stable and Convergent Finite-Difference Scheme for the Phase Field Crystal Equation , 2009, SIAM J. Numer. Anal..

[35]  Cheng Wang,et al.  A Second-Order Energy Stable BDF Numerical Scheme for the Cahn-Hilliard Equation , 2018 .

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

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

[38]  Zhen Guan,et al.  Convergence analysis for second‐order accurate schemes for the periodic nonlocal Allen‐Cahn and Cahn‐Hilliard equations , 2016, 1610.06973.

[39]  Steven M. Wise,et al.  Convergence Analysis of a Second Order Convex Splitting Scheme for the Modified Phase Field Crystal Equation , 2012, SIAM J. Numer. Anal..

[40]  A. A. Wheeler,et al.  Phase-field theory of edges in an anisotropic crystal , 2006, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[41]  R. Kobayashi.  Modeling and numerical simulations of dendritic crystal growth , 1993 .

[42]  Xiaofeng Yang,et al.  Fast, provably unconditionally energy stable, and second-order accurate algorithms for the anisotropic Cahn–Hilliard Model , 2019, Computer Methods in Applied Mechanics and Engineering.

[43]  Wenqiang Feng,et al.  An energy stable fourth order finite difference scheme for the Cahn-Hilliard equation , 2017, J. Comput. Appl. Math..

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

[45]  A. Quarteroni,et al.  Approximation results for orthogonal polynomials in Sobolev spaces , 1982 .

[46]  A. Miranville.  EXISTENCE OF SOLUTIONS FOR A ONE-DIMENSIONAL ALLEN-CAHN EQUATION , 2013 .

[47]  W. Craig Carter,et al.  Orientation-dependent surface tension functions for surface energy minimizing calculations , 2005 .

[48]  Amanda E. Diegel,et al.  Stability and Convergence of a Second Order Mixed Finite Element Method for the Cahn-Hilliard Equation , 2014, 1411.5248.

[49]  Cheng Wang,et al.  Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential , 2017, J. Comput. Phys. X.

[50]  John Lowengrub,et al.  Simulating interfacial anisotropy in thin-film growth using an extended Cahn-Hilliard model. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[51]  Axel Voigt,et al.  Morphological Evolution of Pit-Patterned Si(001) Substrates Driven by Surface-Energy Reduction , 2017, Nanoscale Research Letters.

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

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

[54]  Cheng Wang,et al.  A Linear Energy Stable Scheme for a Thin Film Model Without Slope Selection , 2012, J. Sci. Comput..

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

[56]  Xiaoming Wang,et al.  Long Time Stability of a Classical Efficient Scheme for Two-dimensional Navier-Stokes Equations , 2011, SIAM J. Numer. Anal..