Existence of Weak Solutions to a Power-Law Model for Compressible Non-Newtonian Fluids on 1D Unbounded Domain

Jianing Yang; Siran Li; Yuantu Zhu

REVIEW 4 major objections 4 minor 26 references

Existence of Weak Solutions to a Power-Law Model for Compressible Non-Newtonian Fluids on 1D Unbounded Domain

T0 review · 4 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

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

Pith's one-line read The paper proves that as the power-law index p tends to infinity in a one-dimensional compressible non-Newtonian fluid model, the limiting velocity gradient obeys |∂xu|≤1 and the limiting stress acts through a Lagrange multiplier, yielding

desk verdict A serious extension of Bresch-Burtea-Szlenk to the whole line, but two load-bearing gaps in Section 5 leave the main theorem conditional. read the letter →

arxiv 2608.00689 v1 pith:A7I2EFA7 submitted 2026-08-01 math.AP

Siran Li , Jianing Yang , Yuantu Zhu This is my paper · ORCID

classification math.AP MSC 35Q3076A0535D3035K5535B4035B4576N10

keywords compressiblenon-Newtonianfluidspower-lawmodelweaksolutionssingularlimitsaturationconstraintunboundeddomainp-Laplaciandensitybounds

topics Navier-Stokes Regularity

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 establishes that the p→∞ limit of a one-dimensional compressible power-law non-Newtonian fluid model has weak solutions on the entire real line, not just on a periodic domain. In this limit the fluid becomes infinitely shear-thickening, forcing the velocity gradient to saturate at magnitude one and the shear stress to act as a Lagrange multiplier enforcing that constraint. The authors prove that, given finite-energy initial data with a velocity gradient initially below one, there exist density and velocity functions satisfying the continuity and momentum equations in the distributional sense, with density remaining bounded above and below on every compact set. The proof works by truncating the domain to finite intervals, deriving estimates uniform in both the power-law index and the truncation size, and then passing to limits by compactness and diagonalisation. This matters because it extends the physically motivated saturation limit from the torus to unbounded domains, where the two limits p→∞ and domain-size→∞ do not commute.

What carries the argument

The argument revolves around the quantity σ_{p,k}=μ|∂xu_{p,k}|^{p−2}∂xu_{p,k}−aρ_{p,k}^γ, the stress minus the pressure. A parabolic maximum principle for σ_{p,k} yields the pointwise bound that controls the velocity gradient. The main novelty for Dirichlet boundaries is a boundary maximum principle: at the boundary the momentum equation forces a Neumann condition, and Hopf's lemma rules out a boundary maximum, so the estimate is uniform in p and k. Density upper bounds come from a modified potential function that integrates the momentum flux and subtracts a spatial mean; density lower bounds follow from a Lagrangian flow argument. Compactness is provided by the Aubin–Lions–Simon lemma, and

What would settle it

Construct smooth initial data satisfying the theorem's hypotheses and, for a fixed compact set K, numerically track the backward flow map X_t^{-1}(x) for x∈K across increasing truncations k. If the initial points X_t^{-1}(K) drift out of every fixed compact set or the maximal time for which the truncated density stays positive on K shrinks to zero as k→∞, then the k-independent lower bound (7) would fail, refuting Theorem 1.

Watch

Extended reading notes

Core claim

Theorem 1 states that for initial data satisfying local bounds on the density, a velocity gradient strictly below one, and finite total energy, there exist functions (ρ,u) and a Radon measure τ on (0,T)×R such that |∂xu|≤1, τ=π∂xu, π≥0, and π(1−|∂xu|)=0 almost everywhere, and (ρ,u) solves the continuity and momentum equations in the distributional sense. The density is locally bounded away from zero and infinity on every compact set, with constants depending only on the data and the compact set. The construction first solves the truncated problem on Ωk=(−2k−2,2k+2) with homogeneous Dirichlet boundary conditions, obtains estimates independent of p and k, passes p→∞ for fixed k to get a satura

Load-bearing premise

The local density lower bound on a compact set is proved by tracing particles backward along the flow, and it assumes the initial density at the backward point is controlled by the given compact-set bounds even though that point need not lie in the same compact set; without a uniform bound on the velocity or flow map, the constant could depend on the truncation size.

Editorial extensions

If this is right

The saturated p=∞ model is solvable on the whole real line for any finite-energy data satisfying the stated bounds, so periodicity is not essential for the saturation limit.
The constraint |∂xu|≤1 and the complementary-slackness condition π(1−|∂xu|)=0 hold almost everywhere on the unbounded domain, confirming that the Lagrange-multiplier structure is a genuine feature of the limit.
The density remains locally bounded away from zero and infinity on every compact set for all time, with constants depending only on the initial data and the compact set, so no vacuum or unbounded compression is created in finite time.
The convergence is strong for density in C([0,T];L^r_loc(R)) and for velocity in L^2_loc, while the stress converges only as Radon measures, meaning the limiting equations hold in the distributional sense.
The truncation-plus-diagonalisation strategy provides a template for other singular limits on unbounded domains when uniform estimates can be made independent of the truncation parameter.

Reading between the lines

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

The boundary maximum principle and Hopf-lemma argument do not use periodicity, suggesting the same saturation limit should hold on exterior domains or half-lines with suitable boundary conditions.
The energy-comparison step that identifies τ relies on the one-dimensional ordering of the real line; extending the result to R^d would need a different mechanism to select the limit stress, so the 1D unbounded case is likely the natural limit of this method.
The local density lower bound is the fragile point: if the backward flow map can pull initially low density from outside a compact set into it, the k-independence of the lower-bound constant would fail, so tracking Lagrangian trajectories numerically for large truncations would test the proof's key step.
The theorem assumes |∂xu0|<1 strictly; testing initial data with |∂xu0|=1 on a set of positive measure would show whether the saturation limit is stable at the boundary of the admissible class.

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.

Desk Editor's Note

read the letter

Short version: the paper takes the periodic result of Bresch-Burtea-Szlenk and extends it to the real line by domain truncation, with a new boundary maximum principle and a diagonalization argument. The fixed-k analysis in Appendix A is careful and mostly credible; the authors are honest that the framework is inherited. The problem is the passage k to infinity. Two steps in Section 5 do not hold as written.

First, Lemma 9 claims a uniform-in-k lower bound for the density on any compact K. The lower bound in Proposition 5 is pointwise along the Lagrangian flow: it controls rho(t,X_t(x)) in terms of rho_0(x). To get a bound at a point (t,x) in K you need X_t^{-1}(x) to stay in a set where rho_0 has a uniform lower bound. The flow map is not controlled on K, so the constant may depend on k. The same lemma replaces |Omega_k| by |K| in the potential bound for the upper density estimate; that is not justified because the potential is defined by integrals over Omega_k. So the local uniform density bounds are not established.

Second, the identification of tau = pi du in Step 3 uses a local energy comparison. The passage from (33) to (35) assumes that terms supported on d phi_m vanish as the cutoffs concentrate on dK. But tau is only a Radon measure; (31) gives L^1 bounds on tau_k, which do not pass to an L^1 limit. The product tau u may carry singular mass on dK, and the claimed convergence int tau_k u_k dphi -> int tau u dphi is not justified: weak-* convergence of measures plus strong L^2 convergence of u_k is not enough. The displayed energy identity (33) also does not match the usual local conservation law, so the formal comparison in (35)-(36) is not on solid ground.

These are load-bearing gaps, not cosmetics. The final theorem may be true, and the strategy looks repairable - one could try to prove tau in L^1 via better estimates, or use a different identification argument - but as written the proof of Theorem 1 is incomplete.

The paper deserves a serious referee. The area is active, the extension is natural, and the boundary maximum principle is a genuine new tool. I would send it to review, but the referee should be asked to focus on Section 5 and Lemma 9. I would not rely on the theorem in my own work until the gaps are closed.

Referee Report

4 major / 4 minor

Summary. This paper studies the singular limit p→∞ for a one-dimensional compressible power-law fluid (1)–(2) on the whole real line, and claims Theorem 1: from finite-energy initial data with local density bounds, one obtains weak solutions (ρ,u) and a Radon measure τ satisfying the continuity and momentum equations in D′, the saturation constraint |∂xu|≤1, τ=π∂xu, π≥0, and π(1−|∂xu|)=0 a.e., together with local-in-space density bounds. The strategy is to truncate to intervals Ωk, prove estimates uniform in p and k, pass p→∞ for each fixed k (Theorem 8), then pass k→∞ by a diagonal compactness argument. The final identification of τ as a Lagrange multiplier is carried out in Section 5, Step 3, via a local energy comparison.

Significance. If correct, the result would be a nontrivial extension of the periodic result of Bresch–Burtea–Szlenk [4] to the whole line, and would handle the non-commuting limits p→∞ and k→∞ with τ obtained only as a Radon measure. The paper contains substantial and carefully written a priori estimates: the boundary maximum principle for σ, the modified Basov–Shelukhin potential, the L²(0,T;L∞) stress bound, and the detailed appendix computations are internally consistent and constitute a useful technical core. The compactness framework is clearly laid out. However, as detailed below, the proof of the central identification of τ and of the k-uniform density lower bound has load-bearing gaps, and the main theorem as stated is not established.

major comments (4)

[Lemma 9, §5 (Eq. (28))] The k-uniform local lower bound on ρ is not proven. Proposition 5's lower bound (44) is obtained along Lagrangian trajectories: ρ_{p,k}(t,X_t(x)) is estimated from ρ_{p,k,0}(x). For a point x∈K, the preimage X_t^{-1}(x) need not remain in K, and on the whole interval Ωk the initial lower bound c1(Ωk) is not controlled as k→∞, since ρ0 is only assumed bounded below on compact sets. No estimate controlling the Lagrangian flow on a neighborhood of K is available; the bound (30) for ∂tu actually uses (28), so it cannot supply this control. Consequently the constants in (28) and the conclusion (7) are unsupported.
[§5, Step 3, Eq. (33)] The displayed local energy identity is not the standard conservation law for system (6). A direct computation gives d/dt∫φ(1/2ρu²+a/(γ−1)ργ) + μ∫φτ∂xu = ∫∂xφ [u(1/2ρu² + aγ/(γ−1)ργ) − μτu], with u³ and aγ/(γ−1)ργu in the flux. Equation (33) has instead ρu²∂xφ − aργu∂xφ − μτu∂xφ, which has different powers and signs. Since (35) is obtained by passing to the limit in (33), the energy comparison that identifies τ is based on an incorrect identity.
[§5, Step 3, Eqs. (34)–(35)] The passage from (34) to (35) is not justified. From (31) one only has uniform L¹ bounds on τk; the weak-* limit τ as Radon measures need not be absolutely continuous, so the statement “τ∈L¹_loc” is false. Consequently τu is a priori only a measure, and ∫_{supp∂xφm} τu ∂xφm need not vanish as m→∞: a singular part of τ supported on ∂K can give a nonzero contribution even as supp ∂xφm shrinks. Moreover, the convergence ∫τk uk ∂xφ → ∫τu ∂xφ is not a consequence of τk *⇀ τ and uk→u in L²; for example τk = k 1_{(0,1/k)} dx *⇀ δ0 and uk→0 in L² can produce nonzero limits for ∫uk τk. These two failures break the comparison ∫|τ| ≤ ∫τ∂xu and hence the conclusion τ=π∂xu.
[§5, Step 3, Eq. (36)] The energy equality (36) is obtained by testing the distributional momentum equation with φu. At the level of regularity established in the paper, this operation is not justified: u∈L²(0,T;H¹_loc), ∂tu∈L¹(0,T;H⁻²(K)), ρu has limited regularity, and τ is only a Radon measure. In particular, the term involving τ∂x(φu) requires τ to act on an L² function, which is not available for a singular measure. Since (36) is compared with (35) to deduce τ∂xu=|τ|, this is another load-bearing gap.

minor comments (4)

[§5, Lemma 9] The sentence “This proves (25)” should read “This proves (29)”.
[§5, Step 1] The diagonal extraction should be from nested subsequences; as written, k_j^{(j)} need not be a subsequence of the previously selected k_j^{(m)}.
[Eq. (22)] The exponent in |Ωk|^{2−2/p} should be |Ωk|^{2−4/p} for the L⁴ bound in terms of L^p.
[§1] “in contract to” should be “in contrast to”.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the limiting stress relation τ=π∂xu is identified by an a posteriori energy comparison, not assumed or fitted.

full rationale

The paper's derivation chain is not circular in any of the senses enumerated. No parameter is fitted to data; the initial data and bounds are stated assumptions, not extracted from the target conclusion. The paper does not cite its own authors, and its reliance on Bresch–Burtea–Szlenk [4] is an explicit adaptation with the key estimates reproduced in Appendix A, so that citation is independent external support rather than a self-citation chain. For fixed k, the limiting stress τk is obtained as the weak L2 limit of the p-Laplacian term |∂xu_{p,k}|^{p−2}∂xu_{p,k}; the saturation relation τk=|τk|∂xuk is then derived by comparing the approximate energy inequality (9) with the limiting energy equality (25) obtained by testing the already-weak-converged momentum equation with uk. This is a genuine a posteriori comparison, not a definitional identity or a fitted-input rename. On the whole line, the same pattern is repeated with Radon-measure limits, and even though the passage from (33) to (35) contains a technical gap — the assertion that τ∈L1_loc follows from (31) is not justified because a weak-* limit of L1-bounded measures need not be absolutely continuous — that is a correctness/rigor issue, not a circularity. Likewise, the local density lower bound in Lemma 9 may depend on Lagrangian initial points outside K, but this does not make any result equivalent to its own input. The appended AI statement is immaterial to circularity. Overall, the central claim has independent mathematical content and is not forced by definition or by self-citation.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

No free parameters are fitted. The constants a, μ, γ, T, E0 are inputs. The Lagrange multiplier π is a standard multiplier for the constraint, not an invented entity. The main axioms are standard PDE tools plus an external local existence result.

assumptions (3)

domain assumption Local existence of strong solutions to the truncated system (4) for fixed p and k
Borrowed from Kalousek-Macha-Necasova [15]; needed to start the a priori estimates.
standard math Aubin-Lions-Simon compactness lemma
Used repeatedly to extract subsequences.
standard math The parabolic maximum principle and Hopf lemma are valid for the σ equation
Used in Proposition 5 for the boundary maximum principle; requires sufficient regularity of the solution.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Existence of Weak Solutions to a Power-Law Model for Compressible Non-Newtonian Fluids on 1D Unbounded Domain." pith.science (2026). https://pith.science/paper/A7I2EFA7

@misc{pith2026260800689,
  author       = {Pith},
  title        = {Pith review of: Existence of Weak Solutions to a Power-Law Model for Compressible Non-Newtonian Fluids on 1D Unbounded Domain},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/A7I2EFA7}},
  note         = {Machine review of arXiv:2608.00689}
}

abstract

This paper is concerned with the analysis of a one-dimensional power-law model for compressible fluid dynamics on $\mathbb{R}$, in which the shear stress takes the form $\mu |\partial_{x}u|^{p-2}\partial_{x}u$, where $\mu$ is the viscosity coefficient and $u$ is the velocity. We prove that, in the singular limit $p\rightarrow\infty$, the solutions converge to functions $(\rho,u)$ satisfying $|\partial_{x}u|\leq 1$, $\tau = \pi \partial_{x}u$, $\pi \geq 0$, and $\pi (1 - |\partial_{x}u|) = 0$ a.e. on $\mathbb{R}$. Moreover, we rigorously justify the existence of weak solutions to the limiting equation. The convergence as $p \to \infty$ is obtained via domain truncation and compactness arguments, of which the key challenge is to show that the density remains bounded away from zero and infinity on any compact subset. This extends the recent result of Bresch, Burtea, and Szlenk [Nonlinearity 26 (2026), no. 5, Paper No. 055010.] from one-dimensional periodic domain to the whole real line.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

[4]

Bresch, C

D. Bresch, C. Burtea, M. Szlenk, Discontinuous shear-thickening asymptotic for power-law systems related to compressible flows,Nonlinearity26(2026), no. 5, Paper No. 055010

work page 2026
[1]

Basov, V

I. Basov, V. Shelukhin, Generalized solutions to the equations of compressible Bingham flows,ZAMM Z. Angew. Math. Mech.79(1999), 185–192

work page 1999
[2]

Bothe, J

D. Bothe, J. Pruss,Lp-Theory for a class of non-Newtonian fluids,SIAM J. Math. Anal.39(2007), 379–421

work page 2007
[3]

Berselli, L

L.C. Berselli, L. Diening, M. Růžička, Existence of strong solutions for incompressible fluids with shear dependent viscosities,J. Math. Fluid Mech.12(2010), 101–132. 31

work page 2010
[5]

Chhabra, J.F

R.P. Chhabra, J.F. Richardson, Non-Newtonian Flow and Applied Rheology, 2nd edition, Oxford University Press, 2008, 1–55

work page 2008
[6]

Chhabra, Rheology of Complex Fluids, Springer New York, 2010, 3–34

R.P. Chhabra, Rheology of Complex Fluids, Springer New York, 2010, 3–34

work page 2010
[7]

Diening, Lars and M

L. Diening, Lars and M. Růžička, Strong solutions for generalized Newtonian fluids,J. Math. Fluid Mech.7 (2005), 413–450

work page 2005
[8]

L. Fang, A. Zang, Global existence of strong solutions to the Cauchy problem for a one-dimensional compressible non-Newtonian fluid,J. Math. Fluid Mech.25(2023), no. 1, Paper No. 15, 37 pp

work page 2023

Show all 26 references

[9]

Feireisl, Y.S

E. Feireisl, Y.S. Kwon, A. Novotný, On the long-time behavior of dissipative solutions to models of non-Newtonian compressible fluids,Arch. Ration. Mech. Anal.241(2021), 475–495

2021
[10]

B. Guo, P. Zhu, Partial regularity of suitable weak solutions to the system of the incompressible non-Newtonian fluids,J. Differential Equations178(2002), 281–297

2002
[11]

Hoff, Global existence for 1D, compressible, isentropic Navier–Stokes equations with large initial data,Trans

D. Hoff, Global existence for 1D, compressible, isentropic Navier–Stokes equations with large initial data,Trans. Amer. Math. Soc.303(1987), 169–181

1987
[12]

Hoff, Global solutions of the Navier–Stokes equations for multidimensional compressible flow with discontin- uous initial data,J

D. Hoff, Global solutions of the Navier–Stokes equations for multidimensional compressible flow with discontin- uous initial data,J. Differential Equations120(1995), 215–254

1995
[13]

O. A. Ladyzhenskaya, New equations for the description of the motions of viscous incompressible fluids and global solvability for their boundary value problems,Trudy Mat. Inst. Steklov.102(1967), 85–104

1967
[14]

Jiang, Large-time behavior of solutions to the equations of a one-dimensional viscous polytropic ideal gas in unbounded domains,Comm

S. Jiang, Large-time behavior of solutions to the equations of a one-dimensional viscous polytropic ideal gas in unbounded domains,Comm. Math. Phys.200(1999), no. 1, 181–193

1999
[15]

Kalousek, V

M. Kalousek, V. Mácha, Š. Nečasová, Local-in-time existence of strong solutions to a class of the compressible non-Newtonian Navier–Stokes equations,Math. Ann.384(2022), no. 3-4, 1057–1089

2022
[16]

S. Lian, W. Gao, H. Yuan, C. Cao, Existence of solutions to an initial Dirichlet problem of evolutionalp(x)- Laplace equations,Ann. Inst. H. Poincaré C Anal. Non Linéaire29(2012), 377–399

2012
[17]

Málek, J

J. Málek, J. Nečas, M. Růžička, On the non-Newtonian incompressible fluids,Math. Models Methods Appl. Sci. 3(1993), 35–63

1993
[18]

Málek, J

J. Málek, J. Nečas, M. Růžička, On weak solutions to a class of non-Newtonian incompressible fluids in bounded three-dimensional domains: the casep≥2,Adv. Differential Equations6(2001), 257–302

2001
[19]

Moscariello, M.M

G. Moscariello, M.M. Porzio, On the behavior in time of solutions to motion of non-Newtonian fluids,NoDEA Nonlinear Differential Equations Appl.27(2020), no. 4, Paper No. 42, 18 pp

2020
[20]

Muhammad, L

J. Muhammad, L. Fang, Z. Guo, Global weak solutions to a class of compressible non-Newtonian fluids with vacuum,Math. Methods Appl. Sci.43(2020), 5234–49

2020
[21]

Simon, Compact sets in the spaceLp(0, T;B),Ann

J. Simon, Compact sets in the spaceLp(0, T;B),Ann. Mat. Pura Appl. (4)146(1987), 65–96

1987
[22]

Wolf, Existence of weak solutions to the equations of non-stationary motion of non-Newtonian fluids with shear rate dependent viscosity,J

J. Wolf, Existence of weak solutions to the equations of non-stationary motion of non-Newtonian fluids with shear rate dependent viscosity,J. Math. Fluid Mech.9(2007), 104–138

2007
[23]

H. Yuan, H. Li, J. Qiao, and F. Li, Global existence of strong solutions of Navier–Stokes equations with non- Newtonian potential for one-dimensional isentropic compressible fluids,Acta Math. Sci. Ser. B (Engl. Ed.)32 (2012), 1467–1486

2012
[24]

H. Yuan, X. Si, Z. Feng, Global strong solutions of a class of non-Newtonian fluids with small initial energy,J. Math. Anal. Appl.474(2019), 72–93

2019
[25]

Zhao, On the Cauchy problem and initial traces for the evolutionp-laplacian equations with strongly nonlinear sources,J

J. Zhao, On the Cauchy problem and initial traces for the evolutionp-laplacian equations with strongly nonlinear sources,J. Differential Equations121(1995), 329–383

1995
[26]

J. Zhao, P. Lei, On the Cauchy problem of evolutionp-Laplacian equations with strongly non-linear sources when 1< p <2,Acta Math. Sin. (Engl. Ser.)17(2001), 455–470. 32 Siran Li: School of Mathematical Sciences&CMA-Shanghai, Shanghai Jiao Tong University, No. 6 Science Buildin...

2001

Pith tools

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

[1] [4]

Bresch, C

D. Bresch, C. Burtea, M. Szlenk, Discontinuous shear-thickening asymptotic for power-law systems related to compressible flows,Nonlinearity26(2026), no. 5, Paper No. 055010

work page 2026

[2] [1]

Basov, V

I. Basov, V. Shelukhin, Generalized solutions to the equations of compressible Bingham flows,ZAMM Z. Angew. Math. Mech.79(1999), 185–192

work page 1999

[3] [2]

Bothe, J

D. Bothe, J. Pruss,Lp-Theory for a class of non-Newtonian fluids,SIAM J. Math. Anal.39(2007), 379–421

work page 2007

[4] [3]

Berselli, L

L.C. Berselli, L. Diening, M. Růžička, Existence of strong solutions for incompressible fluids with shear dependent viscosities,J. Math. Fluid Mech.12(2010), 101–132. 31

work page 2010

[5] [5]

Chhabra, J.F

R.P. Chhabra, J.F. Richardson, Non-Newtonian Flow and Applied Rheology, 2nd edition, Oxford University Press, 2008, 1–55

work page 2008

[6] [6]

Chhabra, Rheology of Complex Fluids, Springer New York, 2010, 3–34

R.P. Chhabra, Rheology of Complex Fluids, Springer New York, 2010, 3–34

work page 2010

[7] [7]

Diening, Lars and M

L. Diening, Lars and M. Růžička, Strong solutions for generalized Newtonian fluids,J. Math. Fluid Mech.7 (2005), 413–450

work page 2005

[8] [8]

L. Fang, A. Zang, Global existence of strong solutions to the Cauchy problem for a one-dimensional compressible non-Newtonian fluid,J. Math. Fluid Mech.25(2023), no. 1, Paper No. 15, 37 pp

work page 2023