NMR experiments on a three-dimensional vibrofluidized granular medium

Cooling of a granular gas mixture in microgravity

Granular gases are fascinating non-equilibrium systems with interesting features such as spontaneous clustering and non-Gaussian velocity distributions. Mixtures of different components represent a much more natural composition than monodisperse ensembles but attracted comparably little attention so far. We present the observation and characterization of a mixture of rod-like particles with different sizes and masses in a drop tower experiment. Kinetic energy decay rates during granular cooling and collision rates were determined and Haff's law for homogeneous granular cooling was confirmed. Thereby, energy equipartition between the mixture components and between individual degrees of freedom is violated. Heavier particles keep a slightly higher average kinetic energy than lighter ones. Experimental results are supported by numerical simulations.

Diffusion of intruders in granular suspensions: Enskog theory and random walk interpretation

The Enskog kinetic theory is applied to compute the mean square displacement of impurities or intruders (modeled as smooth inelastic hard spheres) immersed in a granular gas of smooth inelastic hard spheres (grains). Both species (intruders and grains) are surrounded by an interstitial molecular gas (background) that plays the role of a thermal bath. The influence of the latter on the motion of intruders and grains is modeled via a standard viscous drag force supplemented by a stochastic Langevin-like force proportional to the background temperature. We solve the corresponding Enskog-Lorentz kinetic equation by means of the Chapman-Enskog expansion truncated to first order in the gradient of the intruder number density. The integral equation for the diffusion coefficient is solved by considering the first two Sonine approximations. To test these results, we also compute the diffusion coefficient from the numerical solution of the inelastic Enskog equation by means of the direct simulation Monte Carlo method. We find that the first Sonine approximation generally agrees well with the simulation results, although significant discrepancies arise when the intruders become lighter than the grains. Such discrepancies are largely mitigated by the use of the second Sonine approximation, in excellent agreement with computer simulations even for moderately strong inelasticities and/or dissimilar mass and diameter ratios. We invoke a random walk picture of the intruders' motion to shed light on the physics underlying the intricate dependence of the diffusion coefficient on the main system parameters. This approach, recently employed to study the case of an intruder immersed in a granular gas, also proves useful in the present case of a granular suspension. Finally, we discuss the applicability of our model to real systems in the self-diffusion case. We conclude that collisional effects may strongly impact the diffusion coefficient of the grains.

Kinetic Theory and Memory Effects of Homogeneous Inelastic Granular Gases under Nonlinear Drag

We study a dilute granular gas immersed in a thermal bath made of smaller particles with masses not much smaller than the granular ones in this work. Granular particles are assumed to have inelastic and hard interactions, losing energy in collisions as accounted by a constant coefficient of normal restitution. The interaction with the thermal bath is modeled by a nonlinear drag force plus a white-noise stochastic force. The kinetic theory for this system is described by an Enskog-Fokker-Planck equation for the one-particle velocity distribution function. To get explicit results of the temperature aging and steady states, Maxwellian and first Sonine approximations are developed. The latter takes into account the coupling of the excess kurtosis with the temperature. Theoretical predictions are compared with direct simulation Monte Carlo and event-driven molecular dynamics simulations. While good results for the granular temperature are obtained from the Maxwellian approximation, a much better agreement, especially as inelasticity and drag nonlinearity increase, is found when using the first Sonine approximation. The latter approximation is, additionally, crucial to account for memory effects such as Mpemba and Kovacs-like ones.

Visual analysis of density and velocity profiles in dense 3D granular gases

Granular multiparticle ensembles are of interest from fundamental statistical viewpoints as well as for the understanding of collective processes in industry and in nature. Extraction of physical data from optical observations of three-dimensional (3D) granular ensembles poses considerable problems. Particle-based tracking is possible only at low volume fractions, not in clusters. We apply shadow-based and feature-tracking methods to analyze the dynamics of granular gases in a container with vibrating side walls under microgravity. In order to validate the reliability of these optical analysis methods, we perform numerical simulations of ensembles similar to the experiment. The simulation output is graphically rendered to mimic the experimentally obtained images. We validate the output of the optical analysis methods on the basis of this ground truth information. This approach provides insight in two interconnected problems: the confirmation of the accuracy of the simulations and the test of the applicability of the visual analysis. The proposed approach can be used for further investigations of dynamical properties of such media, including the granular Leidenfrost effect, granular cooling, and gas-clustering transitions.

Energy nonequipartition in a collisional model of a confined quasi-two-dimensional granular mixture

A collisional model of a confined quasi-two-dimensional granular mixture is considered to analyze homogeneous steady states. The model includes an effective mechanism to transfer the kinetic energy injected by vibration in the vertical direction to the horizontal degrees of freedom of grains. The set of Enskog kinetic equations for the velocity distribution functions of each component is derived first to analyze the homogeneous state. As in the one-component case, an exact scaling solution is found where the time dependence of the distribution functions occurs entirely through the granular temperature T. As expected, the kinetic partial temperatures T_{i} of each component are different and, hence, energy equipartition is broken down. In the steady state, explicit expressions for the temperature T and the ratio of partial kinetic temperatures T_{i}/T_{j} are obtained by considering Maxwellian distributions defined at the partial temperatures T_{i}. The (scaled) granular temperature and the temperature ratios are given in terms of the coefficients of restitution, the solid volume fraction, the (scaled) parameters of the collisional model, and the ratios of mass, concentration, and diameters. In the case of a binary mixture, the theoretical predictions are exhaustively compared with both direct simulation Monte Carlo and molecular dynamics simulations with a good agreement. The deviations are identified to be originated in the non-Gaussianity of the velocity distributions and on microsegregation patterns, which induce spatial correlations not captured in the Enskog theory.

Velocity Distribution of a Homogeneously Cooling Granular Gas

In contrast to molecular gases, granular gases are characterized by inelastic collisions and require therefore permanent driving to maintain a constant kinetic energy. The kinetic theory of granular gases describes how the average velocity of the particles decreases after the driving is shut off. Moreover, it predicts that the rescaled particle velocity distribution will approach a stationary state with overpopulated high-velocity tails as compared to the Maxwell-Boltzmann distribution. While this fundamental theoretical result was reproduced by numerical simulations, an experimental confirmation is still missing. Using a microgravity experiment that allows the spatially homogeneous excitation of spheres via magnetic fields, we confirm the theoretically predicted exponential decay of the tails of the velocity distribution.

Free Cooling of a Granular Gas of Rodlike Particles in Microgravity

Granular gases as dilute ensembles of particles in random motion are at the basis of elementary structure-forming processes in the Universe, involved in many industrial and natural phenomena, and also excellent models to study fundamental statistical dynamics. The essential difference to molecular gases is the energy dissipation in particle collisions. Its most striking manifestation is the so-called granular cooling, the gradual loss of mechanical energy E(t) in the absence of external excitation. We report an experimental study of homogeneous cooling of three-dimensional granular gases in microgravity. The asymptotic scaling E(t)∝t^{-2} obtained by Haff's minimal model [J. Fluid Mech. 134, 401 (1983)JFLSA70022-112010.1017/S0022112083003419] proves to be robust, despite the violation of several of its central assumptions. The shape anisotropy of the grains influences the characteristic time of energy loss quantitatively but not qualitatively. We compare kinetic energies in the individual degrees of freedom and find a slight predominance of translational motions. In addition, we observe a preferred rod alignment in the flight direction, as known from active matter or animal flocks.

Controlled Viscosity in Dense Granular Materials

We experimentally investigate the fluidization of a granular material subject to mechanical vibrations by monitoring the angular velocity of a vane suspended in the medium and driven by an external motor. On increasing the frequency, we observe a reentrant transition, as a jammed system first enters a fluidized state, where the vane rotates with high constant velocity, and then returns to a frictional state, where the vane velocity is much lower. While the fluidization frequency is material independent, the viscosity recovery frequency shows a clear dependence on the material that we rationalize by relating this frequency to the balance between dissipative and inertial forces in the system. Molecular dynamics simulations well reproduce the experimental data, confirming the suggested theoretical picture.

Heat flux of driven granular mixtures at low density: Stability analysis of the homogeneous steady state

The Navier-Stokes order hydrodynamic equations for a low-density driven granular mixture obtained previously [Khalil and Garzó, Phys. Rev. E 88, 052201 (2013)PLEEE81539-375510.1103/PhysRevE.88.052201] from the Chapman-Enskog solution to the Boltzmann equation are considered further. The four transport coefficients associated with the heat flux are obtained in terms of the mass ratio, the size ratio, composition, coefficients of restitution, and the driven parameters of the model. Their quantitative variation on the control parameters of the system is demonstrated by considering the leading terms in a Sonine polynomial expansion to solve the exact integral equations. As an application of these results, the stability of the homogeneous steady state is studied. In contrast to the results obtained in undriven granular mixtures, the stability analysis of the linearized Navier-Stokes hydrodynamic equations shows that the transversal and longitudinal modes are (linearly) stable with respect to long enough wavelength excitations. This conclusion agrees with a previous analysis made for single granular gases.

Magnetic resonance imaging of granular materials

Magnetic Resonance Imaging (MRI) has become one of the most important tools to screen humans in medicine; virtually every modern hospital is equipped with a Nuclear Magnetic Resonance (NMR) tomograph. The potential of NMR in 3D imaging tasks is by far greater, but there is only "a handful" of MRI studies of particulate matter. The method is expensive, time-consuming, and requires a deep understanding of pulse sequences, signal acquisition, and processing. We give a short introduction into the physical principles of this imaging technique, describe its advantages and limitations for the screening of granular matter, and present a number of examples of different application purposes, from the exploration of granular packing, via the detection of flow and particle diffusion, to real dynamic measurements. Probably, X-ray computed tomography is preferable in most applications, but fast imaging of single slices with modern MRI techniques is unmatched, and the additional opportunity to retrieve spatially resolved flow and diffusion profiles without particle tracking is a unique feature.

Pattern transition, microstructure, and dynamics in a two-dimensional vibrofluidized granular bed

Experiments are conducted in a two-dimensional monolayer vibrofluidized bed of glass beads, with a goal to understand the transition scenario and the underlying microstructure and dynamics in different patterned states. At small shaking accelerations (Γ=Aω^{2}/g<1, where A and ω=2πf are the amplitude and angular frequency of shaking and g is the gravitational acceleration), the particles remain attached to the base of the vibrating container; this is known as the solid bed (SB). With increasing Γ (at large enough shaking amplitude A/d) and/or with increasing A/d (at large enough Γ), the sequence of transitions/bifurcations unfolds as follows: SB ("solid bed") to BB ("bouncing bed") to LS ("Leidenfrost state") to "2-roll convection" to "1-roll convection" and finally to a gas-like state. For a given length of the container, the coarsening of multiple convection rolls leading to the genesis of a "single-roll" structure (dubbed the multiroll transition) and its subsequent transition to a granular gas are two findings of this work. We show that the critical shaking intensity (Γ_{BB}^{LS}) for the BB→LS transition has a power-law dependence on the particle loading (F=h_{0}/d, where h_{0} is the number of particle layers at rest and d is the particle diameter) and the shaking amplitude (A/d). The characteristics of BB and LS states are studied by calculating (i) the coarse-grained density and temperature profiles and (ii) the pair correlation function. It is shown that while the contact network of particles in the BB state represents a hexagonal-packed structure, the contact network within the "floating cluster" of the LS resembles a liquid-like state. An unsteadiness of the Leidenfrost state has been uncovered wherein the interface (between the floating cluster and the dilute collisional layer underneath) and the top of the bed are found to oscillate sinusoidally, with the oscillation frequency closely matching the frequency of external shaking. Therefore, the granular Leidenfrost state is a period-1 wave as is the case for the BB state.

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.

Transport coefficients for driven granular mixtures at low density

The transport coefficients of a granular binary mixture driven by a stochastic bath with friction are determined from the inelastic Boltzmann kinetic equation. A normal solution is obtained via the Chapman-Enskog method for states near homogeneous steady states. The mass, momentum, and heat fluxes are determined to first order in the spatial gradients of the hydrodynamic fields, and the associated transport coefficients are identified. They are given in terms of the solutions of a set of coupled linear integral equations. As in the monocomponent case, since the collisional cooling cannot be compensated for locally by the heat produced by the external driving, the reference distributions (zeroth-order approximations) f(i)((0)) (i=1,2) for each species depend on time through their dependence on the pressure and the temperature. Explicit forms for the diffusion transport coefficients and the shear viscosity coefficient are obtained by assuming the steady-state conditions and by considering the leading terms in a Sonine polynomial expansion. A comparison with previous results obtained for granular Brownian motion and by using a (local) stochastic thermostat is also carried out. The present work extends previous theoretical results derived for monocomponent dense gases [Garzó, Chamorro, and Vega Reyes, Phys. Rev. E 87, 032201 (2013)] to granular mixtures at low density.

Granular gases of rod-shaped grains in microgravity

Granular gases are convenient model systems to investigate the statistical physics of nonequilibrium systems. In the literature, one finds numerous theoretical predictions, but only few experiments. We study a weakly excited dilute gas of rods, confined in a cuboid container in microgravity during a suborbital rocket flight. With respect to a gas of spherical grains at comparable filling fraction, the mean free path is considerably reduced. This guarantees a dominance of grain-grain collisions over grain-wall collisions. No clustering was observed, unlike in similar experiments with spherical grains. Rod positions and orientations were determined and tracked. Translational and rotational velocity distributions are non-Gaussian. Equipartition of kinetic energy between translations and rotations is violated.

Dynamics of the wet granular Leidenfrost phenomenon

By event-driven molecular dynamics simulations, we study the Leidenfrost effect for wet granular matter driven from below. In marked contrast to all earlier studies on other fluids, the dense plug hovering on the hot gas cushion undergoes an undamped oscillation. The location of the Hopf bifurcation leading to this oscillation is strongly dependent on the inelasticity of the grain impacts. The vertical separation into a gas phase with a condensed plug hovering above it is particularly pronounced due to the cohesiveness of the granulate. For sufficiently large system sizes, the Rayleigh-Taylor instability terminates the oscillatory state at late times.

Transport coefficients of driven granular fluids at moderate volume fraction

In a recent publication [Phys. Rev. E 83, 011301 (2011)], Vollmayr-Lee et al. have determined by computer simulations the thermal diffusivity and the longitudinal viscosity coefficients of a driven granular fluid of hard spheres at intermediate volume fractions. Although they compare their simulation results with the predictions of kinetic theory, they use the dilute expressions for the driven system and the modified Sonine approximations for the undriven system. The goal here is to carry out this comparison by proposing a modified Sonine approximation to the Enskog equation for driven granular fluids that leads to a better quantitative agreement.

Particle kinematics in a dilute, three-dimensional, vibration-fluidized granular medium

We report an experimental study of particle kinematics in a three-dimensional system of inelastic spheres fluidized by intense vibration. The motion of particles in the interior of the medium is tracked by high-speed video imaging, yielding a spatially resolved measurement of the velocity distribution. The distribution is wider than a Gaussian and broadens continuously with increasing volume fraction. The deviations from a Gaussian distribution for this boundary-driven system are different in sign and larger in magnitude than predictions for homogeneously driven systems. We also find correlations between velocity components which grow with increasing volume fraction.

Mass transport of impurities in a moderately dense granular gas

Transport coefficients associated with the mass flux of impurities immersed in a moderately dense granular gas of hard disks or spheres described by the inelastic Enskog equation are obtained by means of the Chapman-Enskog expansion. The transport coefficients are determined as the solutions of a set of coupled linear integral equations recently derived for polydisperse granular mixtures [Garzó, Phys. Rev. E 76, 031304 (2007)]. With the objective of obtaining theoretical expressions for the transport coefficients that are sufficiently accurate for highly inelastic collisions, we solve the above integral equations by using the second Sonine approximation. As a complementary route, we numerically solve by means of the direct simulation Monte Carlo method (DSMC) the inelastic Enskog equation to get the kinetic diffusion coefficient D0 for two and three dimensions. We have observed in all our simulations that the disagreement, for arbitrarily large inelasticity, in the values of both solutions (DSMC and second Sonine approximation) is less than 4%. Moreover, we show that the second Sonine approximation to D0 yields a dramatic improvement (up to 50%) over the first Sonine approximation for impurity particles lighter than the surrounding gas and in the range of large inelasticity. The results reported in this paper are of direct application in important problems in granular flows, such as segregation driven by gravity and a thermal gradient. We analyze here the segregation criteria that result from our theoretical expressions of the transport coefficients.

Visualization of collisional substructure in granular shock waves

We study shock wave formation and propagation in an experimental vertically driven quasi-two-dimensional granular gas. We measure the moments of the single particle velocity distribution as a function of space and time. The space-time fields of the velocity moments show acoustic waves with a serrated substructure on the scale of a particle diameter. We show that this substructure is the result of collisional transport in which sequential collisions each transport momentum and energy by one particle diameter.

Linear response and hydrodynamics for granular fluids

A formal derivation of linear hydrodynamics for a granular fluid is given. The linear response to small spatial perturbations of a homogeneous reference state is studied in detail, using methods of nonequilibrium statistical mechanics. A transport matrix for macroscopic excitations in the fluid is defined in terms of the response functions. An expansion in the wave vector to second order allows identification of all phenomenological susceptibilities and transport coefficients through Navier-Stokes order in terms of appropriate time correlation functions. The transport coefficients in this representation are the generalization to granular fluids of the familiar Helfand and Green-Kubo relations for normal fluids. The analysis applies to a variety of collision rules. Important differences in both the analysis and results from those for normal fluids are identified and discussed. A scaling limit is described corresponding to the conditions under which idealized inelastic hard sphere models can apply. Further details and interpretation are provided in the paper following this one, by specialization to the case of smooth, inelastic hard spheres with constant coefficient of restitution.

Enskog theory for polydisperse granular mixtures. I. Navier-Stokes order transport

A hydrodynamic description for an s -component mixture of inelastic, smooth hard disks (two dimensions) or spheres (three dimensions) is derived based on the revised Enskog theory for the single-particle velocity distribution functions. In this first part of the two-part series, the macroscopic balance equations for mass, momentum, and energy are derived. Constitutive equations are calculated from exact expressions for the fluxes by a Chapman-Enskog expansion carried out to first order in spatial gradients, thereby resulting in a Navier-Stokes order theory. Within this context of small gradients, the theory is applicable to a wide range of restitution coefficients and densities. The resulting integral-differential equations for the zeroth- and first-order approximations of the distribution functions are given in exact form. An approximate solution to these equations is required for practical purposes in order to cast the constitutive quantities as algebraic functions of the macroscopic variables; this task is described in the companion paper.

Time-of-flight variant to image mixing of granular media in a 3D fluidized bed

This paper describes a variant of time-of-flight magnetic resonance (MR) imaging that provides a method of measuring the inherent mixing in a fluidized bed without the introduction of tracer particles. The modifications to conventional time-of-flight imaging enable the measurement of the axial mixing of a precisely controlled initial particle distribution, thereby providing measurements suitable for a direct comparison with models of solids mixing in granular systems. The imaging sequence is applied to characterize mixing, over time scales of 25-1000 ms, in a gas-fluidized bed of Myosotis seed particles; mixing over short timescales, inaccessible using conventional tracer techniques, is studied using this technique. The mixing pattern determined by this pulse sequence is used in conjunction with MR velocity images of the motion of the particles to provide new insight into the mechanism of solids mixing in granular systems.

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.

Velocity statistics in dissipative, dense granular media

We use a two-dimensional random-force model to investigate the velocity distributions in driven granular media. In general, the shape of the distribution is found to depend on the degree of dissipation and the packing fraction but, in highly dissipative systems, the velocity distributions have tails close to exponential. We show that these arise from the dynamics of single particles traveling in dilute regions and influenced predominantly by the random force. A self-consistent kinetic theory is developed to describe this behavior.

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.

Flux through a hole from a shaken granular medium

We have measured the flux of grains from a hole in the bottom of a shaken container of grains. We find that the peak velocity of the vibration, v max, controls the flux, i.e., the flux is nearly independent of the frequency and acceleration amplitude for a given value of v max. The flux decreases with increasing peak velocity and then becomes almost constant for the largest values of v max. The data at low peak velocity can be quantitatively described by a simple model, but the crossover to nearly constant flux at larger peak velocity suggests a regime in which the granular density near the container bottom is independent of the energy input to the system.

Mechanisms in the size segregation of a binary granular mixture

A granular mixture of particles of two sizes that is shaken vertically will in most cases segregate. If the larger particles accumulate at the top of the sample, this is called the Brazil-nut effect (BNE); if they accumulate at the bottom, it is called the reverse Brazil-nut effect (RBNE). While this process is of great industrial importance in the handling of bulk solids, it is not well understood. In recent years ten different mechanisms have been suggested to explain when each type of segregation is observed. However, the dependence of the mechanisms on driving conditions and material parameters and hence their relative importance is largely unknown. In this paper we present experiments and simulations where both types of particles are made from the same material and shaken under low air pressure, which reduces the number of mechanisms to be considered to seven. We observe both BNE and RBNE by varying systematically the driving frequency and amplitude, diameter ratio, ratio of total volume of small to large particles, and overall sample volume. All our results can be explained by a combination of three mechanisms: a geometrical mechanism called void filling, transport of particles in sidewall-driven convection rolls, and thermal diffusion, a mechanism predicted by kinetic theory.

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.

Shock wave propagation in vibrofluidized granular materials

Shock wave formation and propagation in vertically vibrated quasi-two-dimensional granular materials are studied by digital high speed photography. Steep density and temperature wave fronts form at the bottom of the granular layer when the layer collides with vibrating plate. Then the fronts propagate upwards through the layer. The temperature front is always in the transition region between the upward and downward granular flows. The effects of driving parameters and particle number on the shock are also explored.

Application of k- and q-space encoding NMR techniques on granular media in a 3D model fluidized bed reactor

A combination of PFG-NMR imaging and velocity encoding methods was applied to investigate the dynamic behavior of a bed of poppy seeds subjected to air flow, representing a model setup for fluidized bed reactors. The particle motion is described both from a statistical point of view, by determining propagators and dispersion coefficients representing an average over the whole bed volume, as well as combined with spatial resolution by generating velocity maps. Velocity images of different horizontal slices in the bed confirm the notion of a toroidal particle flow pattern inside the shallow granular bed. Despite the need of considerable averaging due to the random motion of the relatively few particles in the bed, quantitative velocity images and statistical information about the random particle motion can be obtained from monitoring the fluid component in the seeds by conventional spin-echo techniques.

Transport coefficients for an inelastic gas around uniform shear flow: linear stability analysis

The inelastic Boltzmann equation for a granular gas is applied to spatially inhomogeneous states close to uniform shear flow. A normal solution is obtained via a Chapman-Enskog-like expansion around a local shear flow distribution. The heat and momentum fluxes are determined to first order in the deviations of the hydrodynamic field gradients from their values in the reference state. The corresponding transport coefficients are determined from a set of coupled linear integral equations which are approximately solved by using a kinetic model of the Boltzmann equation. The main new ingredient in this expansion is that the reference state f(0) (zeroth-order approximation) retains all the hydrodynamic orders in the shear rate. In addition, since the collisional cooling cannot be compensated locally for viscous heating, the distribution f(0) depends on time through its dependence on temperature. This means that in general, for a given degree of inelasticity, the complete nonlinear dependence of the transport coefficients on the shear rate requires analysis of the unsteady hydrodynamic behavior. To simplify the analysis, the steady-state conditions have been considered here in order to perform a linear stability analysis of the hydrodynamic equations with respect to the uniform shear flow state. Conditions for instabilities at long wavelengths are identified and discussed.

Collision frequencies and energy flux in a dilute granular gas

Recent experimental study of a granular gas fluidized by vibrations in a low gravity environment has reported that the collision frequency nu(rho) of the particles with the container boundary scales roughly like N(alpha) with alpha=0.6 +/- 0.1, where N is the number of particles. Using numerical simulations, we show that this scaling is observed on a wide range of N, both for nu(rho) and for the particle-particle collision frequency nu(c). Simple scaling arguments show that this behavior is related to the energy flux in the granular gas, from injection by the moving boundary to dissipation by inelastic collisions. We predict in the dilute limit that the collision frequencies scale such as (sqrt)N are in fair agreement with experimental measurements.

Granular Leidenfrost effect: experiment and theory of floating particle clusters

Granular material is vertically vibrated in a 2D container: above a critical shaking strength, and for a sufficient number of beads, a crystalline cluster is elevated and supported by a dilute gaseous layer of fast beads underneath. We call this phenomenon the granular Leidenfrost effect. The experimental observations are explained by a hydrodynamic model featuring three dimensionless control parameters: the energy input S, the number of particle layers F, and the inelasticity of the particle collisions epsilon. The (S,F) phase diagram, in which the Leidenfrost state lies between the purely solid and gas phases, shows accurate agreement between experiment and theory.

Statistics of particle velocities in dense granular flows

We present measurements of the particle velocity distribution in the slow flow of granular material through vertical channels. The velocities of particles adjacent to the smooth, transparent front face of the channel were determined by video imaging and particle tracking. We find that the mean velocity changes sharply in shear layers near the side walls, but remains constant in a substantial core. The velocity distribution is non-Gaussian, is anisotropic, and follows a power law at large velocities. Remarkably, the distribution is identical in the shear layer and the core. We show evidence of spatially correlated motion, and propose a mechanism for the generation of fluctuational motion in the absence of shear.

Transport properties of dense dissipative hard-sphere fluids for arbitrary energy loss models

The revised Enskog approximation for a fluid of hard spheres which lose energy upon collision is discussed for the case that the energy is lost from the normal component of the velocity at collision but is otherwise arbitrary. Granular fluids with a velocity-dependent coefficient of restitution are an important special case covered by this model. A normal solution to the Enskog equation is developed using the Chapman-Enskog expansion. The lowest order solution describes the general homogeneous cooling state and a generating function formalism is introduced for the determination of the distribution function. The first order solution, evaluated in the lowest Sonine approximation, provides estimates for the transport coefficients for the Navier-Stokes hydrodynamic description. All calculations are performed in an arbitrary number of dimensions.

Instabilities in a free granular fluid described by the Enskog equation

A linear stability analysis of the hydrodynamic equations with respect to the homogeneous cooling state is carried out to identify the conditions for stability as functions of the wave vector, the dissipation, and the density. In contrast to previous studies, this description is based on the results derived from the Enskog equation for inelastic hard spheres [V. Garzó and J. W. Dufty, Phys. Rev. E 59, 5895 (1999)], which takes into account the dependence of the transport coefficients on dissipation. As expected, linear stability shows two transversal (shear) modes and a longitudinal "heat") mode to be unstable with respect to long enough wavelength excitations. Comparison with previous results (which neglect the influence of dissipation on transport) shows quantitative discrepancies for strong dissipation.

Crucial role of sidewalls in velocity distributions in quasi-two-dimensional granular gases

Our experiments and three-dimensional molecular dynamics simulations of particles confined to a vertical monolayer by closely spaced frictional walls (sidewalls) yield velocity distributions with non-Gaussian tails and a peak near zero velocity. Simulations with frictionless sidewalls are not peaked. Thus interactions between particles and their containers are an important determinant of the shape of the distribution and should be considered when evaluating experiments on a constrained monolayer of particles.