Analysis of the second-order BDF scheme with variable steps for the molecular beam epitaxial model without slope selection

In this work, we are concerned with the stability and convergence analysis of the second order BDF (BDF2) scheme with variable steps for the molecular beam epitaxial model without slope selection. We first show that the variable-step BDF2 scheme is convex and uniquely solvable under a weak time-step constraint. Then we show that it preserves an energy dissipation law if the adjacent time-step ratios $r_k:=\tau_k/\tau_{k-1}<3.561.$ Moreover, with a novel discrete orthogonal convolution kernels argument and some new discrete convolutional inequalities, the $L^2$ norm stability and rigorous error estimates are established, under the same step-ratios constraint that ensuring the energy stability., i.e., $0<r_k<3.561.$ This is known to be the best result in literature. We finally adopt an adaptive time-stepping strategy to accelerate the computations of the steady state solution and confirm our theoretical findings by numerical examples.

[1]  M. Crouzeix,et al.  The Convergence of Variable-Stepsize, Variable-Formula, Multistep Methods , 1984 .

[2]  Amar.  Effects of crystalline microstructure on epitaxial growth. , 1996, Physical review. B, Condensed matter.

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

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

[5]  Joakim Becker,et al.  A second order backward difference method with variable steps for a parabolic problem , 1998 .

[6]  Jie Shen,et al.  Highly Efficient and Accurate Numerical Schemes for the Epitaxial Thin Film Growth Models by Using the SAV Approach , 2019, J. Sci. Comput..

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

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

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

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

[11]  P. Thiel,et al.  A Little Chemistry Helps the Big Get Bigger , 2010, Science.

[12]  Lu-ming Zhang,et al.  An adaptive BDF2 implicit time-stepping method for the phase field crystal model , 2020, ArXiv.

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

[14]  Qi Wang,et al.  Arbitrarily High-order Unconditionally Energy Stable Schemes for Thermodynamically Consistent Gradient Flow Models , 2019, Comput. Phys. Commun..

[15]  Yuan Ma,et al.  An adaptive time-stepping strategy for solving the phase field crystal model , 2013, J. Comput. Phys..

[16]  Rolf Dieter Grigorieff,et al.  Stability of multistep-methods on variable grids , 1983 .

[17]  Etienne Emmrich,et al.  Stability and error of the variable two-step BDF for semilinear parabolic problems , 2005 .

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

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

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

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

[22]  Joachim Krug,et al.  Coarsening of Surface Structures in Unstable Epitaxial Growth , 1997 .

[23]  Tao Zhou,et al.  On energy stable, maximum-principle preserving, second order BDF scheme with variable steps for the Allen-Cahn equation , 2020, SIAM J. Numer. Anal..

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

[25]  Zhimin Zhang,et al.  Analysis of adaptive BDF2 scheme for diffusion equations , 2019, Math. Comput..

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