The Enskog kinetic theory for moderately dense granular suspensions is considered as a model to determine the Navier-Stokes transport coefficients. The influence of the interstitial gas on solid particles is modeled by a viscous drag force term plus a stochastic Langevin-like term. The suspension model is solved by means of the Chapman-Enskog method conveniently adapted to dissipative dynamics. The momentum and heat fluxes as well as the cooling rate are obtained to first order in the deviations of the hydrodynamic field gradients from their values in the homogeneous steady state. Since the cooling terms (arising from collisional dissipation and viscous friction) cannot be compensated for by the energy gained by grains due to collisions with the interstitial gas, the reference distribution (zeroth-order approximation of the Chapman-Enskog solution) depends on time through its dependence on temperature. On the other hand, to simplify the analysis and given that we are interested in computing transport properties in the first order of deviations from the reference state, the steady-state conditions are considered. This simplification allows us to get explicit expressions for the Navier-Stokes transport coefficients. The present work extends previous results [Garzó et al. 2013, Phys Rev. E 87, 032201] since it incorporates two extra ingredients (an additional density dependence of the zeroth-order solution and the density dependence of the reduced friction coefficient) not accounted for by the previous theoretical attempt. While these two new ingredients do not affect the shear viscosity coefficient, the transport coefficients associated with the heat flux as well as the first-order contribution to the cooling rate are different from those obtained in the previous study. In addition, as expected, the results show that the dependence of the transport coefficients on both inelasticity and density is clearly different from that found in its granular counterpart (no gas phase). Finally, a linear stability analysis of the hydrodynamic equations with respect to the homogeneous steady state is performed. In contrast to the granular case (no gas-phase), no instabilities are found and hence, the homogeneous steady state is (linearly) stable. * Electronic address: ruben@unex.es † Electronic address: vicenteg@unex.es; URL: http://www.unex.es/eweb/fisteor/vicente/ Boltzmann equation for inelastic Maxwell models [10]) compare quite well with molecular dynamics simulations [9] for conditions of practical interest. This good agreement highlights again the good performance of kinetic theory tools in reproducing the transport properties of gas-solid flows.On the other hand, to the best of our knowledge, most of the efforts in kinetic theory of granular suspensions has been mainly focused on non-Newtonian transport properties (which are directly related with the pressure tensor). In particular, much less is known about the energy transport in gas-solid flows. The knowledge of the transport coefficients associated with the heat flux ...

The transport coefficients of a dilute gas of inelastic hard spheres immersed in a gas of elastic hard spheres (molecular gas) are determined. We assume that the number density of the granular gas is much smaller than that of the surrounding molecular gas, so that the latter is not affected by the presence of the granular particles. In this situation, the molecular gas may be treated as a thermostat (or bath) of elastic hard spheres at a fixed temperature. The Boltzmann kinetic equation is the starting point of the present work. The first step is to characterise the reference state in the perturbation scheme, namely the homogeneous state. Theoretical results for the granular temperature and kurtosis obtained in the homogeneous steady state are compared against Monte Carlo simulations showing a good agreement. Then, the Chapman–Enskog method is employed to solve the Boltzmann equation to first order in spatial gradients. In dimensionless form, the Navier–Stokes–Fourier transport coefficients of the granular gas are given in terms of the mass ratio $m/m_g$ ( $m$ and $m_g$ being the masses of a granular and a gas particle, respectively), the (reduced) bath temperature and the coefficient of restitution. Interestingly, previous results derived from a suspension model based on an effective fluid–solid interaction force are recovered in the Brownian limit ( $m/m_g \to \infty$ ). Finally, as an application of the theory, a linear stability analysis of the homogeneous steady state is performed showing that this state is always linearly stable.

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)] to higher densities.

The Boltzmann kinetic equation for low-density granular suspensions under simple shear flow is considered to determine the velocity moments through the fourth degree. The influence of the interstitial gas on solid particles is modeled by a viscous drag force term plus a stochastic Langevin-like term. Two independent but complementary approaches are followed to achieve exact results. First, to keep the structure of the Boltzmann collision operator, the so-called inelastic Maxwell models (IMM) are considered. In this model, since the collision rate is independent of the relative velocity of the two colliding particles, the forms of the collisional moments can be obtained without the knowledge of the velocity distribution function. As a complement of the previous effort, a BGK-type kinetic model adapted to granular gases is solved to get the velocity moments of the velocity distribution function. The analytical predictions of the rheological properties (which are exactly obtained in terms of the coefficient of restitution α and the reduced shear rate a * ) show in general an excellent agreement with event-driven simulations performed for inelastic hard spheres. In particular, both theoretical approaches show clearly that the temperature and non-Newtonian viscosity exhibit an S shape in a plane of stress-strain rate (discontinuous shear thickening effect). With respect to the fourth-degree velocity moments, we find that while those moments have unphysical values for IMM in a certain region of the parameter space of the system, they are well defined functions of both α and a * in the case of the BGK kinetic model. The explicit shear-rate dependence of the fourth-degree moments beyond this critical region is also obtained and compared against available computer simulations. * Electronic address: ruben@unex.es † Electronic address: vicenteg@unex.es;

The Mpemba effect occurs when two samples at different initial temperatures evolve in such a way that the temperatures cross each other during the relaxation toward equilibrium. In this paper, we show the emergence of a Mpemba-like effect in a molecular binary mixture in contact with a thermal reservoir (bath). The interaction between the gaseous particles of the mixture and the thermal reservoir is modeled via a viscous drag force plus a stochastic Langevin-like term. The presence of the external bath couples the time evolution of the total and partial temperatures of each component allowing the appearance of the Mpemba phenomenon, even when the initial temperature differences are of the same order of the temperatures themselves. Analytical results are obtained by considering multitemperature Maxwellian approximations for the velocity distribution functions of each component. The theoretical analysis is carried out for initial states close to and far away (large Mpemba-like effect) from equilibrium. The former situation allows us to develop a simple theory where the time evolution equation for the temperature is linearized around its asymptotic equilibrium solution. This linear theory provides an expression for the crossover time. We also provide a qualitative description of the large Mpemba effect. Our theoretical results agree very well with computer simulations obtained by numerically solving the Enskog kinetic equation by means of the direct simulation Monte Carlo method and by performing molecular dynamics simulations. Finally, preliminary results for driven granular mixtures also show the occurrence of a Mpemba-like effect for inelastic collisions.

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 integro-differential 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 degree 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 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.

The time evolution of a homogeneous bidisperse granular suspension is studied in the context of the Enskog kinetic equation. The influence of the surrounding viscous gas on the solid particles is modeled via a deterministic viscous drag force plus a stochastic Langevin-like term. It is found first that, regardless of the initial conditions, the system reaches (after a transient period lasting a few collisions per particle) a universal unsteady hydrodynamic regime where the distribution function of each species not only depends on the dimensionless velocity (as in the homogeneous cooling state) but also on the instantaneous temperature scaled with respect to the background temperature. To confirm this result, theoretical predictions for the time-dependent partial temperatures are compared against direct simulation Monte Carlo (DSMC) results; the comparison shows an excellent agreement confirming the applicability of hydrodynamics in granular suspensions. Also, in the transient regime, the so-called Mpemba-like effect (namely, when an initially hotter sample cools sooner than the colder one) is analyzed for inelastic collisions. The theoretical analysis of the Mpemba effect is performed for initial states close to and far away from the asymptotic steady state. In both cases, good agreement is found again between theory and DSMC results. As a complement to the previous studies, we determine in this paper the dependence of the steady values of the dynamic properties of the suspension on the parameter space of the system. More specifically, we focus our attention on the temperature ratio T1/T2 and the fourth degree cumulants c1 and c2 (measuring the departure of the velocity distributions f1 and f2 from their Maxwellian forms). While our approximate theoretical expression for T1/T2 agrees very well with computer simulations, some discrepancies are found for the cumulants. Finally, a linear stability analysis of the steady state solution is also carried out showing that the steady state is always linearly stable.

The diffusion transport coefficients of a binary granular suspension where one of the components is present in tracer concentration are determined from the (inelastic) Enskog kinetic equation. The effect of the interstitial gas on the solid particles is accounted for in the kinetic equation through two different terms: (i) a viscous drag force proportional to the particle velocity and (ii) stochastic Langevin-like term defined in terms of the background temperature. The transport coefficients are obtained as the solutions of a set of coupled linear integral equations recently derived for binary granular suspensions with arbitrary concentration [Gómez González et al., “Enskog kinetic theory for multicomponent granular suspensions,” Phys. Rev. E 101, 012904 (2020)]. To achieve analytical expressions for the diffusion coefficients, which can be sufficiently accurate for highly inelastic collisions and/or disparate values of the mass and diameter rations, the above integral equations are approximately solved by considering the so-called second Sonine approximation (two terms in the Sonine polynomial expansion of the distribution function). The theoretical results for the tracer diffusion coefficient D0 (coefficient connecting the mass flux with the gradient of density of tracer particles) are compared with those obtained by numerically solving the Enskog equation by means of the direct simulation Monte Carlo method. Although the first-Sonine approximation to D0 yields, in general, a good agreement with simulation results, we show that the second-Sonine approximation leads to an improvement over the first-Sonine correction, especially when the tracer particles are much lighter than the granular gas. The expressions derived here for the diffusion coefficients are also used for two different applications. First, the stability of the homogeneous steady state is discussed. Second, segregation induced by a thermal gradient is studied. As expected, the results show that the corresponding phase diagrams for segregation clearly differ from those found in previous works when the effect of gas phase on grains is neglected.