Asymptotic Analysis of Shallow Water Moment Equations

✨

This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how a river flows, or how a tsunami wave crashes onto a shore. To do this, scientists use mathematical models.

The Problem: The "Flat" Map vs. The "3D" Reality

The Old Way (Shallow Water Equations - SWE):
Think of the traditional model (SWE) as looking at a river from a drone flying high above. It sees the water's surface and the average speed of the current. It's like saying, "The whole river is moving at 5 miles per hour."

Pros: It's super fast to calculate.
Cons: It's a bit of a lie. In reality, water near the bottom moves slower because of friction with the riverbed, while water near the top moves faster. The "drone view" misses this vertical detail, leading to errors when the flow gets complicated (like in a dam break or a tsunami).

The Better Way (Shallow Water Moment Equations - SWME):
To fix this, scientists created a more detailed model (SWME). Instead of just one speed, this model tries to describe the shape of the water's speed from bottom to top. It uses a stack of variables (called "moments") to build a 3D profile of the flow.

Pros: It's very accurate. It knows the water is slow at the bottom and fast at the top.
Cons: It's computationally expensive. Imagine trying to solve a puzzle with 100 pieces instead of 2. If you want to simulate a massive ocean, this takes forever on a computer.

The Breakthrough: The "Reduced" Model (RSWME)

The authors of this paper asked a clever question: "What if the water is mostly calm and uniform, but we just need to account for tiny wiggles?"

They realized that in many real-world scenarios (like deep water with high friction), the water wants to settle into a smooth, uniform flow. The complex "wiggles" (the extra variables in the SWME) become very small.

The Analogy: The Heavy Suit vs. The Tuxedo
Imagine the full SWME model is a person wearing a heavy, bulky winter suit with 100 pockets (variables). It's great for extreme cold (complex flows), but it's exhausting to run in.
The SWE model is the same person wearing only a t-shirt. It's easy to run, but you freeze if the wind changes.

The authors developed the RSWME (Reduced Shallow Water Moment Equations).

The Idea: They realized that when the "wind" (friction) is strong, the person doesn't actually need all 100 pockets. They can take off the heavy coat and just wear a tuxedo with two essential pockets.
The Magic: They used a mathematical technique (called asymptotic analysis, similar to how physicists simplify complex equations by ignoring tiny, negligible details) to figure out exactly how to calculate those two pockets based on the physics of the situation.
The Result: They created a new model that has the speed of the t-shirt (SWE) but the accuracy of the heavy coat (SWME) for specific conditions.

How They Tested It

They ran three types of simulations to prove their new model works:

The Sharp Wave: A wave with a very steep, jagged edge.
- Result: The new model was much closer to the "perfect" (but slow) model than the old simple model, especially when the water was deep and friction was high.
The Smooth Sine Wave: A gentle, rolling wave.
- Result: The new model was incredibly accurate, improving precision by up to 88% compared to the old simple model, while running just as fast.
The Square Root Profile: A specific, tricky speed pattern.
- Result: The new model captured the complex vertical shape of the water perfectly.

The Big Win

The most exciting part is the speed.

The old detailed model (SWME) took a long time to run, especially as they added more layers of detail (more "moments").
The new model (RSWME) ran up to 77% faster than the detailed model.
It did this without losing the ability to see the vertical details that the simple model missed.

In a Nutshell

The authors took a complex, slow, but accurate model of water flow and found a mathematical shortcut. This shortcut allows computers to run simulations much faster (saving time and energy) while still capturing the important 3D details of how water moves, provided the water is in a "calm" or "friction-heavy" state. It's like finding a secret tunnel that lets you drive a Ferrari at the speed of a bicycle, but with the handling of a race car.

1. Problem Statement

The Shallow Water Equations (SWE) are widely used for modeling free-surface flows but rely on a depth-averaged velocity assumption. This simplification fails to capture complex vertical velocity structures, leading to inaccuracies in scenarios involving non-uniform profiles (e.g., tsunamis, dam breaks, or flows with significant bottom friction).

To address this, the Shallow Water Moment Equations (SWME) were developed. SWME approximate the vertical velocity profile using a polynomial expansion (Legendre polynomials) with coefficients known as "moment variables." While SWME offer higher accuracy than SWE, they introduce a significant computational burden:

Variable Explosion: The system size scales with the number of moments ( $N$ ), requiring $N+2$ variables (height, mean velocity, and $N$ moments) instead of just 2.
Equilibrium Redundancy: In many physical regimes (specifically those with high viscosity and large slip lengths), the system approaches a "constant-velocity equilibrium" where moment variables vanish. However, the full SWME model still requires solving for all $N$ moment variables even when they are negligible, leading to unnecessary computational cost.

The core problem is to derive a reduced model that retains the accuracy of SWME near these equilibria while significantly reducing the number of variables and computational cost.

2. Methodology

The authors employ asymptotic analysis using a technique analogous to the Chapman–Enskog expansion from kinetic theory.

Scaling Assumptions: The analysis focuses on a "viscous slip" regime where both the viscosity ( $\nu$ ) and the slip length ( $\lambda$ ) are large. Specifically, they introduce a small parameter $\epsilon \ll 1$ and scale the parameters as:
$\nu = \frac{\nu_0}{\epsilon}, \quad \lambda = \frac{\lambda_0}{\epsilon}$
This scaling corresponds to a stable equilibrium where the vertical velocity profile becomes uniform (constant velocity), and the moment variables ( $\alpha_j$ ) approach zero.
Asymptotic Expansion: The moment variables $\alpha_i$ are expanded in powers of $\epsilon$ :
$\alpha_i = \alpha_i^{(0)} + \epsilon \alpha_i^{(1)} + \epsilon^2 \alpha_i^{(2)} + O(\epsilon^3)$
By substituting this expansion into the governing SWME equations and performing an order-of-magnitude analysis, the authors derive explicit closure relations. These relations express the higher-order moment variables as functions of the primary variables (fluid height $h$ and mean velocity $u_m$ ) and their spatial derivatives.
Dimensional Reduction: The derived closure relations are substituted back into the original SWME system. This eliminates the need to solve evolution equations for the moments, resulting in a closed system of only two equations (for $h$ and $u_m$ ). This new system is termed the Reduced Shallow Water Moment Equations (RSWME).
Hyperbolicity Analysis: The authors analyze the eigenvalues of the RSWME system matrix to ensure the model remains hyperbolic (real propagation speeds). They identify conditions under which hyperbolicity breaks down and propose a hyperbolic regularisation term (adding a specific $O(\epsilon^4)$ derivative term) to ensure global hyperbolicity.

3. Key Contributions

Derivation of RSWME: The paper derives a closed-form, dimensionally reduced model (RSWME) for arbitrary numbers of moments ( $N$ ).
Universality for $N \geq 2$ : A significant theoretical finding is that for any $N \geq 2$ , the resulting closed RSWME system is identical. This means that increasing the order of the original SWME expansion does not change the reduced system's structure; the complexity of the original model is effectively collapsed into a single 2-variable system for near-equilibrium flows.
Hyperbolic Regularisation: The authors provide a rigorous analysis of the hyperbolicity of the RSWME and introduce a regularisation term to guarantee real eigenvalues, preventing numerical instabilities.
Boussinesq Coefficient Consistency: The derived RSWME advection terms are shown to be consistent with the Boussinesq coefficient approximation, linking the new model to established non-hydrostatic shallow water theories.

4. Results

Numerical tests were conducted comparing the SWE, the full SWME (reference), and the new RSWME across three scenarios:

Wave with Sharp Height Gradient:
- RSWME outperformed SWE in accuracy, particularly for moderate deviations from equilibrium ( $\epsilon = 0.1$ ).
- For large deviations ( $\epsilon = 1$ ), RSWME showed larger errors in moment reconstruction due to sharp gradients in the closure terms, though it still generally outperformed SWE.
Smooth Sine Wave:
- RSWME demonstrated superior accuracy over SWE for both height and velocity.
- Accuracy Improvement: RSWME reduced relative errors by up to 88% compared to SWE.
- The model remained stable and accurate even for larger deviations from equilibrium in smooth flow scenarios.
Square Root Velocity Profile:
- RSWME successfully reconstructed complex vertical velocity profiles that SWE could not capture.
- Computational Efficiency: The RSWME model solved the system for 2 variables, whereas SWME solved for $N+2$ .
- Runtime Reduction: Compared to the full SWME, RSWME reduced computational time by up to 77% (for $N=6$ ). The overhead compared to SWE was negligible (~5%).

5. Significance

The RSWME framework offers a "best of both worlds" solution for shallow water modeling:

Efficiency: It achieves computational costs comparable to the simple SWE (solving only 2 variables) while avoiding the $O(N)$ scaling of the full SWME.
Accuracy: It significantly outperforms the SWE by capturing vertical velocity profile variations, making it suitable for flows with non-uniform vertical structures where SWE fails.
Scalability: The finding that the reduced system is identical for all $N \geq 2$ implies that users can choose high-order SWME for theoretical accuracy without incurring the computational penalty of solving a large system, provided the flow is near the viscous slip equilibrium.

This work provides a robust mathematical foundation for simulating free-surface flows where vertical structure matters but computational resources are limited, bridging the gap between simple depth-averaged models and complex 3D Navier-Stokes simulations.