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 Fully Discrete Finite-Element Numerical Scheme with Second-Order Temporal Accuracy for the Phase-Field Crystal Model

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

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

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

A Linearly Second-Order, Unconditionally Energy Stable Scheme and Its Error Estimates 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

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

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

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

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

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

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

A New Space-Fractional Modified Phase Field Crystal Equation and Its Numerical Algorithm

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

A Second Order Accurate Scalar Auxiliary Variable (SAV) Numerical Method for the Square Phase Field Crystal Equation

Semi-implicit, Unconditionally Energy Stable, Stabilized Finite Element Method Based on Multiscale Enrichment for the Cahn-Hilliard-Navier-Stokes Phase-Field Model

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

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

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

A new Allen–Cahn type two-model phase-field crystal model for fcc ordering and its numerical approximation