A Third Order Exponential Time Differencing Numerical Scheme for No-Slope-Selection Epitaxial Thin Film Model with Energy Stability

In this paper we propose and analyze a (temporally) third order accurate exponential time differencing (ETD) numerical scheme for the no-slope-selection (NSS) equation of the epitaxial thin film growth model, with Fourier pseudo-spectral discretization in space. A linear splitting is applied to the physical model, and an ETD-based multistep approximation is used for time integration of the corresponding equation. In addition, a third order accurate Douglas-Dupont regularization term, in the form of $$-A {\Delta t}^2 \phi _0 (L_N) \Delta _N^2 ( u^{n+1} - u^n)$$ , is added in the numerical scheme. A careful Fourier eigenvalue analysis results in the energy stability in a modified version, and a theoretical justification of the coefficient A becomes available. As a result of this energy stability analysis, a uniform in time bound of the numerical energy is obtained. And also, the optimal rate convergence analysis and error estimate are derived in details, in the $$\ell ^\infty (0,T; H_h^1) \cap \ell ^2 (0,T; H_h^3)$$ norm, with the help of a careful eigenvalue bound estimate, combined with the nonlinear analysis for the NSS model. This convergence estimate is the first such result for a third order accurate scheme for a gradient flow. Some numerical simulation results are presented to demonstrate the efficiency of the numerical scheme and the third order convergence. The long time simulation results for $$\varepsilon =0.02$$ (up to $$T=3 \times 10^5$$ ) have indicated a logarithm law for the energy decay, as well as the power laws for growth of the surface roughness and the mound width. In particular, the power index for the surface roughness and the mound width growth, created by the third order numerical scheme, is more accurate than those produced by certain second order energy stable schemes in the existing literature.

[1]  Richard L. Schwoebel,et al.  Step Motion on Crystal Surfaces. II , 1966 .

[2]  Endre Süli,et al.  An Implicit Midpoint Spectral Approximation of Nonlocal Cahn-Hilliard Equations , 2014, SIAM J. Numer. Anal..

[3]  Cheng Wang,et al.  A second‐order energy stable backward differentiation formula method for the epitaxial thin film equation with slope selection , 2017, 1706.01943.

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

[5]  Weidong Zhao,et al.  Fast High-Order Compact Exponential Time Differencing Runge–Kutta Methods for Second-Order Semilinear Parabolic Equations , 2016, J. Sci. Comput..

[6]  Ole H. Hald,et al.  Convergence of Fourier Methods for Navier-Stokes Equations , 1981 .

[7]  E Weinan,et al.  Convergence of Fourier methods for the Navier-Stokes equations , 1993 .

[8]  M. Hochbruck,et al.  Exponential integrators , 2010, Acta Numerica.

[9]  Zhonghua Qiao,et al.  Error analysis of a finite difference scheme for the epitaxial thin film model with slope selection with an improved convergence constant , 2017 .

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

[11]  Tao Tang,et al.  An Adaptive Time-Stepping Strategy for the Molecular Beam Epitaxy Models , 2011, SIAM J. Sci. Comput..

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

[13]  Qiang Du,et al.  Efficient and stable exponential time differencing Runge-Kutta methods for phase field elastic bending energy models , 2016, J. Comput. Phys..

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

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

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

[17]  Robert V. Kohn,et al.  Energy-Driven Pattern Formation , 2006 .

[18]  Zhengru Zhang,et al.  The stability and convergence of two linearized finite difference schemes for the nonlinear epitaxial growth model , 2012 .

[19]  Leonardo Golubović,et al.  Interfacial Coarsening in Epitaxial Growth Models without Slope Selection , 1997 .

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

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

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

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

[24]  Hui Zhang,et al.  Convergence of a Fast Explicit Operator Splitting Method for the Epitaxial Growth Model with Slope Selection , 2017, SIAM J. Numer. Anal..

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

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

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

[28]  D. Moldovan,et al.  Interfacial coarsening dynamics in epitaxial growth with slope selection , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

[30]  Jian Zhang,et al.  Fast and accurate algorithms for simulating coarsening dynamics of Cahn–Hilliard equations , 2015 .

[31]  P. Jimack,et al.  High Accuracy Benchmark Problems for Allen-Cahn and Cahn-Hilliard Dynamics , 2019, Communications in Computational Physics.

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

[33]  Xiao Li,et al.  Energy stability and error estimates of exponential time differencing schemes for the epitaxial growth model without slope selection , 2018, Math. Comput..

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

[35]  Robert V. Kohn,et al.  Upper bound on the coarsening rate for an epitaxial growth model , 2003 .

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

[37]  Lili Ju,et al.  Compact implicit integration factor methods for a family of semilinear fourth-order parabolic equations , 2014 .

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

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

[40]  E Weinan.  Convergence of spectral methods for Burgers' equation , 1992 .

[41]  Bo Li,et al.  Epitaxial Growth Without Slope Selection: Energetics, Coarsening, and Dynamic Scaling , 2004, J. Nonlinear Sci..

[42]  Zhi-Zhong Sun,et al.  Stability and convergence of second-order schemes for the nonlinear epitaxial growth model without slope selection , 2014, Math. Comput..

[43]  Wenbin Chen,et al.  A Second Order BDF Numerical Scheme with Variable Steps for the Cahn-Hilliard Equation , 2019, SIAM J. Numer. Anal..

[44]  Wenbin Chen,et al.  A mixed finite element method for thin film epitaxy , 2012, Numerische Mathematik.

[45]  F. Hudda,et al.  Atomic View of Surface Self‐Diffusion: Tungsten on Tungsten , 1966 .

[46]  Bo Li,et al.  High-order surface relaxation versus the Ehrlich–Schwoebel effect , 2006 .

[47]  Alexander Ostermann,et al.  Exponential multistep methods of Adams-type , 2011 .

[48]  S. Cox,et al.  Exponential Time Differencing for Stiff Systems , 2002 .

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

[50]  J. M. Keiser,et al.  A New Class of Time Discretization Schemes for the Solution of Nonlinear PDEs , 1998 .

[51]  Wenbin Chen,et al.  A stabilized second order exponential time differencing multistep method for thin film growth model without slope selection , 2019, ESAIM: Mathematical Modelling and Numerical Analysis.

[52]  Junseok Kim,et al.  Effective Time Step Analysis of a Nonlinear Convex Splitting Scheme for the Cahn–Hilliard Equation , 2019, Communications in Computational Physics.