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.

Solution of Two-Dimensional Riemann Problems of Gas Dynamics by Positive Schemes

Solution of Two-Dimensional Riemann Problems of Gas Dynamics by Positive Schemes

Positive Schemes for Solving Multi-Dimensional Hyperbolic Systems of Conservation Laws Ii

Two-dimensional Riemann Problem for Pressureless Gas Dynamics Equations with Functional Solutions

Two-Dimensional Riemann Problems for Chaplygin Gas

An optimization based limiter for enforcing positivity in a semi-implicit discontinuous Galerkin scheme for compressible Navier-Stokes equations

The global smooth symmetric solution to 2-D full compressible Euler system of Chaplygin gases

Positivity-Preserving Analysis of Explicit and Implicit Lax–Friedrichs Schemes for Compressible Euler Equations

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

Conservative difference schemes for one-dimensional flows of polytropic gas

A Complete Global Solution to the Pressure Gradient Equation

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

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

Global entropy solutions to multi-dimensional isentropic gas dynamics with spherical symmetry

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

A Positivity‐preserving Lagrangian Discontinuous Galerkin Scheme with Exact Riemann Solver for Gas‐water Compressible Flows

Extensions and investigations of space‐time generalized Riemann problems numerical schemes for linear systems of conservation laws with source terms

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

Relaxation schemes for entropy dissipative system of viscous conservation laws