🔢 mathematics

Asymptotic preserving scheme for the shallow water equations with non-flat bottom topography and Manning friction term

This paper proposes a class of high-order, computationally efficient asymptotic preserving (AP) schemes for shallow water equations with non-flat topography and Manning friction by utilizing a semi-implicit IMEX-RK time discretization and high-order WENO reconstruction, eliminating the need for the costly penalization terms used in previous methods.

Original authors: Guanlan Huang, Sebastiano Boscarino, Tao Xiong

Published 2026-04-28

📖 3 min read🧠 Deep dive

CC BY 4.0

Original authors: Guanlan Huang, Sebastiano Boscarino, Tao Xiong

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ 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 simulate how water flows in a massive, complex river system—one with jagged rocks on the bottom (topography) and thick, muddy silt that slows the water down (friction).

To do this on a computer, scientists use mathematical "engines" called Shallow Water Equations. However, simulating these equations is like trying to drive a car through two completely different worlds at the same time:

The High-Speed World (Convection): When the water is rushing fast, it behaves like a racing car. It’s unpredictable, can create "shocks" (like sudden waves or rapids), and requires a very fast, reactive steering system.
The Slow-Motion World (Diffusion): When the water is thick with silt or moving very slowly, it behaves like a heavy truck driving through deep mud. It’s sluggish, heavy, and moves in a slow, predictable crawl.

The Problem: The "Gearbox" Struggle

In the past, computer models struggled to switch between these two worlds. If you set the "engine" to be fast enough to catch the racing car waves, the "muddy" part of the simulation would become unstable and crash (like a car spinning out on ice). If you set it to be stable enough for the mud, the racing car part would move in slow motion, making the simulation take forever to finish.

Previous researchers tried to fix this by adding a "fake weight" (a mathematical penalty) to the simulation to keep it stable. But this "weight" made the simulation sluggish and computationally expensive—it was like driving a race car with an anchor tied to the back.

The Solution: The "Smart Automatic Transmission"

The authors of this paper have designed a new, high-order "Smart Automatic Transmission" for these equations.

Instead of using an anchor to stay stable, they developed a method called SI-IMEX (Semi-Implicit Implicit-Explicit). Here is how it works using a metaphor:

The Explicit Part (The Reflexes): For the fast-moving, "racing car" parts of the water, the computer uses quick, direct calculations. This allows it to catch sudden waves and splashes instantly.
The Implicit Part (The Foresight): For the "heavy mud" parts (the friction), the computer doesn't just react to what is happening now; it uses a bit of mathematical "foresight" to calculate what the friction will do in the next moment. This keeps the simulation from "spinning out" when the water gets thick and slow.

Why is this a big deal? (The "Three Wins")

The paper proves their new method achieves three major goals:

Asymptotic Preserving (The Seamless Shift): As the water transitions from a rushing river to a slow, muddy crawl, the math doesn't break. It "preserves" the physics perfectly, shifting gears without the computer crashing.
Asymptotically Accurate (The High-Definition View): Even when the water is moving through the "muddy" phase, the simulation stays incredibly sharp and detailed. It doesn't become a blurry, low-resolution mess.
Well-Balanced (The "Still Water" Test): If the water is supposed to be perfectly still (like a calm lake), the math ensures it stays perfectly still. Many other models accidentally create "fake waves" even when the water should be calm, but this model keeps the lake steady.

The Bottom Line

By removing the "mathematical anchor" used in previous studies and replacing it with this "smart transmission," the researchers have created a way to simulate complex water movements that is faster, more stable, and much more accurate than before. It allows scientists to study everything from ocean currents to river flooding with much higher confidence and less computing power.

Technical Summary: Asymptotic Preserving Scheme for the Shallow Water Equations with Non-Flat Bottom Topography and Manning Friction

1. Problem Statement

The paper addresses the numerical simulation of the Shallow Water Equations (SWEs) subject to two complex physical phenomena: non-flat bottom topography and Manning friction.

The system is characterized by a multiscale nature governed by a parameter $\epsilon$ . When $\epsilon$ is large ( $O(1)$ ), the system is in a convection-dominated (hyperbolic) regime, where shocks and discontinuities are prevalent. As $\epsilon \to 0$ , the system transitions into a friction-dominated (diffusive) regime, governed by a nonlinear, degenerate parabolic equation known as the $p$ -Laplacian.

Standard numerical methods face a dilemma:

Explicit methods are simple but suffer from severe stability constraints (stiffness) as $\epsilon \to 0$ , requiring extremely small time steps.
Implicit methods are stable but computationally expensive.
Previous AP (Asymptotic Preserving) methods (specifically the authors' previous work) utilized a "penalization technique" (adding weighted diffusive terms) to ensure stability, but this increased computational cost and slowed convergence in the intermediate regime.

2. Methodology

The authors propose a new class of high-order Semi-Implicit Implicit-Explicit Runge-Kutta (SI-IMEX-RK) schemes. The core innovation is the removal of the penalization term while maintaining the AP property.

Key components of the methodology include:

Time Discretization (SI-IMEX-RK): Instead of additive penalization, the authors use a semi-implicit framework. They treat the stiff terms (friction and pressure gradients) implicitly and the non-stiff terms (convection) explicitly. To ensure high-order accuracy, they employ Stiffly Accurate (SA) DIRK (Diagonally Implicit Runge-Kutta) schemes.
Spatial Discretization:
- First-order derivatives: Discretized using high-order WENO (Weighted Essentially Non-Oscillatory) reconstruction to accurately capture shocks and prevent oscillations.
- Second-order derivatives: Discretized using a central difference scheme to handle the diffusive limit.
Well-Balanced Property: To maintain the "lake at rest" steady state ( $h+B = \text{const}, m=0$ ), the authors implement a specialized flux-splitting technique and a consistent treatment of the bottom topography gradient.
Implicit Treatment of Momentum Norm: A critical technical finding is that the $L_2$ norm of the momentum ( $|m|$ ) in the Manning friction term must be treated implicitly to preserve the AP property.

3. Key Contributions

Elimination of Penalization: Developed a scheme that achieves the AP property without the need for artificial weighted diffusive terms, thereby reducing computational overhead.
High-Order Accuracy: Proved that the scheme is not only AP (consistent with the limit equation) but also Asymptotically Accurate (AA), meaning it maintains its high-order temporal accuracy even as $\epsilon \to 0$ .
Efficiency in Intermediate Regimes: Demonstrated that the SI-IMEX approach is significantly more efficient than previous additive IMEX methods, particularly when transitioning between convection and diffusion.
Theoretical Rigor: Provided formal proofs for the AP and AA properties using Chapman-Enskog expansion and mathematical induction.

4. Results and Validation

The authors validated their scheme through several numerical experiments:

1D Accuracy Tests: Confirmed third-order accuracy for both linear and nonlinear Manning friction across various scales of $\epsilon$ (from $1$ to $10^{-6}$ ).
AP Property Verification: Showed that as $\epsilon$ decreases, the numerical solution converges to the solution of the $p$ -Laplacian limit equation.
Comparison with Additive IMEX: Experiments proved that the SI-IMEX scheme requires significantly fewer iterations to reach convergence, especially in the intermediate and diffusive regimes.
2D Tests: Demonstrated the scheme's ability to handle complex 2D smooth and discontinuous initial conditions, as well as its ability to capture wave propagation around bottom topography (well-balanced property).
Necessity of Implicit $|m|$ : Numerical tests showed that treating $|m|$ explicitly (SI-S2 scheme) causes the scheme to fail in the diffusive limit, confirming that implicit treatment is mandatory for stability and AP preservation.

5. Significance

This work provides a robust, high-order numerical framework for geophysical fluid dynamics. By solving the stiffness problem without the "brute force" of penalization, the scheme allows for efficient, large-scale simulations of ocean and river flows where bottom topography and friction are dominant. It bridges the gap between hyperbolic and parabolic regimes, making it highly suitable for multiscale environmental modeling.