An adaptive time-stepping strategy for solving the phase field crystal model

In this work, we will propose an adaptive time step method for simulating the dynamics of the phase field crystal (PFC) model. The numerical simulation of the PFC model needs long time to reach steady state, and then large time-stepping method is necessary. Unconditionally energy stable schemes are used to solve the PFC model. The time steps are adaptively determined based on the time derivative of the corresponding energy. It is found that the use of the proposed time step adaptivity cannot only resolve the steady state solution, but also the dynamical development of the solution efficiently and accurately. The numerical experiments demonstrate that the CPU time is significantly saved for long time simulations.

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

[2]  R. Nicolaides,et al.  Numerical analysis of a continuum model of phase transition , 1991 .

[3]  Yinnian He,et al.  On large time-stepping methods for the Cahn--Hilliard equation , 2007 .

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

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

[6]  J. Swift,et al.  Hydrodynamic fluctuations at the convective instability , 1977 .

[7]  Susan E. Minkoff,et al.  A comparison of adaptive time stepping methods for coupled flow and deformation modeling , 2006 .

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

[9]  B. Vollmayr-Lee,et al.  Fast and accurate coarsening simulation with an unconditionally stable time step. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

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

[12]  Axel Voigt,et al.  Nucleation and growth by a phase field crystal (PFC) model , 2007 .

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

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

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

[16]  James A. Warren,et al.  An efficient algorithm for solving the phase field crystal model , 2008, J. Comput. Phys..

[17]  Yunqing Huang,et al.  Moving mesh methods with locally varying time steps , 2004 .

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

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

[20]  R. Li,et al.  Moving Mesh Method with Error-Estimator-Based Monitor and Its Applications to Static Obstacle Problem , 2004, J. Sci. Comput..

[21]  Gustaf Söderlind,et al.  Automatic Control and Adaptive Time-Stepping , 2002, Numerical Algorithms.

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

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

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