Self-diffusion in dense granular shear flows

Rheology of cohesive granular materials across multiple dense-flow regimes

We investigate the dense-flow rheology of cohesive granular materials through discrete element simulations of homogeneous, simple shear flows of frictional, cohesive, spherical particles. Dense shear flows of noncohesive granular materials exhibit three regimes: quasistatic, inertial, and intermediate, which persist for cohesive materials as well. It is found that cohesion results in bifurcation of the inertial regime into two regimes: (a) a new rate-independent regime and (b) an inertial regime. Transition from rate-independent cohesive regime to inertial regime occurs when the kinetic energy supplied by shearing is sufficient to overcome the cohesive energy. Simulations reveal that inhomogeneous shear band forms in the vicinity of this transition, which is more pronounced at lower particle volume fractions. We propose a rheological model for cohesive systems that captures the simulation results across all four regimes.

Effect of inertia on sheared disordered solids: critical scaling of avalanches in two and three dimensions

Molecular dynamics simulations with varying damping are used to examine the effects of inertia and spatial dimension on sheared disordered solids in the athermal quasistatic limit. In all cases the distribution of avalanche sizes follows a power law over at least three orders of magnitude in dissipated energy or stress drop. Scaling exponents are determined using finite-size scaling for systems with 10(3)-10(6) particles. Three distinct universality classes are identified corresponding to overdamped and underdamped limits, as well as a crossover damping that separates the two regimes. For each universality class, the exponent describing the avalanche distributions is the same in two and three dimensions. The spatial extent of plastic deformation is proportional to the energy dissipated in an avalanche. Both rise much more rapidly with system size in the underdamped limit where inertia is important. Inertia also lowers the mean energy of configurations sampled by the system and leads to an excess of large events like that seen in earthquake distributions for individual faults. The distribution of stress values during shear narrows to zero with increasing system size and may provide useful information about the size of elemental events in experimental systems. For overdamped and crossover systems the stress variation scales inversely with the square root of the system size. For underdamped systems the variation is determined by the size of the largest events.

Kinematics of monodisperse and bidisperse granular flows in quasi-two-dimensional bounded heaps

Quasi-two-dimensional bounded heap flow is a useful model for many granular flows in industry and nature. It belongs to a family of free surface flows—inclined chute flow, rotating tumbler flow and unbounded heap flow—but differs from the others in that uniform deposition of particles onto the static bed results in the uniform rise of the heap. The kinematics, however, are only partially understood. We performed discrete element method simulations to study granular flows in quasi-two-dimensional bounded heaps. The experimentally validated computational results show a universal functional form for the streamwise velocity profile for both monodisperse and bidisperse systems when velocities and coordinates are scaled by the local surface velocity and the local flowing layer thickness. This holds true regardless of streamwise location, feed rate, particle size distribution and, most surprisingly, the local particle concentration for bidisperse flows. The local surface velocity decreases linearly in the streamwise direction, while the flowing layer thickness remains nearly constant; both quantities depending only on local flow rate and local mean particle diameter. Additionally, the velocity profile normal to the overall flow, which is important in understanding segregation, can be predicted analytically from the streamwise velocity and matches the simulation results.

Granular convection in microgravity

We investigate the role of gravity on convection in a dense granular shear flow. Using a microgravity-modified Taylor-Couette shear cell under the conditions of parabolic flight microgravity, we demonstrate experimentally that secondary, convective-like flows in a sheared granular material are close to zero in microgravity and enhanced under high-gravity conditions, though the primary flow fields are unaffected by gravity. We suggest that gravity tunes the frictional particle-particle and particle-wall interactions, which have been proposed to drive the secondary flow. In addition, the degree of plastic deformation increases with increasing gravitational forces, supporting the notion that friction is the ultimate cause.

Stratification, segregation, and mixing of granular materials in quasi-two-dimensional bounded heaps

Segregation and mixing of granular mixtures during heap formation has important consequences in industry and agriculture. This research investigates three different final particle configurations of bidisperse granular mixtures--stratified, segregated and mixed--during filling of quasi-two-dimensional silos. We consider a large number and wide range of control parameters, including particle size ratio, flow rate, system size, and heap rise velocity. The boundary between stratified and unstratified states is primarily controlled by the two-dimensional flow rate, with the critical flow rate for the transition depending weakly on particle size ratio and flowing layer length. In contrast, the transition from segregated to mixed states is controlled by the rise velocity of the heap, a control parameter not previously considered. The critical rise velocity for the transition depends strongly on the particle size ratio.

Resolving a paradox of anomalous scalings in the diffusion of granular materials

Granular materials do not perform Brownian motion, yet diffusion can be observed in such systems when agitation causes inelastic collisions between particles. It has been suggested that axial diffusion of granular matter in a rotating drum might be "anomalous" in the sense that the mean squared displacement of particles follows a power law in time with exponent less than unity. Further numerical and experimental studies have been unable to definitively confirm or disprove this observation. We show two possible resolutions to this apparent paradox without the need to appeal to anomalous diffusion. First, we consider the evolution of arbitrary (non-point-source) initial data towards the self-similar intermediate asymptotics of diffusion by deriving an analytical expression for the instantaneous collapse exponent of the macroscopic concentration profiles. Second, we account for the concentration-dependent diffusivity in bidisperse mixtures, and we give an asymptotic argument for the self-similar behavior of such a diffusion process, for which an exact self-similar analytical solution does not exist. The theoretical arguments are verified through numerical simulations of the governing partial differential equations, showing that concentration-dependent diffusivity leads to two intermediate asymptotic regimes: one with an anomalous scaling that matches the experimental observations for naturally polydisperse granular materials, and another with a "normal" diffusive scaling (consistent with a "normal" random walk) at even longer times.

Particle dynamics and effective temperature of jammed granular matter in a slowly sheared three-dimensional Couette cell

We report experimental measurements of particle dynamics on slowly sheared granular matter in a three-dimensional Couette cell. A closely packed ensemble of transparent spherical beads is confined by an external pressure and filled with fluid to match both the density and refractive index of the beads. This allows us to track tracer particles embedded in the system and obtain three-dimensional trajectories [r(t),theta(t),z(t)] as a function of time. We study the probability distribution function of the vertical and radial displacements, finding Gaussian and exponential distributions, respectively. For slow shear rates, the mean-square fluctuations in all three directions are found to be dependent only on the angular displacement of the Couette cell, Delta theta e, (Delta z 2) approximately Delta theta e, (Delta r2) approximately Delta theta e alpha, Delta theta 2 approximately Delta theta e beta, where alpha and beta are constants. With Delta theta e proportional to the time between measurements, the values of the constants, alpha and beta , are found to be subdiffusive and superdiffusive, respectively. ThFe linear relation between (Delta z 2) and angular displacement implies a diffusive process, from which we can calculate an "effective temperature," T eff, in the vertical direction, through a fluctuation-dissipation relation. It is of interest to determine whether these systems can be described by analogous equilibrium statistical mechanics concepts such as "effective temperature" and "compactivity." By studying the dynamics of tracer particles, we find the effective temperature defined by the Stokes-Einstein relation to be independent of the tracer particle characteristic features, such as density and size, and dependent only on the packing density of the system. For slow shear rate, both the diffusivity and mobility of tracer particles are proportional to the shear rate, giving rise to a constant effective temperature, characteristic of the jammed system. We finally discuss the significance of the existence of T eff for a statistical mechanics formulation of granular matter.

Experimental measures of affine and nonaffine deformation in granular shear

Through 2D granular Couette flow experiments, we probe failure and deformation of disordered solids under shear. Shear produces a mean azimuthal flow, smooth affine deformations, and irreversible so-called nonaffine particle displacements. We find that these processes are all of comparable magnitude and depend on the local shear rate. We compute the parameter of Falk and Langer characterizing nonaffine motion, Dmin2, and find that it is reasonably well described in terms of collections of single particles making locally nearly isotropic random steps, delta ri. Distributions for single particle nonaffine displacements, delta ri, satisfy P1(delta ri) proportional, variantexp[-|delta ri/Delta r|alpha] (alpha < or approximately 2).

Superdiffusion and non-Gaussian statistics in a driven-dissipative 2D dusty plasma

Anomalous diffusion and non-Gaussian statistics are detected experimentally in a two-dimensional driven-dissipative system. A single-layer dusty plasma suspension with a Yukawa interaction and frictional dissipation is heated with laser radiation pressure to yield a structure with liquid ordering. Analyzing the time series for mean-square displacement, superdiffusion is detected at a low but statistically significant level over a wide range of temperatures. The probability distribution function fits a Tsallis distribution, yielding q, a measure of nonextensivity for non-Gaussian statistics.

Anomalous diffusion in silo drainage

The silo discharge process is studied by molecular dynamics simulations. The development of the velocity profile and the probability density function for the displacements in the horizontal and vertical axis are obtained. The PDFs obtained at the beginning of the discharge reveal non-Gaussian statistics and superdiffusive behaviors. When the stationary flow is developed, the PDFs at shorter temporal scales are non-Gaussian too. For big orifices a well-defined transition between ballistic and diffusive regime is observed. In the case of a small outlet orifice, no well-defined transition is observed. We use a nonlinear diffusion equation introduced in the framework of non-extensive thermodynamics in order to describe the movements of the grains. The solution of this equation gives a well-defined relationship (gamma = 2/(3 - q)) between the anomalous diffusion exponent gamma and the entropic parameter q introduced by the non-extensive formalism to fit the PDF of the fluctuations.

Asymmetric response of a jammed plastic bead raft

The successful development of an effective temperature would be an important step in the application of statistical mechanics principles to systems driven far from equilibrium. One direction that has shown promise is the use of fluctuation-dissipation relations. However, driven systems break time-reversal symmetry, and understanding the implications of this for fluctuation-dissipation relations is a critical step in developing effective temperatures. Here we study the response function in a driven system of plastic beads as a function of the density in order to elucidate the generality of the use of fluctuation-dissipation relations. We find that even when a linear response is observed, the time scale of the response is dependent on the direction of the applied stress.

Experimental measurement of an effective temperature for jammed granular materials

A densely packed granular system is an example of an out-of-equilibrium system in the jammed state. It has been a longstanding problem to determine whether this class of systems can be described by concepts arising from equilibrium statistical mechanics, such as an effective temperature and compactivity. The measurement of the effective temperature is realized in the laboratory by slowly shearing a closely packed ensemble of spherical beads confined by an external pressure in a Couette geometry. All of the probe particles considered in this study, independent of their characteristic features, equilibrate at the same temperature, given by the packing density of the system.

Diffusion and mobility in a stirred dense granular material

We describe a probe of diffusivity (D) and mobility (B) for a dense 2D granular system. We introduce random motion by stirring, and characterize D by particle tracking. To measure B we measure the force needed to push a particle through the medium at fixed velocity, v, using three sizes of tracer particle. We find simple Brownian diffusion, but B depends strongly on v because the force needed to push a tracer through a sample is nearly independent of v. Data for D/B depend on the tracer particle size.

Transients in sheared granular matter

As dense granular materials are sheared, a shear band and an anisotropic force network form. The approach to steady-state behavior depends on the history of the packing and the existing force and contact network. We present experiments on shearing of dense granular matter in a 2D Couette geometry in which we probe the history and evolution of shear bands by measuring particle trajectories and stresses during transients. We find that when shearing is stopped and restarted in the same direction, steady-state behavior is immediately reached, in agreement with the typical assumption that the system is quasistatic. Although some relaxation of the force network is observed when shearing is stopped, quasistatic behavior is maintained because the contact network remains essentially unchanged. When the direction of shear is reversed, a transient occurs in which stresses initially decrease, changes in the force network reach further into the bulk, and particles far from the wheel become more mobile. This occurs because the force network is fragile to changes transverse to the force network established under previous shear; particles must rearrange before becoming jammed again, thereby providing resistance to shear in the reversed direction. The strong force network is re-established after displacing the shearing surface approximately equal 3d, where d is the mean grain diameter. Steady-state velocity profiles are reached after a shear of < or approximately equal 30 d. Particles immediately outside of the shear band move on average less than 1 diameter before becoming jammed again. We also examine particle rotation during this transient and find that mean particle spin decreases during the transient, which is related to the fact that grains are not interlocked as strongly.