Invariant-region-preserving DG methods for multi-dimensional hyperbolic conservation law systems, with an application to compressible Euler equations
arXiv:1804.08829 · doi:10.1016/j.jcp.2018.03.004
Abstract
An invariant-region-preserving (IRP) limiter for multi-dimensional hyperbolic conservation law systems is introduced, as long as the system admits a global invariant region which is a convex set in the phase space. It is shown that the order of approximation accuracy is not destroyed by the IRP limiter, provided the cell average is away from the boundary of the convex set. Moreover, this limiter is explicit, and easy for computer implementation. A generic algorithm incorporating the IRP limiter is presented for high order finite volume type schemes. For arbitrarily high order discontinuous Galerkin (DG) schemes to hyperbolic conservation law systems, sufficient conditions are obtained for cell averages to remain in the invariant region provided the projected one-dimensional system shares the same invariant region as the full multi-dimensional hyperbolic system {does}. The general results are then applied to both one and two dimensional compressible Euler equations so to obtain high order IRP DG schemes. Numerical experiments are provided to validate the proven properties of the IRP limiter and the performance of IRP DG schemes for compressible Euler equations.
33 pages, 8 tables, 5 figures, accepted for publication in Journal of Computational Physics, 2018
References in corpus (1)
Cited by in corpus (10)
- Positivity-preserving third order DG schemes for Poisson--Nernst--Planck equations
- Is the Classic Convex Decomposition Optimal for Bound-Preserving Schemes in Multiple Dimensions?
- Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part I: The one-dimensional case
- Maximum bound principles for a class of semilinear parabolic equations and exponential time differencing schemes
- Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part II: The multidimensional case
- Geometric Quasilinearization Framework for Analysis and Design of Bound-Preserving Schemes
- On the invariant region for compressible Euler equations with a general equation of state
- Minimum Principle on Specific Entropy and High-Order Accurate Invariant Region Preserving Numerical Methods for Relativistic Hydrodynamics
- Analysis of algebraic flux correction for semi-discrete advection problems
- Numerical Study of Non-uniqueness for 2D Compressible Isentropic Euler Equations