Pith. sign in

REVIEW 3 major objections 5 minor 108 references

A particle simulation predicts how dense rod suspensions thicken with packing and aspect ratio, and how rods align under shear.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.5

2026-07-31 14:42 UTC pith:TTCOY5PI

load-bearing objection Useful DEM tool paper for dense rod suspensions: real algorithmic fixes and systematic η_s(φ,A)/alignment maps, with the expected hydro-closure caveats. the 3 major comments →

arxiv 2607.28206 v1 pith:TTCOY5PI submitted 2026-07-30 cond-mat.soft physics.flu-dyn

Rheology of dense suspensions of granular spherocylinders by particle-based simulation

Alex Dixon , John Hone , Gavin Melaugh , Christopher Ness This is my paper
classification cond-mat.soft physics.flu-dyn PACS 47.57.E-47.57.Qk83.80.Hj45.70.-n
keywords dense suspensionsspherocylindersgranular rodsrheologyparticle-based simulationlubrication forcesorientational ordershear flow
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Dense suspensions of rod-shaped grains matter in rivers, magma, crystallisers and paper making, but their high-packing rheology is hard to measure and incomplete in theory. This paper builds a discrete-element model of hard, non-Brownian spherocylinders under simple shear that includes frictional contacts, blended short-range lubrication, and orientation-dependent drag and torque. The simulation shows a start-up viscosity spike that relaxes into a steady viscosity rising with solids fraction and aspect ratio, while collective alignment along the flow grows until an aspect-ratio-dependent packing where order softens. Those trends match the sparse experiments on granular rods and give a controllable tool for mapping jamming, normal stresses and inhomogeneous flow in this class of materials.

Core claim

For aspect ratios from 1 to 20 the model predicts a transient viscosity spike at shear start-up that gives way to rate-independent steady viscosities increasing systematically with volume fraction and aspect ratio, while the orientational order parameter rises with packing up to an aspect-ratio-dependent critical point and then declines; the same trends corroborate limited experimental rod-suspension rheology.

What carries the argument

A DEM trajectory integrator for spherocylinders that couples Hookean-Coulomb contacts to a discontinuity-free, weighted blend of shaft-shaft, end-shaft and end-end lubrication forces, plus prolate-spheroid single-body hydrodynamics, with a dynamic timestep set by the shortest contact or lubrication timescale.

Load-bearing premise

Hydrodynamic drag, lift, torque and near-field lubrication on the rods are taken from prolate-spheroid formulas and a manually weighted blend of existing lubrication pieces rather than a fully resolved spherocylinder solution.

What would settle it

Measure steady shear viscosity and flow-alignment order versus packing for well-characterised hard rods of aspect ratio near 5–20 and compare the measured η_s(φ) divergence and the packing at which order peaks against the simulation curves in Figures 4 and 6.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Steady viscosity of dense granular rod suspensions is rate-independent and rises with both packing and aspect ratio, with jamming packing falling as rods lengthen.
  • An initial random orientation produces a higher start-up viscosity than the aligned steady state, so unsheared and sheared materials can jam at different packings.
  • Flow alignment strengthens with packing only up to an aspect-ratio-dependent threshold, after which frustration reduces order and inclination angle.
  • The same framework can be used to map friction, polydispersity, normal stresses and pressure-imposed or inhomogeneous flows for rod suspensions.

Where Pith is reading between the lines

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

  • If the hydrodynamic approximations hold, the model supplies the missing dense-regime branch of the Krieger–Dougherty or viscous-number picture for rods, allowing engineering estimates of φ_m(A) without new experiments for every aspect ratio.
  • The start-up spike and later drop in order suggest that process history (how a slurry is first sheared) can control whether a log-jam or crystalliser slurry locks or flows at the same solids fraction.
  • Extending the dynamic-timestep and blended-lubrication scheme to flexible or adhesive rods would open direct comparison with pulp fibres and whisker crystallisation, where bending and sticking are first-order.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The manuscript presents a DEM-style particle-based model for dense, non-Brownian suspensions of frictional spherocylinders in simple shear. Particles evolve under short-range Hookean–Coulomb contacts, a blended multi-mode lubrication force (shaft–shaft, end–shaft, end–end), and single-body drag/lift/Jeffery torques taken from prolate-spheroid analytics, integrated with a configuration-dependent dynamic timestep. For aspect ratios A = 1–20 and volume fractions from dilute up toward jamming, the model yields a start-up viscosity spike that relaxes to a rate-independent steady η_s that rises with φ and A, together with flow-aligned order S(φ) that increases up to an A-dependent threshold and then declines. Force decomposition, stress- versus rate-control equivalence, and low-φ Jeffery-period checks are reported. The authors position the tool as filling the experimental gap for dense granular-rod rheology and as qualitatively consistent with the limited available data.

Significance. Dense suspensions of granular rods are common in natural and industrial flows, yet controlled rheology and microstructure data remain sparse above semi-dilute φ, especially at large A. A tractable particle-based model that retains frictional contacts and near-field hydrodynamics—analogous to the sphere-suspension DEM tools that clarified Boyer-type rheology—is therefore of clear value. The technical contributions (discontinuity-free lubrication blending; dynamic min(t_c, t_lub,F, t_lub,T) timestepping; bidisperse spherocylinder contact detection) are concrete and reusable. If the predicted η_s(φ,A), start-up spike, and non-monotonic S(φ) prove robust to the hydrodynamic closures, the work supplies a practical platform for mapping φ_m(A), normal stresses, and inhomogeneous flows that experiments currently struggle to access.

major comments (3)
  1. [§II.C.3, Eqs. (25)–(26); Fig. 4–6] §II.C.3, Eqs. (25)–(26): The lubrication law is a weighted blend of Yamane SS and Butler–Shaqfeh ES/EE terms with manually chosen correctors a=0.1, b=1.0, c=5.0 fixed after inspecting primitive pair trajectories. The text asserts that weight choice has “minimal impact on the rheological predictions” but provides no sensitivity of the actual Fig. 4–6 curves (η_s(φ,A), start-up spike, or S(φ) peak location). At intermediate φ, Fig. 4(d) shows lubrication still comparable to contacts; for large A, shaft–end transitions are frequent. Because the start-up spike and the alignment critical point are not in the contact-dominated limit, a short sensitivity panel (or appendix) varying a,b,c (and h_min) on the reported observables is needed to secure the claim that the trends are material rather than closure artefacts.
  2. [Abstract; §III.A; §IV] Abstract and §III–IV claim that the model “corroborates the limited experimental rheology data” (citing e.g. Tapia et al. 2017). The manuscript shows only qualitative consistency (η_s rising with φ and A; φ_m decreasing with A). No direct overlay of simulated η_s(φ) against published fibre-suspension data at matched A, nor a quantitative statement of agreement/disagreement on φ_m or the magnitude of η_s, is given. A single comparative panel or table against the densest available non-colloidal rod data would make the corroboration claim falsifiable and load-bearing rather than rhetorical.
  3. [§II.B, Eqs. (3)–(10); Fig. 6(a) inset] §II.B, Eqs. (3)–(10): Single-body drag, lift and Jeffery torque are taken from analytical prolate-spheroid formulae and applied unchanged to spherocylinders. The approximation is acknowledged, but there is no validation (even for an isolated rod’s orbit period or drag coefficient versus known spherocylinder results) before the many-body runs. Because low-φ order oscillations and the fluid stress baseline feed the interpretation of the start-up transient, a brief isolated-particle check against Jeffery/spherocylinder benchmarks would strengthen confidence that orientation-dependent driving is not systematically biased.
minor comments (5)
  1. [§II.D.1, Eq. (28)] §II.D.1, Eq. (28): The fluid contribution σ_F = ηγ̇(1 + 5/2 φ + φ A²/16) is a dilute isotropic expression. The authors note it becomes negligible at high φ; stating explicitly the φ range where σ_F / σ_P drops below a few percent (or omitting it above that threshold) would avoid confusion.
  2. [Table II; §II.D.3] Table II and Fig. 4(b): N is increased with A and φ to keep the box above the minimum interaction length, but finite-size checks (e.g. doubling N at one high-A, high-φ point) are not reported. A sentence on residual box-size sensitivity would help.
  3. [Fig. 4(b); §III.A] Fig. 4(b): Power-law fits for A > 5 are described as “poorer” and φ_m values are only approximate. Either show the fits or replace the verbal φ_m estimates with a clearer operational definition (e.g. φ where η_s exceeds a threshold).
  4. [§II] Notation: both A (aspect ratio) and A in the inertia tensor components appear; ˙γ* is defined late. A short symbol table or earlier definition of ˙γ* would aid readability.
  5. [Throughout] Typos / style: “SIMULA TION”, “F orces”, “RESUL TS” (spaced capitals in headings); “e.g.” spacing; arXiv date “July 31, 2026” is fine for the preprint but should be checked at production.

Circularity Check

0 steps flagged

No circularity: forward DEM outputs from fixed microscopic closures, not fits or self-definitional reductions of η_s(φ,A) or S.

full rationale

The paper’s central claims are numerical outputs of Newtonian particle dynamics (Eqs. 1–2) under fixed contact, lubrication, and single-body hydrodynamic closures. Viscosity is computed from interaction and fluid stress (Eqs. 27–28) after integrating trajectories; order S is the largest eigenvalue of the Q-tensor (Eq. 29). Microscopic parameters (μ, e_n, e_t, h_min, h_max, and lubrication blend correctors a=0.1, b=1.0, c=5.0) are set a priori or by smoothness of primitive pair motions, then held fixed; the text does not fit them to the reported steady η_s(φ,A) or S(φ) curves. Rate-independence checks and stress-control feedback are consistency tests, not inverted predictions. Citations to the group’s sphere-suspension models supply the A=1 baseline and methodological precedent only; they do not force the rod viscosity or alignment results. Approximations (spheroid drag/Jeffery torques on spherocylinders; manually weighted lubrication) are modeling assumptions that may affect correctness, but they are not circular reductions of the claimed outputs to their inputs. No self-definitional loop, fitted-input-as-prediction, or load-bearing uniqueness import is present.

Axiom & Free-Parameter Ledger

8 free parameters · 8 axioms · 2 invented entities

The central claims rest on standard granular-suspension modeling choices plus several hand-set numerical parameters and continuum approximations for non-spherical hydrodynamics. No new physical entities are postulated; the load is carried by DEM contact laws, short-range lubrication closures, and prolate-spheroid single-body hydrodynamics under the assumption that inertia is negligible and long-range hydrodynamics can be dropped in the dense regime.

free parameters (8)
  • lubrication blend correctors a,b,c = a=0.1, b=1.0, c=5.0
    Manual weights (a=0.1, b=1.0, c=5.0) rescale SS/ES/EE contributions after inverse-distance weighting to smooth force transitions; chosen by inspecting primitive pair motions, not derived.
  • h_min (minimum surface gap fraction) = 0.025
    Caps lubrication divergence and sets asperity-like cutoff; directly affects lubrication stress plateau at high φ.
  • h_max (lubrication cutoff) = 0.9
    Maximum separation for computing lubrication; affects neighbor list and near-field stress.
  • friction coefficient μ = 0.5
    Coulomb friction in tangential contact; authors note qualitative alignment is μ-insensitive but defer systematic study; sets contact stress partition.
  • restitution e_n, e_t = 0.5, 0.5
    Enter damping and tangential stiffness via contact-time formulas; modest literature-typical values, not material-specific.
  • normal stiffness k_n = 1.5e6 (sim units)
    Sets hard-particle limit via ˙γ√(ρ_P R^3/k_n)<10^{-4}; large value 1.5×10^6 in simulation units.
  • stress-control gain K_p = ~1e-3 ˙γ
    PI-like update ˙γ←˙γ+(σ_fixed-σ_t)/σ_fixed * K_p with K_p≈10^{-3}˙γ; chosen to avoid unwanted transients.
  • bidisperse size ratio = 1 : 1.4
    1:1.4 radii and lengths to frustrate ordering; authors claim little rheology change vs mono/tri-disperse but still a setup choice.
axioms (8)
  • domain assumption Dense non-Brownian suspension stress is dominated by short-range contacts and lubrication; long-range many-body hydrodynamics may be omitted for tractability at high φ.
    Stated motivation versus Stokesian dynamics cost (§I); underpins entire force inventory.
  • domain assumption Single-body drag, lift, and Jeffery torque on spherocylinders equal those of prolate spheroids of the same aspect ratio.
    §II.B explicitly adopts Zhang/Brenner/Jeffery spheroid formulas in lieu of exact spherocylinder Stokes solutions.
  • domain assumption Pairwise lubrication is adequately represented by weighted Yamane SS plus Butler–Shaqfeh ES/EE normal squeeze terms only (no full resistance matrix).
    §II.C.3; tangential lubrication and higher-order couplings omitted.
  • domain assumption Contacts are Hookean with Coulomb friction and viscous damping set by restitution and contact time (Mahajan/Pournin-style DEM).
    §II.C.2, Eqs. 15–24.
  • domain assumption Particle inertia is negligible when ˙γ*=ρ_F ˙γ R^2/η<10^{-2}, so rate-independent viscous granular rheology (Boyer et al.) applies irrespective of shape.
    §I and §II.D.3; lift retained but argued negligible.
  • ad hoc to paper Fluid-phase shear stress contribution is η˙γ(1+5/2 φ+φ A^2/16) even though derived for dilute isotropic conditions.
    Eq. 28 (Kuhn); authors assert negligibility at high φ but still include it in η_s.
  • standard math Lees–Edwards boundaries plus linear imposed v_f produce simple shear representative of bulk rheology.
    §II.D.3; standard suspension DEM practice.
  • standard math Orientational order is captured by the largest eigenvalue of the nematic Q-tensor averaged over particles.
    §II.D.2, Eq. 29.
invented entities (2)
  • Weighted multi-mode spherocylinder lubrication blender (C_SS, C_ES, C_EE with manual correctors) no independent evidence
    purpose: Remove force discontinuities when pair geometry crosses shaft/end regimes while retaining Yamane/Butler–Shaqfeh asymptotes.
    Not a new physical force; an algorithmic closure. Independent evidence is only internal smoothness tests and a claim of weak rheology sensitivity—no external validation against resolved hydrodynamics.
  • Dynamic per-step timestep from min(t_c, t_lub,F, t_lub,T) no independent evidence
    purpose: Resolve stiff lubrication/contact events without a globally tiny fixed dt.
    Numerical device derived from approximate interaction timescales (Eqs. 19, 30–33); not a physical entity. Correctness judged by stability and velocity-profile diagnostics only.

pith-pipeline@v1.2.0-daily-grok45 · 26468 in / 4791 out tokens · 96840 ms · 2026-07-31T14:42:19.925444+00:00 · methodology

0 comments
read the original abstract

Dense suspensions of rod-shaped granular particles are widespread in nature and manufacturing, where their fluid mechanical properties are often paramount. We have developed a particle-based simulation that models such suspensions under simple shear flow, providing predictions of the viscosity and microstructure for a given solids volume fraction and particle aspect ratio. The model tracks the trajectories of spherocylindrical rods under the action of short-range frictional contact and hydrodynamic forces, inspired by similar tools that have generated new insight into suspensions of granular spheres. It incorporates new schemes for the computation of lubrication forces between spherocylinders and the dynamic determination of the timestep. For aspect ratios up to 20, the model predicts a viscosity spike at shear start-up, giving way to steady state viscosities that increase systematically with volume fraction and aspect ratio. Likewise, particle alignment increases with volume fraction up to an aspect-ratio-dependent critical point. Our model corroborates the limited experimental rheology data available for suspensions of granular rods, and offers a tool for fundamental exploration of the fluid mechanics, microstructure and rheology of this widespread material.

Figures

Figures reproduced from arXiv: 2607.28206 by Alex Dixon, Christopher Ness, Gavin Melaugh, John Hone.

Figure 1
Figure 1. Figure 1: FIG. 1. Particle-based simulation of suspensions of sheared granular spherocylindrical rods. Shown in (a) is a render of a typical [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Results of prerequisite calculations and diagnostic simulations used to determine parameter values. (a) Shown by [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Simulation snapshots of the initial configurations of granular rod suspensions prior to shearing. Shown are systems [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Simulation results showing the viscosity [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Dense suspensions of granular rods under shear, showing example snapshots of the steady state alignment obtained [PITH_FULL_IMAGE:figures/full_fig_p020_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Simulation results showing the alignment of granular rods in suspension under simple shear flow. Shown in (a) are [PITH_FULL_IMAGE:figures/full_fig_p021_6.png] view at source ↗

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

108 extracted references · 4 linked inside Pith

  1. [1]

    Spherocylinders have both a long and short dimension, meaning the minimum distance between their surfaces cannot be determined from the positions of their centres of mass

    Interaction Detection Before introducing the contact and lubrication forces, we first describe our algorithm for detecting whether neigh- bouring particles are close enough to be interacting. Spherocylinders have both a long and short dimension, meaning the minimum distance between their surfaces cannot be determined from the positions of their centres of...

  2. [2]

    Contact Forces Direct contacts between particles are modelled as Hookean with Coulombic friction, with normal and tangential components of the forces computed as follows, based on the work by Mahajanet al.[77]. The total contact force 9 acting on particleiis the sum of the normal ⃗FC,n,ij and tangential ⃗FC,t,ij components contributed by contacting neighb...

  3. [3]

    These forces originate in the viscosity of the fluid in the narrow gap between the particles, opposing their relative motion

    Lubrication forces Lubrication forces arise from short-range hydrodynamic interactions generated when two particles within a fluid move relative to each other in close proximity. These forces originate in the viscosity of the fluid in the narrow gap between the particles, opposing their relative motion. Lubrication forces become particularly important at ...

  4. [4]

    Stress Calculation The stress within the system is calculated as the sum of the stress contributions from pairwise particle interactions σP, and a fluid contributionσ F. The former is a 3×3 tensor constructed from the moment arm ⃗ℓand force ⃗Fvectors, and is defined as [81] σP =− 1 V   NCX ⃗ℓij ⊗ ⃗FC,ij + NLX ⃗ℓik ⊗ ⃗FL,ik,   (27) whereN C is the numb...

  5. [5]

    We quantify the degree of collective ordering of the particles by computing the largest eigenvalueSof the order tensorT, the 13 TABLE I

    Quantifying Alignment To understand the material response to flow, we characterise the microstructural arrangement of the particles borrowing insight from the Q-tensor commonly used to characterise alignment in nematic liquid crystals. We quantify the degree of collective ordering of the particles by computing the largest eigenvalueSof the order tensorT, ...

  6. [6]

    An example snapshot of the simulated system is given in Figure 1(a)

    Simulation Setup and Parameters Particles governed by the above set of forces and torques are placed with initially random positions and orien- tations in a cubic domain consisting of periodic boundaries in thexandydirections and Lees-Edwards boundary conditions [84] across thezaxis. An example snapshot of the simulated system is given in Figure 1(a). The...

    2003

  7. [7]

    The discrete time step dtused for this integration needs to be carefully chosen to resolve all interactions while not adding unnecessarily to the computation time

    Time step size The equations of motion given above are integrated in time to update the particle trajectories. The discrete time step dtused for this integration needs to be carefully chosen to resolve all interactions while not adding unnecessarily to the computation time. To achieve this we therefore need to approximate the characteristic timescale of e...

  8. [8]

    Throughout we prescribe a linear velocity profile by specifying a spatially invariant value of the gradient as ˙γ=∂v f,x /∂z

    System Control Modes As described in Section II B, the fluid velocity profile⃗ vf (z) determines the single-body hydrodynamic forces on each particle. Throughout we prescribe a linear velocity profile by specifying a spatially invariant value of the gradient as ˙γ=∂v f,x /∂z. The simulation can impose this shear rate under two modes. In the first, we simp...

    2022

  9. [9]

    J. E. Butler and B. Snook, Microstructural dynamics and rheology of suspensions of rigid fibers, Annual Review of Fluid Mechanics50, 299 (2018)

    2018

  10. [10]

    N. S. Deshpande and B. T. Crosby, Logjams are not jammed: measurements of log motions in Big Creek, Idaho, arXiv preprint arXiv:1911.01518 (2019)

    Pith/arXiv arXiv 1911

  11. [11]

    Cimarelli, A

    C. Cimarelli, A. Costa, S. Mueller, and H. M. Mader, Rheology of magmas with bimodal crystal size and shape distribu- tions: Insights from analog experiments, Geochemistry, Geophysics, Geosystems12(2011)

    2011

  12. [12]

    Moitra and H

    P. Moitra and H. Gonnermann, Effects of crystal shape-and size-modality on magma rheology, Geochemistry, Geophysics, Geosystems16, 1 (2015)

    2015

  13. [13]

    T. F. Choo, M. A. Mohd Salleh, K. Y. Kok, and K. A. Matori, A review on synthesis of mullite ceramics from industrial wastes, Recycling4, 39 (2019)

    2019

  14. [14]

    Zhang and B

    H. Zhang and B. W. Darvell, Synthesis and characterization of hydroxyapatite whiskers by hydrothermal homogeneous precipitation using acetamide, Acta Biomaterialia6, 3216 (2010)

    2010

  15. [15]

    K. Cao, G. Liu, H. Li, and Z. Huang, Mechanical properties and microstructure of calcium sulfate whisker-reinforced cement-based composites, Materials15, 947 (2022)

    2022

  16. [16]

    Derakhshandeh, R

    B. Derakhshandeh, R. J. Kerekes, S. G. Hatzikiriakos, and C. P. Bennington, Rheology of pulp fibre suspensions: A critical review, Chemical Engineering Science66, 3460 (2011)

    2011

  17. [17]

    Y. Jung, T. Plumb-Reyes, H.-Y. G. Lin, and L. Mahadevan, Entanglement transition in random rod packings, Proceedings of the National Academy of Sciences122, e2401868122 (2025)

    2025

  18. [18]

    G. G. Stokes, On the effect of the internal friction of fluids on the motion of pendulums, Transactions of the Cambridge Philosophical Society9, 8–106 (1851)

  19. [19]

    Einstein, Berichtigung zu meiner arbeit: Eine neue bestimmung der molek¨ uldimensionen, Annalen der Physik34, 591 (1911)

    A. Einstein, Berichtigung zu meiner arbeit: Eine neue bestimmung der molek¨ uldimensionen, Annalen der Physik34, 591 (1911)

    1911

  20. [20]

    Batchelor and J

    G. Batchelor and J. Green, The determination of the bulk stress in a suspension of spherical particles to order c2, Journal of Fluid Mechanics56, 401 (1972)

    1972

  21. [21]

    I. M. Krieger and T. J. Dougherty, A mechanism for non-newtonian flow in suspensions of rigid spheres, Transactions of the Society of Rheology3, 137 (1959)

    1959

  22. [22]

    Boyer, ´E

    F. Boyer, ´E. Guazzelli, and O. Pouliquen, Unifying suspension and granular rheology, Physical Review Letters107, 188301 (2011)

    2011

  23. [23]

    G. B. Jeffery, The motion of ellipsoidal particles immersed in a viscous fluid, Proceedings of the Royal Society of London 102, 161 (1922)

    1922

  24. [24]

    B¨ orzs¨ onyi, B

    T. B¨ orzs¨ onyi, B. Szab´ o, G. T¨ or¨ os, S. Wegner, J. T¨ or¨ ok, E. Somfai, T. Bien, and R. Stannarius, Orientational order and alignment of elongated particles induced by shear, Physical Review Letters108, 228302 (2012)

    2012

  25. [25]

    B¨ orzs¨ onyi, B

    T. B¨ orzs¨ onyi, B. Szab´ o, S. Wegner, K. Harth, J. T¨ or¨ ok, E. Somfai, T. Bien, and R. Stannarius, Shear-induced alignment and dynamics of elongated granular particles, Physical Review E86, 051304 (2012)

    2012

  26. [26]

    Franceschini, E

    A. Franceschini, E. Filippidi, E. Guazzelli, and D. J. Pine, Transverse alignment of fibers in a periodically sheared suspension: an absorbing phase transition with a slowly varying control parameter, Physical Review Letters107, 250603 (2011)

    2011

  27. [27]

    Snook, E

    B. Snook, E. Guazzelli, and J. E. Butler, Vorticity alignment of rigid fibers in an oscillatory shear flow: Role of confinement, Physics of Fluids24(2012). 24

    2012

  28. [28]

    Franceschini, E

    A. Franceschini, E. Filippidi, E. Guazzelli, and D. J. Pine, Dynamics of non-brownian fiber suspensions under periodic shear, Soft Matter10, 6722 (2014)

    2014

  29. [29]

    Ralambotiana, R

    T. Ralambotiana, R. Blanc, and M. Chaouche, Viscosity scaling in suspensions of non-brownian rodlike particles, Physics of Fluids9, 3588 (1997)

    1997

  30. [30]

    M. A. Bibb´ o,Rheology of semiconcentrated fiber suspensions, Ph.D. thesis, Massachusetts Institute of Technology (1987)

    1987

  31. [31]

    M. B. Mackaplow and E. S. Shaqfeh, A numerical study of the rheological properties of suspensions of rigid, non-brownian fibres, Journal of Fluid Mechanics329, 155 (1996)

    1996

  32. [32]

    Chaouche and D

    M. Chaouche and D. L. Koch, Rheology of non-brownian rigid fiber suspensions with adhesive contacts, Journal of Rheology45, 369 (2001)

    2001

  33. [33]

    Tapia, S

    F. Tapia, S. Shaikh, J. E. Butler, O. Pouliquen, and ´E. Guazzelli, Rheology of concentrated suspensions of non-colloidal rigid fibres, Journal of Fluid Mechanics827, R5 (2017)

    2017

  34. [34]

    Snook, L

    B. Snook, L. M. Davidson, J. E. Butler, O. Pouliquen, and E. Guazzelli, Normal stress differences in suspensions of rigid fibres, Journal of Fluid Mechanics758, 486 (2014)

    2014

  35. [35]

    S. R. Strand, S. Kim, and S. J. Karrila, Computation of rheological properties of suspensions of rigid rods: stress growth after inception of steady shear flow, Journal of Non-Newtonian Fluid Mechanics24, 311 (1987)

    1987

  36. [36]

    Qui˜ nones and P

    C. Qui˜ nones and P. D. Olmsted, A smoluchowski equation for a sheared suspension of frictionally interacting rods, arXiv preprint arXiv:2512.19149 (2025)

    Pith/arXiv arXiv 2025

  37. [37]

    Yamane, Y

    Y. Yamane, Y. Kaneda, and M. Doi, Numerical simulation of semi-dilute suspensions of rodlike particles in shear flow, Journal of Non-Newtonian Fluid Mechanics54, 405 (1994)

    1994

  38. [38]

    M. Khan, R. V. More, A. A. Banaei, L. Brandt, and A. M. Ardekani, Rheology of concentrated fiber suspensions with a load-dependent friction coefficient, Physical Review Fluids8, 044301 (2023)

    2023

  39. [39]

    J. F. Brady and G. Bossis, Stokesian dynamics, Annual Review of Fluid Mechanics20, 111 (1988)

    1988

  40. [40]

    T. N. Phung, J. F. Brady, and G. Bossis, Stokesian dynamics simulation of brownian suspensions, Journal of Fluid Mechanics313, 181 (1996)

    1996

  41. [41]

    D. R. Foss and J. F. Brady, Structure, diffusion and rheology of brownian suspensions by stokesian dynamics simulation, Journal of Fluid Mechanics407, 167 (2000)

    2000

  42. [42]

    Singh and P

    A. Singh and P. R. Nott, Normal stresses and microstructure in bounded sheared suspensions via stokesian dynamics simulations, Journal of Fluid Mechanics412, 279 (2000)

    2000

  43. [43]

    Kutteh, Methods for stokesian dynamics simulations of nonspherical particles and chains, Physical Review E69, 011406 (2004)

    R. Kutteh, Methods for stokesian dynamics simulations of nonspherical particles and chains, Physical Review E69, 011406 (2004)

    2004

  44. [44]

    Mari, Shear thickening of suspensions of dimeric particles, Journal of Rheology64, 239 (2020)

    R. Mari, Shear thickening of suspensions of dimeric particles, Journal of Rheology64, 239 (2020)

    2020

  45. [45]

    J. E. Butler and E. S. Shaqfeh, Dynamic simulations of the inhomogeneous sedimentation of rigid fibres, Journal of Fluid Mechanics468, 205 (2002)

    2002

  46. [46]

    X. Fan, N. Phan-Thien, and R. Zheng, A direct simulation of fibre suspensions, Journal of Non-Newtonian Fluid Mechanics 74, 113 (1998)

    1998

  47. [47]

    Durlofsky, J

    L. Durlofsky, J. F. Brady, and G. Bossis, Dynamic simulation of hydrodynamically interacting particles, Journal of Fluid Mechanics180, 21 (1987)

    1987

  48. [48]

    P. A. Cundall and O. D. Strack, A discrete numerical model for granular assemblies, Geotechnique29, 47 (1979)

    1979

  49. [49]

    Plimpton, Fast parallel algorithms for short-range molecular dynamics, Journal of Computational Physics117, 1 (1995)

    S. Plimpton, Fast parallel algorithms for short-range molecular dynamics, Journal of Computational Physics117, 1 (1995)

    1995

  50. [50]

    Yousefian and M

    Z. Yousefian and M. Trulsson, Orientational arrest in dense suspensions of elliptical particles under oscillatory shear flows, Europhysics Letters136, 36002 (2022)

    2022

  51. [51]

    Zhang and N

    Y. Zhang and N. Menon, Hindered stokesian settling of discs and rods, Physical Review Letters134, 168202 (2025). 25

    2025

  52. [52]

    Trulsson, Directional shear jamming of frictionless ellipses, Physical Review E104, 044614 (2021)

    M. Trulsson, Directional shear jamming of frictionless ellipses, Physical Review E104, 044614 (2021)

    2021

  53. [53]

    Bilotto, M

    J. Bilotto, M. Trulsson, and J.-F. Molinari, Shear flow of frictional spheroids: Comparison between elongated and flattened particles, Physical Review E112, 045432 (2025)

    2025

  54. [54]

    Ness, Simulating dense, rate-independent suspension rheology using lammps, Computational Particle Mechanics10, 2031 (2023)

    C. Ness, Simulating dense, rate-independent suspension rheology using lammps, Computational Particle Mechanics10, 2031 (2023)

    2031

  55. [55]

    Dong and M

    J. Dong and M. Trulsson, Unifying viscous and inertial regimes of discontinuous shear thickening suspensions, Journal of Rheology64, 255 (2020)

    2020

  56. [56]

    Ge and L

    Z. Ge and L. Brandt, Implementation note on a minimal hybrid lubrication/granular dynamics model for dense suspen- sions, arXiv preprint arXiv:2005.12755 (2020)

    Pith/arXiv arXiv 2005

  57. [57]

    Gallier, E

    S. Gallier, E. Lemaire, F. Peters, and L. Lobry, Rheology of sheared suspensions of rough frictional particles, Journal of Fluid Mechanics757, 514 (2014)

    2014

  58. [58]

    R. V. More and A. M. Ardekani, Unifying disparate rate-dependent rheological regimes in non-brownian suspensions, Physical Review E103, 062610 (2021)

    2021

  59. [59]

    Y. Guo, C. Wassgren, W. Ketterhagen, B. Hancock, B. James, and J. Curtis, A numerical study of granular shear flows of rod-like particles using the discrete element method, Journal of Fluid Mechanics713, 1 (2012)

    2012

  60. [60]

    A. Pol, R. Artoni, and P. Richard, Unified scaling law for wall friction in laterally confined flows of shape anisotropic particles, Physical Review Fluids8, 084302 (2023)

    2023

  61. [61]

    W. Nan, Y. Wang, Y. Liu, and H. Tang, Dem simulation of the packing of rodlike particles, Advanced Powder Technology 26, 527 (2015)

    2015

  62. [62]

    Y. Guo, C. Wassgren, B. Hancock, W. Ketterhagen, and J. Curtis, Computational study of granular shear flows of dry flexible fibres using the discrete element method, Journal of Fluid Mechanics775, 24 (2015)

    2015

  63. [63]

    Berzi, N

    D. Berzi, N. Thai-Quang, Y. Guo, and J. Curtis, Stresses and orientational order in shearing flows of granular liquid crystals, Physical Review E93, 040901 (2016)

    2016

  64. [64]

    Nath and C

    T. Nath and C. Heussinger, Rheology in dense assemblies of spherocylinders: frictional vs. frictionless, The European Physical Journal E42, 1 (2019)

    2019

  65. [65]

    Berzi, K

    D. Berzi, K. E. Buettner, and J. S. Curtis, Dense shearing flows of soft, frictional cylinders, Soft Matter18, 80 (2022)

    2022

  66. [66]

    Anzivino, C

    C. Anzivino, C. Ness, A. S. Moussa, and A. Zaccone, Shear flow of non-brownian rod-sphere mixtures near jamming, Physical Review E109, L042601 (2024)

    2024

  67. [67]

    Cheal and C

    O. Cheal and C. Ness, Rheology of dense granular suspensions under extensional flow, Journal of Rheology62, 501 (2018)

    2018

  68. [68]

    Seto and G

    R. Seto and G. G. Giusteri, Normal stress differences in dense suspensions, Journal of Fluid Mechanics857, 200 (2018)

    2018

  69. [69]

    Trulsson, B

    M. Trulsson, B. Andreotti, and P. Claudin, Transition from the viscous to inertial regime in dense suspensions, Physical Review Letters109, 118305 (2012)

    2012

  70. [70]

    Gan and A

    J. Gan and A. Yu, Dem simulation of the packing of cylindrical particles, Granular Matter22, 1 (2020)

    2020

  71. [71]

    Lin and S.-P

    A. Lin and S.-P. Han, On the distance between two ellipsoids, SIAM Journal on Optimization13, 298 (2002)

    2002

  72. [72]

    Ch´ eron, F

    V. Ch´ eron, F. Evrard, and B. van Wachem, Drag, lift and torque correlations for axi-symmetric rod-like non-spherical particles in locally linear shear flows, International Journal of Multiphase Flow171, 104692 (2024)

    2024

  73. [73]

    Feng and E

    Z. Feng and E. E. Michaelides, A general and accurate correlation for the drag on spherocylinders, International Journal of Multiphase Flow168, 104579 (2023)

    2023

  74. [74]

    Happel and H

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

    2012

  75. [75]

    Zhang, G

    H. Zhang, G. Ahmadi, F.-G. Fan, and J. B. McLaughlin, Ellipsoidal particles transport and deposition in turbulent channel flows, International Journal of Multiphase Flow27, 971 (2001)

    2001

  76. [76]

    Brenner, The stokes resistance of an arbitrary particle, Chemical Engineering Science18, 1 (1963)

    H. Brenner, The stokes resistance of an arbitrary particle, Chemical Engineering Science18, 1 (1963). 26

    1963

  77. [77]

    Fan and G

    F.-G. Fan and G. Ahmadi, A sublayer model for wall deposition of ellipsoidal particles in turbulent streams, Journal of Aerosol Science26, 813 (1995)

    1995

  78. [78]

    Harper and I.-D

    E. Harper and I.-D. Chang, Maximum dissipation resulting from lift in a slow viscous shear flow, Journal of Fluid Mechanics33, 209 (1968)

    1968

  79. [79]

    J. Wang, S. Yu, S. Luo, G. Xia, and L. Zong, Lift forces on axial symmetry particles rotating in a linear shear flow of a rarefied gas, Physics of Fluids30(2018)

    2018

  80. [80]

    Y. Cui, R. Decker, D. Janocha, C. K¨ orner, and B. Weigand, On constitutive models for the momentum transfer to particles in fluid-dominated two-phase flows, inAdvances in Mechanics of Materials and Structural Analysis: In Honor of Reinhold Kienzler(Springer International Publishing, Cham, 2018) pp. 1–25

    2018

Showing first 80 references.