pith. sign in

arxiv: 2604.11634 · v1 · submitted 2026-04-13 · ⚛️ physics.flu-dyn · cond-mat.soft

From Sedimentation to Suspension: Critical Strain as a Predictor of Particle Resuspension Thresholds

Mohammadreza Mahmoudian , Simon A. Rogers , Parisa Mirbod This is my paper

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

classification ⚛️ physics.flu-dyn cond-mat.soft
keywords viscous resuspensionnon-Brownian suspensionsoscillatory shearrheometryparticle suspensionsedimentationvolume fractionresuspension threshold
0
0 comments X

The pith

Strain determines when sedimented particles resuspend in dense flows

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

Dense suspensions of non-Brownian particles form sediment beds that can be resuspended by applied shear flow. The paper establishes that the transition to suspension is controlled by the total strain experienced rather than by stress levels or flow duration. This critical strain depends on the particle volume fraction through collective motions and interparticle collisions. The finding applies equally to steady shear and oscillatory shear, leading to a predictive model and state diagram for the different regimes of particle behavior.

Core claim

Strain is the key control parameter governing the transition from a sedimented bed to a fully suspended state. This strain-driven onset is mediated by effective interparticle collisions and collective particle motion. A predictive model captures the observed strain thresholds as a function of volume fraction, allowing construction of a state diagram for sedimentation, resuspension, and full suspension regimes under both steady and oscillatory conditions.

What carries the argument

Critical strain threshold that varies with particle volume fraction, determined through bulk rheometry and in situ rheo-microscopy observations.

If this is right

  • Resuspension thresholds can be predicted directly from applied strain for any given particle concentration.
  • The same thresholds hold for both constant-rate and oscillating shear flows.
  • A state diagram can now separate the conditions for settled, partially resuspended, and fully suspended states.

Where Pith is reading between the lines

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

  • If the volume-fraction dependence holds more broadly, it could reduce the need for case-by-case testing in engineering applications.
  • Connections may exist to other strain-controlled transitions in suspensions, such as the onset of yielding or jamming.

Load-bearing premise

The critical strain thresholds are determined primarily by volume fraction without significant effects from particle size distribution or fluid properties outside the tested range.

What would settle it

A set of experiments at fixed volume fraction but with particles of different diameters or in fluids of different viscosities that yield inconsistent critical strains would falsify the primary dependence on volume fraction alone.

Figures

Figures reproduced from arXiv: 2604.11634 by Mohammadreza Mahmoudian, Parisa Mirbod, Simon A. Rogers.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Schematic of the rheometer setup with side-view imaging for particle visualization. (b) Magnitude of complex viscosity [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Transient evolution of viscosity and microstructure during shear startup ( [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Transient stress evolution and strain-governed resuspension dynamics under steady shear. (a) Shear stress plotted as a function of strain [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of resuspension behavior under oscillatory and steady shear across different particle volume fractions ranging from 0.30 [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Evolution of the effective friction coefficient [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. State diagrams and microstructural snapshots in dense suspensions. (a, b) Regime maps for oscillatory and steady shear flows, [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
read the original abstract

Viscous resuspension, the process by which sedimented particles are re-entrained into a fluid under flow, is central to numerous natural and industrial systems, including environmental contaminant transport, riverbed erosion, and biogeochemical cycling. Despite its ubiquity and importance, predicting when and how resuspension occurs remains challenging, particularly under oscillatory shear, where particle interactions are nonlinear, collective, and time-dependent. Here, we examine the resuspension dynamics of dense, non-Brownian suspensions under both steady and oscillatory shear using bulk rheometry and in situ rheo-microscopy over a broad range of particle volume fractions ({\phi}= 0.30 to 0.55). We demonstrate that strain is the key control parameter governing the transition from a sedimented bed to a fully suspended state. This strain-driven onset is mediated by effective interparticle collisions and collective particle motion. We develop a predictive model that captures the observed strain thresholds as a function of volume fraction, allowing for the construction of a new state diagram delineating sedimentation, resuspension, and full suspension regimes. These findings reveal a robust, strain-controlled resuspension mechanism and establish a unified framework for predicting suspension behavior across steady and oscillatory flows, offering new tools for managing particle-laden transport in geophysical, biological, and industrial environments.

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

2 major / 2 minor

Summary. The manuscript examines viscous resuspension of dense non-Brownian suspensions (φ = 0.30–0.55) under steady and oscillatory shear via bulk rheometry and in situ rheo-microscopy. It claims that accumulated strain is the primary control parameter governing the transition from a sedimented bed to a fully suspended state, mediated by interparticle collisions and collective motion, and presents a predictive model for critical strain thresholds as a function of volume fraction that enables a state diagram separating sedimentation, resuspension, and full-suspension regimes.

Significance. If the central claim is supported by independent derivation and reproducible data, the work would offer a unified strain-based framework for resuspension prediction across flow types, with direct relevance to geophysical, industrial, and biological particle transport. The emphasis on collective dynamics and the resulting state diagram represent potentially useful organizing concepts.

major comments (2)
  1. [Abstract / Model development] Abstract and model section: the claim that the model 'captures the observed strain thresholds as a function of volume fraction' is presented without any derivation, governing equations, or fitting procedure; this leaves open whether the relation is derived from force balance or is an empirical fit to the same data used for validation, directly undermining the 'predictive' assertion.
  2. [Experimental methods] Experimental design: the reported strain(φ) relation is obtained at fixed particle radius a and fluid viscosity η (no indication that these were systematically varied); classical viscous resuspension criteria involve dimensionless groups such as ηγ̇a/(Δρga²) that do not drop out when a and η are held constant, so the φ-only model cannot be regarded as general without explicit decoupling of these parameters.
minor comments (2)
  1. [Abstract] The abstract states the volume-fraction range but provides no error bars, number of replicates, or exclusion criteria for the reported thresholds.
  2. [Abstract] Notation for strain (accumulated vs. amplitude) is not defined in the summary statement of the central claim.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which have helped us identify areas where the manuscript can be strengthened. We address each major comment below and describe the revisions that will be made.

read point-by-point responses

  1. Referee: [Abstract / Model development] Abstract and model section: the claim that the model 'captures the observed strain thresholds as a function of volume fraction' is presented without any derivation, governing equations, or fitting procedure; this leaves open whether the relation is derived from force balance or is an empirical fit to the same data used for validation, directly undermining the 'predictive' assertion.

    Authors: We agree that the model development section requires expansion to avoid ambiguity. The critical strain relation is obtained from a scaling analysis of collision-induced vertical displacements: the mean free path between particles scales as a/φ^{1/3}, and the number of collisions per unit strain is proportional to φ, yielding a functional form γ_c(φ) that increases with φ. Prefactors are calibrated to the measured thresholds. In the revised manuscript we will add a dedicated subsection containing the governing equations for the collision rate, the force-balance argument leading to the φ dependence, and the explicit fitting procedure. We will also separate the derivation from the subsequent validation against the full dataset (steady and oscillatory) to clarify the predictive aspect within the viscous regime. revision: yes

  2. Referee: [Experimental methods] Experimental design: the reported strain(φ) relation is obtained at fixed particle radius a and fluid viscosity η (no indication that these were systematically varied); classical viscous resuspension criteria involve dimensionless groups such as ηγ̇a/(Δρga²) that do not drop out when a and η are held constant, so the φ-only model cannot be regarded as general without explicit decoupling of these parameters.

    Authors: We acknowledge that a and η were held fixed to isolate the role of volume fraction over the range φ = 0.30–0.55. In the viscous limit relevant to our experiments, the accumulated strain governs the transition because the resuspension mechanism is kinematic (collision-mediated lift) rather than force-threshold based; the cited dimensionless group influences the time scale but drops out of the strain threshold itself. The revised manuscript will include a new paragraph discussing the relevant dimensionless numbers, scaling arguments for why the φ dependence is expected to persist when a and η vary (provided the flow remains viscous and non-inertial), and an explicit statement that systematic variation of a and η constitutes an important direction for future work. This will qualify the generality of the model without requiring new experiments at this stage. revision: partial

Circularity Check

1 steps flagged

Predictive model reduces to empirical fit of observed strain(φ) thresholds at fixed particle/fluid parameters

specific steps
  1. fitted input called prediction [Abstract]
    "We develop a predictive model that captures the observed strain thresholds as a function of volume fraction, allowing for the construction of a new state diagram delineating sedimentation, resuspension, and full suspension regimes."

    The thresholds are first observed experimentally (at fixed particle size and fluid properties, per the skeptic analysis and lack of parameter variation in the abstract). The model is then constructed to match those same observed thresholds versus φ. This makes the 'prediction' equivalent to a post-hoc fit of the input data rather than an independent derivation; the strain(φ) relation is forced by construction from the measurements used to validate it.

full rationale

The paper's central claim is that a 'predictive model' captures strain thresholds as a function of volume fraction φ, enabling a state diagram. However, the abstract explicitly states the model 'captures the observed strain thresholds as a function of volume fraction' without evidence that particle radius a, density contrast Δρ, or viscosity η were varied. This matches the fitted-input-called-prediction pattern: thresholds are measured experimentally at fixed a, η, Δρ, then the φ-dependence is fitted and relabeled as a general predictor. No first-principles derivation or decoupling of other dimensional groups is shown in the provided text, so the 'prediction' is statistically forced by the input data. The derivation chain is otherwise self-contained against external benchmarks but carries this partial circularity in the model claim.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the experimental observation that strain thresholds vary with volume fraction and on the assumption that a simple function of volume fraction suffices to predict them. No independent derivation or external validation is described in the abstract.

free parameters (1)
  • strain threshold function parameters
    The model is said to capture observed thresholds as a function of volume fraction, implying parameters were adjusted to match the experimental data.
axioms (1)
  • domain assumption Strain is the dominant control parameter for the sedimentation-to-suspension transition
    Invoked directly as the key finding from the rheometry and microscopy observations.

pith-pipeline@v0.9.0 · 5540 in / 1343 out tokens · 69192 ms · 2026-05-10T15:40:39.136457+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

4 extracted references · 4 canonical work pages

  1. [1]

    A gravitational force is generated owing to the density difference between the par- ticles and the silicone oil, ensuring sedimentation in the sus- pension system

    pur- chased from Sigma-Aldrich, which behaves as a Newtonian fluid with a viscosity ofη m= 0.02 Pa·s. A gravitational force is generated owing to the density difference between the par- ticles and the silicone oil, ensuring sedimentation in the sus- pension system. To prepare the suspension, the particles were dispersed in the Newtonian silicone oil using...

    work page 1920

  2. [2]

    Viscous resuspension,

    State diagrams and microstructural snapshots in dense suspensions. (a, b) Regime maps for oscillatory and steady shear flows, respectively, were constructed over a wide range of particle volume fractions (φ=0.30 to 0.55). The diagrams delineate four distinct flow regimes: sedimentation (green), resuspension (purple), fully suspended flow (orange), and jam...

    work page 1979

  3. [3]

    Transport layer structure in intense bed- load,

    bed-load transport by laminar shearing flows,” Journal of Fluid Mechanics 636, 295–319 (2009). 33H. Capart and L. Fraccarollo, “Transport layer structure in intense bed- load,” Geophysical Research Letters38(2011). 34A. G. Kidanemariam and M. Uhlmann, “Interface-resolved direct numerical simulation of the erosion of a sediment bed sheared by laminar chann...

    work page 2009

  4. [4]

    Causes and effects of noise in landscape dynamics,

    con- centration dependence,” Rheologica acta38, 2–13 (1999). 42D. J. Jerolmack, “Causes and effects of noise in landscape dynamics,” Eos, Transactions American Geophysical Union92, 385–386 (2011). 43M. Mahmoudian, F. Goharpey, M. Behzadnasab, and Z. Daneshfar, “Shear thickening and hysteresis in dense suspensions: The effect of particle shape,” Journal of...

    work page 1999