Abstract:In this paper, we develop a fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme for simulating the multicomponent, chemically reacting, compressible Euler equations with complex thermodynamics. The proposed formulation is an extension of the fully conservative, high-order numerical method previously developed by Johnson and Kercher (2020) [14] that maintains pressure equilibrium between adjacent elements. In this first part of our two-part paper, we focus on the one-dimensional case. Our methodology is rooted in the minimum entropy principle satisfied by entropy solutions to the multicomponent, compressible Euler equations, which was proved by Gouasmi et al. (2020) [16] for nonreacting flows. We first show that the minimum entropy principle holds in the reacting case as well. Next, we introduce the ingredients, including a simple linear-scaling limiter, required for the discrete solution to have nonnegative species concentrations, positive density, positive pressure, and bounded entropy. We also discuss how to retain the aforementioned ability to preserve pressure equilibrium between elements. Operator splitting is employed to handle stiff chemical reactions. To guarantee discrete satisfaction of the minimum entropy principle in the reaction step, we develop an entropy-stable discontinuous Galerkin method based on diagonal-norm summation-by-parts operators for solving ordinary differential equations. The developed formulation is used to compute canonical one-dimensional test cases, namely thermal-bubble advection, advection of a low-density Gaussian wave, multicomponent shock-tube flow, and a moving hydrogen-oxygen detonation wave with detailed chemistry. We demonstrate that the formulation can achieve optimal high-order convergence in smooth flows. Furthermore, we find that the enforcement of an entropy bound can considerably reduce the large-scale nonlinear instabilities that emerge when only the positivity property is enforced, to an even greater extent than in the monocomponent, calorically perfect case. Finally, mass, total energy, and atomic elements are shown to be discretely conserved.

A positivity preserving strategy for entropy stable discontinuous Galerkin discretizations of the compressible Euler and Navier-Stokes equations

Entropy stable discontinuous Galerkin methods for the shallow water equations with subcell positivity preservation

Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Navier-Stokes equations

Entropy-stable hybridized discontinuous Galerkin methods for the compressible Euler and Navier-Stokes equations

A general positivity-preserving algorithm for implicit high-order finite volume schemes solving the Euler and Navier-Stokes equations

Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part I: The one-dimensional case

Entropy stabilization and property-preserving limiters for discontinuous Galerkin discretizations of nonlinear hyperbolic equations

A positivity-preserving implicit-explicit scheme with high order polynomial basis for compressible Navier-Stokes equations

Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part II: The multidimensional case

First-order positivity-preserving entropy stable spectral collocation scheme for the 3-D compressible Navier-Stokes equations

Uniformly High-Order Structure-Preserving Discontinuous Galerkin Methods for Euler Equations with Gravitation: Positivity and Well-Balancedness

Positivity-preserving entropy filtering for the ideal magnetohydrodynamics equations

High order entropy stable discontinuous Galerkin spectral element methods through subcell limiting

The Entropy Satisfying Discontinuous Galerkin Method for Fokker–Planck Equations

Positivity-preserving discontinuous spectral element methods for compressible multi-species flows

Outflow Positivity Limiting for Hyperbolic Conservation Laws. Part I: Framework and Recipe

Bound-preserving and entropy stable enriched Galerkin methods for nonlinear hyperbolic equations

High order entropy stable schemes for the quasi-one-dimensional shallow water and compressible Euler equations

Positivity-preserving hybrid DG/FV method for compressible Euler equations with stiff source terms