An adaptive BDF2 implicit time-stepping method for the phase field crystal model

An adaptive BDF2 implicit time-stepping method is analyzed for the phase field crystal model. The suggested method is proved to preserve a modified energy dissipation law at the discrete levels if the time-step ratios $r_k:=\tau_k/\tau_{k-1}<3.561$, a recent zero-stability restriction of variable-step BDF2 scheme for ordinary differential problems. By using the discrete orthogonal convolution kernels and the corresponding convolution inequalities, an optimal $L^2$ norm error estimate is established under the weak step-ratio restriction $0<r_k<3.561$ ensuring the energy stability. This is the first time such error estimate is theoretically proved for a nonlinear parabolic equation. On the basis of ample tests on random time meshes, a useful adaptive time-stepping strategy is suggested to efficiently capture the multi-scale behaviors and to accelerate the numerical simulations.

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

[2]  Yibao Li,et al.  An efficient and stable compact fourth-order finite difference scheme for the phase field crystal equation , 2017 .

[3]  Steven M. Wise,et al.  An Energy Stable BDF2 Fourier Pseudo-Spectral Numerical Scheme for the Square Phase Field Crystal Equation , 2019, Communications in Computational Physics.

[4]  Jie Shen,et al.  Spectral Methods: Algorithms, Analysis and Applications , 2011 .

[5]  Jie Shen,et al.  On the maximum principle preserving schemes for the generalized Allen–Cahn equation , 2016 .

[6]  Tao Zhou,et al.  Analysis of the second-order BDF scheme with variable steps for the molecular beam epitaxial model without slope selection , 2020, Science China Mathematics.

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

[8]  Liquan Mei,et al.  Efficient numerical schemes with unconditional energy stabilities for the modified phase field crystal equation , 2019, Adv. Comput. Math..

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

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

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

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

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

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

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

[16]  Badrinarayan P. Athreya,et al.  Using the phase-field crystal method in the multi-scale modeling of microstructure evolution , 2007 .

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

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

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

[20]  Martin Grant,et al.  Modeling elasticity in crystal growth. , 2001, Physical review letters.

[21]  M. Grant,et al.  Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[23]  E. Asadi,et al.  A Review of Quantitative Phase-Field Crystal Modeling of Solid–Liquid Structures , 2015 .

[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]  Cheng Wang,et al.  Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation , 2009, J. Comput. Phys..

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

[27]  Jia Zhao,et al.  Energy-stable Runge-Kutta schemes for gradient flow models using the energy quadratization approach , 2019, Appl. Math. Lett..

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

[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]  Qi Wang,et al.  Linear second order energy stable schemes for phase field crystal growth models with nonlocal constraints , 2020, Comput. Math. Appl..

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

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

[33]  Xiaoli Li,et al.  Efficient modified stabilized invariant energy quadratization approaches for phase-field crystal equation , 2019, Numerical Algorithms.