🔢 mathematics

Expanding Solutions to Free Boundary 3D Spherically Symmetric Compressible Navier-Stokes-Poisson Equations near the Lane-Emden Stars

This paper establishes the existence of global weak solutions for spherically symmetric compressible Navier-Stokes-Poisson equations modeling viscous polytropic gaseous stars with adiabatic exponents $\gamma \in (\frac{6}{5}, \frac{4}{3}]$ near Lane-Emden stars, while also demonstrating the strong instability of these solutions by proving that the support of any strong solution expands to infinity.

Original authors: Han Cao

Published 2026-06-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Han Cao

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 a giant, glowing ball of gas floating in space, held together by its own gravity. This is a star. In this paper, the author, Han Cao, studies what happens to these stars when they are "sticky" (viscous) and when they are allowed to expand or shrink, with their edges touching the empty vacuum of space.

Here is a breakdown of the paper's story using simple analogies:

1. The Setup: A Sticky, Self-Gravitating Balloon

Think of the star as a massive, self-gravitating balloon.

Gravity: The gas wants to collapse inward, like a deflating balloon.
Pressure: The gas pushes outward, like air trying to escape.
Viscosity (Stickiness): The gas isn't a perfect fluid; it has "thickness" or friction (like honey vs. water). This is the Navier-Stokes part.
The Vacuum: The star doesn't have a hard shell. Its edge is a "free boundary" where the gas density drops to zero, meeting the empty void.

The paper asks: If we start with a star that is almost stable, will it stay that way, or will it explode outward?

2. The "Goldilocks" Zone: The Lane-Emden Star

There is a special, perfectly balanced shape for these stars called the Lane-Emden star. It's like a "Goldilocks" configuration where gravity and pressure are perfectly matched.

The paper focuses on a specific type of gas behavior (controlled by a number called $\gamma$ ).
If the gas behaves in a certain way (specifically when $\gamma$ is between $1.2$ and $1.33$), the paper proves that if you start with a star that is slightly different from this perfect balance, it doesn't just wobble—it runs away.

3. The Main Discovery: The Great Escape

The most exciting result is about instability.

The Claim: If you have a star with a specific amount of "stickiness" (viscosity) and you nudge it slightly away from the perfect Lane-Emden shape, it won't settle back down. Instead, the star will start to expand.
The Metaphor: Imagine a spring-loaded toy that is balanced on a knife-edge. If you push it even a tiny bit, it doesn't fall back; it shoots off. The paper proves that for these specific stars, the "push" of the initial conditions causes the star's edge to accelerate outward forever.
The Result: The paper proves that the radius of the star (how big it is) will grow to infinity. It doesn't just get a little bigger; it expands at a predictable, algebraic rate (like a cube root of time). The star effectively dissolves into the universe.

4. The "Recipe" for the Proof

To prove this, the author had to solve a very difficult math puzzle involving equations that describe how the gas moves and how gravity pulls on it.

The Problem: The edge of the star is tricky. As the gas gets thinner near the edge, the math gets "squeezed" and hard to handle. Also, the center of the star is a singularity (a point where things get weird).
The Solution: The author used a "scaffolding" approach.
1. Approximation: First, they pretended the star had a tiny, hard core (cutting off the very center) and a slightly fuzzy edge to make the math easier.
2. Energy Control: They used a clever "energy accounting" trick. They showed that as long as the star starts with less mass than a critical limit (or is in a specific "safe zone" of energy), the negative energy from gravity can't crush the star. Instead, the internal energy wins, pushing the star outward.
3. The Limit: They then slowly removed the "scaffolding" (letting the core get smaller and the edge get sharper) to show that the solution still exists and behaves the same way.

5. The "Expanding Rate"

The paper doesn't just say "it expands." It calculates how fast.

For the specific type of gas studied, the radius of the star grows roughly like the cube root of time (or a quarter root for a specific critical case).
Analogy: If you watch the star for 8 seconds, it's twice as big as it was at 1 second. If you watch for 27 seconds, it's three times as big. It's a steady, unstoppable growth.

Summary of What the Paper Actually Says

Existence: It proves that "weak solutions" (mathematical descriptions of the star's behavior that allow for some roughness) exist for all time. The star doesn't just disappear or break the math; it keeps evolving.
Instability: It proves that the stable-looking Lane-Emden stars are actually unstable in this specific context. If you are near them, you will drift away, and the star will expand forever.
Conditions: This happens for gases where the "stiffness" ( $\gamma$ $γ$ ) is between $1.2$ and $1.33$.
- If the star is exactly at the critical mass limit ( $\gamma = 4/3$ ), the initial mass must be less than the maximum possible mass for a stable star.
- If the star is slightly less stiff ( $\gamma < 4/3$ ), the initial conditions must be in a specific "invariant set" (a mathematical safe zone) near the Lane-Emden solution.

In a nutshell: This paper mathematically proves that certain sticky, self-gravitating gas stars are like a house of cards on a shaky table. If you start them in a specific, nearly-perfect state, they won't stay put; they will inevitably blow up and expand into the infinite void of space.

Technical Summary: Expanding Solutions to Free Boundary 3D Spherically Symmetric Compressible Navier-Stokes-Poisson Equations near the Lane-Emden Stars

Problem Formulation
This paper investigates the global existence and long-time dynamics of spherically symmetric solutions to the compressible Navier-Stokes-Poisson (NSP) equations with a free boundary in three dimensions. The system models viscous polytropic gaseous stars with the equation of state $P(\rho) = K\rho^\gamma$ , where $\gamma \in (6/5, 4/3]$ . The equations describe the evolution of density $\rho$ and velocity $u$ under self-gravity and viscous stress, with the domain expanding or contracting as a free boundary $a(t)$ where the density vanishes continuously into a vacuum.

The study focuses on two specific regimes:

Mass-critical case ( $\gamma = 4/3$ ): The initial mass is required to be strictly less than the mass of the Lane-Emden star (Chandrasekhar mass).
Sub-critical case ( $\gamma \in (6/5, 4/3)$ ): The initial data must belong to a specific invariant set $I$ defined by variational constraints, which includes initial data near the Lane-Emden stars but allows for configurations not covered by previous literature.

Methodology
The authors employ a combination of variational methods, a priori energy estimates, and compactness arguments to construct global weak solutions.

Variational Control of Gravitational Energy:
A central challenge in the range $\gamma \in (6/5, 4/3]$ is controlling the negative gravitational energy, which can dominate the internal energy. Drawing on the variational methods established in [7], the authors utilize the properties of Lane-Emden solutions to derive lower bounds for the functional $S_\mu(\rho)$ .
- For $\gamma = 4/3$ , the condition $M(\rho_0) < M_{ch}$ ensures the gravitational energy is controlled by the internal energy.
- For $\gamma \in (6/5, 4/3)$ , the authors define an invariant set $I$ where the functional $Q(\rho)$ remains positive. They demonstrate that this set includes initial data near Lane-Emden stars that were previously excluded in works like [25], specifically those with small kinetic energy but specific density profiles.
Approximation and Regularization:
To handle the vacuum boundary and the singularity at the origin, the authors consider an approximate problem on a domain $[\xi, a(t)]$ with $\xi > 0$ . They establish the existence of local strong solutions for this approximate problem and then extend them globally.
A Priori Estimates:
The proof relies on deriving uniform estimates independent of the regularization parameter $\xi$ :
- Basic Energy Estimate: Establishes the conservation of total energy minus the dissipation due to viscosity.
- Higher-Order Estimates: The authors derive higher space-time integrability for the density ( $\rho^{2\gamma}$ ) and weighted $L^2$ estimates for derivatives of velocity and density in Lagrangian coordinates.
- Boundary and Origin Separation: The spatial domain is split into a region near the origin and a region near the free boundary.
  - Near the free boundary, higher regularity allows for strong convergence via the Arzela-Ascoli theorem.
  - Near the origin, where regularity is lower, the authors adapt techniques from [22] (originally for Navier-Stokes on $\mathbb{R}^n$ ) to prove the strong convergence of the density. This involves proving that the weak limit of the pressure term $\rho^\gamma$ equals the pressure of the weak limit, utilizing the convexity of the pressure function and specific integral inequalities.
Compactness and Limit Passage:
By combining the strong convergence near the boundary and the origin, the authors extract a subsequence that converges to a global weak solution satisfying the distributional formulation of the NSP equations.

Key Results

Global Existence of Weak Solutions (Theorem 2.1 & 2.2):
The paper proves the existence of global weak solutions for $\gamma \in (6/5, 4/3]$ under the specified mass or invariant set conditions.
- Theorem 2.1: Establishes the existence of a global weak solution $(\rho, u, a)$ with non-negative density and bounded energy.
- Theorem 2.2: Provides higher-order regularity results. Specifically, it establishes interior regularity ( $H^1$ for density, $H^2$ for velocity) away from the origin and boundary, and specific regularity near the free boundary, including bounds on the effective viscous flux $\sigma$ .
Expansion of the Free Boundary (Theorem 2.3):
Assuming the existence of strong solutions, the paper proves that the support of the solution expands to infinity at an algebraic rate. This indicates the strong instability of Lane-Emden solutions in this regime.
- For $\gamma \in (6/5, 4/3)$ , the radius $a(t)$ satisfies $a(t) \approx (t / (\eta + 4/3\varepsilon))^{1/3}$ .
- For $\gamma = 4/3$ , the expansion rate is bounded by $(1+t)^{1/4} \lesssim a(t) \lesssim (t / (\eta + 4/3\varepsilon))^{1/3}$ .
- The proof utilizes a Virial-type argument involving the functional $H(t) = \int \rho r^4 dr$ and its time derivatives.
Density-Dependent Viscosity (Remark 2.3):
The authors extend the expansion rate analysis to cases with density-dependent viscosity $\varepsilon(\rho), \eta(\rho) \propto \rho^\alpha$ . They show that the expansion rate depends on the exponent $\alpha$ :
- If $0 \le \alpha < 2/3$ , the rate is $a(t) \gtrsim (t/\eta)^{1/(3(1-\alpha))}$ .
- If $2/3 \le \alpha \le \gamma$ , the rate is linear, $a(t) \gtrsim t$ .
- This highlights that the vanishing of bulk viscosity significantly alters the expansion dynamics when $\alpha$ is small.

Significance and Claims
The paper claims to resolve the global existence of weak solutions for the free-boundary NSP equations in the critical and sub-critical adiabatic exponent ranges ( $\gamma \in (6/5, 4/3]$ ) with constant viscosity. A primary contribution is the refinement of the initial data set for $\gamma \in (6/5, 4/3)$ ; the authors show that their invariant set $I$ is broader than previous constraints (e.g., [25]), allowing for initial data near Lane-Emden stars that were previously considered out of reach.

Furthermore, the paper establishes the algebraic expansion of the free boundary, providing a rigorous mathematical confirmation of the instability of Lane-Emden stars in this regime. The authors note that while the expansion rate is proven for strong solutions, the global existence of such strong solutions remains an open problem, making the instability result conditional. The work bridges the gap between the stability results for $\gamma > 4/3$ and the instability results for $\gamma < 4/3$ in the context of viscous gaseous stars.