No Arabic abstract
Cohesive granular media flowing down an inclined plane are studied by discrete element simulations. Previous work on cohesionless granular media demonstrated that within the steady flow regime where gravitational energy is balanced by dissipation arising from intergrain forces, the velocity profile in the flow direction scales with depth in a manner consistent with the predictions of Bagnold. Here we demonstrate that this Bagnold scaling does not hold for the analogous steady-flows in cohesive granular media. We develop a generalization of the Bagnold constitutive relation to account for our observation and speculate as to the underlying physical mechanisms responsible for the different constitutive laws for cohesive and noncohesive granular media.
The way granular materials response to an applied shear stress is of the utmost relevance to both human activities and natural environment. One of the their most intriguing and less understood behavior, is the stick-instability, whose most dramatic manifestation are earthquakes, ultimately governed by the dynamics of rocks and debris jammed within the fault gauge. Many of the features of earthquakes, i.e. intermittency, broad times and energy scale involved, are mimicked by a very simple experimental set-up, where small beads of glass under load are slowly sheared by an elastic medium. Analyzing data from long lasting experiments, we identify a critical dynamical regime, that can be related to known theoretical models used for crackling-noise phenomena. In particular, we focus on the average shape of the slip velocity, observing a breakdown of scaling: while small slips show a self-similar shape, large does not, in a way that suggests the presence of subtle inertial effects within the granular system. In order to characterise the crossover between the two regimes, we investigate the frictional response of the system, which we trat as a stochastic quantity. Computing different averages, we evidence a weakening effect, whose Stribeck threshold velocity can be related to the aforementioned breaking of scaling.
For vertical velocity field $v_{rm z} (r,z;R)$ of granular flow through an aperture of radius $R$, we propose a size scaling form $v_{rm z}(r,z;R)=v_{rm z} (0,0;R)f (r/R_{rm r}, z/R_{rm z})$ in the region above the aperture. The length scales $R_{rm r}=R- 0.5 d$ and $R_{rm z}=R+k_2 d$, where $k_2$ is a parameter to be determined and $d$ is the diameter of granule. The effective acceleration, which is derived from $v_{rm z}$, follows also a size scaling form $a_{rm eff} = v_{rm z}^2(0,0;R)R_{rm z}^{-1} theta (r/R_{rm r}, z/R_{rm z})$. For granular flow under gravity $g$, there is a boundary condition $a_{rm eff} (0,0;R)=-g$ which gives rise to $v_{rm z} (0,0;R)= sqrt{ lambda g R_{rm z}}$ with $lambda=-1/theta (0,0)$. Using the size scaling form of vertical velocity field and its boundary condition, we can obtain the flow rate $W =C_2 rho sqrt{g } R_{rm r}^{D-1} R_{rm z}^{1/2} $, which agrees with the Beverloo law when $R gg d$. The vertical velocity fields $v_z (r,z;R)$ in three-dimensional (3D) and two-dimensional (2D) hoppers have been simulated using the discrete element method (DEM) and GPU program. Simulation data confirm the size scaling form of $v_{rm z} (r,z;R)$ and the $R$-dependence of $v_{rm z} (0,0;R)$.
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.
Based on discrete element method simulations, we propose a new form of the constitution equation for granular flows independent of packing fraction. Rescaling the stress ratio $mu$ by a power of dimensionless temperature $Theta$ makes the data from a wide set of flow geometries collapse to a master curve depending only on the inertial number $I$. The basic power-law structure appears robust to varying particle properties (e.g. surface friction) in both 2D and 3D systems. We show how this rheology fits and extends frameworks such as kinetic theory and the Nonlocal Granular Fluidity model.
We study experimentally the fracture mechanisms of a model cohesive granular medium consisting of glass beads held together by solidified polymer bridges. The elastic response of this material can be controlled by changing the cross-linking of the polymer phase, for example. Here we show that its fracture toughness can be tuned over an order of magnitude by adjusting the stiffness and size of the polymer bridges. We extract a well-defined fracture energy from fracture testing under a range of material preparations. This energy is found to scale linearly with the cross-sectional area of the bridges. Finally, X-ray microcomputed tomography shows that crack propagation is driven by adhesive failure of about one polymer bridge per bead located at the interface, along with microcracks in the vicinity of the failure plane. Our findings provide insight to the fracture mechanisms of this model material, and the mechanical properties of disordered cohesive granular media in general.