Enskog kinetic theory of rheology for a moderately dense inertial suspension

Diffusion of intruders in granular suspensions: Enskog theory and random walk interpretation

The Enskog kinetic theory is applied to compute the mean square displacement of impurities or intruders (modeled as smooth inelastic hard spheres) immersed in a granular gas of smooth inelastic hard spheres (grains). Both species (intruders and grains) are surrounded by an interstitial molecular gas (background) that plays the role of a thermal bath. The influence of the latter on the motion of intruders and grains is modeled via a standard viscous drag force supplemented by a stochastic Langevin-like force proportional to the background temperature. We solve the corresponding Enskog-Lorentz kinetic equation by means of the Chapman-Enskog expansion truncated to first order in the gradient of the intruder number density. The integral equation for the diffusion coefficient is solved by considering the first two Sonine approximations. To test these results, we also compute the diffusion coefficient from the numerical solution of the inelastic Enskog equation by means of the direct simulation Monte Carlo method. We find that the first Sonine approximation generally agrees well with the simulation results, although significant discrepancies arise when the intruders become lighter than the grains. Such discrepancies are largely mitigated by the use of the second Sonine approximation, in excellent agreement with computer simulations even for moderately strong inelasticities and/or dissimilar mass and diameter ratios. We invoke a random walk picture of the intruders' motion to shed light on the physics underlying the intricate dependence of the diffusion coefficient on the main system parameters. This approach, recently employed to study the case of an intruder immersed in a granular gas, also proves useful in the present case of a granular suspension. Finally, we discuss the applicability of our model to real systems in the self-diffusion case. We conclude that collisional effects may strongly impact the diffusion coefficient of the grains.

Homogeneous cooling and heating states of dilute soft-core gases under nonlinear drag

The temperature evolution of dilute soft inertial gas-solid suspensions is theoretically analyzed when the gas particles are influenced by a nonlinear drag force from a background fluid. The kinetic theory is extended to this system, and the time evolutions of the temperature and the kurtosis of the velocity distribution are derived. Molecular dynamics simulations are also performed to check the validity of the theory, and they show good agreement with the theoretical predictions.

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.