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, J. Sci. Comput., 44 (2010), pp.136--155; Carpenter et al., SIAM J. Sci. Comput., 36 (2014), pp.~B835-B867). 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 Chen and Shu (J. Comput. Phys., 345 (2017), pp.~427--461) 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. (SIAM J. Numer. Anal., 51 (2013), pp.~1371--1391), while entropy stable numerical fluxes are used at element interfaces. Time discretization is performed with strong-stability preserving Runge-Kutta schemes. We use this framework for the discretization of two systems in one space-dimension: a $2\times 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.

References in corpus (3)

• Split Form Nodal Discontinuous Galerkin Schemes with Summation-By-Parts Property for the Compressible Euler Equations
• A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes
• Shallow water equations: Split-form, entropy stable, well-balanced, and positivity preserving numerical methods

Cited by in corpus (12)

• Subcell limiting strategies for discontinuous Galerkin spectral element methods
• Entropy stable, robust and high-order DGSEM for the compressible multicomponent Euler equations
• High order entropy preserving ADER scheme
• An Entropy Stable Nodal Discontinuous Galerkin Method for the resistive MHD Equations. Part II: Subcell Finite Volume Shock Capturing
• An entropy stable high-order discontinuous Galerkin spectral element method for the Baer-Nunziato two-phase flow model
• Entropy Stable Discontinuous Galerkin Methods for Balance Laws in Non-Conservative Form: Applications to the Euler Equations with Gravity
• A positivity preserving strategy for entropy stable discontinuous Galerkin discretizations of the compressible Euler and Navier-Stokes equations
• Entropy-Stable Gauss Collocation Methods for Ideal Magneto-Hydrodynamics
• A discontinuous Galerkin spectral element method for a nonconservative compressible multicomponent flow model
• Generalized upwind summation-by-parts operators and their application to nodal discontinuous Galerkin methods
• An entropy stable discontinuous Galerkin method for the two-layer shallow water equations on curvilinear meshes
• Acoustic Propagation/Refraction Through Diffuse Interface Models