Data for Hong et al. (2017) and Tao et al. (2021)

• "Clogging of soft particles in 2D hoppers"
  X Hong, M Kohne, M Morrell, H Wang, & ER Weeks, Phys. Rev. E 96, 062605 (2017).
  • Abstract / PDF / journal / arXiv:1512.02500
  • Publicity: Emory news story
• "Soft particle clogging in two-dimensional hoppers"
  R Tao, M Wilson, & ER Weeks, Phys. Rev. E 104, 044909 (2021)
  • Abstract / PDF / journal / arXiv:2104.01698

Data for figures

Data from Hong 2017 Fig. 5: data-fig5.txt. Format of data: text file, columns = (panel of figure, gravity in g, w/d, Pclog, symbol). Symbol explains what symbol of the figure the data corresponds to, which should be consistent with the value of gravity.

Data from Hong 2017 Fig. 6: data-fig6.txt. Format of data: text file, columns = (delta/d, w/d(Pclog=1/2), symbol).

Data from Tao 2021 Fig. 4: sigmoid_dat.txt. Format of data: text file, columns = (delta/d, w/d(Pclog=1/2), width s of sigmoid fit, particle type).

Representative simulation data

The data are in the format [x,y,radius,time,ID]. (Each particle has a unique ID number.) The simulation code is available on GitHub (written in IDL).

filename	w/d	gravity	final state
w20g01h.0000.txt (in a zip file)	1.0	1e-2	clogged
w36g01h.0000.txt (in a zip file)	1.8	1e-2	clogged
w36g01h.0001.txt (in a zip file)	1.8	1e-2	did not clog
w60g01h.0000.txt (in a zip file)	3.0	1e-2	did not clog
w50g001h.0000.txt (in a zip file)	2.5	1e-3	did not clog
w50g001h.0001.txt (in a zip file)	2.5	1e-3	clogged

Arch data

These data are specifically from Tao et al, PRE 2021.

Starting with 800 soft particles, let them flow out of a hopper, see if they clog. If they clog, save the arch information and final state information:

• archall-new.txt -- all arches except for three big runs
• arches-w125g01.txt -- big run with w/d = 1.25, g = 0.01, 1579 simulations. Clogged every time (P_clog = 1).
• final-w125g01.txt -- final state of the 1579 simulations
• arches-w16g01.txt -- big run with w/d = 1.6, g = 0.01, 1697 simulations. Clogged 1456 times, drained out 241 times (P_clog = 0.858).
• final-w16g01.txt -- final state of the 1456 simulations that clogged
• arches-w18g01.txt -- big run with w/d = 1.8, g = 0.01, 1454 simulations. Clogged 778 times, drained out 676 times (P_clog = 0.535).
• final-w18g01.txt -- final state of the 778 simulations that clogged

Format of 'arch' data files:

• 1st column = number of particles in arch N
• columns 2 through 3N+1: data (x,y,radius) for each particle in arch. Particles should be listed in order of increasing x.
• columns 3N+2, 3N+3: (x,y) of closest point on left hopper wall to the arch particle touching that wall. That is, from this (x,y) position to the (x,y) of the arch particle should be perpendicular to the wall.
• columns 3N+4, 3N+5: (x,y) of closest point on right hopper wall. As a side effect, for a given value of w/d, plotting these columns and the previous two give you various points on two lines, which indicates where the hopper wall is.
• column 47: how many particles left in hopper when clogged (includes arch particles)
• column 49: width in terms of diameter, that is, w/d
• column 50: gravity value

Format of 'final' data files: [x, y, radius, final time, simulation number]. For example, the file final-w18g01.txt starts with the entry [0,0,0,0,0] indicating that simulation #0 did not clog. Then rows 2 through 130 have the form [x,y,r,2442,1] indicating that simulation #1 clogged at time 2442, with final particles in the positions given. Row 131 has the entry [0,0,0,0,2] indicating that simulation #2 did not clog.

If you take a histogram of column 47 of the three arches-w* files, that should reproduce Fig. 5(b) of Tao PRE 2021. Likewise Fig. 5(a) is derivable from that data although less directly.