Pith. sign in

REVIEW 35 references

Formation of Implosion Singularities in 3D Compressible Navier-Stokes-Korteweg Equation

T0 review · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read For alpha less than 1/2, the 3D Navier-Stokes-Korteweg system has smooth initial data whose solutions form a finite-time implosion: the density diverges at a point and the effective velocity blows up.

arxiv 2608.10554 v1 pith:MLV73W2H submitted 2026-08-11 math.AP

Xiangdi Huang , Yongteng Gu This is my paper · ORCID
classification math.AP
keywords alphasmoothsolutionsblowupcitecompressibledatadensity
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper studies a fluid model where friction (viscosity) and surface-tension-like forces (capillarity) depend on fluid density through power laws. The main result is that, when the power alpha is small (alpha less than 1/2), there exist smooth starting states, with no vacuum anywhere, where the fluid nonetheless crushes itself into a single point in a finite time: the density becomes infinite at the origin. The construction starts from known self-similar imploding solutions of the simpler compressible Euler equations, which already describe a fluid squeezing itself. In specially chosen coordinates that zoom in as the collapse time approaches, the extra viscosity and capillarity terms appear multiplied by a rapidly decaying exponential factor. The authors show that these small extra terms do not destroy the implosion. They use weighted energy estimates to control the solution, decompose the linearized problem into a finite set of unstable directions and a stable remainder, and carefully select the initial data so the unstable part stays controlled. The selected initial data form a finite-codimensional family. At the singular time, the density becomes infinite at the origin. The effective velocity of the reformulated system also blows up, although the physical velocity of the original equations may or may not blow up because of a possible cancellation with the density gradient.
Extended reading notes

Core claim

Theorem 1.1 states that for (gamma, alpha) in one of two regimes (P1: 1<gamma<1+2√3 and 0<alpha<1/2 with f1>0, or P2: gamma in (1,1+2√3)\J and 0<alpha<alpha*(gamma)), there exist C^infty initial data with inf rho0 > rho0 such that the solution of the reformulated system (1.12) exists on [0,T) and satisfies lim_{t->T-} rho(t,0)=+infinity and lim sup_{|x|<=r}|v(t,x)|=+infinity for every r>0, with the rescaled solution converging to the prescribed self-similar Euler profile. Corollary 1.1 transfers this to the original NSK system for density, while noting that blowup of the physical velocity is not established.

Load-bearing premise

The construction inherits, without proof, the existence and detailed properties of smooth self-similar Euler profiles: Lemma 5.3 (existence of admissible scaling parameters Lambda and profiles) and Lemma 5.4 (positivity, decay, and repulsivity estimates (5.11)-(5.15)). It also relies on the spectral decomposition of the truncated linearized operator L in Lemma 5.7 (finite-dimensional unstable subspace, spectral gap c_g, stable semigroup decay), which is imported from [10]. These cited results are load-bearing for the bootstrap and unstable-mode selection; if they failed, the construction collapses. Location: Appendix 5.2 (Lemmas 5.3-5.4) and 5.3 (Lemma 5.7).

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Assumptions & free parameters 8 free parameters · 5 assumptions · 0 invented entities

The proof introduces no new physical entities. The effective velocity is a mathematical change of variables, not a postulated mechanism. The main external assumptions are the imported Euler-profile existence/repulsivity and the spectral splitting of the truncated operator, both taken from prior literature.

free parameters (8)
  • sigma0 = unspecified small constant
    Controls the size of the initial perturbation; chosen small to close the bootstrap. Appears throughout the parameter hierarchy (3.6).
  • sigma1 = sigma1 = sigma0^(5/4)
    Smaller than sigma0; used for the lower density bound and bootstrap improvements. The relation is imposed in Remark 3.1.
  • tau0 = large, e^{-f1 tau0} << 1
    Initial self-similar time; taken large so the exponentially decaying dissipative terms are perturbative. See (3.13).
  • R0 = determined by sigma0 ~ R0^(1-Lambda)
    Cutoff radius for the weight function phi; chosen so that profile bounds (3.9) hold. See (3.10).
  • R1 = large, satisfying (3.1)-(3.4)
    Truncation radius for the linearized operator L; chosen to make the exterior profile small. See Lemma 3.1.
  • K, m, J = large integers, 1/K << 1/m << eta, J >= 2m, K >= 6
    Sobolev order and truncation parameters for the spectral splitting; chosen to satisfy the hierarchy (3.6)-(3.8) and Lemma 5.5.
  • eta = small exponent in weight phi, 1/K << eta << sigma_g
    Weight exponent in (2.2); smallness ensures weighted Gagliardo-Nirenberg and positivity of transport terms. See (3.6).
  • varpi (with sigma_g = 25/12 varpi) = small, sigma_g << f1
    Decay rate for the unstable modes; related to the spectral gap c_g = sigma_g. The ratio sigma_g/varpi = 25/12 is fixed in Section 3.7.
assumptions (5)
  • domain assumption Existence of smooth self-similar Euler profiles (Lemma 5.3)
    The paper cites [4,29,30,32] for the existence of admissible scaling parameters Lambda and smooth spherically symmetric profiles of (5.8). The blowup construction uses these profiles as the leading-order singularity.
  • domain assumption Positivity, decay, and repulsivity of the profiles (Lemma 5.4, (5.11)-(5.15))
    The bootstrap and high-order energy estimates use the lower bound Sbar >= C^{-1}<r>^{-(Lambda-1)}, the decay (5.12), and the repulsivity conditions (5.13)-(5.15). These are imported from [10,32] without proof.
  • domain assumption Spectral decomposition of the truncated linearized operator L (Lemmas 5.5-5.7)
    The finite-dimensional unstable subspace X_u, the spectral gap c_g, the stable semigroup decay, and the smoothness of X_u are cited from [10]. The unstable-mode selection in Proposition 3.7 depends entirely on this.
  • standard math Weighted Gagliardo-Nirenberg inequalities (Lemma 5.1)
    Used repeatedly for pointwise decay estimates; stated without proof as a standard analytic tool.
  • standard math No continuous retraction of a ball onto its boundary (Brouwer fixed-point theorem)
    Used in Proposition 3.7 to prove the existence of initial unstable coefficients that keep the trajectory inside the moving neighborhood.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Formation of Implosion Singularities in 3D Compressible Navier-Stokes-Korteweg Equation." pith.science (2026). https://pith.science/paper/MLV73W2H

@misc{pith2026260810554,
  author       = {Pith},
  title        = {Pith review of: Formation of Implosion Singularities in 3D Compressible Navier-Stokes-Korteweg Equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MLV73W2H}},
  note         = {Machine review of arXiv:2608.10554}
}
abstract

Previous works of Gu-Huang-Meng-Zhou~\cite{Gu-Huang-Meng-Zhou} and Huang-Lei-Zhou~\cite{Huang-Lei-Zhou} established global strong solutions away from vacuum for arbitrarily large initial data when $\alpha$ lies in a suitable range. In contrast, we show that, for a class of small positive exponents $\alpha$ $(\alpha<\frac{1}{2}$), there exist smooth initial data with density uniformly separated from vacuum whose corresponding solutions develop finite-time implosion singularities. Our construction is based on smooth self-similar imploding profiles of the compressible Euler equations. After reformulating the system in self-similar coordinates, the viscous and capillary effects appear as exponentially decaying perturbations. We control the resulting non-autonomous system through weighted high-order energy estimates, a stable-unstable decomposition of the linearized operator, and a finite-dimensional selection of the unstable components. The constructed solutions remain smooth before the singular time and converge, after rescaling, to the prescribed imploding profile. In particular, at the blowup time $T$, the density becomes infinite at the origin, while the effective velocity $u + d \alpha \rho^{\alpha-2} \nabla \rho$ is unbounded in every neighborhood of the origin. These results complement the aforementioned global existence theory and exhibit a distinct finite-time blowup mechanism for the small-$\alpha$ regime, where the effective bulk-viscosity structure may no longer be positive.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

35 extracted references · 28 canonical work pages

  1. [1]

    Bourgain, Global well-posedness of the defocusing critical nonlinear Schr¨ odinger equa- tion in the radial case, J

    J. Bourgain, Global well-posedness of the defocusing critical nonlinear Schr¨ odinger equa- tion in the radial case, J. Amer. Math. Soc. 12 no. 1 (1999), 145–171

    work page 1999

  2. [2]

    Bresch, B

    D. Bresch, B. Desjardins, C.-K. Lin, On some compressible fluid models: Korteweg, Lu- brication, and Shallow water systems, Comm. Partial Differential Equations 28 no. 2-3 (2003), 843–868

    work page 2003

  3. [3]

    Bresch, B

    D. Bresch, B. Desjardins, Existence of global weak solutions for the 2D viscous shallow- water equations and convergence to the quasi-geostrophic model, Comm. Math. Phys. 238 no. 1-2 (2003), 211–223

    work page 2003

  4. [4]

    Buckmaster, G

    T. Buckmaster, G. Cao-Labora, and J. G´ omez-Serrano, Smooth imploding solutions for 3D compressible fluids, Forum Math. Pi 13 (2025), Paper No. e6, 139pp

    work page 2025

  5. [5]

    Brevik and Ø

    I. Brevik and Ø. Grøn, Universe models with negative bulk viscosity, Astrophys. Space Sci. 347 (2013), 399–404

    work page 2013

  6. [6]

    Buckmaster, S

    T. Buckmaster, S. Shkoller, and V. Vicol, Formation of shocks for 2D isentropic compress- ible Euler, Comm. Pure Appl. Math. 75 (2022), 2069–2120

    work page 2022

  7. [7]

    Buckmaster, S

    T. Buckmaster, S. Shkoller, and V. Vicol, Shock formation and vorticity creation for 3D Euler, Comm. Pure Appl. Math. 76 (2023), 1965–2072

    work page 2023

  8. [8]

    Buckmaster, S

    T. Buckmaster, S. Shkoller, and V. Vicol, Formation of point shocks for 3D compressible Euler, Comm. Pure Appl. Math. 76 (2023), 2073–2191

    work page 2023

Show all 35 references
  1. [9]

    Cort´ azar, M

    C. Cort´ azar, M. del Pino, and M. Musso, Green’s function and infinite-time bubbling in the critical nonlinear heat equation, J. Eur. Math. Soc. (JEMS) 22 no. 1 (2020), 283–344

    2020

  2. [10]

    Cao-Labora, J

    G. Cao-Labora, J. G´ omez-Serrano, J. Shi, and G. Staffilani, Non-radial implosion for compressible Euler and Navier–Stokes inT 3 andR 3, Camb. J. Math. 13(4) (2026), 753- 885

    2026

  3. [11]

    G.-Q. Chen, L. Liu, S. Zhu, Development of Implosions of Solutions to the Three- Dimensional Degenerate Compressible Navier–Stokes Equations, arXiv:2603.10141 (2026)

    2026

  4. [12]

    Colliander, M

    J. Colliander, M. Keel, G. Staffilani, H. Takaoka, and T. Tao, Global well-posedness and scattering for the energy-critical nonlinear Schr¨ odinger equation inR3, Ann. of Math. (2) 167 no. 3 (2008), 767–865

    2008

  5. [13]

    Daskalopoulos, M

    P. Daskalopoulos, M. del Pino, and N. Sesum, Type II ancient compact solutions to the Yamabe flow, J. Reine Angew. Math. 738 (2018), 1–71

    2018

  6. [14]

    Duyckaerts, C

    T. Duyckaerts, C. Kenig, and F. Merle, Universality of the blow-up profile for small radial type II blow-up solutions of the energy-critical wave equation, J. Eur. Math. Soc. (JEMS) 13 no. 3 (2011), 533–599

    2011

  7. [15]

    Duyckaerts, C

    T. Duyckaerts, C. Kenig, and F. Merle, Universality of the blow-up profile for small type II blow-up solutions of the energy-critical wave equation: the nonradial case, J. Eur. Math. Soc. (JEMS) 14 no. 5 (2012), 1389–1454

    2012

  8. [16]

    J. Dunn, J. Serrin, On the thermomechanics of interstitial working, Arch. Rational Mech. Anal. 88 no. 2 (1985), 95–133

    1985

  9. [17]

    M. G. Grillakis, On nonlinear Schr¨ odinger equations, Comm. Partial Differential Equations 25 no. 9–10 (2000), 1827–1844

    2000

  10. [18]

    Y. Gu, X. Huang, Global spherically symmetric classical solutions for arbitrary large ini- tial data of the multi-dimensional non-isentropic compressible Navier–Stokes equations, arXiv:2512.24799 (2025). 89

    2025

  11. [19]

    Y. Gu, X. Huang, W. Meng, H. Zhou, Global-in-time strong solutions for the 2D and 3D generalized compressible Navier–Stokes–Korteweg system with arbitrarily large initial data, arXiv:2603.11762 (2026)

    2026 arXiv

  12. [20]

    Huang, Y

    X. Huang, Y. Gu, M. Lei, Global strong solutions with large initial data for the Cauchy problem of the multi-dimensional compressible Navier–Stokes–Korteweg system, arXiv:2602.10700 (2026)

    2026

  13. [21]

    Huang, M

    X. Huang, M. Lei, H. Zhou, On the Cauchy problem for the multi-dimensional compressible Navier–Stokes–Korteweg system: Global strong solutions with arbitrarily large initial data, arXiv:2604.23652 (2026)

    2026 arXiv

  14. [22]

    Huang, W

    X. Huang, W. Meng, X. 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 spherical symmetry, arXiv:2512.15029 (2025)

    2025

  15. [23]

    Huang, W

    X. Huang, W. Meng, X. Zhang, Global regularity of the multi-dimensional compressible Navier–Stokes–Korteweg system, arXiv:2602.00455 (2026)

    2026

  16. [24]

    Korteweg, Sur la forme que prennent les ´ equations du mouvement des fluides si l’on tient compte des forces capillaires par des variations de densit´ e, Arch

    D. Korteweg, Sur la forme que prennent les ´ equations du mouvement des fluides si l’on tient compte des forces capillaires par des variations de densit´ e, Arch. N´ eer. Sci. Exactes S´ er. II 6 (1901), 1–24

    1901

  17. [25]

    C. E. Kenig and F. Merle, Global well-posedness, scattering and blow-up for the energy- critical focusing non-linear wave equation, Acta Math. 201 no. 2 (2008), 147–212

    2008

  18. [26]

    Krieger, W

    J. Krieger, W. Schlag, and D. Tataru, Slow blow-up solutions for theH 1(R3)-critical focusing semilinear wave equation, Duke Math. J. 147 no. 1 (2009), 1–53

    2009

  19. [27]

    P. D. Lax, Development of singularities of solutions of nonlinear hyperbolic partial differ- ential equations, J. Math. Phys. 5 (1964), 611–613

    1964

  20. [28]

    J. Li, Z. Xin, Global existence of weak solutions to the barotropic compressible Navier– Stokes flows with degenerate viscosities, arXiv:1504.06826 (2015)

    2015 arXiv

  21. [29]

    Merle, P

    F. Merle, P. Rapha¨ el, I. Rodnianski, and J. Szeftel, On the implosion of a compressible fluid I: Smooth self-similar inviscid profiles, Ann. of Math. 196 (2022), 567–778

    2022

  22. [30]

    Merle, P

    F. Merle, P. Rapha¨ el, I. Rodnianski, and J. Szeftel, On the implosion of a compressible fluid II: Singularity formation, Ann. of Math. 196 (2022), 779–889

    2022

  23. [31]

    T. C. Sideris, Formation of singularities in three-dimensional compressible fluids, Comm. Math. Phys. 101 (1985), 475–485

    1985

  24. [32]

    F. Shao, S. Wang, D. Wei, Z. Zhang, Blow-up of the 3-D compressible Navier–Stokes equations for monatomic gases, arXiv:2501.15701v2 (2025)

    2025 arXiv

  25. [33]

    Tankeshwar, Generalized negative bulk viscosity in liquids, J

    K. Tankeshwar, Generalized negative bulk viscosity in liquids, J. Phys. (1994), 9295–9300

    1994

  26. [34]

    Vasseur, C

    A. Vasseur, C. Yu, Existence of global weak solutions for 3D degenerate compressible Navier–Stokes equations, Invent. Math. 206(3) (2016), 935–974

    2016

  27. [35]

    Zhang, Spherically symmetric strong solution of compressible flow with large data and density-dependent viscosities, J

    X. Zhang, Spherically symmetric strong solution of compressible flow with large data and density-dependent viscosities, J. Math. Anal. Appl. 549(2) (2025), Paper no. 129488. 90

    2025

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.