Entropy stable DGSEM for nonlinear hyperbolic systems in nonconservative form with application to two-phase flows
Abstract
In this work, we consider the discretization of nonlinear hyperbolic systems in nonconservative form with the high-order discontinuous Galerkin spectral element method (DGSEM) based on collocation of quadrature and interpolation points (Kopriva and Gassner (2010) [39]). We present a general framework for the design of such schemes that satisfy a semi-discrete entropy inequality for a given convex entropy function at any approximation order. The framework is closely related to the one introduced for conservation laws by Fisher and Carpenter (2013) [27] and Gassner (2013) [32] and relies on the modification of the integral over discretization elements where we replace the physical fluxes by entropy conservative numerical fluxes from Castro et al. (2013) [12], while entropy stable numerical fluxes are used at element interfaces. Time discretization is performed with strong-stability preserving Runge-Kutta schemes. The present method is first validated on the Burgers and Euler equations in nonconservative form. We then use this framework for the discretization of two-phase flow systems in one space-dimension: a 2 × 2 system with a nonconservative product associated to a linearly-degenerate field for which the DGSEM fails to capture the physically relevant solution, and the isentropic Baer-Nunziato model. For the latter, we derive conditions on the numerical parameters of the discrete scheme to further keep positivity of the partial densities and a maximum principle on the void fractions. Numerical experiments support the conclusions of the present analysis and highlight stability and robustness of the present schemes.
- Publication:
-
Journal of Computational Physics
- Pub Date:
- April 2019
- DOI:
- 10.1016/j.jcp.2018.12.035
- arXiv:
- arXiv:1806.10880
- Bibcode:
- 2019JCoPh.382....1R
- Keywords:
-
- Nonconservative hyperbolic systems;
- Entropy stable schemes;
- Discontinuous Galerkin method;
- Summation-by-parts;
- Two-phase flows;
- Mathematics - Numerical Analysis;
- 65M12;
- 65M70;
- 76T10
- E-Print:
- doi:10.1016/j.jcp.2018.12.035