Abstract:We present a framework for devising discretely energy-stable schemes for general dissipative systems based on a generalized auxiliary variable. The auxiliary variable, a scalar number, can be defined in terms of the energy functional by a general class of functions, not limited to the square root function adopted in previous approaches. The current method has another remarkable property: the computed values for the generalized auxiliary variable are guaranteed to be positive on the discrete level, regardless of the time step sizes or the external forces. This property of guaranteed positivity is not available in previous approaches. A unified procedure for treating the dissipative governing equations and the generalized auxiliary variable on the discrete level has been presented. The discrete energy stability of the proposed numerical scheme and the positivity of the computed auxiliary variable have been proved for general dissipative systems. The current method, termed gPAV (generalized Positive Auxiliary Variable), requires only the solution of linear algebraic equations within a time step. With appropriate choice of the operator in the algorithm, the resultant linear algebraic systems upon discretization involve only constant and time-independent coefficient matrices, which only need to be computed once and can be pre-computed. Several specific dissipative systems are studied in relative detail using the gPAV framework. Ample numerical experiments are presented to demonstrate the performance of the method, and the robustness of the scheme at large time step sizes.

gPAV-based unconditionally energy-stable schemes for the Cahn–Hilliard equation: Stability and error analysis

A new decoupled unconditionally stable scheme and its optimal error analysis for the Cahn-Hilliard-Navier-Stokes equations

An efficient data assimilation based unconditionally stable scheme for Cahn–Hilliard equation

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

Error Analysis with Polynomial Dependence on in SAV Methods for the Cahn-Hilliard Equation

Highly efficient, robust and unconditionally energy stable second order schemes for approximating the Cahn-Hilliard-Brinkman system

An exponential-free Runge--Kutta framework for developing third-order unconditionally energy stable schemes for the Cahn--Hilliard equation

Global error estimates of high-order fully decoupled schemes for the Cahn-Hilliard-Navier-Stokes model of Two-Phase Incompressible Flows

A General Class of Linear Unconditionally Energy Stable Schemes for the Gradient Flows

A Roadmap for Discretely Energy-Stable Schemes for Dissipative Systems Based on a Generalized Auxiliary Variable with Guaranteed Positivity

High-order energy stable variable-step schemes for the time-fractional Cahn-Hilliard model

A high-order accurate unconditionally stable bound-preserving numerical scheme for the Cahn-Hilliard-Navier-Stokes equations

High-efficiency and positivity-preserving stabilized SAV methods for gradient flows

Variable-step numerical schemes and energy dissipation laws for time fractional Cahn–Hilliard model

Efficiently high-order time-stepping R-GSAV schemes for the Navier–Stokes-Poisson-Nernst–Planck equations

An energy stable fourth order finite difference scheme for the Cahn-Hilliard equation

Second-Order Decoupled Linear Energy-Law Preserving gPAV Numerical Schemes for Two-Phase Flows in Superposed Free Flow and Porous Media

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

Second-order energy-stable scheme and superconvergence for the finite difference method on non-uniform grids for the viscous Cahn–Hilliard equation

Parallel Energy-Stable Solver for a Coupled Allen--Cahn and Cahn--Hilliard System

Efficient, decoupled, and second-order unconditionally energy stable numerical schemes for the coupled Cahn–Hilliard system in copolymer/homopolymer mixtures