{"id":"e88b04e8-d70c-4830-a357-7124e25bf1d0","arxiv_id":"2608.00689","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Existence of weak solutions to the p=∞ limiting system of 1D compressible power-law fluids is established on ℝ by passing p→∞ on truncated intervals and then letting the intervals grow.","lead":"This paper proves that a shear-thickening power-law fluid model, as the power-law exponent tends to infinity, has weak solutions on the whole real line, with the velocity gradient limited in magnitude to at most 1 and a stress playing the role of a Lagrange multiplier for that constraint. The result extends a recent periodic-domain theorem to the unbounded line via domain truncation and compactness.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Section 5 Step 3: local energy comparison used to identify τ=π∂xu on R is not justified — τ is only a Radon measure, so the boundary terms in (33)→(35) need not vanish, and the displayed local identity is not the standard conservation law.","rationale":"The reader's verdict is CONDITIONAL, and the rationale explicitly cites the boundary terms in the local energy comparison when the stress is a measure as a structural gap. The single most load-bearing issue is therefore this local energy comparison in Section 5, Step 3: it is the only argument identifying the limiting stress as a Lagrange multiplier on R. The density-lower-bound gap in Lemma 9 is also real and affects the regularity claim (7), but the measure/boundary-term issue directly undermines the central identity τ=π∂xu, without which the theorem's main conclusion fails even if density bounds are repaired. Since the reader already flagged this issue and chose CONDITIONAL, the present stress-test does not change the verdict; it sharpens the concern by noting that the local energy identity (33) itself appears inconsistent with the standard conservation-law calculation, in addition to the measure-convergence problem.","tokens_in":28557,"tokens_out":18013,"duration_ms":217834,"concrete_test":"Independently re-derive the local energy identity for a smooth solution with cutoff φ and compare it with (33); if the correct identity contains the u^3 and aγ/(γ−1)ρ^γu flux terms, the displayed inequality (33) is incorrect. Then test the measure limit: take τ_k=k1_{(0,1/k)} so τ_k⇀δ_0 as measures, and g_k→0 in L^2 but not uniformly; show ∫τ_k g_k does not tend to 0. Apply this with g_k=u_k∂xφ_m on supp∂xφ_m→∂K to demonstrate that the vanishing boundary-term claim requires uniform convergence of u_k∂xφ_m or absolute continuity of τ, neither of which is established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The identification of the limiting stress as τ=π∂xu on the whole line is carried out in Section 5, Step 3, via a local energy inequality. Two load-bearing transitions are not justified. First, the local energy identity (33) does not match the standard local conservation law: for smooth solutions one obtains d/dt∫φ(1/2ρu^2+a/(γ−1)ρ^γ)+μ∫φτ∂xu = −∫∂xφ[u(1/2ρu^2 + aγ/(γ−1)ρ^γ) − μuτ], whose flux contains u^3 and aγ/(γ−1)ρ^γu; the right-hand side in (33) (ρu^2∂xφ − aρ^γu∂xφ − μτu∂xφ) has different powers and signs. Second, the passage from (33) to (35) sends supp ∂xφ_m → ∂K and asserts the boundary terms vanish using 'τ∈L^1_loc'. But (31) only gives uniform L^1 bounds on τ_k; the weak-* limit is a Radon measure and need not be absolutely continuous. The product τu can have singular mass on ∂K, so ∫τu∂xφ_m need not tend to 0. Moreover, with only weak-* convergence τ_k⇀τ as measures, strong L^2 convergence of u_k is insufficient to justify ∫τ_k u_k∂xφ → ∫τu∂xφ; a counterexample such as τ_k=k1_{(0,1/k)} shows the product integral can fail to converge under merely L^2-convergent multipliers. Without (35), the comparison ∫|τ|≤∫τ∂xu cannot be made, and τ=π∂xu is unproved.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":28987,"tokens_out":15928,"duration_ms":197739,"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":[{"comment":"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.","section":"Lemma 9, §5 (Eq. (28))"},{"comment":"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.","section":"§5, Step 3, Eq. (33)"},{"comment":"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.","section":"§5, Step 3, Eqs. (34)–(35)"},{"comment":"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.","section":"§5, Step 3, Eq. (36)"}],"minor_comments":[{"comment":"The sentence “This proves (25)” should read “This proves (29)”.","section":"§5, Lemma 9"},{"comment":"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)}.","section":"§5, Step 1"},{"comment":"The exponent in |Ωk|^{2−2/p} should be |Ωk|^{2−4/p} for the L⁴ bound in terms of L^p.","section":"Eq. (22)"},{"comment":"“in contract to” should be “in contrast to”.","section":"§1"}],"recommendation":"reject","confidential_remarks":"The paper contains genuine technical content in the a priori estimates, but the main theorem is not established. The global identification of the Lagrange multiplier is the heart of the extension to R, and the measure-theoretic issues there are not cosmetic. I do not see how to fix them within the current framework without either additional assumptions (e.g., a global-in-space lower bound on ρ0 or uniform equi-integrability of τk) or a substantially different argument. I therefore recommend rejection, though the authors may wish to resubmit a revised version with the local energy argument replaced and the density lower bound reworked."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nFirst, 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.\n\nSecond, 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.\n\nThese 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.\n\nThe 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.","headline":"A serious extension of Bresch-Burtea-Szlenk to the whole line, but two load-bearing gaps in Section 5 leave the main theorem conditional.","tokens_in":29430,"tokens_out":5958,"would_cite":true,"duration_ms":74488,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q30","76A05","35D30","35K55","35B40","35B45","76N10"],"pacs":[],"model":"deepseek-v4-flash","headline":"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","keywords":["compressible non-Newtonian fluids","power-law model","weak solutions","singular limit","saturation constraint","unbounded domain","p-Laplacian","density bounds"],"falsifier":"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.","tokens_in":28478,"feed_emoji":"🌊","tokens_out":4877,"duration_ms":58783,"temperature":0.7,"pith_summary":"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.","feed_headline":"Power-law fluid limit solved on the whole real line","feed_subtitle":"Saturation constraint |∂xu|≤1 and Lagrange-multiplier stress survive without periodic boundaries.","key_machinery":"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","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the periodic-domain result and the σ maximum-principle framework that this paper adapts to Dirichlet boundaries and unbounded domains.","marker":"[4]"},{"why":"Introduces the potential function used to bound the density from above in the truncated problem.","marker":"[1]"},{"why":"Provides the Aubin–Lions–Simon compactness lemma used to pass to the limit in both p and k.","marker":"[21]"},{"why":"Gives the local-in-time existence of strong solutions used to start the approximate-solution construction on each truncation.","marker":"[15]"},{"why":"Supplies ideas from one-dimensional compressible Navier–Stokes analysis used for the time-derivative estimates.","marker":"[11]"},{"why":"Supplies additional ideas for estimating time derivatives of the velocity, adapted here to the power-law context.","marker":"[12]"}],"fun_headline_variants":["Power-law fluid existence extended to whole real line","Weak solutions on R for compressible non-Newtonian fluids","Saturation limit p→∞ solved for power-law fluids on R","Whole-line weak solutions for power-law fluids with saturation"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Power-law fluid existence extended to whole real line","Weak solutions on R for compressible non-Newtonian fluids","Saturation limit p→∞ solved for power-law fluids on R","Whole-line weak solutions for power-law fluids with saturation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000362,"raw_usage":{"total_tokens":1817,"prompt_tokens":795,"completion_tokens":1022,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":965}},"tokens_in":539,"tokens_out":1022,"duration_ms":10474,"temperature":1.0,"reasoning_tokens":965,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T00:33:41.632553+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Bresch, C","cited_arxiv_id":null,"evidence_quote":"Supplies the periodic-domain result and the σ maximum-principle framework that this paper adapts to Dirichlet boundaries and unbounded domains."},{"cited_title":"Basov, V","cited_arxiv_id":null,"evidence_quote":"Introduces the potential function used to bound the density from above in the truncated problem."},{"cited_title":"Simon, Compact sets in the spaceLp(0, T;B),Ann","cited_arxiv_id":null,"evidence_quote":"Provides the Aubin–Lions–Simon compactness lemma used to pass to the limit in both p and k."},{"cited_title":"Kalousek, V","cited_arxiv_id":null,"evidence_quote":"Gives the local-in-time existence of strong solutions used to start the approximate-solution construction on each truncation."},{"cited_title":"Hoff, Global existence for 1D, compressible, isentropic Navier–Stokes equations with large initial data,Trans","cited_arxiv_id":null,"evidence_quote":"Supplies ideas from one-dimensional compressible Navier–Stokes analysis used for the time-derivative estimates."},{"cited_title":"Hoff, Global solutions of the Navier–Stokes equations for multidimensional compressible flow with discontin- uous initial data,J","cited_arxiv_id":null,"evidence_quote":"Supplies additional ideas for estimating time derivatives of the velocity, adapted here to the power-law context."}],"review_version":1}