No Arabic abstract
The effects of movement of the side walls of a confined granular packing are studied by discrete element, molecular dynamics simulations. The dynamical evolution of the stress is studied as a function of wall movement both in the direction of gravity as well as opposite to it. For all wall velocities explored, the stress in the final state of the system after wall movement is fundamentally different from the original state obtained by pouring particles into the container and letting them settle under the influence of gravity. The original packing possesses a hydrostatic-like region at the top of the container which crosses over to a depth-independent stress. As the walls are moved in the direction opposite to gravity, the saturation stress first reaches a minimum value independent of the wall velocity, then increases to a steady-state value dependent on the wall-velocity. After wall movement ceases and the packing reaches equilibrium, the stress profile fits the classic Janssen form for high wall velocities, while it has some deviations for low wall velocities. The wall movement greatly increases the number of particle-wall and particle-particle forces at the Coulomb criterion. Varying the wall velocity has only small effects on the particle structure of the final packing so long as the walls travel a similar distance.
We report numerical results of effective attractive forces on the packing properties of two-dimensional elongated grains. In deposits of non-cohesive rods in 2D, the topology of the packing is mainly dominated by the formation of ordered structures of aligned rods. Elongated particles tend to align horizontally and the stress is mainly transmitted from top to bottom, revealing an asymmetric distribution of local stress. However, for deposits of cohesive particles, the preferred horizontal orientation disappears. Very elongated particles with strong attractive forces form extremely loose structures, characterized by an orientation distribution, which tends to a uniform behavior when increasing the Bond number. As a result of these changes, the pressure distribution in the deposits changes qualitatively. The isotropic part of the local stress is notably enhanced with respect to the deviatoric part, which is related to the gravity direction. Consequently, the lateral stress transmission is dominated by the enhanced disorder and leads to a faster pressure saturation with depth.
For packings of hard but not perfectly rigid particles, the length scales that govern the packing geometry and the contact forces are well separated. This separation of length scales is explored in the force network ensemble, where one studies the space of allowed force configurations for a given, frozen contact geometry. Here we review results of this approach, which yields nontrivial predictions for the effect of packing dimension and anisotropy on the contact force distribution $P(f)$, the response to overall shear and point forcing, all of which can be studied in great numerical detail. Moreover, there are emerging analytical approaches that very effectively capture, for example, the form of force distributions.
The structure and stresses of static granular packs in cylindrical containers are studied using large-scale discrete element molecular dynamics simulations in three dimensions. We generate packings by both pouring and sedimentation and examine how the final state depends on the method of construction. The vertical stress becomes depth-independent for deep piles and we compare these stress depth-profiles to the classical Janssen theory. The majority of the tangential forces for particle-wall contacts are found to be close to the Coulomb failure criterion, in agreement with the theory of Janssen, while particle-particle contacts in the bulk are far from the Coulomb criterion. In addition, we show that a linear hydrostatic-like region at the top of the packings unexplained by the Janssen theory arises because most of the particle-wall tangential forces in this region are far from the Coulomb yield criterion. The distributions of particle-particle and particle-wall contact forces $P(f)$ exhibit exponential-like decay at large forces in agreement with previous studies.
We present a multiscale simulation algorithm for amorphous materials, which we illustrate and validate in a canonical case of dense granular flow. Our algorithm is based on the recently proposed Spot Model, where particles in a dense random packing undergo chain-like collective displacements in response to diffusing spots of influence, carrying a slight excess of interstitial free volume. We reconstruct the microscopic dynamics of particles from the coarse grained dynamics of spots by introducing a localized particle relaxation step after each spot-induced block displacement, simply to enforce packing constraints with a (fairly arbitrary) soft-core repulsion. To test the model, we study to what extent it can describe the dynamics of up to 135,000 frictional, viscoelastic spheres in granular drainage simulated by the discrete-element method (DEM). With only five fitting parameters (the radius, volume, diffusivity, drift velocity, and injection rate of spots), we find that the spot simulations are able to largely reproduce not only the mean flow and diffusion, but also some subtle statistics of the flowing packings, such as spatial velocity correlations and many-body structural correlations. The spot simulations run over 100 times faster than DEM and demonstrate the possibility of multiscale modeling for amorphous materials, whenever a suitable model can be devised for the coarse-grained spot dynamics.
We consider a granular gas under the action of gravity, fluidized by a vibrating base. We show that a horizontal temperature gradient, here induced by limiting dissipative lateral walls (DLW), leads always to a granular thermal convection (DLW-TC) that is essentially different from ordinary buoyancy-driven convection (BD-TC). In an experiment where BD-TC is inhibited, by reducing gravity with an inclined plane, we always observe a DLW-TC cell next to each lateral wall. Such a cell squeezes towards the nearest wall as the gravity and/or the number of grains increase. Molecular dynamics simulations reproduce the experimental results and indicate that at large gravity or number of grains the DLW-TC is barely detectable.