🔢 mathematics

Local well-posedness of strong solutions to the compressible Navier-Stokes equations with degenerate viscosities and far field vacuum in 3D exterior domains

This paper establishes the local well-posedness of strong solutions to the isentropic compressible Navier-Stokes equations with degenerate viscosities and far-field vacuum in 3D exterior domains under Navier-slip boundary conditions, utilizing a method that simultaneously addresses boundary and vacuum challenges while decoupling the viscosity exponent from the adiabatic coefficient.

Original authors: Jiaxu Li, Boqiang Lü, Bing Yuan

Published 2026-02-05

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Jiaxu Li, Boqiang Lü, Bing Yuan

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 vast, open ocean (the "exterior domain") where a fluid is flowing. This fluid is special: it's a gas that can get incredibly thin, eventually turning into a complete vacuum (empty space) as you move further out to the horizon. Furthermore, the fluid's "thickness" or resistance to flow (viscosity) isn't constant; it changes depending on how dense the gas is. When the gas gets very thin, it becomes almost frictionless, which makes the math describing its movement extremely difficult to solve.

This paper is like a team of mathematicians building a new, more robust bridge to cross a very tricky river. Here is the breakdown of their journey:

1. The Problem: A Slippery Slope with Empty Space

The scientists are trying to predict how this gas moves using the Navier-Stokes equations (the standard rulebook for fluid motion). However, they are facing two major hurdles:

The "Far-Field Vacuum": As you go far away from the center, the gas density drops to zero. In math terms, this is like a wall that disappears. When density hits zero, the equations usually break down because you can't divide by zero or multiply by nothing.
The "Slippery" Boundary: The gas is flowing near a solid object (like a rock in the river). The paper uses a "Navier-slip" condition, which means the gas doesn't stick perfectly to the rock; it slides a little bit. This sliding adds a layer of complexity, especially when the gas is getting thin near the edges.

Previous attempts to solve this were like trying to walk on a tightrope that only existed if you followed a very specific, narrow path (a strict relationship between the gas's heat properties and its thinning rate). If the gas didn't fit that narrow path, the math failed.

2. The Solution: A New Toolkit

The authors (Li, Lü, and Yuan) developed a new mathematical method that acts like a universal key. Instead of forcing the gas to fit a narrow path, they reorganized the equations to handle the "slippery" nature of the viscosity directly.

The "Magic" Reformulation: They took the messy equations and divided them by the density in a clever way. This transformed the problem into a form where the "vacuum" (zero density) doesn't cause the math to explode.
The "Anchor" Trick: One of the biggest challenges was estimating the speed of the gas near the boundary. Because the gas gets thin, the math usually loses its grip. The authors discovered a new way to prove that the gas density stays "thick enough" right next to the boundary, acting like an anchor that keeps the calculations stable. This allowed them to handle the sliding (slip) condition without needing the gas to follow those old, restrictive rules.

3. The Result: Local Well-Posedness

The paper proves "Local Well-Posedness." In plain English, this means:

Existence: If you start with a specific, realistic setup (a certain amount of gas moving at a certain speed), a solution does exist.
Uniqueness: There is only one correct way that gas will behave for a short period of time. You won't get two different answers for the same starting point.
Regularity: The solution is "strong," meaning it's smooth and well-behaved, not jagged or chaotic, for a finite amount of time.

Think of it as successfully predicting the weather for the next hour with high confidence, even though the atmosphere is getting thinner and thinner at the edges of the map.

4. The "Blow-Up" Warning

The paper also includes a "safety alarm" (Theorem 1.2). It says: "We can predict the flow for a while, but if the flow starts to get too wild (specifically, if the twisting of the fluid or the thinning of the gas near the boundary gets too extreme), the prediction will eventually fail." This tells us exactly what to watch out for if the solution stops working.

Summary Analogy

Imagine trying to predict how a crowd of people moves through a giant, open plaza that fades into a foggy void at the edges.

Old Method: You could only predict the crowd if everyone moved at a very specific speed relative to how crowded the area was. If the crowd got too thin, your prediction broke.
This Paper's Method: The authors invented a new way to track the crowd that works even when the edges are foggy and empty. They proved that for a short time, you can predict exactly where everyone will be, regardless of how the crowd thins out, as long as you know where they started. They also identified the exact moment the crowd might get so chaotic that your prediction would fail.

In short: This paper removes a major mathematical roadblock, allowing scientists to model compressible gases in open spaces with changing viscosity much more freely and accurately than before, without needing to force the gas into a narrow, artificial box.

Technical Summary: Local Well-Posedness of Strong Solutions to Compressible Navier-Stokes Equations with Degenerate Viscosities and Far-Field Vacuum in 3D Exterior Domains

Problem Statement
This paper investigates the local well-posedness of strong solutions to the three-dimensional isentropic compressible Navier-Stokes equations in a general exterior domain $\Omega = \mathbb{R}^3 \setminus \bar{D}$ , where $D$ is a bounded smooth domain. The system is characterized by density-dependent (degenerate) viscosities, specifically $\mu(\rho) = \mu\rho^\delta$ and $\lambda(\rho) = \lambda\rho^\delta$ with $0 < \delta < 1$ . The study addresses the initial-boundary value problem (IBVP) where the initial density approaches a far-field vacuum ( $\rho \to 0$ as $|x| \to \infty$ ) and the fluid satisfies Navier-slip boundary conditions on $\partial\Omega$ :
$u \cdot n = 0, \quad \text{curl} u \times n = -A(x)u \quad \text{on } \partial\Omega,$
where $A$ is a symmetric matrix. The primary challenge lies in the double degeneracy of the momentum equation when vacuum appears and the complexities introduced by the boundary terms in an exterior domain.

Methodology
The authors adopt a reformulation strategy motivated by Li-Li [17], dividing the momentum equation by $\rho^\delta$ to handle the degeneracy. The system is transformed into:
$\rho^{1-\delta}u_t + \rho^{1-\delta}u \cdot \nabla u - Lu + \frac{a\gamma}{\gamma-\delta}\nabla \rho^{\gamma-\delta} = \delta \nabla \log \rho \cdot S(u),$
where $L$ is the Lamé operator.

To establish the local existence of strong solutions, the authors employ the following technical approaches:

Approximation in Annular Domains: The problem is first solved in a bounded annular domain $\Omega_R = \Omega \cap B_R$ to approximate the exterior domain. A priori estimates are derived that are independent of the radius $R$ .
Handling Boundary Terms: A significant portion of the analysis is dedicated to managing the boundary terms arising from the Navier-slip condition. Unlike the Cauchy problem, the boundary term in the energy estimates involves the matrix $A$ , which is not necessarily positive definite. The authors utilize a new observation that the density $\rho$ retains a positive lower bound near the boundary $\partial\Omega$ (specifically in $\Omega_0 = \Omega \cap B_{2R_0}$ ), allowing for the control of negative powers of $\rho$ in boundary integrals.
Weighted Estimates: The authors derive crucial weighted estimates for terms such as $\nabla \rho^{\frac{\delta-1}{2}}$ and $\rho^{\frac{\delta-1}{2}} \text{div} u$ . They overcome the difficulty of deriving elliptic estimates in annular domains (where the first Betti number is non-zero) by utilizing a cut-off function method and a specific div-curl lemma adapted for the exterior domain structure.
Blow-up Criterion: To address the extension of local solutions to global time, a blow-up criterion is established. This involves controlling the $L^1(0,T; L^\infty)$ norm of the deformation tensor $D(u)$ and specific norms of the density gradient near the boundary.

Key Contributions and Results
The paper establishes the following main results:

Theorem 1.1 (Local Well-Posedness): Under the conditions $\gamma > 1$ and $0 < \delta < 1$ , if the initial data $(\rho_0, u_0)$ satisfies specific regularity and compatibility conditions (including $\rho_0^{\frac{\gamma-1+\delta}{2}} \in D^1_0 \cap D^2$ and a compatibility condition $L(u_0) = \rho_0^{\frac{1-\delta}{2}}g$ ), there exists a unique strong solution on a time interval $[0, T_0]$ . The solution satisfies a specific set of regularity properties, including $\rho^{\frac{\delta-1}{2}} \in L^\infty(0, T_0; L^6(\Omega_0))$ , which ensures the density remains bounded away from zero near the boundary.
Independence of Parameters: A critical contribution is that the selection of the viscosity exponent $\delta$ is independent of the adiabatic exponent $\gamma$ . This contrasts with previous works (e.g., Xin-Zhu [33], Liu-Zhong [23]) which required a restrictive relationship between $\delta$ and $\gamma$ (specifically $\delta \to 1$ as $\gamma \to 1$ ).
Relaxed Compatibility Conditions: The paper requires only one compatibility condition (1.13), removing two additional restrictive constraints found in earlier literature for exterior domains.
Theorem 1.2 (Blow-up Criterion): A criterion for the finite-time blow-up of strong solutions is provided. If the maximal existence time $T^*$ is finite, then the sum of the $L^1(0,T; L^\infty)$ norm of $D(u)$ , the $L^\infty(0,T; L^6(\Omega_0))$ norm of $\rho^{\frac{\delta-1}{2}}$ , and the $L^\infty(0,T; D^{1,2})$ norm of $\nabla \rho^{\frac{\delta-1}{2}}$ must diverge.

Significance
The paper claims to extend the theory of compressible Navier-Stokes equations with degenerate viscosities from the Cauchy problem to more general boundary value problems in exterior domains. By decoupling the constraints on $\delta$ and $\gamma$ , the results encompass a wider range of physical scenarios, including Maxwellian molecules, without requiring the specific limiting relations previously necessary. The development of new technical methods to handle boundary terms in the presence of far-field vacuum and non-positive definite boundary matrices represents a significant advancement over existing theories for IBVPs in exterior domains. The work demonstrates that strong solutions exist locally even when the initial density vanishes at infinity, provided the density remains strictly positive in a neighborhood of the boundary.