Arbitrarily High-order Unconditionally Energy Stable Schemes for Thermodynamically Consistent Gradient Flow Models

Abstract

We present a systematical approach to developing arbitrarily high order, unconditionally energy stable numerical schemes for thermodynamically consistent gradient flow models that satisfy energy dissipation laws. Utilizing the energy quadratization (EQ) method, We formulate the gradient flow model into an equivalent form with a corresponding quadratic free energy functional. Based on the equivalent form with a quadratic energy, we propose two classes of energy stable numerical approximations. In the first approach, we use a prediction-correction strategy to improve the accuracy of linear numerical schemes. In the second approach, we adopt the Gaussian collocation method to discretize the equivalent form with a quadratic energy, arriving at an arbitrarily high-order scheme for gradient flow models. Schemes derived using both approaches are proved rigorously to be unconditionally energy stable. The proposed schemes are then implemented in four gradient flow models numerically to demonstrate their accuracy and effectiveness. Detailed numerical comparisons among these schemes are carried out as well. These numerical strategies are rather general so that they can be readily generalized to solve any thermodynamically consistent PDE models.

References in corpus (7)

• Modeling Elasticity in Crystal Growth
• Modeling Elastic and Plastic Deformations in Non-Equilibrium Processing Using Phase Field Crystals
• Numerical Approximations for a three components Cahn-Hilliard phase-field Model based on the Invariant Energy Quadratization method
• A second order in time, uniquely solvable, unconditionally stable numerical scheme for Cahn-Hilliard-Navier-Stokes equation
• Mass Conservative and Energy Stable Finite Difference Methods for the Quasi-incompressible Navier-Stokes-Cahn-Hilliard system: Primitive Variable and Projection-Type Schemes
• An energy-stable convex splitting for the phase-field crystal equation
• A new class of efficient and robust energy stable schemes for gradient flows

Cited by in corpus (7)

• Arbitrarily High-order Linear Schemes for Gradient Flow Models
• An energy stable $C^0$ finite element scheme for a quasi-incompressible phase-field model of moving contact line with variable density
• A general class of linear unconditionally energy stable schemes for the gradient flows
• Supplementary Variable Method for Developing Structure-Preserving Numerical Approximations to Thermodynamically Consistent Partial Differential Equations
• Energy stable Runge-Kutta discontinuous Galerkin schemes for fourth order gradient flows
• Scalar auxiliary variable approach for conservative/dissipative partial differential equations with unbounded energy
• Energy stable arbitrary order ETD-MS method for gradient flows with Lipschitz nonlinearity