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+\alpha)/2$, 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\geq 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 collisionsComment: 20 pages; 6 figures; change of title; few minor changes; accepted for publication in PR

Self-diffusion in granular gases

The coefficient of self-diffusion for a homogeneously cooling granular gas changes significantly if the impact-velocity dependence of the restitution coefficient $\epsilon$ is taken into account. For the case of a constant $\epsilon$ the particles spread logarithmically slow with time, whereas the velocity dependent coefficient yields a power law time-dependence. The impact of the difference in these time dependences on the properties of a freely cooling granular gas is discussed.Comment: 6 pages, no figure

Uniform shear flow in dissipative gases. Computer simulations of inelastic hard spheres and (frictional) elastic hard spheres

In the preceding paper (cond-mat/0405252), we have conjectured that the main transport properties of a dilute gas of inelastic hard spheres (IHS) can be satisfactorily captured by an equivalent gas of elastic hard spheres (EHS), 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\leq \alpha\leq 0.95$ and two values of the imposed shear rate $a$. 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. The EHS results are in general hardly distinguishable from the IHS ones if $\alpha\gtrsim 0.7$, so that the distinct signature of the IHS gas (higher anisotropy and overpopulation) only manifests itself at relatively high dissipationsComment: 23 pages; 18 figures; Figs. 2 and 9 include new simulations; two new figures added; few minor changes; accepted for publication in PR

Critical Behavior of a Heavy Particle in a Granular Fluid

Behavior analogous to a second order phase transition is observed for the homogeneous cooling state of a heavy impurity particle in a granular fluid. The order parameter $\phi$ is the ratio of impurity mean square velocity to that of the fluid, with a conjugate field $h$ proportional to the mass ratio. A parameter $\beta$, measuring the fluid cooling rate relative to the impurity--fluid collision rate, is the analogue of the inverse temperature. For $\beta <1$ the fluid is ``normal'' with $\phi =0$ at $h=0$, as in the case of a system with elastic collisions. For $\beta >1$ an ``ordered'' state with $\phi \neq 0$ occurs at $h=0$, representing an extreme breakdown of equipartition. Critical slowing and qualitative changes in the velocity distribution function for the impurity particle near the transition are notedComment: 4 pages (4 figures included

Boundary effects in a quasi-two-dimensional driven granular fluid

The effect of a confining boundary on the spatial variations in granular temperature of a driven quasi-2d layer of particles is investigated experimentally. The radial drop in the relative granular temperature ΔT/T, exhibits a maximum at intermediate particle numbers which coincides with a crossover from kinetic to collisional transport of energy. It is also found that at low particle numbers, the distributions of radial velocities are increasingly asymmetric as one approaches the boundary. The radial and tangential granular temperatures split, and in the tails of the radial velocity distribution there is a higher population of fast moving particles travelling away rather than towards the boundary

Velocity Tails for Inelastic Maxwell Models

We study the velocity distribution function for inelastic Maxwell models, characterized by a Boltzmann equation with constant collision rate, independent of the energy of the colliding particles. By means of a nonlinear analysis of the Boltzmann equation, we find that the velocity distribution function decays algebraically for large velocities, with exponents that are analytically calculated.Comment: 4 pages, 2 figure

Diffusion in a Granular Fluid - Theory

Many important properties of granular fluids can be represented by a system of hard spheres with inelastic collisions. Traditional methods of nonequilibrium statistical mechanics are effective for analysis and description of the inelastic case as well. This is illustrated here for diffusion of an impurity particle in a fluid undergoing homogeneous cooling. An appropriate scaling of the Liouville equation is described such that the homogeneous cooling ensemble and associated time correlation functions map to those of a stationary state. In this form the familiar methods of linear response can be applied, leading to Green - Kubo and Einstein representations of diffusion in terms of the velocity and mean square displacement correlation functions. These correlation functions are evaluated approximately using a cumulant expansion and from kinetic theory, providing the diffusion coefficient as a function of the density and the restitution coefficients. Comparisons with results from molecular dynamics simulation are given in the following companion paper

Hydrodynamics of driven granular gases

Hydrodynamic equations for granular gases driven by the Fokker-Planck operator are derived. Transport coefficients appeared in Navier-Stokes order change from the values of a free cooling state to those of a steady state.Comment: 5 pages, 3 figure

Transport Coefficients for Granular Media from Molecular Dynamics Simulations

Under many conditions, macroscopic grains flow like a fluid; kinetic theory pred icts continuum equations of motion for this granular fluid. In order to test the theory, we perform event driven molecular simulations of a two-dimensional gas of inelastic hard disks, driven by contact with a heat bath. Even for strong dissipation, high densities, and small numbers of particles, we find that continuum theory describes the system well. With a bath that heats the gas homogeneously, strong velocity correlations produce a slightly smaller energy loss due to inelastic collisions than that predicted by kinetic theory. With an inhomogeneous heat bath, thermal or velocity gradients are induced. Determination of the resulting fluxes allows calculation of the thermal conductivity and shear viscosity, which are compared to the predictions of granular kinetic theory, and which can be used in continuum modeling of granular flows. The shear viscosity is close to the prediction of kinetic theory, while the thermal conductivity can be overestimated by a factor of 2; in each case, transport is lowered with increasing inelasticity.Comment: 14 pages, 17 figures, 39 references, submitted to PRE feb 199

Gaussian Kinetic Model for Granular Gases

A kinetic model for the Boltzmann equation is proposed and explored as a practical means to investigate the properties of a dilute granular gas. It is shown that all spatially homogeneous initial distributions approach a universal "homogeneous cooling solution" after a few collisions. The homogeneous cooling solution (HCS) is studied in some detail and the exact solution is compared with known results for the hard sphere Boltzmann equation. It is shown that all qualitative features of the HCS, including the nature of over population at large velocities, are reproduced semi-quantitatively by the kinetic model. It is also shown that all the transport coefficients are in excellent agreement with those from the Boltzmann equation. Also, the model is specialized to one having a velocity independent collision frequency and the resulting HCS and transport coefficients are compared to known results for the Maxwell Model. The potential of the model for the study of more complex spatially inhomogeneous states is discussed.Comment: to be submitted to Phys. Rev.