Self-diffusion of grains in a two-dimensional vibrofluidized bed

Positron Emission Particle Tracking of Granular Flows

Positron emission particle tracking (PEPT) is a noninvasive technique capable of imaging the three-dimensional dynamics of a wide variety of powders, particles, grains, and/or fluids. The PEPT technique can track the motion of particles with high temporal and spatial resolution and can be used to study various phenomena in systems spanning a broad range of scales, geometries, and physical states. We provide an introduction to the PEPT technique, an overview of its fundamental principles and operation, and a brief review of its application to a diverse range of scientific and industrial systems.

Nuclear magnetic resonance measurements of velocity distributions in an ultrasonically vibrated granular bed

We report the results of nuclear magnetic resonance imaging experiments on granular beds of mustard grains fluidized by vertical vibration at ultrasonic frequencies. The variation of both granular temperature and packing fraction with height was measured within the three-dimensional cell for a range of vibration frequencies, amplitudes and numbers of grains. Small increases in vibration frequency were found--contrary to the predictions of classical 'hard-sphere' expressions for the energy flux through a vibrating boundary--to result in dramatic reductions in granular temperature. Numerical simulations of the grain-wall interactions, using experimentally determined Hertzian contact stiffness coefficients, showed that energy flux drops significantly as the vibration period approaches the grain-wall contact time. The experiments thus demonstrate the need for new models for 'soft-sphere' boundary conditions at ultrasonic frequencies.

Long-time tails and cage effect in driven granular fluids

We study the velocity autocorrelation function of a driven granular fluid in the stationary state in three dimensions. As the critical volume fraction of the glass transition in the corresponding elastic system is approached, we observe pronounced cage effects in the velocity autocorrelation function as well as a strong decrease of the diffusion constant, depending on the inelasticity. At moderate densities the velocity autocorrelation function is shown to decay algebraically in time, like t(-3/2), if momentum is conserved locally, and like t(-1), if momentum is not conserved by the driving. A simple scaling argument supports the observed long-time tails.

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.

NMR measurements and hydrodynamic simulations of phase-resolved velocity distributions within a three-dimensional vibrofluidized granular bed

We report the results of nuclear magnetic resonance imaging experiments on vertically vibrated granular beds of mustard grains. A novel spin-echo velocity profiling technique was developed that allows granular temperature, mean velocity and packing fraction distributions within the three-dimensional cell to be measured as a function of both vertical position and vibration phase. Bimodal velocity distributions were observed at certain portions of the vibration cycle, and in general the ability to acquire time-resolved data demonstrated the significant distortions to the velocity distributions and the systematic errors in calculated temperature distributions that may arise with time-averaged measurements. The experimental behaviour was compared with predictions from a time-varying one-dimensional hydrodynamic model using the experimental parameters as input to the code. In both cases, damping of longitudinal sound waves was linked to significant volume heating effects, which contrasts with the usual heat transport mechanism (i.e. diffusion from the boundaries) currently assumed in most steady-state models. This leads to a new explanation for the counterintuitive upturn in granular temperature in vibrofluidized granular beds, based on amplification and damping of sound waves in the high-altitude region.

Fluidization of a vertically vibrated two-dimensional hard sphere packing: a granular meltdown

We report measurements of the fluidization process in vertically vibrated two-dimensional granular packings. An initially close packed granular bed is exposed to sinusoidal container oscillations with gradually increasing amplitude. At first the particles close to the free surface become mobile. When a critical value of the forcing strength is reached the remaining crystal suddenly breaks up and the bed fluidizes completely. This transition leads to discontinuous changes in the density distribution and in the root mean square displacement of the individual particles. Likewise the vertical center of mass coordinate increases by leaps and bounds at the transition. It turns out that the maximum container velocity v0 is the crucial driving parameter determining the state of a fully fluidized system. For particles of various sizes the transition to full fluidization occurs at the same value of v 2 0/gd, where d is the particle diameter and g is the gravitational acceleration. A discontinuous fluidization transition is only observed when the particles are highly elastic.

NMR experiments on a three-dimensional vibrofluidized granular medium

A three-dimensional granular system fluidized by vertical container vibrations was studied using pulsed field gradient NMR coupled with one-dimensional magnetic resonance imaging. The system consisted of mustard seeds vibrated vertically at 50 Hz, and the number of layers N(l)<or=4 was sufficiently low to achieve a nearly time-independent granular fluid. Using NMR, the vertical profiles of density and granular temperature were directly measured, along with the distributions of vertical and horizontal grain velocities. The velocity distributions showed modest deviations from Maxwell-Boltzmann statistics, except for the vertical velocity distribution near the sample bottom, which was highly skewed and non-Gaussian. Data taken for three values of N(l) and two dimensionless accelerations Gamma=15,18 were fitted to a hydrodynamic theory, which successfully models the density and temperature profiles away from the vibrating container bottom. A temperature inversion near the free upper surface is observed, in agreement with predictions based on the hydrodynamic parameter micro which is nonzero only in inelastic systems.

Self-diffusion in dense granular shear flows

Diffusivity is a key quantity in describing velocity fluctuations in granular materials. These fluctuations are the basis of many thermodynamic and hydrodynamic models which aim to provide a statistical description of granular systems. We present experimental results on diffusivity in dense, granular shear flows in a two-dimensional Couette geometry. We find that self-diffusivities D are proportional to the local shear rate gamma; with diffusivities along the direction of the mean flow approximately twice as large as those in the perpendicular direction. The magnitude of the diffusivity is D approximately gamma;a(2), where a is the particle radius. However, the gradient in shear rate, coupling to the mean flow, and strong drag at the moving boundary lead to particle displacements that can appear subdiffusive or superdiffusive. In particular, diffusion appears to be superdiffusive along the mean flow direction due to Taylor dispersion effects and subdiffusive along the perpendicular direction due to the gradient in shear rate. The anisotropic force network leads to an additional anisotropy in the diffusivity that is a property of dense systems and has no obvious analog in rapid flows. Specifically, the diffusivity is suppressed along the direction of the strong force network. A simple random walk simulation reproduces the key features of the data, such as the apparent superdiffusive and subdiffusive behavior arising from the mean velocity field, confirming the underlying diffusive motion. The additional anisotropy is not observed in the simulation since the strong force network is not included. Examples of correlated motion, such as transient vortices, and Lévy flights are also observed. Although correlated motion creates velocity fields which are qualitatively different from collisional Brownian motion and can introduce nondiffusive effects, on average the system appears simply diffusive.

Collision statistics of driven granular materials

We present an experimental investigation of the statistical properties of spherical granular particles on an inclined plane that are excited by an oscillating side wall. The data is obtained by high-speed imaging and particle tracking techniques. We identify all particles in the system and link their positions to form trajectories over long times. Thus, we identify particle collisions to measure the effective coefficient of restitution and find a broad distribution of values for the same impact angles. We find that the energy inelasticity can take on values greater than one, which implies that the rotational degrees of freedom play an important role in energy transfer. We also measure the distance and the time between collision events in order to directly determine the distribution of path lengths and the free times. These distributions are shown to deviate from expected theoretical forms for elastic spheres, demonstrating the inherent clustering in this system. We describe the data with a two-parameter fitting function and use it to calculate the mean free path and collision time. We find that the ratio of these values is consistent with the average velocity. The velocity distributions are observed to be strongly non-Gaussian and do not demonstrate any apparent universal behavior. We report the scaling of the second moment, which corresponds to the granular temperature, and higher order moments as a function of distance from the driving wall. Additionally, we measure long-time correlation functions in both space and in the velocities to probe diffusion in a dissipative gas.

Transport theory of granular swarms

The transport of trace granular gas (swarm) in a carrier granular fluid is studied by means of the Boltzmann-Lorentz kinetic equation. Time-dependent perturbation theory is used to follow the evolution of the granular swarm from an arbitrary initial distribution. A nonhydrodynamic extension of the diffusion equation is derived, with transport coefficients that are time dependent and implicitly depend on the wave vector. Transport coefficients of any order are obtained as velocity moments of the solutions of the corresponding kinetic equations derived from the Boltzmann-Lorentz equation. For the special case of the initial distribution of swarm particles, transport coefficients are identified as time derivatives of the moments of the number density. Finally the granular particle transport theory is extended by the introduction of the concept of non-particle-conserving collisions.

Coexistence of two granular temperatures in binary vibrofluidized beds

An investigation into the granular temperature distributions of a binary vibrofluidized granular bed has been conducted using positron emission particle tracking. By repeating each experiment with the tracer selected in turn from the two size components, the granular temperature and packing fraction distributions for each phase were determined. It was found that, for a range of size fractions, the granular temperature of the larger particles was higher than that of the smaller diameter grains, a result which was supported by a simple theoretical analysis based on the steady state energy equation.

Measurements of grain motion in a dense, three-dimensional granular fluid

We have used an NMR technique to measure the short-time, three-dimensional displacement of grains in a system of mustard seeds vibrated vertically at 15 g. The technique averages over a time interval in which the grains move ballistically, giving a direct measurement of the granular temperature profile. The dense, lower portion of the sample is well described by a recent hydrodynamic theory for inelastic hard spheres. Near the free upper surface the mean free path is longer than the particle diameter and the hydrodynamic description fails.

Numerical solution of the Smoluchowski equation for a vibrofluidized granular bed

A stochastic approach, similar to that used to describe Brownian motion, was used to model the displacement probability of grains in a three-dimensional vibrofluidized granular bed. As neither an analytical description nor measurements of the diffusion coefficients were available, the governing partial differential equation, namely, the Smoluchowski equation, was solved numerically using an iterative procedure, modifying the granular temperature profile at each step. The results of this stochastic model were compared to experimental measurements of the displacement probability density made using positron emission particle tracking. The results indicate that methods based on hard elastic systems such as the Smoluchowski equation are appropriate to granular systems, particularly over timescales greater than the mean collision time.

Granular temperature profiles in three-dimensional vibrofluidized granular beds

The motion of grains in a three-dimensional vibrofluidized granular bed has been measured using the technique of positron emission particle tracking, to provide three-dimensional packing fraction and granular temperature distributions. The mean square fluctuation velocity about the mean was calculated through analysis of the short time mean squared displacement behavior, allowing measurement of the granular temperature at packing fractions of up to eta approximately 0.15. The scaling relationship between the granular temperature, the number of layers of grains, and the base velocity was determined. Deviations between the observed scaling exponents and those predicted by recent theories are attributed to the influence of dissipative grain-sidewall collisions.

Convection in highly fluidized three-dimensional granular beds

Free, buoyancy-driven convection has been observed experimentally in three-dimensional highly fluidized granular flows for the first time. Positron emission particle tracking was used to determine the position of a tracer grain in a vibrofluidized bed, from which packing fraction distributions as well as the velocity fields could be determined. The convection rolls, although small compared to the magnitude of velocity fluctuations (<5%), were consistently observed for a range of grain numbers and shaker amplitudes. Density variations are a signature of free convection and, with negative temperature gradients also present, were interpreted as the mechanism by which the convection rolls were initiated.

Self-diffusion in a gas-fluidized bed of fine powder

We have investigated the self-diffusion in a stable gas-fluidized bed of fine powder. Two regimes have been observed: for gas velocities v(g) above the minimum fluidization velocity v(m) and below a critical gas velocity v(c) smaller than the minimum bubbling velocity v(b) the powder does not mix. Experimental measurements show the existence of yield stresses in this regime which are responsible for the static behavior of the bed. For v(g)>v(c) the yield stress vanishes; the bed behaves like a fluid and displays a diffusive dynamics. In this region we have found that the diffusion coefficient D increases with gas velocity until the bed expansion approaches its maximum value.

Statistical mechanics of fluidized granular media: short-range velocity correlations

A statistical mechanical study of fluidized granular media is presented. Using a special energy injection mechanism, homogeneous fluidized stationary states are obtained. Molecular dynamics simulations and theoretical analysis of the inelastic hard-disk model show that there is a large asymmetry in the two-particle distribution function between pairs that approach and separate. Large velocity correlations appear in the postcollisional states due to the dissipative character of the collision rule. These correlations can be well-characterized by a state dependent pair correlation function at contact. It is also found that velocity correlations are present for pairs that are about to collide. Particles arrive at collisions with a higher probability that their velocities are parallel rather than antiparallel. These dynamical correlations lead to a decrease of the pressure and of the collision frequency as compared to their Enskog values. A phenomenological modified equation of state is presented.

Single-particle motion in three-dimensional vibrofluidized granular beds

A technique to probe the interior of three-dimensional dynamic granular systems is presented. Positron emission particle tracking (PEPT) allows a single tracer particle to be followed around a three dimensional vibrofluidized granular bed for periods up to six hours. At present the technique is able to resolve the position of the grains to +/-4 mm, with an average temporal resolution of about 7 ms. Packing fraction profiles are calculated by making use of the ergodicity of the system, and granular temperature profiles are obtained, in the dilute case, from the short time behavior of the mean squared displacement. At longer times, the mean squared displacement shows a range of behavior which can be explained by the presence of strong gradients in the packing fraction. Convection currents were observed, but were sufficiently small in magnitude to be ignored during the analysis of grain motion. The system was modeled using the Smoluchowski equation, which was solved numerically, and the results compared with the experimentally determined displacement probability density functions. Good agreement between experiment and numerical results was achieved using Brownian motion relationships modified to accommodate differences between granular systems and thermal systems.

Nonlinear analysis of the shearing instability in granular gases

It is known that a finite-size homogeneous granular fluid develops a hydrodynamiclike instability when dissipation crosses a threshold value. This instability is analyzed in terms of modified hydrodynamic equations: first, a source term is added to the energy equation which accounts for the energy dissipation at collisions and the phenomenological Fourier law is generalized according to previous results. Second, a rescaled time formalism is introduced that maps the homogeneous cooling state into a nonequilibrium steady state. A nonlinear stability analysis of the resulting equations is done which predicts the appearance of flow patterns. A stable modulation of density and temperature is produced that does not lead to clustering. Also a global decrease of the temperature is obtained, giving rise to a decrease of the collision frequency and dissipation rate. Good agreement with molecular dynamics simulations of inelastic hard disks is found for low dissipation.