Model for the atomic-scale structure of the homogeneous cooling state of granular fluids

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.

Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. II. Stability analysis

Conditions for the stability under linear perturbations around the homogeneous cooling state are studied for dilute granular gases of inelastic and rough hard disks or spheres with constant coefficients of normal (α) and tangential (β) restitution. After a formally exact linear stability analysis of the Navier-Stokes-Fourier hydrodynamic equations in terms of the translational (d_{t}) and rotational (d_{r}) degrees of freedom, the transport coefficients derived in the companion paper [A. Megías and A. Santos, "Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. I. Transport coefficients" Phys. Rev. E 104, 034901 (2021)10.1103/PhysRevE.104.034901] are employed. Known results for hard spheres [Garzó, Santos, and Kremer, Phys. Rev. E 97, 052901 (2018)10.1103/PhysRevE.97.052901] are recovered by setting d_{t}=d_{r}=3, while novel results for hard disks (d_{t}=2, d_{r}=1) are obtained. In the latter case, a high-inelasticity peculiar region in the (α,β) parameter space is found, inside which the critical wave number associated with the longitudinal modes diverges. Comparison with event-driven molecular dynamics simulations for dilute systems of hard disks at α=0.2 shows that this theoretical region of absolute instability may be an artifact of the extrapolation to high inelasticity of the approximations made in the derivation of the transport coefficients, although it signals a shrinking of the conditions for stability. In the case of moderate inelasticity (α=0.7), however, a good agreement between the theoretical predictions and the simulation results is found.

Transport coefficients of solid particles immersed in a viscous gas

Transport properties of a suspension of solid particles in a viscous gas are studied. The dissipation in such systems arises from two sources: inelasticity in particle collisions and viscous dissipation due to the effect of the gas phase on the particles. Here we consider a simplified case in which the mean relative velocity between the gas and solid phases is taken to be zero, such that "thermal drag" is the only remaining gas-solid interaction. Unlike the previous, more general, treatment of the drag force [Garzó et al., J. Fluid Mech. 712, 129 (2012)]JFLSA70022-112010.1017/jfm.2012.404, here we take into account contributions to the (scaled) transport coefficients η^{*} (shear viscosity), κ^{*} (thermal conductivity), and μ^{*} (Dufour-like coefficient) coming from the temperature dependence of the (dimensionless) friction coefficient γ^{*} characterizing the amplitude of the drag force. At moderate densities, the thermal drag model (which is based on the Enskog kinetic equation) is solved by means of the Chapman-Enskog method and the Navier-Stokes transport coefficients are determined in terms of the coefficient of restitution, the solid volume fraction, and the friction coefficient. The results indicate that the effect of the gas phase on η^{*} and μ^{*} is non-negligible (especially in the case of relatively dilute systems) while the form of κ^{*} is the same as the one obtained in the dry granular limit. Finally, as an application of these results, a linear stability analysis of the hydrodynamic equations is carried out to analyze the conditions for stability of the homogeneous cooling state. A comparison with direct numerical simulations shows a good agreement for conditions of practical interest.

Anomalous self-diffusion in a freely evolving granular gas near the shearing instability

The self-diffusion coefficient of a granular gas in the homogeneous cooling state is analyzed near the shearing instability. Using mode-coupling theory, it is shown that the coefficient diverges logarithmically as the instability is approached, due to the coupling of the diffusion process with the shear modes. The divergent behavior, which is peculiar in granular gases and disappears in the elastic limit, does not depend on any other transport coefficient. The theoretical prediction is confirmed by molecular dynamics simulation results for two-dimensional systems.

Quantifying non-ergodic dynamics of force-free granular gases

We demonstrate how non-ergodicity arises in simple mechanistic systems such as force free, dissipative granular gases. This behaviour results from the strong non-stationarity of the process mirrored in the continuous decay of the gas temperature.

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.

Memory effect in uniformly heated granular gases

We evidence a Kovacs-like memory effect in a uniformly driven granular gas. A system of inelastic hard particles, in the low density limit, can reach a nonequilibrium steady state when properly forced. By following a certain protocol for the drive time dependence, we prepare the gas in a state where the granular temperature coincides with its long time value. The temperature subsequently does not remain constant but exhibits a nonmonotonic evolution with either a maximum or a minimum, depending on the dissipation and on the protocol. We present a theoretical analysis of this memory effect at Boltzmann-Fokker-Planck equation level and show that when dissipation exceeds a threshold, the response can be called anomalous. We find excellent agreement between the analytical predictions and direct Monte Carlo simulations.

Kovacs-like memory effect in driven granular gases

While memory effects have been reported for dense enough disordered systems such as glasses, we show here by a combination of analytical and simulation techniques that they are also intrinsic to the dynamics of dilute granular gases. By means of a certain driving protocol, we prepare the gas in a state where the granular temperature T coincides with its long time limit. However, T does not subsequently remain constant but exhibits a nonmonotonic evolution before reaching its nonequilibrium steady value. The corresponding so-called Kovacs hump displays a normal behavior for weak dissipation (as observed in molecular systems) but is reversed under strong dissipation, where it, thus, becomes anomalous.

Shearing instability of a dilute granular mixture

The shearing instability of a dilute granular mixture composed of smooth inelastic hard spheres or disks is investigated. By using the Navier-Stokes hydrodynamic equations, it is shown that the scaled transversal velocity mode exhibits a divergent behavior, similarly to what happens in one-component systems. The theoretical prediction for the critical size is compared with direct Monte Carlo simulations of the Boltzmann equations describing the system, and a good agreement is found. The total energy fluctuations in the vicinity of the transition are shown to scale with the second moment of the distribution. The scaling distribution function is the same as found in other equilibrium and nonequilibrium phase transitions, suggesting the existence of some kind of universality.

Universal reference state in a driven homogeneous granular gas

We study the dynamics of a homogeneous granular gas heated by a stochastic thermostat, in the low density limit. It is found that, before reaching the stationary regime, the system quickly "forgets" the initial condition and then evolves through a universal state that does not only depend on the dimensionless velocity, but also on the instantaneous temperature, suitably renormalized by its steady state value. We find excellent agreement between the theoretical predictions at the Boltzmann equation level for the one-particle distribution function and the direct Monte Carlo simulations. We conclude that at variance with the homogeneous cooling phenomenology, the velocity statistics should not be envisioned as a single-parameter, but as a two-parameter scaling form, keeping track of the distance to stationarity.

Fluctuating hydrodynamics for dilute granular gases: a Monte Carlo study

We investigate hydrodynamic noise in a dilute granular gas during the homogeneous cooling state, by means of a proper application of the direct simulation Monte Carlo (DSMC) algorithm. The DSMC includes a source of randomization which is not present in molecular dynamics (MD) for inelastic hard disks. Notwithstanding this difference, a fair quantitative agreement is found, including a violation of the fluctuation-dissipation relation for the noise amplitude of the same order observed in MD. This study suggests that deterministic collision dynamics is not an essential ingredient to reproduce, up to a good degree of approximation, hydrodynamic fluctuations in dilute granular gases.

Spatial correlations in sheared isothermal liquids: from elastic particles to granular particles

Spatial correlations in sheared isothermal liquids for both elastic and granular cases are theoretically investigated. Using the generalized fluctuating hydrodynamics, correlation functions for both the microscopic scale and the macroscopic scale are obtained. We find the existence of long-range correlations obeying power laws. The validity of our theoretical predictions has been verified from molecular-dynamics simulation.

Shear state of freely evolving granular gases

Hydrodynamic equations are used to identify the final state reached by a freely evolving granular gas above but close to its shear instability. The theory predicts the formation of a two bands shear state with a steady density profile. There is a modulation between temperature and density profiles as a consequence of the energy balance, the density fluctuations remaining small, without producing clustering. Moreover, the time dependence of the velocity field can be scaled out with the squared root of the average temperature of the system. The latter follows the Haff law, but with an effective cooling rate that is smaller than that of the free homogeneous state. The theoretical predictions are compared with numerical results for inelastic hard disks obtained by using the direct Monte Carlo simulation method, and a good agreement is obtained for low inelasticity.

Segregation induced by inelasticity in a vibrofluidized granular mixture

We investigate the segregation of a dense binary mixture of granular particles that only differ in their restitution coefficient. The mixture is vertically vibrated in the presence of gravity. We find a partial segregation of the species, where most dissipative particles submerge in the less dissipative ones. The segregation occurs even if one type of the particles is elastic. In order to have a complete description of the system, we study the structure of the fluid at microscopic scale (few particle diameters). The density and temperature pair distribution functions show strong enhancements with respect to the equilibrium ones at the same density. In particular, there is an increase in the probability that the more inelastic particles group together in pairs (microsegregation). Microscopically the segregation is buoyancy driven, by the appearance of a dense and cold region around the more inelastic particles.

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.

Linear response and hydrodynamics for granular fluids

A formal derivation of linear hydrodynamics for a granular fluid is given. The linear response to small spatial perturbations of a homogeneous reference state is studied in detail, using methods of nonequilibrium statistical mechanics. A transport matrix for macroscopic excitations in the fluid is defined in terms of the response functions. An expansion in the wave vector to second order allows identification of all phenomenological susceptibilities and transport coefficients through Navier-Stokes order in terms of appropriate time correlation functions. The transport coefficients in this representation are the generalization to granular fluids of the familiar Helfand and Green-Kubo relations for normal fluids. The analysis applies to a variety of collision rules. Important differences in both the analysis and results from those for normal fluids are identified and discussed. A scaling limit is described corresponding to the conditions under which idealized inelastic hard sphere models can apply. Further details and interpretation are provided in the paper following this one, by specialization to the case of smooth, inelastic hard spheres with constant coefficient of restitution.

Long-time tails in freely cooling granular gases

The long-time behavior of the current autocorrelation functions for the velocity, the shear stress, and the heat flux is investigated in freely cooling granular gases. It is found that the correlation functions for the velocity and the shear stress have the long-time tails obeying tau(-d/2), while the correlation function for heat flux decays as tau(-(d+2)/2) exp(-zeta*tau) with the dimensionless cooling rate zeta*, the spatial dimension d, and the scaled time tau in terms of the collision frequency. The result of our numerical simulation of the freely cooling granular gases is consistent with the theoretical prediction.

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 part 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.

Transport properties of dense dissipative hard-sphere fluids for arbitrary energy loss models

The revised Enskog approximation for a fluid of hard spheres which lose energy upon collision is discussed for the case that the energy is lost from the normal component of the velocity at collision but is otherwise arbitrary. Granular fluids with a velocity-dependent coefficient of restitution are an important special case covered by this model. A normal solution to the Enskog equation is developed using the Chapman-Enskog expansion. The lowest order solution describes the general homogeneous cooling state and a generating function formalism is introduced for the determination of the distribution function. The first order solution, evaluated in the lowest Sonine approximation, provides estimates for the transport coefficients for the Navier-Stokes hydrodynamic description. All calculations are performed in an arbitrary number of dimensions.

Instabilities in a free granular fluid described by the Enskog equation

A linear stability analysis of the hydrodynamic equations with respect to the homogeneous cooling state is carried out to identify the conditions for stability as functions of the wave vector, the dissipation, and the density. In contrast to previous studies, this description is based on the results derived from the Enskog equation for inelastic hard spheres [V. Garzó and J. W. Dufty, Phys. Rev. E 59, 5895 (1999)], which takes into account the dependence of the transport coefficients on dissipation. As expected, linear stability shows two transversal (shear) modes and a longitudinal "heat") mode to be unstable with respect to long enough wavelength excitations. Comparison with previous results (which neglect the influence of dissipation on transport) shows quantitative discrepancies for strong dissipation.

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.

Rheology of dense polydisperse granular fluids under shear

The solution of the Enskog equation for the one-body velocity distribution of a moderately dense arbitrary mixture of inelastic hard spheres undergoing planar shear flow is described. A generalization of the Grad moment method, implemented by means of a novel generating function technique, is used so as to avoid any assumptions concerning the size of the shear rate. The result is illustrated by using it to calculate the pressure, normal stresses, and shear viscosity of a model polydisperse granular fluid in which grain size, mass, and coefficient of restitution vary among the grains. The results are compared to a numerical solution of the Enskog equation as well as molecular-dynamics simulations. Most bulk properties are well described by the Enskog theory and it is shown that the generalized moment method is more accurate than the simple (Grad) moment method. However, the description of the distribution of temperatures in the mixture predicted by Enskog theory does not compare well to simulation, even at relatively modest densities.

Simulation study of the Green-Kubo relations for dilute granular gases

The Green-Kubo relations for dilute granular gases are employed to compute their transport coefficients by means of the direct simulation Monte Carlo method. This requires not only to follow the dynamics of the system, but also to identify some modified fluxes appearing in the time-correlation functions. The results are compared with those obtained from the Boltzmann equation by means of the Chapman-Enskog procedure in the first Sonine approximation. A good agreement is found for the shear viscosity over a wide range of inelasticities. Nevertheless, for the two transport coefficients associated with the heat flux, significant discrepancies appear for strong inelasticity. Their origin is discussed, showing that they are partially due to the presence of velocity correlations in the homogeneous cooling state of a dilute granular fluid.

Steady-state representation of the homogeneous cooling state of a granular gas

The properties of a dilute granular gas in the homogeneous cooling state are mapped to those of a stationary state by means of a change in the time scale that does not involve any internal property of the system. The new representation is closely related with a general property of the granular temperature in the long time limit. The physical and practical implications of the mapping are discussed. In particular, simulation results obtained by the direct simulation Monte Carlo method applied to the scaled dynamics are reported. This includes ensemble averages and also the velocity autocorrelation function, as well as the self-diffusion coefficient obtained from the latter by means of the Green-Kubo representation. In all cases, the obtained results are compared with theoretical predictions.

Kinetic theory and hydrodynamics of dense, reacting fluids far from equilibrium

The kinetic theory for a fluid of hard spheres which undergo endothermic and/or exothermic reactions with mass transfer is developed. The exact balance equations for concentration, density, velocity, and temperature are derived. The Enskog approximation is discussed and used as the basis for the derivation, via the Chapman-Enskog procedure, of the Navier-Stokes reaction equations under various assumptions about the speed of the chemical reactions. It is shown that the phenomenological description consisting of a reaction-diffusion equation with a convective coupling to the Navier-Stokes equations is of limited applicability.

Validity of the Boltzmann equation to describe low-density granular systems

The departure of a granular gas in the instable region of parameters from the initial homogeneous cooling state is studied. Results from molecular dynamics and from direct Monte Carlo simulation of the Boltzmann equation are compared. The results indicate that the Boltzmann equation accurately predicts the low-density limit of the system. The relevant role played by the parallelization of the velocities as time proceeds and the dependence of this effect on the density is analyzed in detail.

Atomic-scale structure of hard-core fluids under shear flow

The effect of velocity correlations on the equal-time density autocorrelation function, e.g., the pair distribution function (PDF), of a hard-sphere fluid undergoing shear flow is investigated. The PDF at contact is calculated within the Enskog approximation and is shown to be in good agreement with molecular dynamics simulations for shear rates below the shear-induced ordering transition. These calculations are used to construct a nonequilibrium generalized mean-spherical approximation for the PDF at finite separations, which is also found to agree well with the simulation data.

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.

Diffusion in a granular fluid. II. Simulation

The linear-response description for impurity diffusion in a granular fluid undergoing homogeneous cooling is developed in the preceding paper. The formally exact Einstein and Green-Kubo expressions for the self-diffusion coefficient are evaluated there from an approximation to the velocity autocorrelation function. These results are compared here to those from molecular-dynamics simulations over a wide range of density and inelasticity, for the particular case of self-diffusion. It is found that the approximate theory is in good agreement with simulation data up to moderate densities and degrees of inelasticity. At higher density, the effects of inelasticity are stronger, leading to a significant enhancement of the diffusion coefficient over its value for elastic collisions. Possible explanations associated with an unstable long wavelength shear mode are explored, including the effects of strong fluctuations and mode coupling.

Nonequilibrium phase transition for a heavy particle in a granular fluid

It is shown that the homogeneous cooling state (HCS) for a heavy impurity particle in a granular fluid supports two distinct phases. The order parameter straight phi;(s) is the mean square velocity of the impurity particle relative to that of a fluid particle, and the control parameter xi* is the fluid cooling rate relative to the impurity collision rate. For xi*<1 there is a "normal" phase for which straight phi;(s) scales as the fluid/impurity mass ratio, just as for a system with elastic collisions. For xi*>1 an "ordered" phase occurs in which straight phi;(s) is finite even for vanishingly small mass ratio, representing an extreme violation of energy equipartition. The phenomenon can be described in terms of a Landau-like free energy for a second order phase transition. The dynamics leading to the HCS is studied in detail using an asymptotic analysis of the Enskog-Lorentz kinetic equation near each phase and the critical domain. Critical slowing is observed with a divergent relaxation time at the critical point. The stationary velocity distributions are determined in each case, showing a crossover from Maxwellian in the normal phase to an exponential quartic function of the velocity that is sharply peaked about the nonzero straight phi;(s) for the ordered phase. It is shown that the diffusion coefficient in the normal phase diverges at the critical point and remains so in the ordered phase. This is interpreted as a transition from diffusive to ballistic dynamics between the normal and ordered phases.