{"id":"d15bf4de-8465-4a3f-89b2-6a4e744cf4f6","arxiv_id":"2501.10631","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A new flow rule, v = v∞[1−exp(−l/L)]^{1/2} with v∞ ∝ μ_r^{3/2}[(tanθ−μ_r)h]^{4/3}, describes undeveloped and fully developed dense granular flows down rough inclines in the tested data.","lead":"Researchers measured how fast dense sand, spheres, and their mixtures flow down a rough inclined chute and found a new relationship between speed, friction, incline, thickness, and travel distance. The resulting formula also predicts how far a flow must travel before reaching full speed, which could help estimate the destructive energy of avalanches.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The unified rule hinges on v0 ∝ sqrt(l), an untested rigid-block assumption for undeveloped granular flows.","rationale":"The reader's CONDITIONAL verdict is appropriate. The central claim is an empirical correlation with fitted exponents, and the weakest link is the rigid-block model for the undeveloped stage. This assumption is physically questionable, but the paper's multi-position data can test it directly. No internal inconsistency or obvious contradiction was found. The question of whether Eq. (7) is the exact spatial saturation law is secondary because the collapse plot is the empirical evidence, and the proposed check on Eq. (5) would also partially test the functional form. I therefore agree with the reader's weakest_assumption and keep the verdict unchanged.","tokens_in":7001,"tokens_out":9628,"duration_ms":103548,"concrete_test":"Use the existing measurements at l = 0.95, 1.3, and 1.7 m (the paper states these were collected) for runs in the undeveloped regime (l ≪ L). For each run, test the rigid-block prediction v0^2 = c0^2 (tanθ−μ_r) l: fit v0^2 = a + b l and check (i) a ≈ 0 and (ii) b ≈ constant c0^2 across materials, or equivalently v(l2)/v(l1) ≈ sqrt(l2/l1). A nonzero intercept a from entrance velocity or a deviating l-exponent would refute Eq. (5), requiring replacement of L in Eq. (8); a clean pass would strongly support the unified rule.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation of the development length L = (c∞^2/c0^2) μ_r^3[(tanθ−μ_r)h]^{5/3} h rests on Eq. (5): v0 = c0 sqrt(l(tanθ−μ_r)). This treats the undeveloped flow as a rigid block that starts from rest and slides with kinetic friction μ_r. Real hopper-fed chute flows enter through an opening gate with a finite velocity and develop a shear profile over a rough, no-slip base; the free-surface velocity is not that of a single sliding block. Any additive entrance velocity or h-dependent shear contribution changes the l-dependence of v0^2, so the cancellation leading to L and the exponents in Eq. (8) are not protected. The paper never validates Eq. (5) against its own multi-position measurements, and the DEM evidence cited [18] is for temporal, not spatial, saturation. Thus the central unification inherits an unverified physical assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports new laboratory chute measurements of the free-surface velocity of dense granular flows of spheres, two natural sands, spheres-sand mixtures, and a zirconia-based mixture down rough inclines. It argues that the standard flow rules, in which the Froude number scales with h/h_s or (tanθ/μ_r)^2 h/h_s, fail because they are not universal across granular materials and because many measurements are made before the flow is fully developed. The authors propose a unified flow rule v = v∞ [1 − exp(−l/L)]^{1/2}, with v∞ = c∞ μ_r^{3/2} [(tanθ − μ_r) h]^{4/3} and L = (c∞^2/c0^2) μ_r^3 [(tanθ − μ_r) h]^{5/3} h, derived from a dimensional-analysis argument motivated by a recent unification of aeolian and fluvial sediment transport. The rule is compared with the new data and with data from Ref. [6] and is reported to collapse the combined measurements for undeveloped and fully developed flows.","tokens_in":7307,"tokens_out":5221,"duration_ms":51729,"significance":"If the proposed rule holds, it would be a substantial step forward for dense granular chute flows: it replaces the two-parameter h/h_s-based flow rules with a rule built around the single material parameter μ_r, it provides a first measurements-based estimate of the development length L, and it explicitly addresses the pre-saturation regime that previous flow rules ignored. The paper also supplies new experimental data and a simple, falsifiable functional form that can be tested independently. The main strengths are the breadth of tested materials and the direct comparison with earlier experiments. However, the claimed universality is presently only partial: the exponents p and q in Eq. (6) are determined by fitting to the same combined data used to demonstrate the collapse, and the constants c∞ and c0 are fitted separately for nonmixtures and mixtures. The central physical assumption for the undeveloped stage, Eq. (5), is also not directly validated against the multi-position measurements that the paper itself reports.","major_comments":[{"comment":"Equation (5), v0 = c0 sqrt(l(tanθ − μ_r)), is the load-bearing assumption for the undeveloped stage, but it is not tested against the paper's own multi-position measurements. Treating the entire undeveloped flow as a rigid block sliding with kinetic friction μ_r is a strong assertion; dense granular flows generally develop internal velocity profiles and basal slip distributions, and the free-surface velocity need not follow a single sliding-friction law. Because Eq. (5) enters directly into Eq. (8) through L = l v∞^2 / v0^2, the predicted development length and the spatial saturation form are not grounded if Eq. (5) fails. I ask the authors to validate Eq. (5) directly using their measurements at l = 0.95 m, 1.3 m, and 1.7 m, or to justify the rigid-block approximation with additional evidence.","section":"Eq. (5)"},{"comment":"The exponents p ≈ 3/2 and q ≈ 4/3 in Eq. (6) are determined by weighted-least-squares fitting and Kendall rank maximization using the same combined data that are then displayed as a data collapse in Fig. 5, and c∞ and c0 are fitted separately for nonmixtures and mixtures. Consequently, the collapse in Fig. 5 is not an independent test of the functional form; it is a demonstration that a two-exponent, two-constant family can parameterize the data. The claim of universality would be strengthened substantially by validating Eq. (6) on a holdout subset of the data or against an independent data set, rather than by showing agreement with the same data used for calibration.","section":"Eq. (6) and Fig. 5"},{"comment":"The translation from the temporal saturation observed in previous DEM simulations (Ref. [18]) to the spatial saturation form in Eq. (7) is not derived rigorously. The step 'l = ∫ v0(t′) dt′ ≈ v0 t/2 and L ∝ v∞ T' assumes that v0(t′) is roughly constant during the undeveloped stage, which is inconsistent with the accelerating rigid-block motion described by Eq. (5). No direct evidence is provided that the spatial velocity profile follows [1 − exp(−l/L)]^{1/2}; the authors themselves note that 'other choices are also possible.' Since Eq. (7) is one of the central predictions, this gap in the derivation should be addressed, either by a more careful asymptotic argument or by a direct test against measured spatial profiles.","section":"Eq. (7) and preceding derivation"},{"comment":"The paper acknowledges that c∞ and c0 are different for nonmixtures and mixtures (c∞ = 11, c0 = 1.5 versus c∞ = 15, c0 = 1.7). The proposed rule therefore contains two material-class-dependent free constants, and the universality across granular materials is only partial. The text attributes the difference to segregation in the mixtures, but as written the model does not predict this difference. The authors should either state clearly that the rule is class-specific with separate calibrated constants, or provide a physical prediction for c∞ and c0 in terms of measurable mixture properties.","section":"Fig. 5 caption and text after Eq. (8)"}],"minor_comments":[{"comment":"In the sentence introducing Eq. (7), the expression 'l = ∫ v0(t′)t′' appears to be a typographical error; it should read 'l = ∫ v0(t′) dt′.'","section":"Eq. (7) derivation"},{"comment":"Figures 4 and 5 do not show measurement uncertainties. Since the quantitative claims rest on the quality of the data collapse, reporting run-to-run scatter or error bars would make the comparison more convincing.","section":"Figs. 4 and 5"},{"comment":"The numerical values c∞ = 11, c0 = 1.5 and c∞ = 15, c0 = 1.7 appear only in the Fig. 5 caption; they should be listed in the main text or in Table I, together with their uncertainties.","section":"Fig. 5 caption"},{"comment":"Reference [6] is cited as 'Physics Review E'; the correct journal name is 'Physical Review E.'","section":"Reference [6]"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the scope of the journal and the experimental dataset is valuable. My main concern is the gap between the strong universal claims and the fact that the central scaling exponents and constants are fitted on the same data used for validation, and that the undeveloped-stage assumption in Eq. (5) is not tested. These issues are fixable with additional analysis and a more cautious framing, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper gives a compact empirical flow rule for granular chute flows: v = v∞[1 − exp(−l/L)]^{1/2}, with v∞ ∝ μ_r^{3/2}[(tanθ − μ_r)h]^{4/3} and L ∝ μ_r^3[(tanθ − μ_r)h]^{5/3} h. The scaling is new, the data collapse in Fig. 5 is real, and the authors correctly show that the two standard rules (Pouliquen, Börzsönyi–Ecke) fail across materials and in the undeveloped regime. They deserve credit for combining new measurements with published data and for being transparent about the fitting procedure.\n\nThe genuinely new pieces are the μ_r^{3/2} dependence, the (tanθ − μ_r)^{4/3} h^{4/3} form for v∞, and the first measurement-based estimate of the development length L. The paper also makes a fair point that hs-based rules diverge at tanθ = μ_2 and that the existence of a phase transition there is not established.\n\nThe soft spots are real but not fatal. Most important is the early-stage velocity v0 = c0 sqrt(l(tanθ − μ_r)). This treats the undeveloped flow as a rigid sliding block, and the paper never directly validates it against its own multi-position measurements (l = 1.7, 1.3, 0.95 m). That assumption feeds directly into the derivation of L, so if it fails, the L scaling is not protected. The stress-test note is right about this. Second, the exponents p and q are fitted on the same data used to demonstrate the collapse, so the rule is not an independent prediction. Third, c∞ and c0 differ between nonmixtures and mixtures (11 vs 15, 1.5 vs 1.7), so the “unified” rule is only partially unified. These are addressable concerns, not contradictions.\n\nMy bottom line: this is a solid empirical Letter that deserves a serious referee. I would send it out, asking for a direct test of Eq. (5) using the multi-position data, an out-of-sample check of the fitted scalings, and some physical discussion of why the mixture constants differ. The Supplement (not available to me) should contain the uncertainty analysis and the comparison with µ(I)-rheology predictions; those details matter for judging how much weight to put on the exponent values.\n\nIf I needed a compact flow rule for chute-flow speed or runout estimates, I would cite this with caveats. It is a correlation with plausible dimensional motivation, not yet a derived law.","headline":"Useful empirical flow rule with fitted exponents; the early-stage rigid-block assumption is the main soft spot before calling it unified.","tokens_in":7877,"tokens_out":3164,"would_cite":true,"duration_ms":34244,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a single formula, v = v∞[1 − exp(−l/L)]^(1/2), describes the free-surface velocity of dense granular flows down rough inclines from release to full development, with v∞ and L set only by the dynamic friction…","keywords":["granular chute flows","free-surface velocity","flow rule","dynamic friction coefficient","development length","dense granular flow","avalanche","rough incline"],"falsifier":"Measure the free-surface velocity at several downstream positions along a chute long enough to include both undeveloped and fully developed states, for fixed material, thickness, and angle. If $v^2$ plotted against $l$ does not follow $v_\\infty^2[1 - \\exp(-l/L)]$ with $L = (c_\\infty^2/c_0^2) \\mu_r^3[(\\tan\\theta - \\mu_r)h]^{5/3} h$, or if the fitted $L$ does not scale as $[(\\tan\\theta - \\mu_r)h]^{5/3} h$, the rule is falsified. A more direct test of the weakest assumption is to image the velocity profile across the depth early in the flow: if the surface moves faster than the base, the rigid-block law is invalid.","tokens_in":6782,"feed_emoji":"🏔️","tokens_out":11233,"duration_ms":85426,"temperature":0.7,"pith_summary":"Dry, dense granular flows down rough inclines—the laboratory cousin of avalanches—take a finite distance to accelerate, and on many chutes they never reach a steady speed. The paper argues that existing flow rules miss this because they assume full development and because they are not universal across granular materials. It proposes a single flow rule, $v = v_\\infty[1 - \\exp(-l/L)]^{1/2}$, that describes both the undeveloped and fully developed regimes using only the dynamic friction coefficient $\\mu_r$, the excess driving stress $(\\tan\\theta - \\mu_r)h$, the thickness $h$, and the travel distance $l$. If correct, flow speed can be predicted without the stop thickness $h_s(\\theta)$, and the characteristic development length $L$ can be estimated from measurements for the first time. This matters because knowing how quickly a flow accelerates controls the energy and destructiveness of natural avalanches.","feed_headline":"One flow rule captures granular chute flows from start to steady state","feed_subtitle":"One formula fixes avalanche speed from material friction and slope geometry alone.","key_machinery":"The central machinery is the spatial saturation function $v = v_\\infty[1 - \\exp(-l/L)]^{1/2}$, obtained by merging the undeveloped velocity law $v_0 = c_0 \\sqrt{l(\\tan\\theta - \\mu_r)}$ with the fully developed velocity $v_\\infty = c_\\infty \\mu_r^{3/2}[(\\tan\\theta - \\mu_r)h]^{4/3}$. The development length $L = l v_\\infty^2 / v_0^2$ follows algebraically from the two limits. The dimensional analysis for $v_\\infty$ borrows from sediment transport: the steady-state velocity is a power law of the excess shear stress $(\\tan\\theta - \\mu_r)h$ times a dimensionless environment scale $\\mu_r$, with exponents $p = 3/2$ and $q = 4/3$ fixed by maximizing the Kendall rank correlation across the data. This machinery replaces the old $h/h_s$ scaling with a rule that contains only one material parameter, $\\mu_r$, and explicitly accounts for the finite chute length.","core_discovery":"The central claim is that the free-surface velocity $v$ of a dense granular flow down a rough incline, as a function of travel distance $l$, follows a square-root exponential saturation: $v = v_\\infty[1 - \\exp(-l/L)]^{1/2}$, where the fully developed velocity is $v_\\infty = c_\\infty \\mu_r^{3/2}[(\\tan\\theta - \\mu_r)h]^{4/3}$ and the development length is $L = (c_\\infty^2/c_0^2) \\mu_r^3[(\\tan\\theta - \\mu_r)h]^{5/3} h$. Here $\\mu_r$ is the dynamic friction coefficient, $h$ is the flow thickness, $\\theta$ is the inclination angle, and $c_0$ and $c_\\infty$ are proportionality constants that take one set of values for nonmixtures and another for shape-bidisperse mixtures. The rule merges an undeveloped-stage velocity $v_0 = c_0 \\sqrt{l(\\tan\\theta - \\mu_r)}$, the law of a rigid block sliding with friction $\\mu_r$, with a fully developed velocity that depends only on the excess shear stress and the environment scale $\\mu_r$, and not on the stop thickness $h_s(\\theta)$. The paper shows that this one formula collapses combined measurements of spheres, two sands, three spheres-sand mixtures, and a zirconia mixture equally well in the undeveloped and fully developed states.","pith_inferences":["Because the constants $c_0$ and $c_\\infty$ differ between nonmixtures and mixtures, an untested extension is whether they vary continuously with mixing ratio or grain-shape distribution; the paper leaves this open.","If the null hypothesis that no phase transition occurs at $\\mu_2$ holds, the flow rule should remain valid above the stop-thickness divergence, which could be tested in long chutes or with simulations at higher inclination angles.","The rigid-block assumption implies a testable prediction for the undeveloped stage: the surface velocity should scale as $\\sqrt{l}$ with a prefactor independent of $h$ and $\\theta$ at fixed $\\mu_r$, which high-speed imaging of the velocity profile could confirm or refute.","Connecting $L$ to natural avalanche runout gives an independent geological check: long-runout avalanches should show spatial velocity profiles that approach the exponential saturation, not linear-in-distance acceleration."],"forward_implications":["For short chutes or thick flows, where the development length $L$ exceeds the travel length $l$, the flow speed stays below the fully developed value; the rule provides a measurements-based estimate of $L$ for the first time.","The Froude number $v/\\sqrt{gh}$ is predicted to be controlled by $\\mu_r$, the excess shear stress, and the ratio $l/L$, not by the stop thickness ratio $h/h_s$.","The stop thickness $h_s(\\theta)$ and the second friction coefficient $\\mu_2$ drop out of the velocity prediction, so measuring only the dynamic friction coefficient $\\mu_r$ is sufficient.","Compared with the $\\mu(I)$-rheology estimate $L \\propto (\\tan\\theta - \\mu_r)h^3$, the data imply a substantially stronger dependence on $\\tan\\theta - \\mu_r$ and on $\\mu_r$, indicating that the standard rheology misses observed trends."],"supporting_citations":[{"why":"Supplies the previous free-surface velocity data for spheres and sand and the stop-thickness method that the new rule must reproduce.","marker":"[6]"},{"why":"Provides DEM simulation evidence that the velocity saturates as $v \\approx v_\\infty[1 - \\exp(-t/T)]$, the time-domain analogue of the spatial saturation function.","marker":"[18]"},{"why":"Borrows the insight that steady-state transport velocity depends only on excess shear stress and an environment scale, the basis of the dimensional analysis for $v_\\infty$.","marker":"[19]"},{"why":"One of the standard flow rules (Eq. 1) that the paper tests and finds inconsistent; supplies earlier chute data and the $h/h_s$ scaling baseline.","marker":"[4]"},{"why":"Two-dimensional DEM simulations showing steady uniform flows exist above $\\mu_2$, supporting the null hypothesis that $h_s$ need not enter $v_\\infty$.","marker":"[25]"}],"fun_headline_variants":["One rule unifies granular flow start and steady state","Granular flows obey same law from onset to full speed","Single formula describes avalanche development and steady state","Unified flow rule collapses granular chute data from start to finish","One equation captures granular flow speed from onset to saturation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that during the undeveloped stage the entire flow accelerates like a rigid block sliding on the incline with friction coefficient $\\mu_r$, giving $v_0 = c_0 \\sqrt{l(\\tan\\theta - \\mu_r)}$; if internal shear and basal slip develop early, this law fails and the predicted development length $L$ loses its foundation.","fun_headline_variants_meta":{"raw":{"variants":["One rule unifies granular flow start and steady state","Granular flows obey same law from onset to full speed","Single formula describes avalanche development and steady state","Unified flow rule collapses granular chute data from start to finish","One equation captures granular flow speed from onset to saturation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000541,"raw_usage":{"total_tokens":2669,"prompt_tokens":1100,"completion_tokens":1569,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":716,"completion_tokens_details":{"reasoning_tokens":1491}},"tokens_in":716,"tokens_out":1569,"duration_ms":11015,"temperature":1.0,"reasoning_tokens":1491,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:02:13.059691+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the free-surface velocity at several downstream positions along a chute long enough to include both undeveloped and fully developed states, for fixed material, thickness, and angle. If $v^2$ plotted against $l$ does not follow $v_\\infty^2[1 - \\exp(-l/L)]$ with $L = (c_\\infty^2/c_0^2) \\mu_r^3[(\\tan\\theta - \\mu_r)h]^{5/3} h$, or if the fitted $L$ does not scale as $[(\\tan\\theta - \\mu_r)h]^{5/3} h$, the rule is falsified. A more direct test of the weakest assumption is to image the velocity profile across the depth early in the flow: if the surface moves faster than the base, the rigid-block law is invalid.","supporting_citations":[{"cited_title":"Ancey, Dry granular flows down an inclined channel: Experimental investigations on the frictional-collisional regime, Physical Review E 65, 011304 (2001)","cited_arxiv_id":null,"evidence_quote":"Supplies the previous free-surface velocity data for spheres and sand and the stop-thickness method that the new rule must reproduce."},{"cited_title":"B¨ orzs¨ onyi and R","cited_arxiv_id":null,"evidence_quote":"Provides DEM simulation evidence that the velocity saturates as $v \\approx v_\\infty[1 - \\exp(-t/T)]$, the time-domain analogue of the spatial saturation function."},{"cited_title":"Parez and E","cited_arxiv_id":null,"evidence_quote":"Borrows the insight that steady-state transport velocity depends only on excess shear stress and an environment scale, the basis of the dimensional analysis for $v_\\infty$."},{"cited_title":"Petley, Global patterns of loss of life from landslides, Geology 40, 927 (2012)","cited_arxiv_id":null,"evidence_quote":"One of the standard flow rules (Eq. 1) that the paper tests and finds inconsistent; supplies earlier chute data and the $h/h_s$ scaling baseline."},{"cited_title":"Delannay, M","cited_arxiv_id":null,"evidence_quote":"Two-dimensional DEM simulations showing steady uniform flows exist above $\\mu_2$, supporting the null hypothesis that $h_s$ need not enter $v_\\infty$."}],"review_version":1}