REVIEW 3 major objections 4 minor 29 references
Weak solution for granular model
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A simplified granular flow model coupling threshold rheology with dilatancy is shown to have a weak solution, with velocity and pressure uniquely determined.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [§3.3, Proposition 2 (approximate existence)] 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.
- [§3.6, Corollary 1 and the paragraph after (33)] 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.
- [§4.1, Theorem 2 and Remark 6] 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.
minor comments (4)
- [§3.6, notation in (31) and (33)] 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.
- [Abstract and Introduction] 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.
- [§4.2, asymptotic derivation] 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.
- [§3.3, Proposition 2 uniqueness] 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.
Circularity Check
No significant circularity: the existence proof is self-contained; [14] is contextual and [2] is an external Galerkin reference.
full rationale
Walking the derivation of Theorem 1: the approximate system (17)-(20) is asserted to have solutions by a Galerkin argument referenced to [2] (external authors, not a self-citation) combined with the energy identity (24); the uniform estimates in Proposition 3 follow from that identity; weak compactness in Proposition 4 and the convex-analysis Lemmas 1-2 identify the nonlinear limits. None of these steps is equivalent to the conclusion. In particular, the energy inequality (15) in Definition 1 is not assumed in order to force the rheology: it is derived from the approximate energy estimates and passed to the limit, and Proposition 1 only shows that the inequality |sigma| <= p together with the energy inequality forces the exact flow rule sigma:Su = 2p|Su| by comparing the energy identity with (15). The nonlinear viscosity is explicitly declared 'essentially technical' in Remark 3, so this is a transparent modeling assumption rather than a hidden ansatz. The self-citation [14] provides the physical model and stability context, but no proof step in Theorem 1 relies on [14] for a mathematical fact. The proof of Proposition 2 is admittedly sketched with 'see [2]', and Theorem 2 is proved through formal estimates, but these are exposition/rigor gaps rather than circular reductions: they do not make the claimed result equal to its inputs. Score 1 reflects only the minor contextual self-citation, which is not load-bearing.
Assumptions & free parameters
free parameters (2)
- normalized coefficients in (9)-(11) =
alpha0 = beta0 = nu0 = phi0 rho0 = 1
- regularization strength xi in (37)-(40) =
xi > 0, arbitrary
assumptions (8)
- standard math Korn's inequality in W^{1,3} for bounded domains
- standard math Aubin-Lions-Simon compactness lemma
- standard math Zhikov lemma on weak convergence of fluxes (Lemma 2)
- standard math Convex and concave weak lower and upper semicontinuity (Lemma 1)
- standard math Existence and uniqueness for the Galerkin approximate system (17)-(20)
- domain assumption Drucker-Prager rheology and Roux-Radjai dilatancy law as constitutive relations
- ad hoc to paper Non-linear viscosity 2|Du|Du is sufficient for the L^3 framework
- ad hoc to paper Artificial regularization xi H and modified derivative D_t in the volume-fraction model
Cite this review
Pith. "Pith review of Weak solution for granular model." pith.science (2026). https://pith.science/paper/S46WVIUN
@misc{pith2026250517588,
author = {Pith},
title = {Pith review of: Weak solution for granular model},
year = {2026},
howpublished = {\url{https://pith.science/paper/S46WVIUN}},
note = {Machine review of arXiv:2505.17588}
}
read the original abstract
This article is devoted to questions concerning the existence of solutions for partial differential equation problems modeling granular flows. The models studied take into account the complex threshold rheology of these flows, as well as the dilatance effects. It is the coupling of these two physical phenomena that ensures stability and the existence of dissipated energy. The key point of the article is to understand how this energy can ensure the existence of a weak solution. We first establish a complete result on a simplified model, then demonstrate how it can be extended to more general cases. This work represents a real breakthrough in the mathematical analysis of this type of models for complex flows.
Reference graph
Works this paper leans on
-
[2]
A. Abbatiello, T. Los, J. M´ alek, and O. Sou vcek. Three-dimensional flows of pore pressure-activated bingham fluids. Mathematical Models and Methods in Applied Sciences , 29(11):2089–2125, 2019
work page 2019
-
[1]
A. Abbatiello, M. Bul ´ ıvcek, T. Los, J. M´ alek, and O. Sou vcek. On unsteady flows of pore pressure- activated granular materials. Zeitschrift f¨ ur angewandte Mathematik und Physik, 72(1):6, 2021
work page 2021
-
[3]
B. Andreotti, Y. Forterre, and O. Pouliquen. Les milieux granulaires-entre fluide et solide: Entre fluide et solide . EDP sciences, 2012
work page 2012
- [4]
- [5]
- [6]
-
[7]
A. Beck. First-order methods in optimization . SIAM, 2017
work page 2017
-
[8]
F. Bouchut, E. D. Fernandez-Nieto, A. Mangeney, and G. Narb ona-Reina. A two-phase shallow debris flow model with energy balance. ESAIM: Mathematical Modelling and Numerical Analysis , 49(1):101– 140, 2015
work page 2015
Show all 29 references
-
[9]
Bouchut, E
F. Bouchut, E. D. Fern´ andez-Nieto, A. Mangeney, and G. Narbona-Reina. A two-phase two-layer model for fluidized granular flows with dilatancy effects. J. Fluid Mech. , 801:166–221, 2016
2016
-
[10]
Bouchut, E
F. Bouchut, E. D. Fern´ andez-Nieto, A. Mangeney, G. Narbo na-Reina, et al. Dilatancy in dry granular flows with a compressible µ(I)-rheology. Journal of Computational Physics , 429:110013, 2021
2021
-
[11]
Boyer and P
F. Boyer and P. Fabrie. Mathematical Tools for the Study of the Incompressible Navi er-Stokes Equations andRelated Models, volume 183. Springer Science & Business Media, 2012
2012
-
[12]
E. C. Breard, L. Fullard, J. Dufek, M. Tennenbaum, A. Fernan dez Nieves, and J. F. Dietiker. Inves- tigating the rheology of fluidized and non-fluidized gas-particle beds : implications for the dynamics of geophysical flows and substrate entrainment. Granul. Matter , 24(1):34, 2022
2022
-
[13]
H. Brezis. Analyse fonctionnelle. Th´ eorie et applications, 1983
1983
-
[14]
Chupin and T
L. Chupin and T. Dubois. Non-isochoric stable granular models ta king into account fluidisation by pore gas pressure. Journal of Fluid Mechanics , 979:A14, 2024
2024
-
[15]
Chupin, T
L. Chupin, T. Dubois, M. Phan, and O. Roche. Pressure-depen dent threshold in a granular flow: Numerical modeling and experimental validation. J. Non-Newton. Fluid Mech. , 291:104529, 2021
2021
-
[16]
Chupin and J
L. Chupin and J. Mathe. Existence theorem for homogeneous in compressible navier–stokes equation with variable rheology. European Journal of Mechanics-B/Fluids , 61:135–143, 2017
2017
-
[17]
Fang and Z
L. Fang and Z. Guo. Global weak solutions to a three-dimensiona l compressible non-newtonian fluid. Communications in Mathematical Sciences , 20(6):1703–1733, 2022
2022
-
[18]
Hutter and K
K. Hutter and K. Rajagopal. On flows of granular materials. Continuum Mechanics and Thermody- namics, 6:81–139, 1994
1994
-
[19]
P. Jop, Y. Forterre, and O. Pouliquen. A constitutive law for de nse granular flows. Nature, 441:727–730, 2006
2006
-
[20]
Malek, J
J. Malek, J. Necas, M. Rokyta, and M. Ruzicka. Weak and Measure-Valued Solutions to Evolutionary PDEs. Chapman and Hall, 1996
1996
-
[21]
G. MiDi. On dense granular flows. Eur. Phys. J. E , 14:341–365, 2004
2004
-
[22]
Montell` a, J
E. Montell` a, J. Chauchat, C. Bonamy, D. Weij, G. Keetels, and T. Hsu. Numerical investigation of mode failures in submerged granular columns, flow, 3, e28, 2023
2023
-
[23]
E. P. Montell` a, J. Chauchat, B. Chareyre, C. Bonamy, and T.- J. Hsu. A two-fluid model for immersed granular avalanches with dilatancy effects. Journal of Fluid Mechanics , 925:A13, 2021
2021
-
[24]
Radjai and F
F. Radjai and F. Dubois. Discrete-element modeling of granular materials . Wiley-Iste, 2011
2011
-
[25]
J. A. Robinson, D. J. Holland, and L. Fullard. Complex behavior in c ompressible nonisochoric granular flows. Phys. Rev. Fluid , 8(1):014304, 2023
2023
-
[26]
O. Roche. Depositional processes and gas pore pressure in py roclastic flows: an experimental perspec- tive. Bull. Volcanol., 74(8):1807–1820, 2012
2012
-
[27]
Roux and F
S. Roux and F. Radja ¨ ı. Texture-dependent rigid-plastic beha vior. In Physics of dry granular media , pages 229–236. Springer, 1998. 19
1998
-
[28]
Schaeffer, T
D. Schaeffer, T. Barker, D. Tsuji, P. Gremaud, M. Shearer, a nd J. Gray. Constitutive relations for compressible granular flow in the inertial regime. J. Fluid Mech. , 874:926–951, 2019
2019
-
[29]
V. Zhikov. On the weak convergence of fluxes to a flux. In Doklady Mathematics , volume 81-1, 2010. 20
2010
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.