Hydrodynamics of driven granular gases

Mpemba effect in inertial suspensions

The Mpemba effect (a counterintuitive thermal relaxation process where an initially hotter system may cool down to the steady state sooner than an initially colder system) is studied in terms of a model of inertial suspensions under shear. The relaxation to a common steady state of a suspension initially prepared in a quasiequilibrium state is compared with that of a suspension initially prepared in a nonequilibrium sheared state. Two classes of Mpemba effect are identified, the normal and the anomalous one. The former is generic, in the sense that the kinetic temperature starting from a cold nonequilibrium sheared state is overtaken by the one starting from a hot quasiequilibrium state, due to the absence of initial viscous heating in the latter, resulting in a faster initial cooling. The anomalous Mpemba effect is opposite to the normal one since, despite the initial slower cooling of the nonequilibrium sheared state, it can eventually overtake an initially colder quasiequilibrium state. The theoretical results based on kinetic theory agree with those obtained from event-driven simulations for inelastic hard spheres. It is also confirmed the existence of the inverse Mpemba effect, which is a peculiar heating process, in these suspensions. More particularly, we find the existence of a mixed process in which both heating and cooling can be observed during relaxation.

Enskog kinetic theory of rheology for a moderately dense inertial suspension

The Enskog kinetic theory for moderately dense inertial suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the background fluid on suspended particles is modeled via a viscous drag force plus a Langevin-like term defined in terms of the background temperature. In a previous paper [Hayakawa et al., Phys. Rev. E 96, 042903 (2017)10.1103/PhysRevE.96.042903], Grad's moment method with the aid of a linear shear-rate expansion was employed to obtain a theory which gave good agreement with the results of event-driven Langevin simulations of hard spheres for low densities and/or small shear rates. Nevertheless, the previous approach had a limitation of not being applicable to the high-shear-rate and high-density regime. Thus, in the present paper, we extend the previous work and develop Grad's theory including higher-order terms in the shear rate. This improves significantly the theoretical predictions, a quantitative agreement between theory and simulation being found in the high-density region (volume fractions smaller than or equal to 0.4).

Kinetic theory of shear thickening for a moderately dense gas-solid suspension: From discontinuous thickening to continuous thickening

The Enskog kinetic theory for moderately dense gas-solid suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the environmental fluid on solid particles is modeled via a viscous drag force plus a stochastic Langevin-like term. The Enskog equation is solved by means of two independent but complementary routes: (i) Grad's moment method and (ii) event-driven Langevin simulation of hard spheres. Both approaches clearly show that the flow curve (stress-strain rate relation) depends significantly on the volume fraction of the solid particles. In particular, as the density increases, there is a transition from the discontinuous shear thickening (observed in dilute gases) to the continuous shear thickening for denser systems. The comparison between theory and simulations indicates that while the theoretical predictions for the kinetic temperature agree well with simulations for densities φ≲0.5, the agreement for the other rheological quantities (the viscosity, the stress ratio, and the normal stress differences) is limited to more moderate densities (φ≲0.3) if the inelasticity during collisions between particles is not large.

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.

Theory of thermostatted inhomogeneous granular fluids: A self-consistent density functional description

The authors present a study of the nonequilibrium statistical properties of a one dimensional hard-rod fluid dissipating energy via inelastic collisions and subject to the action of a Gaussian heat bath, simulating an external driving mechanism. They show that the description of the fluid based on the one-particle phase-space reduced distribution function, in principle necessary because of the presence of velocity dependent collisional dissipation, can be contracted to a simpler description in configurational space. Indeed, by means of a multiple-time-scale method the authors derive a self-consistent governing equation for the particle density distribution function. This equation is similar to the dynamic density functional equation employed in the study of colloids, but contains additional terms taking into account the inelastic nature of the fluid. Such terms cannot be derived from a Liapunov generating functional and contribute not only to the relaxational properties, but also to the nonequilibrium steady state properties. A validation of the theory against molecular dynamics simulations is presented in a series of cases, and good agreement is found.

Local heat flux and energy loss in a two-dimensional vibrated granular gas

We performed event-driven simulations of a two-dimensional granular gas between two vibrating walls and directly measured the local heat flux and local energy dissipation in the stationary state. Describing the local heat flux as a function of the coordinate in the direction perpendicular to the driving walls, we test a generalization of Fourier's law, q(x)=-kappa inverted delta T(x)+mu inverted delta rho(x), by relating the local heat flux to the local gradients of the temperature and density. This ansatz accounts for the fact that heat flux can also be generated by density gradients, not only by temperature gradients. Assuming the transport coefficients kappa and mu to be independent of x, we check the validity of this assumption and test the generalized Fourier law in the simulations. Both kappa and mu are determined for different system parameters, in particular, for a wide range of coefficients of restitution. We also compare our numerical results to existing hydrodynamic theories. Agreement is found for kappa for very small inelasticities only, i.e., when the gradients are small. Beyond this region, kappa and mu exhibit a striking nonmonotonic behavior. This may hint that hydrodynamics to Navier-Stokes order cannot be applied to moderately inelastic vibrated systems.