Assessing a Hydrodynamic Description for Instabilities in Highly Dissipative, Freely Cooling Granular Gases

An intriguing phenomenon displayed by granular flows and predicted by kinetic-theory-based models is the instability known as particle clustering, which refers to the tendency of dissipative grains to form transient, loose regions of relatively high concentration. In this work, we assess a modified-Sonine approximation recently proposed [Garzo et al., Physica A 376, 94 (2007)] for a granular gas via an examination of system stability. In particular, we determine the critical length scale associated with the onset of two types of instabilities -vortices and clusters- via stability analyses of the Navier-Stokes-order hydrodynamic equations by using the expressions of the transport coefficients obtained from both the standard and the modified-Sonine approximations. We examine the impact of both Sonine approximations over a range of solids fraction phi <0.2 for small restitution coefficients e=0.25--0.4, where the standard and modified theories exhibit discrepancies. The theoretical predictions for the critical length scales are compared to molecular dynamics (MD) simulations, of which a small percentage were not considered due to inelastic collapse. Results show excellent quantitative agreement between MD and the modified-Sonine theory, while the standard theory loses accuracy for this highly dissipative parameter space. The modified theory also remedies a (highdissipation) qualitative mismatch between the standard theory and MD for the instability that forms more readily. Furthermore, the evolution of cluster size is briefly examined via MD, indicating that domain-size clusters may remain stable or halve in size, depending on system parameters.

Enskog Theory for Polydisperse Granular Mixtures. III. Comparison of dense and dilute transport coefficients and equations of state for a binary mixture

The objective of this study is to assess the impact of a dense-phase treatment on the hydrodynamic description of granular, binary mixtures relative to a previous dilute-phase treatment. Two theories were considered for this purpose. The first, proposed by Garzo and Dufty (GD) [Phys. Fluids {bf 14}, 146 (2002)], is based on the Boltzmann equation which does not incorporate finite-volume effects, thereby limiting its use to dilute flows. The second, proposed by Garzo, Hrenya and Dufty (GHD) [Phys. Rev. E {bf 76}, 31303 and 031304 (2007)], is derived from the Enskog equation which does account for finite-volume effects; accordingly this theory can be applied to moderately dense systems as well. To demonstrate the significance of the dense-phase treatment relative to its dilute counterpart, the ratio of dense (GHD) to dilute (GD) predictions of all relevant transport coefficients and equations of state are plotted over a range of physical parameters (volume fraction, coefficients of restitution, material density ratio, diameter ratio, and mixture composition). These plots reveal the deviation between the two treatments, which can become quite large ($>$100%) even at moderate values of the physical parameters. Such information will be useful when choosing which theory is most applicable to a given situation, since the dilute theory offers relative simplicity and the dense theory offers improved accuracy. It is also important to note that several corrections to original GHD expressions are presented here in the form of a complete, self-contained set of relevant equations.

Mass transport of impurities in a moderately dense granular gas

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.

Segregation by thermal diffusion of an intruder in a moderately dense granular fluid

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.

Brazil-nut effect versus reverse Brazil-nut effect in a moderately dense granular fluid

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 note on the violation of the Einstein relation in a driven moderately dense granular gas

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.

Mass transport in a strongly sheared binary mixture of Maxwell molecules

Transport coefficients associated with the mass flux of a binary mixture of Maxwell molecules under uniform shear flow are exactly determined from the Boltzmann kinetic equation. A normal solution is obtained via a Chapman--Enskog-like expansion around a local shear flow distribution that retains all the hydrodynamics orders in the shear rate. In the first order of the expansion the mass flux is proportional to the gradients of mole fraction, pressure, and temperature but, due to the anisotropy induced in the system by the shear flow, mutual diffusion, pressure diffusion and thermal diffusion tensors are identified instead of the conventional scalar coefficients. These tensors are obtained in terms of the shear rate and the parameters of the mixture (particle masses, concentrations, and force constants). The description is made both in the absence and in the presence of an external thermostat introduced in computer simulations to compensate for the viscous heating. As expected, the analysis shows that there is not a simple relationship between the results with and without the thermostat. The dependence of the three diffusion tensors on the shear rate is illustrated in the tracer limit case, the results showing that the deviation of the generalized transport coefficients from their equilibrium forms is in general quite important. Finally, the generalized transport coefficients associated with the momentum and heat transport are evaluated from a model kinetic equation of the Boltzmann equation.

Shear-rate dependent transport coefficients for inelastic Maxwell models

The Boltzmann equation for d-dimensional inelastic Maxwell models is considered to analyze transport properties in spatially inhomogeneous states close to the simple shear flow. A normal solution is obtained via a Chapman--Enskog--like expansion around a local shear flow distribution f^{(0)} that retains all the hydrodynamic orders in the shear rate. The constitutive equations for the heat and momentum fluxes are obtained to first order in the deviations of the hydrodynamic field gradients from their values in the reference state and the corresponding generalized transport coefficients are {em exactly} determined in terms of the coefficient of restitution alpha and the shear rate a. Since f^{(0)} applies for arbitrary values of the shear rate and is not restricted to weak dissipation, the transport coefficients turn out to be nonlinear functions of both parameters a and alpha. A comparison with previous results obtained for inelastic hard spheres from a kinetic model of the Boltzmann equation is also carried out.

Kinetic Theory for Binary Granular Mixtures at Low-Density

Many features of granular media can be modelled as a fluid of hard spheres with {em inelastic} collisions. Under rapid flow conditions, the macroscopic behavior of grains can be described through hydrodynamic equations. At low-density, a fundamental basis for the derivation of the hydrodynamic equations and explicit expressions for the transport coefficients appearing in them is provided by the Boltzmann kinetic theory conveniently modified to account for inelastic binary collisions. The goal of this chapter is to give an overview of the recent advances made for binary granular gases by using kinetic theory tools. Some of the results presented here cover aspects such as transport properties, energy nonequipartition, instabilities, segregation or mixing, non-Newtonian behavior, .... In addition, comparison of the analytical results with those obtained from Monte Carlo and molecular dynamics simulations is also carried out, showing the reliability of kinetic theory to describe granular flows even for strong dissipation.