Mobility and diffusion of intruders in granular suspensions: Einstein relation

The Enskog kinetic equation is considered to determine the diffusion D and mobility λ transport coefficients of intruders immersed in a granular gas of inelastic hard spheres (grains). Intruders and grains are in contact with a thermal bath, which plays the role of a background gas. As usual, the influence of the latter on the dynamics of intruders and grains is accounted for via a viscous drag force plus a stochastic Langevin-like term proportional to the background temperature T b . In this case, the starting kinetic equations are the Enskog and Enskog–Lorentz equations for grains and intruders, respectively, with the addition of Fokker–Planck terms to each one of the above equations. The transport coefficients λ and D are determined by solving the Enskog–Lorentz kinetic equation by means of the Chapman–Enskog method adapted to inelastic collisions. As for elastic collisions, both transport coefficients are given in terms of the solutions of two integral equations, which are approximately solved up to the second order in a Sonine polynomial expansion. Theoretical results are compared against numerical solutions of the inelastic Enskog equation by means of the direct simulation Monte Carlo method. In general, good agreement between theory and simulations is found, especially in the case of the second Sonine approximation. Knowledge of the coefficients λ and D allows us to assess the departure of the (conventional) Einstein relation ϵ = D / ( T b λ ) from 1. As expected from previous results for driven granular gases, it is shown that when the bath temperature T b is replaced by the intruder temperature T 0 in the Einstein relation, the origin of the deviation of ε from 1 is only due to the non-Maxwellian behavior of the reference state of intruders (measured by the cumulant c 0 ). Since the magnitude of c 0 is in general very small, deviations of the (modified) Einstein relation ϵ 0 = D / ( T 0 λ ) from 1 cannot be detected in computer simulations of dilute granular gases. This conclusion agrees well with previous computer simulation results.