12,292 research outputs found
On discretely entropy conservative and entropy stable discontinuous Galerkin methods
High order methods based on diagonal-norm summation by parts operators can be
shown to satisfy a discrete conservation or dissipation of entropy for
nonlinear systems of hyperbolic PDEs. These methods can also be interpreted as
nodal discontinuous Galerkin methods with diagonal mass matrices. In this work,
we describe how use flux differencing, quadrature-based projections, and
SBP-like operators to construct discretely entropy conservative schemes for DG
methods under more arbitrary choices of volume and surface quadrature rules.
The resulting methods are semi-discretely entropy conservative or entropy
stable with respect to the volume quadrature rule used. Numerical experiments
confirm the stability and high order accuracy of the proposed methods for the
compressible Euler equations in one and two dimensions
An Entropy Stable Nodal Discontinuous Galerkin Method for the Two Dimensional Shallow Water Equations on Unstructured Curvilinear Meshes with Discontinuous Bathymetry
We design an arbitrary high-order accurate nodal discontinuous Galerkin
spectral element approximation for the nonlinear two dimensional shallow water
equations with non-constant, possibly discontinuous, bathymetry on
unstructured, possibly curved, quadrilateral meshes. The scheme is derived from
an equivalent flux differencing formulation of the split form of the equations.
We prove that this discretisation exactly preserves the local mass and
momentum. Furthermore, combined with a special numerical interface flux
function, the method exactly preserves the mathematical entropy, which is the
total energy for the shallow water equations. By adding a specific form of
interface dissipation to the baseline entropy conserving scheme we create a
provably entropy stable scheme. That is, the numerical scheme discretely
satisfies the second law of thermodynamics. Finally, with a particular
discretisation of the bathymetry source term we prove that the numerical
approximation is well-balanced. We provide numerical examples that verify the
theoretical findings and furthermore provide an application of the scheme for a
partial break of a curved dam test problem
Entropy stable DGSEM for nonlinear hyperbolic systems in nonconservative form with application to two-phase flows
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
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
On the convergence of a shock capturing discontinuous Galerkin method for nonlinear hyperbolic systems of conservation laws
In this paper, we present a shock capturing discontinuous Galerkin (SC-DG)
method for nonlinear systems of conservation laws in several space dimensions
and analyze its stability and convergence. The scheme is realized as a
space-time formulation in terms of entropy variables using an entropy stable
numerical flux. While being similar to the method proposed in [14], our
approach is new in that we do not use streamline diffusion (SD) stabilization.
It is proved that an artificial-viscosity-based nonlinear shock capturing
mechanism is sufficient to ensure both entropy stability and entropy
consistency, and consequently we establish convergence to an entropy
measure-valued (emv) solution. The result is valid for general systems and
arbitrary order discontinuous Galerkin method.Comment: Comments: Affiliations added Comments: Numerical results added,
shortened proo
- …