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 high-efficiency second-order numerical scheme for time-fractional phase field models by using extended SAV method

The Exponential SAV Approach for the Time-Fractional Allen–Cahn and Cahn–Hilliard Phase-Field Models

A Space-Time Fractional Phase-Field Model with Tunable Sharpness and Decay Behavior and Its Efficient Numerical Simulation

An efficient second order stabilized scheme for the two dimensional time fractional Allen-Cahn equation

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

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

High-order energy stable variable-step schemes for the time-fractional Cahn-Hilliard model

On Energy Dissipation Theory and Numerical Stability for Time-Fractional Phase-Field Equations

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

A Numerical Algorithm Based on Quadratic Finite Element for Two-Dimensional Nonlinear Time Fractional Thermal Diffusion Model

A Fractional Phase-Field Model for Two-Phase Flows with Tunable Sharpness: Algorithms and Simulations

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

Two Second-Order and Linear Numerical Schemes for the Multi-Dimensional Nonlinear Time-Fractional Schrödinger Equation

Variable-step numerical schemes and energy dissipation laws for time fractional Cahn–Hilliard model

A numerical scheme for time-fractional Allen–Cahn equation with application in phase separation

A Fourier Spectral Method for Fractional-in-space Cahn–Hilliard Equation

The Adaptive Second-Order Stable Schemes for the Fractional Allen–Cahn Equation

A Second-order Energy Stable Fourier Spectral Method for the Fractional Cahn-Hilliard Equation

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

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

A fully-decoupled second-order-in-time and unconditionally energy stable scheme for a phase-field model of two phase flow with variable density