An entropy stable discontinuous Galerkin method for the shallow water equations on curvilinear meshes with wet/dry fronts accelerated by GPUs

Abstract

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 $N=7$). We provide numerical examples that challenge the shock capturing and positivity properties of our scheme to verify our theoretical findings.

References in corpus (5)

• Split Form Nodal Discontinuous Galerkin Schemes with Summation-By-Parts Property for the Compressible Euler Equations
• OCCA: A unified approach to multi-threading languages
• GPU performance analysis of a nodal discontinuous Galerkin method for acoustic and elastic models
• Shallow water equations: Split-form, entropy stable, well-balanced, and positivity preserving numerical methods
• A continuous/discontinuous Galerkin solution of the shallow water equations with dynamic viscosity, high-order wetting and drying, and implicit time integration

Cited by in corpus (7)

• A provably entropy stable subcell shock capturing approach for high order split form DG for the compressible Euler Equations
• A staggered semi-implicit hybrid finite volume / finite element scheme for the shallow water equations at all Froude numbers
• High-order entropy stable discontinuous Galerkin methods for the shallow water equations: curved triangular meshes and GPU acceleration
• Symmetry-preserving finite-difference discretizations of arbitrary order on structured curvilinear staggered grids
• Efficient entropy stable Gauss collocation methods
• Semi-Analytical Solutions of Shallow Water Waves with Idealised Bottom Topographies
• Entropy stable discontinuous Galerkin methods for nonlinear conservation laws on networks and multi-dimensional domains