Scalar auxiliary variable approach for conservative/dissipative partial differential equations with unbounded energy

by   Tomoya Kemmochi, et al.

In this paper, we present a novel investigation of the so-called SAV approach, which is a framework to construct linearly implicit geometric numerical integrators for partial differential equations with variational structure. SAV approach was originally proposed for the gradient flows that have lower-bounded nonlinear potentials such as the Allen-Cahn and Cahn-Hilliard equations, and this assumption on the energy was essential. In this paper, we propose a novel approach to address gradient flows with unbounded energy such as the KdV equation by a decomposition of energy functionals. Further, we will show that the equation of the SAV approach, which is a system of equations with scalar auxiliary variables, is expressed as another gradient system that inherits the variational structure of the original system. This expression allows us to construct novel higher-order integrators by a certain class of Runge-Kutta methods. We will propose second and fourth order schemes for conservative systems in our framework and present several numerical examples.


page 1

page 2

page 3

page 4


Highly efficient exponential scalar auxiliary variable approaches with relaxation (RE-SAV) for gradient flows

For the past few years, scalar auxiliary variable (SAV) and SAV-type app...

A highly efficient and accurate exponential semi-implicit scalar auxiliary variable (ESI-SAV) approach for dissipative system

The scalar auxiliary variable (SAV) approach is a very popular and effic...

Local sensitivity analysis of the `Membrane shape equation' derived from the Helfrich energy

The Helfrich energy is commonly used to model the elastic bending energy...

An algebraic multigrid method based on an auxiliary topology with edge matrices

This paper introduces a novel approach to algebraic multigrid methods fo...

Data driven gradient flows

We present a framework enabling variational data assimilation for gradie...

Geometric Quasilinearization Framework for Analysis and Design of Bound-Preserving Schemes

Solutions to many partial differential equations satisfy certain bounds ...

Local Characteristic Decomposition Based Central-Upwind Scheme

We propose novel less diffusive schemes for conservative one- and two-di...