Asymptotic behaviour of time stepping methods for phase field models

DOI10.1007/S10915-020-01391-XMaRDI QIDQ2027928zbMATH OpenOpenAlexFDO

Authors Xinyu Cheng, Dong Li, Brian Wetton, Keith Promislow

Publication date 28 May 2021

Published in Journal of Scientific Computing (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/1905.02299

zbMATH Keywords

Allen-Cahn equationCahn-Hilliard equationenergy stabilityphase field modeltime stepping

Mathematics Subject Classification ID

Asymptotic behavior of solutions to PDEs (35B40) KdV equations (Korteweg-de Vries equations) (35Q53) Numerical solutions to abstract evolution equations (65J08) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06)

Abstract: Adaptive time stepping methods for metastable dynamics of the Allen Cahn and Cahn Hilliard equations are investigated in the spatially continuous, semi-discrete setting. We analyse the performance of a number of first and second order methods, formally predicting step sizes required to satisfy specified local truncation error $sigma$ in the limit of small order parameter $epsilonightarrow0$ during meta-stable dynamics. The formal predictions are made under stability assumptions that include the preservation of the asymptotic structure of the diffuse interface, a concept we call profile fidelity. In this setting, definite statements about the relative behaviour of time stepping methods can be made. Some methods, including all so-called energy stable methods but also some fully implicit methods, require asymptotically more time steps than others.The formal analysis is confirmed in computational studies. We observe that some provably energy stable methods popular in the literature perform worse than some more standard schemes. We show further that when Backward Euler is applied to meta-stable Allen Cahn dynamics, the energy decay and profile fidelity properties for these discretizations are preserved for much larger time steps than previous analysis would suggest. The results are established asymptotically for general interfaces, with a rigorous proof for radial interfaces. It is shown analytically and computationally that for most reaction terms, Eyre type time stepping performs asymptotically worse due to loss of profile fidelity.

Recommendations

• Second order schemes and time-step adaptivity for Allen-Cahn and Cahn-Hilliard models
• Parameter-Free Time Adaptivity Based on Energy Evolution for the Cahn-Hilliard Equation
• On large time-stepping methods for the Cahn-Hilliard equation
• Error estimates for time discretizations of Cahn-Hilliard and Allen-Cahn phase-field models for two-phase incompressible flows
• On effective numerical methods for phase-field models

Cites work

• A diffuse-interface method for simulating two-phase flows of complex fluids
• A second-order energy stable BDF numerical scheme for the Cahn-Hilliard equation
• Convergence analysis of the fast subspace descent method for convex optimization problems
• Convergence of the Cahn-Hilliard equation to the Hele-Shaw model
• Finite Element Approximation of the Cahn--Hilliard Equation with Degenerate Mobility
• Free energy of a nonuniform system. I: Interfacial free energy
• Front migration in the nonlinear Cahn-Hilliard equation
• High accuracy solutions to energy gradient flows from material science models
• scientific article; zbMATH DE number 3554969 (Why is no real title available?)
• scientific article; zbMATH DE number 940566 (Why is no real title available?)
• Implicit-Explicit Methods for Time-Dependent Partial Differential Equations
• Numerical Analysis of a Continuum Model of Phase Transition
• Numerical analysis of a model for phase separation of a multi-component alloy
• Numerical analysis of the Cahn-Hilliard equation and approximation for the Hele-Shaw problem
• On second order semi-implicit Fourier spectral methods for 2D Cahn-Hilliard equations
• On the stability and accuracy of partially and fully implicit schemes for phase field modeling
• On the stabilization size of semi-implicit Fourier-spectral methods for 3D Cahn-Hilliard equations
• On the unconditionally gradient stable scheme for the Cahn-Hilliard equation and its implementation with Fourier method
• Phase transitions, minimal surfaces and a conjecture of de Giorgi
• Preconditioned steepest descent methods for some nonlinear elliptic equations involving p-Laplacian terms
• Spectral and dynamical stability of nonlinear waves
• Stability and convergence of a second-order mixed finite element method for the Cahn-Hilliard equation
• Stabilized linear semi-implicit schemes for the nonlocal Cahn-Hilliard equation
• The scalar auxiliary variable (SAV) approach for gradient flows
• Unconditionally stable finite difference, nonlinear multigrid simulation of the Cahn-Hilliard-Hele-Shaw system of equations

Cited in

(7)

• Solving phase-field models in the tensor train format to generate microstructures of bicontinuous composites
• Analysis of time-stepping methods for the monodomain model
• A linear second-order in time unconditionally energy stable finite element scheme for a Cahn-Hilliard phase-field model for two-phase incompressible flow of variable densities
• Asymptotic stability of many numerical schemes for phase-field modeling
• Second order energy stable semi-implicit schemes for the 2D Allen-Cahn equation
• Benchmark computation of morphological complexity in the functionalized Cahn-Hilliard gradient flow
• Energy stable semi-implicit schemes for the 2D Allen-Cahn and fractional Cahn-Hilliard equations

Describes a project that uses

Uses Software

• RODAS

This page was built for publication: Asymptotic behaviour of time stepping methods for phase field models

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2027928)