Fate of Boltzmann's Breathers: Stokes Hypothesis and Anomalous Thermalization.

Boltzmann showed that in spite of momentum and energy redistribution through collisions, a rarefied gas confined in a isotropic harmonic trapping potential does not reach equilibrium; it evolves instead into a breathing mode where density, velocity, and temperature oscillate. This counterintuitive prediction is upheld by cold atoms experiments. Yet, are the breathers eternal solutions of the dynamics even in an idealized and isolated system? We show by a combination of hydrodynamic arguments and molecular dynamics simulations that an original dissipative mechanism is at work, where the minute and often neglected bulk viscosity eventually thermalizes the system, which thus reaches equilibrium.

Self-diffusion in a quasi-two-dimensional gas of hard spheres.

A quasi-two-dimensional system of hard spheres strongly confined between two parallel plates is considered. The attention is focused on the macroscopic self-diffusion process observed when the system is seen from above or from below. The transport equation, and the associated self-diffusion coefficient, are derived from a Boltzmann-Lorentz kinetic equation, valid in the dilute limit. Since the equilibrium state of the system is inhomogeneous, this requires the use of a modified Chapman-Enskog expansion that distinguishes between equilibrium and nonequilibrium gradients of the density of labeled particles. The self-diffusion coefficient is obtained as a function of the separation between the two confining plates. The theoretical predictions are compared with molecular dynamics simulation results and a good agreement is found.

Understanding an instability in vibrated granular monolayers.

We investigate the dynamics of an ensemble of smooth inelastic hard spheres confined between two horizontal plates separated by a distance smaller than twice the diameter of the particles, in such a way that the system is quasi-two-dimensional. The bottom wall is vibrating and, therefore, it injects energy into the system in the vertical direction and a stationary state is reached. It is found that if the size of the plates is small enough, the stationary state is homogeneous. Otherwise, a cluster of particles is developed. The instability is understood by using some effective hydrodynamic equations in the horizontal plane. Moreover, the theoretical prediction for the size of the system above which it becomes unstable agrees very well with molecular dynamics simulation results without any fitting parameter.

Kinetic equation and nonequilibrium entropy for a quasi-two-dimensional gas.

A kinetic equation for a dilute gas of hard spheres confined between two parallel plates separated a distance smaller than two particle diameters is derived. It is a Boltzmann-like equation, which incorporates the effect of the confinement on the particle collisions. A function S(t) is constructed by adding to the Boltzmann expression a confinement contribution. Then it is shown that for the solutions of the kinetic equation, S(t) increases monotonically in time, until the system reaches a stationary inhomogeneous state, when S becomes the equilibrium entropy of the confined system as derived from equilibrium statistical mechanics. From the entropy, other equilibrium properties are obtained, and molecular dynamics simulations are used to verify some of the theoretical predictions.

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.

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.

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.

Internal energy fluctuations of a granular gas under steady uniform shear flow.

The stochastic properties of the total internal energy of a dilute granular gas in the steady uniform shear flow state are investigated. A recent theory formulated for fluctuations about the homogeneous cooling state is extended by analogy with molecular systems. The theoretical predictions are compared with molecular dynamics simulation results. Good agreement is found in the limit of weak inelasticity, while systematic and relevant discrepancies are observed when the inelasticity increases. The origin of this behavior is discussed.

Rheological effects in the linear response and spontaneous fluctuations of a sheared granular gas.

The decay of a small homogeneous perturbation in the temperature of a dilute granular gas in the steady uniform shear flow state is investigated. Using kinetic theory based on the inelastic Boltzmann equation, a closed equation for the decay of the perturbation is derived. The equation involves the generalized shear viscosity of the gas in the time-dependent shear flow state, and therefore, it predicts relevant rheological effects beyond the quasielastic limit. Good agreement is found when comparing the theory with molecular dynamics simulation results. Moreover, the Onsager postulate on the regression of fluctuations is fulfilled.

Breakdown of hydrodynamics in the inelastic Maxwell model of granular gases.

Both the right and left eigenfunctions and eigenvalues of the linearized homogeneous Boltzmann equation for inelastic Maxwell molecules corresponding to the hydrodynamic modes are calculated. Also, some nonhydrodynamic modes are identified. It is shown that below a critical value of the parameter characterizing the inelasticity, one of the kinetic modes decays slower than one of the hydrodynamic ones. As a consequence, a closed hydrodynamic description does not exist in that regime. Some implications of this behavior on the formally computed Navier-Stokes transport coefficients are discussed.

Fluctuating hydrodynamics for dilute granular gases.

Starting from the kinetic equations for the fluctuations and correlations of a dilute gas of inelastic hard spheres or disks, a Boltzmann-Langevin equation for the one-particle distribution function of the homogeneous cooling state is constructed. This equation is the linear Boltzmann equation with a fluctuating white-noise term. Balance equations for the fluctuating hydrodynamic fields are derived. New fluctuating forces appear as compared with the elastic limit. The particular case of the transverse velocity field is investigated in detail. Its fluctuations can be described by means of a Langevin equation, but exhibiting two main differences with the Landau-Lifshitz theory: the noise is not white, and its second moment is not determined by the shear viscosity. This shows that the fluctuation-dissipation relations for molecular fluids do not straightforwardly carry over to inelastic gases. The theoretical predictions are shown to be in good agreement with molecular-dynamics simulation results.

Mesoscopic theory of critical fluctuations in isolated granular gases.

Fluctuating hydrodynamics is used to describe the total energy fluctuations of a freely evolving gas of inelastic hard spheres near the threshold of the clustering instability. They are shown to be governed only by vorticity fluctuations that also lead to a renormalization of the average total energy. The theory predicts a power-law divergent behavior of the scaled second moment of the fluctuations, and a scaling property of their probability distribution, both in agreement with simulations results. A more quantitative comparison between theory and simulation for the critical amplitudes and the form of the scaling function is also carried out.

Scaling and universality of critical fluctuations in granular gases.

The total energy fluctuations of a low-density granular gas in the homogeneous cooling state near the threshold of the clustering instability are studied by means of molecular dynamics simulations. The relative dispersion of the fluctuations is shown to exhibit a power-law divergent behavior. Moreover, the probability distribution of the fluctuations presents data collapse as the system approaches the instability, for different values of the inelasticity. The function describing the collapse turns out to be the symmetric of the one found in several molecular equilibrium and nonequilibrium systems.

Energy fluctuations in the homogeneous cooling state of granular gases.

Starting from the hierarchy of equations for microscopic densities in phase space, a general theory for fluctuations and correlations in a dilute granular gas of hard particles is developed. Then, the particular case of the homogeneous cooling state is addressed. Explicit expressions for some distributions describing the presence of velocity correlations and their dynamics are obtained. These correlations are inherent to the dissipative dynamics of the collisions. The implications for the behavior of the total energy of the system are analyzed and the results are expressed in terms of a fluctuation-dissipation theorem. The theoretical predictions are shown to be in agreement with results obtained by molecular dynamics simulations, which also indicate that energy fluctuations are well described by a Gaussian distribution.