Granular fluid thermostated by a bath of elastic hard spheres

Rheology of granular particles immersed in a molecular gas under uniform shear flow

Non-Newtonian transport properties of a dilute gas of inelastic hard spheres immersed in a molecular gas are determined. We assume that the granular gas is sufficiently rarefied, and hence the state of the molecular gas is not disturbed by the presence of the solid particles. In this situation, one can treat the molecular gas as a bath (or thermostat) of elastic hard spheres at a given temperature. Moreover, in spite of the fact that the number density of grains is quite small, we take into account their inelastic collisions among themselves in its kinetic equation. The system (granular gas plus a bath of elastic hard spheres) is subjected to a simple (or uniform) shear flow. In the low-density regime, the rheological properties of the granular gas are determined by solving the Boltzmann kinetic equation by means of Grad's moment method. These properties turn out to be highly nonlinear functions of the shear rate and the remaining parameters of the system. Our results show that the kinetic granular temperature and the non-Newtonian viscosity present a discontinuous shear thickening effect for sufficiently high values of the mass ratio m/m_{g} (m and m_{g} being the mass of grains and gas particles, respectively). This effect becomes more pronounced as the mass ratio m/m_{g} increases. In particular, in the Brownian limit (m/m_{g}→∞) the expressions of the non-Newtonian transport properties derived here are consistent with those previously obtained by considering a coarse-grained approach where the effect of gas phase on grains is through an effective force. Theoretical results are compared against computer simulations, showing an excellent agreement.

Exact Results for Non-Newtonian Transport Properties in Sheared Granular Suspensions: Inelastic Maxwell Models and BGK-Type Kinetic Model

The Boltzmann kinetic equation for dilute granular suspensions under simple (or uniform) shear flow (USF) is considered to determine the non-Newtonian transport properties of the system. In contrast to previous attempts based on a coarse-grained description, our suspension model accounts for the real collisions between grains and particles of the surrounding molecular gas. The latter is modeled as a bath (or thermostat) of elastic hard spheres at a given temperature. Two independent but complementary approaches are followed to reach exact expressions for the rheological properties. First, the Boltzmann equation for the so-called inelastic Maxwell models (IMM) is considered. The fact that the collision rate of IMM is independent of the relative velocity of the colliding spheres allows us to exactly compute the collisional moments of the Boltzmann operator without the knowledge of the distribution function. Thanks to this property, the transport properties of the sheared granular suspension can be exactly determined. As a second approach, a Bhatnagar-Gross-Krook (BGK)-type kinetic model adapted to granular suspensions is solved to compute the velocity moments and the velocity distribution function of the system. The theoretical results (which are given in terms of the coefficient of restitution, the reduced shear rate, the reduced background temperature, and the diameter and mass ratios) show, in general, a good agreement with the approximate analytical results derived for inelastic hard spheres (IHS) by means of Grad's moment method and with computer simulations performed in the Brownian limiting case (m/mg→∞, where mg and m are the masses of the particles of the molecular and granular gases, respectively). In addition, as expected, the IMM and BGK results show that the temperature and non-Newtonian viscosity exhibit an S shape in a plane of stress-strain rate (discontinuous shear thickening, DST). The DST effect becomes more pronounced as the mass ratio m/mg increases.

Kovacs-like memory effect in driven granular gases

While memory effects have been reported for dense enough disordered systems such as glasses, we show here by a combination of analytical and simulation techniques that they are also intrinsic to the dynamics of dilute granular gases. By means of a certain driving protocol, we prepare the gas in a state where the granular temperature T coincides with its long time limit. However, T does not subsequently remain constant but exhibits a nonmonotonic evolution before reaching its nonequilibrium steady value. The corresponding so-called Kovacs hump displays a normal behavior for weak dissipation (as observed in molecular systems) but is reversed under strong dissipation, where it, thus, becomes anomalous.

Velocity-correlation distributions in thermal and granular hard-core gases

Collision statistics of hard-core systems (thermal and dissipative) is investigated through the velocity-correlation distributions after n collisions of a tagged hard-core particle: These quantities provide information on the velocity correlations for a given number of collisions. We obtain exact results for arbitrary dimension for the velocity-correlation distribution after the first collision as well as for the velocity-correlation function after an infinite number of collisions. For Gaussian velocity distributions, we show that the decay of the first-collision velocity-correlation distribution for negative argument is always exponential in any dimension; the decay rate is then a function of the mass and the coefficient of restitution. For granular gases, where deviations from Gaussian are relevant, expressions including Sonine corrections are also derived for the velocity-correlation distribution and a comparison with a direct simulation Monte Carlo (DSMC) shows accurate agreement with theoretical results. We emphasize that these quantities can be easily obtained in simulations and most likely also in experiments: therefore they could be an efficient probe of the local environment and of the degree of inelasticity of the collisions.

Angular velocity distribution of a granular planar rotator in a thermalized bath

The kinetics of a granular planar rotator with a fixed center undergoing inelastic collisions with bath particles is analyzed both numerically and analytically by means of the Boltzmann equation. The angular velocity distribution evolves from quasi-Gaussian in the Brownian limit to an algebraic decay in the limit of an infinitely light particle. In addition, we compare this model to that of a planar rotator with a free center and discuss the prospects for experimental confirmation of these results.

Thermalization of an anisotropic granular particle

We investigate the dynamics of a needle in a two-dimensional bath composed of thermalized point particles. Collisions between the needle and points are inelastic and characterized by a normal restitution coefficient alpha<1. By using the Enskog-Boltzmann equation, we obtain analytical expressions for the translational and rotational granular temperatures of the needle and show that these are, in general, different from the bath temperature. The translational temperature always exceeds the rotational one, though the difference decreases with increasing moment of inertia. The predictions of the theory are in very good agreement with numerical simulations of the model.

Hydrodynamics of driven granular gases

Hydrodynamic equations for granular gases driven by the Fokker-Planck operator are derived. Transport coefficients appeared in Navier-Stokes order change from the values of a free-cooling state to those of a steady state. The mismatch between the granular temperature and the constant temperature of the heat bath produces a nontrivial steady state.