Abstract:SIAM Journal on Numerical Analysis, Volume 62, Issue 6, Page 2667-2697, December 2024. We develop a two-stage, second-order, global-in-time energy stable implicit-explicit Runge–Kutta (IMEX RK(2, 2)) scheme for the phase field crystal equation with an [math] time step constraint, and without the global Lipschitz assumption. A linear stabilization term is introduced to the system with Fourier pseudo-spectral spatial discretization, and a novel compact reformulation is devised by rewriting the IMEX RK(2, 2) scheme as an approximation to the variation-of-constants formula. Under the assumption that all stage solutions are a priori bounded in the [math] norm, we first demonstrate that the original energy obtained by this second-order scheme is nonincreasing for any time step with a sufficiently large stabilization parameter. To justify the a priori [math] bound assumption, we establish a uniform-in-time [math] estimate for all stage solutions, subject to an [math] time step constraint. This results in a uniform-in-time bound for all stage solutions through discrete Sobolev embedding from [math] to [math]. Consequently, we achieve an [math] stabilization parameter, ensuring global-in-time energy stability. Additionally, we provide an optimal rate convergence analysis and error estimate for the IMEX RK(2, 2) scheme in the [math] norm. The global-in-time energy stability represents a novel achievement for a two-stage, second-order accurate scheme for a gradient flow without the global Lipschitz assumption. Numerical experiments substantiate the second-order accuracy and energy stability of the proposed scheme.

Efficient and energy stable numerical schemes for the two-mode phase field crystal equation

Efficient Numerical Schemes with Unconditional Energy Stabilities for the Modified Phase Field Crystal Equation

Efficient Second-Order Unconditionally Stable Numerical Schemes for the Modified Phase Field Crystal Model with Long-Range Interaction

Efficient Unconditionally Stable Numerical Schemes for a Modified Phase Field Crystal Model with a Strong Nonlinear Vacancy Potential

A Second-Order, Uniquely Solvable, Energy Stable BDF Numerical Scheme for the Phase Field Crystal Model

Parallel Energy-Stable Phase Field Crystal Simulations Based on Domain Decomposition Methods.

An Efficient and Stable Compact Fourth-Order Finite Difference Scheme for the Phase Field Crystal Equation

An Efficient Linear Second Order Unconditionally Stable Direct Discretization Method for the Phase-Field Crystal Equation on Surfaces

A Linearly Second-Order, Unconditionally Energy Stable Scheme and Its Error Estimates for the Modified Phase Field Crystal Equation

High-order energy stable schemes of incommensurate phase-field crystal model

Efficient numerical methods for computing the stationary states of phase field crystal models

Efficient and unconditionally energy stable exponential-SAV schemes for the phase field crystal equation

A Second Order Unconditionally Stable Scheme for the Modified Phase Field Crystal Model with Elastic Interaction and Stochastic Noise Effect

A Second-Order, Global-in-Time Energy Stable Implicit-Explicit Runge–Kutta Scheme for the Phase Field Crystal Equation

New third-order convex splitting methods and analysis for the phase field crystal equation

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

An Unconditionally Energy Stable Method for the Anisotropic Phase-Field Crystal Model in Two Dimension

High-oder Symplectic FDTD Scheme for Solving Time-Dependent Schrdinger Equation

An efficient numerical method for the anisotropic phase field dendritic crystal growth model

Global-in-time energy stability analysis for the exponential time differencing Runge-Kutta scheme for the phase field crystal equation