Search

We perform discrete element simulations of freely cooling, dense granular materials, previously sheared at a constant rate. Particles are identical, frictional spheres interacting via linear springs and dashpots and the solid volume fraction is constant and equal to 60% during both shearing and cooling. We measure the average and the distributions of contacts per particle and the anisotropy of the contact network. We observe that the granular material, at the beginning of cooling, can be shear-jammed, fragile or unjammed. The initial state determines the subsequent evolution of the dense assembly into either an anisotropic solid, an isotropic or an anisotropic fluid, respectively. While anisotropic solids and isotropic fluids rapidly reach an apparent final steady configuration, the microstructure continues to evolve for anisotropic fluids. We explain this with the presence of vortices in the flow field that counteract the randomizing and structure-annihilating effect of collisions. We notice, in accordance with previous findings, that the initial fraction of mechanically stable particles permits to distinguish between shear-jammed, fragile or unjammed states and, therefore, determine beforehand the fate of the freely evolving granular materials. We also find that the fraction of mechanically stable particles is in a one-to-one relation with the average number of contacts per particle. The latter is, therefore, a variable that must be incorporated in continuum models of granular materials, even in the case of unjammed states, where it was widely accepted that the solid volume fraction was sufficient to describe the geometry of the system.

In this paper, the authors analyse numerical and experimental results concerning either dry or saturated granular flows under steady, simple shear conditions. A new constitutive model is introduced, on the base of the mixture theory, according to which granular and liquid phases are considered separately. The constitutive relationship refers to the representative elementary volume and assumes that the mean values of all kinematic variables, of both granular and liquid phases, coincide. For the two phases, a parallel scheme is chosen. As regards the granular contribution, the authors employ an already conceived constitutive model where the critical state concept and the kinetic theory of granular gases are merged, and in which the granular temperature plays the role of state variable for the material. Under saturated conditions, the new model accounts for granular–liquid coupling effects. In fact, the liquid viscosity is a function of granular concentration, whereas the evolution of granular temperature is influenced by the liquid molecular viscosity. The model is validated against numerical results and critically discussed. For sufficiently small values of concentration, the transition from a Newtonian to a Bagnoldian regime, for increasing values of strain rate, is correctly reproduced.

Simple homogeneous shear flows of frictionless, deformable particles are studied by particle simulations at large shear rates and for differently soft, deformable particles. The particle stiffness sets a time-scale that can be used to scale the physical quantities; thus the dimensionless shear rate, i.e. the inertial number (inversely proportional to pressure), can alternatively be expressed as inversely proportional to the square root of the particle stiffness. Asymptotic scaling relations for the field variables pressure, shear stress and granular temperature are inferred from simulations in both fluid and solid regimes, corresponding to unjammed and jammed conditions. Then the limit cases are merged to unique constitutive relations that cover also the transition zone in proximity of jamming. By exploiting the diverging behavior of the scaling laws at the jamming density, we arrive at continuous and differentiable phenomenological constitutive relations for the stresses and the granular temperature as functions of the volume fraction, shear rate, particle stiffness and distance from jamming. In contrast to steady shear flows of hard particles the (shear) stress ratio does not collapse as a function of the inertial number, indicating the need for an additional control parameter. In the range of particle stiffnesses investigated, in the solid regime, only the pressure is rate independent, whereas the shear stress exhibits a slight shear rate- and stiffness-dependency., Comment: 37 pages, 14 figures, submitted to Soft Matter

The numerical simulation of rapid landslides is quite complex mainly because constitutive models capable of simulating the mechanical behaviour of granular materials in the pre-collapse and post-collapse regimes are still missing. The goal of this paper is to introduce a constitutive model capable of capturing the response of dry granular flows from quasi-static to dynamic conditions, in particular when the material experiences a sort of solid-to-fluid phase transition. An ideal assembly of identical spheres under simple shear conditions is considered. In the constitutive model, void ratio and granular temperature have been chosen as state variables, and both shear and normal stresses are computed as the sum of two contributions: the quasi-static one and the collisional one. The former is determined by using a perfect elasto-plastic model including the critical state concept, while the latter is derived from the kinetic theory of granular gases. The evolution of the granular temperature, fundamentally governing the material phase transition, is obtained by imposing the kinetic fluctuating energy balance. The constitutive relationship has been integrated, under both constant pressure and constant volume conditions, and the influence of shear strain rate, initial void ratio and normal pressure on the mechanical response has been investigated.

We use previous results from discrete element simulations of simple shear flows of rigid, identical spheres in the collisional regime to show that the volume fractiondependence of the stresses is singular at the shear rigidity. Here, we identify the shear rigidity, which is a decreasing function of the interparticle friction, as the maximum volume fraction beyond which a random collisional assembly of grains cannot be sheared without developing force chains that span the entire domain. In the framework of extended kinetic theory, i.e., kinetic theory that accounts for the decreasing in the collisional dissipation due to the breaking of molecular chaos at volume fractions larger than 0.49, we also show that the volume fraction-dependence of the correlation length (measure of the velocity correlation) is singular at random close packing, independent of the interparticle friction. The difference in the singularities ensures that the ratio of the shear stress to the pressure at shear rigidity is different from zero even in the case of frictionless spheres: we identify that with the yield stress ratio of granular materials, and we show that the theoretical predictions, once the different singularities are inserted into the functions of extended kinetic theory, are in excellent agreement with the results of numerical simulations.

This work focuses on the mechanical response of dry granular materials under steady, simple shear conditions. In particular, the goal is to obtain a complete rheology able to describe the material behavior within the entire range of concentrations for which the flow can be considered dense. The total stress is assumed to be the linear sum of a frictional and a kinetic component. The frictional and the kinetic contributions are modeled in the context of the critical state theory and the kinetic theory of dense granular gases, respectively; in the latter, the correlated motion among the particles, which is likely to occur at high concentration, is also included. In accordance with recent findings on disordered granular packings, the frictional component of stresses is assumed to vanish when the concentration is below the random loose packing. According to this approach, four nondimensional quantities govern steady, simple shear flows: the concentration, the shear to normal stress ratio, the ratio of the time scales associated with the motion perpendicular and parallel to the flow, and the ratio between the particle stiffness and the normal stress. The present theory allows us to reproduce, in a notable way, both numerical simulations on simple shear flows of disks and physical experiments on incline flows of glass spheres taken from the literature.