Linear response and hydrodynamics for granular fluids

Stability analysis of the homogeneous hydrodynamics of a model for a confined granular gas

The linear hydrodynamic stability of a model for confined quasi-two-dimensional granular gases is analyzed. The system exhibits homogeneous hydrodynamics, i.e., there are macroscopic evolution equations for homogeneous states. The stability analysis is carried out around all these states and not only the homogeneous steady state reached eventually by the system. It is shown that in some cases the linear analysis is not enough to reach a definite conclusion on the stability, and molecular dynamics simulation results are presented to elucidate these cases. The analysis shows the relevance of nonlinear hydrodynamic contributions to describe the behavior of spontaneous fluctuations occurring in the system, that lead even to the transitory formation of clusters of particles. The conclusion is that the system is always stable. The relevance of the results for describing the instabilities of confined granular gases observed experimentally is discussed.

Divergence of Viscosity in Jammed Granular Materials: A Theoretical Approach

A theory for jammed granular materials is developed with the aid of a nonequilibrium steady-state distribution function. The approximate nonequilibrium steady-state distribution function is explicitly given in the weak dissipation regime by means of the relaxation time. The theory quantitatively agrees with the results of the molecular dynamics simulation on the critical behavior of the viscosity below the jamming point without introducing any fitting parameter.

Hydrodynamics for a model of a confined quasi-two-dimensional granular gas

The hydrodynamic equations for a model of a confined quasi-two-dimensional gas of smooth inelastic hard spheres are derived from the Boltzmann equation for the model, using a generalization of the Chapman-Enskog method. The heat and momentum fluxes are calculated to Navier-Stokes order, and the associated transport coefficients are explicitly determined as functions of the coefficient of normal restitution and the velocity parameter involved in the definition of the model. Also an Euler transport term contributing to the energy transport equation is considered. This term arises from the gradient expansion of the rate of change of the temperature due to the inelasticity of collisions, and it vanishes for elastic systems. The hydrodynamic equations are particularized for the relevant case of a system in the homogeneous steady state. The relationship with previous works is analyzed.

Power-law decay of the velocity autocorrelation function of a granular fluid in the homogeneous cooling state

The hydrodynamic part of the velocity autocorrelation function of a granular fluid in the homogeneous cooling state has been calculated by using mode-coupling theory for a finite system with periodic boundary conditions. The existence of the shearing instability, leading to a divergent behavior of the velocity flow fluctuations, is taken into account. A time region in which the velocity autocorrelation function exhibits a power-law decay, when time is measured by the number of collisions per particle, has been been identified. Also the explicit form of the exponential asymptotic long time decay has been obtained. The theoretical prediction for the power-law decay is compared with molecular dynamics simulation results, and a good agreement is found, after taking into account finite size corrections. The effects of approaching the shearing instability are also explored.

Homogeneous hydrodynamics of a collisional model of confined granular gases

The hydrodynamic equation governing the homogeneous time evolution of the temperature in a model of confined granular gas is studied by means of the Enskog equation. The existence of a normal solution of the kinetic equation is assumed as a condition for hydrodynamics. Dimensional analysis implies a scaling of the distribution function that is used to determine it in the first Sonine approximation, with a coefficient that evolves in time through its dependence on the temperature. The theoretical predictions are compared to numerical results obtained by the direct simulation Monte Carlo method and a good agreement is found. The relevance of the normal homogeneous distribution function to derive inhomogeneous hydrodynamic equations, for instance using the Champan-Enskog algorithm, is indicated.

Homogeneous steady state of a confined granular gas

The nonequilibrium statistical mechanics and kinetic theory for a model of a confined quasi-two-dimensional gas of inelastic hard spheres is presented. The dynamics of the particles includes an effective mechanism to transfer the energy injected in the vertical direction to the horizontal degrees of freedom. The Enskog approximation is formulated and used as the basis to investigate the temperature and the distribution function of the steady state eventually reached by the system. An exact scaling of the distribution function of the system having implications on the form of its moments is pointed out. The theoretical predictions are compared with numerical results obtained by a particle simulation method, and a good agreement is found.

Nonequilibrium identities and response theory for dissipative particles

We derive some nonequilibrium identities such as the integral fluctuation theorem and the Jarzynski equality starting from a nonequilibrium state for dissipative classical systems. Thanks to the existence of the integral fluctuation theorem we can naturally introduce an entropy-like quantity for dissipative classical systems in far from equilibrium states. We also derive the generalized Green-Kubo formula as a nonlinear response theory for a steady dynamics around a nonequilibrium state. We numerically verify the validity of the derived formulas for sheared frictionless granular particles.

Linear hydrodynamics for driven granular gases

We study the dynamics of a granular gas heated by a stochastic thermostat. From a Boltzmann description, we derive the hydrodynamic equations for small perturbations around the stationary state that is reached in the long time limit. Transport coefficients are identified as Green-Kubo formulas obtaining explicit expressions as a function of the inelasticity and the spatial dimension.

Structural responses of quasi-two-dimensional colloidal fluids to excitations elicited by nonequilibrium perturbations

We investigate the response of a dense monodisperse quasi-two-dimensional colloid suspension when a particle is dragged by a constant velocity optical trap. Consistent with microrheological studies of other geometries, the perturbation induces a leading density wave and trailing wake. We also use a hybrid version of Stokesian dynamics simulations to parse direct colloid-colloid and hydrodynamic interactions. We go on to analyze the underlying individual particle-particle collisions in the experimental images. The displacements of particles occur in chains reminiscent of stress propagation in sheared granular materials. From these data, we can reconstruct steady-state dipolar-like flow patterns that were predicted for dilute suspensions and previously observed in granular analogs to our system. The decay of this field differs, however, from point Stokeslet calculations, indicating that the nonzero size of the colloids is important. Moreover, there is a pronounced angular dependence that corresponds to the surrounding colloid structure, which develops in response to the perturbation. Put together, our results show that the response of the complex fluid is highly anisotropic owing to the fact that the effects of the perturbation propagate through the structured medium via chains of colloid-colloid collisions.

Generalized Green-Kubo relation and integral fluctuation theorem for driven dissipative systems without microscopic time reversibility

We derive the generalized Green-Kubo relation and an integral form of the fluctuation theorem that apply to driven dissipative systems, in which time-reversal symmetry or local detailed balance is broken. Uniformly sheared granular systems and driven inelastic Lorentz-gas model are considered as examples. It is discussed how statistical mechanical theory dealing with nonequilibrium steady-state properties can be constructed for such systems for which equilibrium state does not exist.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter =(1-e);{12} , where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the "conduction length." In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k;{23} in the limit k-->0 , if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in . The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t;{-52} in two dimensions and t;{-154} in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t;{-3} in two dimensions and t;{-72} in three dimensions.

Vibrated granular gas confined by a piston

The steady state of a vibrated granular gas confined by a movable piston on top is discussed. Particular attention is given to the hydrodynamic boundary conditions to be used when solving the inelastic Navier-Stokes equations. The relevance of an exact general condition relating the grain fluxes approaching and moving away from each of the walls is emphasized. It is shown how it can be used to get a consistent hydrodynamic description of the boundaries. The obtained expressions for the fields do not contain any undetermined parameter. Comparison of the theoretical predictions with molecular-dynamics simulation results is carried out, and a good agreement is observed for low density and not too large inelasticity. A practical way of introducing small finite density corrections to the dilute limit theory is proposed to improve the accuracy of the theory.

Dynamics of a dilute sheared inelastic fluid. II. The effect of correlations

The effect of correlations on the viscosity of a dilute sheared inelastic fluid is analyzed using the ring-kinetic equation for the two-particle correlation function. The leading-order contribution to the stress in an expansion in =(1-e);{12} is calculated, and it is shown that the leading-order viscosity is identical to that obtained from the Green-Kubo formula, provided the stress autocorrelation function in a sheared steady state is used in the Green-Kubo formula. A systemmatic extension of this to higher orders is also formulated, and the higher-order contributions to the stress from the ring-kinetic equation are determined in terms of the terms in the Chapman-Enskog solution for the Boltzmann equation. The series is resummed analytically to obtain a renormalized stress equation. The most dominant contributions to the two-particle correlation function are products of the eigenvectors of the conserved hydrodynamic modes of the two correlated particles. In Part I, it was shown that the long-time tails of the velocity autocorrelation function are not present in a sheared fluid. Using those results, we show that correlations do not cause a divergence in the transport coefficients; the viscosity is not divergent in two dimensions, and the Burnett coefficients are not divergent in three dimensions. The equations for three-particle and higher correlations are analyzed diagrammatically. It is found that the contributions due to the three-particle and higher correlation functions to the renormalized viscosity are smaller than those due to the two-particle distribution function in the limit -->0 . This implies that the most dominant correlation effects are due to the two-particle correlations.

Transport coefficients for the hard-sphere granular fluid

In the preceding paper, linear response methods have been applied to obtain formally exact expressions for the parameters of Navier-Stokes order hydrodynamics. The analysis there is general, applying to both normal and granular fluids with a wide range of collision rules. Those results are specialized here to the case of smooth, inelastic, hard spheres with constant coefficient of normal restitution, for further elaboration. Explicit expressions for the cooling rate, pressure, and transport coefficients are given and compared with the corresponding expressions for a system of elastic hard spheres. The scope of the results for further analytical explorations and possible numerical evaluation is discussed.