612 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
Construction of Modern Robust Nodal Discontinuous Galerkin Spectral Element Methods for the Compressible Navier-Stokes Equations
Discontinuous Galerkin (DG) methods have a long history in computational
physics and engineering to approximate solutions of partial differential
equations due to their high-order accuracy and geometric flexibility. However,
DG is not perfect and there remain some issues. Concerning robustness, DG has
undergone an extensive transformation over the past seven years into its modern
form that provides statements on solution boundedness for linear and nonlinear
problems.
This chapter takes a constructive approach to introduce a modern incarnation
of the DG spectral element method for the compressible Navier-Stokes equations
in a three-dimensional curvilinear context. The groundwork of the numerical
scheme comes from classic principles of spectral methods including polynomial
approximations and Gauss-type quadratures. We identify aliasing as one
underlying cause of the robustness issues for classical DG spectral methods.
Removing said aliasing errors requires a particular differentiation matrix and
careful discretization of the advective flux terms in the governing equations.Comment: 85 pages, 2 figures, book chapte
A free-energy stable nodal discontinuous Galerkin approximation with summation-by-parts property for the Cahn-Hilliard equation
We present a nodal Discontinuous Galerkin (DG) scheme for the Cahn-Hilliard
equation that satisfies the summation-by-parts simultaneous-approximation-term
(SBP-SAT) property. The latter permits us to show that the discrete free-energy
is bounded, and as a result, the scheme is provably stable. The scheme and the
stability proof are presented for general curvilinear three-dimensional
hexahedral meshes. We use the Bassi-Rebay 1 (BR1) scheme to compute interface
fluxes, and an IMplicit-EXplicit (IMEX) scheme to integrate in time. Lastly, we
test the theoretical findings numerically and present examples for two and
three-dimensional problems
An entropy stable discontinuous Galerkin method for the shallow water equations on curvilinear meshes with wet/dry fronts accelerated by GPUs
We extend the entropy stable high order nodal discontinuous Galerkin spectral
element approximation for the non-linear two dimensional shallow water
equations presented by Wintermeyer et al. [N. Wintermeyer, A. R. Winters, G. J.
Gassner, and D. A. Kopriva. An entropy stable nodal discontinuous Galerkin
method for the two dimensional shallow water equations on unstructured
curvilinear meshes with discontinuous bathymetry. Journal of Computational
Physics, 340:200-242, 2017] with a shock capturing technique and a positivity
preservation capability to handle dry areas. The scheme preserves the entropy
inequality, is well-balanced and works on unstructured, possibly curved,
quadrilateral meshes. For the shock capturing, we introduce an artificial
viscosity to the equations and prove that the numerical scheme remains entropy
stable. We add a positivity preserving limiter to guarantee non-negative water
heights as long as the mean water height is non-negative. We prove that
non-negative mean water heights are guaranteed under a certain additional time
step restriction for the entropy stable numerical interface flux. We implement
the method on GPU architectures using the abstract language OCCA, a unified
approach to multi-threading languages. We show that the entropy stable scheme
is well suited to GPUs as the necessary extra calculations do not negatively
impact the runtime up to reasonably high polynomial degrees (around ). We
provide numerical examples that challenge the shock capturing and positivity
properties of our scheme to verify our theoretical findings
- …