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 stable and convergent finite-difference model which conserves the positivity and the dissipativity of Gibbs’ free energy for a nonlinear combustion equation

Low Mach number lattice Boltzmann model for turbulent combustion: flow in confined geometries

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

A positive and bounded convergent scheme for general space-fractional diffusion-reaction systems with inertial times

Numerical simulation of a fractional Davey–Stewartson model via a conservative scheme: preservation of mass, energy and momenta

A numerical method with properties of consistency in the energy domain for a class of dissipative nonlinear wave equations with applications to a Dirichlet boundary-value problem

Numerical Analysis for a Non-isothermal Incompressible Navier–Stokes–Allen–Cahn System

Quenching Study of Two-Dimensional Fractional Reaction–diffusion Equation from Combustion Process

Entropy conservative discretization of compressible Euler equations with an arbitrary equation of state

A staggered pressure correction numerical scheme to compute a travelling reactive interface in a partially premixed mixture

A Finite Difference Method for a Nonlinear Diffusive Model

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

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

An Entropy Satisfying Conservative Method for the Fokker-Planck Equation of the Finitely Extensible Nonlinear Elastic Dumbbell Model

Supplementary Variable Method for Developing Structure-Preserving Numerical Approximations to Thermodynamically Consistent Partial Differential Equations

Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes/Cahn-Hilliard system

Entropy stable modeling of non-isothermal multi-component diffuse-interface two-phase flows with realistic equations of state

Asymptotically entropy-conservative and kinetic-energy preserving numerical fluxes for compressible Euler equations

Property-preserving numerical approximation of a Cahn-Hilliard-Navier-Stokes model with variable density and degenerate mobility

Optimal temporal decay rates of solutions for combustion of compressible fluids

A dual-pairing summation-by-parts finite difference framework for nonlinear conservation laws