Abstract:The nonlinear shallow water equations (SWEs) are widely used to model the unsteady water flows in rivers and coastal areas. In this work, we present a novel class of locally conservative, entropy stable and well-balanced discontinuous Galerkin (DG) methods for the nonlinear shallow water equation with a non-flat bottom topography. The major novelty of our work is the use of velocity field as an independent solution unknown in the DG scheme, which is closely related to the entropy variable approach to entropy stable schemes for system of conservation laws proposed by Tadmor (in: Tezduyar, Hughes T (eds) Proceedings of the winter annual meeting of the American Society of Mechanical Engineering 1986) back in 1986, where recall that velocity is part of the entropy variable for the shallow water equations. Due to the use of velocity as an independent solution unknown, no specific numerical quadrature rules are needed to achieve entropy stability of our scheme on general unstructured meshes in two dimensions. The proposed DG semi-discretization is then carefully combined with the classical explicit strong stability preserving Runge–Kutta (SSP–RK) time integrators (Gottlieb et al. in SIAM Rev. 43, 89–112, 2001) to yield a locally conservative, well-balanced, and positivity preserving fully discrete scheme. Here the positivity preservation property is enforced with the help of a simple scaling limiter. In the fully discrete scheme, we re-introduce discharge as an auxiliary unknown variable. In doing so, standard slope limiting procedures can be applied on the conservative variables (water height and discharge) without violating the local conservation property. Here we apply a characteristic-wise TVB limiter (Cockburn and Shu in J Comput Phys 141:199–224, 1998) on the conservative variables using the Fu-Shu troubled cell indicator (Fu and Shu in J Comput Phys 347:305–327, 2017) in each inner stage of the Runge–Kutta time stepping to suppress numerical oscillations. This fully discrete can be readily applied to various SWEs simulations without dry areas where the water height is close to zero . The case with dry areas need further special attention, where the velocity approximation can be unphysically large near cells with a small water height, which may eventually crashes the simulation if no special treatment is used near these cells. Here we propose a simple wetting/drying treatment for the velocity update without violating the local conservation property to enhance the robustness of the overall scheme. One- and two-dimensional numerical experiments are presented to demonstrate the performance of the proposed methods.

The application of DG in shallow water equations

A High-Order Velocity-Based Discontinuous Galerkin Scheme for the Shallow Water Equations: Local Conservation, Entropy Stability, Well-Balanced Property, and Positivity Preservation

Well-balanced path-conservative discontinuous Galerkin methods with equilibrium preserving space for two-layer shallow water equations

An Oscillation-free Discontinuous Galerkin Method for Shallow Water Equations

A Path-Conservative ADER Discontinuous Galerkin Method for Non-Conservative Hyperbolic Systems: Applications to Shallow Water Equations

A high-order well-balanced positivity-preserving moving mesh DG method for the shallow water equations with non-flat bottom topography

An entropy stable discontinuous Galerkin method for the two-layer shallow water equations on curvilinear meshes

An Implicit Quadrature-Free Modal Discontinuous Galerkin (DG) Scheme for Shallow Water Equations

Discontinuous Galerkin method for the diffusive-viscous wave equation

A high-order domain preserving DG method for the two-layer shallow water equations

Ensuring 'well-balanced' shallow water flows via a discontinuous Galerkin finite element method: issues at lowest order

Shallow Water Equations: Magic and Limitations

Conservation and stability in a discontinuous Galerkin method for the vector invariant spherical shallow water equations

Discontinuous Galerkin method for solving wave equations in two-phase and viscoelastic media

A high order semi-Lagrangian discontinuous Galerkin method for the two-dimensional incompressible Euler equations and the guiding center Vlasov model without operator splitting

p-adaptive discontinuous Galerkin method for the shallow water equations on heterogeneous computing architectures

Runge-Kutta Discontinuous Local Evolution Galerkin Methods for the Shallow Water Equations on the Cubed-Sphere Grid

High resolution finite volume method for solving shallow water equations

A well-balanced positivity-preserving quasi-Lagrange moving mesh DG method for the shallow water equations