Hydrodynamic Modes for Granular Gases

The eigenfunctions and eigenvalues of the linearized Boltzmann equation for inelastic hard spheres (d=3) or disks (d=2) corresponding to d+2 hydrodynamic modes, are calculated in the long wavelength limit for a granular gas. The transport coefficients are identified and found to agree with those from the Chapman-Enskog solution. The dominance of hydrodynamic modes at long times and long wavelengths is studied via an exactly solvable kinetic model. A collisional continuum is bounded away from the hydrodynamic spectrum, assuring a hydrodynamic description at long times. The bound is closely related to the power law decay of the velocity distribution in the reference homogeneous cooling state

Kinetic Theory of Response Functions for the Hard Sphere Granular Fluid

The response functions for small spatial perturbations of a homogeneous granular fluid have been described recently. In appropriate dimensionless variables, they have the form of stationary state time correlation functions. Here, these functions are expressed in terms of reduced single particle functions that are expected to obey a linear kinetic equation. The functional assumption required for such a kinetic equation, and a Markov approximation for its implementation are discussed. If, in addition, static velocity correlations are neglected, a granular fluid version of the linearized Enskog kinetic theory is obtained. The derivation makes no a priori limitation on the density, space and time scale, nor degree of inelasticity. As an illustration, recently derived Helfand and Green-Kubo expressions for the Navier-Stokes order transport coefficients are evaluated with this kinetic theory. The results are in agreement with those obtained from the Chapman-Enskog solution to the nonlinear Enskog kinetic equation.Comment: Submitted to J. Stat. Mec

Aging to non-Newtonian hydrodynamics in a granular gas

The evolution to the steady state of a granular gas subject to simple shear flow is analyzed by means of computer simulations. It is found that, regardless of its initial preparation, the system reaches (after a transient period lasting a few collisions per particle) a non-Newtonian (unsteady) hydrodynamic regime, even at strong dissipation and for states where the time scale associated with inelastic cooling is shorter than the one associated with the irreversible fluxes. Comparison with a simplified rheological model shows a good agreement.Comment: 6 pages, 4 figures; v2: improved version to be published in EP

Enskog Theory for Polydisperse Granular Mixtures. I. Navier-Stokes order Transport

A hydrodynamic description for an $s$-component mixture of inelastic, smooth hard disks (two dimensions) or spheres (three dimensions) is derived based on the revised Enskog theory for the single-particle velocity distribution functions. In this first portion of the two-part series, the macroscopic balance equations for mass, momentum, and energy are derived. Constitutive equations are calculated from exact expressions for the fluxes by a Chapman-Enskog expansion carried out to first order in spatial gradients, thereby resulting in a Navier-Stokes order theory. Within this context of small gradients, the theory is applicable to a wide range of restitution coefficients and densities. The resulting integral-differential equations for the zeroth- and first-order approximations of the distribution functions are given in exact form. An approximate solution to these equations is required for practical purposes in order to cast the constitutive quantities as algebraic functions of the macroscopic variables; this task is described in the companion paper.Comment: 36 pages, to be published in Phys. Rev.

Nonlinear viscosity and velocity distribution function in a simple longitudinal flow

A compressible flow characterized by a velocity field $u_x(x,t)=ax/(1+at)$ is analyzed by means of the Boltzmann equation and the Bhatnagar-Gross-Krook kinetic model. The sign of the control parameter (the longitudinal deformation rate $a$) distinguishes between an expansion ($a>0$) and a condensation ($a<0$) phenomenon. The temperature is a decreasing function of time in the former case, while it is an increasing function in the latter. The non-Newtonian behavior of the gas is described by a dimensionless nonlinear viscosity $\eta^*(a^*)$, that depends on the dimensionless longitudinal rate $a^*$. The Chapman-Enskog expansion of $\eta^*$ in powers of $a^*$ is seen to be only asymptotic (except in the case of Maxwell molecules). The velocity distribution function is also studied. At any value of $a^*$, it exhibits an algebraic high-velocity tail that is responsible for the divergence of velocity moments. For sufficiently negative $a^*$, moments of degree four and higher may diverge, while for positive $a^*$ the divergence occurs in moments of degree equal to or larger than eight.Comment: 18 pages (Revtex), including 5 figures (eps). Analysis of the heat flux plus other minor changes added. Revised version accepted for publication in PR

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.Comment: 5 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

Gaussian density fluctuations and Mode Coupling Theory for supercooled liquids

The equations of motion for the density modes of a fluid, derived from Newton's equations, are written as a linear generalized Langevin equation. The constraint imposed by the fluctuation-dissipation theorem is used to derive an exact form for the memory function. The resulting equations, solved under the assumption that the noise, and consequently density fluctuations, of the liquid are gaussian distributed, are equivalent to the random-phase-approximation for the static structure factor and to the well known ideal mode coupling theory (MCT) equations for the dynamics. This finding suggests that MCT is the canonical mean-field theory of the fluid dynamics.Comment: 4 pages, REVTE

Quantum Reciprocity Conjecture for the Non-Equilibrium Steady State

By considering the lack of history dependence in the non-equilibrium steady state of a quantum system we are led to conjecture that in such a system, there is a set of quantum mechanical observables whose retarded response functions are insensitive to the arrow of time, and which consequently satisfy a quantum analog of the Onsager reciprocity relations. Systems which satisfy this conjecture can be described by an effective Free energy functional. We demonstrate that the conjecture holds in a resonant level model of a multi-lead quantum dot.Comment: References revised to take account of related work on Onsager reciprocity in mesoscopics by Christen, and in hydrodynamics by Mclennan, Dufty and Rub

Hydrodynamic modes, Green-Kubo relations, and velocity correlations in dilute granular gases

It is shown that the hydrodynamic modes of a dilute granular gas of inelastic hard spheres can be identified, and calculated in the long wavelength limit. Assuming they dominate at long times, formal expressions for the Navier-Stokes transport coefficients are derived. They can be expressed in a form that generalizes the Green-Kubo relations for molecular systems, and it is shown that they can also be evaluated by means of $N$-particle simulation methods. The form of the hydrodynamic modes to zeroth order in the gradients is used to detect the presence of inherent velocity correlations in the homogeneous cooling state, even in the low density limit. They manifest themselves in the fluctuations of the total energy of the system. The theoretical predictions are shown to be in agreement with molecular dynamics simulations. Relevant related questions deserving further attention are pointed out