An energy-stable convex splitting for the phase-field crystal equation

Abstract

The phase-field crystal equation, a parabolic, sixth-order and nonlinear partial differential equation, has generated considerable interest as a possible solution to problems arising in molecular dynamics. This is because the phase-field crystal model can capture atomic-scale effects at time-scales that are orders of magnitude larger than what molecular dynamics simulations can afford presently. Nonetheless, solving this equation is not a trivial task, as a non-increasing free energy and mass conservation need to be verified for the numerical solution to be valid. This work focuses on these issues, and proposes a novel algorithm that guarantees mass conservation, unconditional energy stability and is second-order accurate in time. This is achieved through a convex-concave splitting of the nonlinearity present in the equation, along with the use of a stabilization term that bounds possible increases in free energy. We present numerical results that validate our mathematical proofs, and show two and three dimensional simulations involving crystal growth that showcase the robustness of the method.

43 pages, 13 figures

References in corpus (2)

• Phase-field crystal modeling of equilibrium bcc-liquid interfaces
• An energy-stable convex splitting for the phase-field crystal equation

Cited by in corpus (9)

• An energy-stable convex splitting for the phase-field crystal equation
• Arbitrarily High-order Unconditionally Energy Stable Schemes for Thermodynamically Consistent Gradient Flow Models
• Coarse-grained modeling of crystals by the amplitude expansion of the phase-field crystal model: an overview
• Energy Exchange Analysis in Droplet Dynamics via the Navier-Stokes-Cahn-Hilliard Model
• Unconditionally stable, second-order accurate schemes for solid state phase transformations driven by mechano-chemical spinodal decomposition
• On the Thermodynamics of the Swift-Hohenberg Theory
• Refined isogeometric analysis for generalized Hermitian eigenproblems
• Optimal spectral approximation of $2n$-order differential operators by mixed isogeometric analysis
• Error analysis of SAV finite element method to phase field crystal model