Global existence of classical solutions for the multi-dimensional compressible Navier-Stokes-Poisson equations on solid balls for arbitrary spherically symmetric large initial data
Pith reviewed 2026-05-10 17:40 UTC · model grok-4.3
The pith
The compressible Navier-Stokes-Poisson equations admit global classical solutions for arbitrarily large spherically symmetric initial data on solid balls.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We establish the global existence of spherically symmetric classical solutions to the compressible Navier-Stokes-Poisson equations for both gaseous stars and plasmas with arbitrarily large initial data on solid balls, when the viscosity coefficients satisfy the BD-type entropy equality μ=ρ^α, λ=(α-1)ρ^α with N=2 and α in (1/2,1] or N=3 and α in (5/6,1]. The proof controls the central singularity by bounding the growth of density and gravitational potential, then uses the BD entropy together with the coupling between effective velocity and velocity to obtain L^∞ estimates that deliver upper and lower bounds on density.
What carries the argument
The BD-type entropy equality from the density-dependent viscosities μ=ρ^α and λ=(α-1)ρ^α, which supplies structural control, combined with spherical symmetry to handle the central singularity and derive density bounds via effective-velocity coupling.
If this is right
- Classical solutions exist for all positive time without breakdown.
- Density remains strictly positive and bounded from above and below.
- The gravitational potential remains controlled for all time.
- The result applies simultaneously to gaseous-star and plasma models.
- The same control works in both two and three space dimensions under the given viscosity ranges.
Where Pith is reading between the lines
- The same entropy-plus-symmetry strategy might extend to other radially symmetric domains or to the non-symmetric case with small perturbations.
- Numerical schemes that preserve the BD entropy could reliably simulate large-data regimes without artificial blow-up.
- The central-singularity estimates suggest that similar bounds could be obtained for related systems such as compressible MHD with gravity.
Load-bearing premise
The initial data must be spherically symmetric on the solid ball and the viscosities must obey exactly μ=ρ^α, λ=(α-1)ρ^α for α in the stated intervals.
What would settle it
A concrete spherically symmetric initial datum of arbitrary size on the ball for which the density reaches zero or infinity in finite time, violating the derived L^∞ bounds.
read the original abstract
Whether the 3D compressible Navier-Stokes-Poisson equations admit global classical solutions for general large initial data has long been a challenging open problem. In this paper, we provide an affirmative answer to this question under spherical symmetry on solid balls . Specifically, we consider the initial-boundary value problem for the multi-dimensional compressible equations with density-dependent viscosity coefficients satisfying the BD-type entropy equality, namely, assuming $\mu=\rho^{\alpha},\ \lambda=(\alpha-1)\rho^{\alpha}$ with $N=2, \alpha\in (\frac{1}{2},1]$ and $N=3, \alpha\in (\frac{5}{6},1]$, we establish the global existence of spherically symmetric classical solutions to the compressible Navier-Stokes-Poisson equations for both gaseous stars and plasmas with arbitrarily large initial data on solid balls. Our key observation lies in successfully handling the singularity at the center of the ball. By controlling the growth orders of the density and the gravitational potential at the central singularity, leveraging the structural advantages of the BD entropy and spherical symmetry, and fully exploiting the coupling between the effective velocity and the velocity, we establish $L^\infty$ estimates for the key quantities, which in turn yield upper and lower bound estimates for the density. This can be regarded as the first result on the existence of global classical solutions for arbitrarily large initial data to the compressible Navier-Stokes-Poisson equations in a truly multi-dimensional domain with high-dimensional features.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proves global-in-time existence of classical spherically symmetric solutions to the compressible Navier-Stokes-Poisson system on the solid ball in dimensions N=2 and N=3. The viscosity coefficients are required to satisfy the BD entropy relation μ=ρ^α, λ=(α-1)ρ^α with α∈(1/2,1] (N=2) or α∈(5/6,1] (N=3). The result holds for arbitrarily large initial data, for both the gaseous-star (attractive) and plasma (repulsive) cases, by reducing via spherical symmetry, controlling the origin singularity through growth-order estimates on density and gravitational/electrostatic potential, and closing L^∞ bounds via the effective-velocity formulation.
Significance. If the central estimates are correct, the result is a substantial advance: it supplies the first global classical existence theorem for large data in a genuinely multi-dimensional domain (solid ball) for the NSP system. The structural exploitation of the BD entropy, spherical symmetry, and effective-velocity coupling to obtain uniform bounds without smallness assumptions is a clear technical strength and may influence related problems with singular coefficients or geometric singularities.
major comments (2)
- §3 (a priori estimates): the passage from the BD entropy identity to the pointwise growth-order control on ρ and Φ near r=0 must be verified in detail; the argument appears to rely on integrating the effective-velocity equation against a carefully chosen test function, but it is unclear whether the resulting differential inequality closes uniformly up to the origin without an additional logarithmic correction or cutoff.
- Theorem 1.2 and the continuation criterion (presumably §4): once the L^∞ bound on density is obtained, the paper must confirm that the standard Beale-Kato-Majda-type criterion for the NSP system is satisfied; the current sketch leaves open whether the time-integral of ||∇u||_∞ remains finite solely from the effective-velocity estimates.
minor comments (2)
- The precise statement of the initial-boundary conditions on the solid ball (no-flux or Dirichlet for velocity, etc.) should be written explicitly in the introduction rather than deferred to §2.
- Notation: the effective velocity w = u + ∇(something) is introduced without a numbered equation; adding an equation label would improve readability when it is later coupled to the density equation.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation of our work and the constructive comments, which help strengthen the presentation. We address the two major comments point by point below, providing clarifications and indicating the revisions made to the manuscript.
read point-by-point responses
-
Referee: §3 (a priori estimates): the passage from the BD entropy identity to the pointwise growth-order control on ρ and Φ near r=0 must be verified in detail; the argument appears to rely on integrating the effective-velocity equation against a carefully chosen test function, but it is unclear whether the resulting differential inequality closes uniformly up to the origin without an additional logarithmic correction or cutoff.
Authors: We appreciate the referee drawing attention to the details of the origin estimates. In the original manuscript, Lemma 3.2 derives the growth-order bounds by testing the effective-velocity equation with a radial cutoff function that is identically 1 away from a small neighborhood of r=0 and transitions smoothly; the BD entropy identity supplies the necessary integrability to absorb the commutator terms without logarithmic factors. The resulting differential inequality for the weighted integrals of ρ and |∇Φ| closes uniformly because the effective velocity damps the singular contributions at the origin under the given range of α. To remove any ambiguity, we have expanded the proof in the revised Section 3 with an explicit computation of all boundary terms at r=0 and a remark explaining why no additional cutoff or log correction is required. revision: yes
-
Referee: Theorem 1.2 and the continuation criterion (presumably §4): once the L^∞ bound on density is obtained, the paper must confirm that the standard Beale-Kato-Majda-type criterion for the NSP system is satisfied; the current sketch leaves open whether the time-integral of ||∇u||_∞ remains finite solely from the effective-velocity estimates.
Authors: We thank the referee for this observation on the continuation argument. After establishing the uniform L^∞ bound on density via the effective-velocity formulation, the revised Section 4 now contains a dedicated paragraph that verifies the Beale-Kato-Majda criterion explicitly. Using the relation between u and the effective velocity together with the already-obtained bounds on density and the Poisson potential, we show that ∫_0^T ||∇u(t)||_∞ dt remains finite; the key step is an L^∞ estimate for ∇u that is integrable in time by virtue of the damping provided by the BD structure and the spherical symmetry. This confirms that the local solution extends to all time, completing the proof of Theorem 1.2. revision: yes
Circularity Check
No significant circularity detected in the proof strategy
full rationale
The paper is a pure existence proof for classical solutions of the compressible Navier-Stokes-Poisson system under spherical symmetry on a solid ball. It assumes explicit BD-type viscosity coefficients μ=ρ^α, λ=(α-1)ρ^α with given ranges for α, reduces via symmetry, and derives a priori L^∞ bounds on density and potential by controlling growth orders at the origin together with effective-velocity coupling. All steps are forward estimates from the PDE structure and entropy identity; no quantity is defined in terms of itself, no fitted parameter is relabeled as a prediction, and no load-bearing step collapses to a self-citation or ansatz smuggled from prior work by the same authors. The result is therefore independent of its own inputs.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Viscosity coefficients satisfy μ=ρ^α, λ=(α-1)ρ^α with the stated ranges of α for N=2,3, yielding the BD entropy equality.
- domain assumption Initial data and solution remain spherically symmetric on the solid ball domain.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
assuming μ=ρ^α, λ=(α−1)ρ^α with N=2, α∈(1/2,1] and N=3, α∈(5/6,1], … effective velocity w=u+r^{N−1}(ρ^α)_x … weighted estimates … higher integrability … Sobolev embedding
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
spherically symmetric classical solutions … on solid balls … singularity at the center … Alexander duality not invoked
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
-
[4]
P. Bourguignon, H. Brezis: Remarks on the Euler equation, J. Funct. Anal. 15 (4) (1974) 341–363
work page 1974
-
[5]
Q. Duan, H.-L. Li: Global existence of weak solution for the compressible Navier- Stokes-Poisson system for gaseous stars. J. Differ. Equ. 259, 5302–5330 (2015) 60
work page 2015
-
[6]
D. Hoff: Global solutions of the Navier-Stokes equations for multidimensional com- pressible flow with discontinuous initial data. J. Differ. Equ. 120(1), 215-254 (1995)
work page 1995
-
[7]
X. D. Huang, W. L. Meng, X. Y. Zhang: On global classical and weak solutions with arbitrary large initial data to the multi-dimensional viscous Saint-Venant system and compressible Navier-Stokes equations subject to the BD entropy condition under spher- ical symmetry, arXiv:2512.15029 (2025)
-
[8]
J. Li, Z. P. Xin: Global existence of weak solutions to the barotropic compressible Navier-Stokes flows with degenerate viscosities. arXiv:1504.06826 (2015)
work page Pith review arXiv 2015
-
[9]
T. Luo, Z. Xin, H. Zeng: On nonlinear asymptotic stability of the Lane-Emden solutions for the viscous gaseous star problem. Adv. Math. 291, 90–182 (2016)
work page 2016
-
[10]
T. Luo, Z. Xin, H. Zeng: Nonlinear asymptotic stability of the Lane-Emden solutions for the viscous gaseous star problem with degenerate density dependent viscosities. Commun. Math. Phys. 347, 657–702 (2016)
work page 2016
- [11]
-
[12]
Z. Tan, Y. Zhang: Strong solutions of the coupled Navier–Stokes–Poisson equations for isentropic compressible fluids. Acta Math. Sci. 30(4), 1280–1290 (2010)
work page 2010
-
[13]
A. F. Vasseur, C. Yu: Existence of global weak solutions for 3D degenerate compressible Navier-Stokes equations. Invent. Math. 206(3), 935–974 (2016)
work page 2016
-
[14]
X. Y. Zhang: Spherically symmetric strong solution of compressible flow with large data and density-dependent viscosities. J. Math. Anal. Appl. 549(2), Paper no. 129488 (2025)
work page 2025
- [15]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.