Order parameter description of stationary partially fluidized shear granular flows

Relaxation-type nonlocal inertial-number rheology for dry granular flows

We propose a constitutive model to describe the nonlocality, hysteresis, and several flow features of dry granular materials. Taking the well-known inertial number I as a measure of sheared-induced local fluidization, we derive a relaxation model for I according to the evolution of microstructure during avalanche and dissipation processes. The model yields a nonmonotonic flow law for a homogeneous flow, accounting for hysteretic solid-fluid transition and intermittency in quasistatic flows. For an inhomogeneous flow, the model predicts a generalized Bagnold shear stress revealing the interplay of two microscopic nonlocal mechanisms: collisions among correlated structures and the diffusion of fluidization within the structures. In describing a uniform flow down an incline, the model reproduces the hysteretic starting and stopping heights and the Pouliquen flow rule for mean velocity. Moreover, a dimensionless parameter reflecting the nonlocal effect on the flow is discovered, which controls the transition between Bagnold and creeping flow dynamics.

Microscopic Description of the Granular Fluidity Field in Nonlocal Flow Modeling

A recent granular rheology based on an implicit "granular fluidity" field has been shown to quantitatively predict many nonlocal phenomena. However, the physical nature of the field has not been identified. Here, the granular fluidity is found to be a kinematic variable given by the velocity fluctuation and packing fraction. This is verified with many discrete element simulations, which show that the operational fluidity definition, solutions of the fluidity model, and the proposed microscopic formula all agree. Kinetic theoretical and Eyring-like explanations shed insight into the obtained form.

Unified rheology of vibro-fluidized dry granular media: From slow dense flows to fast gas-like regimes

Granular media take on great importance in industry and geophysics, posing a severe challenge to materials science. Their response properties elude known soft rheological models, even when the yield-stress discontinuity is blurred by vibro-fluidization. Here we propose a broad rheological scenario where average stress sums up a frictional contribution, generalizing conventional μ(I)-rheology, and a kinetic collisional term dominating at fast fluidization. Our conjecture fairly describes a wide series of experiments in a vibrofluidized vane setup, whose phenomenology includes velocity weakening, shear thinning, a discontinuous thinning transition, and gaseous shear thickening. The employed setup gives access to dynamic fluctuations, which exhibit a broad range of timescales. In the slow dense regime the frequency of cage-opening increases with stress and enhances, with respect to μ(I)-rheology, the decrease of viscosity. Diffusivity is exponential in the shear stress in both thinning and thickening regimes, with a huge growth near the transition.

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.

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.

Friction and the oscillatory motion of granular flows

This contribution reports on numerical simulations of two-dimensional granular flows on erodible beds. The broad aim is to investigate whether simple flows of model granular matter exhibit spontaneous oscillatory motion in generic flow conditions, and in this case, whether the frictional properties of the contacts between grains may affect the existence or the characteristics of this oscillatory motion. The analysis of different series of simulations shows that the flow develops an oscillatory motion with a well-defined frequency which increases like the inverse of the velocity's square root. We show that the oscillation is essentially a surface phenomenon. The amplitude of the oscillation is higher for lower volume fractions and can thus be related to the flow velocity and grains' friction properties. The study of the influence of the periodic geometry of the simulation cell shows no significant effect. These results are discussed in relation to sonic sands.

Effect of noise on solid-to-liquid transition in small granular systems under shear

The effect of noise on the solid-to-liquid transition of a dense granular assembly under planar shear is studied numerically using soft-particle molecular dynamics simulations in two dimensions. We focus on small systems in a thin planar Couette cell, examining the bistable region while increasing shear, with varying amounts of random noise, and determine statistics of the shear required for fluidization. In the absence of noise, the threshold value of the shear stress depends on the preparation of the system and has a broad distribution. However, adding force fluctuations both lowers the mean threshold value of the shear stress and decreases its variability. This behavior is interpreted as thermoactivated escape through a fluctuating barrier.

Spreading of a granular droplet

The influence of controlled vibrations on the granular rheology is investigated in a specifically designed experiment in which a granular film spreads under the action of horizontal vibrations. A nonlinear diffusion equation is derived theoretically that describes the evolution of the deposit shape. A self-similar parabolic shape (the "granular droplet") and a spreading dynamics are predicted that both agree quantitatively with the experimental results. The theoretical analysis is used to extract effective friction coefficients between the base and the granular layer under sustained and controlled vibrations. A shear thickening regime characteristic of dense granular flows is evidenced at low vibration energy, both for glass beads and natural sand. Conversely, shear thinning is observed at high agitation.

Hydrodynamics of fluid-solid coexistence in dense shear granular flow

We consider dense rapid shear flow of inelastically colliding hard disks. Navier-Stokes granular hydrodynamics is applied accounting for the recent finding that shear viscosity diverges at a lower density than the rest of the constitutive relations. New interpolation formulas for constitutive relations between dilute and dense cases are proposed and justified in molecular dynamics (MD) simulations. A linear stability analysis of the uniform shear flow is performed and the full phase diagram is presented. It is shown that when the inelasticity of particle collision becomes large enough, the uniform sheared flow gives way to a two-phase flow, where a dense "solidlike" striped cluster is surrounded by two fluid layers. The results of the analysis are verified in event-driven MD simulations, and a good agreement is observed.

Influences of the interstitial liquid on segregation patterns of granular slurries in a rotating drum

Granular mixtures immersed in a liquid (slurries) show segregation dynamics which are quantitatively and qualitatively different from those of dry systems. The principal mechanisms of the segregation dynamics in slurries, as well as the relevant material parameters that must be taken into account in a dynamic description are not sufficiently understood so far. We investigate experimentally the influence of the viscosity of the interstitial liquid on the coarsening of axial segregation patterns in a horizontally rotating mixer. It is found that not only the characteristic time scales but also fundamental structural features of these patterns are influenced by the viscous properties of the liquid component.

Towards a theoretical picture of dense granular flows down inclines

Unlike most fluids, granular materials include coexisting solid, liquid or gaseous regions, which produce a rich variety of complex flows. Dense flows down inclines preserve this complexity but remain simple enough for detailed analysis. In this review we survey recent advances in this rapidly evolving area of granular flow, with the aim of providing an organized, synthetic review of phenomena and a characterization of the state of understanding. The perspective that we adopt is influenced by the hope of obtaining a theory for dense, inclined flows that is based on assumptions that can be tested in physical experiments and numerical simulations, and that uses input parameters that can be independently measured. We focus on dense granular flows over three kinds of inclined surfaces: flat-frictional, bumpy-frictional and erodible. The wealth of information generated by experiments and numerical simulations for these flows has led to meaningful tests of relatively simple existing theories.

Rapid two-dimensional imaging of bubbles and slugs in a three-dimensional, gas-solid, two-phase flow system using ultrafast magnetic resonance

Ultrafast magnetic resonance has been applied to measure the geometry of bubbles and slugs in a three-dimensional gas-solid two-phase flow. A bed of particles of diameter 0.5 mm were fluidized with gas velocities in the range of 0.08-0.26 m/s. Bubbles were imaged in transverse as well as vertical planes with an acquisition time of down to 25 ms and a spatial resolution down to 1.7 mm. Owing to the ultrafast character of these measurements, it is not only possible to evaluate correlations, e.g., for the bubble diameter, but also evaluate models of complex hydrodynamic phenomena, such as the splitting and coalescence of bubbles.

Coarsening of axial segregation patterns of slurries in a horizontally rotating drum

Segregation structures of granular mixtures in rotating drums represent classical examples of pattern formation in granular material. We investigate the coarsening of axial segregation patterns of slurries in a long horizontally rotating cylinder. The dynamics and the three-dimensional geometry of the segregation structures are analyzed with optical methods and nuclear magnetic resonance imaging. Previous studies have mainly considered global statistical features of the pattern dynamics. In order to get insight into driving mechanisms for the coarsening process, we focus on the details of the dissolution of individual bands. We treat the coarsening as a consequence of interactions of adjacent bands in the pattern, which are determined by their geometrical relations. In addition to initially homogeneous mixtures, which evolve to spontaneously formed patterns, we study the evolution of specially prepared simple initial states. The role of the three-dimensional geometry of the axial core in the dissolution process of segregation bands is demonstrated. Relations between geometry and dynamic processes are established, which may help to find the correct microscopic models for the coarsening mechanism.

Shear-induced crystallization of a dense rapid granular flow: hydrodynamics beyond the melting point

We investigate shear-induced crystallization in a very dense flow of monodisperse inelastic hard spheres. We consider a steady plane Couette flow under constant pressure and neglect gravity. We assume that the granular density is greater than the melting point of the equilibrium phase diagram of elastic hard spheres. We employ a Navier-Stokes hydrodynamics with constitutive relations all of which (except the shear viscosity) diverge at the crystal-packing density, while the shear viscosity diverges at a smaller density. The phase diagram of the steady flow is described by three parameters: an effective Mach number, a scaled energy loss parameter, and an integer number m: the number of half-oscillations in a mechanical analogy that appears in this problem. In a steady shear flow the viscous heating is balanced by energy dissipation via inelastic collisions. This balance can have different forms, producing either a uniform shear flow or a variety of more complicated, nonlinear density, velocity, and temperature profiles. In particular, the model predicts a variety of multilayer two-phase steady shear flows with sharp interphase boundaries. Such a flow may include a few zero-shear (solidlike) layers, each of which moving as a whole, separated by fluidlike regions. As we are dealing with a hard sphere model, the granulate is fluidized within the "solid" layers: the granular temperature is nonzero there, and there is energy flow through the boundaries of the solid layers. A linear stability analysis of the uniform steady shear flow is performed, and a plausible bifurcation diagram of the system, for a fixed m, is suggested. The problem of selection of m remains open.

Transverse instability of avalanches in granular flows down an incline

Avalanche experiments on an erodible substrate are analyzed using the "partial fluidization" model of dense granular flows. The model identifies a family of propagating solitonlike avalanches with shape and velocity controlled by the inclination angle and the depth of the substrate. At high inclination angles, the solitons display a transverse instability, followed by coarsening and fingering similar to recent experimental observation. A primary cause for the transverse instability is directly related to the dependence of the soliton velocity on the granular mass trapped in the avalanche.

Real-time measurement of bubbling phenomena in a three-dimensional gas-fluidized bed using ultrafast magnetic resonance imaging

Ultrafast magnetic resonance imaging has been applied for the first time to measure simultaneously both the rise velocities and coalescence of bubbles, and the dynamics of the solid phase in a gas-solid two-phase flow. Here, we consider the hydrodynamics within a gas-fluidized bed of particles of diameter 0.5 mm contained within a column of internal diameter 50 mm; gas velocities in the range of 0.18-0.54 m/s were studied. The data are of sufficient temporal and spatial resolution that bubble size and the evolution of bubble size and velocity following coalescence events are determined.

Objective decomposition of the stress tensor in granular flows

A model for the stress tensor in granular flows [Volfson, Tsimring, and Aranson, Phys. Rev. Lett. 90, 254301 (2003)] is correctly generalized to an objective form that is independent of the coordinate system. The objective representation correctly models the isotropic and anisotropic parts of the stress tensor, whereas the original model for stress tensor components is dependent on the coordinate system. This general objective form of the model also relaxes the assumption in the original model that the principal axes of the granular stress tensor be coaxial with that of the "fluid" stress tensor. This generalization expands the applicability of the model to a wider class of granular flows. The objective representation is also useful in analyzing other models based on additive decomposition of the stress tensor in granular flows.

Stick-slip dynamics of a granular layer under shear

Stick-slip regime of shear granular flows is studied theoretically and numerically. Numerical experiments are carried out for a thin Couette cell using soft-particle molecular dynamics code in two dimensions. We apply order parameter theory of partially fluidized granular flows and find a good agreement with simulations and experiments by Nasuno et al.

Partially fluidized shear granular flows: continuum theory and molecular dynamics simulations

The continuum theory of partially fluidized shear granular flows is tested and calibrated using two-dimensional soft particle molecular dynamics simulations. The theory is based on the relaxational dynamics of the order parameter that describes the transition between static and flowing regimes of granular material. We define the order parameter as a fraction of static contacts among all contacts between particles. We also propose and verify by direct simulations the constitutive relation based on the splitting of the shear stress tensor into a"fluid part" proportional to the strain rate tensor, and a remaining "solid part." The ratio of these two parts is a function of the order parameter. The rheology of the fluid component agrees well with the kinetic theory of granular fluids even in the dense regime. Based on the hysteretic bifurcation diagram for a thin shear granular layer obtained in simulations, we construct the "free energy" for the order parameter. The theory calibrated using numerical experiments with the thin granular layer is applied to the surface-driven stationary two-dimensional granular flows in a thick granular layer under gravity.