Mass transport of impurities in a moderately dense granular gas

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.

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.

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.

Instabilities in granular binary mixtures at moderate densities

A linear stability analysis of the Navier-Stokes (NS) granular hydrodynamic equations is performed to determine the critical length scale for the onset of vortices and clusters instabilities in granular dense binary mixtures. In contrast to previous attempts, our results (which are based on the solution to the inelastic Enskog equation to NS order) are not restricted to nearly elastic systems since they take into account the complete nonlinear dependence of the NS transport coefficients on the coefficients of restitution α(ij). The theoretical predictions for the critical length scales are compared to molecular dynamics (MD) simulations in flows of strong dissipation (α(ij) ≥ 0.7) and moderate solid volume fractions (ϕ ≤ 0.2). We find excellent agreement between MD and kinetic theory for the onset of velocity vortices, indicating the applicability of NS hydrodynamics to polydisperse flows even for strong inelasticity, finite density, and particle dissimilarity.

Thermal segregation of intruders in the Fourier state of a granular gas

A low density binary mixture of granular gases is considered within the Boltzmann kinetic theory. One component, the intruders, is taken to be dilute with respect to the other, and thermal segregation of the two species is described for a special solution to the Boltzmann equation. This solution has a macroscopic hydrodynamic representation with a constant temperature gradient and is referred to as the Fourier state. The thermal diffusion factor characterizing conditions for segregation is calculated without the usual restriction to Navier-Stokes hydrodynamics. Integral equations for the coefficients in this hydrodynamic description are calculated approximately within a Sonine polynomial expansion. Molecular dynamics simulations are reported, confirming the existence of this idealized Fourier state. Good agreement is found for the predicted and simulated thermal diffusion coefficient, while only qualitative agreement is found for the temperature ratio.

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.

Segregation by thermal diffusion in moderately dense granular mixtures

A theory based on a solution of the inelastic Enskog equation that goes beyond the weak dissipation limit is used to determine the thermal diffusion factor of a binary granular mixture under gravity. The Enskog equation that aims to describe moderate densities neglects velocity correlations but retains spatial correlations arising from volume exclusion effects. As expected, the thermal diffusion factor provides a segregation criterion that shows the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE) by varying the parameters of the system (masses, sizes, composition, density and coefficients of restitution). The form of the phase diagrams for the BNE/RBNE transition is illustrated in detail in the tracer limit case, showing that the phase diagrams depend sensitively on the value of gravity relative to the thermal gradient. Two specific situations are considered: i) absence of gravity, and ii) homogeneous temperature. In the latter case, after some approximations, our results are consistent with previous theoretical results derived from the Enskog equation. Our results also indicate that the influence of dissipation on thermal diffusion is more important in the absence of gravity than in the opposite limit. The present analysis, which is based on a preliminary short report of the author (Phys. Rev. E 78, 020301(R) (2008)), extends previous theoretical results derived in the dilute limit case.