🔢 mathematics

Global regularity of the multi-dimensional compressible Navier-Stokes-Korteweg system with large initial data

This paper establishes the global existence of strong solutions to the multi-dimensional compressible Navier-Stokes-Korteweg system with arbitrarily large initial data on the torus by introducing a novel modified Nash-Moser iteration that links the effective velocity to the density's lower bound.

Original authors: Xiangdi Huang, Weili Meng, Xueyao Zhang

Published 2026-02-09

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Xiangdi Huang, Weili Meng, Xueyao Zhang

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.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 a fluid, like a gas or a liquid, that isn't just flowing smoothly but is also trying to pull itself together or push itself apart based on how tightly packed its molecules are. This is the world of the Navier-Stokes-Korteweg system. Think of it as a complex dance where the fluid has three main personalities:

Density: How crowded the dancers are.
Velocity: How fast they are moving.
Capillarity: A special "surface tension" force that acts like a rubber band, trying to smooth out sharp changes in how crowded the dancers are.

For a long time, mathematicians could predict how this dance would go if the dancers started off moving slowly and gently (small initial data). But what if the dancers started off sprinting, jumping, and spinning wildly (large initial data)? Could the system hold together forever, or would it eventually collapse into chaos (a mathematical "blow-up")?

This paper, by Huang, Meng, and Zhang, says: Yes, it can hold together. They proved that even if the fluid starts off in a chaotic, high-energy state, it will continue to exist and behave predictably for all time, provided the fluid follows certain physical rules (specifically regarding how pressure relates to density).

Here is how they solved this puzzle, broken down into simple steps:

1. The Problem: The "Crunch" and the "Stretch"

In fluid dynamics, there are two main fears:

The Crunch: The density gets so high it becomes infinite (like a black hole forming in the fluid).
The Stretch: The density gets so low it hits zero (creating a vacuum), which breaks the math equations.

The authors needed to prove that the density stays in a "Goldilocks zone"—not too high, not too low—no matter how wild the start was.

2. The Secret Weapon: The "Effective Velocity"

The authors introduced a clever trick. Instead of looking at the fluid's speed ( $u$ ) and its density ( $\rho$ ) separately, they combined them into a new character called the "Effective Velocity" ( $v$ ).

Analogy: Imagine you are walking on a moving walkway at an airport. Your total speed is your walking speed plus the speed of the walkway. The authors realized that if they looked at the "walkway speed" (the density gradient) combined with your walking speed, the math became much easier to handle. This new character, $v$ , acts like a stabilizer for the whole system.

3. Step One: Keeping the Crowd from Exploding (Upper Bound)

First, they had to prove the density wouldn't get infinitely high.

The Method: They used a technique called De Giorgi iteration.
The Analogy: Imagine a crowd of people in a room. If the room gets too crowded, the pressure builds up. The authors showed that even if the crowd starts huge, the "pressure" (density) has a natural limit. They used the "Effective Velocity" to show that the crowd can't squeeze into a space smaller than a certain size. It's like proving that no matter how hard you push a spring, it has a maximum compression point before it stops.

4. Step Two: Keeping the Crowd from Vanishing (Lower Bound)

This was the hardest part. They had to prove the density wouldn't drop to zero (vacuum).

The Old Problem: Previous attempts to solve this relied on a mathematical estimate that worked for small crowds but failed when the crowd was huge. It was like trying to measure a mountain with a ruler meant for a hill; the ruler broke.
The New Solution: The authors invented a modified Nash-Moser iteration.
The Analogy: Think of this as a "feedback loop" or a "self-correcting thermostat." They discovered a special relationship between the "Effective Velocity" and the "Inverse Density" (how empty the room is).
- They proved that if the room starts to get too empty, the "Effective Velocity" (the stabilizer) automatically kicks in to push the density back up.
- They showed that the speed of this stabilizer grows only as fast as the square root of the logarithm of how empty the room is.
- Why this matters: Logarithms grow very slowly. This means the stabilizer is incredibly efficient. Even if the room gets 1,000 times emptier, the stabilizer doesn't need to go 1,000 times faster; it only needs to go a tiny bit faster. This prevents the system from ever hitting zero.

5. The Result: A Forever Dance

By proving the density stays within a safe range (not too high, not too low), they could then prove that the fluid's movement remains smooth and predictable forever.

The Claim: This is the first time anyone has proven that this specific type of fluid (compressible, with capillarity) can handle arbitrarily large starting chaos in three dimensions without breaking down.

Summary

The paper is a mathematical victory lap. The authors took a fluid system that was known to be unstable with wild starting conditions and proved it is actually robust. They did this by:

Combining speed and density into a single "Effective Velocity" character.
Using a "thermostat" logic to prove the fluid can never become a vacuum.
Using a "pressure valve" logic to prove the fluid can never become infinitely dense.

They didn't just say "it works"; they built a mathematical bridge that connects the chaotic start to a stable, eternal future for the fluid.

Technical Summary: Global Regularity of the Multi-dimensional Compressible Navier-Stokes-Korteweg System with Large Initial Data

Problem Statement
This paper addresses the global existence of strong solutions for the multi-dimensional compressible Navier–Stokes–Korteweg (NSK) system with density-dependent viscosity and capillarity coefficients. The system, derived by Dunn and Serrin, models capillarity in compressible fluids and reduces to the quantum Navier–Stokes system under specific coefficient choices ( $\mu(\varrho) = \nu_0 \varrho$ , $\lambda(\varrho) = 0$ , $\kappa(\varrho) = \kappa_0^2/\varrho$ with $\nu_0 = \kappa_0$ ).

The primary mathematical challenge lies in establishing global regularity for arbitrarily large initial data in dimensions $N=2$ and $N=3$ , specifically when the fluid is away from vacuum. While global weak solutions are well-studied, the existence of global strong solutions for large data in higher dimensions remains an open problem. The authors focus on the torus $T^N$ with periodic initial data $(\varrho_0, u_0)$ satisfying $0 < \underline{\varrho} \le \varrho_0 \le \overline{\varrho}$ .

Methodology
The proof strategy relies on transforming the original system into a parabolic system using the effective velocity $v = u + \nabla \log \varrho$ . This transformation, introduced by Jüngel, converts the NSK system into:
$\begin{cases} \varrho_t + \text{div}(\varrho v) - \Delta \varrho = 0, \\ \varrho v_t + \varrho(v - \nabla \log \varrho) \cdot \nabla v + \nabla P = \text{div}(\varrho \nabla v). \end{cases}$
The analysis proceeds in three main steps:

Upper Bound for Density:
The authors first establish an upper bound for the density $\varrho$ . This is achieved by deriving higher integrability estimates for the effective velocity $v$ (specifically $\|\varrho^{1/(q+2)}v\|_{L^{q+2}}$ ) and applying the De Giorgi iteration technique. This step adapts methods previously used for 2D problems to the 3D case, requiring specific constraints on the adiabatic exponent $\gamma$ ( $\gamma < 7/3$ for $N=3$ ).
Positive Lower Bound for Density:
The most critical and novel contribution is obtaining a positive lower bound for the density (preventing vacuum formation). Previous approaches (e.g., Haspot) relied on estimates of $\|\varrho^{1/(p+2)}v\|_{L^{p+2}}$ that were independent of $p$ to pass to the limit $p \to \infty$ . The authors demonstrate that such uniform estimates are unattainable. Instead, they develop an original modified Nash-Moser type iteration.
- They establish a critical logarithmic control relation between the effective velocity and the inverse density:
  $\|v\|_{L^\infty(T^N \times (0,T))} \le C \sqrt{\log(e + \|\varrho^{-1}\|_{L^\infty(T^N \times (0,T)))}}.$
- This estimate allows them to close the bootstrap argument for the lower bound of the density without requiring the uniform $L^\infty$ bound on the effective velocity that previous methods failed to achieve.
Higher-Order Regularity:
With both upper and lower bounds for the density established, the authors derive higher-order energy estimates for $\varrho$ and $v$ . These estimates ensure that the local strong solution can be extended globally in time.

Key Results
The main theorem (Theorem 1.1) establishes the global existence and uniqueness of strong solutions $(\varrho, u)$ for the system on $T^N$ ( $N=2, 3$ ) with large initial data, provided:

For $N=2$ : $\gamma \in [1, \infty)$ .
For $N=3$ : $\gamma \in [1, 7/3)$ .

The solution satisfies:

Uniform bounds: $(C(T))^{-1} \le \varrho(x,t) \le C(T)$ .
Regularity: $\varrho \in C([0,T]; H^3) \cap L^2(0,T; H^4)$ and $u \in C([0,T]; H^2) \cap L^2(0,T; H^3)$ .

Significance and Claims
The paper claims to present the first existence result of global strong solutions for compressible fluid dynamics equations with physical significance in general three-dimensional domains with arbitrarily large initial data.

Contrast with Incompressible Fluids: The authors note that for incompressible fluids, the only known result for global strong solutions in 3D with large data dates back to Cao and Titi (2007).
Overcoming Previous Limitations: The work resolves the difficulty encountered in previous studies (such as Haspot [25]) where the dependence of estimates on the parameter $p$ prevented the derivation of a positive lower bound for the density in the limit $p \to \infty$ .
Scope: The result applies to both isothermal ( $\gamma=1$ ) and non-isothermal cases. The restriction on $\gamma$ in 3D is attributed to the technical requirements of estimating the upper bound of the density via De Giorgi iteration.

The authors maintain a modest tone regarding the scope, noting that the restriction on $\gamma$ in 3D is a consequence of the current method for bounding the density, and the 2D case allows for any $\gamma \ge 1$ due to stronger energy estimates.