Abstract:In the recent breakthrough work \cite{xu2023lack}, a rigorous numerical analysis was conducted on the numerical solution of a scalar ODE containing a cubic polynomial derived from the Allen-Cahn equation. It was found that only the implicit Euler method converge to the correct steady state for any given initial value $u_0$ under the unique solvability and energy stability. But all the other commonly used second-order numerical schemes exhibit sensitivity to initial conditions and may converge to an incorrect equilibrium state as $t_n\to\infty$. This indicates that energy stability may not be decisive for the long-term qualitative correctness of numerical solutions. We found that using another fundamental property of the solution, namely monotonicity instead of energy stability, is sufficient to ensure that many common numerical schemes converge to the correct equilibrium state. This leads us to introduce the critical step size constant $h^*=h^*(u_0,\epsilon)$ that ensures the monotonicity and unique solvability of the numerical solutions, where the scaling parameter $\epsilon \in(0,1)$. We prove that the implicit Euler scheme $h^*=h^*(\epsilon)$, which is independent of $u_0$ and only depends on $\epsilon$. Hence regardless of the initial value taken, the simulation can be guaranteed to be correct when $h<h^*$. But for various other numerical methods, no mater how small the step size $h$ is in advance, there will always be initial values that cause simulation errors. In fact, for these numerical methods, we prove that $\inf_{u_0\in \mathbb{R}}h^*(u_0,\epsilon)=0$. Various numerical experiments are used to confirm the theoretical analysis.

Asymptotic Behaviour of Time Stepping Methods for Phase Field Models

Asymptotic stability of many numerical schemes for phase-field modeling

On the Stability and Accuracy of Partially and Fully Implicit Schemes for Phase Field Modeling

An Adaptive Time-Stepping Procedure Based On The Scaled Boundary Finite Element Method For Elastodynamics

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

On Efficient Second Order Stabilized Semi-implicit Schemes for the Cahn–Hilliard Phase-Field Equation

Stability and Robustness of Time-discretization Schemes for the Allen-Cahn Equation via Bifurcation and Perturbation Analysis

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

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

Simple maximum-principle preserving time-stepping methods for time-fractional Allen-Cahn equation

Explicit Local Time-Stepping for the Inhomogeneous Wave Equation with Optimal Convergence

Energy stability and convergence of variable-step L1 scheme for the time fractional Swift-Hohenberg model

Highly efficient and energy dissipative schemes for the time fractional Allen-Cahn equation

Energy stability of variable-step L1-type schemes for time-fractional Cahn-Hilliard model

Finite element solution of nonlocal Cahn-Hilliard equations with feedback control time step size adaptivity

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

Stability Analysis and Error Estimate of the Explicit Single-Step Time-Marching Discontinuous Galerkin Methods with Stage-Dependent Numerical Flux Parameters for a Linear Hyperbolic Equation in One Dimension

Asymptotically compatible energy and dissipation law of the nonuniform L2-$1_σ$ scheme for time fractional Allen-Cahn model

An efficient and robust staggered scheme based on adaptive time field for phase field fracture model

Energetic Spectral-Element Time Marching Methods for Phase-Field Nonlinear Gradient Systems