Steady-state velocity distributions of an oscillated granular gas

Tuning the distance to equipartition by controlling the collision rate in a driven granular gas experiment

In a granular gas experiment of magnetized particles confined in a thin layer, the rate of dissipative collisions is tuned by adjusting the amplitude of an external magnetic field. The velocity statistics are analyzed using the dynamic and static structure factors of transverse velocity modes. Using the fluctuating hydrodynamics theory, we measure the deviation from kinetic energy equipartition in this out-of-equilibrium system as a function of the dissipative collision rate. When the collision rate is decreased, the distance to equipartition becomes smaller, meaning that the dynamical properties of this granular gas approach by analogy those of a molecular gas in thermal equilibrium.

Cluster-Induced Deagglomeration in Dilute Gravity-Driven Gas-Solid Flows of Cohesive Grains

Clustering is often presumed to lead to enhanced agglomeration between cohesive grains due to the reduced relative velocities of particles within a cluster. Our discrete-particle simulations on gravity-driven, gas-solid flows of cohesive grains exhibit the opposite trend, revealing a new mechanism we coin "cluster-induced deagglomeration." Specifically, we examine relatively dilute gas-solid flows and isolate agglomerates of cohesive origin from overall heterogeneities in the system, i.e., agglomerates of cohesive origin and clusters of hydrodynamic origin. We observe enhanced clustering with an increasing system size (as is the norm for noncohesive systems) as well as reduced agglomeration. The reduced agglomeration is traced to the increased collisional impact velocities of particles at the surface of a cluster; i.e., higher levels of clustering lead to larger relative velocities between the clustered and nonclustered regions, thereby serving as an additional source of granular temperature. This physical picture is further evidenced by a theoretical model based on a balance between the generation and breakage rates of agglomerates. Finally, cluster-induced deagglomeration also provides an explanation for a surprising saturation of agglomeration levels in gravity-driven, gas-solid systems with increasing levels of cohesion, as opposed to the monotonically increasing behavior seen in free-evolving or driven granular systems in the absence of gravity. Namely, higher cohesion leads to more energy dissipation, which is associated with competing effects: enhanced agglomeration and enhanced clustering, the latter of which results in more cluster-induced deagglomeration.

Velocity Distribution of a Homogeneously Driven Two-Dimensional Granular Gas

The theory of homogeneously driven granular gases of hard particles predicts that the stationary state is characterized by a velocity distribution function with overpopulated high-energy tails as compared to the exponential decay valid for molecular gases. While this fundamental theoretical result was confirmed by numerous numerical simulations, an experimental confirmation is still missing. Using self-rotating active granular particles, we find a power-law decay of the velocity distribution whose exponent agrees well with the theoretic prediction.

Entropic force on granular chains self-extracting from one-dimensional confinement

The entropic forces on the self-retracting granular chains, which are confined in channels with different widths, are determined. The time dependence of the length of chain remaining in the channel Lin(t) is measured. The entropic force is treated as the only parameter in fitting the solution of the nonlinear equation of motion of Lin(t) to the experimental data. The dependence of the entropic force on the width of the confining channel can be expressed as a power-law with an exponent of 1.3, which is consistent with the previous theoretical predictions for the entropy loss due to confinement.

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.

Experimental measurements of the collapse of a two-dimensional granular gas under gravity

We experimentally measure the decay of a quasi-two-dimensional granular gas under gravity. A granular gas is created by vibrofluidization, after which the energy input is halted, and the time-dependent statistical properties of the decaying gas are measured with video particle tracking. There are two distinct cooling stages separated by a high temperature settling shock. In the final stage, the temperature of a fluid packet decreases as a power law T proportional, variant(t{c}-t);{alpha} just before the system collapses to a static state. The measured value of alpha ranges from 3.3 to 6.1 depending on the height, significantly higher than the exponent of 2 found in theoretical work on this problem [D. Volfson, B. Meerson, and L. S. Tsimring, Phys Rev. E 73, 61305 (2006)]. We also address the question of whether the collapse occurs simultaneously at different heights in the system.

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.

Orientational correlation and velocity distributions in uniform shear flow of a dilute granular gas

Using particle simulations of the uniform shear flow of a rough dilute granular gas, we show that the translational and rotational velocities are strongly correlated in direction, but there is no orientational correlation-induced singularity at perfectly smooth (beta=-1) and rough (beta=1) limits for elastic collisions (e=1); both the translational and rotational velocity distribution functions remain close to a Gaussian for these two limiting cases. Away from these two limits, the orientational as well as spatial velocity correlations are responsible for the emergence of non-Gaussian high-velocity tails. The tails of both distribution functions follow stretched exponentials, with the exponents depending on normal (e) and tangential (beta) restitution coefficients.

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.

Velocity distribution and the effect of wall roughness in granular Poiseuille flow

From event-driven simulations of a gravity-driven channel flow of inelastic hard disks, we show that the velocity distribution function remains close to a Gaussian for a wide range densities (even when the Knudsen number is of order 1) if the walls are smooth and the particle collisions are nearly elastic. For dense flows, a transition from a Gaussian to a power-law distribution for the high-velocity tails occurs with increasing dissipation in the center of the channel, irrespective of wall roughness. For a rough wall, the near-wall distribution functions are distinctly different from those in the bulk, even in the quasielastic limit.

Energy distribution and effective temperatures in a driven dissipative model

We investigate nonequilibrium behavior of driven dissipative systems, using a model we recently presented [Phys. Rev. Lett., 93, 240601 (2004)]. We solve the non-Boltzmann steady state energy distribution and the temporal evolution to it, and find its high energy tail to behave exponentially. We demonstrate that various measures of effective temperatures generally differ. We discuss infinite hierarchies of effective temperatures defined from moments of the nonexponential energy distribution, and relate them to the "configurational temperature," measured directly from instantaneous particle locations without any kinetic information. We calculate the "granular temperature," characterizing the average energy in the system, two different "fluctuation temperatures," scaling fluctuation-dissipation relations, and the "entropic temperature," defined from differentiating the entropy with respect to energy.

Impact of high-energy tails on granular gas properties

The velocity distribution function of granular gases in the homogeneous cooling state as well as some heated granular gases decays for large velocities as f proportional to exp(-const x v). That is, its high-energy tail is overpopulated as compared with the Maxwell distribution. At the present time, there is no theory to describe the influence of the tail on the kinetic characteristics of granular gases. We develop an approach to quantify the overpopulated tail and analyze its impact on granular gas properties, in particular on the cooling coefficient. We observe and explain anomalously slow relaxation of the velocity distribution function to its steady state.

Free cooling of the one-dimensional wet granular gas

The free cooling behavior of a wet granular gas is studied in one dimension. We employ a particularly simple model system in which the interaction of wet grains is characterized by a fixed energy loss assigned to each collision. Macroscopic laws of energy dissipation and cluster formation are studied on the basis of numerical simulations and mean-field analytical calculations. We find a number of remarkable scaling properties which may shed light on earlier unexplained results for related systems.

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.

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.

Equilibrium statistics of an inelastically bouncing ball, subject to gravity and a random force

We consider a particle moving on the half line x > 0 and subject to a constant force in the -x direction plus a delta-correlated random force. At x = 0 the particle is reflected inelastically. The velocities just after and before the reflection satisfy v(f) = -r v(i), where r is the coefficient of restitution. This simple model is of interest in connection with studies of driven granular matter in a gravitational field. With an exact analytical approach and simulations we study the steady-state distribution function P(x,v).

System of elastic hard spheres which mimics the transport properties of a granular gas

The prototype model of a fluidized granular system is a gas of inelastic hard spheres (IHS) with a constant coefficient of normal restitution alpha. Using a kinetic theory description we investigate the two basic ingredients that a model of elastic hard spheres (EHS) must have in order to mimic the most relevant transport properties of the underlying IHS gas. First, the EHS gas is assumed to be subject to the action of an effective drag force with a friction constant equal to half the cooling rate of the IHS gas, the latter being evaluated in the local equilibrium approximation for simplicity. Second, the collision rate of the EHS gas is reduced by a factor (1/2)(1+alpha), relative to that of the IHS gas. Comparison between the respective Navier-Stokes transport coefficients shows that the EHS model reproduces almost perfectly the self-diffusion coefficient and reasonably well the two transport coefficients defining the heat flux, the shear viscosity being reproduced within a deviation less than 14% (for alpha > or = 0.5). Moreover, the EHS model is seen to agree with the fundamental collision integrals of inelastic mixtures and dense gases. The approximate equivalence between IHS and EHS is used to propose kinetic models for inelastic collisions as simple extensions of known kinetic models for elastic collisions.

Uniform shear flow in dissipative gases: computer simulations of inelastic hard spheres and frictional elastic hard spheres

In the preceding paper, we have conjectured that the main transport properties of a dilute gas of inelastic hard spheres (IHSs) can be satisfactorily captured by an equivalent gas of elastic hard spheres (EHSs), provided that the latter are under the action of an effective drag force and their collision rate is reduced by a factor (1+alpha)/2 (where alpha is the constant coefficient of normal restitution). In this paper we test the above expectation in a paradigmatic nonequilibrium state, namely, the simple or uniform shear flow, by performing Monte Carlo computer simulations of the Boltzmann equation for both classes of dissipative gases with a dissipation range 0.5 < or = alpha < or = 0.95 and two values of the imposed shear rate a. It is observed that the evolution toward the steady state proceeds in two stages: a short kinetic stage (strongly dependent on the initial preparation of the system) followed by a slower hydrodynamic regime that becomes increasingly less dependent on the initial state. Once conveniently scaled, the intrinsic quantities in the hydrodynamic regime depend on time, at a given value of alpha, only through the reduced shear rate a*(t) is proportional to a/square root(T(t)), until a steady state, independent of the imposed shear rate and of the initial preparation, is reached. The distortion of the steady-state velocity distribution from the local equilibrium state is measured by the shear stress, the normal stress differences, the cooling rate, the fourth and sixth cumulants, and the shape of the distribution itself. In particular, the simulation results seem to be consistent with an exponential overpopulation of the high-velocity tail. These properties are common to both the IHS and EHS systems. In addition, the EHS results are in general hardly distinguishable from the IHS ones if alpha approximately > 0.7, so that the distinct signature of the IHS gas (higher anisotropy and overpopulation) only manifests itself at relatively high dissipations.

Simulations of vibrated granular medium with impact-velocity-dependent restitution coefficient

We report numerical simulations of strongly vibrated granular materials designed to mimic recent experiments performed in both the presence and the absence of gravity. The coefficient of restitution used here depends on the impact velocity by taking into account both the viscoelastic and plastic deformations of particles, occurring at low and high velocities, respectively. We show that this model with impact-velocity-dependent restitution coefficient reproduces results that agree with experiments. We measure the scaling exponents of the granular temperature, collision frequency, impulse, and pressure with the vibrating piston velocity as the particle number increases. As the system changes from a homogeneous gas state at low density to a clustered state at high density, these exponents are all found to decrease continuously with increasing particle number. All these results differ significantly from classical inelastic hard sphere kinetic theory and previous simulations, both based on a constant restitution coefficient.

Velocity distribution in a viscous granular gas

We investigate the velocity relaxation of a viscous one-dimensional granular gas in which neither energy nor momentum is conserved in a collision. Of interest is the distribution of velocities in the gas as it cools, and the time dependence of the relaxation behavior. A Boltzmann equation of instantaneous binary collisions leads to a two-peaked distribution, as do numerical simulations of grains on a line. Of particular note is that in the presence of friction there is no inelastic collapse, so there is no need to invoke additional assumptions such as the quasielastic limit.

Local equation of state and velocity distributions of a driven granular gas

We present event-driven simulations of a granular gas of inelastic hard disks with incomplete normal restitution in two dimensions between vibrating walls (without gravity). We measure hydrodynamic quantities such as the stress tensor, density and temperature profiles, as well as velocity distributions. Relating the local pressure to the local temperature and local density, we construct a local constitutive equation. For strong inelasticities the local constitutive relation depends on global system parameters, like the volume fraction and the aspect ratio. For moderate inelasticities the constitutive relation is approximately independent of the system parameters and can hence be regarded as a local equation of state, even though the system is highly inhomogeneous with heterogeneous temperature and density profiles arising as a consequence of energy injection. With respect to local velocity distributions we find that they do not scale with the square root of the local granular temperature. Moreover the high-velocity tails are different for the distribution of the x and the y components of the velocity, and even depend on the position in the sample, the global volume fraction, and the coefficient of restitution.

Role of friction in pattern formation in oscillated granular layers

Particles in granular flows are often modeled as frictionless (smooth) inelastic spheres; however, there exist no frictionless grains, just as there are no elastic grains. Our molecular dynamics simulations reveal that friction is essential for realistic modeling of vertically oscillated granular layers: simulations of frictionless particles yield patterns with an onset at a container acceleration about 30% smaller than that observed in experiments and simulations with friction. More importantly, even though square and hexagonal patterns form for a wide range of the oscillation parameters in experiments and in our simulations of frictional inelastic particles, only stripe patterns form in the simulations without friction, even if the inelasticity is increased to obtain as much dissipation as in frictional particles. We also consider the effect of particle friction on the shock wave that forms each time the granular layer strikes the container. While a shock wave still forms for frictionless particles, the spatial and temporal dependence of the hydrodynamic fields differ for the cases with and without friction.