Universal reference state in a driven homogeneous granular gas

Kinetic glass transition in granular gases and nonlinear molecular fluids

In this paper, we investigate, both analytically and numerically, the emergence of a kinetic glass transition in two different model systems: a uniformly heated granular gas and a molecular fluid with nonlinear drag. Despite the profound differences between these two physical systems, their behavior in thermal cycles share strong similarities, which stem from the relaxation time diverging algebraically at low temperatures for both systems. When the driving intensity--for the granular gas-or the bath temperature-for the molecular fluid-is decreased to sufficiently low values, the kinetic temperature of both systems becomes "frozen" at a value that depends on the cooling rate through a power law with the same exponent. Interestingly, this frozen glassy state is universal in the following sense: for a suitable rescaling of the relevant variables, its velocity distribution function becomes independent of the cooling rate. Upon reheating, i.e., when either the driving intensity or the bath temperature is increased from this frozen state, hysteresis cycles arise and the apparent heat capacity displays a maximum. The numerical results obtained from the simulations are well described by a perturbative approach.

Enskog kinetic theory of binary granular suspensions: Heat flux and stability analysis of the homogeneous steady state

The Enskog kinetic theory of multicomponent granular suspensions employed previously [Gómez González, Khalil, and Garzó, Phys. Rev. E 101, 012904 (2020)2470-004510.1103/PhysRevE.101.012904] is considered further to determine the four transport coefficients associated with the heat flux. These transport coefficients are obtained by solving the Enskog equation by means of the application of the Chapman-Enskog method around the local version of the homogeneous state. Explicit forms of the heat flux transport coefficients are provided in steady-state conditions by considering the so-called second Sonine approximation to the distribution function of each species. Their quantitative variation on the control parameters of the mixture (masses and diameters, coefficients of restitution, concentration, volume fraction, and the background temperature) is demonstrated and the results show that in general the dependence of the heat flux transport coefficients on inelasticity is clearly different from that found in the absence of the gas phase (dry granular mixtures). As an application of the general results, the stability of the homogeneous steady state is analyzed by solving the linearized Navier-Stokes hydrodynamic equations. The linear stability analysis (which holds for wavelengths long compared with the mean free path) shows that the transversal and longitudinal modes are always stable with respect to long-enough wavelength excitations. This conclusion agrees with previous results derived for monocomponent and (dilute) bidisperse granular suspensions but contrasts with the instabilities found in previous works in dry (no gas phase) granular mixtures.

Optimal Control of Uniformly Heated Granular Fluids in Linear Response

We present a detailed analytical investigation of the optimal control of uniformly heated granular gases in the linear regime. The intensity of the stochastic driving is therefore assumed to be bounded between two values that are close, which limits the possible values of the granular temperature to a correspondingly small interval. Specifically, we are interested in minimising the connection time between the non-equilibrium steady states (NESSs) for two different values of the granular temperature by controlling the time dependence of the driving intensity. The closeness of the initial and target NESSs make it possible to linearise the evolution equations and rigorously-from a mathematical point of view-prove that the optimal controls are of bang-bang type, with only one switching in the first Sonine approximation. We also look into the dependence of the optimal connection time on the bounds of the driving intensity. Moreover, the limits of validity of the linear regime are investigated.

Strong nonexponential relaxation and memory effects in a fluid with nonlinear drag

We analyze the dynamical evolution of a fluid with nonlinear drag, for which binary collisions are elastic, described at the kinetic level by the Enskog-Fokker-Planck equation. This model system, rooted in the theory of nonlinear Brownian motion, displays a really complex behavior when quenched to low temperatures. Its glassy response is controlled by a long-lived nonequilibrium state, independent of the degree of nonlinearity and also of the Brownian-Brownian collisions rate. The latter property entails that this behavior persists in the collisionless case, where the fluid is described by the nonlinear Fokker-Planck equation. The observed response, which includes nonexponential, algebraic, relaxation, and strong memory effects, presents scaling properties: the time evolution of the temperature-for both relaxation and memory effects-falls onto a master curve, regardless of the details of the experiment. To account for the observed behavior in simulations, it is necessary to develop an extended Sonine approximation for the kinetic equation-which considers not only the fourth cumulant but also the sixth one.

Linear response in the uniformly heated granular gas

We analyze the linear response properties of the uniformly heated granular gas. The intensity of the stochastic driving fixes the value of the granular temperature in the nonequilibrium steady state reached by the system. Here, we investigate two specific situations. First, we look into the "direct" relaxation of the system after a single (small) jump of the driving intensity. This study is carried out by two different methods. Not only do we linearize the evolution equations around the steady state, but we also derive generalized out-of-equilibrium fluctuation-dissipation relations for the relevant response functions. Second, we investigate the behavior of the system in a more complex experiment, specifically a Kovacs-like protocol with two jumps in the driving. The emergence of an anomalous Kovacs response is explained in terms of the properties of the direct relaxation function: it is the second mode changing sign at the critical value of the inelasticity that demarcates anomalous from normal behavior. The analytical results are compared with numerical simulations of the kinetic equation, and a good agreement is found.

Enskog kinetic theory for multicomponent granular suspensions

The Navier-Stokes transport coefficients of multicomponent granular suspensions at moderate densities are obtained in the context of the (inelastic) Enskog kinetic theory. The suspension is modeled as an ensemble of solid particles where the influence of the interstitial gas on grains is via a viscous drag force plus a stochastic Langevin-like term defined in terms of a background temperature. In the absence of spatial gradients, it is shown first that the system reaches a homogeneous steady state where the energy lost by inelastic collisions and viscous friction is compensated for by the energy injected by the stochastic force. Once the homogeneous steady state is characterized, a normal solution to the set of Enskog equations is obtained by means of the Chapman-Enskog expansion around the local version of the homogeneous state. To first order in spatial gradients, the Chapman-Enskog solution allows us to identify the Navier-Stokes transport coefficients associated with the mass, momentum, and heat fluxes. In addition, the first-order contributions to the partial temperatures and the cooling rate are also calculated. Explicit forms for the diffusion coefficients, the shear and bulk viscosities, and the first-order contributions to the partial temperatures and the cooling rate are obtained in steady-state conditions by retaining the leading terms in a Sonine polynomial expansion. The results show that the dependence of the transport coefficients on inelasticity is clearly different from that found in its granular counterpart (no gas phase). The present work extends previous theoretical results for dilute multicomponent granular suspensions [Khalil and Garzó, Phys. Rev. E 88, 052201 (2013)10.1103/PhysRevE.88.052201] to higher densities.

When the Hotter Cools More Quickly: Mpemba Effect in Granular Fluids

Under certain conditions, two samples of fluid at different initial temperatures present a counterintuitive behavior known as the Mpemba effect: it is the hotter system that cools sooner. Here, we show that the Mpemba effect is present in granular fluids, both in uniformly heated and in freely cooling systems. In both cases, the system remains homogeneous, and no phase transition is present. Analytical quantitative predictions are given for how differently the system must be initially prepared to observe the Mpemba effect, the theoretical predictions being confirmed by both molecular dynamics and Monte Carlo simulations. Possible implications of our analysis for other systems are also discussed.

Global stability and H theorem in lattice models with nonconservative interactions

In kinetic theory, a system is usually described by its one-particle distribution function f(r,v,t), such that f(r,v,t)drdv is the fraction of particles with positions and velocities in the intervals (r,r+dr) and (v,v+dv), respectively. Therein, global stability and the possible existence of an associated Lyapunov function or H theorem are open problems when nonconservative interactions are present, as in granular fluids. Here, we address this issue in the framework of a lattice model for granularlike velocity fields. For a quite general driving mechanism, including both boundary and bulk driving, we show that the steady state reached by the system in the long-time limit is globally stable. This is done by proving analytically that a certain H functional is nonincreasing in the long-time limit. Moreover, for a quite general energy injection mechanism, we are able to demonstrate that the proposed H functional is nonincreasing for all times. Also, we put forward a proof that clearly illustrates why the "classical" Boltzmann functional H_{B}[f]=∫drdvf(r,v,t)lnf(r,v,t) is inadequate for systems with nonconservative interactions. This is done not only for the simplified kinetic description that holds in the lattice models analyzed here but also for a general kinetic equation, like Boltzmann's or Enskog's.

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.

Statistics of the dissipated energy in driven diffusive systems

Understanding the physics of non-equilibrium systems remains one of the major open questions in statistical physics. This problem can be partially handled by investigating macroscopic fluctuations of key magnitudes that characterise the non-equilibrium behaviour of the system of interest; their statistics, associated structures and microscopic origin. During the last years, some new general and powerful methods have appeared to delve into fluctuating behaviour that have drastically changed the way to address this problem in the realm of diffusive systems: macroscopic fluctuation theory (MFT) and a set of advanced computational techniques that make it possible to measure the probability of rare events. Notwithstanding, a satisfactory theory is still lacking in a particular case of intrinsically non-equilibrium systems, namely those in which energy is not conserved but dissipated continuously in the bulk of the system (e.g. granular media). In this work, we put forward the dissipated energy as a relevant quantity in this case and analyse in a pedagogical way its fluctuations, by making use of a suitable generalisation of macroscopic fluctuation theory to driven dissipative media.

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.

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.

Memory effects in the relaxation of a confined granular gas

The accuracy of a model to describe the horizontal dynamics of a confined quasi-two-dimensional system of inelastic hard spheres is discussed by comparing its predictions for the relaxation of the temperature in a homogenous system with molecular dynamics simulation results for the original system. A reasonably good agreement is found. Next the model is used to investigate the peculiarities of the nonlinear evolution of the temperature when the parameter controlling the energy injection is instantaneously changed while the system was relaxing. This can be considered as a nonequilibrium generalization of the Kovacs effect. It is shown that, in the low-density limit, the effect can be accurately described by using a simple kinetic theory based on the first Sonine approximation for the one-particle distribution function. Some possible experimental implications are indicated.

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.

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.

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.