Hydrodynamic correlation functions of a driven granular fluid in steady state

Pseudo hard-sphere viscosities from equilibrium Molecular Dynamics

Transport coefficients like shear, bulk and longitudinal viscosities are sensitive to the intermolecular interaction potential and finite size effects when are numerically determined. For the hard-sphere (HS) fluid, such transport properties are determined almost exclusively with computer simulations. However, their systematic determination and analysis throughout shear stress correlation functions and the Green-Kubo formalism can not be done due to discontinuous nature of the interaction potential. Here, we use the pseudo hard-sphere (PHS) potential to determine pressure correlation functions as a function of volume fraction in order to compute mentioned viscosities. Simulation results are compared to available event-driven molecular dynamics of the HS fluid and also used to propose empirical corrections for the Chapman-Enskog zero density limit of shear viscosity. Moreover, we show that PHS potential is a reliable representation of the HS fluid and can be used to compute transport coefficients. The molecular simulation results of the present work are valuable for further exploration of HS-type fluids or extend the approach to compute transport properties of hard-colloid suspensions.

Wave spectroscopy in a driven granular material

Driven granular media constitute model systems in out-of-equilibrium statistical physics. By assimilating the motions of granular particles to those of atoms, by analogy, one can obtain macroscopic equivalent of phase transitions. Here, we study fluid-like and crystal-like two-dimensional states in a driven granular material. In our experimental device, a tunable magnetic field induces and controls remote interactions between the granular particles. We use high-speed video recordings to analyse the velocity fluctuations of individual particles in stationary regime. Using statistical averaging, we find that the particles self-organize into collective excitations characterized by dispersion relations in the frequency-wavenumber space. These findings thus reveal that mechanical waves analogous to condensed matter phonons propagate in driven granular media. When the magnetic coupling is weak, the waves are longitudinal, as expected for a fluid-like phase. When the coupling is stronger, both longitudinal and transverse waves propagate, which is typically seen in solid-like phases. We model the dispersion relations using the spatial distribution of particles and their interaction potential. Finally, we infer the elastic parameters of the granular assembly from equivalent sound velocities, thus realizing the spectroscopy of a granular material.

Reticulate collisional structure in boundary-driven granular gases

We report a peculiar head-on collision network between two vibrating boundaries in experiments performed during a parabolic flight and in a laboratory using horizontal vibration. This structure is a new ordering, which is due to an orientation correlation between the relative position and velocity of any particle pair. It weakens the collision frequency and produces a long-range boundary effect. Moreover, we find the molecular chaos assumption is violated in a larger portion of the phase space. Using an anisotropic distribution model, we modify angular integration results and compare them to the results of the kinetic theory.

Heat flux of driven granular mixtures at low density: Stability analysis of the homogeneous steady state

The Navier-Stokes order hydrodynamic equations for a low-density driven granular mixture obtained previously [Khalil and Garzó, Phys. Rev. E 88, 052201 (2013)PLEEE81539-375510.1103/PhysRevE.88.052201] from the Chapman-Enskog solution to the Boltzmann equation are considered further. The four transport coefficients associated with the heat flux are obtained in terms of the mass ratio, the size ratio, composition, coefficients of restitution, and the driven parameters of the model. Their quantitative variation on the control parameters of the system is demonstrated by considering the leading terms in a Sonine polynomial expansion to solve the exact integral equations. As an application of these results, the stability of the homogeneous steady state is studied. In contrast to the results obtained in undriven granular mixtures, the stability analysis of the linearized Navier-Stokes hydrodynamic equations shows that the transversal and longitudinal modes are (linearly) stable with respect to long enough wavelength excitations. This conclusion agrees with a previous analysis made for single granular gases.

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.

Transport coefficients for driven granular mixtures at low density

The transport coefficients of a granular binary mixture driven by a stochastic bath with friction are determined from the inelastic Boltzmann kinetic equation. A normal solution is obtained via the Chapman-Enskog method for states near homogeneous steady states. The mass, momentum, and heat fluxes are determined to first order in the spatial gradients of the hydrodynamic fields, and the associated transport coefficients are identified. They are given in terms of the solutions of a set of coupled linear integral equations. As in the monocomponent case, since the collisional cooling cannot be compensated for locally by the heat produced by the external driving, the reference distributions (zeroth-order approximations) f(i)((0)) (i=1,2) for each species depend on time through their dependence on the pressure and the temperature. Explicit forms for the diffusion transport coefficients and the shear viscosity coefficient are obtained by assuming the steady-state conditions and by considering the leading terms in a Sonine polynomial expansion. A comparison with previous results obtained for granular Brownian motion and by using a (local) stochastic thermostat is also carried out. The present work extends previous theoretical results derived for monocomponent dense gases [Garzó, Chamorro, and Vega Reyes, Phys. Rev. E 87, 032201 (2013)] to granular mixtures at low density.

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.

Hydrodynamic modes in a confined granular fluid

Confined granular fluids, placed in a shallow box that is vibrated vertically, can achieve homogeneous stationary states due to energy injection mechanisms that take place throughout the system. These states can be stable even at high densities and inelasticities allowing for a detailed analysis of the hydrodynamic modes that govern the dynamics of granular fluids. By analyzing the decay of the time correlation functions it is shown that there is a crossover from a quasielastic regime in which energy evolves as a slow mode to an inelastic regime with energy slaved to the other conserved fields. The two regimes have well differentiated transport properties and in the inelastic regime the dynamics can be described by a reduced hydrodynamics with modified longitudinal viscosity and sound speed. The crossover between the two regimes takes place at a wave vector that is proportional to the inelasticity. A two-dimensional granular model, with collisions that mimic the energy transfers that take place in a confined system, is studied by means of microscopic simulations. The results show excellent agreement with the theoretical framework and allow validation of hydrodynamiclike models.

Longitudinal viscosity of two-dimensional Yukawa liquids

The longitudinal viscosity η(l) is obtained for a two-dimensional (2D) liquid using a Green-Kubo method with a molecular dynamics simulation. The interparticle potential used has the Debye-Hückel or Yukawa form, which models a 2D dusty plasma. The longitudinal η(l) and shear η(s) viscosities are found to have values that match very closely, with only negligible differences for the entire range of temperatures that is considered. For a 2D Yukawa liquid, the bulk viscosity η(b) is determined to be either negligibly small or not a meaningful transport coefficient.

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.

Segregation of an intruder in a heated granular dense gas

A recent segregation criterion [Phys. Rev. E 78, 020301(R) (2008)] based on the thermal diffusion factor Λ of an intruder in a heated granular gas described by the inelastic Enskog equation is revisited. The sign of Λ provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE). The present theory incorporates two extra ingredients not accounted for by the previous theoretical attempt. First, the theory is based upon the second Sonine approximation to the transport coefficients of the mass flux of the intruder. Second, the dependence of the temperature ratio (intruder temperature over that of the host granular gas) on the solid volume fraction is taken into account in the first and second Sonine approximations. In order to check the accuracy of the Sonine approximation considered, the Enskog equation is also numerically solved by means of the direct simulation Monte Carlo method to get the kinetic diffusion coefficient D(0). The comparison between theory and simulation shows that the second Sonine approximation to D(0) yields an improvement over the first Sonine approximation when the intruder is lighter than the gas particles in the range of large inelasticity. With respect to the form of the phase diagrams for the BNE-RBNE transition, the kinetic theory results for the factor Λ indicate that while the form of these diagrams depends sensitively on the order of the Sonine approximation considered when gravity is absent, no significant differences between both Sonine solutions appear in the opposite limit (gravity dominates the thermal gradient). In the former case (no gravity), the first Sonine approximation overestimates both the RBNE region and the influence of dissipation on thermal diffusion segregation.

Slow dynamics and precursors of the glass transition in granular fluids

We use event driven simulations to analyze glassy dynamics as a function of density and energy dissipation in a two-dimensional bidisperse granular fluid under stationary conditions. Clear signatures of a glass transition are identified, such as an increase of relaxation times over several orders of magnitude. As the inelasticity is increased, the glass transition is shifted to higher densities, and the precursors of the transition become less and less pronounced, in agreement with a recent mode-coupling theory. We analyze the long-time tails of the velocity autocorrelation and discuss its consequences for the nonexistence of the diffusion constant in two dimensions.

Transport coefficients of driven granular fluids at moderate volume fraction

In a recent publication [Phys. Rev. E 83, 011301 (2011)], Vollmayr-Lee et al. have determined by computer simulations the thermal diffusivity and the longitudinal viscosity coefficients of a driven granular fluid of hard spheres at intermediate volume fractions. Although they compare their simulation results with the predictions of kinetic theory, they use the dilute expressions for the driven system and the modified Sonine approximations for the undriven system. The goal here is to carry out this comparison by proposing a modified Sonine approximation to the Enskog equation for driven granular fluids that leads to a better quantitative agreement.