Hydrodynamics and transport coefficients for dilute granular gases

Translational and rotational non-Gaussianities in homogeneous freely evolving granular gases

The importance of roughness in the modeling of granular gases has been increasingly considered in recent years. In this paper, a freely evolving homogeneous granular gas of inelastic and rough hard disks or spheres is studied under the assumptions of the Boltzmann kinetic equation. The homogeneous cooling state is studied from a theoretical point of view using a Sonine approximation, in contrast to a previous Maxwellian approach. A general theoretical description is done in terms of d_{t} translational and d_{r} rotational degrees of freedom, which accounts for the cases of spheres (d_{t}=d_{r}=3) and disks (d_{t}=2, d_{r}=1) within a unified framework. The non-Gaussianities of the velocity distribution function of this state are determined by means of the first nontrivial cumulants and by the derivation of non-Maxwellian high-velocity tails. The results are validated by computer simulations using direct simulation Monte Carlo and event-driven molecular dynamics algorithms.

Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. I. Transport coefficients

The transport coefficients for dilute granular gases of inelastic and rough hard disks or spheres with constant coefficients of normal (α) and tangential (β) restitution are obtained in a unified framework as functions of the number of translational (d_{t}) and rotational (d_{r}) degrees of freedom. The derivation is carried out by means of the Chapman-Enskog method with a Sonine-like approximation in which, in contrast to previous approaches, the reference distribution function for angular velocities does not need to be specified. The well-known case of purely smooth d-dimensional particles is recovered by setting d_{t}=d and formally taking the limit d_{r}→0. In addition, previous results [G. M. Kremer, A. Santos, and V. Garzó, Phys. Rev. E 90, 022205 (2014)10.1103/PhysRevE.90.022205] for hard spheres are reobtained by taking d_{t}=d_{r}=3, while novel results for hard-disk gases are derived with the choice d_{t}=2, d_{r}=1. The singular quasismooth limit (β→-1) and the conservative Pidduck's gas (α=β=1) are also obtained and discussed.

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

Kinetic theory of shear thickening for a moderately dense gas-solid suspension: From discontinuous thickening to continuous thickening

The Enskog kinetic theory for moderately dense gas-solid suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the environmental fluid on solid particles is modeled via a viscous drag force plus a stochastic Langevin-like term. The Enskog equation is solved by means of two independent but complementary routes: (i) Grad's moment method and (ii) event-driven Langevin simulation of hard spheres. Both approaches clearly show that the flow curve (stress-strain rate relation) depends significantly on the volume fraction of the solid particles. In particular, as the density increases, there is a transition from the discontinuous shear thickening (observed in dilute gases) to the continuous shear thickening for denser systems. The comparison between theory and simulations indicates that while the theoretical predictions for the kinetic temperature agree well with simulations for densities φ≲0.5, the agreement for the other rheological quantities (the viscosity, the stress ratio, and the normal stress differences) is limited to more moderate densities (φ≲0.3) if the inelasticity during collisions between particles is not large.

Velocity Distribution of a Homogeneously Driven Two-Dimensional Granular Gas

The theory of homogeneously driven granular gases of hard particles predicts that the stationary state is characterized by a velocity distribution function with overpopulated high-energy tails as compared to the exponential decay valid for molecular gases. While this fundamental theoretical result was confirmed by numerous numerical simulations, an experimental confirmation is still missing. Using self-rotating active granular particles, we find a power-law decay of the velocity distribution whose exponent agrees well with the theoretic prediction.

Structure of velocity distributions in shock waves in granular gases with extension to molecular gases

Velocity distributions in normal shock waves obtained in dilute granular flows are studied. These distributions cannot be described by a simple functional shape and are believed to be bimodal. Our results show that these distributions are not strictly bimodal but a trimodal distribution is shown to be sufficient. The usual Mott-Smith bimodal description of these distributions, developed for molecular gases, and based on the coexistence of two subpopulations (a supersonic and a subsonic population) in the shock front, can be modified by adding a third subpopulation. Our experiments show that this additional population results from collisions between the supersonic and subsonic subpopulations. We propose a simple approach incorporating the role of this third intermediate population to model the measured probability distributions and apply it to granular shocks as well as shocks in molecular gases.

Kinetic theory for dilute cohesive granular gases with a square well potential

We develop the kinetic theory of dilute cohesive granular gases in which the attractive part is described by a square well potential. We derive the hydrodynamic equations from the kinetic theory with the microscopic expressions for the dissipation rate and the transport coefficients. We check the validity of our theory by performing the direct simulation Monte Carlo.

Declustering in a granular gas as a finite-size effect

The existence of dense clusters has been shown to be a transient phenomenon for realistic models of granular collisions, where the coefficient of restitution depends on the impact velocity. We report direct numerical simulations that elucidate the conditions for the disappearance of structures. We find that upon cluster formation the granular temperature and the convective kinetic energy couple and both follow Haff's law. Furthermore, we show that clusters will eventually dissolve in all finite-size systems. We find the strong power law t'∝L(12) for the dependency of the declustering time on system size. Our results imply that only in systems close to the initial critical system size both clustering and declustering transitions are observable.

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.

Stochastic behavior of the coefficient of normal restitution

We consider the collision of a rough sphere with a plane by detailed analysis of the collision geometry. Using stochastic methods, the effective coefficient of restitution may be described as a fluctuating quantity whose probability density follows an asymmetric Laplace distribution. This result agrees with recent experiments by Montaine et al. [Phys. Rev. E 84, 041306 (2011)].

Long-time tails in freely cooling granular gases

The long-time behavior of the current autocorrelation functions for the velocity, the shear stress, and the heat flux is investigated in freely cooling granular gases. It is found that the correlation functions for the velocity and the shear stress have the long-time tails obeying tau(-d/2), while the correlation function for heat flux decays as tau(-(d+2)/2) exp(-zeta*tau) with the dimensionless cooling rate zeta*, the spatial dimension d, and the scaled time tau in terms of the collision frequency. The result of our numerical simulation of the freely cooling granular gases is consistent with the theoretical prediction.

Self-diffusion in granular gases: Green-Kubo versus Chapman-Enskog

We study the diffusion of tracers (self-diffusion) in a homogeneously cooling gas of dissipative particles, using the Green-Kubo relation and the Chapman-Enskog approach. The dissipative particle collisions are described by the coefficient of restitution epsilon which for realistic material properties depends on the impact velocity. First, we consider self-diffusion using a constant coefficient of restitution, epsilon=const, as frequently used to simplify the analysis. Second, self-diffusion is studied for a simplified (stepwise) dependence of epsilon on the impact velocity. Finally, diffusion is considered for gases of realistic viscoelastic particles. We find that for epsilon=const both methods lead to the same result for the self-diffusion coefficient. For the case of impact-velocity dependent coefficients of restitution, the Green-Kubo method is, however, either restrictive or too complicated for practical application, therefore we compute the diffusion coefficient using the Chapman-Enskog method. We conclude that in application to granular gases, the Chapman-Enskog approach is preferable for deriving kinetic coefficients.

Clustering of atoms in a model with multiple thermostats

We propose a model for a one-dimensional chain of interacting particles in an external periodic potential. In this model the particles have a complex structure treated in a mean-field fashion: particle collisions are inelastic and also each particle is considered as having its own thermostat. We derived the Fokker-Planck equation for this model and demonstrated that the model has a truly equilibrium ground state. When an external dc force is applied to the atoms, the model exhibits a hysteresis even at high temperatures due to the clustering of atoms with the same velocity. Another effect of clustering is phase separation in the steady state when the system splits into regions of immobile atoms ("traffic jams") and regions of running atoms.

Transient structures in a granular gas

A force-free granular gas is considered with an impact-velocity-dependent coefficient of restitution as it follows from the model of viscoelastic particles. We analyze structure formation in this system by means of three independent methods: molecular dynamics, numerical solution of the hydrodynamic equations, and linear stability analysis of these equations. All these approaches indicate that structure formation occurs in force-free granular gases only as a transient process.