Segregation of an intruder in a heated granular dense 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.

Generalized transport coefficients for inelastic Maxwell mixtures under shear flow

The Boltzmann equation framework for inelastic Maxwell models is considered to determine the transport coefficients associated with the mass, momentum, and heat fluxes of a granular binary mixture in spatially inhomogeneous states close to the simple shear flow. The Boltzmann equation is solved by means of a Chapman-Enskog-type expansion around the (local) shear flow distributions f(r)(0) for each species that retain all the hydrodynamic orders in the shear rate. Due to the anisotropy induced by the shear flow, tensorial quantities are required to describe the transport processes instead of the conventional scalar coefficients. These tensors are given in terms of the solutions of a set of coupled equations, which can be analytically solved as functions of the shear rate a, the coefficients of restitution α(rs), and the parameters of the mixture (masses, diameters, and composition). Since the reference distribution functions f(r)(0) apply for arbitrary values of the shear rate and are not restricted to weak dissipation, the corresponding generalized coefficients turn out to be nonlinear functions of both a and α(rs). The dependence of the relevant elements of the three diffusion tensors on both the shear rate and dissipation is illustrated in the tracer limit case, the results showing that the deviation of the generalized transport coefficients from their forms for vanishing shear rates is in general significant. A comparison with the previous results obtained analytically for inelastic hard spheres by using Grad's moment method is carried out, showing a good agreement over a wide range of values for the coefficients of restitution. Finally, as an application of the theoretical expressions derived here for the transport coefficients, thermal diffusion segregation of an intruder immersed in a granular gas is also studied.

Stress partition and microstructure in size-segregating granular flows

When a granular mixture involving grains of different sizes is shaken, sheared, mixed, or left to flow, grains tend to separate by sizes in a process known as size segregation. In this study, we explore the size segregation mechanism in granular chute flows in terms of the pressure distribution and granular microstructure. Therefore, two-dimensional discrete numerical simulations of bidisperse granular chute flows are systematically analyzed. Based on the theoretical models of J. M. N. T. Gray and A. R. Thornton [Proc. R. Soc. A 461, 1447] and K. M. Hill and D. S. Tan [J. Fluid Mech. 756, 54 (2014)], we explore the stress partition in the phases of small and large grains, discriminating between contact stresses and kinetic stresses. Our results support both gravity-induced and shear-gradient-induced segregation mechanisms. However, we show that the contact stress partition is extremely sensitive to the definition of the partial stress tensors and, more specifically, to the way mixed contacts (i.e., involving a small grain and a large grain) are handled, making conclusions on gravity-induced segregation uncertain. By contrast, the computation of the partial kinetic stress tensors is robust. The kinetic pressure partition exhibits a deviation from continuum mixture theory of a significantly higher amplitude than the contact pressure and displays a clear dependence on the flow dynamics. Finally, using a simple approximation for the contact partial stress tensors, we investigate how the contact stress partition relates to the flow microstructure and suggest that the latter may provide an interesting proxy for studying gravity-induced segregation.

Hydrodynamic granular segregation induced by boundary heating and shear

Segregation induced by a thermal gradient of an impurity in a driven low-density granular gas is studied. The system is enclosed between two parallel walls from which we input thermal energy to the gas. We study here steady states occurring when the inelastic cooling is exactly balanced by some external energy input (stochastic force or viscous heating), resulting in a uniform heat flux. A segregation criterion based on Navier-Stokes granular hydrodynamics is written in terms of the tracer diffusion transport coefficients, whose dependence on the parameters of the system (masses, sizes, and coefficients of restitution) is explicitly determined from a solution of the inelastic Boltzmann equation. The theoretical predictions are validated by means of Monte Carlo and molecular dynamics simulations, showing that Navier-Stokes hydrodynamics produces accurate segregation criteria even under strong shearing and/or inelasticity.

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.

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.