24 research outputs found
A Provably Stable Discontinuous Galerkin Spectral Element Approximation for Moving Hexahedral Meshes
We design a novel provably stable discontinuous Galerkin spectral element
(DGSEM) approximation to solve systems of conservation laws on moving domains.
To incorporate the motion of the domain, we use an arbitrary
Lagrangian-Eulerian formulation to map the governing equations to a fixed
reference domain. The approximation is made stable by a discretization of a
skew-symmetric formulation of the problem. We prove that the discrete
approximation is stable, conservative and, for constant coefficient problems,
maintains the free-stream preservation property. We also provide details on how
to add the new skew-symmetric ALE approximation to an existing discontinuous
Galerkin spectral element code. Lastly, we provide numerical support of the
theoretical results
A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes
We are interested in the approximation of a steady hyperbolic problem. In
some cases, the solution can satisfy an additional conservation relation, at
least when it is smooth. This is the case of an entropy. In this paper, we
show, starting from the discretisation of the original PDE, how to construct a
scheme that is consistent with the original PDE and the additional conservation
relation. Since one interesting example is given by the systems endowed by an
entropy, we provide one explicit solution, and show that the accuracy of the
new scheme is at most degraded by one order. In the case of a discontinuous
Galerkin scheme and a Residual distribution scheme, we show how not to degrade
the accuracy. This improves the recent results obtained in [1, 2, 3, 4] in the
sense that no particular constraints are set on quadrature formula and that a
priori maximum accuracy can still be achieved. We study the behavior of the
method on a non linear scalar problem. However, the method is not restricted to
scalar 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
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
Affordable, Entropy Conserving and Entropy Stable Flux Functions for the Ideal MHD Equations
In this work, we design an entropy stable, finite volume approximation for
the ideal magnetohydrodynamics (MHD) equations. The method is novel as we
design an affordable analytical expression of the numerical interface flux
function that discretely preserves the entropy of the system. To guarantee the
discrete conservation of entropy requires the addition of a particular source
term to the ideal MHD system. Exact entropy conserving schemes cannot dissipate
energy at shocks, thus to compute accurate solutions to problems that may
develop shocks, we determine a dissipation term to guarantee entropy stability
for the numerical scheme. Numerical tests are performed to demonstrate the
theoretical findings of entropy conservation and robustness.Comment: arXiv admin note: substantial text overlap with arXiv:1509.06902;
text overlap with arXiv:1007.2606 by other author