Formulation of Entropy-Stable schemes for the multicomponent compressible Euler equations

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 spurious pressure oscillations typically observed with conservative schemes (Karni, 1994). 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 and Sturtevant, 1987; Quirk and Karni, 1996).