Asymptotic stability of many numerical schemes for phase-field modeling

Abstract

In the recent breakthrough work 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 u0 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 tn∞. 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*(u0,ε) that ensures the monotonicity and unique solvability of the numerical solutions, where the scaling parameter ε ∈(0,1). We prove that the implicit Euler scheme h*=h*(ε), which is independent of u0 and only depends on ε. 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 ∈fu0∈ Rh*(u0,ε)=0. Various numerical experiments are used to confirm the theoretical analysis.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…