🔢 mathematics

Asymptotic-Preserving and Well-Balanced Linearly Implicit IMEX Schemes for the Anelastic Limit of the Isentropic Euler Equations with Gravity

This paper presents and analyzes a class of higher-order, linearly implicit IMEX Runge-Kutta schemes that are both asymptotic-preserving and well-balanced for solving the isentropic Euler equations with gravity, ensuring accuracy in the zero-Mach-number anelastic limit through a novel combination of steady-state penalization, finite-volume reconstruction, and source term discretization.

Original authors: Marco Artiano, Hendrik Ranocha, Saurav Samantaray

Published 2026-04-14

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Marco Artiano, Hendrik Ranocha, Saurav Samantaray

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 the weather or the flow of air around an airplane wing using a computer. The math behind this (the Euler equations) is incredibly complex because it has to handle two very different speeds at once:

The Fast Stuff: Sound waves moving at supersonic speeds (like a sonic boom).
The Slow Stuff: Gentle breezes or slow-moving clouds (low-speed flows).

The Problem: The "Speed Trap"

When you try to simulate slow-moving air (like a gentle breeze), standard computer methods get stuck in a speed trap.

Think of it like trying to film a slow-motion video of a snail, but your camera is forced to take pictures at the speed of a hummingbird's wings. To keep the math stable, the computer has to take tiny, tiny time steps. It's like trying to walk across a room by taking steps the size of a grain of sand. You'll eventually get there, but it will take forever, and the computer will overheat.

Furthermore, in the real world, gravity creates a "steady state" (like the air pressure in a room staying constant). If your computer simulation is slightly off, it might think the air is suddenly rushing up or down when it should be still. This is like a car with a broken cruise control that thinks it's going downhill when it's actually on flat ground, causing it to speed up uncontrollably.

The Solution: A "Smart" New Algorithm

The authors of this paper, Marco, Hendrik, and Saurav, have built a new type of mathematical "engine" (a numerical scheme) to solve these problems. They call it an Asymptotic-Preserving (AP) and Well-Balanced (WB) scheme.

Here is how they did it, using some everyday analogies:

1. The "Asymptotic-Preserving" Feature (The Chameleon)

Imagine a chameleon that changes its skin to match its environment perfectly.

In the fast world (high speed): The algorithm acts like a standard, robust solver for fast-moving air.
In the slow world (low speed): As the speed drops, the algorithm automatically transforms itself into a different, more efficient engine designed specifically for slow flows.
The Magic: It doesn't need to be told to switch modes. It just knows to change its behavior so that the computer doesn't get stuck taking tiny steps. It remains stable and fast whether the air is screaming past or gently drifting.

2. The "Well-Balanced" Feature (The Perfect Scale)

Imagine a perfectly balanced scale. On one side, you have the weight of the air (pressure). On the other, you have the pull of gravity. In a calm room, these forces cancel each other out perfectly.

Old Schemes: They were like a scale with a tiny, invisible defect. Even if you put nothing on it, the scale would slowly tip over, thinking there was a ghost weight. This caused the simulation to create fake winds and errors.
The New Scheme: The authors built a "perfect scale." They used a clever trick called Penalization.
- Instead of fighting the complex, non-linear math of gravity, they added a "helper" equation that is linear (simpler) but has the exact same answer.
- Think of it like this: If you want to balance a heavy, wobbly rock, you don't try to hold it with your bare hands. You put it on a simple, flat board that you know balances perfectly. The new algorithm uses this "flat board" (the linear balance) to keep the simulation steady, ensuring that if the air is supposed to be still, it stays still.

3. The "Linearly Implicit" Trick (The Shortcut)

Usually, solving these complex equations requires a "Newton Solver," which is like trying to find the exit of a maze by trying every single path until you find the right one. It's accurate but very slow and computationally expensive.

The Innovation: The authors found a way to turn the maze into a straight hallway. By using their "linear balance" trick, they made the math linearly implicit.
The Result: Instead of wandering the maze, the computer can just walk straight to the answer. It's much faster and doesn't require the heavy computing power of the old methods.

Why Does This Matter?

This new method is a game-changer for simulating:

Weather patterns: Where slow winds and gravity interact.
Astrophysics: Modeling stars and gas clouds.
Engineering: Designing more efficient aircraft or ventilation systems.

In a nutshell: The authors created a "smart" computer simulation that doesn't get stuck when things move slowly, doesn't get confused by gravity, and solves the math much faster than before. It's like upgrading from a bicycle with a flat tire to a high-speed maglev train that can handle any terrain without slowing down.

1. Problem Statement

The paper addresses the numerical approximation of the compressible isentropic Euler equations in the presence of a gravitational potential, specifically focusing on the low-Mach number regime (where the Mach number $\varepsilon \to 0$ ).

Governing Equations: The system models fluid flow where the Mach number is scaled by a small parameter $\varepsilon$ . As $\varepsilon \to 0$ , the system converges to an anelastic incompressible limit system, where density variations are driven by the gravitational potential (hydrostatic balance) rather than acoustic waves.
Challenges:
1. Stiffness: Explicit time-stepping schemes require time steps proportional to $\varepsilon$ (CFL condition), making them computationally prohibitive as $\varepsilon \to 0$ .
2. Asymptotic Preserving (AP): Standard schemes often fail to correctly capture the transition to the incompressible limit or become unstable without $\varepsilon$ -dependent constraints.
3. Well-Balancing (WB): In low-Mach regimes, flows are often small perturbations of a hydrostatic steady state. Standard schemes often fail to preserve this equilibrium exactly, leading to spurious oscillations or the amplification of small perturbations, which destroys accuracy.
4. Non-linearity: Implicit treatment of the full non-linear flux terms usually requires expensive Newton solvers.

2. Methodology

The authors propose a novel class of higher-order, linearly implicit, Asymptotic-Preserving (AP), and Well-Balanced (WB) IMEX Runge-Kutta (RK) schemes. The methodology is built on three core pillars:

A. Linear Penalization and Reformulation

To avoid non-linear solvers while maintaining the AP property, the authors introduce a linear penalization technique:

They derive a linear hydrostatic balance equation that shares the same steady-state solution ( $\rho_{eq}$ ) as the original non-linear balance equation.
They add and subtract this linear balance term in the momentum equation.
Key Insight: The non-linear terms are treated explicitly, while the linearized stiff terms (involving the pressure gradient and the penalization source) are treated implicitly. This transforms the implicit step into a linear elliptic problem for the density, avoiding Newton iterations.

B. Asymptotic-Preserving (AP) Time Discretization

First-Order Scheme: A semi-implicit scheme is constructed where mass flux is fully implicit, and momentum fluxes are split. The penalization term ensures that as $\varepsilon \to 0$ , the scheme consistently discretizes the anelastic limit system.
Higher-Order Extension: The first-order scheme is extended to arbitrary orders using Additive IMEX Runge-Kutta (IMEX-RK) frameworks (specifically Type-A and Type-CK schemes).
Strong AP Property: The paper proves that for Type-A IMEX-RK schemes, the method achieves strong asymptotic consistency (converging to the limit system) even for ill-prepared initial data (data not initially satisfying the hydrostatic balance) within two time steps.

C. Well-Balanced (WB) Spatial Discretization

To preserve the hydrostatic equilibrium exactly, the authors employ a specialized Finite Volume (FV) approach:

Hydrostatic Reconstruction: Instead of reconstructing primitive variables directly, they reconstruct the variable $\omega = (p e^{-H}, u)$ , where $H$ is a function of the gravitational potential. This ensures that if the flow is in equilibrium ( $u=0$ ), the reconstructed states at cell interfaces exactly satisfy the steady-state condition.
Source Term Discretization: The gravitational source term and the penalization source term are discretized using a specific balance-preserving formulation (based on the exponential form of the hydrostatic solution).
Fluxes: The explicit part uses a Rusanov-type flux, while the implicit part uses central differences.

3. Key Contributions

Linearly Implicit AP-WB Scheme: The primary contribution is a scheme that is both AP and WB but remains linearly implicit. Unlike previous AP-WB schemes that may require solving non-linear systems, this approach only requires solving a linear elliptic system for density, significantly reducing computational cost.
Theoretical Analysis:
- Proof of Asymptotic Consistency: The scheme converges to the correct anelastic limit system as $\varepsilon \to 0$ .
- Proof of Well-Balancing: The scheme exactly preserves the hydrostatic steady state for any grid resolution.
- Analysis of Ill-Prepared Data: Demonstration that Type-A schemes recover the well-prepared state (and thus the divergence-free condition) after two time steps, even from ill-prepared initial conditions.
Higher-Order Accuracy: The framework allows for arbitrary high-order time accuracy using standard IMEX-RK coefficients (e.g., DP2-A(2,4,2)), maintaining stability independent of $\varepsilon$ .

4. Numerical Results

The authors validate the scheme through several numerical case studies:

Well-Balancing Test: The scheme maintains the hydrostatic steady state with machine precision ( $L_2$ errors $\approx 10^{-30}$ ) across a wide range of Mach numbers ( $\varepsilon \in [1, 10^{-4}]$ ) and different gravitational potentials.
Perturbation of Steady State:
- For small perturbations ( $\zeta = \varepsilon^2$ ), the proposed WB-AP scheme correctly damps the perturbation.
- In contrast, a non-well-balanced (NWB) AP scheme amplifies the perturbation or generates spurious oscillations, leading to solution blow-up for small $\varepsilon$ .
Asymptotic Order of Convergence (AOC): The scheme demonstrates uniform second-order convergence in time and space as $\varepsilon \to 0$ , with time steps independent of $\varepsilon$ .
2D Vortex Benchmark: The scheme simulates a stationary vortex with minimal dissipation. The kinetic energy decay and vorticity profiles remain identical across different Mach numbers ( $\varepsilon = 10^{-1}$ to $10^{-4}$ ), proving uniform low-Mach performance.
Riemann Problem: The scheme accurately captures shock structures in a Riemann problem, showing robustness in the compressible regime as well.

5. Significance

This work bridges a critical gap in computational fluid dynamics for atmospheric and astrophysical flows:

Efficiency: By achieving linear implicitness, it removes the bottleneck of non-linear solvers, making high-order AP schemes feasible for large-scale simulations.
Robustness: The combination of AP and WB properties ensures that the scheme works seamlessly from the fully compressible regime ( $\varepsilon \sim 1$ ) to the incompressible anelastic limit ( $\varepsilon \to 0$ ) without changing the algorithm or requiring mesh refinement.
Generality: The penalization approach provides a template for designing similar schemes for more complex equations of state (e.g., including potential temperature variations), which is identified as future work.

In summary, the paper presents a mathematically rigorous and computationally efficient numerical framework that solves the long-standing challenge of simulating low-Mach number flows with gravity while maintaining accuracy, stability, and equilibrium preservation.