Fast and Accurate Coarsening Simulation with an Unconditionally Stable Time Step

Abstract

We present Cahn-Hilliard and Allen-Cahn numerical integration algorithms that are unconditionally stable and so provide significantly faster accuracy-controlled simulation. Our stability analysis is based on Eyre's theorem and unconditional von Neumann stability analysis, both of which we present. Numerical tests confirm the accuracy of the von Neumann approach, which is straightforward and should be widely applicable in phase-field modeling. We show that accuracy can be controlled with an unbounded time step Delta-t that grows with time t as Delta-t ~ t^alpha. We develop a classification scheme for the step exponent alpha and demonstrate that a class of simple linear algorithms gives alpha=1/3. For this class the speed up relative to a fixed time step grows with the linear size of the system as N/log N, and we estimate conservatively that an 8192^2 lattice can be integrated 300 times faster than with the Euler method.

14 pages, 6 figures

References in corpus (2)

• Theory of Phase Ordering Kinetics
• Efficient Computation of Dendritic Microstructures using Adaptive Mesh Refinement

Cited by in corpus (19)

• An Efficient Algorithm for Solving the Phase Field Crystal Model
• Consistent Hydrodynamics for Phase Field Crystals
• Eighth-order phase-field-crystal model for two-dimensional crystallization
• Simulation of incompressible two-phase flow in porous media with large timesteps
• Maximally-fast coarsening algorithms
• Controlling the accuracy of unconditionally stable algorithms in Cahn-Hilliard Equation
• Vacancy diffusion in colloidal crystals as determined by dynamical density-functional theory and the phase-field-crystal model
• Glassy phases and driven response of the phase-field-crystal model with random pinning
• Development of Knife-Edge Ridges on Ion-Bombarded Surfaces
• Low Mach Asymptotic Preserving Scheme for the Euler-Korteweg Model
• Theory of Microphase Separation in Elastomers
• Long wavelength properties of phase field crystal models with second order dynamics
• Efficient and stable time integration of Cahn-Hilliard equations: explicit, implicit and explicit iterative schemes
• A Method Based on Total Variation for Network Modularity Optimization using the MBO Scheme
• Navigating the Complex Landscape of Shock Filter Cahn-Hilliard Equation: From Regularized to Entropy Solutions
• The Vortex Kinetics of Conserved and Non-conserved O(n) Models
• Parallel energy-stable solver for a coupled Allen-Cahn and Cahn-Hilliard system
• A Method for Finding Structured Sparse Solutions to Non-negative Least Squares Problems with Applications
• Stable Numerical Integration of an Epitaxial Growth Model with Slope Selection