Pith. sign in

REVIEW 3 major objections 5 minor 35 references

Hindered stokesian settling of discs and rods

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper claims that the hindered settling of discs and rods collapses onto the sphere curve after a horizontal shift equal to the particle volume relative to a sphere, so mean settling speed is set by volume, not orientation.

desk verdict New flat-particle settling data and a plausible volume-based collapse, but the exponents rest on three points and a possibly biased Stokes normalization. read the letter →

arxiv 2411.14363 v2 pith:U32KJRUL submitted 2024-11-21 cond-mat.soft

Yating Zhang , Narayanan Menon This is my paper · ORCID
classification cond-mat.soft
keywords hinderedsettlingStokesflowanisotropicparticlesdiscsrodssedimentationRichardson-Zakiaspectratio
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper measures how suspensions of discs and rods settle in viscous fluid at low Reynolds number. It claims that all six particle shapes follow the same hindered settling function as spheres once each data set is shifted horizontally by a factor S(A), and that this factor is just the volume of a sphere of the same largest dimension divided by the particle volume. The result would matter because it turns a many-body hydrodynamic problem into a volume-ratio correction: for any nonpolar axisymmetric particle, the mean settling curve could be predicted from sphere data without accounting for orientation dynamics. The paper also argues that this supports backflow, proportional to particle volume, as the dominant hindering mechanism.

What carries the argument

The argument is carried by three objects: the hindered settling function H=U/U_S, the normalized interparticle separation a* = a/L = k(A) $φ^{{-1/3}}$, and the horizontal shift factor S(A) found by total least squares collapse of each shape's H(a*) onto the sphere Richardson-Zaki curve (1-φ)^4.5. S(A) absorbs all shape dependence, and its measured $A^{{-1}}$ and $A^{2}$ scalings match the ratio of sphere volume to particle volume. Physically, the paper invokes Batchelor's dilute-limit calculation, which attributes -5.5 of the -6.55 φ slope to upward backflow proportional to particle volume; the successful volume-only collapse suggests this backflow term dominates even at finite φ.

What would settle it

Measure the orientation distribution of rods and discs during settling and recompute H using a single-particle Stokes velocity matched to that distribution; if the required shift S(A) changes with the orientation state, or if a cylinder with aspect ratio one does not land at the predicted S=2/3, the volume-only collapse is an artifact of the normalization.

Watch

Extended reading notes

Core claim

The central discovery is that, over aspect ratios from thin flat particles to long rods, the hindered settling function H(φ)=U(φ)/U_S is the same universal curve as for spheres once plotted against normalized interparticle separation and shifted by a shape-dependent factor S(A). For flat particles S(A) follows $A^{{-1}}$ and for rods $A^{2}$, with prefactor 2/3, which is exactly the ratio of the volume of a sphere of diameter L to the volume of the particle. The collapse includes previously published fiber data and holds up to the semi-dilute regime where the mean separation is comparable to the particle size. The authors conclude that the dominant contribution to hindered settling is volume-proportional backflow, and that orientational degrees of freedom do not change the mean settling velocity.

Load-bearing premise

The load-bearing premise is that the single-particle Stokes velocity used to normalize each suspension is the correct reference for the particles' orientation state while settling, since rods were normalized by the vertically oriented velocity and flat particles by an average over initial orientations.

Editorial extensions

If this is right

  • The mean hindered settling of any nonpolar axisymmetric particle, such as a clay platelet, a paper fiber, or a red blood cell, can be predicted from the sphere master curve using only the ratio of the sphere volume to the particle volume.
  • Existing empirical sphere fits, such as Richardson-Zaki with n≈4.5, can serve as the universal reference curve for shaped particles after the S(A) shift.
  • At equal normalized interparticle separation, flat particles and rods hinder settling less than spheres, with hindering beginning only when the typical separation falls below one largest dimension.
  • The dilute-limit backflow contribution identified by Batchelor for spheres remains the dominant hindering mechanism in the semi-dilute regime for anisotropic shapes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the volume-only rule holds more generally, it provides a cheap predictive protocol: calibrate once with spheres and compute particle volume to estimate settling in any new suspension.
  • A natural theoretical target would be a Batchelor-style dilute-limit calculation for anisotropic particles; the data suggest the backflow coefficient is shape-independent, while only the volume normalization changes.
  • The rule may extend to non-axisymmetric shapes such as rectangular platelets or ellipsoidal grains, but the present data only directly support effectively axisymmetric particles.
  • Orientation may still matter for quantities the paper did not measure, such as the sharpness of the settling interface and velocity fluctuations, even if it drops out of the mean settling speed.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This manuscript presents sedimentation experiments on monodisperse suspensions of three flat particle types (hexagons with aspect ratio A = 0.15, 0.07 and a disc with A = 0.05) and three rod types (A = 3, 10, 19), all in the Stokes regime. The authors extract the mean settling velocity from interface tracking, with PIV used at low volume fraction, and normalize by the measured single-particle Stokes velocity to form the hindered settling function H = U/U_S. They plot H against the normalized interparticle separation a* = a/L and compare with the literature sphere Richardson–Zaki curve. They find that the shaped-particle data collapse onto the sphere curve after a horizontal shift S(A), and that S(A) is proportional to A^-1 for flat particles and A^2 for rods, which they identify with the ratio of the sphere volume to the particle volume at fixed largest dimension. From this they argue that backflow proportional to particle volume is the dominant mechanism of hindered settling for nonpolar axisymmetric particles.

Significance. If the proposed volume-ratio collapse is correct, the paper provides a simple, useful rule for estimating hindered settling of rods and platelets and extends the classical sphere-based empirical framework to anisotropic shapes. The experimental work has clear strengths: the particles are well characterized and monodisperse, Reynolds numbers are very small, the settling interface and PIV measurements provide complementary determinations of U, finite-size checks are reported, and the comparison target is an independently compiled sphere Richardson–Zaki curve, so the collapse is not internally circular. The principal risk is that the fitted shift factors S(A), which carry the entire volume-scaling claim, are sensitive to the normalization of H by a single-particle Stokes velocity measured in a specific orientation state and to the conversion from volume fraction to interparticle separation. These systematic uncertainties are acknowledged in the End Matter but not quantified. Because the multi-shape collapse is the central result, the scaling law is not yet established at the level claimed.

major comments (3)
  1. [Normalization paragraph after Fig. 2; Table I; End Matter] The central claim depends on the normalization H = U/U_S, but U_S is measured for vertically oriented rods and as an average over initial orientations for the flat particles. In a settling semi-dilute suspension, steric and hydrodynamic interactions constrain particle rotations, so the orientation distribution may differ from the one used in the single-particle measurement. Since the settling velocity of an anisotropic particle depends on orientation, this multiplies H by a factor f(phi,A) that is not necessarily constant. The shift factor S(A) is then obtained by minimizing the deviation of the shifted H(a*) from the sphere R-Z curve, and a vertical rescaling of a data set can be partly absorbed as a different horizontal shift. The End Matter notes that an error in U_S is a systematic error in H for a given particle type, but it does not quantify the size, sign, or phi-dependence of the effect. The authors should measure or tightly bound the orientation distribution in the settling suspension (the particles are visible in the images) and report the range of single-particle settling velocities over orientations. Without this, the fitted S(A), and hence the A^-1 and A^2 scaling, may be biased.
  2. [Fig. 3 and End Matter] The horizontal shift S(A) is not a directly measured quantity: it is the one-parameter horizontal translation that best fits each data set to the sphere R-Z curve. This procedure can absorb systematic errors in the conversion from volume fraction to interparticle separation, a* = k(A) phi^{-1/3}. The factor k(A) is taken from the literature, and an error in k(A) or in phi (both acknowledged as systematic in the End Matter) is largely indistinguishable from a shift in S. The paper reports no residuals, no confidence intervals for S(A), and no propagation of the stated uncertainties in U, U_S, phi, and k(A). It also includes literature data [22,29] without stating whether the same definitions of L, a*, and U_S were used. To make the collapse convincing, the authors should tabulate S(A) with uncertainties, show residual plots of the shifted data against the R-Z curve, and demonstrate that plausible variations in k(A) and phi do not change the fitted power-law exponents.
  3. [Fig. 3(c); Table I] The flat-particle branch of the volume-ratio law rests on only three aspect ratios (A = 0.15, 0.07, 0.05), spanning a factor of about three in A, and two of these three points are hexagons rather than discs. For a regular hexagonal prism with circumdiameter L, the geometric volume is (3*sqrt(3)/8) A L^3, whereas the disc volume used in the text is (pi/4) A L^3; at fixed A the sphere-to-particle volume ratio differs by about 20%. If the scaling is truly volume-only, the hexagon points should be placed using their actual volume, not the disc formula, and the statement that the flat shapes are 'effectively axisymmetric' should be checked against the single-particle Stokes velocity. Moreover, Fig. 3(c) shows no error bars on S(A), so the A^-1 slope for the flat branch is not yet supported with stated precision. The authors should report fitted prefactors and confidence intervals, and ideally add at least one more disc aspect ratio to separate a volume effect from a shape effect.
minor comments (5)
  1. [Table I] Table I contains typographical spacing errors in the Reynolds number entries (e.g., '6 .33 × 10−5' and '4 .96 × 10−4'); these should be corrected.
  2. [Fig. 3(c)] The power-law fits should be stated with numerical prefactors and uncertainties; the text says the prefactors are 'of order unity' without giving values, and the abstract's phrase 'exactly the ratio' is stronger than the data support.
  3. [Fig. 3(b)] The text says the collapse is achieved for 'all six of our particle shapes' but also includes two literature data sets; clarify that the comparison in Fig. 3(b) involves eight sets and state the aspect ratios and normalization used for Refs. [22,29].
  4. [Abstract and conclusion] The interpretive statement that backflow is the 'dominant contribution' and that the effect 'emerges from terms simply proportional to volume' goes beyond what an empirical horizontal-shift collapse can demonstrate; it would be better framed as a consequence of the volume-ratio scaling rather than as an independently established mechanism.
  5. [Fig. 1 and Table I] The paper should clarify the definition of 'effectively axisymmetric' for hexagons and report how the hexagon's single-particle Stokes velocity compares with a disc of the same circumdiameter and volume.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the collapse is tested against an external sphere benchmark and the volume-ratio scaling is an empirical fit, not a by-construction reduction.

full rationale

The central claim is that the hindered-settling curves for discs and rods collapse onto the sphere Richardson-Zaki curve after a horizontal shift S(A), and that S(A) follows the sphere-to-shape volume ratio. This is not circular. The sphere curve is an external literature benchmark (Brzinski & Durian; Richardson & Zaki), and S(A) is obtained by a total-least-squares fit that is permitted to fail; the paper itself hedges with 'To the extent that these eight data sets collapse within experimental error.' The volume-ratio scaling of S(A) is read off the fitted S(A) values in Fig. 3(c) and compared to independently computed geometric volume ratios; it was not an input to the fit. Thus the 'prediction' for other nonpolar axisymmetric shapes is an extrapolation of an empirical scaling law, not a parameter forced by construction. The End Matter statement that 'an error in Us is a systematic error in H for a given particle type' identifies a legitimate systematic-error concern in the normalization of flat-particle data (orientation-averaged US), but this is an accuracy issue, not circularity. Self-citations (Refs. 17 and 19) are background physics on disk pairs and lattices and are not load-bearing for the central argument. No equation is shown to be equal to another by definition, and no fitted parameter is renamed as an independent prediction.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the geometric normalization a* = k(A) phi^{-1/3}, the external sphere R-Z reference curve, the assumption of Stokes/non-Brownian conditions, and an interpretive backflow-dominance hypothesis. The fitted shift factors S(A) are free parameters; their volume-ratio scaling is inferred from the same fits, with prefactors left unconstrained.

free parameters (2)
  • S(A) shift factor for each shape = Six values (plus two from literature), exact numbers not reported
    Chosen by total least squares to collapse each shape's H(a*) data onto the Richardson-Zaki sphere curve with n=4.5. These fitted shifts are the basis for the claimed volume-ratio scaling.
  • Power-law prefactors for S(A) versus A = Not reported, described only as 'of order unity'
    The scaling S ~ A^-1 for discs and S ~ A^2 for rods is inferred from the fitted S(A) values. The prefactor is not predicted by the volume-ratio argument, only bounded as O(1), so it is an effective free parameter in the predictive rule.
assumptions (4)
  • domain assumption The suspension is in the Stokes regime (Re << 1) and non-Brownian (Pe >> 1)
    Ensures the comparison to sphere hindered-settling data is valid; stated in the introduction and quantified in Table I.
  • domain assumption Mean interparticle separation follows a* = k(A) phi^{-1/3} with k(A) from the literature for random distributions
    Used to convert measured volume fraction to the normalized separation plotted in Fig. 3(a); cited to Refs. [30-34], but specific k(A) values are not given in the main text.
  • domain assumption The Richardson-Zaki form with n=4.5 is a valid representation of sphere hindered settling in this regime
    Used as the reference curve for the collapse; taken from the Brzinski-Durian compilation of sphere sedimentation data.
  • ad hoc to paper The dominant hindering contribution is backflow, proportional to particle volume
    Inferred from the observed collapse, not derived from a Batchelor-style calculation for anisotropic particles. The paper explicitly notes no such calculation exists.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Hindered stokesian settling of discs and rods." pith.science (2026). https://pith.science/paper/U32KJRUL

@misc{pith2026241114363,
  author       = {Pith},
  title        = {Pith review of: Hindered stokesian settling of discs and rods},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/U32KJRUL}},
  note         = {Machine review of arXiv:2411.14363}
}
read the original abstract

We report measurements of the mean settling velocities for suspensions of discs and rods in the stokes regime for a number of particle aspect ratios. All these shapes display ''hindered settling'', namely, a decrease in settling speed as the solid volume fraction is increased. A comparison of our data to spheres reveals that discs and rods show less hindering than spheres at the same relative interparticle separation. The data for all six of our particle shapes may be scaled to collapse on that of spheres, with a scaling factor that depends only on the volume of the particle relative to a sphere. Despite the orientational degrees of freedom available with nonspherical particles, it thus appears that the dominant contribution to the hindered settling emerges from terms that are simply proportional to the volume of the sedimenting particles, enabling prediction of hindered settling of other nonpolar axisymmetric shapes.

Figures

Figures reproduced from arXiv: 2411.14363 by the authors.

Figure 1
Figure 1. FIG. 1: Data for discs and rods are compared to those of spheres held at the same normalized interparticle sepa [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: (a) Sample of [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: (a) Hindered settling function [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

35 extracted references · 34 canonical work pages

  1. [1]

    P. M. Chaikin, Thermodynamics and hydrodynamics of hard spheres: The role of gravity, in Soft and fragile matter: nonequilibrium dynamics, metastability and flow , edited by M. E. Cates and M. R. Evans (Institute of Physics Publishing, London, 2000) pp. 315–348

    work page 2000

  2. [2]

    R. E. Caflisch and J. H. Luke, Variance in the sedimen- tation speed of a suspension, Physics of Fluids 28, 759 (1985)

    work page 1985

  3. [3]

    Nicolai, B

    H. Nicolai, B. Herzhaft, E. Hinch, L. Oger, and E. Guazzelli, Particle velocity fluctuations and hy- 6 drodynamic self-diffusion of sedimenting non-brownian spheres, Physics of Fluids 7, 12 (1995)

    work page 1995

  4. [4]

    J.-Z. Xue, E. Herbolzheimer, M. Rutgers, W. Russel, and P. Chaikin, Diffusion, dispersion, and settling of hard spheres, Physical review letters 69, 1715 (1992)

    work page 1992

  5. [5]

    P. J. Mucha, S.-Y. Tee, D. A. Weitz, B. I. Shraiman, and M. P. Brenner, A model for velocity fluctuations in sedimentation, Journal of fluid mechanics501, 71 (2004)

    work page 2004

  6. [6]

    R. H. Davis and A. Acrivos, Sedimentation of noncol- loidal particles at low reynolds numbers, Annual Review of Fluid Mechanics 17, 91 (1985)

    work page 1985

  7. [7]

    Ramaswamy, Issues in the statistical mechanics of steady sedimentation, Advances in Physics 50, 297 (2001)

    S. Ramaswamy, Issues in the statistical mechanics of steady sedimentation, Advances in Physics 50, 297 (2001)

    work page 2001

  8. [8]

    Brzinski III and D

    T. Brzinski III and D. Durian, Observation of two branches in the hindered settling function at low reynolds number, Physical Review Fluids 3, 124303 (2018)

    work page 2018

Show all 35 references
  1. [9]

    Richardson and W

    J. Richardson and W. Zaki, Sedimentation and fluidisa- tion: Part i, Chemical Engineering Research and Design 75, S82 (1997)

    1997

  2. [10]

    Batchelor, Sedimentation in a dilute dispersion of spheres, Journal of fluid mechanics 52, 245 (1972)

    G. Batchelor, Sedimentation in a dilute dispersion of spheres, Journal of fluid mechanics 52, 245 (1972)

    1972

  3. [11]

    G. I. Taylor, Motion of axisymmetric bodies in viscous fluids, in Problems of Hydrodynamics and Continuum Mechanics, edited by L. I. Sedov (SIAM, Philadelphia,

  4. [12]

    Batchelor, Slender-body theory for particles of arbi- trary cross-section in stokes flow, Journal of Fluid Me- chanics 44, 419 (1970)

    G. Batchelor, Slender-body theory for particles of arbi- trary cross-section in stokes flow, Journal of Fluid Me- chanics 44, 419 (1970)

    1970

  5. [13]

    M. B. Mackaplow, A study of the transport properties and sedimentation characteristics of fiber suspensions , Ph.D. thesis, Stanford University, Stanford, CA (1996)

    1996

  6. [14]

    M. B. Mackaplow and E. S. Shaqfeh, A numerical study of the sedimentation of fibre suspensions, Journal of Fluid Mechanics 376, 149 (1998)

    1998

  7. [15]

    Happel and H

    J. Happel and H. Brenner, Low Reynolds number hydro- dynamics: with special applications to particulate media , Vol. 1 (Springer Science & Business Media, 2012)

    2012

  8. [16]

    Happel and R

    J. Happel and R. Pfeffer, The motion of two spheres fol- lowing each other in a viscous fluid, AIChE Journal 6, 129 (1960)

    1960

  9. [17]

    Chajwa, N

    R. Chajwa, N. Menon, and S. Ramaswamy, Kepler or- bits in pairs of disks settling in a viscous fluid, Physical Review Letters 122, 224501 (2019)

    2019

  10. [18]

    S. Jung, S. Spagnolie, K. Parikh, M. Shelley, and A.- K. Tornberg, Periodic sedimentation in a stokesian fluid, Physical Review E 74, 035302 (2006)

    2006

  11. [19]

    Chajwa, N

    R. Chajwa, N. Menon, S. Ramaswamy, and R. Govin- darajan, Waves, algebraic growth, and clumping in sed- imenting disk arrays, Physical Review X 10, 041016 (2020)

    2020

  12. [20]

    J. M. Crowley, Viscosity-induced instability of a one- dimensional lattice of falling spheres, Journal of Fluid Mechanics 45, 151 (1971)

    1971

  13. [21]

    Herzhaft, ´E

    B. Herzhaft, ´E. Guazzelli, M. B. Mackaplow, and E. S. Shaqfeh, Experimental investigation of the sedimentation of a dilute fiber suspension, Physical review letters 77, 290 (1996)

    1996

  14. [22]

    Herzhaft and ´E

    B. Herzhaft and ´E. Guazzelli, Experimental study of the sedimentation of dilute and semi-dilute suspensions of fibres, Journal of Fluid Mechanics 384, 133 (1999)

    1999

  15. [23]

    Guazzelli and J

    ´E. Guazzelli and J. Hinch, Fluctuations and instability in sedimentation, Annual review of fluid mechanics 43, 97 (2011)

    2011

  16. [24]

    Metzger, E

    B. Metzger, E. Guazzelli, and J. E. Butler, Large-scale streamers in the sedimentation of a dilute fiber suspen- sion, Physical review letters 95, 164506 (2005)

    2005

  17. [25]

    Metzger, J

    B. Metzger, J. E. Butler, and ´E. Guazzelli, Experimen- tal investigation of the instability of a sedimenting sus- pension of fibres, Journal of Fluid Mechanics 575, 307 (2007)

    2007

  18. [26]

    See Supplemental Material at URL-will-be-inserted by- publisher for more details of data, methods and analysis

  19. [27]

    G. H. Golub, Some modified matrix eigenvalue problems, SIAM review 15, 318 (1973)

    1973

  20. [28]

    G. H. Golub and C. F. Van Loan, An analysis of the total least squares problem, SIAM journal on numerical analysis 17, 883 (1980)

    1980

  21. [29]

    M. A. Turney, M. K. Cheung, R. L. Powell, and M. J. McCarthy, Hindered settling of rod-like particles mea- sured with magnetic resonance imaging, AIChE Journal 41, 251 (1995)

    1995

  22. [30]

    Hertz, ¨Uber den gegenseitigen durchschnittlichen ab- stand von punkten, die mit bekannter mittlerer dichte im raume angeordnet sind, Mathematische Annalen 67, 387 (1909)

    P. Hertz, ¨Uber den gegenseitigen durchschnittlichen ab- stand von punkten, die mit bekannter mittlerer dichte im raume angeordnet sind, Mathematische Annalen 67, 387 (1909)

    1909

  23. [31]

    Chandrasekhar, Stochastic problems in physics and astronomy, Reviews of modern physics 15, 1 (1943)

    S. Chandrasekhar, Stochastic problems in physics and astronomy, Reviews of modern physics 15, 1 (1943)

    1943

  24. [32]

    G. J. Vroege and H. N. Lekkerkerker, Phase transitions in lyotropic colloidal and polymer liquid crystals, Reports on Progress in Physics 55, 1241 (1992)

    1992

  25. [33]

    Eppenga and D

    R. Eppenga and D. Frenkel, Monte carlo study of the isotropic and nematic phases of infinitely thin hard platelets, Molecular physics 52, 1303 (1984)

    1984

  26. [34]

    Veerman and D

    J. Veerman and D. Frenkel, Phase behavior of disklike hard-core mesogens, Physical Review A 45, 5632 (1992)

    1992

  27. [35]

    Torquato, Nearest-neighbor statistics for packings of hard spheres and disks, Physical Review E 51, 3170 (1995)

    S. Torquato, Nearest-neighbor statistics for packings of hard spheres and disks, Physical Review E 51, 3170 (1995)

    1995

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.