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

Abstract

In this paper, we develop a fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme for simulating the chemically reacting, compressible Euler equations with complex thermodynamics. The proposed formulation is an extension of the conservative, high-order numerical method previously developed by Johnson and Kercher [J. Comput. Phys., 423 (2020), 109826] 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. [ESAIM: Math. Model. Numer. Anal., 54 (2020), 373--389] for nonreacting flows. We first show that the minimum entropy principle holds in the reacting case as well. Next, we introduce the ingredients required for the 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 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, multicomponent shock-tube flow, and a moving hydrogen-oxygen detonation wave with detailed chemistry. 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.

References in corpus (16)

• Split Form Nodal Discontinuous Galerkin Schemes with Summation-By-Parts Property for the Compressible Euler Equations
• On discretely entropy conservative and entropy stable discontinuous Galerkin methods
• Invariant domain preserving discretization-independent schemes and convex limiting for hyperbolic systems
• An optimization-based approach for high-order accurate discretization of conservation laws with discontinuous solutions
• Multiscale modeling strategy and general theory of non-equilibrium plasma assisted ignition and combustion
• Fast estimation from above of the maximum wave speed in the Riemann problem for the Euler equations
• Positivity-Preserving Entropy-Based Adaptive Filtering for Discontinuous Spectral Element Methods
• A Conservative Discontinuous Galerkin Discretization for the Chemically Reacting Navier-Stokes Equations
• Sparse invariant domain preserving discontinuous Galerkin methods with subcell convex limiting
• Entropy-Bounded Discontinuous Galerkin Scheme for Euler Equations
• Invariant-region-preserving DG methods for multi-dimensional hyperbolic conservation law systems, with an application to compressible Euler equations
• Entropy stable, robust and high-order DGSEM for the compressible multicomponent Euler equations
• Formulation of Entropy-Stable schemes for the multicomponent compressible Euler equations
• A minimum entropy principle in the compressible multicomponent Euler equations
• Second Order Unconditional Positive Preserving Schemes for Non-equilibrium Reactive Flows with Mass and Mole Balance
• An Efficient Algorithm for the Multicomponent Compressible Navier-Stokes Equations in Low- and High-Mach Number Regimes

Cited by in corpus (1)

• Finite-difference compatible entropy-conserving schemes for the compressible Euler equations