REVIEW 2 major objections 5 minor 20 references
Global existence and optimal time-decay rates of the compressible Navier-Stokes equations with density-dependent viscosities
T0 review · 2 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read This paper establishes global classical solutions and optimal decay rates for the 3D isentropic compressible Navier–Stokes equations with density-dependent viscosities when initial perturbations are small in L1∩L2 and α is close to 1…
desk verdict The main theorem is not proved as written: the bootstrap criterion (3.1) demands L∞ smallness at t=0, which the L1∩L2 assumption does not imply, and the linear estimates (3.42) and (3.51) silently drop high-derivative factors. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the Green function G(t,x) of the linearized system, decomposed into a low-frequency heat-like part G_L, a high-frequency regular part G_HR that can bear one derivative and decays exponentially, and a singular high-frequency part G_HS that behaves like $e^{{-ct}}$δ(x) and appears only in the continuity equation. The argument runs through the Duhamel formula V(t)=G(t)*V_0+∫_0^t G(t-s)*N(V)(s)ds, and the bootstrap is a time-decay regularity criterion: assuming ∥(ϱ,u)∥_{L^∞}≤1/2(1+t)^{-3/2} and ∥∇(ϱ,u)∥_{L^∞}≤η(1+t)^{-2}, the paper proves conditional energy estimates needing only $L^{2}$ smallness, then closes the criterion with an auxiliary $L^{{4/3}}$ estimate and the Green function's L^p decay. The viscosity constraint |α-1|≤δ enters through H_α(ϱ)=(1+ϱ)^{α-1}-1, whose L^∞ norm is O(|α-1|), allowing the high-order velocity dissipation terms in Lemma 3.4 to be absorbed.
What would settle it
Take initial data V_0=(ϱ_0,u_0) as a bump of height M and width w with M $w^{{3/2}}$=C_0: the $L^{2}$ norm is C_0 and the $L^{1}$ norm is $C_0^{2}$/M, so for large M the data satisfy (1.7) while violating the pointwise bound at t=0. A check of the linear estimate (3.42) or (3.51) with the factor ∥∇^2V_0∥^{3/4} kept would settle whether the claimed decay holds without an additional high-frequency assumption.
Extended reading notes
Core claim
On its own terms, the central discovery is Theorem 1.1: with (ϱ_0,u_0)∈$H^{4}$, ∥(ϱ_0,u_0)∥_{$L^{1}$∩$L^{2}$}≤C_0 and |α-1|≤δ, the Cauchy problem has a unique global classical solution (ϱ,u)∈C([0,∞);$H^{4}$) with ∇u∈$L^{2}$([0,∞);$H^{4}$), and the optimal uniform decay (1.9)–(1.11) holds. The point of the theorem is that no smallness is imposed on any Sobolev norm of derivatives; the $H^{4}$ norm of the initial data may be arbitrarily large. The mechanism is a bootstrap that starts from the time-decay regularity criterion (3.1), proves conditional energy estimates that need only $L^{2}$ smallness, and then closes the criterion through the Duhamel representation using low-frequency heat-like decay of the Green function and exponential high-frequency estimates. The condition |α-1|≤δ is essential to the closure, since the linearized viscosity perturbation H_α(ϱ)=(1+ϱ)^{α-1}-1 is of size O(|α-1|) and can be absorbed into the velocity dissipation.
Load-bearing premise
The proof's bootstrap needs the initial data to satisfy the pointwise bounds ∥(ϱ_0,u_0)∥_{L^∞}≤1/2 and ∥∇(ϱ_0,u_0)∥_{L^∞}≤η, but smallness of the $L^{1}$∩$L^{2}$ norm alone does not force those bounds when higher-order derivatives are large.
Editorial extensions
If this is right
- Global classical solutions exist for all time for small L^1∩L^2 perturbations of a constant-density state even when the H^4 norms of initial derivatives are large.
- The density and velocity approach the equilibrium at heat-like rates: ∥(ϱ,u)∥_{L^∞}≤C(1+t)^{-3/2} and ∥∇(ϱ,u)∥_{L^∞}≤C(1+t)^{-2}.
- The H^k norms for k=0,1,2 decay at the rates ∥∇^k(ϱ,u)∥_{L^2}≤C(1+t)^{-3/4-k/2}.
- The L^p interpolation bounds (1.11) hold for 2≤q≤∞ with k'=0,1.
- The theorem also gives uniqueness of the global classical solution, so the asymptotic state is determined by the low-frequency part of the data even when high-frequency derivatives are large.
Reading between the lines
- An extension the authors do not pursue: the auxiliary L^{4/3} estimate in Lemma 3.5 could likely be replaced by any L^q with q<2, since it only supplies a small algebraic decay surplus; the final rates would be unchanged and the interpolation might simplify.
- If the mechanism is robust, the same Green-function-plus-regularity-criterion bootstrap should transfer to other systems whose linearization has a heat-like low-frequency block, such as viscous shallow water equations, giving analogous no-small-derivative decay results.
- Because the bootstrap requires the criterion (3.1) at t=0, the advertised regime of arbitrarily large derivatives is fully valid only if a short-time estimate can force the initial pointwise bounds; supplying such an estimate would complete the theorem exactly as stated.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the three-dimensional isentropic compressible Navier-Stokes equations with density-dependent viscosities μ=ρ^α, λ=ρ^α near the constant state. Theorem 1.1 claims global existence of classical solutions and optimal decay rates for H^4 initial data under only smallness of the L1∩L2 norm of the perturbation and |α−1| small, with no smallness imposed on derivatives (Remark 1.2). The proof combines a Green's function decomposition (low-frequency heat-like part, regular high-frequency part, and singular part) with conditional energy estimates under a time-decay bootstrap assumption (3.1), an auxiliary L^{4/3} estimate, and a continuity argument in Section 4. The advertised novelty is that the Sobolev norms of derivatives of the initial data may be arbitrarily large.
Significance. If valid, the result would be a substantial improvement over the classical Matsumura-Nishida and Guo-Wang frameworks and over the recent Luo-Yang result for density-dependent viscosities, since it would remove derivative smallness. The paper contains a careful Green's function decomposition, detailed conditional energy estimates, and an interesting L^{4/3} auxiliary estimate. However, the bootstrap cannot be initialized under the stated hypotheses, and the linear estimates silently discard uncontrolled high-derivative factors. These issues concern the central claim, so the theorem is not established as written.
major comments (2)
- [§3.2 and §4.1, Eq. (3.1), (3.28), (1.7)] The bootstrap assumption (3.1) is required to hold on [0,T], hence in particular at t=0 it demands ∥(ϱ0,u0)∥_{L∞} ≤ 1/2 and ∥∇(ϱ0,u0)∥_{L∞} ≤ η. The hypotheses of Theorem 1.1 only assume (1.7), i.e. smallness in L1∩L2, and Remark 1.2 explicitly allows the Sobolev norms of derivatives to be arbitrarily large. Smallness in L1∩L2 does not imply either of the required L∞ bounds: for a bump ϱ0(x)=M φ(x/w) with fixed profile φ, one has ∥ϱ0∥_{L2}=M w^{3/2} and ∥ϱ0∥_{L1}=M w^3; with M=C0 w^{-3/2}, the L2 norm equals C0 and the L1 norm tends to 0 as w→0, while ∥ϱ0∥_{L∞}=C0 w^{-3/2} and ∥∇ϱ0∥_{L∞}∼C0 w^{-5/2} diverge. Consequently the quantity T* defined in §4.1 can be 0, and the continuity argument cannot be started. The proof needs either a separate short-time argument establishing (3.1) under (1.7) or a change of the main hypotheses; as written, the advertised result is not proven.
- [Eq. (3.42) and (3.51)] In the linear estimate M1 in (3.42), the proof bounds ∥GHS(t)∗V0∥_{L∞} ≤ e^{-Ct}∥V0∥_{L∞} ≤ e^{-Ct}∥V0∥_{L2}^{1/4}∥∇²V0∥_{L2}^{3/4} and then replaces the right side by C C0^{1/4}(1+t)^{-3/2}. The factor ∥∇²V0∥_{L2}^{3/4} is dropped without justification. Similarly, in (3.51) the terms ∥∇V0∥_{L2} and ∥∇V0∥_{L∞} are interpolated as ∥V0∥_{L2}^{1/2}∥∇²V0∥_{L2}^{1/2} and ∥V0∥_{L2}^{1/6}∥∇³V0∥_{L2}^{5/6}, and the high-derivative factors are then replaced by constants C0^{1/2} and C0^{1/6}. Under the hypotheses of Theorem 1.1 these high-derivative norms are not controlled by C0; they may be arbitrarily large. These dropped factors are exactly what would make the linear terms small, so the closing of the bootstrap in Lemmas 3.6 and 3.7 is not justified.
minor comments (5)
- [Theorem 1.1, Eq. (1.11)] The statement says 'For 2 ≤ q ≤ ∞' but the estimate uses L^p; the index p is not defined and the relation between q and p is not given.
- [§3.1, Lemma 3.5] In the K2 estimate, the displayed powers of C0 do not follow from (3.8): since (3.8) gives ∥V∥_{L2} ≤ C C0^{1/2}, the product ∥V∥_{L2}∥∇V∥_{L2} is O(C0^{11/10}) rather than O(C0^{8/5}). The exponents should be recomputed or clarified.
- [§2 and §3] There are several typos: 'Naiver-Stokes' appears in §1 and in references [3,16]; 'Gargliardo-Nirenberg' appears in Lemma 3.5; 'The the general solution' appears in §2; and in (3.17) the sign of the (λ+μ)∇divu term differs from (1.6).
- [§4.2] The exponents α1 and α2 in (4.2) are asserted to be 'verified in Section 3.2', but no explicit verification is given; this should be supplied.
- [Remark 4.1] The remark says that only optimal decay for ∥∇^k u∥_{L^2} is obtained, whereas Theorem 1.1 and §4.2 state decay for ∇^k(ϱ,u); the text should be reconciled.
Circularity Check
No significant circularity: the paper's bootstrap closure uses independent Green's-function and energy estimates, and its decay conclusions are not re-imported as inputs; the flagged initialization gap for the regularity criterion at t=0 is a proof gap, not a circular step.
full rationale
The derivation chain for Theorem 1.1 is a standard bootstrap: the time-decay regularity criterion (3.1) is assumed on [0,T], conditional energy estimates (Lemmas 3.3–3.4, (3.8), (3.16)) and an L^{4/3} auxiliary estimate (Lemma 3.5, (3.31)) are derived from it together with Lemma 2.2's Green's-function estimates, and these close the criterion with improved constants 1/4 and eta/2 (Lemmas 3.6–3.7, (3.40), (3.49)). The decay rates in (1.9) are not fed back into the linear part: (3.42) and (3.51) bound the linear terms by ||V0||_{L1} with the intrinsic (1+t)^{-3/2} and (1+t)^{-2} rates, and the nonlinear integrals are controlled by the energy estimates. No parameter is fitted; C0 and delta are smallness constants, and the H^k decay rates (1.10) in Section 4.2 are derived after the bootstrap from the same estimates, not used as inputs. The Green's-function estimates are cited to external works [3,7,8,16]; local well-posedness is cited to [20]. The only self-citations ([5], [6]) appear in the introduction's literature review and are not load-bearing. Flagged but not circular: the criterion (3.1) is imposed at t=0 ('||(rho,u)||_{L∞} <= 1/2(1+t)^{-3/2}, ||grad(rho,u)||_{L∞} <= eta(1+t)^{-2}'), while the hypothesis (1.7) with 'the Sobolev norms of the spatial derivatives of the initial data may be arbitrarily large' (Remark 1.2) does not imply it; estimates (3.42) and (3.51) drop the factors ||grad^2 V0||^{3/4}_{L2} and ||grad^3 V0||^{5/6}_{L2}, so T* defined in Section 4.1 may be 0 and the continuity argument may never start. This is a missing argument for initialization, not an equivalence of inputs and outputs, so it does not raise the circularity score.
Assumptions & free parameters
free parameters (3)
- C0 =
sufficiently small, not quantified
- delta =
sufficiently small, depends on μ, λ, γ
- eta =
sufficiently small, not quantified
assumptions (4)
- domain assumption Green's function estimates in Lemma 2.2: low-frequency part behaves like the heat kernel, high-frequency regular part decays exponentially, and the singular part satisfies Lp bounds.
- domain assumption Local well-posedness for the density-dependent viscosity system in H^4 is taken from Zhang-Zhao [20].
- ad hoc to paper The bootstrap regularity criterion (3.1) is assumed to hold on [0,T] before being improved.
- standard math Sobolev embedding, Gagliardo-Nirenberg interpolation, commutator estimates from Lemma 3.2, and Gronwall's inequality.
Cite this review
Pith. "Pith review of Global existence and optimal time-decay rates of the compressible Navier-Stokes equations with density-dependent viscosities." pith.science (2026). https://pith.science/paper/DDL6AORZ
@misc{pith2026250716436,
author = {Pith},
title = {Pith review of: Global existence and optimal time-decay rates of the compressible Navier-Stokes equations with density-dependent viscosities},
year = {2026},
howpublished = {\url{https://pith.science/paper/DDL6AORZ}},
note = {Machine review of arXiv:2507.16436}
}
abstract
This paper is devoted to studying the Cauchy problem for the three-dimensional isentropic compressible Navier-Stokes equations with density-dependent viscosities given by $\mu=\rho^\alpha,\lambda=\rho^\alpha(\alpha>0)$. We establish the global existence and optimal decay rates of classical solutions under the assumptions of small initial data in $L^1(\mathbb{R}^3)\cap L^2(\mathbb{R}^3)$ and the viscosity constraint $|\alpha-1|\ll 1$. The key idea of our proof lies in the combination of Green's function method, energy method and a time-decay regularity criterion. In contrast to previous works, the Sobolev norms of the spatial derivatives of the initial data may be arbitrarily large in our analysis
Reference graph
Works this paper leans on
-
[10]
Z. Luo, W. Yang, Optimal decay rates for the compressible Navier–Stokes equations with density–dependent viscosities, J. Differ. Equ. 441, 113483 (2025)
work page 2025
-
[4]
Y. Guo, Y. Wang, Decay of dissipative equations and negative Sobolev space,Commun. Partial Differ. Equ. 37, 2165-2208 (2012)
work page 2012
- [1]
- [2]
-
[3]
Q. Chen, C. Miao, Z. Zhang, Global well-posedness for compressible Naiver-Stokes equations with highly oscillating initial velocity, Commun. Pure Appl. Math. 9, 1173- 1224 (2010)
work page 2010
-
[5]
X. Huang, J. Li, R. Zhang, Global large strong solutions to the compressible Navier- Stokes equations with density-dependent viscosities, case I: isentropic flows. arXiv: 2408.04305v1
- [6]
-
[7]
D. Hoff, K. Zumbrun, Multi-dimensional diffusion waves for the Navier-Stokes equations of compressible flow, Indiana Univ. Math. J. 44, 603-676 (1995)
work page 1995
Show all 20 references
-
[8]
D. Hoff, K. Zumbrun, Pointwise decay estimates for multidimensional Navier-Stokes diffusion waves, Z. Angew. Math. Phys. 48, 597-614 (1997)
1997
-
[9]
Li, The Green’s function of the Navier-Stokes equations for gas dynamics in R3, Comm
D. Li, The Green’s function of the Navier-Stokes equations for gas dynamics in R3, Comm. Math. Phys. 257, 579-619 (2005)
2005
-
[11]
J. Li, Z. Xin, Global existence of weak solutions to the barotropic compressible Navier- Stokes flows with degenerate viscosities. arXiv: 1504.06826v2
-
[12]
T. Liu, W. Wang, The pointwise estimates of diffusion waves for the Navier-Stokes equations in odd multi-dimensions. Comm. Math. Phys. 196, 145-173 (1997)
1997
-
[13]
Matsumura, T
A. Matsumura, T. Nishita, The initial value problem for the equations of motion of viscous and heat-conductive gases, J. Math. Kyoto Univ. 20 (1), 67-104 (1980)
1980
-
[14]
Majda, A.-L
A. Majda, A.-L. Bertozzi, Vorticity and Incompressible Flow, Cambridge University Press, Cambridge, 2002
2002
-
[15]
Mellet, A
A. Mellet, A. Vasseur, On the barotropic compressible Navier-Stokes equations, Com- mun. Partial Differ. Equ. 32 (1-3), 431-452 (2007). 27
2007
-
[16]
J. Qi, W. Wang, Global solutions to the Cauchy problem of viscous shallow water equations in some classes of large data, Discrete Contin. Dyn. Syst. 44 (11), 3572–3598 (2024)
2024
-
[17]
Vasseur, C
A. Vasseur, C. Yu, Existence of global weak solutions for 3D degenerate compressible Navier-Stokes equations, Invent. Math. 206 (3), 935-974 (2016)
2016
-
[18]
Yu, Global existence of strong solutions to the 3d isentropic compressible Navier–stokes equations with density-dependent viscosities, Math
H. Yu, Global existence of strong solutions to the 3d isentropic compressible Navier–stokes equations with density-dependent viscosities, Math. Models Methods Appl. Sci. 46, 10123-10136 (2023)
2023
-
[19]
Zeng, Asymptotic behavior of compressible isentropic viscous 1-D flow, Comm
Y. Zeng, Asymptotic behavior of compressible isentropic viscous 1-D flow, Comm. Pure Appl. Math. 47, 1053-1082 (1994)
1994
-
[20]
Zhang, J
P. Zhang, J. Zhao, The existence of local solutions for the compressible Navier-Stokes equations with the density-dependent viscosities, Commun. Math Sci. 12, 1277–1302 (2014). 28
2014
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.