Friction-mediated flow and jamming in a two-dimensional silo with two exit orifices
Pith reviewed 2026-05-25 02:08 UTC · model grok-4.3
The pith
Interparticle friction nonmonotonically tunes the frequency of unjamming events at a small silo orifice when a larger one flows steadily nearby.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The interparticle friction coefficient significantly influences the flow and jamming behavior of granular materials exiting through the orifice of a two-dimensional silo in the presence of another orifice located in its vicinity. The fluctuations emanating from a continuous flow through a larger orifice result in an intermittent flow through the smaller orifice consisting of sequential jamming and flowing events. The mean time duration of jammed and flow events, respectively, increase and decrease monotonically with increasing interparticle friction coefficient, while the frequency of unjamming instances shows a nonmonotonic behavior comprising an increase followed by a decrease. The overall
What carries the argument
Fluctuation-mediated unjamming at the smaller orifice, whose strength is set by interorifice distance and whose effectiveness is modulated by the interparticle friction coefficient.
If this is right
- Higher friction lengthens jammed intervals and shortens flowing intervals at the small orifice.
- Unjamming frequency reaches a maximum at intermediate friction before the system trends toward either permanent jamming or permanent flow.
- The location of the frequency maximum shifts with interorifice separation because that separation sets the amplitude of arriving fluctuations.
- Probability distributions of jammed times remain exponential and flowing times remain power-law tailed across different friction values and separations.
Where Pith is reading between the lines
- The nonmonotonic peak may mark the friction value at which the jamming threshold and the fluctuation amplitude are optimally matched.
- The same mechanism could produce analogous intermittent regimes in three-dimensional silos or when more than two orifices are present.
- Power-law tails in flowing times hint that the small orifice sits near a critical point whose location is shifted by friction.
Load-bearing premise
The strength of fluctuations that reach the smaller orifice is determined by the distance to the larger orifice and that distance systematically controls how unjamming frequency depends on friction.
What would settle it
Recording the unjamming frequency versus friction coefficient at fixed interorifice distance and observing a strictly monotonic curve instead of a rise followed by a fall.
Figures
read the original abstract
We show that the interparticle friction coefficient significantly influences the flow and jamming behavior of granular materials exiting through the orifice of a two-dimensional silo in the presence of another orifice located in its vicinity. The fluctuations emanating from a continuous flow through a larger orifice results in an intermittent flow through the smaller orifice consisting of sequential jamming and flowing events. The mean time duration of jammed and flow events, respectively, increase and decrease monotonically with increasing interparticle friction coefficient. The frequency of unjamming instances ($n_{u}$), however, shows a nonmonotonic behavior comprising an increase followed by a decrease with increasing friction coefficient. The decrease on either side of the maximum, then, represents a system moving progressively towards a permanently jammed or a permanently flowing state. The overall behavior shows a systematic dependence on the interorifice distance which determines the strength of the fluctuations reaching the smaller orifice leading to unjamming instances. The probability distributions of jamming and flowing times are nearly similar for different combinations of friction coefficients and interorifice distances studied and, respectively, exhibit exponential and power-law tails.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports discrete-element simulations of a 2D silo with two orifices of unequal size. It claims that interparticle friction μ controls jamming/flow at the smaller orifice via fluctuations from the larger orifice: mean jammed-event duration increases monotonically with μ while mean flowing-event duration decreases; unjamming frequency n_u is non-monotonic (rises then falls) with μ, the location and height of the peak depending on inter-orifice distance; jamming-time distributions are exponential and flowing-time distributions have power-law tails, both largely insensitive to the parameter combinations studied.
Significance. If the non-monotonic n_u(μ) trend and its distance dependence survive statistical scrutiny, the work supplies a concrete, fluctuation-mediated mechanism by which friction and geometry together tune the transition between permanently jammed, intermittent, and continuously flowing regimes in multi-orifice silos. The result is directly relevant to industrial hopper design and to the broader literature on force-chain fluctuations and unjamming.
major comments (2)
- [Results (unjamming-frequency plots)] Results section (plots of n_u versus μ for several inter-orifice distances): no error bars, no reported sample sizes or number of independent runs, and no statement of how post-hoc binning or run selection was performed. Without these, it is impossible to judge whether the reported non-monotonic peak is statistically robust or an artifact of finite sampling, directly undermining the central claim that n_u exhibits a maximum at intermediate friction.
- [Methods/Results (fluctuation-strength paragraph)] Methods and Results: the claim that inter-orifice distance “determines the strength of the fluctuations reaching the smaller orifice” is asserted but not quantified; no auxiliary measurement (e.g., force or velocity fluctuation spectra at the smaller orifice as a function of distance) is provided to support the causal link between distance and fluctuation amplitude.
minor comments (2)
- [Figure captions] Figure captions and axis labels should explicitly state the number of particles, the integration time step, and the precise definition (force or velocity threshold) used to classify a state as jammed versus flowing.
- [Abstract and distribution subsection] The abstract states that the probability distributions “are nearly similar for different combinations”; the text should report a quantitative measure (e.g., Kolmogorov–Smirnov distances) rather than a qualitative statement.
Simulated Author's Rebuttal
We thank the referee for the careful review and valuable suggestions. We address the major comments point-by-point below and will revise the manuscript accordingly to improve statistical reporting and quantification of fluctuations.
read point-by-point responses
-
Referee: Results section (plots of n_u versus μ for several inter-orifice distances): no error bars, no reported sample sizes or number of independent runs, and no statement of how post-hoc binning or run selection was performed. Without these, it is impossible to judge whether the reported non-monotonic peak is statistically robust or an artifact of finite sampling, directly undermining the central claim that n_u exhibits a maximum at intermediate friction.
Authors: We agree that error bars, sample sizes, and methodological details are necessary to establish robustness. In the revised manuscript we will add error bars (standard error of the mean) to all n_u(μ) curves, report that each data point is averaged over 20 independent runs of total duration 10^5 time units each, and explicitly state that n_u is computed as the total number of unjamming events divided by the aggregate flowing time across runs. The non-monotonic trend remains statistically significant under these controls. revision: yes
-
Referee: Methods and Results: the claim that inter-orifice distance “determines the strength of the fluctuations reaching the smaller orifice” is asserted but not quantified; no auxiliary measurement (e.g., force or velocity fluctuation spectra at the smaller orifice as a function of distance) is provided to support the causal link between distance and fluctuation amplitude.
Authors: We accept that a direct quantification strengthens the mechanistic interpretation. The revised manuscript will include a new paragraph and supplementary figure reporting the root-mean-square velocity fluctuations measured in a fixed region above the smaller orifice; these amplitudes decay monotonically with increasing inter-orifice separation, consistent with the observed distance dependence of the n_u peak. revision: yes
Circularity Check
No significant circularity
full rationale
The paper reports direct observations from discrete element method (DEM) simulations of granular flow through a two-orifice 2D silo. It defines jammed and flowing states from time series of particle velocities or discharge rates, then measures mean durations and unjamming frequency n_u as functions of friction coefficient μ and interorifice distance. These quantities are computed from the raw simulation trajectories without any intermediate fitting, rescaling, or analytical model whose parameters are tuned to the same data. Probability distributions of event durations are likewise extracted directly from the histograms of the same trajectories. No equations, uniqueness theorems, or self-citations are invoked to derive the reported trends; the central claims are therefore the simulation outputs themselves and do not reduce to their own inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Granular particles interact through contact forces that include friction
Reference graph
Works this paper leans on
-
[1]
K. To, P. Y. Lai, and H. K. Pak, Phys. Rev. Lett. 86, 71 (2001)
2001
-
[2]
Tang and R
J. Tang and R. P. Behringer, Chaos 21, 041107 (2011)
2011
-
[3]
Tewari, M
S. Tewari, M. Dichter, and B. Chakraborty, Soft Matter 9, 5016 (2013)
2013
-
[4]
Zuriguel, Pap
I. Zuriguel, Pap. Phys. 6, 060014 (2014)
2014
-
[5]
Zuriguel, L
I. Zuriguel, L. A. Pugnaloni, A. Garcimartı’n, and D. Maza, Phys. Rev. E 68, 030301(R) (2003)
2003
-
[6]
Janda, I
A. Janda, I. Zuriguel, A. Garcimart´ ın, A. Pugnaloni, and D. Maza, Eur. Phys. Lett. 84, 44002 (2008)
2008
-
[7]
Janda, D
A. Janda, D. Maza, A. Garcimart´ ın, E. Kolb, J. Lanuze, and E. Cl´ ement, Eur. Phys. Lett.87, 24002 (2009)
2009
-
[8]
Kunte, P
A. Kunte, P. Doshi, and A. V. Orpe, Phys. Rev. E (Rapid Comm.) 90, 020201(R) (2014)
2014
-
[9]
Mondal and M
S. Mondal and M. M. Sharma, Granular Matter 16, 125 (2014)
2014
-
[10]
Kamath, A
S. Kamath, A. Kunte, P. Doshi, and A. V. Orpe, Phys. Rev. E 90, 062206 (2014)
2014
-
[11]
Zuriguel, A
I. Zuriguel, A. Janda, A. Garcimart´ ın, C. Lozano, R. Ar´ evalo, and D. Maza, Phys. Rev. Lett.107, 278001 (2011)
2011
-
[12]
Zuriguel, D
I. Zuriguel, D. R. Parisi, R. C. Hidalgo, C. Lozano, A. Janda, P. A. Gago, J. P. Peralta, L. M. Ferrer, L. A. Pugnaloni, E. Cl´ ement, D. Maza, I. Pagonabarraga, and A. Garcimart´ ın, Sci. Rep.4, 7324 (2014)
2014
-
[13]
Zuriguel, A
I. Zuriguel, A. Janda, R. Ar´ evalo, D. Maza, and A. Garcimart´ ın, EPJ Web. Conf.140, 01002 (2017)
2017
-
[14]
R. C. Hidalgo, A. G. ni Arana, A. Hern´ andez-Puerta, and I. Pagonabarraga, Phys. Rev. E 97, 012611 (2018)
2018
-
[15]
Nicolas, A
A. Nicolas, A. Garcimart´ ın, and I. Zuriguel, Phys. Rev. Lett. 120, 198002 (2018)
2018
-
[16]
L. A. Fullard, E. C. P. Breard, C. E. Davies, A. J. R. Godfrey, M. Fukuoka, A. Wade, J. Dufek, and G. Lube, Proc. R. Soc. A 475, 20180462 (2019)
2019
-
[17]
Maiti, G
R. Maiti, G. Das, and P. K. Das, Phys. Fluids 29, 103303 (2017)
2017
-
[18]
C. H. Rycroft, A. V. Orpe, and A. Kudrolli, Phys. Rev. E 80, 031305 (2009)
2009
-
[19]
http://lammps.sandia.gov/
-
[20]
J. W. Landry, G. S. Grest, L. E. Silbert, and S. J. Plimp- ton, Phys. Rev. E 67, 041303 (2003)
2003
-
[21]
K. A. Reddy, Y. Forterre, and O. Pouliquen, Phys. Rev. Lett. 106, 108301 (2011)
2011
-
[22]
Nichol and M
K. Nichol and M. van Hecke, Phys. Rev. E 85, 061309 (2012)
2012
-
[23]
Kamrin and G
K. Kamrin and G. Koval, Phys. Rev. Lett. 108, 178301 (2012)
2012
-
[24]
G. A. Patterson, P. I. Fierens, F. S. Jimka, P. G. K ˙.onig, A. Garcimart´ ın, I. Zuriguel, L. A. Pugnaloni, and D. R. Parisi, Phys. Rev. Lett. 119, 248301 (2017)
2017
-
[25]
C. C. Thomas and D. J. Durian, Phys. Rev. E 87, 052201 (2013)
2013
-
[26]
C. C. Thomas and D. J. Durian, Phys. Rev. Lett. 114, 178001 (2015)
2015
-
[27]
C. C. Thomas and D. J. Durian, Phys. Rev. E 94, 022901 (2016)
2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.