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.

Energy-preserving schemes for conservative PDEs based on periodic quasi-interpolation methods

A Linearly-Fitted Conservative (dissipative) Scheme for Efficiently Solving Conservative (dissipative) Nonlinear Wave PDEs

A variant of the discrete gradient method for the solution of the semilinear wave equation under different boundary conditions

General Local Energy-Preserving Integrators for Solving Multi-Symplectic Hamiltonian PDEs.

Arbitrarily high-order energy-preserving schemes for the Camassa-Holm equation

High Resolution Schemes for Conservation Laws and Convection-Diffusion Equations with Varying Time and Space Grids

NUMERICAL ANALYSIS OF AN ENERGY-CONSERVATION SCHEME FOR TWO-DIMENSIONAL HAMILTONIAN WAVE EQUATIONS WITH NEUMANN BOUNDARY CONDITIONS

The Onsager principle and structure preserving numerical schemes

An Energy-Preserving Algorithm for Nonlinear Hamiltonian Wave Equations with Neumann Boundary Conditions

A particle-in-Fourier method with semi-discrete energy conservation for non-periodic boundary conditions

Supplementary Variable Method for Developing Structure-Preserving Numerical Approximations to Thermodynamically Consistent Partial Differential Equations

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

Numerical solutions of the nonlinear wave equations with energy-preserving sixth-order finite difference schemes

The Formulation and Analysis of Energy-Preserving Schemes for Solving High-Dimensional Nonlinear Klein–Gordon Equations

Meshless Conservative Scheme for Multivariate Nonlinear Hamiltonian PDEs.

Energy-dissipative spectral renormalization exponential integrator method for gradient flow problems

An Asymptotic Preserving and Energy Stable Scheme for the Euler-Poisson System in the Quasineutral Limit

Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations

Fully discrete energy-dissipative and conservative discrete gradient particle methods for a class of continuity equations

Conservative Multiquadric Quasi-Interpolation Method for Hamiltonian Wave Equations