pith. sign in

arxiv: 2604.14936 · v1 · submitted 2026-04-16 · ❄️ cond-mat.soft

Effect of sub-critical fluid shear flow on granular bed strength

Dong Wang , Sophie Bodek , Nicholas T. Ouellette , Mark D. Shattuck , Corey S. O'Hern This is my paper

Pith reviewed 2026-05-10 10:10 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords granular bedsub-critical flowstrength anisotropysurface grainsdislodgable fractionShields numberDEM simulationfluid shear
0
0 comments X

The pith

Granular beds strengthen in the direction of sub-critical fluid flow due to fewer dislodgable surface grains.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors use discrete element simulations to examine how fluid flows weaker than needed to entrain grains affect granular bed strength. They discover that strength increases along the flow direction and decreases against it. This directional change correlates directly with the share of surface grains that a fluid force can move in each direction. The anisotropy does not last indefinitely but decays on a timescale governed by the Shields number. Friction between particles is unnecessary for the effect but influences its strength.

Core claim

The strength of a granular bed in a particular direction after sub-critical fluid shear flow is highly correlated with the fraction of surface grains that can be dislodged by a fluid force applied in that direction. The sub-critical grain motion does not cause significant bed compaction. The anisotropic bed strength persists only over a finite time scale set by the Shields number. Inter-particle static friction is not required for the anisotropy, but its magnitude depends on the friction coefficient.

What carries the argument

The fraction of surface grains dislodgable in a given direction, serving as the control on directional bed strength following sub-critical conditioning flow.

If this is right

  • Bed strength anisotropy arises without inter-particle static friction.
  • Changes in inter-particle friction alter the degree of anisotropy.
  • Anisotropy duration is determined by the Shields number.
  • History of the bed surface fabric must be considered when predicting strength.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If flow direction varies over time in real environments, the bed may retain a memory that affects erosion patterns.
  • Monitoring the mobility of surface grains could offer a practical method to estimate bed strength in field studies.
  • Future models might incorporate this surface fraction to simulate evolving bed properties under different flows.

Load-bearing premise

The model fluid flows in the DEM simulations accurately capture sub-critical grain motion without significant bed compaction, making the dislodgable fraction the primary control on strength.

What would settle it

Measure the dislodgable surface grain fraction in different directions after sub-critical flow and verify if it quantitatively predicts the directional strengths observed in tests.

Figures

Figures reproduced from arXiv: 2604.14936 by Corey S. O'Hern, Dong Wang, Mark D. Shattuck, Nicholas T. Ouellette, Sophie Bodek.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Top panel: Sketch of the experimental setup of fluid-driven sediment beds. The blue-shaded regions contain water. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Schematic diagrams of the DEM simulations of fluid-driven granular beds in (a) 2D and (b) 3D. In (a) and (b), the left [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Mean grain velocity normalized by the fluid speed at the bed surface, [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The ratio [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Normalized time-averaged mean surface grain velocity [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. (a) The fraction of mobile surface grains [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. (a) The fraction [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Change in the average local packing fraction ∆ [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. (a) Illustration of the calculation of the critical force [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Fraction [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Fraction [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Change in the average local packing fraction ∆ [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. (a) The mean value [PITH_FULL_IMAGE:figures/full_fig_p014_13.png] view at source ↗
read the original abstract

Interactions between fluids and granular materials are prevalent on the Earth's surface. In the case of fluid flow over a sediment bed, the fluid imparts a shear stress to the granular materials. When the applied shear stress is above a critical value, the grains become entrained in the fluid flow. Prior experimental studies have shown that granular beds subjected to a sub-critical fluid flow can strengthen in the same direction as the sub-critical flow. In contrast, granular beds can become weaker in the direction opposite to the sub-critical fluid flow. To investigate the grain-scale mechanisms that control directional strengthening and weakening, we perform discrete element method (DEM) simulations of granular beds subjected to model fluid flows in two (2D) and three (3D) dimensions with varied inter-particle static friction coefficients and conditioning flow speeds. In these studies, the sub-critical grain motion does not cause significant bed compaction. Instead, we find that the strength of a granular bed in a particular direction is highly correlated with the fraction of {\it surface} grains that can be dislodged by a fluid force applied in that direction. Further, the anisotropic bed strength only persists over a finite time scale that is set by the Shields number. We also show that inter-particle static friction is not required for bed strength anisotropy, but varying the friction affects the magnitude of the anisotropy. This research enhances the grain-scale understanding of erosion of granular beds caused by fluid flows and underscores the importance of tracking the history of the fabric of the bed surface since it couples strongly to bed strength.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript uses 2D and 3D DEM simulations of granular beds subjected to sub-critical model fluid shear flows, with variations in inter-particle friction and conditioning flow speeds. It reports that directional bed strength correlates strongly with the fraction of surface grains that can be dislodged by a fluid force applied in that direction, that the resulting anisotropy decays over a finite time scale controlled by the Shields number, and that friction is not required for anisotropy but modulates its magnitude. Sub-critical motion is stated to produce no significant compaction.

Significance. If the reported correlation and time-scale dependence hold under validated fluid modeling, the work supplies a concrete grain-scale mechanism connecting surface fabric to directional erosion resistance, extending prior experimental observations of sub-critical strengthening/weakening. The parameter variations and explicit separation of compaction from fabric effects are positive features that could inform continuum sediment-transport models.

major comments (3)
  1. [Methods] Methods (simulation setup): The 'model fluid flows' used to compute dislodgement forces and sub-critical grain motion are not specified in sufficient detail (e.g., whether drag is Stokesian, whether lift or pressure-gradient terms are included, or how the flow is imposed without resolving the fluid). This is load-bearing because the central claim equates the measured dislodgable fraction to the control on strength; without explicit validation against known sub-critical hydrodynamic benchmarks or experimental grain trajectories, the fraction could be an artifact of the force model.
  2. [Results] Results (correlation analysis): The abstract and results state that strength 'is highly correlated' with the dislodgable fraction, yet no quantitative metric (Pearson coefficient, R², or regression slope with uncertainty), number of independent realizations, or error bars on the fraction itself are provided. This weakens the claim that the fraction is the primary control rather than a correlated proxy.
  3. [Results] Results (compaction check): The statement that sub-critical motion 'does not cause significant bed compaction' is central to ruling out density changes as the source of anisotropy, but no quantitative metrics (bed-height change, porosity evolution, or coordination-number shift) with statistical comparison to the initial state are referenced. These data are required to substantiate the claim.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'model fluid flows' should be replaced by a concise description of the force law or a reference to the methods section so readers can immediately assess the hydrodynamic fidelity.
  2. [Introduction] Notation: Ensure the Shields number is defined explicitly (including the grain diameter and fluid density used) at first use, and confirm consistency between 2D and 3D definitions.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed review. The comments identify areas where additional detail and quantification will strengthen the manuscript. We address each major comment below and will revise accordingly.

read point-by-point responses

  1. Referee: [Methods] Methods (simulation setup): The 'model fluid flows' used to compute dislodgement forces and sub-critical grain motion are not specified in sufficient detail (e.g., whether drag is Stokesian, whether lift or pressure-gradient terms are included, or how the flow is imposed without resolving the fluid). This is load-bearing because the central claim equates the measured dislodgable fraction to the control on strength; without explicit validation against known sub-critical hydrodynamic benchmarks or experimental grain trajectories, the fraction could be an artifact of the force model.

    Authors: We agree that the fluid force model requires more explicit specification to support the central claim. Section 2.2 of the manuscript describes the model as a linear shear profile with a local drag force on each grain, but we will expand this section in revision to include: (i) the explicit Stokesian drag formula with the coefficient calibrated to experimental terminal settling velocities, (ii) confirmation that lift and pressure-gradient forces are omitted under the sub-critical approximation, and (iii) the precise numerical implementation of the imposed flow (fixed far-field velocity gradient without full fluid resolution). We will also add a dedicated validation subsection comparing simulated incipient-motion thresholds to the Shields curve and selected experimental grain trajectories from the literature. These additions will demonstrate that the dislodgable fraction is not an artifact of the force model. revision: yes

  2. Referee: [Results] Results (correlation analysis): The abstract and results state that strength 'is highly correlated' with the dislodgable fraction, yet no quantitative metric (Pearson coefficient, R², or regression slope with uncertainty), number of independent realizations, or error bars on the fraction itself are provided. This weakens the claim that the fraction is the primary control rather than a correlated proxy.

    Authors: The referee correctly notes the absence of quantitative correlation metrics. In the revised manuscript we will add these statistics: Pearson correlation coefficients (with p-values), R² values and regression slopes with uncertainties from linear fits, the number of independent realizations per parameter combination (typically 8–12), and error bars on the dislodgable fraction computed as the standard deviation across realizations. These will be presented both in the text and in a new or augmented figure. The added metrics will allow readers to assess the strength of the correlation and the degree to which the dislodgable fraction serves as the primary control. revision: yes

  3. Referee: [Results] Results (compaction check): The statement that sub-critical motion 'does not cause significant bed compaction' is central to ruling out density changes as the source of anisotropy, but no quantitative metrics (bed-height change, porosity evolution, or coordination-number shift) with statistical comparison to the initial state are referenced. These data are required to substantiate the claim.

    Authors: We acknowledge that the compaction claim was stated qualitatively and requires quantitative support. In the revision we will insert a new subsection (or expanded paragraph) reporting: (i) time series of bed height with mean change and standard deviation, (ii) porosity profiles versus depth before and after conditioning, and (iii) average coordination number with statistical comparison (e.g., paired t-test or mean difference with 95% confidence intervals) to the initial state. These metrics will show variations below 1% and no statistically significant compaction, thereby confirming that the observed anisotropy arises from surface fabric rather than bulk density changes. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained via independent DEM measurements

full rationale

The paper's central result—that directional bed strength correlates with the fraction of surface grains dislodgable by a directional fluid force, with anisotropy decaying on a Shields-number timescale—arises from direct DEM simulations in 2D and 3D. Strength is measured via separate shear tests, while the dislodgable fraction is computed by applying model fluid forces to surface grains; these are distinct observables, not defined in terms of each other. No parameters are fitted to the target correlation, no self-citations provide load-bearing uniqueness theorems, and the model fluid flows are treated as an input assumption whose validity is external to the derivation chain. Varying friction and flow speeds further tests the correlation without reducing it to a tautology. The absence of compaction is reported as an observation, not a definitional constraint.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of DEM modeling of sub-critical flows and the assumption that surface grain dislodgability directly governs macroscopic directional strength without compaction effects.

free parameters (2)
  • inter-particle static friction coefficient
    Varied across simulations to test its effect on the magnitude of strength anisotropy.
  • conditioning flow speeds
    Varied to investigate dependence of bed strengthening on sub-critical flow intensity.
axioms (1)
  • domain assumption Discrete element method (DEM) with model fluid flows accurately represents sub-critical grain dynamics and interactions in granular beds.
    Invoked as the basis for all simulation results on grain motion and bed strength.

pith-pipeline@v0.9.0 · 5590 in / 1499 out tokens · 34378 ms · 2026-05-10T10:10:08.528975+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

35 extracted references · 35 canonical work pages

  1. [1]

    Nearest neighbor grains may also prevent test grainifrom moving a distanceσ l and cause small deviations betweenF 0 g,c,i andF g,c,i in 3D

    We observe an excellent match betweenF 0 g,c,1 andF g,c,1 for all anglesξ g. Nearest neighbor grains may also prevent test grainifrom moving a distanceσ l and cause small deviations betweenF 0 g,c,i andF g,c,i in 3D. 11 The excellent agreement betweenF g,c from MD simulations over a wide range ofµand the analytical expressions in Eqs. 9 and 10 in 2D and 3...

    work page

  2. [2]

    D. J. Jerolmack and K. E. Daniels, Viewing earth’s surface as a soft-matter landscape, Nat. Rev. Phys.1, 716 (2019)

    work page 2019

  3. [3]

    V. J. Galay, Causes of river bed degradation, Water Resour. Res.19, 1057 (1983)

    work page 1983

  4. [4]

    J. M. Buffington and D. R. Montgomery, A systematic analysis of eight decades of incipient motion studies, with special reference to gravel-bedded rivers, Water Resour. Res.33, 1993 (1997)

    work page 1993

  5. [5]

    P¨ ahtz, A

    T. P¨ ahtz, A. H. Clark, M. Valyrakis, and O. Dur´ an, The physics of sediment transport initiation, cessation, and entrainment across aeolian and fluvial environments, Rev. Geophys.58, e2019RG000679 (2020)

    work page 2020

  6. [6]

    Dey, Fluvial hydrodynamics: hydrodynamic and sediment transport phenomena, Fluvial Hydrodynamics: Hydrody- namic and Sediment Transport Phenomena , 1 (2014)

    S. Dey, Fluvial hydrodynamics: hydrodynamic and sediment transport phenomena, Fluvial Hydrodynamics: Hydrody- namic and Sediment Transport Phenomena , 1 (2014)

    work page 2014

  7. [7]

    Poesen, Soil erosion in the anthropocene: Research needs, Earth Surf

    J. Poesen, Soil erosion in the anthropocene: Research needs, Earth Surf. Process. Landf.43, 64 (2018)

    work page 2018

  8. [8]

    X. Wu, N. Bi, J. Xu, J. A. Nittrouer, Z. Yang, Y. Saito, and H. Wang, Stepwise morphological evolution of the active yellow river (huanghe) delta lobe (1976–2013): Dominant roles of riverine discharge and sediment grain size, Geomorphology292, 115 (2017)

    work page 1976

  9. [9]

    Vellinga, Beach and dune erosion during storm surges, Coast

    P. Vellinga, Beach and dune erosion during storm surges, Coast. Eng.6, 361 (1982)

    work page 1982

  10. [10]

    R. A. Feagin, D. J. Sherman, and W. E. Grant, Coastal erosion, global sea-level rise, and the loss of sand dune plant habitats, Front. Ecol. Environ.3, 359 (2005)

    work page 2005

  11. [11]

    A. A. Berhe, R. T. Barnes, J. Six, and E. Mar´ ın-Spiotta, Role of soil erosion in biogeochemical cycling of essential elements: Carbon, nitrogen, and phosphorus, Annu. Rev. Earth Planet. Sci.46, 521 (2018)

    work page 2018

  12. [12]

    R. G. Hilton and A. J. West, Mountains, erosion and the carbon cycle, Nat. Rev. Earth Environ.1, 284 (2020)

    work page 2020

  13. [13]

    D. Im, H. Kang, K.-H. Kim, and S.-U. Choi, Changes of river morphology and physical fish habitat following weir removal, Ecol. Eng.37, 883 (2011)

    work page 2011

  14. [14]

    Paphitis and M

    D. Paphitis and M. B. Collins, Sand grain threshold, in relation to bed ‘stress history’: an experimental study, Sedimen- tology52, 827 (2005)

    work page 2005

  15. [15]

    Haynes and G

    H. Haynes and G. Pender, Stress history effects on graded bed stability, J. Hydraul. Eng.133, 343 (2007)

    work page 2007

  16. [16]

    M. M. Piedra, H. Haynes, and T. B. Hoey, The spatial distribution of coarse surface grains and the stability of gravel river beds, Sedimentology59, 1014 (2012)

    work page 2012

  17. [17]

    Mao, The effects of flood history on sediment transport in gravel-bed rivers, Geomorphology322, 196 (2018)

    L. Mao, The effects of flood history on sediment transport in gravel-bed rivers, Geomorphology322, 196 (2018)

    work page 2018

  18. [18]

    C. C. Masteller, N. J. Finnegan, J. M. Turowski, E. M. Yager, and D. Rickenmann, History-dependent threshold for motion revealed by continuous bedload transport measurements in a steep mountain stream, Geophys. Res. Lett.46, 2583 (2019)

    work page 2019

  19. [19]

    C. C. Masteller, J. P. L. Johnson, D. Rickenmann, and J. M. Turowski, Modeling memory in gravel-bed rivers: a flow- history-dependent relation for evolving thresholds of motion, Earth Surface Dynamics13, 593 (2025)

    work page 2025

  20. [20]

    W. E. Dietrich, J. W. Kirchner, H. Ikeda, and F. Iseya, Sediment supply and the development of the coarse surface layer in gravel-bedded rivers, Nature340, 215 (1989)

    work page 1989

  21. [21]

    Ferdowsi, C

    B. Ferdowsi, C. P. Ortiz, M. Houssais, and D. J. Jerolmack, River-bed armouring as a granular segregation phenomenon, Nat. Commun.8, 1363 (2017)

    work page 2017

  22. [22]

    Allen and A

    B. Allen and A. Kudrolli, Granular bed consolidation, creep, and armoring under subcritical fluid flow, Phys. Rev. Fluids 3, 074305 (2018)

    work page 2018

  23. [23]

    Charru, H

    F. Charru, H. Mouilleron, and O. Eiff, Erosion and deposition of particles on a bed sheared by a viscous flow, J. Fluid Mech.519, 55–80 (2004)

    work page 2004

  24. [24]

    T. S. Majmudar and R. P. Behringer, Contact force measurements and stress-induced anisotropy in granular materials, Nature435, 1079 (2005)

    work page 2005

  25. [25]

    R. P. Behringer and B. Chakraborty, The physics of jamming for granular materials: a review, Rep. Prog. Phys.82, 012601 (2018)

    work page 2018

  26. [26]

    Y. Wang, Y. Wang, and J. Zhang, Connecting shear localization with the long-range correlated polarized stress fields in granular materials, Nat. Commun.11, 4349 (2020). 16

    work page 2020

  27. [27]

    Galanis, M

    M. Galanis, M. D. Shattuck, C. S. O’Hern, and N. T. Ouellette, Directional strengthening and weakening in hydrodynam- ically sheared granular beds, Phys. Rev. Fluids7, 013802 (2022)

    work page 2022

  28. [28]

    Bodek, D

    S. Bodek, D. Wang, M. D. Shattuck, C. S. O’Hern, and N. T. Ouellette, Anisotropic stress history effects in erodible sediment beds, J. Geophys. Res. Earth Surf.131, e2025JF008561 (2026)

    work page 2026

  29. [29]

    A. H. Clark, M. D. Shattuck, N. T. Ouellette, and C. S. O’Hern, Onset and cessation of motion in hydrodynamically sheared granular beds, Phys. Rev. E92, 042202 (2015)

    work page 2015

  30. [30]

    A. H. Clark, M. D. Shattuck, N. T. Ouellette, and C. S. O’Hern, Role of grain dynamics in determining the onset of sediment transport, Phys. Rev. Fluids2, 034305 (2017)

    work page 2017

  31. [31]

    P. A. Cundall and O. D. L. Strack, A discrete numerical model for granular assemblies, G´ eotechnique29, 47 (1979)

    work page 1979

  32. [32]

    X. An, R. Yang, R. Zou, and A. Yu, Effect of vibration condition and inter-particle frictions on the packing of uniform spheres, Powder Technol.188, 102 (2008)

    work page 2008

  33. [33]

    P. A. Gago, D. Maza, and L. A. Pugnaloni, Ergodic-nonergodic transition in tapped granular systems: The role of persistent contacts, Papers Phys.8, 080001 (2016)

    work page 2016

  34. [34]

    L. E. Silbert, Jamming of frictional spheres and random loose packing, Soft Matter6, 2918 (2010)

    work page 2010

  35. [35]

    Kudrolli, D

    A. Kudrolli, D. Scheff, and B. Allen, Critical shear rate and torque stability condition for a particle resting on a surface in a fluid flow, J. Fluid Mech.808, 397 (2016)

    work page 2016