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

by   Kailiang Wu, et al.

Solutions to many partial differential equations satisfy certain bounds or constraints. For example, the density and pressure are positive for equations of fluid dynamics, and in the relativistic case the fluid velocity is upper bounded by the speed of light, etc. As widely realized, it is crucial to develop bound-preserving numerical methods that preserve such intrinsic constraints. Exploring provably bound-preserving schemes has attracted much attention and is actively studied in recent years. This is however still a challenging task for many systems especially those involving nonlinear constraints. Based on some key insights from geometry, we systematically propose an innovative and general framework, referred to as geometric quasilinearization (GQL), which paves a new effective way for studying bound-preserving problems with nonlinear constraints. The essential idea of GQL is to equivalently transfer all nonlinear constraints into linear ones, through properly introducing some free auxiliary variables. We establish the fundamental principle and general theory of GQL via the geometric properties of convex regions, and propose three simple effective methods for constructing GQL. We apply the GQL approach to a variety of partial differential equations, and demonstrate its effectiveness and remarkable advantages for studying bound-preserving schemes, by diverse challenging examples and applications which cannot be easily handled by direct or traditional approaches.


Energy plus maximum bound preserving Runge-Kutta methods for the Allen-Cahn equation

It is difficult to design high order numerical schemes which could prese...

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

In this paper, we present a novel investigation of the so-called SAV app...

Non-Dissipative and Structure-Preserving Emulators via Spherical Optimization

Approximating a function with a finite series, e.g., involving polynomia...

A new Lagrange multiplier approach for constructing structure-preserving schemes, II. bound preserving

In the second part of this series, we use the Lagrange multiplier approa...

Monolithic Multigrid for Magnetohydrodynamics

The magnetohydrodynamics (MHD) equations model a wide range of plasma ph...