Abstract:We introduce a high-order space-time approximation of the Shallow Water Equations with sources that is invariant-domain preserving (IDP) and well-balanced with respect to rest states. The employed time-stepping technique is a novel explicit Runge-Kutta (ERK) approach which is an extension of the class of ERK-IDP methods introduced by Ern and Guermond (SIAM J. Sci. Comput. 44(5), A3366--A3392, 2022) for systems of non-linear conservation equations. The resulting method is then numerically illustrated through verification and validation.
What problem does this paper attempt to address?
This paper attempts to solve several key problems in the numerical simulation of the Shallow Water Equations (SWE), specifically including:
1. **High - order accuracy**: Many existing methods can only achieve relatively low accuracy in time and space. In order to improve the accuracy of numerical solutions, it is very important to maintain high - order accuracy especially in long - time simulations.
2. **Well - balanced**: For rest states, that is, when the fluid is at rest, the numerical method should be able to accurately maintain this equilibrium state. This requires that the method can handle terrain and other external source terms without introducing unnecessary perturbations.
3. **Invariant - domain - preserving (IDP)**: Ensure that the numerical solution is always within a physically reasonable range, for example, the water depth \(h\) is always positive, avoiding negative values or unrealistic results.
4. **Robustness with respect to dry states**: When dealing with the wet - dry interface area, the numerical method needs to have good robustness to prevent problems such as division by zero from occurring.
### Specific problem description
The paper aims to develop a new explicit Runge - Kutta (ERK) method for solving the shallow water equations including terrain and external source terms. This method should have the following characteristics:
- **High - order accuracy**: It has high - order accuracy in both time and space.
- **Well - balanced**: It remains balanced for rest states, even in the presence of terrain changes or other external source terms.
- **Invariant - domain - preserving**: Ensure that the numerical solution always remains within a physically reasonable range.
- **Robustness with respect to dry states**: It can stably handle the wet - dry interface area.
### Mathematical model
The shallow water equations can be expressed as:
\[
\frac{\partial u}{\partial t}+\nabla\cdot f(u) = b(u, z(x))
\]
where:
- \(u:=(h, q)^T\in\mathbb{R}^{d + 1}\), \(h(x, t)\) is the water depth, and \(q(x, t)\in\mathbb{R}^d\) is the average momentum vector.
- \(z(x)\) is a known terrain function.
- \(f(u):=\left(q, q\otimes v+\frac{1}{2}gh^2I_d\right)^T\), where \(v := h^{-1}q\) is the average velocity vector, \(g\) is the gravitational acceleration constant, and \(I_d\in\mathbb{R}^{d\times d}\) is the identity matrix.
- \(b(u, z(x)):=\left(0, -gh\nabla z\right)^T\) is the source term.
### Solution
The paper proposes a new high - order explicit Runge - Kutta method, combined with convex limiting technology, to ensure that the method maintains high - order accuracy while also having well - balanced and invariant - domain - preserving properties. In addition, this method further improves the robustness with respect to dry states by introducing modified auxiliary states and local boundary control.
### Conclusion
Through these improvements, the paper provides an efficient, accurate, and robust numerical solution method for the shallow water equations suitable for complex terrains and external source terms, especially suitable for long - time simulations and multi - physics field extension applications.