Low Mach Number Limit and Convergence Rates for a Compressible Two-Fluid Model with Algebraic Pressure Closure

This paper rigorously establishes the low Mach number limit and derives explicit convergence rates for a three-dimensional viscous compressible two-fluid model with algebraic pressure closure, proving that its well-prepared strong solutions converge to the incompressible Navier–Stokes equations as the Mach number tends to zero.

The Gist

Imagine you are watching a busy highway. Sometimes, the cars are packed so tightly together that they can't move freely; they bump into each other, compress, and expand like a springy rubber band. This is a compressible fluid (like air or gas). Other times, the cars are spaced out, moving smoothly without ever really squishing together. This is an incompressible fluid (like water).

In physics, the "Mach number" is like a speedometer that tells us how "squishy" the fluid is behaving. A high Mach number means the fluid is acting like a springy gas (compressible). A low Mach number means it's acting like a smooth, incompressible liquid.

The Big Problem: The "Black Box" Pressure

For decades, scientists have known how to predict what happens when a gas slows down and starts acting like a liquid. They have a simple rulebook (equations) for this.

However, this paper tackles a much trickier scenario: Two different fluids mixed together (like oil and water, or gas and liquid) that are moving at the same speed but have different densities.

The real headache here is the pressure. In simple models, pressure is like a clear instruction: "If you squeeze me this much, I push back that much." But in this specific two-fluid model, the pressure is a black box. You can't just look at the ingredients and calculate the pressure directly. Instead, the pressure is determined by a hidden, complex relationship between the two fluids. It's like trying to guess the temperature of a soup by tasting it, but the recipe is written in a code you can't see.

Because this "code" is hidden (mathematically called an implicit algebraic closure), proving what happens when the fluids slow down (the low Mach number limit) is incredibly difficult. It's like trying to solve a puzzle where the pieces keep changing shape as you try to fit them together.

The Solution: A Mathematical Tightrope Walk

The authors of this paper (Yang Li, Mária Lukáčová-Medviďová, and Ewelina Zatorska) managed to walk this tightrope. They proved that even with this messy, hidden pressure rule, if you start with the fluids in a "well-prepared" state (meaning they aren't chaotic at the start), the system behaves beautifully as it slows down.

Here is what they did, using some analogies:

1. The "Uniform" Safety Net:
   Usually, when you try to slow down a complex system, the math might blow up or become unstable. The authors built a "safety net" (mathematically called uniform estimates). They proved that no matter how small the "squishiness" gets (as the Mach number approaches zero), the solution stays stable and doesn't explode. They showed that the two fluids will eventually settle down to a constant, calm state, and the velocity will match the smooth flow of water (the incompressible Navier-Stokes equations).

2. The "Relative Energy" Scale:
   To measure how fast the system converges to the smooth liquid state, they used a tool called a relative energy argument. Imagine you have a heavy, wobbly cart (the compressible gas) and a smooth, sleek sled (the ideal liquid). The authors didn't just say, "They look similar." They put the cart and the sled on a giant, sensitive scale. They measured the exact "energy difference" between the two.

   They found that as the "squishiness" (Mach number) gets smaller, the cart doesn't just eventually look like the sled; it transforms into the sled at a predictable speed. They calculated exactly how much error remains at every step.

The Results: Speed and Precision

The paper delivers two major victories:

• Existence: They proved that a solution actually exists for a specific amount of time, regardless of how small the Mach number is. It's not just a theoretical guess; the math holds up.
• Convergence Rates: This is the "cool factor." They didn't just say "it converges." They gave a receipt. They proved that the error in the density is proportional to the square of the Mach number ($\epsilon^2$), and the error in the speed is proportional to the Mach number ($\epsilon$).
  • Analogy: If you reduce the "squishiness" by half, the error in the density drops by a factor of four, and the error in the speed drops by half. This gives engineers and scientists a precise way to predict how good their approximations will be.

Why Does This Matter?

This isn't just abstract math. This kind of model is used in real-world engineering, like:

• Oil and Gas Pipelines: Where liquid oil and gas bubbles flow together.
• Rocket Engines: Where fuel and oxidizer mix.
• Medical Devices: Like blood flow with gas bubbles.

Previously, scientists had to make huge simplifications to study these systems. This paper removes the need for those simplifications, proving that even with the most complex, "hidden" pressure rules, nature still follows a predictable path toward smooth, liquid-like flow. They turned a chaotic, black-box problem into a clear, quantifiable journey.

Technical Summary

Here is a detailed technical summary of the paper "Low Mach Number Limit and Convergence Rates for a Compressible Two-Fluid Model with Algebraic Pressure Closure" by Yang Li, Mária Lukáčová-Medviďová, and Ewelina Zatorska.

1. Problem Statement

The paper investigates the low Mach number limit for a viscous compressible two-fluid model in a three-dimensional periodic domain ($\mathbb{T}^3$). The system describes two immiscible, isentropic fluids sharing a common velocity field $u$ and coupled through an algebraic pressure closure.

The Model:
The governing equations (in conservative variables $R = \alpha_+ \rho_+$ and $Q = \alpha_- \rho_-$) are:
$$\begin{cases} \partial_t R + \text{div}(Ru) = 0, \\ \partial_t Q + \text{div}(Qu) = 0, \\ \partial_t [(R+Q)u] + \text{div}[(R+Q)u \otimes u] + \nabla p = \mu \Delta u + (\mu+\lambda)\nabla \text{div} u, \end{cases}$$
subject to the algebraic pressure closure:
$$p = (\rho_+)^{\gamma_+} = (\rho_-)^{\gamma_-}, \quad \gamma_\pm > 1.$$
Crucially, the pressure $p$ is not an explicit function of $R$ and $Q$. Instead, it is defined implicitly via the relation $Q = (1 - R/Z)Z^\gamma$, where $Z = \rho_+$ is determined implicitly as a function $Z(R, Q)$.

The Challenge:
The primary mathematical difficulty lies in the implicit structure of the pressure law. Unlike models with explicit pressure laws (e.g., $P(R,Q) = R^\gamma + Q^\beta$), the derivatives of the pressure with respect to the densities involve the derivatives of the implicitly defined map $Z(R,Q)$. This makes the derivation of uniform estimates and the analysis of the singular limit (as the Mach number $\epsilon \to 0$) significantly more delicate.

The goal is to prove that, for well-prepared initial data, the solutions to the rescaled compressible system converge to the solution of the incompressible Navier-Stokes equations as $\epsilon \to 0$, and to establish explicit convergence rates.

2. Methodology

The authors employ a combination of uniform high-order energy estimates and a relative energy argument, adapted specifically to handle the implicit pressure structure.

A. Non-dimensional Scaling

The system is rescaled using the Mach number $\epsilon$:
$$R \to R_\epsilon, \quad Q \to Q_\epsilon, \quad u \to \epsilon u_\epsilon, \quad \mu \to \epsilon \mu_\epsilon.$$
The rescaled momentum equation contains a singular pressure term $\nabla p(Z_\epsilon)/\epsilon^2$. The expected limit system is the incompressible Navier-Stokes equations:
$$\partial_t u + u \cdot \nabla u + \nabla \Pi = \frac{\mu}{2} \Delta u, \quad \text{div } u = 0.$$

B. Existence and Uniform Regularity (Theorem 2.1)

To prove existence on a time interval independent of $\epsilon$, the authors use a fixed-point argument (contraction mapping principle) on a suitable set of functions $A_{T^*, \epsilon}$.

1. Linearization: They define a map $\chi$ where the nonlinear terms (convection and pressure gradients) are evaluated using a "frozen" background state $(r_\epsilon, q_\epsilon, v_\epsilon)$.
2. Symmetrization: They construct a symmetrized energy functional $F_s$ based on the structure of the linearized system. This functional accounts for the implicit derivatives $\partial_R Z$ and $\partial_Q Z$.
3. Uniform Estimates: They derive high-order Sobolev estimates ($H^s$ with $s \ge 4$) for the linearized system. A key technical step involves carefully estimating commutators arising from the implicit pressure term $Z(R,Q)$. They utilize Sobolev embedding and specific bounds on the derivatives of the implicit function $Z$ to close the estimates.
4. Contraction: They show that for a sufficiently small time $T^*$ and small $\epsilon$, the map $\chi$ is a contraction, ensuring the existence of a unique classical solution $(R_\epsilon, Q_\epsilon, u_\epsilon)$ with uniform bounds independent of $\epsilon$.

C. Convergence Rates (Theorem 2.3)

To obtain quantitative rates, the authors utilize a relative energy inequality.

1. Relative Energy: They define a relative energy functional measuring the distance between the compressible solution $(R_\epsilon, Q_\epsilon, u_\epsilon)$ and the incompressible limit $(1, 1, u)$.
2. Target Test Function: They test the momentum equation of the compressible system against the incompressible velocity $u$.
3. Remainder Estimates: The proof involves estimating the remainder terms generated by the difference between the compressible and incompressible systems. The implicit pressure structure requires careful handling of the terms involving $\nabla R_\epsilon$ and $\nabla Q_\epsilon$ scaled by $1/\epsilon^2$.
4. Gronwall's Inequality: By combining the energy inequality with the uniform bounds, they derive a differential inequality that yields the convergence rates via Gronwall's lemma.

3. Key Contributions and Results

Main Theorems

1. Uniform Existence (Theorem 2.1): For well-prepared initial data (where initial perturbations are of order $\epsilon^2$ for densities and $\epsilon$ for velocity), the system admits a unique classical solution on a time interval $[0, T^*]$ independent of $\epsilon$.

   • Convergence: As $\epsilon \to 0$, $R_\epsilon \to 1$ and $Q_\epsilon \to 1$ strongly in $L^\infty(0, T^*; H^s)$, and $u_\epsilon \to u$ strongly in $C([0, T^*]; H^{s-\delta'})$, where $u$ solves the incompressible Navier-Stokes equations.

2. Convergence Rates (Theorem 2.3): Under an additional smallness assumption on the initial energy (specifically related to the acoustic energy), the authors establish explicit rates:

   • Densities: $\|R_\epsilon - 1\|_{H^s}^2 + \|Q_\epsilon - 1\|_{H^s}^2 \le C \epsilon^2$.
   • Velocity: $\|u_\epsilon - u\|_{L^2}^2 + \int_0^t \|u_\epsilon - u\|_{H^1}^2 d\tau \le C \epsilon$.
   • Divergence: $\|\text{div } u_\epsilon\|_{H^{s-1}} \le C \epsilon$.

Technical Innovations

• Implicit Pressure Handling: The paper successfully navigates the non-explicit nature of the pressure law. The authors derive precise bounds for the derivatives of the implicit map $Z(R,Q)$ and incorporate them into the energy estimates, avoiding the need for explicit algebraic inversion which is impossible here.
• Relaxed Assumptions: Unlike previous works on liquid-gas models (e.g., Yao et al. [24]) which required a uniform bound on the ratio of initial densities ($Q_0/R_0 \le 1$), this result does not require such a bound. The convergence rates hold under the general assumption $\gamma_\pm > 1$.
• Optimality of Rates: The derived rates ($O(\epsilon)$ for velocity, $O(\epsilon^2)$ for densities) are consistent with standard low Mach number limits for compressible fluids but are the first to be established for this specific algebraic closure model in the strong solution setting.

4. Significance

• Rigorous Justification: This work provides the first rigorous justification of the low Mach number limit for compressible two-fluid models with algebraic pressure closure in the framework of strong solutions.
• Quantitative Analysis: By establishing explicit convergence rates, the paper bridges the gap between theoretical asymptotic analysis and practical numerical simulations, offering error bounds for approximating compressible two-phase flows with incompressible models.
• Broader Applicability: The methodology developed for handling implicit pressure laws can be extended to other complex multiphase flow models where pressure is defined through complex equations of state or algebraic constraints, not just simple explicit laws.
• Comparison to Literature: The results improve upon existing literature by removing restrictive assumptions on initial density ratios and providing explicit rates where previous works (like Lebot [11]) only identified the limit system without rates or required stronger constraints on adiabatic constants ($\gamma \ge 2$).

In summary, the paper resolves a significant open problem in mathematical fluid dynamics by proving that the complex, implicitly coupled two-fluid system behaves like the standard incompressible Navier-Stokes equations in the low Mach limit, with precise quantitative error estimates.