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

Thomas Kent,Onno Bokhove

DOI: https://doi.org/10.48550/arXiv.2006.03370

2020-06-05

Abstract:The discontinuous Galerkin finite element method (DGFEM) developed by Rhebergen et al. (2008) offers a robust method for solving systems of nonconservative hyperbolic partial differential equations but, as we show here, does not satisfactorily deal with topography in shallow water flows at lowest order (so-called DG0, or equivalently finite volume). In particular, numerical solutions of the space-DG0 discretised one-dimensional shallow water equations over varying topography are not truly `well-balanced'. A numerical scheme is well-balanced if trivial steady states are satisfied in the numerical solution; in the case of the shallow water equations, initialised rest flow should remain at rest for all times. Whilst the free-surface height and momentum remain constant and zero, respectively, suggesting that the scheme is indeed well-balanced, the fluid depth and topography evolve in time. This is both undesirable and unphysical, leading to incorrect numerical solutions for the fluid depth, and is thus a concern from a predictive modelling perspective. We expose this unsatisfactory issue, both analytically and numerically, and indicate a solution that combines the DGFEM formulation for nonconservative products with a fast and stable well-balanced finite-volume method. This combined scheme bypasses the offending issue and successfully integrates nonconservative hyperbolic shallow water-type models with varying topography at lowest order. We briefly discuss implications for the definition of a well-balanced scheme, and highlight applications when higher-order schemes may not be desired, which give further value to our finding beyond its exposure alone.

Computational Physics,Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to ensure that the numerical scheme can be "well - balanced" at the lowest order (i.e., DG0 or the equivalent finite - volume method) when using the Discontinuous Galerkin Finite Element Method (DGFEM) to solve the shallow - water equations. Specifically, the paper points out that the DGFEM scheme proposed by Rhebergen et al. cannot well maintain the stationary flow state at the lowest order (DG0) when dealing with shallow - water flows with varying topography. This means that although the free - surface height \(h + b\) and momentum \(hu\) remain constant and zero, indicating that the scheme seems to be balanced, the fluid depth \(h\) and topography \(b\) evolve over time, which is an undesirable and unphysical situation.

### Specific problem description
1. **Background**:
- The Shallow Water Equations (SWEs) are a set of hyperbolic partial differential equations (PDEs) used to describe shallow - water flows.
- DGFEM is a powerful numerical method suitable for solving non - conservative hyperbolic PDEs systems.

2. **Core of the problem**:
- At the lowest order (DG0), the DGFEM scheme of Rhebergen et al. cannot correctly handle topographic changes, resulting in the evolution of the fluid depth \(h\) and topography \(b\) over time.
- This evolution is undesirable because it leads to incorrect numerical solutions and affects the accuracy of prediction models.

3. **Specific manifestations**:
- Although the free - surface height \(h + b\) and momentum \(hu\) remain constant and zero, the fluid depth \(h\) and topography \(b\) change over time.
- In numerical simulations, this situation is manifested as the "diffusion" of the topography, which in turn leads to changes in the fluid depth.

4. **Solution**:
- The paper proposes a scheme that combines DGFEM and a fast - stable balanced finite - volume method to overcome this problem.
- Through this combined scheme, non - conservative hyperbolic shallow - water models with varying topography can be successfully solved at the lowest order.

### Mathematical formulation
- **Shallow - water equations**:
\[
\begin{aligned}
&\frac{\partial h}{\partial t}+\frac{\partial (hu)}{\partial x} = 0,\\
&\frac{\partial (hu)}{\partial t}+\frac{\partial \left(hu^{2}+\frac{1}{2}gh^{2}\right)}{\partial x}=-gh\frac{\partial b}{\partial x},\\
&\frac{\partial b}{\partial t} = 0.
\end{aligned}
\]

- **Non - conservative form**:
\[
\frac{\partial \mathbf{U}}{\partial t}+\frac{\partial \mathbf{F}(\mathbf{U})}{\partial x}+\mathbf{G}(\mathbf{U})\frac{\partial \mathbf{U}}{\partial x}=0,
\]
where,
\[
\mathbf{U}=\begin{pmatrix}h\\hu\\b\end{pmatrix},\quad \mathbf{F}(\mathbf{U})=\begin{pmatrix}hu\\hu^{2}+\frac{1}{2}gh^{2}\\0\end{pmatrix},