Energy Stable and Structure-Preserving Schemes for the Stochastic Galerkin Shallow Water Equations

Dihan Dai,Yekaterina Epshteyn,Akil Narayan

2023-10-10

Abstract:The shallow water flow model is widely used to describe water flows in rivers, lakes, and coastal areas. Accounting for uncertainty in the corresponding transport-dominated nonlinear PDE models presents theoretical and numerical challenges that motivate the central advances of this paper. Starting with a spatially one-dimensional hyperbolicity-preserving, positivity-preserving stochastic Galerkin formulation of the parametric/uncertain shallow water equations, we derive an entropy-entropy flux pair for the system. We exploit this entropy-entropy flux pair to construct structure-preserving second-order energy conservative, and first- and second-order energy stable finite volume schemes for the stochastic Galerkin shallow water system. The performance of the methods is illustrated on several numerical experiments.

Numerical Analysis

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to solve several key challenges in the numerical simulation of Stochastic Shallow Water Equations (SSWE). Specifically, the authors focus on how to develop numerical methods that can maintain energy stability and structure preservation to deal with the following problems:

1. **Uncertainty handling**: There are uncertain factors in the shallow - water equation model, such as incomplete information on initial conditions or bottom topography functions. These uncertainties are usually modeled as random variables or processes, causing the solution \(U\) to become a random quantity. Therefore, numerical methods that can effectively handle this uncertainty need to be developed.

2. **Numerical stability**: Traditional numerical methods face many challenges when dealing with nonlinear hyperbolic conservation laws and balance laws, such as the discontinuity of solutions, the non - uniqueness of weak solutions, and the difficulty in implementing implicit time - integration solvers. These problems become more complicated after the introduction of uncertainty.

3. **Physical property preservation**: Numerical schemes should be able to preserve some important physical properties, such as the positivity of water depth, the hyperbolicity of the system, and the dynamic properties of approximate equilibrium states (the so - called "well - balanced" property). These properties are crucial for ensuring the physical meaning of numerical solutions.

4. **Energy stability and conservation**: In order to ensure the long - term stability and accuracy of numerical solutions, numerical schemes that can maintain energy conservation and energy stability at the discrete level need to be developed.

### Main contributions

In response to the above problems, the main contributions of this paper include:

- **Derivation of entropy - entropy flux pairs**: The authors derived the entropy - entropy flux pairs of the one - dimensional hyperbolicity - preserving and positivity - preserving stochastic Galerkin shallow - water equation system in space, which provides a theoretical basis for proposing the entropy - acceptability criterion.

- **Energy - conserving and stable finite - volume schemes**: Using the entropy - entropy flux pairs, the authors designed second - order energy - conserving and first - order and second - order energy - stable finite - volume schemes, and these schemes also have the "well - balanced" property. These schemes are stochastic extensions of deterministic schemes.

- **Numerical experiment verification**: The simulation capabilities of the new scheme are demonstrated through a series of numerical experiments. To the best of the authors' knowledge, this is the first proposed stochastic Galerkin shallow - water equation scheme that has energy stability, the "well - balanced" property, and simultaneously preserves positivity and hyperbolicity.

Through these contributions, this paper provides an effective method for dealing with the shallow - water equation model with uncertainty and ensures the stability and physical consistency of numerical solutions.