Formulation of Entropy-Stable schemes for the multicomponent compressible Euler equations
arXiv:1904.00972 · doi:10.1016/j.cma.2020.112912
Abstract
In this work, Entropy-Stable (ES) schemes are formulated for the multicomponent compressible Euler equations. Entropy-conservative (EC) and ES fluxes are derived. Particular attention is paid to the limit case of zero partial densities where the structure required by ES schemes breaks down (the entropy variables are no longer defined). It is shown that while an EC flux is well-defined in this limit, a well-defined upwind ES flux requires appropriately averaged partial densities in the dissipation matrix. A similar result holds for the high-order TecNO reconstruction. However, this does not prevent the numerical solution from developing negative partial densities or internal energy. Numerical experiments were performed on one-dimensional and two-dimensional interface and shock-interface problems. The present scheme exactly preserves stationary interfaces. On moving interfaces, it produces pressure anomalies typically observed with conservative schemes [Karni, \textit{J. Comput. Phys.}, 112 (1994) 1]. We find that these anomalies, which are not present in the single-component case, violate neither entropy stability nor a minimum principle of the specific entropy. Finally, we show that the scheme is able to reproduce the physical mechanisms of the two-dimensional shock-bubble interaction problem [Haas \& Sturtevant, J. Fluid Mech. 181 (1987) 41, Quirk \& Karni, J. Fluid Mech. 318 (1996) 129].
Revised manuscript
References in corpus (11)
- On admissibility criteria for weak solutions of the Euler equations
- Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier-Stokes equations
- A counterexample to well-posedness of entropy solutions to the compressible Euler system
- ENO reconstruction and ENO interpolation are stable
- Fast estimation from above of the maximum wave speed in the Riemann problem for the Euler equations
- A Novel High-Order, Entropy Stable, 3D AMR MHD Solver with Guaranteed Positive Pressure
- Analysis and entropy stability of the line-based discontinuous Galerkin method
- A minimum entropy principle in the compressible multicomponent Euler equations
- Entropy Conservative Schemes and the Receding Flow Problem
- Entropy-stable hybridized discontinuous Galerkin methods for the compressible Euler and Navier-Stokes equations
- On entropy stable temporal fluxes
Cited by in corpus (8)
- Entropy stable, robust and high-order DGSEM for the compressible multicomponent Euler equations
- High-Order Methods for Hypersonic Flows with Strong Shocks and Real Chemistry
- Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part I: The one-dimensional case
- High-order accurate entropy stable adaptive moving mesh finite difference schemes for (multi-component) compressible Euler equations with the stiffened equation of state
- Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part II: The multidimensional case
- A Discontinuous Galerkin Solver in the FLASH Multi-Physics Framework
- Finite-difference compatible entropy-conserving schemes for the compressible Euler equations
- A Kinetic Scheme Based On Positivity Preservation For Multi-component Euler Equations