pith. sign in

arxiv: 2507.00684 · v1 · submitted 2025-07-01 · ❄️ cond-mat.soft · cond-mat.dis-nn· cond-mat.stat-mech

Structural Order Drives Diffusion in a Granular Packing

David Luce , Adrien Gans , S\'ebastien Kiesgen de Richter , Nicolas Vandewalle This is my paper

Pith reviewed 2026-05-19 06:50 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.dis-nncond-mat.stat-mech
keywords granular flowsilo dischargecrystallizationhexatic orderdiffusion lengthstructural ordervelocity profilesbidisperse particles
0
0 comments X

The pith

Crystallization enhances the diffusion length that governs velocity profiles in granular silo flows.

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

The paper investigates flows of bidisperse granular particles in a quasi-two-dimensional silo and shows that greater local crystalline order produces a larger diffusion length b. This length b sets the spatial scale over which velocities vary across the flowing region. Adjusting the mass fraction of the two particle sizes tunes the amount of ordering and produces a clear correlation between b and the hexatic order parameter. Pressure gradients further stabilize orientational order even without full crystallization, causing b to grow with height. The work therefore connects particle-level arrangement and pressure to the overall transport behavior.

Core claim

Crystallization significantly enhances the diffusion length b, a key parameter controlling the velocity profiles within the flowing medium. A strong correlation exists between b and the hexatic order parameter <|ψ6|>t. Pressure gradients within the silo promote the stabilization of orientational order even in the absence of crystallization, intrinsically increasing b with height. These findings establish a direct link between microstructural order, pressure, and transport properties in granular silo flows.

What carries the argument

The diffusion length b extracted via kinematic modeling of high-speed images, which scales the velocity profiles and correlates with the time-averaged hexatic order parameter.

If this is right

  • Velocity profiles broaden as the degree of crystalline order rises.
  • Local structural organization directly controls macroscopic flow behavior.
  • b increases with height because pressure stabilizes orientational order without requiring full crystallization.
  • Particle size distribution can tune flow properties through its effect on ordering.

Where Pith is reading between the lines

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

  • Adjusting grain size ratios could provide a practical way to modify discharge rates in industrial silos.
  • The same structure-transport link may appear in other confined granular systems such as hoppers or shear bands.
  • Experiments that hold order fixed while varying pressure independently could separate the two contributions to b.

Load-bearing premise

Kinematic modeling from high-speed imaging extracts the diffusion length b without significant bias from resolution limits or particle tracking errors.

What would settle it

Direct measurement of velocity fields in packings with controlled order but constant b, or inconsistent b values when the same flows are reanalyzed with different tracking algorithms, would break the reported correlation.

Figures

Figures reproduced from arXiv: 2507.00684 by Adrien Gans, David Luce, Nicolas Vandewalle, S\'ebastien Kiesgen de Richter.

Figure 1
Figure 1. Figure 1: FIG. 1. Four snapshots of the bidisperse granular media in [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Probability Distribution Function (PDF) of Voronoi [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Normalized vertical velocity [PITH_FULL_IMAGE:figures/full_fig_p002_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Time-average field of the normalized vertical velocity field [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Evolution of the normalized parameter [PITH_FULL_IMAGE:figures/full_fig_p003_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Correlation between the local modulus of the hexatic [PITH_FULL_IMAGE:figures/full_fig_p004_6.png] view at source ↗
read the original abstract

We investigate how structural ordering, i.e. crystallization, affects the flow of bidisperse granular materials in a quasi-two-dimensional silo. By systematically varying the mass fraction of two particle sizes, we finely tune the degree of local order. Using high-speed imaging and kinematic modeling, we show that crystallization significantly enhances the diffusion length $b$, a key parameter controlling the velocity profiles within the flowing medium. We reveal a strong correlation between $b$ and the hexatic order parameter $\left<|\psi_6|\right>_t$, highlighting the role of local structural organization in governing macroscopic flow behavior. Furthermore, we demonstrate that pressure gradients within the silo promote the stabilization of orientational order even in the absence of crystallization, thus intrinsically increasing $b$ with height. These findings establish a direct link between microstructural order, pressure, and transport properties in granular silo flows.

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

1 major / 2 minor

Summary. The paper investigates the effect of structural ordering (crystallization) on the flow behavior of bidisperse granular materials in a quasi-two-dimensional silo. By systematically varying the mass fraction of two particle sizes to tune the degree of local order, and employing high-speed imaging combined with kinematic modeling, the authors report that crystallization enhances the diffusion length b that controls velocity profiles. They demonstrate a strong correlation between b and the time-averaged hexatic order parameter <|ψ6|>t, and further show that pressure gradients can stabilize orientational order even without full crystallization, causing b to increase with height.

Significance. If the central correlation holds after addressing potential methodological biases, the work would establish a direct connection between local microstructural order and a key macroscopic transport parameter in granular silo flows. This has potential implications for granular rheology and silo discharge modeling. The systematic composition variation and use of an independently measured order parameter (rather than a derived one) are positive features supporting the claim.

major comments (1)
  1. [Kinematic modeling and results on b extraction] Kinematic modeling section: The diffusion length b is extracted by fitting a standard exponential form (v(x) ~ exp(-x/b)) to particle velocity data from high-speed imaging. The manuscript does not report stratified goodness-of-fit statistics (e.g., R² or residual distributions) separated by regions of high vs. low <|ψ6|>t, nor any cross-validation against model-free particle image velocimetry. Because coherent motions in crystalline regions can stabilize fits while noisier trajectories in disordered regions can bias b downward, this differential sensitivity is load-bearing for the reported correlation strength and must be quantified to rule out artifact.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'significantly enhances' is used without accompanying quantitative measure (e.g., factor by which b increases or correlation coefficient value); adding the actual correlation strength or effect size would improve clarity.
  2. [Results on pressure and order] The description of pressure-gradient effects on order stabilization is intriguing but would benefit from an explicit statement of how height-dependent pressure is measured or estimated in the silo geometry.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive comments and positive assessment of the significance of our work. We address the major comment below and will incorporate the requested analysis into the revised manuscript.

read point-by-point responses

  1. Referee: [Kinematic modeling and results on b extraction] Kinematic modeling section: The diffusion length b is extracted by fitting a standard exponential form (v(x) ~ exp(-x/b)) to particle velocity data from high-speed imaging. The manuscript does not report stratified goodness-of-fit statistics (e.g., R² or residual distributions) separated by regions of high vs. low <|ψ6|>t, nor any cross-validation against model-free particle image velocimetry. Because coherent motions in crystalline regions can stabilize fits while noisier trajectories in disordered regions can bias b downward, this differential sensitivity is load-bearing for the reported correlation strength and must be quantified to rule out artifact.

    Authors: We agree that providing stratified goodness-of-fit statistics and cross-validation is important to confirm that the extracted diffusion length b is not biased by differences in trajectory coherence between crystalline and disordered regions. In the revised manuscript we will report R² values together with residual distributions for velocity-profile fits, stratified explicitly by regions of high versus low time-averaged hexatic order parameter. We will also add a direct comparison of the kinematic-model results against model-free particle-image-velocimetry fields to demonstrate that the reported correlation between b and <|ψ6|>t remains robust under these checks. revision: yes

Circularity Check

0 steps flagged

No circularity: empirical correlation between independently extracted b and hexatic order

full rationale

The paper extracts the diffusion length b via kinematic fits to high-speed imaging velocity data and computes the hexatic order parameter <|ψ6|>t directly from particle position configurations. The reported correlation is between these two separately measured quantities. No equation or step defines b in terms of the order parameter (or vice versa), no fitted parameter is relabeled as a prediction, and no self-citation chain or uniqueness theorem is invoked to force the central claim. The derivation chain remains self-contained against external data benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard experimental assumptions in granular physics plus one key modeling choice; no new entities are postulated.

free parameters (1)
  • diffusion length b
    Fitted parameter in the kinematic model that controls the shape of velocity profiles; its value is extracted from data rather than predicted a priori.
axioms (1)
  • domain assumption The hexatic order parameter accurately quantifies local structural ordering in bidisperse granular packings.
    Invoked when correlating <|ψ6|>t with b to claim that order governs diffusion.

pith-pipeline@v0.9.0 · 5691 in / 1276 out tokens · 33489 ms · 2026-05-19T06:50:42.220477+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

22 extracted references · 22 canonical work pages

  1. [1]

    Torquato and F

    S. Torquato and F. H. Stillinger, Jammed hard-particle packings: From kepler to bernal and beyond, Reviews of modern physics 82, 2633 (2010)

    work page 2010

  2. [2]

    Watanabe and H

    K. Watanabe and H. Tanaka, Direct observation of medium-range crystalline order in granular liquids near the glass transition, Physical review letters 100, 158002 (2008)

    work page 2008

  3. [3]

    Kondic, A

    L. Kondic, A. Goullet, C. O’Hern, M. Kramar, K. Mis- chaikow, and R. Behringer, Topology of force networks in compressed granular media, Europhysics Letters 97, 54001 (2012)

    work page 2012

  4. [4]

    P. M. Reis, R. A. Ingale, and M. D. Shattuck, Crystalliza- tion of a quasi-two-dimensional granular fluid, Physical review letters 96, 258001 (2006)

    work page 2006

  5. [5]

    Benyamine, M

    M. Benyamine, M. Djermane, B. Dalloz-Dubrujeaud, and P. Aussillous, Discharge flow of a bidisperse granular me- dia from a silo, Physical Review E 90, 032201 (2014). 5

    work page 2014

  6. [6]

    Downs, N

    J. Downs, N. Smith, K. Mandadapu, J. Garrahan, and M. Smith, Topographic control of order in quasi-2d gran- ular phase transitions, Physical Review Letters 127, 268002 (2021)

    work page 2021

  7. [7]

    J. Bai, J. Li, G. Hong, J. Pan, and H. Fei, Mesoscopic evolution and kinetic properties of dense granular flow crystallization under continuous shear induction, Powder Technology 426, 118615 (2023)

    work page 2023

  8. [8]

    C. M. Carlevaro, R. Kozlowski, and L. A. Pugnaloni, Flow rate in 2d silo discharge of binary granular mixtures: the role of ordering in monosized systems, Frontiers in Soft Matter 4, 1340744 (2024)

    work page 2024

  9. [9]

    G. Hagen, ¨Uber den Druck und die Bewegung des trocknen Sandes, Bericht ¨ uber die zur Bekanntmachung geeigneten Verhandlungen der K¨ oniglich Preussischen Akademie der Wissenschaften zu Berlin , 35 (1852)

    work page

  10. [10]

    W. A. Beverloo, H. A. Leniger, and J. Van de Velde, The flow of granular solids through orifices, Chemical engineering science 15, 260 (1961)

    work page 1961

  11. [11]

    Janda, I

    A. Janda, I. Zuriguel, and D. Maza, Flow rate of parti- cles through apertures obtained from self-similar density and velocity profiles, Physical review letters 108, 248001 (2012)

    work page 2012

  12. [12]

    Litwiniszyn, An application of the random walk argu- ment to the mechanics of granular media, Rheology and Soil Mechanics: Symposium Grenoble, April 1–8, 1964 , 82 (1964)

    J. Litwiniszyn, An application of the random walk argu- ment to the mechanics of granular media, Rheology and Soil Mechanics: Symposium Grenoble, April 1–8, 1964 , 82 (1964)

    work page 1964

  13. [13]

    Mullins, Experimental evidence for the stochastic the- ory of particle flow under gravity, Powder Technology 9, 29 (1974)

    W. Mullins, Experimental evidence for the stochastic the- ory of particle flow under gravity, Powder Technology 9, 29 (1974)

    work page 1974

  14. [14]

    Nedderman and U

    R. Nedderman and U. T¨ uz¨ un, A kinematic model for the flow of granular materials, Powder Technology 22, 243 (1979)

    work page 1979

  15. [15]

    Medina, J

    A. Medina, J. Cordova, E. Luna, and C. Trevino, Velocity field measurements in granular gravity flow in a near 2d silo, Physics Letters A 250, 111 (1998)

    work page 1998

  16. [16]

    J. Choi, A. Kudrolli, and M. Z. Bazant, Velocity pro- file of granular flows inside silos and hoppers, Journal of Physics: Condensed Matter 17, S2533 (2005)

    work page 2005

  17. [17]

    Benyamine, P

    M. Benyamine, P. Aussillous, and B. Dalloz-Dubrujeaud, Discharge flow of a granular media from a silo: effect of the packing fraction and of the hopper angle, EPJ Web of Conferences 140, 03043 (2017)

    work page 2017

  18. [18]

    Perge, M

    C. Perge, M. A. Aguirre, P. A. Gago, L. A. Pugnaloni, D. Le Tourneau, and J.-C. G´ eminard, Evolution of pres- sure profiles during the discharge of a silo, Physical Re- view E—Statistical, Nonlinear, and Soft Matter Physics 85, 021303 (2012)

    work page 2012

  19. [19]

    Samadani, A

    A. Samadani, A. Pradhan, and A. Kudrolli, Size segre- gation of granular matter in silo discharges, Physical Re- view E 60, 7203 (1999)

    work page 1999

  20. [20]

    Q. Chen, R. Li, W. Xiu, V. Zivkovic, and H. Yang, Mea- surement of granular temperature and velocity profile of granular flow in silos, Powder Technology 392, 123 (2021)

    work page 2021

  21. [21]

    Zuriguel, D

    I. Zuriguel, D. Maza, A. Janda, R. C. Hidalgo, and A. Garcimart´ ın, Velocity fluctuations inside two and three dimensional silos, Granular Matter 21, 1 (2019)

    work page 2019

  22. [22]

    Walker, An approximate theory for pressures and arching in hoppers, Chemical Engineering Science 21, 975 (1966)

    D. Walker, An approximate theory for pressures and arching in hoppers, Chemical Engineering Science 21, 975 (1966)

    work page 1966