  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.
