Abstract:An adaptive BDF2 implicit time-stepping method is analyzed for the phase field crystal model. The suggested method is proved to preserve a modified energy dissipation law at the discrete levels if the time-step ratios $r_k:=\tau_k/\tau_{k-1}<3.561$, a recent zero-stability restriction of variable-step BDF2 scheme for ordinary differential problems. By using the discrete orthogonal convolution kernels and the corresponding convolution inequalities, an optimal $L^2$ norm error estimate is established under the weak step-ratio restriction $0<r_k<3.561$ ensuring the energy stability. This is the first time such error estimate is theoretically proved for a nonlinear parabolic equation. On the basis of ample tests on random time meshes, a useful adaptive time-stepping strategy is suggested to efficiently capture the multi-scale behaviors and to accelerate the numerical simulations.

What problem does this paper attempt to address?

This paper attempts to solve several key problems in the numerical simulation of the Phase Field Crystal (PFC) model: 1. **Discrete Preservation of the Energy Dissipation Law**: The paper proposes an adaptive BDF2 implicit time - stepping method. This method can preserve a modified energy dissipation law at the discrete level, that is, under the condition that the time - step ratio $r_k := \frac{\tau_k}{\tau_{k - 1}}<3.561$, the numerical method can preserve the energy dissipation characteristics. This is achieved by using discrete orthogonal convolution kernels and the corresponding convolution inequalities. 2. **Optimal $L^2$ - Norm Error Estimation**: For the first time, the paper theoretically proves that for nonlinear parabolic equations, the adaptive BDF2 method has an optimal $L^2$ - norm error estimation under the weak time - step ratio limit $0 < r_k < 3.561$. This provides theoretical support for the numerical stability of nonlinear problems. 3. **Effective Capture of Multiscale Behaviors**: The paper proposes an effective adaptive time - stepping strategy, which can efficiently capture multiscale behaviors on a random time grid and accelerate numerical simulations. This method is particularly useful when dealing with problems at different time scales and can significantly reduce the computational cost. ### Main Contributions: 1. **Establishment of the Discrete Energy Dissipation Law**: Under the condition that the ratio of adjacent time steps satisfies $0 < r_k < 3.561$, the BDF2 implicit method preserves a modified energy dissipation law at the discrete time level, which ensures the maximum - norm stability of the numerical solution. 2. **Proof of $L^2$ - Norm Convergence**: Under the condition $0 < r_k < 3.561$, the convergence of the BDF2 implicit method in the $L^2$ - norm is proved, and if the condition $0 < r_k < 1+\sqrt{2}$ is further satisfied, second - order accuracy can be achieved. 3. **Extensive Numerical Experiments**: Through a large number of numerical experiments and comparison with the Crank - Nicolson method, the effectiveness of the BDF2 time - stepping method is verified, especially when combined with the adaptive time - stepping strategy. ### Mathematical Formulas: - **Free - Energy Functional**: \[ E[\Phi]=\int_{\Omega}\left(\frac{1}{4}\Phi^{4}+\frac{1}{2}\Phi\left(-\epsilon+(1 + \Delta)^{2}\right)\Phi\right)dx \] - **Phase - Field Crystal Equation**: \[ \partial_{t}\Phi=\Delta\mu\quad\text{with}\quad\mu = \frac{\delta E}{\delta\Phi}=\Phi^{3}-\epsilon\Phi+(1 + \Delta)^{2}\Phi \] - **Discrete Energy Dissipation Law**: \[ \frac{dE}{dt}=\left(\frac{\delta E}{\delta\Phi},\partial_{t}\Phi\right)=(\mu,\Delta\mu)=-\|\nabla\mu\|^{2}\leq0 \] - **BDF2 Formula**: \[ D_{2}v_{n}:=\sum_{k = 1}^{n}b(n)_{n - k}\nabla_{\tau}v_{k} \] where the discrete convolution kernel $b(n)_{n - k}$ is defined as: \[ b(n)_0:=\frac{1 + 2r_n}{\tau_n(1 + r_n)},\quad b(n)_1:=-\frac{r_n^2}{\tau_n(1 + r_n)},\quad b(n)_j := 0

An adaptive BDF2 implicit time-stepping method for the phase field crystal model

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

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

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

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

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

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

A Scalable Implicit Solver for Phase Field Crystal Simulations

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

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

A Second-Order, Global-in-Time Energy Stable Implicit-Explicit Runge–Kutta Scheme for the Phase Field Crystal Equation

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

An efficient and fast adaptive numerical method for a novel phase-field model of crystal growth

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

Phase-Field Simulations of Crystal Growth with Adaptive Mesh Refinement

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

An Adaptive Block Bregman Proximal Gradient Method for Computing Stationary States of Multicomponent Phase-Field Crystal Model

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

Asymptotic Behaviour of Time Stepping Methods for Phase Field Models

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