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.

Step-by-step solving virtual element schemes based on scalar auxiliary variable with relaxation for Allen–Cahn-type gradient flows

A weighted scalar auxiliary variable method for solving gradient flows: bridging the nonlinear energy-based and Lagrange multiplier approaches

The scalar auxiliary variable (SAV) approach for gradient flows

A New Scalar Auxiliary Variable Approach for Gradient Flows

A New High-Order Compact Scheme Of Unstructured Finite Volume Method

Efficiently high-order time-stepping R-GSAV schemes for the Navier–Stokes-Poisson-Nernst–Planck equations

A novel high-order linearly implicit and energy-stable additive Runge-Kutta methods for gradient flow models

Improved scalar auxiliary variable schemes for original energy stability of gradient flows

A General Class of Linear Unconditionally Energy Stable Schemes for the Gradient Flows

High-efficiency and positivity-preserving stabilized SAV methods for gradient flows

A convergent stochastic scalar auxiliary variable method

Quad-SAV scheme for gradient systems

Second-order linear adaptive time-stepping schemes for the fractional Allen–Cahn equation

A novel energy-optimal scalar auxiliary variable (EOP-SAV) approach for gradient flows

Improving the accuracy and consistency of the energy quadratization method with an energy-optimized technique

Global constraints preserving SAV schemes for gradient flows

Unconditionally energy stable ESAV-VEM schemes with variable time steps for the time fractional Allen-Cahn equation

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

Adaptive linear second-order energy stable schemes for time-fractional Allen-Cahn equation with volume constraint

High order spatial discretization for variational time implicit schemes: Wasserstein gradient flows and reaction-diffusion systems

Supplementary Variable Method for Developing Structure-Preserving Numerical Approximations to Thermodynamically Consistent Partial Differential Equations