{"id":"67a73241-f655-4194-9c67-3a2e8f63d36e","arxiv_id":"2507.15410","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Letting the power-law exponent tend to infinity in compressible fluid equations yields a limit system with the unilateral constraint |shear rate| <= 1; proven fully in 1D, conditionally in multiple dimensions.","lead":"This mathematics paper studies what happens to the equations for compressible power-law fluids when the power-law exponent is sent to infinity. The authors prove that, in several regimes, the limit is a 'thick fluid' model in which the shear rate is capped by a hard constraint, extending earlier incompressible results.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1D fixed-p global-existence step is cited, not proved; without it Theorems 1.1 and 1.5 lack a solution sequence to pass to the p→∞ limit.","rationale":"The reader's weakest_assumption identifies exactly this gap: the fixed-p global weak solutions with uniform bounds are required for the compactness passages, and in the 1D case the upgrade from local to global existence is skipped. I agree that this is the most load-bearing condition. The other concerns—the stress maximum-principle sign, the conditional multi-D theorem, and the Section 4 sketch—are real but less central: the stress bound is plausibly fixable with a sign-aware maximum-principle argument; the multi-D result is explicitly conditional and Theorem 1.10 offers a regularized construction; Section 4 is declared a sketch. In contrast, the 1D fixed-p global existence underlies the main unconditional claim and no substitute construction is provided. My concern therefore reinforces the CONDITIONAL verdict rather than changing it; I recommend keeping the reader's verdict as CONDITIONAL, so the verdict adjustment is UNCHANGED.","tokens_in":24811,"tokens_out":15712,"duration_ms":182494,"concrete_test":"Attempt the missing continuation argument for fixed p: start from the local strong solution supplied by [24] and prove that the bounds of Propositions 2.1–2.4 (uniform density lower/upper bounds, ∫ρ|dot u|^2≤C(T), and ∫|∂_x u|^p≤E0) prevent finite-time blow-up and yield a global weak solution for each p≥2. Concretely, derive a blow-up criterion: if T*<T is the maximal existence time, identify the quantity whose control would allow restarting the local solution after an interval whose length depends only on E0,c1,c2,T; check whether the existing bounds control that quantity. If no such continuation theorem can be produced or located, Theorem 1.1 must be restated as conditional and Theorem 1.5 remains unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 starts by invoking [24] for local-in-time strong solutions, says 'we skip the details', and then states Theorem 1.1: global weak solutions for fixed p≥2 with uniform-in-p bounds. Theorem 1.5 then needs a subsequence of such global solutions to pass to the limit. Since [24] is local-in-time, no continuation or global-existence argument is supplied. Propositions 2.1 and 2.4 give a priori estimates on smooth solutions—density away from vacuum, ∫ρ|dot u|^2≤C(T), and ∫|∂_x u|^p≤E0—but a priori estimates do not by themselves construct a global weak solution on [0,T]. If fixed-p global solutions with these bounds do not exist, the convergence statements in (1.6), Proposition 2.6, and the strong-form identification in Section 2.2.3 have no sequence to act on, and Theorem 1.5 collapses. This is the load-bearing condition of the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the limit p -> infinity for several compressible power-law fluid systems. In the one-dimensional non-stationary case, the authors claim that, under initial data satisfying |partial_x u_{p,0}| <= 1, the compressible power-law system (1.3) converges to a 'thick' compressible system with the constraint |partial_x u| <= 1, the stress relation tau = pi partial_x u, and the complementarity condition pi(1 - |partial_x u|) = 0. In the multi-dimensional semi-stationary Stokes case, they establish a conditional compactness result toward a variational inequality, and they construct solutions to a regularized variant of the limit problem. The proofs combine a stress maximum principle, uniform density bounds, Hoff-type estimates, and monotonicity arguments.","tokens_in":24935,"tokens_out":12786,"duration_ms":134316,"significance":"If made fully rigorous, the 1D result is a natural and nontrivial extension of the known p -> infinity thick-fluid limits to the compressible setting, with motivation from shear-thickening suspensions in which large pressures may require compressibility. The paper is honest about the multi-dimensional theorem being conditional, and the appendix contains a complete proof of a useful continuity property. No free parameters or fitted data enter the analysis; the limit systems are derived from a priori estimates. The compactness chain is plausible and contains interesting ideas. However, the missing fixed-p global existence step is load-bearing for the main theorem, and the ordering of one compactness argument appears circular as written.","major_comments":[{"comment":"Theorem 1.1 asserts global weak solutions for fixed p (large enough) with uniform-in-p bounds, and this is the starting point for the entire limit passage. The proof only cites [24] for local-in-time strong solutions and then states 'we skip the details'; the a priori estimates in Propositions 2.1 and 2.4 concern smooth solutions and do not by themselves produce a global weak solution on [0,T]. A continuation argument, or a reference to a global existence theorem for (1.3)-(1.5), is missing. Since Theorem 1.5 and the convergence statements in (1.6) require a sequence of such global solutions, this gap is load-bearing.","section":"Section 2, Theorem 1.1"},{"comment":"The strong density convergence argument uses the bound |partial_x u| <= 1 ('Note that we know that |partial_x u| <= 1...') before Proposition 2.6 in Section 2.2.2 establishes this bound. As written the proof is circular: the density compactness step appears to rely on the very constraint that is only derived afterward. Please reorder the proof so that the L^infty bound on partial_x u is proved first, or explain which terms in Section 2.2.1 can be handled without it.","section":"Section 2.2.1"},{"comment":"The notion of solution used in Theorem 1.9 is not aligned with Definition 1.8. Definition 1.8 defines variational solutions through the energy inequality, the continuity equation, and a variational momentum inequality, but the proof of Theorem 1.9 tests the PDE (1.8) by psi u_p and psi u, which requires the PDE formulation. Please state precisely which assumptions are needed. In addition, Theorem 1.10 relies on the sentence 'from the proof it is clear' that the global existence result of [18] applies to the modified system (3.8), with the deviatoric part removed and the singular term changed. This is a nontrivial transfer of a cited theorem and is the basis for the approximating sequence; it needs a verification or a direct construction.","section":"Section 3, Theorems 1.9 and 1.10"}],"minor_comments":[{"comment":"The convergence 'partial_t u_p -> partial_t u weakly-* in L^2(0,T;L^infty(T))' appears to be a typo; the available estimate on partial_t u_p is in L^2((0,T)xT), so the convergence should be stated in L^2((0,T)xT) unless an additional L^infty bound is proved.","section":"Theorem 1.1, display (1.6)"},{"comment":"The level-set proof bounds the energy by the entropy expression integral(rho_0 log rho_0 - rho_0), which is not the quantity supplied by the energy inequality; the available uniform bound is E_0. Since any fixed finite constant would suffice for the contradiction, this is not fatal, but the displayed estimate should be corrected.","section":"Proposition 2.6"},{"comment":"The variational inequality satisfied by global weak solutions at fixed p is asserted without derivation. A short derivation from the weak formulation and the convexity of the stress potential would improve the readability and would clarify the role of the constraint |partial_x v| <= 1.","section":"Section 2.2.2"},{"comment":"Remark 1.6 ends with the sentence 'Note also that the formulation (1.7) with the regularity found on (rho,u) given in Theorem 1.1 is equivalent to the problem' and then stops; the sentence is incomplete.","section":"Remark 1.6"},{"comment":"There are several minor typographical and grammatical issues, for example 'stated with the work' should be 'starting with the work' in the abstract, and 'Concerinig' should be 'Concerning' in the introduction; these should be corrected in a revised version.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript contains a good idea and the eventual theorems would be a strong fit for the journal. My main concern is that the missing fixed-p global existence in the 1D section is an honest gap rather than a purely presentational issue. If the authors supply a complete continuation or compactness argument, or clearly state the 1D result as conditional, I would be willing to see a revised version. The multi-dimensional existence is already honestly labeled as open, so I would not require that gap to be filled."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a real extension of the thick-fluid p→∞ theory to the compressible setting, and the 1D part is mostly sound. The main gap is the fixed-p global existence step in Section 2, which is cited rather than proved; that is fixable but it is load-bearing for the limit. The multi-D results are more conditional than the 1D theorem, and Section 4 is a sketch.\n\nWhat's new: the 1D nonstationary theorem gives the expected limit system with |∂xu|≤1 and the complementarity relation, and it goes through a coherent compactness chain: stress maximum principle, density bounds, Hoff-type estimate on ρ u̇, strong convergence of density, and a monotonicity argument to identify the stress. The maximum principle on σp = μ|∂xu|^{p−2}∂xu − aρ^γ is an elegant piece; the sign term the reader flagged is compressed in presentation but is actually non-positive by monotonicity of H^{−1}, so the argument survives scrutiny. The multidimensional semi-stationary compactness argument yields the |Du|≤1 constraint and strong density convergence under an energy bound, and the paper is honest that existence of the approximating sequence is open. The singular-stress extension in Section 4 is plausible but only sketched.\n\nSoft spots, in order of importance. First, Section 2 says 'we skip the details' for local existence and then states Theorem 1.1: global weak solutions for fixed p with uniform-in-p bounds. The cited [24] is local-in-time, and no continuation argument is supplied. The a priori estimates in Propositions 2.1 and 2.4 are exactly what a continuation would need, but the step is not written. This matters because Theorem 1.5 needs a sequence of such global solutions; if the step fails, the limit statement has nothing to act on. I think it is likely fixable by standard continuation, but the paper should say so. Second, Theorem 1.10 leans on an unproved adaptation of [18] (dropping the convective term, replacing D_d u by Du). A referee should ask for that proof or a precise statement. Third, Theorem 1.9 is explicitly conditional, which is fine but limits its status. Section 4 should either be expanded or labeled as formal. The citation pattern is fine; the self-citations are to prior published compactness lemmas and singular-limit results, not to prop up this argument.\n\nWho this is for: people working on p→∞ limits, thick fluids, compressible non-Newtonian flows. The 1D result is likely to be cited. My recommendation: send it to referees, with the expectation of heavy revision rather than acceptance as is.","headline":"A genuinely new p→∞ thick-fluid limit for 1D compressible power-law flows, but the fixed-p global existence step is cited rather than proved and the multi-D results lean on conditional or sketched arguments; deserves refereeing.","tokens_in":25528,"tokens_out":6069,"would_cite":true,"duration_ms":61478,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q35","76A05","76N10","35B25","49J40"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that the p→∞ limit of compressible power-law shear-thickening fluids is a thick fluid whose shear rate is capped at a maximum value.","keywords":["discontinuous shear thickening","power-law fluids","compressible Navier-Stokes","thick fluids","maximum admissible shear rate","p-Laplacian asymptotics","variational inequalities","unilateral constraint"],"falsifier":"Run a high-p direct numerical simulation of the periodic one-dimensional system (1.3) with data satisfying |∂x u_{p,0}| ≤ 1: if the density develops a finite-time blow-up while the energy stays bounded, the uniform estimates of Proposition 2.1 are false and the convergence to (1.7) cannot hold; the paper's claim conversely predicts that such simulations stay bounded in density and converge to the constrained system with |∂x u| ≤ 1.","tokens_in":24553,"feed_emoji":"🌊","tokens_out":10352,"duration_ms":101260,"temperature":0.7,"pith_summary":"The paper proves that the power-law model for compressible shear-thickening fluids, in which the stress grows like |D u|^{p−2} D u, has a well-defined limit as the exponent p tends to infinity, and that the limit is a 'thick' compressible fluid whose shear rate is capped at a maximum value. In the one-dimensional non-stationary case, the paper shows that any sequence of weak solutions with initial data satisfying |∂x u_{p,0}| ≤ 1 converges to a solution of a PDE in which the velocity gradient is constrained by |∂x u| ≤ 1 and the viscous stress takes the form τ = π ∂x u with a Lagrange multiplier π ≥ 0 that vanishes unless |∂x u| = 1. The same limiting mechanism is obtained for the multi-dimensional semi-stationary Stokes system, where the limit is described by a variational inequality with a constraint on the symmetric strain tensor, provided the approximating fixed-p solutions exist with the stated uniform bounds. The paper also shows that the same limit arises from a singular shear-rate-dependent stress law in one dimension. A reader should care because this gives a mathematical justification, in the compressible setting, for the unilateral strain constraint that models discontinuous shear thickening in suspensions such as cornstarch.","feed_headline":"Compressible power-law flow gains a shear-rate cap as p→∞","feed_subtitle":"The limit caps the shear rate and turns the stress into a Lagrange multiplier, a direct PDE model for thick fluids.","key_machinery":"The argument is carried by p-uniform a priori estimates built on two devices. The first is a maximum-principle bound for the generalized Cauchy stress σ_p = µ|∂x u_p|^{p−2}∂x u_p − a ρ_p^γ, which controls the positive part of the strain rate and, via a transport argument, yields the lower density bound; the upper density bound uses the Basov–Shelukhin identity on ρ_p^µ exp(−ψ_p). The second device is the monotonicity of the p-Laplacian-type stress tensor S_p(Du)=|Du|^{p−2}Du, which gives strong convergence of the density through a Grönwall-type inequality on the convexity gap X_p=∫(ρ_p^γ−ρ^γ−$γρ^{{γ−1}}$(ρ_p−ρ)), and then forces the limit relation |τ|=τ∂x u, from which τ=π∂x u with π≥0 and π(1−|∂x u|)=0 follows.","core_discovery":"The central claim is Theorem 1.5: for the one-dimensional periodic compressible power-law system (1.3), initial data with bounded energy, density bounded away from zero and infinity, and |∂x u_{p,0}| ≤ 1, produce a subsequence converging strongly enough to pass to the limit in the equations, and the limit (ρ,u) satisfies the continuity equation, the momentum equation with τ = |∂x u|^{p−2}∂x u replaced by τ = π ∂x u, the inequality |∂x u| ≤ 1, and the complementarity condition π(1 − |∂x u|) = 0. In other words, the p→∞ limit of the compressible power-law system is exactly the thick compressible fluid with maximum admissible shear rate. In the multi-dimensional semi-stationary Stokes case (Theorem 1.9), the paper reaches the same constraint |D u| ≤ 1 and characterizes the limit by the variational inequality (1.12), assuming the approximating fixed-p weak solutions exist; the construction is made unconditional in Theorem 1.10 by adding a regularizing term for div u and using the singular-stress existence theory. Section 4 shows that the one-dimensional singular stress law ∂x u/√(1 − |∂x u|²) leads to the same limiting system.","pith_inferences":["If global weak solutions of the multi-dimensional semi-stationary system (1.8) with the uniform energy bound can be constructed, Theorem 1.9 would upgrade from a conditional compactness result to a fully constructive convergence theorem; the one-dimensional proof suggests the missing ingredient is a fixed-p existence theory with bounds independent of p.","The complementarity condition π(1−|∂x u|)=0 is structurally the same as the 'free/congested' density constraint π(1−ρ)=0 studied elsewhere, and the two regimes might be combined in a single two-constraint model for jamming flows.","The one-dimensional mechanism, based on maximum principles and transport estimates, is likely to extend to bounded intervals with no-slip boundary conditions if the Poincaré–Wirtinger step is replaced by a suitable estimate on the mean velocity, which would make the model applicable to squeezing-flow experiments.","A testable quantitative prediction is that in the limiting model the shear stress saturates at the maximal shear rate and becomes independent of further increases in applied pressure; rheometric experiments on dense suspensions could look for such a plateau."],"forward_implications":["The p→∞ limit of the compressible power-law model is a well-defined PDE with a unilateral constraint, so the maximum-admissible-shear-rate condition can be derived from the power-law model rather than assumed.","In the limit, the fluid cannot shear faster than the threshold |∂x u| = 1 (or |D u| ≤ 1 in multiple dimensions), and the stress is a Lagrange multiplier that is active only at the threshold, giving a precise rheological characterization of the thickened state.","The same limit obtained from the singular stress law in Section 4 indicates that the asymptotic constrained model is insensitive to whether the viscosity diverges through a power law or through a singular square-root-type term.","The variational-inequality formulation in multiple dimensions provides a framework for studying existence, uniqueness, and numerical approximation of the thick compressible fluid without needing strong solutions."],"supporting_citations":[{"why":"Supplies the local-in-time existence of strong solutions for the fixed-p system that the paper refers to for the existence step underlying the approximation.","marker":"[24]"},{"why":"Provides the Basov–Shelukhin identity used to control the density from above uniformly in p, a key ingredient in the uniform estimates.","marker":"[3]"},{"why":"Introduces the material-time-derivative estimate used in Proposition 2.4 to bound ρ_p|u_p dot|² uniformly and thereby control the stress.","marker":"[21, 22]"},{"why":"Supplies the method for passing from the variational limit to the strong form τ=π∂x u; the paper's proof of Theorem 1.5 is inspired by this result.","marker":"[14]"},{"why":"Gives the variational and quasi-variational inequality framework with gradient constraints used to characterize the limit in the non-stationary and semi-stationary systems.","marker":"[35, 36]"},{"why":"Provides global existence of variational solutions for the non-Newtonian compressible system with singular stress, used in the proof of Theorem 1.10 to construct solutions of the multi-dimensional limit problem.","marker":"[18]"}],"fun_headline_variants":["Power-law limit caps shear rate in compressible flows","Shear-thickening limit yields thick fluid constraint","Compressible power-law converges to maximum shear-rate cap","Discontinuous thickening limit enforces shear-rate cap in fluids"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that global weak solutions of the fixed-p compressible power-law systems exist with the uniform-in-p bounds stated in Theorem 1.1; the paper cites a local-in-time existence result and does not supply the full upgrade to global weak solutions, and for the multi-dimensional semi-stationary system existence is explicitly left open, so if those fixed-p solutions or their uniform bounds fail, the compactness passages do not run.","fun_headline_variants_meta":{"raw":{"variants":["Power-law limit caps shear rate in compressible flows","Shear-thickening limit yields thick fluid constraint","Compressible power-law converges to maximum shear-rate cap","Discontinuous thickening limit enforces shear-rate cap in fluids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00054,"raw_usage":{"total_tokens":2646,"prompt_tokens":1056,"completion_tokens":1590,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":672,"completion_tokens_details":{"reasoning_tokens":1525}},"tokens_in":672,"tokens_out":1590,"duration_ms":12052,"temperature":1.0,"reasoning_tokens":1525,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:34:32.369160+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a high-p direct numerical simulation of the periodic one-dimensional system (1.3) with data satisfying |∂x u_{p,0}| ≤ 1: if the density develops a finite-time blow-up while the energy stays bounded, the uniform estimates of Proposition 2.1 are false and the convergence to (1.7) cannot hold; the paper's claim conversely predicts that such simulations stay bounded in density and converge to the constrained system with |∂x u| ≤ 1.","supporting_citations":[{"cited_title":"Kalousek, V","cited_arxiv_id":null,"evidence_quote":"Supplies the local-in-time existence of strong solutions for the fixed-p system that the paper refers to for the existence step underlying the approximation."},{"cited_title":"Basov and V","cited_arxiv_id":null,"evidence_quote":"Provides the Basov–Shelukhin identity used to control the density from above uniformly in p, a key ingredient in the uniform estimates."},{"cited_title":"De Pascale, L","cited_arxiv_id":null,"evidence_quote":"Supplies the method for passing from the variational limit to the strong form τ=π∂x u; the paper's proof of Theorem 1.5 is inspired by this result."},{"cited_title":"Feireisl, X","cited_arxiv_id":null,"evidence_quote":"Provides global existence of variational solutions for the non-Newtonian compressible system with singular stress, used in the proof of Theorem 1.10 to construct solutions of the multi-dimensional limit problem."}],"review_version":1}