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.

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

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

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

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

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

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

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

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

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

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

Efficient Numerical Scheme for the Anisotropic 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

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

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

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

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

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

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

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

A Second-Order Crank-Nicolson Leap-Frog Scheme for the Modified Phase Field Crystal Model with Long-Range Interaction

Efficient Linearly And Unconditionally Energy Stable Schemes For The Phase Field Model Of Solid-State Dewetting Problems