{"id":"16932619-d36a-4301-881f-bff0d8785449","arxiv_id":"2505.17588","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"The paper proves existence and uniqueness of weak solutions for a simplified granular flow model with threshold rheology and dilatancy.","lead":"Mathematicians prove that a simplified model of dense granular flow, combining a pressure-dependent yield threshold with volume-change (dilatancy), admits a weak solution with unique velocity and pressure. The proof uses a new energy-based weak formulation that recovers the threshold rheology from an energy inequality, a step toward well-posedness for realistic granular models.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 1 relies on Proposition 2, which asserts existence for the ε-approximate system, but the proof is only an energy estimate plus a reference to [2]; the Galerkin passage for the coupled nonlinear system is not demonstrated, leaving the central limit argument without a proven starting point.","rationale":"Read in good faith: the mathematical structure of the proof is coherent and the L3 framework is internally consistent. The convex-analysis lemmas and the identification steps in §3.6 appear sound, and Proposition 1 correctly recovers the rheology from the weak formulation. The reader's main concern about the artificial viscosity is explicitly acknowledged by the authors in Remark 3; the paper does not claim Theorem 1 for the inviscid model, so that point concerns scope and physical relevance rather than correctness of the stated theorem. The most concrete, load-bearing issue is Proposition 2: the approximate system is the starting point of the entire limiting argument, and its existence proof is deferred to a standard-reference claim without the needed details. This does not show the theorem is false; it shows the manuscript is incomplete as written and justifies a conditional verdict. I therefore uphold the reader's CONDITIONAL verdict.","tokens_in":20674,"tokens_out":30876,"duration_ms":229725,"concrete_test":"Write out a complete Galerkin proof for Proposition 2: take finite-dimensional subspaces for (u,p), obtain the Galerkin ODE system, prove the Galerkin analogue of the energy estimate (24), derive bounds on ∂tu_N and ε∂tp_N uniform in the Galerkin dimension N, and pass N→∞ in (22)-(23), explicitly identifying the weak limit of ∫ p_N Su_N/(|Su_N|+ε) : Sϕ. If the identification requires strong convergence of Su_N or p_N beyond what Aubin-Lions gives from the stated bounds, Proposition 2 fails and Theorem 1 is not established; if the passage closes, the gap is expository and the theorem stands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proof of Theorem 1 proceeds by constructing approximate solutions (17)-(20) for each ε>0. Proposition 2 claims existence and uniqueness of these solutions, but its 'proof' consists of stating that a Galerkin method applies 'see [2]' and then deriving the energy identity (24). This is load-bearing: without a rigorously constructed sequence (uε,pε,σε), the uniform estimates in Proposition 3, the compactness in Proposition 4, and the identification of weak limits in §3.6 have no object to act on. The missing part is not cosmetic: for fixed ε the system couples a p-Laplacian-type momentum equation to an ε-heat equation for pε and to the non-smooth rheology σε = pε Suε/(|Suε|+ε). Passing to the limit in the Galerkin dimension requires, at minimum, uniform (in the Galerkin index) bounds on ∂tuε and ε∂tpε that must be derived before the approximate solution is known to exist; the paper only proves such bounds for an already-existing solution in Proposition 3. In particular, identifying the weak limit of ∫ pε Suε/(|Suε|+ε) : Sϕ needs strong compactness of Suε or pε, which the energy estimate (24) alone does not provide. The reference to [2] may cover a similar model, but the exact adaptation to this system with the p=3 viscosity is not shown. This concern is independent of the separate scope question about whether the technical viscosity is physical.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes to prove existence and uniqueness of weak solutions for a simplified dense granular flow model that couples a pressure-dependent threshold rheology (Drucker-Prager type) with a dilatancy law. In Section 3, the authors define a weak formulation (Definition 1) that incorporates an energy inequality, show in Proposition 1 that this weak formulation recovers the full rheological relation, and then state Theorem 1: for initial data in L² and forcing in L^{3/2}(0,T;W^{-1,3/2}), there exists a weak solution, with u and p uniquely determined. The proof constructs a sequence of approximate problems (17)--(20) involving a nonlinear viscosity |Du|Du, a regularized yield law σε = pε Suε/(|Suε|+ε), and an ε-heat equation for the pressure. After establishing uniform estimates (Proposition 3) and compactness (Proposition 4), the authors pass to the limit using convexity arguments, Zhikov's lemma, and careful weak-strong identifications. Section 4 sketches an extension to a model with variable volume fraction (Theorem 2) and discusses the relation of the simplified system to the µ(I)-rheology via an asymptotic expansion.","tokens_in":21025,"tokens_out":18727,"duration_ms":138407,"significance":"If the gaps noted below are closed, the main theorem would be a valuable contribution: the proof strategy is innovative in exploiting the coupling between rheology and dilatancy to obtain dissipated energy, and the weak-strong limit identification via convexity and Zhikov’s lemma is elegant. The paper is largely self-contained for the simplified model and gives a complete uniqueness proof. However, two load-bearing issues remain: the existence proof of the ε-approximate system is only sketched by reference to a Galerkin method, and the identification V(p)=√p is not fully established on the zero set of p. The result is also conditional on the technical nonlinear viscosity |Du|Du, which is not part of the original inviscid granular model; the authors acknowledge this in Remark 3, but the abstract’s claim of a “real breakthrough” should be tempered accordingly. Overall, the paper is promising and the main ideas are sound, but the proof as written is not yet complete.","major_comments":[{"comment":"The proof of Proposition 2 is only a sketch: it derives the energy identity (24) and states that existence follows from a classical Galerkin method, citing [2]. This is load-bearing because Theorem 1’s limit passage in §3.5–3.6 requires a fully constructed sequence (uε,pε,σε). The Galerkin passage for the coupled system (22)–(23) needs to be demonstrated: one must show uniform bounds at the Galerkin level, pass to the limit in the finite-dimensional approximations, and identify the nonlinear term σε = pε Suε/(|Suε|+ε). A bare reference to [2] is not sufficient unless the applicability to this specific coupled system is spelled out. Please either provide the full Galerkin argument or state precisely which theorem in [2] covers the present system and why.","section":"§3.3, Proposition 2 (approximate existence)"},{"comment":"The conclusion “equality B = 0 implies V(p) = √p” is not justified on the set where p = 0. Corollary 1 gives V(p) ≤ √p and pV(p) ≥ p√p; since p ≥ 0, these imply V(p) = √p only on {p > 0}. On {p = 0} they only give V(p) ≤ 0. The term B = ∫(pV(p) − pV(p)) vanishes identically on {p = 0}, because both products contain the factor p, so B = 0 imposes no constraint on V(p) there. Since Vε(0) = 0 but Vε(x) = x/ε for x < 0, the weak limit V(p) could in principle carry a negative contribution supported on {p = 0} (e.g. if pε ≈ −εχE on a set E of positive measure). Without V(p) = √p a.e., equation (29) is not the target equation (14). An additional argument is needed, for instance a sign or monotonicity estimate showing that such negative contributions vanish, or a modification of the approximation so that Vε is bounded below so that V(p) is forced to be nonnegative.","section":"§3.6, Corollary 1 and the paragraph after (33)"},{"comment":"Theorem 2 is stated as a theorem, but its proof is only sketched in “Ideas for the proof” and relies on formal energy estimates for regular solutions (Propositions 5 and 6). The passage to weak solutions for the ξ-regularized system (37)–(40) is not carried out; in particular, the strong compactness of φε and the identification of the nonlinear terms (φmax − φ)√p and Dt(φ,u) are not demonstrated. Since Remark 6 already notes that the argument breaks down when ξ = 0, the status of Theorem 2 should be clarified: either supply the missing compactness and limit-identification details, or state Theorem 2 as a formal extension rather than a proved result.","section":"§4.1, Theorem 2 and Remark 6"}],"minor_comments":[{"comment":"The notation for the weak limits is ambiguous: the weak limit of pεVε(pε) and the product of p with the weak limit of Vε(pε) are both written as “pV(p)”, differing only by an overline that is lost in the displayed equations. Please use distinct symbols, e.g. \\overline{pV(p)} and p\\overline{V(p)}, to avoid confusion in the argument.","section":"§3.6, notation in (31) and (33)"},{"comment":"The phrase “This work represents a real breakthrough” is stronger than the mathematical content warrants, given that the main theorem requires the technical nonlinear viscosity (Remark 3) and that the extension Theorem 2 is only sketched. I recommend tempering this claim and stating explicitly in the abstract that the existence result is proved for a regularized simplified model.","section":"Abstract and Introduction"},{"comment":"The derivation of the simplified system (55)–(57) from the µ(I)-rheology model is purely formal. The text should state clearly that this is a formal asymptotic reduction, not a rigorous singular-limit result, to avoid giving the impression that Theorem 1 applies to the original inviscid model.","section":"§4.2, asymptotic derivation"},{"comment":"The last sentence of the proof of Proposition 2 refers to Subsection 3.7 for uniqueness, but Subsection 3.7 uses Proposition 1, which is specific to the limiting rheology (10), not to the regularized relation (18). Either give a direct uniqueness argument for the approximate system or omit the uniqueness claim from Proposition 2, since existence is all that is needed for the main theorem.","section":"§3.3, Proposition 2 uniqueness"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of math.AP and the central idea is interesting. However, the two main gaps—the missing Galerkin proof for Proposition 2 and the incomplete identification of V(p) on {p=0}—are load-bearing and should be fixed before publication. The third issue, the sketched proof of Theorem 2, is also significant, although one might accept it as a roadmap if the paper explicitly marks it as such. I encourage the editor to request a major revision rather than rejection, as the gaps appear fixable and the core methodology is sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the paper. The genuinely new thing is Definition 1: putting the threshold rheology into an energy inequality instead of the constitutive relation, then recovering σ:Su = 2p|Su| from the weak form via Proposition 1. That is a clever and useful trick, and the existence proof built around it — Galerkin approximations, uniform L^3 estimates, convexity/Zhikov identification of weak limits, pressure positivity — is mostly careful and correct. The uniqueness argument via monotonicity is also fine. Credit where due: this is the first existence/uniqueness theory for a Drucker-Prager-type granular model with dilatancy coupling, and the authors are honest that the nonlinear viscosity is a technical device (Remark 3).\n\nThe soft spots are real but not fatal. First, Proposition 2 — existence for the ε-approximate system — is asserted on a citation to [2] plus a single energy identity. That is load-bearing: Theorem 1 needs a rigorous sequence (uε,pε,σε), and the paper does not show the Galerkin construction or the uniform-in-dimension estimates required before passing to a solution of the approximating system. It may be standard, but it is not demonstrated here. A referee should ask for a complete proof or a very precise adaptation of [2]. Second, the scope is overstated. The abstract says 'real breakthrough' and 'granular model', but the theorem is for a simplified system with a technical viscosity; the physically original model has no viscosity, and the authors say the L^3 framework depends on it. The extension to the full model (Theorem 2) only works with a positive ξ regularization, and the limit ξ→0 is explicitly left open. These are honest limitations — they are flagged — but they should be reflected in the abstract and introduction. Third, the regularity framework is tight: products like p|Su| make sense only in L^{3/2}, and any weakening of the viscous term would break the argument. That is fine as a theorem, but it narrows the claim.\n\nBottom line: worth a serious referee. With a filled-in Proposition 2 and a toned-down abstract, it becomes a solid contribution to the analysis of non-isochoric granular models. I would send it to review rather than desk reject, but the referee should be told to verify the Galerkin step.","headline":"Genuinely new weak formulation and careful existence proof for a simplified granular model, but the approximate-solution step is cited rather than proved and the abstract overstates the scope.","tokens_in":21535,"tokens_out":2353,"would_cite":false,"duration_ms":28298,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35D30","35Q35","76T25"],"pacs":[],"model":"deepseek-v4-flash","headline":"A simplified granular flow model coupling threshold rheology with dilatancy is shown to have a weak solution, with velocity and pressure uniquely determined.","keywords":["granular flow","threshold rheology","dilatancy","weak solution","existence and uniqueness","pressure positivity","non-Newtonian viscosity","mu-I rheology"],"falsifier":"Find a weak solution of (9)–(12) in the sense of Definition 1 whose pressure is negative on a set of positive measure. Because the energy inequality together with $|\\sigma|\\le p$ forces $\\sigma:Su = 2p|Su|$ almost everywhere, such a solution would contradict Proposition 1 and therefore Theorem 1.","tokens_in":20471,"feed_emoji":"🧱","tokens_out":6823,"duration_ms":56715,"temperature":0.7,"pith_summary":"This paper establishes an existence and uniqueness theory for a PDE model of dense granular flow that combines threshold rheology with dilatancy. The central result, Theorem 1, shows that for a square-integrable initial velocity and a force field in a natural dual space, the simplified model (9)–(11) has a weak solution on a bounded three-dimensional domain, and that the velocity and pressure are unique. The mechanism is that the coupling of the pressure-dependent yield threshold with the dilatancy law $\\mathrm{div}\\,u = 2|Su| - \\sqrt{p}$ produces a dissipated-energy inequality strong enough to control the pressure, keep it non-negative, and identify all nonlinear weak limits. The same strategy extends to a model with evolving volume fraction, provided a small regularizing term is added; the original inviscid granular model is not covered.","feed_headline":"Existence proved for granular flows with threshold rheology","feed_subtitle":"Rheology plus dilatancy yields existence, uniqueness, and positive pressure via energy dissipation.","key_machinery":"The load-bearing object is the energy inequality obtained by testing the momentum equation with $u$ and the dilatancy equation with $p$. The nonlinear viscous term $\\mathrm{div}(2|Du|Du)$ makes the strain-rate tensor live in $L^3$, giving enough regularity to give meaning to products such as $p|Su|$ and to run the compactness argument. The threshold relation is encoded through the regularized stress $\\sigma_\\varepsilon = p_\\varepsilon Su_\\varepsilon/(|Su_\\varepsilon|+\\varepsilon)$, while the square-root pressure law is approximated by concave functions $V_\\varepsilon$; a convexity lemma and a flux-identification lemma then pin down the weak limits, and the dilatancy constraint $\\mathrm{div}\\,u = 2|Su|-\\sqrt{p}$ is recovered exactly.","core_discovery":"For the simplified system (9)–(12), the authors prove that a weak solution $(u,p,\\sigma)$ exists whenever $u_{\\rm init}\\in L^2(\\Omega)$ and $f\\in L^{3/2}(0,T;W^{-1,3/2}(\\Omega))$, and that $u$ and $p$ are unique. The weak formulation (13)–(15) is deliberately built so that the full threshold rheology $\\sigma:Su = 2p|Su|$, $|\\sigma|\\le p$ is recovered from the energy inequality alone, even though it is not imposed directly in the definition of weak solution. Existence is obtained by passing to the limit in a regularized system, with a convex-analysis and monotonicity argument identifying the weak limits of $|Du_\\varepsilon|Du_\\varepsilon$, $V_\\varepsilon(p_\\varepsilon)$, and $|Su_\\varepsilon|$. A notable part of the proof is that the pressure is shown to be non-negative almost everywhere, a property the model needs physically.","pith_inferences":["The technical nonlinear viscosity $\\mathrm{div}(2|Du|Du)$ is the real price of the theorem: the physical model from the paper's starting point has no such term, and the proof gives no indication how to replace the resulting $L^3$ regularity with a weaker mechanism.","The pressure-as-Lagrange-multiplier viewpoint noted in Remark 4 points to a variational route: treating the dilatancy relation as a constraint in a minimization problem could yield the stress relation without the viscous regularization, if the constraint can be handled in the non-smooth regime.","The same energy-dissipation strategy is likely transferable to other threshold rheologies, such as $\\mu(J)$-rheology for immersed granular flows, but the dilation closure and stability conditions would need to be verified case by case.","A concrete testable extension is to build a numerical scheme that preserves the discrete analogue of the energy inequality (15) and the bounds on $\\varphi$; if such a scheme exists, the theorem's structural assumptions are strong enough for practical simulation."],"forward_implications":["If the theorem is correct, the simplified granular model is well-posed for the stated data class, giving numerical discretizations a target problem with existence and uniqueness guarantees.","The proof delivers non-negative pressure as a theorem, not an assumption, which is physically natural and rarely proved rigorously for such models.","The same energy argument works for a volume-fraction model when a small regularizing term $\\xi(\\Delta\\varphi - \\varphi\\sqrt{p})$ is added, and the paper shows that the simplified system is the leading-order approximation of the full $\\mu(I)$-rheology model in a specific regime.","The uniqueness result is restricted to velocity and pressure; the stress is not claimed unique, which suggests that numerical methods should treat stress as a derived quantity.","The inviscid limit and the case $\\xi=0$ are explicitly left open, so the theorem marks a first step rather than a full theory of the original physical model."],"supporting_citations":[{"why":"Supplies the full granular model and the observation that the rheology-dilatancy coupling gives stability and energy dissipation.","marker":"[14]"},{"why":"Source of the dilatancy law relating volume change to shear and pressure, which becomes the constraint $\\mathrm{div}\\,u = 2|Su|-\\sqrt{p}$.","marker":"[27]"},{"why":"Provides the stability conditions and the equilibrium volume-fraction closure used in the model design and in the generalized system.","marker":"[28]"},{"why":"Gives the flux-identification lemma used to identify the weak limit of $|Du_\\varepsilon|Du_\\varepsilon$ and $|Du_\\varepsilon|^3$.","marker":"[29]"},{"why":"Supplies the compactness and time-derivative tools used to extract strongly convergent subsequences and to justify duality pairings.","marker":"[11]"},{"why":"Provides the Korn-type inequality in $W^{1,3}$ used to control the velocity gradient by $Du$ in the energy estimate.","marker":"[20]"},{"why":"Supplies the standard regularization strategy for the stress and the Galerkin approximation framework for the approximate system.","marker":"[2]"}],"fun_headline_variants":["Granular weak solutions exist via energy dissipation","Threshold rheology and dilatancy enable weak solutions","Existence and uniqueness for granular flow equations","Energy dissipation proves granular solution existence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof hinges on the nonlinear viscosity term $\\mathrm{div}(2|Du|Du)$, which is added for technical reasons and is absent from the original physical granular model; without it, the $L^3$ regularity of the strain rate and the identification of the nonlinear terms fail.","fun_headline_variants_meta":{"raw":{"variants":["Granular weak solutions exist via energy dissipation","Threshold rheology and dilatancy enable weak solutions","Existence and uniqueness for granular flow equations","Energy dissipation proves granular solution existence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000703,"raw_usage":{"total_tokens":3116,"prompt_tokens":837,"completion_tokens":2279,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":453,"completion_tokens_details":{"reasoning_tokens":2237}},"tokens_in":453,"tokens_out":2279,"duration_ms":14580,"temperature":1.0,"reasoning_tokens":2237,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:45:30.544397+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Find a weak solution of (9)–(12) in the sense of Definition 1 whose pressure is negative on a set of positive measure. Because the energy inequality together with $|\\sigma|\\le p$ forces $\\sigma:Su = 2p|Su|$ almost everywhere, such a solution would contradict Proposition 1 and therefore Theorem 1.","supporting_citations":[{"cited_title":"Chupin and T","cited_arxiv_id":null,"evidence_quote":"Supplies the full granular model and the observation that the rheology-dilatancy coupling gives stability and energy dissipation."},{"cited_title":"Roux and F","cited_arxiv_id":null,"evidence_quote":"Source of the dilatancy law relating volume change to shear and pressure, which becomes the constraint $\\mathrm{div}\\,u = 2|Su|-\\sqrt{p}$."},{"cited_title":"Schaeffer, T","cited_arxiv_id":null,"evidence_quote":"Provides the stability conditions and the equilibrium volume-fraction closure used in the model design and in the generalized system."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the flux-identification lemma used to identify the weak limit of $|Du_\\varepsilon|Du_\\varepsilon$ and $|Du_\\varepsilon|^3$."},{"cited_title":"Boyer and P","cited_arxiv_id":null,"evidence_quote":"Supplies the compactness and time-derivative tools used to extract strongly convergent subsequences and to justify duality pairings."},{"cited_title":"Malek, J","cited_arxiv_id":null,"evidence_quote":"Provides the Korn-type inequality in $W^{1,3}$ used to control the velocity gradient by $Du$ in the energy estimate."},{"cited_title":"Abbatiello, T","cited_arxiv_id":null,"evidence_quote":"Supplies the standard regularization strategy for the stress and the Galerkin approximation framework for the approximate system."}],"review_version":1}