Abstract:In this paper we propose and analyze a second order accurate (in time) numerical scheme for the square phase field crystal equation, a gradient flow modeling crystal dynamics at the atomic scale in space but on diffusive scales in time. Its primary difference with the standard phase field crystal model is an introduction of the 4-Laplacian term in the free energy potential, which in turn leads to a much higher degree of nonlinearity. To make the numerical scheme linear while preserving the nonlinear energy stability, we make use of the scalar auxiliary variable (SAV) approach, in which a second order backward differentiation formula is applied in the temporal stencil. Meanwhile, a direct application of the SAV method faces certain difficulties, due to the involvement of the 4-Laplacian term, combined with a derivation of the lower bound of the nonlinear energy functional. In the proposed numerical method, an appropriate decomposition for the physical energy functional is formulated, so that the nonlinear energy part has a well-established global lower bound, and the rest terms lead to constant-coefficient diffusion terms with positive eigenvalues. In turn, the numerical scheme could be very efficiently implemented by constant-coefficient Poisson-like type solvers (via FFT), and energy stability is established by introducing an auxiliary variable, and an optimal rate convergence analysis is provided for the proposed SAV method. A few numerical experiments are also presented, which confirm the efficiency and accuracy of the proposed scheme.

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

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

Linear and Unconditionally Energy Stable Schemes for the Modified Phase Field Crystal Equation

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

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

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

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

Efficient Fully Discrete Finite-Element Numerical Scheme with Second-Order Temporal Accuracy for the Phase-Field Crystal Model

A New Lagrange Multiplier Method for the Mass-Conserved Allen–Cahn Type Square Phase-Field Crystal Model

A high-efficiency second-order numerical scheme for time-fractional phase field models by using extended SAV method

A Second-Order, Linear, \(\Boldsymbol{l^\infty}\)-Convergent, and Energy Stable Scheme for the Phase Field Crystal Equation.

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

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

Stability and convergence analysis for a new phase field crystal model with a nonlocal Lagrange multiplier

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

Novel Linear Decoupled and Unconditionally Energy Stable Numerical Methods for the Modified Phase Field Crystal Model

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

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

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

Two Finite Difference Schemes for the Phase Field Crystal Equation

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