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

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.

arXiv comment50 pages, 20 figures, 1 table

References in corpus (8)

• Preserving energy resp. dissipation in numerical PDEs using the "Average Vector Field" method
• Numerical Approximation of Incompressible Navier-Stokes Equations Based on an Auxiliary Energy Variable
• Multiphase Flows of N Immiscible Incompressible Fluids: A Reduction-Consistent and Thermodynamically-Consistent Formulation and Associated Algorithm
• Efficient linear, Stabilized, second-order time marching schemes for an anisotropic phase field dendritic crystal growth model
• An Unconditionally Energy-Stable Scheme Based on an Implicit Auxiliary Energy Variable for Incompressible Two-Phase Flows with Different Densities Involving Only Precomputable Coefficient Matrices
• Partitioned AVF methods
• Multiphase Allen-Cahn and Cahn-Hilliard Models and Their Discretizations with the Effect of Pairwise Surface Tensions
• Study of a chemo-repulsion model with quadratic production. Part I: Analysis of the continuous problem and time-discrete numerical schemes

Cited by in corpus (14)

• On Computing the Hyperparameter of Extreme Learning Machines: Algorithm and Application to Computational PDEs, and Comparison with Classical and High-Order Finite Elements
• A highly efficient and accurate exponential semi-implicit scalar auxiliary variable (ESI-SAV) approach for dissipative system
• Numerical Approximation of Partial Differential Equations by a Variable Projection Method with Artificial Neural Networks
• gPAV-Based Unconditionally Energy-Stable Schemes for the Cahn-Hilliard Equation: Stability and Error Analysis
• A gPAV-Based Unconditionally Energy-Stable Scheme for Incompressible Flows with Outflow/Open Boundaries
• The generalized scalar auxiliary variable approach (G-SAV) for gradient flows
• An Energy-Stable Scheme for Incompressible Navier-Stokes Equations with Periodically Updated Coefficient Matrix
• A novel Lagrange Multiplier approach with relaxation for gradient flows
• Supplementary Variable Method for Developing Structure-Preserving Numerical Approximations to Thermodynamically Consistent Partial Differential Equations
• Stability and convergence of relaxed scalar auxiliary variable schemes for Cahn-Hilliard systems with bounded mass source
• A hybrid phase field method for fluid-structure interactions in viscous fluids
• Scalar auxiliary variable approach for conservative/dissipative partial differential equations with unbounded energy
• A convergent stochastic scalar auxiliary variable method
• Thermodynamically Consistent Algorithms for Models of Diblock Copolymer Solutions Interacting with Electric and Magnetic Fields