{"id":"8c9a1112-fa1a-41f1-b8f9-d51e807840fe","arxiv_id":"2607.17380","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A finite-strain elasto-viscoplastic framework is implemented in the material point method, with analytic and numerical comparisons showing modified Cam-clay can reproduce Drucker-Prager behavior without excessive volume expansion.","lead":"This paper builds a computational framework for yield-stress fluids that can switch between different yield surfaces, and shows how one model (modified Cam-clay) can replace another (Drucker-Prager) without the usual volume-expansion artifacts. It may help engineers simulate mud, snow, and granular flows more reliably.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"MCC-to-DP retrieval rests on unproven critical-state limit in §4.1; no analytic reduction or parameter-convergence evidence, and overstress states are not on the critical-state line.","rationale":"The reader and I identify the same load-bearing concern: the MCC-to-DP equivalence is supported only by the §4.1 heuristic about the critical-state limit, not by a derivation or a systematic parameter-convergence study. The manuscript has independent support—analytic velocity profiles, numerical benchmarks, and open-source code—and no internal inconsistency forces rejection. But the headline claim's utility depends on a limit that is not established. The proposed check would either prove the reduction or reveal the extra conditions (small t_v, stress-rate bound) needed, determining whether the verdict should remain conditional. Since the reader's verdict is already CONDITIONAL and my concern does not move it further, the verdict should remain unchanged.","tokens_in":16858,"tokens_out":13474,"duration_ms":156990,"concrete_test":"Derive the asymptotic reduction claimed in §4.1 directly from Eqs. (9)–(10) and the associative flow rule: take p_c^0→0, ξ→∞ along the distinguished path p_c^0 ξ = const (or another that keeps the hardening rate finite), and show that the effective yield stress at mean pressure p tends to μp and tr(l^P)→0, so the MCC return mapping reduces to the DP-with-von-Mises-potential update (Eq. 28). If the reduction is found to require additionally t_v→0 or a bound on dp/dt, the 'right conditions' in the abstract are narrower than stated and the claim needs qualification.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4.1 proposes that with p_c^0 small and ξ large, plastic states will reside around the critical state line determined by µ, and uses this to conclude that MCC retrieves the DP solution. This is the single load-bearing step for the abstract's central claim. It is not derived: the yield function (9) collapses to a point as p_c^0→0 and the hardening law (10) becomes infinitely stiff as ξ→∞, and no distinguished-limit analysis shows that the composition converges to the DP cone uniformly over the pressures occurring in free-surface flows. In the inclined-plane benchmark the stress ratio is fixed by gravity at τ/p = tan θ (θ=40°, ratio 0.839), while the critical-state slope is µ = tan30° (0.577); therefore the condition cannot be that stress states lie on the critical-state line. The operationally required mechanism is that the ellipse apex (p_c/2, μp_c/2) tracks the current mean stress, i.e. p_c ≈ 2p, via the hardening law. The paper gives no timescale or error estimate for this tracking, and the numerical evidence covers a narrow parameter window (t_v=10^-4 s, one geometry) with no convergence study in (p_c^0, ξ). Figs. 6b/c show departure when p_c or β increases, but no validity region is quantified. Hence the central claim is currently an observation for selected parameters, not a demonstrated property of the model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a finite-strain overstress elasto-viscoplastic framework accommodating von Mises, Drucker–Prager, and modified Cam–clay yield surfaces, implemented in the material point method. It derives analytic steady-state velocity profiles for inclined-plane flow for Peric and Duvaut–Lions overstress models, notes a connection to the Bagnold profile for s=2, introduces a volume-correction technique for the Drucker–Prager cone-tip projection, and proposes that modified Cam–clay with small initial compressive strength and large hardening can retrieve Drucker–Prager behavior while avoiding spurious volume increase. Numerical benchmarks cover inclined-plane profiles, dam breaks, stretching flows, and energy/resolution convergence.","tokens_in":17231,"tokens_out":6249,"duration_ms":63717,"significance":"If substantiated, the proposed use of modified Cam–clay as a volume-stable regularization of Drucker–Prager would be practically useful for free-surface elasto-viscoplastic flow simulations. The paper’s strengths include analytic benchmark profiles that are used as independent references, an open-source implementation, clear numerical comparisons, and a useful discussion of the Drucker–Prager volume-expansion problem. However, the central MCC-to-DP equivalence currently rests on an unproven heuristic supported only by a narrow set of numerical experiments; the claimed retrieval is therefore not yet established as a property of the model. The paper is transparent about this limitation, but the abstract and conclusions state the result more strongly than the evidence supports.","major_comments":[{"comment":"The central claim that modified Cam–clay ‘retrieves’ the Drucker–Prager solution rests on the heuristic that with p0_c small and ξ large, plastic states reside around the critical state line. This is not derived. For the inclined-plane benchmark the stress state satisfies τ/p = tan θ (0.839 for θ=40°), while the critical-state slope is µ=tan30°=0.577; overstress states are necessarily above the critical-state line, so the stated mechanism cannot be ‘residing on the critical state line.’ The operative mechanism must instead be that the ellipse apex tracks the current mean stress (pc ≈ 2p) through Eq. (10). No timescale, error estimate, or distinguished limit is given for this tracking. Please replace the heuristic with an analysis, or at least a quantified convergence study in (p0_c, ξ), and state its regime of validity. As written, the abstract’s claim is an observation for selected para","section":"§4.1 and Abstract"},{"comment":"The numerical evidence for DP–MCC equivalence covers a narrow window: one flow geometry (infinite inclined plane), one set of elastic parameters, and no systematic variation of (p0_c, ξ). Figs. 6b/c show departure when pc or β is increased, but no quantitative criterion for acceptable agreement is given. In the dam-break comparison (§4.2), the statement that DP, DP-corr., and MCC show ‘no significant difference’ is not quantified (e.g., by runout distance, final shape, or material-distribution metric). Please add a quantitative agreement measure and a parameter-space scan that delineates the ‘right conditions’ claimed in the abstract.","section":"Fig. 6 and §4.2"}],"minor_comments":[{"comment":"The condition stated after Eq. (19), ‘we must assume µ > tan(θ) for flow’, appears to have the wrong inequality. For the non-cohesive DP model, the overstress ratio is τ/τ_y = tanθ/µ (Eq. (25)); sustained flow requires µ < tanθ, and the simulations use µ=tan30° with θ=40° (µ<tanθ). If µ>tanθ, hy from Eq. (18) lies above the free surface and the shear stress never reaches yield. This typo should be corrected; it does not affect the simulations but obscures the benchmark setup.","section":"Eq. (19)"},{"comment":"The initial compressive strength is defined as p0_c in Eq. (10) but denoted pc in §4.1 and Fig. 6. Please make the notation consistent or explicitly state that pc in §4.1 is the initial value.","section":"§2.1/§4.1"},{"comment":"The caption contains a garbled phrase: ‘as described in from Section 4.3’ should be ‘as described in Section 4.3’.","section":"Fig. 10 caption"},{"comment":"The discussion states that Cam–clay models ‘can also be used (with other parameter choices) to produce very different results’ than Drucker–Prager. This is acknowledged, but it would help to give one concrete example (e.g., a parameter regime where MCC behaves as a capped compressible model rather than as a regularized DP cone) so that the reader can calibrate the scope of the claimed equivalence.","section":"§5"},{"comment":"The statement that corrected Drucker–Prager and modified Cam–clay give ‘almost identical volumetric behavior’ is visual. Reporting the maximum relative volume difference, or plotting the two curves with a shared quantitative scale, would make the claim precise.","section":"Fig. 11"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You'll get value from this paper if you work on EVP models or MPM. The analytic steady-state velocity profiles for Peric and Duvaut-Lions with von Mises and Drucker-Prager yield are the real contribution; they're derived cleanly, match the simulations, and the connection to Bagnold is a nice payoff. The framework itself is not revolutionary—overstress viscoplasticity with three standard yield surfaces—but the unified finite-strain treatment is well put together, and the open-source implementation is a plus.\n\nThe central claim, that modified Cam-clay can be used to retrieve Drucker-Prager solutions without the cone-tip volume artifact, is more fragile. The paper proposes in §4.1 that with small pc and large ξ, plastic states reside around the critical state line. That heuristic is not derived, and the stress-test note gets this right: in the inclined-plane benchmark, τ/p = tan(40°) ≈ 0.84, while the critical-state slope μ = tan(30°) ≈ 0.58, so the states are not on the critical state line. The actual mechanism is that the hardening law drives the ellipse apex to track the mean stress (pc ≈ 2p), making the yield stress at the operating pressure approximate μp. The paper gives no error estimate or convergence study in (pc0, ξ) for that tracking, and the numerical evidence covers a narrow window (one geometry, tv = 10^-4 s). That doesn't kill the paper—the claim is explicitly 'under the right conditions'—but the abstract reaches a bit beyond what is demonstrated. A referee should ask for a validity region or a distinguished-limit argument.\n\nMinor points: there is a sign error in §2.3—'we must assume µ > tan(θ) for flow' should be μ < tanθ. The derivations that follow use the correct sign, so it's a typo, but it will trip readers. Reproducibility would be improved by pinning commit hashes and per-figure data, but the code DOI is already better than most.\n\nOverall, the paper is solid, honest about its scope, and the analytic benchmarks are genuinely new. It deserves a serious referee, but the MCC-to-DP claim needs tightening before publication.","headline":"Useful EVP comparison with a genuinely new analytic benchmark, but the MCC-to-DP claim rests on a heuristic that needs tightening.","tokens_in":17683,"tokens_out":5836,"would_cite":true,"duration_ms":58991,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Under the right parameter regime, modified Cam-clay reproduces Drucker-Prager flow without the volume blow-up.","keywords":["Elasto-viscoplastic flows","Free-surface flows","Yield surfaces","Overstress model","Modified Cam-clay","Drucker-Prager","Finite strain","Material point method"],"falsifier":"Run an inclined-shear-flow simulation with the modified Cam–clay model using a small hardening parameter ξ or an initial pc comparable to the gravitational confining pressure. If the steady-state velocity profile still matches the Drucker–Prager analytic solution while volume stays constant, the claim is robust; if the profile departs or the material compacts noticeably, the circumvention fails outside the tuned regime. The paper's own figures already show such departure when pc and β are increased, so the decisive test is to map the boundary of the matching regime in (pc, ξ) space.","tokens_in":16763,"feed_emoji":"🧱","tokens_out":6860,"duration_ms":69579,"temperature":0.7,"pith_summary":"The paper sets out a single finite-strain, overstress-type elasto-viscoplastic framework that can host any of three yield surfaces—pressure-independent, frictional, and capped—and compares how each behaves in free-surface flows. Its central claim is that the capped modified Cam–clay model, given a small initial compressive strength and a sufficiently fast hardening law, drives plastic states to the critical-state line and therefore reproduces the frictional Drucker–Prager solution while eliminating the troublesome volume expansion that Drucker–Prager produces under large deformation. A sympathetic reader would care because this turns a geotechnical constitutive model into a practical, volume-stable substitute for Drucker–Prager in yield-stress fluid simulations, removing the need for a special cone-tip projection and its spurious dilatancy. The paper also derives analytic steady-state velocity profiles for the two overstress variants and shows that one of them contains the familiar granular-flow velocity profile as a special case. The equivalence claim is demonstrated numerically for shear flow, dam break, and stretching flows, not proven from the constitutive equations.","feed_headline":"Cam-clay sidesteps Drucker-Prager flow blow-up","feed_subtitle":"A capped yield surface produces the same frictional flow solution while avoiding the excessive volume gain of large deformations.","key_machinery":"The load-bearing object is the modified Cam–clay yield surface, an ellipse in the pressure–shear-stress plane defined by τy(p) = µ sqrt((p − pt)(pc − p)), where pc is an isotropic compressive strength that hardens as pc(εP_V) = p0_c exp(−ξ εP_V) with accumulated plastic volumetric strain εP_V. A sufficiently stiff hardening law and small initial pc force stress states toward the apex of the ellipse, which lies on the critical state line of slope µ—exactly where Drucker–Prager plastic states would sit. Surrounding that is the paper's finite-strain overstress framework with a non-associative deviatoric plastic potential, which lets all three yield surfaces share one elastic predictor–plastic c","core_discovery":"On the paper's own terms, the central discovery is that the singularity of the Drucker–Prager yield surface—the cone tip that forces a special, volume-increasing return mapping—can be sidestepped without changing the observable flow solution. By replacing Drucker–Prager with the modified Cam–clay yield surface and choosing a large hardening parameter ξ and a small initial compressive strength pc, the plastic stress states remain close to the critical state line set by the friction slope µ. In that regime the model behaves as an effective Drucker–Prager fluid, matching the analytic inclined-plane velocity profile and the dam-break runout in the paper's tests, while the capped yield surface pe","pith_inferences":["The equivalence between Cam-clay and Drucker-Prager is a tuned-regime statement, not a material equivalence: outside the small-pc, large-ξ regime, Cam-clay intentionally produces different, more compressible flows, so practitioners should treat parameter selection as part of the model definition.","A natural next test is to calibrate the hardening law ξ against measured compaction in a collapsing column or debris flow; the paper's framework predicts that the volume-change history discriminates between corrected Drucker-Prager and Cam-clay even when final runout matches.","Because the constitutive framework is not tied to the particular particle method used here, the same yield-surface comparison could be reproduced in any large-deformation continuum solver, which would test whether the Drucker-Prager-to-Cam-clay equivalence is numerical-scheme independent.","For yield-stress fluids whose elastic response matters at startup, the framework suggests that a capped surface may regularize both the solid-fluid transition and the pressure singularity simultaneously, a combination not offered by the uncapped frictional surface."],"forward_implications":["If the claimed equivalence holds, modified Cam–clay can replace Drucker–Prager in elasto-viscoplastic free-surface flow simulations, giving the same velocity profiles and runout without the special cone-tip volume correction.","The uncorrected Drucker–Prager model is unsafe for large-deformation flows: in the paper's stretching test it gains roughly 50% volume, while corrected Drucker–Prager and modified Cam–clay stay close in volume.","The cone-tip correction technique and the modified Cam-clay model give almost identical volumetric behavior in the tested flows, so the correction is a viable alternative when a Cam-clay parameter set is not available.","The analytic velocity profiles provide ready verification targets: the projection-type overstress variant with exponent s = 2 reduces to a familiar granular-flow profile, linking the framework to granular rheology.","Yield surface choice, not just yield stress, controls plug formation and compressibility in the simulations."],"fun_headline_variants":["Capped yield surface bypasses Drucker-Prager singularity","Cam-clay cap matches Drucker-Prager flow without blow-up","Modified Cam-clay dodges Drucker-Prager cone tip","Capped surface avoids Drucker-Prager volume-gain singularity"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole circumvention rests on the heuristic that a small initial compressive strength pc and a large hardening parameter ξ keep plastic states near the critical state line; the paper verifies this numerically for the studied flows and parameter sets but does not derive it from the constitutive equations.","fun_headline_variants_meta":{"raw":{"variants":["Capped yield surface bypasses Drucker-Prager singularity","Cam-clay cap matches Drucker-Prager flow without blow-up","Modified Cam-clay dodges Drucker-Prager cone tip","Capped surface avoids Drucker-Prager volume-gain singularity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000721,"raw_usage":{"total_tokens":3093,"prompt_tokens":783,"completion_tokens":2310,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":527,"completion_tokens_details":{"reasoning_tokens":2235}},"tokens_in":527,"tokens_out":2310,"duration_ms":14714,"temperature":1.0,"reasoning_tokens":2235,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T18:07:42.385363+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run an inclined-shear-flow simulation with the modified Cam–clay model using a small hardening parameter ξ or an initial pc comparable to the gravitational confining pressure. If the steady-state velocity profile still matches the Drucker–Prager analytic solution while volume stays constant, the claim is robust; if the profile departs or the material compacts noticeably, the circumvention fails outside the tuned regime. The paper's own figures already show such departure when pc and β are increased, so the decisive test is to map the boundary of the matching regime in (pc, ξ) space.","supporting_citations":[],"review_version":1}