A discontinuous Galerkin spectral element method for a nonconservative compressible multicomponent flow model

TL;DR

This paper develops a high-order, stable discontinuous Galerkin spectral element method for simulating nonconservative multicomponent compressible flows with curved unstructured meshes, ensuring accuracy, robustness, and entropy stability.

Contribution

It introduces a novel DGSEM discretization for multicomponent flows, including a root-free HLLC solver and stability proofs, extending previous frameworks to complex geometries.

Findings

The scheme is high-order accurate and free-stream preserving.

It is robust and entropy stable across material interfaces.

Numerical experiments confirm the scheme's effectiveness in 1D and 2D.

Abstract

In this work, we propose an accurate, robust, and stable discretization of the gamma-based compressible multicomponent model by Shyue [J. Comput. Phys., 142 (1998), 208-242] where each component follows a stiffened gas equation of state (EOS). We here extend the framework proposed in Renac [J. Comput. Phys. 382 (2019), 1-26] and Coquel et al. [J. Comput. Phys. 431 (2021) 110135] for the discretization of hyperbolic systems, with both fluxes and nonconservative products, to unstructured meshes with curved elements in multiple space dimensions. The framework relies on the discontinuous Galerkin spectral element method (DGSEM) using collocation of quadrature and interpolation points. We modify the integrals over discretization elements where we replace the physical fluxes and nonconservative products by two-point numerical fluctuations. The contributions of this work are threefold. First, we analyze the semi-discrete DGSEM discretization and prove that the scheme is high-order accurate, free-stream preserving, and entropy stable when excluding material interfaces. Second, we design a three-point scheme with a HLLC solver that does not require a root-finding algorithm for approximating the nonconservative products. The scheme is proved to be robust and entropy stable for convex entropies, preserves uniform states across material interfaces, satisfies a discrete minimum principle on the specific entropy and maximum principles on the EOS parameters. Third, the HLLC solver is applied at interfaces in the DGSEM scheme, while we consider two kinds of fluctuations in the integrals over discretization elements: material interface preserving and entropy conservative. Time integration is performed using SSP Runge-Kutta schemes. The high-order accuracy, nonlinear stability, and robustness of the present scheme are assessed through several numerical experiments in one and two space dimensions.