🔢 mathematics

Low-Mach-number limit of a compressible two-phase flow system with algebraic closure

This paper rigorously establishes the low-Mach-number limit of a compressible isentropic two-phase flow system with equal pressure and single velocity, proving that weak solutions converge to an incompressible non-homogeneous fluid system where partial densities become constant and volume fractions are transported by the divergence-free velocity field, utilizing a novel relative entropy functional and a specific comparison of density norms.

Original authors: Cassandre Lebot

Published 2026-03-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Cassandre Lebot

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 watching a high-speed video of a chaotic storm inside a jar. Inside this jar, you have two different types of "air" (let's call them Red Mist and Blue Mist) swirling together. They are compressible, meaning they can be squished tight or spread out thin, and they are moving very fast.

This paper is about what happens when we slow down the video so much that the "Red" and "Blue" mists stop acting like squishy, compressible balloons and start behaving like incompressible water.

Here is the breakdown of the story, the problem, and the clever solution the author found.

1. The Setup: Two Fluids, One Dance

In the real world, when you mix two fluids (like oil and water, or steam and air), they usually have their own speeds and pressures. But in this specific mathematical model, the author makes a simplifying assumption: The two fluids are "best friends."

They move at the exact same speed (single velocity).
They share the exact same pressure (equal pressure).

Think of them as two dancers holding hands, spinning at the same speed. Even though they are different materials (one might be heavy, one light), they move as a single unit.

2. The Problem: The "Squishy" vs. The "Stiff"

The paper looks at the Low Mach Number Limit.

High Mach Number: The fluids are moving fast, like a jet plane. They are "squishy" (compressible). If you push them, they compress, and the pressure waves travel fast.
Low Mach Number: The fluids are moving slowly, like a river. They act "stiff" (incompressible). If you push them, they don't squish; they just flow around.

The Challenge: Mathematically, going from "squishy" to "stiff" is like trying to turn a rubber ball into a steel ball. As the speed slows down, the equations become "singular"—they break or become infinitely difficult to solve because the pressure term explodes.

The author wants to prove rigorously that if you start with a fast, squishy mix and slow it down, it smoothly transforms into a slow, stiff, incompressible flow without the math breaking.

3. The Magic Tool: The "Relative Entropy" Scale

To prove this, the author invents a special measuring tool called Relative Entropy.

Imagine you have a Target State (the perfect, slow, incompressible flow you want to reach) and a Current State (the messy, fast, compressible flow you have right now).

The author creates a "Scorecard" (the Relative Entropy) that measures the distance between the Current State and the Target State.

If the score is 0, the two states are identical.
If the score is High, they are very different.

The goal is to show that as the speed slows down (the "Mach number" goes to zero), this Scorecard naturally drops to zero.

4. The Twist: The "Partial Densities" Puzzle

Here is where it gets tricky. In a single fluid, you just track one density. But here, we have Red Mist and Blue Mist.

The Total Density (Red + Blue) changes as they mix.
The Partial Densities (just Red, just Blue) are what we really care about.

The author discovered a clever trick: The "Volume Fraction" is the key.
Think of the jar as being filled with a certain percentage of Red and Blue.

As the flow slows down, the amount of Red and Blue in the jar stays constant (mass is conserved).
However, the density of the individual mists (how tightly packed the Red molecules are) settles into a constant value.
The Volume Fraction (the percentage of the jar filled with Red vs. Blue) doesn't become constant; instead, it gets carried along by the flow like a leaf on a river.

The author proves that even though the math is messy, the "Red" and "Blue" densities settle down to constant values, while the "Red/Blue ratio" just flows along with the water.

5. The Secret Weapon: A New Inequality

The hardest part of the proof was comparing the "squishy" state to the "stiff" state. The author had to compare two different ways of measuring the difference between the fluids:

The L1 Norm: A measure of the total "amount" of difference (like counting the total number of mismatched pixels).
The L2 Norm: A measure of the "intensity" of the difference (like measuring how bright the mismatched pixels are).

Usually, these are hard to compare. The author developed a new mathematical inequality (a variant of the Csiszár-Kullback-Pinsker inequality) that acts like a translator. It allows them to say: "If the total amount of difference is small, then the intensity of the difference must also be small."

This translator allowed the author to use the "Scorecard" (Entropy) to prove that the fluids must converge to the target state.

6. The Conclusion: The Smooth Transition

The paper concludes that:

As the fluid slows down, the Red and Blue densities stop fluctuating and become constant.
The velocity of the fluid becomes perfectly smooth and "incompressible" (it doesn't expand or contract).
The mixing ratio (how much Red vs. Blue is in each spot) is simply transported by the flow, just like a dye spreading in water.

In simple terms: The author proved that if you take a chaotic, compressible mix of two fluids and slow it down, it doesn't explode or break. Instead, it gracefully transforms into a smooth, incompressible flow where the two fluids move together, maintaining their distinct identities but flowing as one perfect, non-squishy unit.

This is a big deal for engineers and physicists because it gives them a solid mathematical foundation to model things like fuel sprays in engines or blood flow in veins, knowing that their simplified "slow-flow" models are actually correct limits of the complex "fast-flow" reality.

1. Problem Statement

The paper addresses the rigorous mathematical justification of the low Mach number limit for a compressible isentropic two-phase flow system with an algebraic closure.

The System: The model describes two fluid phases (denoted $+$ and $-$) sharing a single common velocity field ( $u$ ) and an equal pressure ( $p_+ = p_-$ ). The phases interact via barotropic pressure laws $p_i(\rho_i) = \rho_i^{\gamma_i}$ with $\gamma_i \geq 2$ . The system is governed by conservation of mass for each phase and a single momentum equation for the mixture.
The Closure: The pressure equality $p_+(\rho_+) = p_-(\rho_-)$ creates a non-trivial algebraic coupling between the partial densities and volume fractions ( $\alpha_\pm$ ). This distinguishes the system from models with differential closures or distinct velocities.
The Limit: The goal is to analyze the behavior of the system as the Mach number $\varepsilon \to 0$ $ε \to 0$ . Formally, this limit should yield an incompressible non-homogeneous Navier-Stokes system where:
- The partial densities $\rho_\pm$ become constant in space and time.
- The volume fractions $\alpha_\pm$ evolve via a transport equation driven by the divergence-free velocity field.
- The mixture density $\rho = \alpha_+\rho_+ + \alpha_-\rho_-$ remains variable and is transported by the flow.

2. Methodology

The author employs the Relative Entropy Method, a powerful tool for proving the convergence of weak solutions of compressible systems to their incompressible limits. The methodology involves several novel technical steps:

A. Construction of the Relative Entropy Functional

The core of the analysis is a tailored relative entropy functional $E = E_1 + E_2$ , measuring the distance between a weak solution of the compressible system $(\rho^\varepsilon_\pm, \alpha^\varepsilon_\pm, u^\varepsilon)$ and a strong solution of the incompressible limit system $(\rho_\pm, \alpha_\pm, u)$ .

$E_1$ (Standard Kinetic and Internal Energy):
- Contains the kinetic energy difference and the standard relative internal energy (based on the convex function $H(\rho|r)$ ) scaled by $1/\varepsilon^2$ .
- Controls the velocity difference and the deviation of densities from their constant limit states.
$E_2$ (Logarithmic Comparison for Partial Masses):
- Novelty: A new term introduced specifically for the two-phase structure:
  $E_2 = \int_\Omega R^\varepsilon_+ \ln_+ \left(\frac{R^\varepsilon_+}{R_+}\right) + \int_\Omega R^\varepsilon_- \ln_+ \left(\frac{R^\varepsilon_-}{R_-}\right)$
  where $R^\varepsilon_\pm = \alpha^\varepsilon_\pm \rho^\varepsilon_\pm$ are the partial mass densities.
- Purpose: This term is designed to control the $L^1$ norm of the difference in partial masses, which is crucial because the algebraic closure links the volume fractions to the densities in a way that standard entropy estimates cannot fully capture.

B. Key Estimates and Structural Exploitation

To close the Gronwall argument required for convergence, the author derives specific inequalities exploiting the algebraic structure of the system:

$L^1$ vs. $L^2$ Comparison for Partial Densities:
- The author proves a lemma showing that the difference in volume fractions $|\alpha^\varepsilon_\pm - \alpha_\pm|$ can be bounded by the differences in partial masses $|R^\varepsilon_\pm - R_\pm|$ .
- This relies on the implicit function theorem applied to the pressure equality constraint $p_+(\rho_+) = p_-(\rho_-)$ .
- This allows the control of the $L^2$ norm of $R^\varepsilon_\pm - R_\pm$ using the entropy $E_1$ and the $L^1$ norm of the difference.
Refined Csiszár-Kullback-Pinsker (CKP) Inequality:
- A variant of the CKP inequality is established to bound the $L^1$ norm $\|f - g\|_{L^1}$ directly by the integral $\int f \ln_+(f/g)$ .
- Significance: Standard CKP inequalities bound the square of the $L^1$ norm. The author argues that a square bound is insufficient for the specific differential inequalities arising in this two-phase context; a direct $L^1$ bound is necessary to apply Gronwall's lemma effectively.

C. Convergence Analysis

The proof proceeds by deriving a differential inequality for the total relative entropy $E(t)$ .

The remainder terms in the entropy inequality are estimated using the structural bounds derived above.
Under the assumption of well-prepared initial data (where the initial energy deviation from the limit is $O(\varepsilon^2)$ ), the Gronwall lemma implies that $E(t) \to 0$ as $\varepsilon \to 0$ .

3. Key Contributions

Rigorous Justification of the Limit: The paper provides the first rigorous proof of the low Mach number limit for a two-phase system with algebraic closure and single velocity, a regime that had previously only been treated formally or in simplified settings.
Novel Relative Entropy Functional: The introduction of the $E_2$ functional involving logarithmic terms of partial masses is a significant methodological innovation. It bridges the gap between the compressible partial densities and the incompressible volume fractions.
Structural Analysis of Algebraic Closures: The paper demonstrates how to handle the non-trivial coupling induced by the pressure equality constraint. The derived estimates relating $L^1$ and $L^2$ norms of partial masses are specific to the bi-fluid algebraic structure.
Refined Inequalities: The adaptation of the Csiszár-Kullback-Pinsker inequality to provide direct $L^1$ control (rather than $L^1$ squared) is a technical contribution that enables the closure of the energy estimates.

4. Main Results

Under the assumptions of well-prepared initial data and sufficient regularity for the limit solution, the main theorem (Theorem 1.6) establishes the following convergence as $\varepsilon \to 0$ :

Partial Masses: $R^\varepsilon_\pm \to R_\pm$ strongly in $L^\infty(0, T; L^2(\Omega))$ .
Partial Densities: $\rho^\varepsilon_\pm \to \rho_\pm$ strongly in $L^\infty(0, T; L^2(\Omega))$ , where $\rho_\pm$ are constant states.
Volume Fractions: $\alpha^\varepsilon_\pm \to \alpha_\pm$ weakly- $*$ in $L^\infty(0, T; L^2(\Omega))$ . The limit volume fractions satisfy the transport equation $\partial_t \alpha_\pm + u \cdot \nabla \alpha_\pm = 0$ .
Velocity: $u^\varepsilon \to u$ strongly in $L^2(0, T; L^2(\Omega))$ .
Incompressibility: $\text{div}(u^\varepsilon) \rightharpoonup 0$ weakly in $L^2$ , confirming the limit velocity is divergence-free.

The limit system is identified as the non-homogeneous incompressible Navier-Stokes equations with a variable density $\rho = \alpha_+ \rho_+ + \alpha_- \rho_-$ determined by the transported volume fractions.

5. Significance

Mathematical Physics: This work fills a critical gap in the theory of compressible multiphase flows. While single-phase low Mach limits are well-understood, the multi-phase case with algebraic closures presents unique analytical challenges due to the coupling of phases.
Engineering Applications: Two-phase flows with equal pressure and single velocity are common in industrial applications (e.g., bubbly flows, spray dynamics) where mechanical relaxation is fast. This paper validates the use of incompressible models for such regimes, providing a theoretical foundation for numerical simulations that switch between compressible and incompressible formulations.
Analytical Framework: The techniques developed, particularly the specific relative entropy functional and the $L^1/L^2$ comparison lemmas, offer a new toolkit for analyzing other singular limits in complex coupled fluid systems.