Transport coefficients for driven granular mixtures at low density

Enskog kinetic theory of binary granular suspensions: Heat flux and stability analysis of the homogeneous steady state

The Enskog kinetic theory of multicomponent granular suspensions employed previously [Gómez González, Khalil, and Garzó, Phys. Rev. E 101, 012904 (2020)2470-004510.1103/PhysRevE.101.012904] is considered further to determine the four transport coefficients associated with the heat flux. These transport coefficients are obtained by solving the Enskog equation by means of the application of the Chapman-Enskog method around the local version of the homogeneous state. Explicit forms of the heat flux transport coefficients are provided in steady-state conditions by considering the so-called second Sonine approximation to the distribution function of each species. Their quantitative variation on the control parameters of the mixture (masses and diameters, coefficients of restitution, concentration, volume fraction, and the background temperature) is demonstrated and the results show that in general the dependence of the heat flux transport coefficients on inelasticity is clearly different from that found in the absence of the gas phase (dry granular mixtures). As an application of the general results, the stability of the homogeneous steady state is analyzed by solving the linearized Navier-Stokes hydrodynamic equations. The linear stability analysis (which holds for wavelengths long compared with the mean free path) shows that the transversal and longitudinal modes are always stable with respect to long-enough wavelength excitations. This conclusion agrees with previous results derived for monocomponent and (dilute) bidisperse granular suspensions but contrasts with the instabilities found in previous works in dry (no gas phase) granular mixtures.

Kinetic Theory of Polydisperse Granular Mixtures: Influence of the Partial Temperatures on Transport Properties-A Review

It is well-recognized that granular media under rapid flow conditions can be modeled as a gas of hard spheres with inelastic collisions. At moderate densities, a fundamental basis for the determination of the granular hydrodynamics is provided by the Enskog kinetic equation conveniently adapted to account for inelastic collisions. A surprising result (compared to its molecular gas counterpart) for granular mixtures is the failure of the energy equipartition, even in homogeneous states. This means that the partial temperatures Ti (measuring the mean kinetic energy of each species) are different to the (total) granular temperature T. The goal of this paper is to provide an overview on the effect of different partial temperatures on the transport properties of the mixture. Our analysis addresses first the impact of energy nonequipartition on transport which is only due to the inelastic character of collisions. This effect (which is absent for elastic collisions) is shown to be significant in important problems in granular mixtures such as thermal diffusion segregation. Then, an independent source of energy nonequipartition due to the existence of a divergence of the flow velocity is studied. This effect (which was already analyzed in several pioneering works on dense hard-sphere molecular mixtures) affects to the bulk viscosity coefficient. Analytical (approximate) results are compared against Monte Carlo and molecular dynamics simulations, showing the reliability of kinetic theory for describing granular flows.

Energy nonequipartition in a collisional model of a confined quasi-two-dimensional granular mixture

A collisional model of a confined quasi-two-dimensional granular mixture is considered to analyze homogeneous steady states. The model includes an effective mechanism to transfer the kinetic energy injected by vibration in the vertical direction to the horizontal degrees of freedom of grains. The set of Enskog kinetic equations for the velocity distribution functions of each component is derived first to analyze the homogeneous state. As in the one-component case, an exact scaling solution is found where the time dependence of the distribution functions occurs entirely through the granular temperature T. As expected, the kinetic partial temperatures T_{i} of each component are different and, hence, energy equipartition is broken down. In the steady state, explicit expressions for the temperature T and the ratio of partial kinetic temperatures T_{i}/T_{j} are obtained by considering Maxwellian distributions defined at the partial temperatures T_{i}. The (scaled) granular temperature and the temperature ratios are given in terms of the coefficients of restitution, the solid volume fraction, the (scaled) parameters of the collisional model, and the ratios of mass, concentration, and diameters. In the case of a binary mixture, the theoretical predictions are exhaustively compared with both direct simulation Monte Carlo and molecular dynamics simulations with a good agreement. The deviations are identified to be originated in the non-Gaussianity of the velocity distributions and on microsegregation patterns, which induce spatial correlations not captured in the Enskog theory.

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).

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.

Influence of the first-order contributions to the partial temperatures on transport properties in polydisperse dense granular mixtures

The Chapman-Enskog solution to the Enskog kinetic equation of polydisperse granular mixtures is revisited to determine the first-order contributions ϖ_{i} to the partial temperatures. As expected, these quantities (which were neglected in previous attempts) are given in terms of the solution to a set of coupled integrodifferential equations analogous to those for elastic collisions. The solubility condition for this set of equations is confirmed and the coefficients ϖ_{i} are calculated by using the leading terms in a Sonine polynomial expansion. These coefficients are given as explicit functions of the sizes, masses, composition, density, and coefficients of restitution of the mixture. Within the context of small gradients, the results apply for arbitrary degrees of inelasticity and are not restricted to specific values of the parameters of the mixture. In the case of elastic collisions, previous expressions of ϖ_{i} for ordinary binary mixtures are recovered. Finally, the impact of the first-order coefficients ϖ_{i} on the bulk viscosity η_{b} and on the first-order contribution ζ^{(1)} to the cooling rate is assessed. It is shown that the effect of ϖ_{i} on η_{b} and ζ^{(1)} is not negligible, specially for disparate mass ratios and strong inelasticity.

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.

Shear-rate-dependent transport coefficients in granular suspensions

A recent model for monodisperse granular suspensions is used to analyze transport properties in spatially inhomogeneous states close to the simple (or uniform) shear flow. The kinetic equation is based on the inelastic Boltzmann (for low-density gases) with the presence of a viscous drag force that models the influence of the interstitial gas phase on the dynamics of grains. A normal solution is obtained via a Chapman-Enskog-like expansion around a (local) shear flow distribution which retains all the hydrodynamic orders in the shear rate. To first order in the expansion, the transport coefficients characterizing momentum and heat transport around shear flow are given in terms of the solutions of a set of coupled linear integral equations which are approximately solved by using a kinetic model of the Boltzmann equation. To simplify the analysis, the steady-state conditions when viscous heating is compensated by the cooling terms arising from viscous friction and collisional dissipation are considered to get the explicit forms of the set of generalized transport coefficients. The shear-rate dependence of some of the transport coefficients of the set is illustrated for several values of the coefficient of restitution.

Hydrodynamic Burnett equations for inelastic Maxwell models of granular gases

The hydrodynamic Burnett equations and the associated transport coefficients are exactly evaluated for generalized inelastic Maxwell models. In those models, the one-particle distribution function obeys the inelastic Boltzmann equation, with a velocity-independent collision rate proportional to the γ power of the temperature. The pressure tensor and the heat flux are obtained to second order in the spatial gradients of the hydrodynamic fields with explicit expressions for all the Burnett transport coefficients as functions of γ, the coefficient of normal restitution, and the dimensionality of the system. Some transport coefficients that are related in a simple way in the elastic limit become decoupled in the inelastic case. As a byproduct, existing results in the literature for three-dimensional elastic systems are recovered, and a generalization to any dimension of the system is given. The structure of the present results is used to estimate the Burnett coefficients for inelastic hard spheres.