Abstract:Computing equilibrium shapes of crystals (ESC) is a challenging problem in materials science that involves minimizing an orientation-dependent (i.e., anisotropic) surface energy functional subject to a prescribed mass constraint. The highly nonlinear and singular anisotropic terms in the problem make it very challenging from both the analytical and numerical aspects. Especially, when the strength of anisotropy is very strong (i.e., strongly anisotropic cases), the ESC will form some singular, sharp corners even if the surface energy function is smooth. Traditional numerical approaches, such as the $H^{-1}$ gradient flow, are unable to produce true sharp corners due to the necessary addition of a high-order regularization term that penalizes sharp corners and rounds them off. In this paper, we propose a new numerical method based on the Davis-Yin splitting (DYS) optimization algorithm to predict the ESC instead of using gradient flow approaches. We discretize the infinite-dimensional phase-field energy functional in the absence of regularization terms and transform it into a finite-dimensional constraint minimization problem. The resulting optimization problem is solved using the DYS method which automatically guarantees the mass-conservation and bound-preserving properties. We also prove the global convergence of the proposed algorithm. These desired properties are numerically observed. In particular, the proposed method can produce real sharp corners with satisfactory accuracy. Finally, we present numerous numerical results to demonstrate that the ESC can be well simulated under different types of anisotropic surface energies, which also confirms the effectiveness and efficiency of the proposed method.

Efficient and unconditionally energy stable exponential-SAV schemes for the phase field crystal equation

A Second Order Accurate Scalar Auxiliary Variable (SAV) Numerical Method for the Square Phase Field Crystal Equation

High-order energy stable schemes of incommensurate phase-field crystal model

A Second-Order, Global-in-Time Energy Stable Implicit-Explicit Runge–Kutta Scheme for the Phase Field Crystal Equation

Parallel Energy-Stable Phase Field Crystal Simulations Based on Domain Decomposition Methods.

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

Unconditionally energy stable ESAV-VEM schemes with variable time steps for the time fractional Allen-Cahn equation

New third-order convex splitting methods and analysis for the phase field crystal equation

Global-in-time energy stability analysis for the exponential time differencing Runge-Kutta scheme for the phase field crystal equation

Efficiently linear and unconditionally energy-stable time-marching schemes with energy relaxation for the phase-field surfactant model

Efficient and accurate exponential SAV algorithms with relaxation for dissipative system

An Unconditionally Energy Stable Method for the Anisotropic Phase-Field Crystal Model in Two Dimension

An efficiently linear and totally decoupled variant of SAV approach for the binary phase-field surfactant fluid model

The high-order exponential semi-implicit scalar auxiliary variable approach for nonlocal Cahn-Hilliard equation

A fully-decoupled second-order-in-time and unconditionally energy stable scheme for a phase-field model of two phase flow with variable density

An operator-splitting optimization approach for phase-field simulation of equilibrium shapes of crystals

A unified framework on the original energy laws of three effective classes of Runge-Kutta methods for phase field crystal type models

Arbitrary high-order linearly implicit energy-preserving algorithms for Hamiltonian PDEs

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

A Second-Order Crank-Nicolson Leap-Frog Scheme for the Modified Phase Field Crystal Model with Long-Range Interaction

The stabilized exponential-SAV approach preserving maximum bound principle for nonlocal Allen-Cahn equation