Constitutive relations for steady, dense granular flows

Numerical model for solid-like and fluid-like behavior of granular flows

We propose a constitutive model for both the solid-like and fluid-like behavior of granular materials by decomposing the stress tensor into quasi-static and collisional components. A hypoplastic model is adopted for the solid-like behavior in the quasi-static regime, while the viscous and dilatant behavior in the fluid-like regime is represented by a modified μ ( I ) rheology model. This model effectively captures the transition between solid-like and fluid-like flows. Performance and validation of the proposed model are demonstrated through numerical simulations of element tests.

Merging fluid and solid granular behavior

Simple homogeneous shear flows of frictionless, deformable particles are studied by particle simulations at large shear rates and for differently soft, deformable particles. 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 I (inversely proportional to pressure), can alternatively be expressed as inversely proportional to the square root of 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 the 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 stresses and 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 pressure is rate independent, whereas shear stress exhibits a slight shear rate- and stiffness-dependency.

Enstrophy cascades in two-dimensional dense granular flows

Employing two-dimensional molecular dynamics simulations of dense granular materials under simple shear deformations, we investigate vortex structures of particle rearrangements. Introducing vorticity fields as a measure of local spinning motions of the particles, we observe their heterogeneous distributions, where statistics of vorticity fields exhibit the highly non-Gaussian behavior and typical domain sizes of vorticity fields significantly increase if the system is yielding under quasistatic deformations. In such dense granular flows, a power-law decay of vorticity spectra can be observed at mesoscopic scale, implying anomalous local structures of kinetic energy dissipation. We explain the power-law decay, or enstrophy cascades in dense granular materials, by a dimensional analysis, where the dependence of vorticity spectra not only on the wave number, but also on the shear rate, is well explained. From our dimensional analyses, the scaling of granular temperature and rotational kinetic energy is also predicted.

Anomalous energy cascades in dense granular materials yielding under simple shear deformations

By using molecular dynamics (MD) simulations of dense granular particles in two dimensions, we study turbulent-like structures of their non-affine velocities under simple shear deformations. We find that the spectrum of non-affine velocities, introduced as an analog of the energy spectrum for turbulent flows, exhibits the power-law decay if the system is yielding in a quasi-static regime, where large-scale collective motions and inelastic interactions of granular particles are crucial for the anomalous cascade of kinetic energy. Based on hydrodynamic equations of dense granular materials, which include both kinetic and contact contributions in constitutive relations, we derive a theoretical expression for the spectrum, where a good agreement between the result of MD simulations and theoretical prediction is established over a wide range of length scales.

Depth-averaged analytic solutions for free-surface granular flows impacting rigid walls down inclines

In the present paper, flows of granular materials impacting wall-like obstacles down inclines are described by depth-averaged analytic solutions. Particular attention is paid to extending the existing depth-averaged equations initially developed for frictionless and incompressible fluids down a horizontal plane. The effects of the gravitational acceleration along the slope, and of the retarding acceleration caused by friction as well, are systematically taken into account. The analytic solutions are then used to revisit existing data on rigid walls impacted by granular flows. This approach allows establishing a complete phase diagram for granular flow-wall interaction.

Hydrodynamic instabilities in shear flows of dry cohesive granular particles

We extend the dynamic van der Waals model introduced by A. Onuki [Phys. Rev. Lett., 2005, 94, 054501] to the description of cohesive granular flows under a plane shear to study their hydrodynamic instabilities. By numerically solving the dynamic van der Waals model, we observed various heterogeneous structures of density fields in steady states, where the viscous heating is balanced with the energy dissipation caused by inelastic collisions. Based on the linear stability analysis, we found that the spatial structures are determined by the mean volume fraction, the applied shear rate, and the inelasticity, where the instability is triggered if the system is thermodynamically unstable, i.e. the pressure, p, and the volume fraction, ϕ, satisfy ∂p/∂ϕ < 0.

Steady shearing flows of deformable, inelastic spheres

We extend models for granular flows based on the kinetic theory beyond the critical volume fraction at which a rate-independent contribution to the stresses develops. This involves the incorporation of a measure of the duration of the particle interaction before and after this volume fraction. At volume fractions less than the critical, the stress components contain contributions from momentum exchanged in collisions that are influenced by the particle elasticity. At volume fractions greater than the critical, the stress components contain both static contributions from particle elasticity and dynamic contributions from the momentum transfer associated with the release of elastic energy by the breaking of force chains. A simple expression for the duration of a collision before and after the critical volume fraction permits a smooth transition between the two regimes and predictions for the components of the stress in steady, homogeneous shearing that are in good agreement with the results of numerical simulations. Application of the theory to steady, inhomogeneous flows reproduces the features of such flows seen in numerical simulations and physical experiments.

Applying GSH to a wide range of experiments in granular media

Granular solid hydrodynamics (GSH) is a continuum-mechanical theory for granular media, whose wide range of applicability is shown in this paper. Simple, frequently analytic solutions are related to classic observations at different shear rates, including: i) static stress distribution, clogging; ii) elasto-plastic motion: loading and unloading, approach to the critical state, angle of stability and repose; iii) rapid dense flow: the μ-rheology, Bagnold scaling and the stress minimum; iv) elastic waves, compaction, wide and narrow shear band. Less conventional experiments have also been considered: shear jamming, creep flow, visco-elastic behavior and non-local fluidization. With all these phenomena ordered, related, explained and accounted for, though frequently qualitatively, we believe that GSH may be taken as a unifying framework, providing the appropriate macroscopic vocabulary and mindset that help one coming to terms with the breadth of granular physics.