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.

On a positive-preserving, energy-stable numerical scheme to mass-action kinetics with detailed balance

Explicit exponential convergence to equilibrium for nonlinear reaction-diffusion systems with detailed balance condition

A Second-Order Accurate, Operator Splitting Scheme for Reaction-Diffusion Systems in an Energetic Variational Formulation

Existence of Positive Steady States for Mass Conserving and Mass-Action Chemical Reaction Networks with a Single Terminal-Linkage Class

Efficient bound preserving and asymptotic preserving semi-implicit schemes for the fast reaction-diffusion system

Convergence analysis of a positivity-preserving numerical scheme for the Cahn-Hilliard-Stokes system with Flory-Huggins energy potential

A second-order accurate, original energy dissipative numerical scheme for chemotaxis and its convergence analysis

Persistence of Delayed Complex Balanced Chemical Reaction Networks

A Decomposition-Based Approach for Deriving Positive Steady States of a Class of Chemical Reaction Networks with Non-Mass-Action Kinetics

A Second Order Accurate, Positivity Preserving Numerical Method for the Poisson–Nernst–Planck System and Its Convergence Analysis

Global stability of first order endotactic reaction systems

Second Order Unconditional Positive Preserving Schemes for Non-equilibrium Reactive Flows with Mass and Mole Balance

Unconditionally positivity preserving and energy dissipative schemes for Poisson--Nernst--Planck equations

A Graphic Formulation of Nonisothermal Chemical Reaction Systems and the Analysis of Detailed Balanced Networks.

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

An overview of computational methods for chemical equilibrium and kinetic calculations for geochemical and reactive transport modeling

A Kinetic Scheme Based On Positivity Preservation For Multi-component Euler Equations

Trend to equilibrium for reaction-diffusion systems arising from complex balanced chemical reaction networks

A Balanced Shaping Method for Stabilizing Chemical Reaction Networks

Sufficient conditions for linear stability of complex-balanced equilibria in generalized mass-action systems