Search

It is well-known that granular mixtures that differ in size or shape segregate when sheared. In the past, two mechanisms have been proposed to describe this effect, and it is unclear if both exist. To settle this question, we consider a bidisperse mixture of spheroids of equal volume in a rotating drum, where the two mechanisms are predicted to act in opposite directions. We present the first evidence that there are two \emph{distinct} segregation mechanisms driven by relative \emph{over-stress}. Additionally, we showed that for non-spherical particles, these two mechanisms can act in different directions leading to a competition between the effects of the two. As a result, the segregation intensity varies non-monotonically as a function of $AR$, and at specific points, the segregation direction changes for both prolate and oblate spheroids, explaining the surprising segregation reversal previously reported. Consistent with previous results, we found that the kinetic mechanism is dominant for (almost) spherical particles. Furthermore, for moderate aspect ratios, the kinetic mechanism is responsible for the spherical particles segregation to the periphery of the drum, and the gravity mechanism plays only a minor role. Whereas, at the extreme values of $AR$, the gravity mechanism notably increases and overtakes its kinetic counterpart.

Very recently, To et al.~have experimentally explored granular flow in a cylindrical silo, with a bottom wall that rotates horizontally with respect to the lateral wall \cite{Kiwing2019}. Here, we numerically reproduce their experimental findings, in particular, the peculiar behavior of the mass flow rate $Q$ as a function of the frequency of rotation $f$. Namely, we find that for small outlet diameters $D$ the flow rate increased with $f$, while for larger $D$ a non-monotonic behavior is confirmed. Furthermore, using a coarse-graining technique, we compute the macroscopic density, momentum, and the stress tensor fields. These results show conclusively that changes in the discharge process are directly related to changes in the flow pattern from funnel flow to mass flow. Moreover, by decomposing the mass flux (linear momentum field) at the orifice into two main factors: macroscopic velocity and density fields, we obtain that the non-monotonic behavior of the linear momentum is caused by density changes rather than by changes in the macroscopic velocity. In addition, by analyzing the spatial distribution of the kinetic stress, we find that for small orifices increasing rotational shear enhances the mean kinetic pressure $\langle p^k \rangle$ and the system dilatancy. This reduces the stability of the arches, and, consequently, the volumetric flow rate increases monotonically. For large orifices, however, we detected that $\langle p^k \rangle$ changes non-monotonically, which might explain the non-monotonic behavior of $Q$ when varying the rotational shear.

We analyze the failure process of a two-component system with widely different fracture strength in the framework of a fiber bundle model with localized load sharing. A fraction 0\leq \alpha \leq 1 of the bundle is strong and it is represented by unbreakable fibers, while fibers of the weak component have randomly distributed failure strength. Computer simulations revealed that there exists a critical composition \alpha_c which separates two qualitatively different behaviors: below the critical point the failure of the bundle is brittle characterized by an abrupt damage growth within the breakable part of the system. Above \alpha_c, however, the macroscopic response becomes ductile providing stability during the entire breaking process. The transition occurs at an astonishingly low fraction of strong fibers which can have importance for applications. We show that in the ductile phase the size distribution of breaking bursts has a power law functional form with an exponent \mu=2 followed by an exponential cutoff. In the brittle phase the power law also prevails but with a higher exponent \mu=9/2. The transition between the two phases shows analogies to continuous phase transitions. Analyzing the microstructure of the damage, it was found that at the beginning of the fracture process cracks nucleate randomly, while later on growth and coalescence of cracks dominate which give rise to power law distributed crack sizes., Comment: 8 pages, 9 figures

A loosely packed bed of sand sits precariously on the fence between mechanically stable and flowing states. This has especially strong implications for animals or vehicles needing to navigate sandy environments, which can sink and become stuck in a "dry quicksand" if their weight exceeds the yield stress of this fragile matter. While it is known that the contact stresses in these systems are loaded by gravity, very little is known about the sinking dynamics of objects into loose granular systems under gravitational accelerations different from the Earth's (g). A fundamental understanding of how objects sink in different gravitational environments is not only necessary for successful planetary navigation and engineering, but it can also improve our understanding of celestial impact dynamics and crater geomorphology. Here we perform and explain the first systematic experiments of the sink dynamics of objects into granular media in different gravitational accelerations. By using an accelerating experimental apparatus, we explore gravitational conditions ranging from 0.4g to 1.2g. With the aid of discrete element modeling simulations, we reproduce these results and extend this range to include objects as small as asteroids and as large as Jupiter. Surprisingly, we find that the final sink depth is independent of the gravitational acceleration, an observation with immediate relevance to the design of future extraterrestrial structures land-roving spacecraft. Using a phenomenological equation of motion that includes a gravity-loaded frictional term, we are able to quantitatively explain the experimental and simulation results.

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., Comment: 6 pages, 6 figures

A two-dimensional fracture model where the interaction among elements is modeled by an anisotropic stress-transfer function is presented. The influence of anisotropy on the macroscopic properties of the samples is clarified, by interpolating between several limiting cases of load sharing. Furthermore, the critical stress and the distribution of failure avalanches are obtained numerically for different values of the anisotropy parameter $\alpha$ and as a function of the interaction exponent $\gamma$. From numerical results, one can certainly conclude that the anisotropy does not change the crossover point $\gamma_c=2$ in 2D. Hence, in the limit of infinite system size, the crossover value $\gamma_c=2$ between local and global load sharing is the same as the one obtained in the isotropic case. In the case of finite systems, however, for $\gamma\le2$, the global load sharing behavior is approached very slowly.

We study the effect of strong heterogeneities on the fracture of disordered materials using a fiber bundle model. The bundle is composed of two subsets of fibers, i.e. a fraction 0<\alpha<1 of fibers is unbreakable, while the remaining 1-\alpha fraction is characterized by a distribution of breaking thresholds. Assuming global load sharing, we show analytically that there exists a critical fraction of the components \alpha_c which separates two qualitatively different regimes of the system: below \alpha_c the burst size distribution is a power law with the usual exponent \tau=5/2, while above \alpha_c the exponent switches to a lower value \tau=9/4 and a cutoff function occurs with a diverging characteristic size. Analyzing the macroscopic response of the system we demonstrate that the transition is conditioned to disorder distributions where the constitutive curve has a single maximum and an inflexion point defining a novel universality class of breakdown phenomena.