Pith. sign in

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 →

arxiv 2507.16436 v1 pith:DDL6AORZ submitted 2025-07-22 math.AP

classification math.AP MSC 76N1035B45
keywords compressibleNavier-Stokesequationsdensity-dependentviscosityglobalwell-posednessoptimaldecayratesGreen'sfunctionmethodregularitycriterionCauchyproblemclassicalsolutions
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 claims that the three-dimensional isentropic compressible Navier–Stokes equations with density-dependent viscosities μ=λ=ρ^α admit a unique global classical solution whenever the initial perturbation of a constant density state is small in $L^{1}$∩$L^{2}$ and the viscosity exponent satisfies |α-1|≤δ, no matter how large the $H^{4}$ Sobolev norms of spatial derivatives of the data are. The solution is shown to decay like a heat kernel, with ∥(ϱ,u)∥_{L^∞}≤C(1+t)^{-3/2}, ∥∇(ϱ,u)∥_{L^∞}≤C(1+t)^{-2}, and ∥∇^k(ϱ,u)∥_{$L^{2}$}≤C(1+t)^{-3/4-k/2} for k=0,1,2, together with matching L^p rates. This sharpens earlier decay results that required smallness of the whole Sobolev norm of the initial data; if correct, it means the long-time behavior is controlled by the low-frequency, integrable part of the data alone, while high-frequency oscillations may be large. The proof combines a Green's function decomposition of the linearized flow with a time-decay regularity criterion and conditional energy estimates, and the near-constant-viscosity condition α≈1 is used to absorb higher-order dissipation terms.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
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.

Referee Report

2 major / 5 minor

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)
  1. [§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.
  2. [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)
  1. [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.
  2. [§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.
  3. [§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. [§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.
  5. [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

0 steps flagged · score 1.0 of 10

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 3 free parameters · 4 assumptions · 0 invented entities

The proof relies on standard harmonic analysis and PDE tools, on Green's function estimates cited from previous papers, and on a bootstrap criterion that is not implied by the theorem's smallness assumptions. There are no invented particles or fitted physical constants. The critical missing support is an initial L∞ or high-frequency control for the bootstrap.

free parameters (3)
  • C0 = sufficiently small, not quantified
    Smallness threshold for the initial L1∩L2 norm in Theorem 1.1. The proof chooses it small enough to close the bootstrap, but this choice cannot compensate for missing initial L∞ control.
  • delta = sufficiently small, depends on μ, λ, γ
    Restricts |α−1| ≤ δ. It is chosen small enough to absorb Hα(ϱ) viscosity error terms in Lemma 3.4 and the bootstrap estimates; no quantitative value is derived.
  • eta = sufficiently small, not quantified
    Constant in the gradient part of the regularity criterion (3.1). It appears in Proposition 3.1 and is chosen small enough for the closing argument; no numerical value is given.
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.
    The paper sketches the frequency analysis leading to these estimates but refers to [3,7,8,16] for proof details. The entire bootstrap in Section 3.2 relies on these decay rates.
  • domain assumption Local well-posedness for the density-dependent viscosity system in H^4 is taken from Zhang-Zhao [20].
    Used in Section 4.1 to extend the solution beyond T*; no proof is reproduced.
  • ad hoc to paper The bootstrap regularity criterion (3.1) is assumed to hold on [0,T] before being improved.
    This is the conditional hypothesis in Proposition 3.1. The gap is that at t=0 it is not implied by (1.7), so the assumption is not guaranteed to be available.
  • standard math Sobolev embedding, Gagliardo-Nirenberg interpolation, commutator estimates from Lemma 3.2, and Gronwall's inequality.
    These are standard tools invoked throughout Lemmas 3.3, 3.4, 3.5, 3.6, and 3.7 without proof.

how reviews work

0 comments
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

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 18 canonical work pages

  1. [10]

    Z. Luo, W. Yang, Optimal decay rates for the compressible Navier–Stokes equations with density–dependent viscosities, J. Differ. Equ. 441, 113483 (2025)

  2. [4]

    Y. Guo, Y. Wang, Decay of dissipative equations and negative Sobolev space,Commun. Partial Differ. Equ. 37, 2165-2208 (2012)

  3. [1]

    Adams, J

    A. Adams, J. Fournier, Sobolev Spaces, second edition. Elsevier/Academic Press, Am- sterdam, 2003. 26

  4. [2]

    Bresch, B

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

  5. [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)

  6. [5]

    Global large strong solutions to the compressible Navier-Stokes equations with density-dependent viscosities, case I: isentropic flows

    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

  7. [6]

    Huang, J

    X. Huang, J. Li, Z. Xin, Global well-posedness of classical solutions with large os- cillations and vacuum to the three-dimensional isentropic compressible Navier-Stokes equations, Commun. Pure Appl. Math. 65 (4), 549-585 (2012)

  8. [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)

Show all 20 references
  1. [8]

    D. Hoff, K. Zumbrun, Pointwise decay estimates for multidimensional Navier-Stokes diffusion waves, Z. Angew. Math. Phys. 48, 597-614 (1997)

  2. [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)

  3. [11]

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

  4. [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)

  5. [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)

  6. [14]

    Majda, A.-L

    A. Majda, A.-L. Bertozzi, Vorticity and Incompressible Flow, Cambridge University Press, Cambridge, 2002

  7. [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

  8. [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)

  9. [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)

  10. [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)

  11. [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)

  12. [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

Pith tools

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