No Arabic abstract
Transport coefficients associated with the mass flux of impurities immersed in a moderately dense granular gas of hard disks or spheres described by the inelastic Enskog equation are obtained by means of the Chapman-Enskog expansion. The transport coefficients are determined as the solutions of a set of coupled linear integral equations recently derived for polydisperse granular mixtures [V. Garzo, J. W. Dufty and C. M. Hrenya, Phys. Rev. E {bf 76}, 031304 (2007)]. With the objective of obtaining theoretical expressions for the transport coefficients that are sufficiently accurate for highly inelastic collisions, we solve the above integral equations by using the second Sonine approximation. As a complementary route, we numerically solve by means of the direct simulation Monte Carlo method (DSMC) the inelastic Enskog equation to get the kinetic diffusion coefficient $D_0$ for two and three dimensions. We have observed in all our simulations that the disagreement, for arbitrarily large inelasticity, in the values of both solutions (DSMC and second Sonine approximation) is less than 4%. Moreover, we show that the second Sonine approximation to $D_0$ yields a dramatic improvement (up to 50%) over the first Sonine approximation for impurity particles lighter than the surrounding gas and in the range of large inelasticity. The results reported in this paper are of direct application in important problems in granular flows, such as segregation driven by gravity and a thermal gradient. We analyze here the segregation criteria that result from our theoretical expressions of the transport coefficients.
The Einstein relation for a driven moderately dense granular gas in $d$-dimensions is analyzed in the context of the Enskog kinetic equation. The Enskog equation neglects velocity correlations but retains spatial correlations arising from volume exclusion effects. As expected, there is a breakdown of the Einstein relation $epsilon=D/(T_0mu) eq 1$ relating diffusion $D$ and mobility $mu$, $T_0$ being the temperature of the impurity. The kinetic theory results also show that the violation of the Einstein relation is only due to the strong non-Maxwellian behavior of the reference state of the impurity particles. The deviation of $epsilon$ from unity becomes more significant as the solid volume fraction and the inelasticity increase, especially when the system is driven by the action of a Gaussian thermostat. This conclusion qualitatively agrees with some recent simulations of dense gases [Puglisi {em et al.}, 2007 {em J. Stat. Mech.} P08016], although the deviations observed in computer simulations are more important than those obtained here from the Enskog kinetic theory. Possible reasons for the quantitative discrepancies between theory and simulations are discussed.
A solution of the inelastic Enskog equation that goes beyond the weak dissipation limit and applies for moderate densities is used to determine the thermal diffusion factor of an intruder immersed in a dense granular gas under gravity. This 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, density and coefficients of restitution). The form of the phase-diagrams for the BNE/RBNE transition depends sensitively on the value of gravity relative to the thermal gradient, so that it is possible to switch between both states for given values of the parameters of the system. Two specific limits are considered with detail: (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 extends previous theoretical results derived in the dilute limit case [V. Garzo, Europhys. Lett. {bf 75}, 521 (2006)] and is consistent with the findings of some recent experimental results.
A new segregation criterion based on the inelastic Enskog kinetic equation is derived to show the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE) by varying the different parameters of the system. In contrast to previous theoretical attempts the approach is not limited to the near-elastic case, takes into account the influence of both thermal gradients and gravity and applies for moderate densities. The form of the phase-diagrams for the BNE/RBNE transition depends sensitively on the value of gravity relative to the thermal gradient, so that it is possible to switch between both states for given values of the mass and size ratios, the coefficients of restitution and the solid volume fraction. In particular, the influence of collisional dissipation on segregation becomes more important when the thermal gradient dominates over gravity than in the opposite limit. The present analysis extends previous results derived in the dilute limit case and is consistent with the findings of some recent experimental results.
A granular gas may be modeled as a set of hard-spheres undergoing inelastic collisions; its microscopic dynamics is thus strongly irreversible. As pointed out in several experimental works bearing on turbulent flows or granular materials, the power injected in a dissipative system to sustain a steady-state over an asymptotically large time window is a central observable. We describe an analytic approach allowing us to determine the full distribution of the power injected in a granular gas within a steady-state resulting from subjecting each particle independently either to a random force (stochastic thermostat) or to a deterministic force proportional to its velocity (Gaussian thermostat). We provide an analysis of our results in the light of the relevance, for other types of systems, of the injected power to fluctuation relations.
We consider a dilute gas of inelastic hard spheres enclosed in a slab under the action of gravity along the longitudinal direction. In addition, the gas is subject to a white-noise stochastic force that mimics the effect of external vibrations customarily used in experiments to compensate for the collisional cooling. The system is described by means of a kinetic model of the inelastic Boltzmann equation and its steady-state solution is derived through second order in gravity. This solution differs from the Navier-Stokes description in that the hydrostatic pressure is not uniform, normal stress differences are present, a component of the heat flux normal to the thermal gradient exists, and the temperature profile includes a positive quadratic term. As in the elastic case, this new term is responsible for a bimodal shape of the temperature profile. The results show that, except for high inelasticities, the effect of inelasticity on the profiles is to slightly decrease the quantitative deviations from the Navier-Stokes results.