No Arabic abstract
We carry out a direct comparison of experimental and numerical realizations of the exact same granular system as it undergoes shear jamming. We adjust the numerical methods used to optimally represent the experimental settings and outcomes up to microscopic contact force dynamics. Measures presented here range form microscopic, through mesoscopic to system-wide characteristics of the system. Topological properties of the mesoscopic force networks provide a key link between micro and macro scales. We report two main findings: the number of particles in the packing that have at least two contacts is a good predictor for the mechanical state of the system, regardless of strain history and packing density. All measures explored in both experiments and numerics, including stress tensor derived measures and contact numbers depend in a universal manner on the fraction of non-rattler particles, $f_{NR}$. The force network topology also tends to show this universality, yet the shape of the master curve depends much more on the details of the numerical simulations. In particular we show that adding force noise to the numerical data set can significantly alter the topological features in the data. We conclude that both $f_{NR}$ and topological metrics are useful measures to consider when quantifying the state of a granular system.
A remarkable feature of static granular matter is the distribution of force along intricate networks. Even regular inter-particle contact networks produce wildly inhomogeneous force networks where certain chains of particles carry forces far larger than the mean. In this paper, we briefly review past theoretical approaches to understanding the geometry of force networks. We then investigate the structure of experimentally-obtained granular force networks using a simple algorithm to obtain corresponding graphs. We compare our observations with the results of geometric models, including random bond percolation, which show similar spatial distributions without enforcing vector force balance. Our findings suggest that some aspects of the mean geometry of granular force networks may be captured by these simple descriptions.
We have made experimental observations of the force networks within a two-dimensional granular silo similar to the classical system of Janssen. Models like that of Janssen predict that pressure within a silo saturates with depth as the result of vertical forces being redirected to the walls of the silo where they can then be carried by friction. By averaging ensembles of experimentally-obtained force networks in different ways, we compare the observed behavior with various predictions for granular silos. We identify several differences between the mean behavior in our system and that predicted by Janssen-like models: We find that the redirection parameter describing how the force network transfers vertical forces to the walls varies with depth. We find that changes in the preparation of the material can cause the pressure within the silo to either saturate or to continue building with depth. Most strikingly, we observe a non-linear response to overloads applied to the top of the material in the silo. For larger overloads we observe the previously reported giant overshoot effect where overload pressure decays only after an initial increase [G. Ovarlez et al., Phys. Rev. E 67, 060302(R) (2003)]. For smaller overloads we find that additional pressure propagates to great depth. This effect depends on the particle stiffness, as given for instance by the Youngs modulus, E, of the material from which the particles are made. Important measures include E, the unscreened hydrostatic pressure, and the applied load. These experiments suggest that when the load and the particle weight are comparable, particle elasticity acts to stabilize the force network, allowing non-linear network effects to be seen in the mean behavior.
Large-scale three dimensional molecular dynamics simulations of hopper flow are presented. The flow rate of the system is controlled by the width of the aperture at the bottom. As the steady-state flow rate is reduced, the force distribution $P(f)$ changes only slightly, while there is a large change in the impulse distribution $P(i)$. In both cases, the distributions show an increase in small forces or impulses as the systems approach jamming, the opposite of that seen in previous Lennard-Jones simulations. This occurs dynamically as well for a hopper that transitions from a flowing to a jammed state over time. The final jammed $P(f)$ is quite distinct from a poured packing $P(f)$ in the same geometry. The change in $P(i)$ is a much stronger indicator of the approach to jamming. The formation of a peak or plateau in $P(f)$ at the average force is not a general feature of the approach to jamming.
We experimentally investigate the fluidization of a granular material subject to mechanical vibrations by monitoring the angular velocity of a vane suspended in the medium and driven by an external motor. On increasing the frequency we observe a re-entrant transition, as a jammed system first enters a fluidized state, where the vane rotates with high constant velocity, and then returns to a frictional state, where the vane velocity is much lower. While the fluidization frequency is material independent, the viscosity recovery frequency shows a clear dependence on the material, that we rationalize by relating this frequency to the balance between dissipative and inertial forces in the system. Molecular dynamics simulations well reproduce the experimental data, confirming the suggested theoretical picture.
Combining X-ray tomography with simultaneous shear force measurement, we investigate shear-induced granular avalanches using spherical particles with different surface roughness. We find that systems consisting of particles with large surface roughness display quasi-periodic avalanches interrupted by crackling-like small ones. In contrast, systems consisting of particles with small roughness display no detectable avalanches. The stress drop of quasi-periodic avalanche shows a linear relation with the correlation length of particle non-affine displacement, suggesting that roughness enhances inter-particle locking and hence particle-level dynamic correlation length. However, the nonaffine displacement is two orders of magnitude smaller than particle size, indicating that stress is mainly released on the length scale of roughness. The correlation length of non-affine displacements abruptly increases when a quasi-periodic avalanche occurs, suggesting that quasi-periodic avalanches can be interpreted as a spinodal nucleation event in a first-order phase transition.