Velocity distributions and correlations in homogeneously heated granular media

Velocity distribution of a driven inelastic one-component Maxwell gas

The nature of the velocity distribution of a driven granular gas, though well studied, is unknown as to whether it is universal or not, and, if universal, what it is. We determine the tails of the steady state velocity distribution of a driven inelastic Maxwell gas, which is a simple model of a granular gas where the rate of collision between particles is independent of the separation as well as the relative velocity. We show that the steady state velocity distribution is nonuniversal and depends strongly on the nature of driving. The asymptotic behavior of the velocity distribution is shown to be identical to that of a noninteracting model where the collisions between particles are ignored. For diffusive driving, where collisions with the wall are modeled by an additive noise, the tails of the velocity distribution is universal only if the noise distribution decays faster than exponential.

Driven inelastic Maxwell gas in one dimension

A lattice version of the driven inelastic Maxwell gas is studied in one dimension with periodic boundary conditions. Each site i of the lattice is assigned with a scalar "velocity," v_{i}. Nearest neighbors on the lattice interact, with a rate τ_{c}^{-1}, according to an inelastic collision rule. External driving, occurring with a rate τ_{w}^{-1}, sustains a steady state in the system. A set of closed coupled equations for the evolution of the variance and the two-point correlation is found. Steady-state values of the variance, as well as spatial correlation functions, are calculated. It is shown exactly that the correlation function decays exponentially with distance, and the correlation length for a large system is determined. Furthermore, the spatiotemporal correlation C(x,t)=〈v_{i}(0)v_{i+x}(t)〉 can also be obtained. We find that there is an interior region -x^{*}<x<x^{*}, where C(x,t) has a time-dependent form, whereas in the exterior region |x|>x^{*}, the correlation function remains the same as the initial form. C(x,t) exhibits second-order discontinuity at the transition points x=±x^{*}, and these transition points move away from the x=0 with a constant speed.

Transition to a labyrinthine phase in a driven granular medium

Labyrinthine patterns arise in two-dimensional physical systems submitted to competing interactions, in fields ranging from solid-state physics to hydrodynamics. For systems of interacting particles, labyrinthine and stripe phases were studied in the context of colloidal particles confined into a monolayer, both numerically by means of Monte Carlo simulations and experimentally using superparamagnetic particles. Here we report an experimental observation of a labyrinthine phase in an out-of-equilibrium system constituted of macroscopic particles. Once sufficiently magnetized, they organize into short chains of particles in contact and randomly orientated. We characterize the transition from a granular gas state towards a solid labyrinthine phase, as a function of the ratio of the interaction strength to the kinetic agitation. The spatial local structure is analyzed by means of accurate particle tracking. Moreover, we explain the formation of these chains using a simple model.

Velocity statistics of dynamic spinners in out-of-equilibrium magnetic suspensions

We report on the velocity statistics of an out-of-equilibrium magnetic suspension in a spinner phase confined at a liquid interface. The suspension is energized by a uniaxial alternating magnetic field applied parallel to the interface. In a certain range of the magnetic field parameters the system spontaneously undergoes a transition into a dynamic spinner phase (ensemble of hydrodynamically coupled magnetic micro-rotors) comprised of two subsystems: self-assembled spinning chains and a gas of rotating single particles. Both subsystems coexist in a dynamic equilibrium via continuous exchange of the particles. Spinners excite surface flows that significantly increase particle velocity correlations in the system. For both subsystems the velocity distributions are strongly non-Maxwellian with nearly exponential high-energy tails, P(v) ∼ exp(-|v/v0|). The kurtosis, the measure of the deviation from the Gaussian statistics, is influenced by the frequency of the external magnetic field. We show that in the single-particle gas the dissipation is mostly collisional, whereas the viscous damping dominates over collisional dissipation for the self-assembled spinners. The dissipation increases with the frequency of the applied magnetic field. Our results provide insights into non-trivial dissipation mechanisms determining self-assembly processes in out-of-equilibrium magnetic suspensions.

Effect of Coulomb friction on orientational correlation and velocity distribution functions in a sheared dilute granular gas

From particle simulations of a sheared frictional granular gas, we show that the Coulomb friction can have dramatic effects on orientational correlation as well as on both the translational and angular velocity distribution functions even in the Boltzmann (dilute) limit. The dependence of orientational correlation on friction coefficient (μ) is found to be nonmonotonic, and the Coulomb friction plays a dual role of enhancing or diminishing the orientational correlation, depending on the value of the tangential restitution coefficient (which characterizes the roughness of particles). From the sticking limit (i.e., with no sliding contact) of rough particles, decreasing the Coulomb friction is found to reduce the density and spatial velocity correlations which, together with diminished orientational correlation for small enough μ, are responsible for the transition from non-gaussian to gaussian distribution functions in the double limit of small friction (μ→0) and nearly elastic particles (e→1). This double limit in fact corresponds to perfectly smooth particles, and hence the maxwellian (gaussian) is indeed a solution of the Boltzmann equation for a frictional granular gas in the limit of elastic collisions and zero Coulomb friction at any roughness. The high-velocity tails of both distribution functions seem to follow stretched exponentials even in the presence of Coulomb friction, and the related velocity exponents deviate strongly from a gaussian with increasing friction.

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.

Boltzmann equation for dissipative gases in homogeneous states with nonlinear friction

Combining analytical and numerical methods, we study within the framework of the homogeneous nonlinear Boltzmann equation a broad class of models relevant for the dynamics of dissipative fluids, including granular gases. We use the method presented in a previous paper [J. Stat. Phys. 124, 549 (2006)] and extend our results to a different heating mechanism: namely, a deterministic nonlinear friction force. We derive analytically the high-energy tail of the velocity distribution and compare the theoretical predictions with high-precision numerical simulations. Stretched exponential forms are obtained when the nonequilibrium steady state is stable. We derive subleading corrections and emphasize their relevance. In marginal stability cases, power-law behaviors arise, with exponents obtained as the roots of transcendental equations. We also consider some simple Bhatnagar-Gross-Krook models, driven by similar heating devices, to test the robustness of our predictions.

Spatial force correlations in granular shear flow. II. Theoretical implications

Numerical simulations are used to test the kinetic theory constitutive relations of inertial granular shear flow. These predictions are shown to be accurate in the dilute regime, where only binary collisions are relevant, but underestimate the measured value in the dense regime, where force networks of size xi are present. The discrepancy in the dense regime is due to non-collisional forces that we measure directly in our simulations and arise from elastic deformations of the force networks. We model the non-collisional stress by summing over all paths that elastic waves travel through force networks. This results in an analytical theory that successfully predicts the stress tensor over the entire inertial regime without any adjustable parameters.

Forcing independent velocity distributions in an experimental granular fluid

We present experimental results on the velocity statistics of a granular fluid with an effective stochastic thermostat, in a quasi-two-dimensional configuration. We find the base state, as measured by the single particle velocity distribution P(c) in the central high-probability regions, to be well described by P(c)=f{MB}[1+a2S2(c2)] : It deviates from a Maxwell-Boltzmann f{MB} by a second order Sonine polynomial S2(c2) with a single adjustable parameter a2. We find a2 to be a function of the filling fraction and independent of the driving over a wide range of frequencies and accelerations. Moreover, there is a consistent overpopulation in the distribution's tails, which scale as P proportional, variant exp(-A x c{3/2}) . To our knowledge, this is the first time that Sonine deviations have been measured in an experimental system.

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.

Dynamics of a tracer granular particle as a nonequilibrium Markov process

The dynamics of a tracer particle in a stationary driven granular gas is investigated. We show how to transform the linear Boltzmann equation, describing the dynamics of the tracer into a master equation for a continuous Markov process. The transition rates depend on the stationary velocity distribution of the gas. When the gas has a Gaussian velocity probability distribution function (PDF), the stationary velocity PDF of the tracer is Gaussian with a lower temperature and satisfies detailed balance for any value of the restitution coefficient alpha. As soon as the velocity PDF of the gas departs from the Gaussian form, detailed balance is violated. This nonequilibrium state can be characterized in terms of a Lebowitz-Spohn action functional W(tau) defined over trajectories of time duration tau. We discuss the properties of this functional and of a similar functional W(tau), which differs from the first for a term that is nonextensive in time. On the one hand, we show that in numerical experiments (i.e., at finite times tau), the two functionals have different fluctuations and W always satisfies an Evans-Searles-like symmetry. On the other hand, we cannot observe the verification of the Lebowitz-Spohn-Gallavotti-Cohen (LS-GC) relation, which is expected for W(tau) at very large times tau. We give an argument for the possible failure of the LS-GC relation in this situation. We also suggest practical recipes for measuring W(tau) and W(tau) in experiments.

Velocity distributions in dilute granular systems

We investigate the idea that velocity distributions in granular gases are determined mainly by eta, the coefficient of restitution and q, which measures the relative importance of heating (or energy input) to collisions. To this end, we study by numerical simulation the properties of inelastic gases as functions of eta, concentration phi, and particle number N with various heating mechanisms. For a wide range of parameters, we find Gaussian velocity distributions for uniform heating and non-Gaussian velocity distributions for boundary heating. Comparison between these results and velocity distributions obtained by other heating mechanisms and for a simple model of a granular gas without spatial degrees of freedom, shows that uniform and boundary heating can be understood as different limits of q, with q>>1 and q < or approximately 1 respectively. We review the literature for evidence of the role of q in the recent experiments.

Velocity distributions of granular gases with drag and with long-range interactions

We study velocity statistics of electrostatically driven granular gases. For two different experiments, (i) nonmagnetic particles in a viscous fluid and (ii) magnetic particles in air, the velocity distribution is non-Maxwellian, and its high-energy tail is exponential, P(upsilon) approximately exp(-/upsilon/). This behavior is consistent with the kinetic theory of driven dissipative particles. For particles immersed in a fluid, viscous damping is responsible for the exponential tail, while for magnetic particles, long-range interactions cause the exponential tail. We conclude that velocity statistics of dissipative gases are sensitive to the fluid environment and to the form of the particle interaction.

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.

Inelastic hard rods in a periodic potential

A simple model of inelastic hard rods subject to a one-dimensional array of identical wells is introduced. The energy loss due to inelastic collisions is balanced by the work supplied by an external stochastic heat bath. We explore the effect of the spatial nonuniformity on the steady states of the system. The spatial variations of the density, granular temperature, and pressure induced by the gradient of the external potential are investigated and compared with the analogous variations in an elastic system. Finally, we study the clustering process by considering the relaxation of the system starting from a uniform homogeneous state.

Velocity distributions in dissipative granular gases

Motivated by recent experiments reporting non-Gaussian velocity distributions in driven dilute granular materials, we study by numerical simulation the properties of 2D inelastic gases. We find theoretically that the form of the observed velocity distribution is governed primarily by the coefficient of restitution eta and q=N(H)/N(C), the ratio between the average number of heatings and the average number of collisions in the gas. The differences in distributions we find between uniform and boundary heating can then be understood as different limits of q, for q>>1 and q less, similar 1, respectively.

Steady-state velocity distributions of an oscillated granular gas

We use a three-dimensional molecular dynamics simulation to study the single particle distribution function of a dilute granular gas driven by a vertically oscillating plate at high accelerations (15g-90g). We find that the density and the temperature fields are essentially time-invariant above a height of about 40 particle diameters, where typically 20% of the grains are contained. These grains form the nonequilibrium steady-state granular gas with a Knudsen number unity or greater. In the steady-state region, the probability distribution function of the horizontal velocity c(x) (scaled by the local horizontal temperature) is found to be nearly independent of height, even though the hydrodynamic fields vary with height. We find that the high energy tails of the distribution function are described by a stretched exponential approximately exp(-Bcalphax), where alpha depends on the restitution coefficient e and falls in the range 1.2<alpha<1.6. However, alpha does not vary significantly for a wide range of friction coefficient values. We find that the distribution function of a frictionless inelastic hard sphere model can be made similar to that of a frictional model by adjusting e. However, there is no single value of e that mimics the frictional model over a range of heights.

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.

Shear viscosity for a heated granular binary mixture at low density

The shear viscosity for a heated granular binary mixture of smooth hard spheres at low density is analyzed. The mixture is heated by the action of an external driving force (Gaussian thermostat) that exactly compensates for cooling effects associated with the dissipation of collisions. The study is made from the Boltzmann kinetic theory, which is solved by using two complementary approaches. First, a normal solution of the Boltzmann equation via the Chapman-Enskog method is obtained up to first order in the spatial gradients. The mass, heat, and momentum fluxes are determined and the corresponding transport coefficients identified. As in the free cooling case [V. Garzó and J. W. Dufty, Phys. Fluids 14, 1476 (2002)], practical evaluation requires a Sonine polynomial approximation, and here it is mainly illustrated in the case of the shear viscosity. Second, to check the accuracy of the Chapman-Enskog results, the Boltzmann equation is numerically solved by means of the direct simulation Monte Carlo method. The simulation is performed for a system under uniform shear flow, using the Gaussian thermostat to control inelastic cooling. The comparison shows an excellent agreement between theory and simulation over a wide range of values of the restitution coefficients and the parameters of the mixture (masses, concentrations, and sizes).

Velocity fluctuations in electrostatically driven granular media

We study experimentally the particle velocity fluctuations in an electrostatically driven dilute granular gas. The velocity distributions have strong deviations from a Maxwellian form over a wide range of parameters. We have found that the tails of the distribution functions are consistent with a stretched exponential law with typical exponents of the order 3/2. Molecular dynamic simulations shows qualitative agreement with experimental data. Our results suggest that this non-Gaussian behavior is typical of most inelastic gases with both short- and long-range interactions.

Forcing and velocity correlations in a vibrated granular monolayer

The role of forcing on the dynamics of a vertically shaken granular monolayer is investigated. Using a flat plate, surprising negative velocity correlations are measured. A mechanism for this anticorrelation is proposed with support from both experimental results and molecular dynamics simulations. Using a rough plate, velocity correlations are positive, and the velocity distribution evolves from a Gaussian at very low densities to a broader distribution at high densities. These results are interpreted as a balance between stochastic forcing, interparticle collisions, and friction with the plate.

Steady-state properties of a mean-field model of driven inelastic mixtures

We investigate a Maxwell model of inelastic granular mixture under the influence of a stochastic driving and obtain its steady-state properties in the context of classical kinetic theory. The model is studied analytically by computing the moments up to the eighth order and approximating the distributions by means of a Sonine polynomial expansion method. The main findings concern the existence of two different granular temperatures, one for each species, and the characterization of the distribution functions, whose tails are in general more populated than those of an elastic system. These analytical results are tested against Monte Carlo numerical simulations of the model and are in general in good agreement. The simulations, however, reveal the presence of pronounced non-Gaussian tails in the case of an infinite temperature bath, which are not well reproduced by the Sonine method.

Randomly driven granular fluids: collisional statistics and short scale structure

We present a molecular-dynamics and kinetic theory study of granular material, modeled by inelastic hard disks, fluidized by a random driving force. The focus is on collisional averages and short-distance correlations in the nonequilibrium steady state, in order to analyze in a quantitative manner the breakdown of molecular chaos, i.e., factorization of the two-particle distribution function, f((2))(x(1),x(2)) approximately chif((1))(x(1))f((1))(x(2)) in a product of single-particle ones, where x(i)=[r(i),v(i)] with i=1,2 and chi represents the position correlation. We have found that molecular chaos is only violated in a small region of the two-particle phase space [x(1),x(2)], where there is a predominance of grazing collisions. The size of this singular region grows with increasing inelasticity. The existence of particle- and noise-induced recollisions magnifies the departure from mean-field behavior. The implications of this breakdown in several physical quantities are explored.