Nonequilibrium phase transition for a heavy particle in a granular fluid

Kinetic Theory of Polydisperse Granular Mixtures: Influence of the Partial Temperatures on Transport Properties-A Review

It is well-recognized that granular media under rapid flow conditions can be modeled as a gas of hard spheres with inelastic collisions. At moderate densities, a fundamental basis for the determination of the granular hydrodynamics is provided by the Enskog kinetic equation conveniently adapted to account for inelastic collisions. A surprising result (compared to its molecular gas counterpart) for granular mixtures is the failure of the energy equipartition, even in homogeneous states. This means that the partial temperatures Ti (measuring the mean kinetic energy of each species) are different to the (total) granular temperature T. The goal of this paper is to provide an overview on the effect of different partial temperatures on the transport properties of the mixture. Our analysis addresses first the impact of energy nonequipartition on transport which is only due to the inelastic character of collisions. This effect (which is absent for elastic collisions) is shown to be significant in important problems in granular mixtures such as thermal diffusion segregation. Then, an independent source of energy nonequipartition due to the existence of a divergence of the flow velocity is studied. This effect (which was already analyzed in several pioneering works on dense hard-sphere molecular mixtures) affects to the bulk viscosity coefficient. Analytical (approximate) results are compared against Monte Carlo and molecular dynamics simulations, showing the reliability of kinetic theory for describing granular flows.

Intruder-induced asymmetry in compartmentalized granular gases

The clustering behavior of a compartmentalized monodisperse granular gas with the addition of one heavy intruding particle is investigated experimentally. Depending on the number of particles, the presence of a heavy intruder leads to three population states: a homogeneous state, an expelled clustering state, and a fully clustering state. These states are found to be consistent with the clustering of a purely monodisperse granular gas in an asymmetric compartmentalized structure. We obtain an exact relation between the size of an intruder and the elevation of the compartment bottom. This relation quantifies the particle-expelling ability of a heavy intruder, and suggests that the one-intruder system is a type of asymmetric system with an intruder-size-related asymmetrical index ξ. Under the framework of the flux model, a ξ-associated Ξ function is proposed to quantitatively reproduce the experimental results.

Intermediate regimes in granular Brownian motion: superdiffusion and subdiffusion

Brownian motion in a granular gas in a homogeneous cooling state is studied theoretically and by means of molecular dynamics. We use the simplest first-principles model for the impact-velocity dependent restitution coefficient, as it follows for the model of viscoelastic spheres. We reveal that for a wide range of initial conditions the ratio of granular temperatures of Brownian and bath particles demonstrates complicated nonmonotonic behavior, which results in a transition between different regimes of Brownian dynamics: It starts from the ballistic motion, switches later to a superballistic one, and turns at still later times into subdiffusion; eventually normal diffusion is achieved. Our theory agrees very well with the molecular dynamics results, although extreme computational costs prevented us from detecting the final diffusion regime. Qualitatively, the reported intermediate diffusion regimes are generic for granular gases with any realistic dependence of the restitution coefficient on the impact velocity.

Impurity in a sheared inelastic Maxwell gas

The Boltzmann equation for inelastic Maxwell models is considered in order to investigate the dynamics of an impurity (or intruder) immersed in a granular gas driven by a uniform shear flow. The analysis is based on an exact solution of the Boltzmann equation for a granular binary mixture. It applies for conditions arbitrarily far from equilibrium (arbitrary values of the shear rate a) and for arbitrary values of the parameters of the mixture (particle masses m(i), mole fractions x(i), and coefficients of restitution α(ij)). In the tracer limit where the mole fraction of the intruder species vanishes, a nonequilibrium phase transition takes place. We thereby identify ordered phases where the intruder bears a finite contribution to the properties of the mixture, in a region of parameter space that is worked out in detail. These findings extend previous results obtained for ordinary Maxwell gases, and further show that dissipation leads to new ordered phases.

Intruder-induced change in condensation temperature of granular gases

The process from a gaseous state to a clustering state for a compartmentalized monodisperse granular gas is accompanied by a drop in the granular temperature to a condensation point. We show experimentally that adding an intruder generally results in a decrease in the condensation point, and a heavier intruder makes this decrease more pronounced. However, once the Brazil nut effect (the intruder on the top of clustering grains) occurs, the condensation point will rise. Through the balance of particle fluxes and the hydrodynamic balance of driving forces, we analytically calculated the condensation point for the monodisperse gases and the intruder-fluid mixtures. The analytical results match the experimental data.

Dynamics of a hard sphere granular impurity

An impurity particle coupling to its host fluid via inelastic hard sphere collisions is considered. It is shown that the exact equation for its distribution function can be mapped onto that for an impurity with elastic collisions and an effective mass. The application of this result to the Enskog-Lorentz kinetic equation leads to several conclusions: (1) every solution in the elastic case is equivalent to a class of solutions in the granular case; (2) for an equilibrium host fluid the granular impurity approaches equilibrium at a different temperature, with a dominant diffusive mode at long times; (3) for a granular host fluid in its scaling state, the granular impurity approaches the corresponding scaling solution.

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.

Self-diffusion in granular gases: Green-Kubo versus Chapman-Enskog

We study the diffusion of tracers (self-diffusion) in a homogeneously cooling gas of dissipative particles, using the Green-Kubo relation and the Chapman-Enskog approach. The dissipative particle collisions are described by the coefficient of restitution epsilon which for realistic material properties depends on the impact velocity. First, we consider self-diffusion using a constant coefficient of restitution, epsilon=const, as frequently used to simplify the analysis. Second, self-diffusion is studied for a simplified (stepwise) dependence of epsilon on the impact velocity. Finally, diffusion is considered for gases of realistic viscoelastic particles. We find that for epsilon=const both methods lead to the same result for the self-diffusion coefficient. For the case of impact-velocity dependent coefficients of restitution, the Green-Kubo method is, however, either restrictive or too complicated for practical application, therefore we compute the diffusion coefficient using the Chapman-Enskog method. We conclude that in application to granular gases, the Chapman-Enskog approach is preferable for deriving kinetic coefficients.

Boltzmann equation for a granular capped rectangle in a thermalized bath of hard disks

By using the Boltzmann approach, we study the steady-state dynamics of a granular capped rectangle placed in a two-dimensional bath of thermalized hard disks. Hard core collisions are assumed elastic between disks and inelastic between the capped rectangle and the disks, with a normal coefficient of restitution alpha < 1. Assuming a Gaussian ansatz for the probability distribution functions, we obtain analytical expressions for the granular temperatures. We show the absence of equipartition and investigate both the role of the anisotropy of the capped rectangle and of the relative ratio of the bath particles to the linear sizes of the capped rectangle. In addition, we investigate a model of a capped rectangle with two normal coefficients of restitution for collisions along the straight and curved surfaces of the capped rectangle. In this case one observes equipartition for a nontrivial ratio of the normal coefficient of restitutions.

Diffusion of impurities in a granular gas

Diffusion of impurities in a granular gas undergoing homogeneous cooling state is studied. The results are obtained by solving the Boltzmann-Lorentz equation by means of the Chapman-Enskog method. In the first order in the density gradient of impurities, the diffusion coefficient D is determined as the solution of a linear integral equation which is approximately solved by making an expansion in Sonine polynomials. In this paper, we evaluate D up to the second order in the Sonine expansion and get explicit expressions for D in terms of the coefficients of restitution for the impurity-gas and gas-gas collisions as well as the ratios of mass and particle sizes. To check the reliability of the Sonine polynomial solution, analytical results are compared with those obtained from numerical solutions of the Boltzmann equation by means of the direct simulation Monte Carlo method. In the simulations, the diffusion coefficient is measured via the mean-square displacement of impurities. The comparison between theory and simulation shows in general an excellent agreement, except for the cases in which the gas particles are much heavier and/or much larger than impurities. In these cases, the second Sonine approximation to D improves significantly the qualitative predictions made from the first Sonine approximation. A discussion on the convergence of the Sonine polynomial expansion is also carried out.

Granular fluid thermostated by a bath of elastic hard spheres

The homogeneous steady state of a fluid of inelastic hard spheres immersed in a bath of elastic hard spheres kept at equilibrium is analyzed by means of the first Sonine approximation to the (spatially homogeneous) Enskog-Boltzmann equation. The temperature of the granular fluid relative to the bath temperature and the kurtosis of the granular distribution function are obtained as functions of the coefficient of restitution, the mass ratio, and a dimensionless parameter beta measuring the cooling rate relative to the friction constant. Comparison with recent results obtained from an iterative numerical solution of the Enskog-Boltzmann equation [Biben et al., Physica A 310, 308 (2002)] shows an excellent agreement. Several limiting cases are also considered. In particular, when the granular particles are much heavier than the bath particles (but have a comparable size and number density), it is shown that the bath acts as a white noise external driving. In the general case, the Sonine approximation predicts the lack of a steady state if the control parameter beta is larger than a certain critical value beta(c) that depends on the coefficient of restitution and the mass ratio. However, this phenomenon appears outside the expected domain of applicability of the approximation.

Kinetic temperatures for a granular mixture

An isolated mixture of smooth, inelastic hard spheres supports a homogeneous cooling state with different kinetic temperatures for each species. This phenomenon is explored here by molecular dynamics simulation of a two component fluid, with comparison to predictions of the Enskog kinetic theory. The ratio of kinetic temperatures is studied for two values of the restitution coefficient alpha=0.95 and 0.80, as a function of mass ratio, size ratio, composition, and density. Good agreement between theory and simulation is found for the lower densities and higher restitution coefficient; significant disagreement is observed otherwise. The phenomenon of different temperatures is also discussed for driven systems, as occurs in recent experiments. Differences between the freely cooling state and driven steady states are illustrated.

Tracer diffusion in granular shear flows

Tracer diffusion in a granular gas in simple shear flow is analyzed. The analysis is made from a perturbation solution of the Boltzmann kinetic equation through first order in the gradient of the mole fraction of tracer particles. The reference state (zeroth-order approximation) corresponds to a Sonine solution of the Boltzmann equation, which holds for arbitrary values of the restitution coefficients. Due to the anisotropy induced in the system by the shear flow, the mass flux defines a diffusion tensor D(ij) instead of a scalar diffusion coefficient. The elements of this tensor are given in terms of the restitution coefficients and mass and size ratios. The dependence of the diffusion tensor on the parameters of the problem is illustrated in the three-dimensional case. The results show that the influence of dissipation on the elements D(ij) is in general quite important, even for moderate values of the restitution coefficients. In the case of self-diffusion (mechanically equivalent particles), the trends observed in recent molecular-dynamics simulations are similar to those obtained here from the Boltzmann kinetic theory.

Diffusion in a granular fluid. I. 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.