A Conservative Discontinuous Galerkin Discretization for the Chemically Reacting Navier-Stokes Equations

Abstract

We present a detailed description and verification of a discontinuous Galerkin finite element method (DG) for the multi-component chemically reacting compressible Navier-Stokes equations that retains the desirable properties of DG, namely discrete conservation and high-order accuracy in smooth regions of the flow. Pressure equilibrium between adjacent elements is maintained through the consistent evaluation of the thermodynamics model and the resulting weak form, as well as the proper choice of nodal basis. As such, the discretization does not generate unphysical pressure oscillations in smooth regions of the flow or at material interfaces where the temperature is continuous. Additionally, we present an hp-adaptive DG method for solving systems of ordinary differential equations, DGODE, which is used to resolve the temporal evolution of the species concentrations due to stiff chemical reactions. The coupled solver is applied to several challenging test problems including multi-component shocked flows as well as chemically reacting detonations, deflagrations, and shear flows with detailed kinetics. We demonstrate that the discretization does not produce unphysical pressure oscillations and, when applicable, we verify that it maintains discrete conservation. The solver is also shown to reproduce the expected temperature and species profiles throughout a detonation as well as the expected two-dimensional cellular detonation structure. We also demonstrate that the solver can produce accurate, high-order, approximations of temperature and species profiles without artificial stabilization for the case of a one-dimensional pre-mixed flame. Finally, high-order solutions of two- and three-dimensional multi-component chemically reacting shear flows, computed without any additional stabilization, are presented.

arXiv commentDirect expansion of arXiv:1910.10544 Replacement: fixes bad figure reference

References in corpus (6)

• pyJac: analytical Jacobian generator for chemical kinetics
• An optimization-based approach for high-order accurate discretization of conservation laws with discontinuous solutions
• Sparse invariant domain preserving discontinuous Galerkin methods with subcell convex limiting
• A five-equation model for the simulation of miscible and viscous compressible fluids
• An Adaptive Moving Finite Element Method for Steady Low Mach Number Compressible Combustion Problems
• A Conservative Discontinuous Galerkin Discretization for the Total Energy Formulation of the Reacting Navier Stokes Equations

Cited by in corpus (4)

• The Moving Discontinuous Galerkin Finite Element Method with Interface Condition Enforcement for Compressible Viscous Flows
• Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part I: The one-dimensional case
• Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part II: The multidimensional case
• Effect of Inflow Turbulence on Premixed Combustion in a Cavity Flameholder