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.

Phase-field Modeling and Linearly Energy-Stable Runge–Kutta Algorithm of Colloidal Crystals on Curved Surfaces

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

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

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

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

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

Simple and practical method for the simulations of two-component PFC models for binary colloidal crystals on curved surfaces

A unified framework on the original energy laws of three effective classes of Runge-Kutta methods for phase field crystal type models

Unconditionally stable method and numerical solution of the hyperbolic phase-field crystal equation

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

Phase-field Computations of Anisotropic Ice Crystal Growth on a Spherical Surface

Phase-Field Simulations of Crystal Growth with Adaptive Mesh Refinement

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

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

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

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

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

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

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