Velocity distribution in granular gases of viscoelastic particles

Granular gases under resetting

We investigate the granular temperatures in force-free granular gases under exponential resetting. When a resetting event occurs, the granular temperature attains its initial value, whereas it decreases because of the inelastic collisions between the resetting events. We develop a theory and perform computer simulations for granular gas cooling in the presence of Poissonian resetting events. We also investigate the probability density function to quantify the distribution of granular temperatures. Our theory may help us to understand the behavior of nonperiodically driven granular systems.

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

Temperature distribution in driven granular mixtures does not depend on mechanism of energy dissipation

We study analytically and numerically the distribution of granular temperatures in granular mixtures for different dissipation mechanisms of inelastic inter-particle collisions. Both driven and force-free systems are analyzed. We demonstrate that the simplified model of a constant restitution coefficient fails to predict even qualitatively a granular temperature distribution in a homogeneous cooling state. At the same time we reveal for driven systems a stunning result - the distribution of temperatures in granular mixtures is universal. That is, it does not depend on a particular dissipation mechanism of inter-particles collisions, provided the size distributions of particles is steep enough. The results of the analytic theory are compared with simulation results obtained by the direct simulation Monte Carlo (DSMC). The agreement between the theory and simulations is perfect. The reported results may have important consequences for fundamental science as well as for numerous application, e.g. for the experimental modelling in a lab of natural processes.

Electrification in granular gases leads to constrained fractal growth

The empirical observation of aggregation of dielectric particles under the influence of electrostatic forces lies at the origin of the theory of electricity. The growth of clusters formed of small grains underpins a range of phenomena from the early stages of planetesimal formation to aerosols. However, the collective effects of Coulomb forces on the nonequilibrium dynamics and aggregation process in a granular gas - a model representative of the above physical processes - have so far evaded theoretical scrutiny. Here, we establish a hydrodynamic description of aggregating granular gases that exchange charges upon collisions and interact via the long-ranged Coulomb forces. We analytically derive the governing equations for the evolution of granular temperature, charge variance, and number density for homogeneous and quasi-monodisperse aggregation. We find that, once the aggregates are formed, the granular temperature of the cluster population, the charge variance of the cluster population and the number density of the cluster population evolve in such a way that their non-dimensional combination obeys a physical constraint of nearly constant dimensionless ratio of characteristic electrostatic to kinetic energy. This constraint on the collective evolution of charged clusters is confirmed both by  our theory and our detailed molecular dynamics simulations. The inhomogeneous aggregation of monomers and clusters in their mutual electrostatic field proceeds in a fractal manner. Our theoretical framework is extendable to more precise charge exchange mechanisms, a current focus of extensive experimentation. Furthermore, it illustrates the collective role of long-ranged interactions in dissipative gases and can lead to novel designing principles in particulate systems.

Increasing temperature of cooling granular gases

The kinetic energy of a force-free granular gas decays monotonously due to inelastic collisions of the particles. For a homogeneous granular gas of identical particles, the corresponding decay of granular temperature is quantified by Haff’s law. Here, we report that for a granular gas of aggregating particles, the granular temperature does not necessarily decay but may even increase. Surprisingly, the increase of temperature is accompanied by the continuous loss of total gas energy. This stunning effect arises from a subtle interplay between decaying kinetic energy and gradual reduction of the number of degrees of freedom associated with the particles’ dynamics. We derive a set of kinetic equations of Smoluchowski type for the concentrations of aggregates of different sizes and their energies. We find scaling solutions to these equations and a condition for the aggregation mechanism predicting growth of temperature. Numerical direct simulation Monte Carlo results confirm the theoretical predictions.

Granular gases—dilute systems composed of dissipatively colliding particles—exhibit anomalous dynamics and numerous surprising phenomena. Here, Brilliantov et al. show that the aggregation mechanism can induce increase of the gas temperature despite the fact that the total kinetic energy decreases.

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.

Granular rotor as a probe for a nonequilibrium bath

This study numerically and analytically investigates the dynamics of a rotor under viscous or dry friction as a nonequilibrium probe of a granular gas. In order to demonstrate the role of the rotor as a probe for a nonequilibrium bath, the molecular dynamics (MD) simulation of the rotor is performed under viscous or dry friction surrounded by a steady granular gas under gravity. A one-to-one map between the velocity distribution function (VDF) of the granular gas and the angular distribution function for the rotor is theoretically derived. The MD simulation demonstrates that the one-to-one map accurately infers the local VDF of the granular gas from the angular VDF of the rotor, and vice versa.

Quantifying non-ergodic dynamics of force-free granular gases

We demonstrate how non-ergodicity arises in simple mechanistic systems such as force free, dissipative granular gases. This behaviour results from the strong non-stationarity of the process mirrored in the continuous decay of the gas temperature.

Collision dynamics of particle clusters in a two-dimensional granular gas

In a granular gas, inelastic collisions produce an instability in which the constituent particles cluster heterogeneously. These clusters then interact with each other, further decreasing their kinetic energy. We report experiments of the free collisions of dense clusters of particles in a two-dimensional geometry. The particles are composed of solid CO(2), which float nearly frictionlessly on a hot surface due to sublimated vapor. After two dense clusters of ≈100 particles collide, there are two distinct stages of evolution. First, the translational kinetic energy rapidly decreases by over 90% as a "jamming front" sweeps across each cluster. Subsequently, the kinetic energy decreases more slowly as the particles approach the container boundaries. In this regime, the measured velocity distributions are non-Gaussian with long tails. Finally, we compare our experiments to computer simulations of colliding, two-dimensional, granular clusters composed of circular, viscoelastic particles with friction.

Velocity distribution function and effective restitution coefficient for a granular gas of viscoelastic particles

We perform large-scale event-driven molecular dynamics (MD) simulations for granular gases of particles interacting with the impact-velocity-dependent restitution coefficient ε(v(imp)). We use ε(v(imp)) as it follows from the simplest first-principles collision model of viscoelastic spheres. Both cases of force-free and uniformly heated gases are studied. We formulate a simplified model of an effective constant restitution coefficient ε(eff), which depends on a current granular temperature, and we compute ε(eff) using the kinetic theory. We develop a theory of the velocity distribution function for driven gases of viscoelastic particles and analyze the evolution of granular temperature and of the Sonine coefficients, which characterize the form of the velocity distribution function. We observe that for not large dissipation the simulation results are in an excellent agreement with the theory for both the homogeneous cooling state and uniformly heated gases. At the same time, a noticeable discrepancy between the theory and MD results for the Sonine coefficients is detected for large dissipation. We analyze the accuracy of the simplified model based on the effective restitution coefficient ε(eff), and we conclude that this model can accurately describe granular temperature. It provides also an acceptable accuracy for the velocity distribution function for small dissipation, but it fails when dissipation is large.

Intermediate regimes in granular Brownian motion: superdiffusion and subdiffusion

Brownian motion in a granular gas in a homogeneous cooling state is studied theoretically and by means of molecular dynamics. We use the simplest first-principles model for the impact-velocity dependent restitution coefficient, as it follows for the model of viscoelastic spheres. We reveal that for a wide range of initial conditions the ratio of granular temperatures of Brownian and bath particles demonstrates complicated nonmonotonic behavior, which results in a transition between different regimes of Brownian dynamics: It starts from the ballistic motion, switches later to a superballistic one, and turns at still later times into subdiffusion; eventually normal diffusion is achieved. Our theory agrees very well with the molecular dynamics results, although extreme computational costs prevented us from detecting the final diffusion regime. Qualitatively, the reported intermediate diffusion regimes are generic for granular gases with any realistic dependence of the restitution coefficient on the impact velocity.

Diffusivity and weak clustering in a quasi-two-dimensional granular gas

We present results from a detailed simulation of a quasi-two-dimensional dissipative granular gas, kept in a noncondensed steady state via vertical shaking over a rough substrate. This gas shows a weak power-law decay in the tails of its pair distribution functions, indicating clustering. This clustering depends monotonically on the dissipation coefficient and disappears when the sphere-sphere collisions are conservative. Clustering is also sensitive to the packing fraction. This gas also displays the standard nonequilibrium characteristics of similar systems, including non-Maxwellian velocity distributions. The diffusion coefficients are calculated over all the conditions of the simulations, and it is found that diluted gases are more diffusive for smaller restitution coefficients.

Equivalence of the freely cooling granular gas to the sticky gas

A freely cooling granular gas with a velocity-dependent restitution coefficient is studied in one dimension. The restitution coefficient becomes near elastic when the relative velocity of the colliding particles is less than a velocity scale delta . Different statistical quantities, namely, density distribution, occupied and empty cluster length distributions, and spatial density and velocity correlation functions, are obtained using event driven molecular dynamic simulations. We compare these with the corresponding quantities of the sticky gas (inelastic gas with zero coefficient of restitution). We find that in the inhomogeneous cooling regime, for times smaller than a crossover time t{1} , where t{1} approximately delta;{-1} , the behavior of the granular gas is equivalent to that of the sticky gas. When delta-->0 , then t{1}-->infinity and, hence, the results support an earlier claim that the freely cooling inelastic gas is described by the inviscid Burgers equation. For a real granular gas with finite delta , the existence of the time scale t{1} shows that, for large times, the granular gas is not described by the inviscid Burgers equation.

Impact of high-energy tails on granular gas properties

The velocity distribution function of granular gases in the homogeneous cooling state as well as some heated granular gases decays for large velocities as f proportional to exp(-const x v). That is, its high-energy tail is overpopulated as compared with the Maxwell distribution. At the present time, there is no theory to describe the influence of the tail on the kinetic characteristics of granular gases. We develop an approach to quantify the overpopulated tail and analyze its impact on granular gas properties, in particular on the cooling coefficient. We observe and explain anomalously slow relaxation of the velocity distribution function to its steady state.

Granular clustering: self-consistent analysis for general coefficients of restitution

We study the equilibrium behavior of one-dimensional granular clusters and one-particle granular gases for a variety of velocity-dependent coefficients of restitution r. We obtain equations describing the long-time behavior for the cluster's pressure, rms velocity, and granular interspacing. We show that for extremely long times, clusters with velocity-dependent coefficients of restitution are unstable and dissolve into homogeneous, quasielastic gases, but clusters with velocity-independent r are permanent. This is in accordance with hydrodynamic studies pointing to the transient nature of density instabilities for granular gases with velocity-dependent r.

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.

Simulations of vibrated granular medium with impact-velocity-dependent restitution coefficient

We report numerical simulations of strongly vibrated granular materials designed to mimic recent experiments performed in both the presence and the absence of gravity. The coefficient of restitution used here depends on the impact velocity by taking into account both the viscoelastic and plastic deformations of particles, occurring at low and high velocities, respectively. We show that this model with impact-velocity-dependent restitution coefficient reproduces results that agree with experiments. We measure the scaling exponents of the granular temperature, collision frequency, impulse, and pressure with the vibrating piston velocity as the particle number increases. As the system changes from a homogeneous gas state at low density to a clustered state at high density, these exponents are all found to decrease continuously with increasing particle number. All these results differ significantly from classical inelastic hard sphere kinetic theory and previous simulations, both based on a constant restitution coefficient.

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.

Hydrodynamics and transport coefficients for dilute granular gases

The hydrodynamics of granular gases of viscoelastic particles, whose collision is described by an impact-velocity dependent coefficient of restitution, which is developed using the Chapman-Enskog approach. We derive the hydrodynamic equations and the according transport coefficients with the assumption that the shape of the velocity distribution function follows adiabatically the decaying temperature. We show numerically that this approximation is justified up to intermediate dissipation. The transport coefficients and the coefficient of cooling are expressed in terms of the elastic and dissipative parameters of the particle material and by the gas parameters. The dependence of these coefficients on temperature differs qualitatively from that obtained with the simplifying assumption of a constant coefficient of restitution which was used in previous studies. The approach formulated for gases of viscoelastic particles may be applied also for other impact-velocity dependencies of the restitution coefficient.

Hydrodynamics of granular gases of viscoelastic particles

Our study examines the long-time behaviour of a force-free granular gas of viscoelastic particles, for which the coefficient of restitution depends on the impact velocity, as it follows from the solution of the impact problem for viscoelastic spheres. Starting from the Boltzmann equation, we derived the hydrodynamic equations and obtained microscopic expressions for the transport coefficients in terms of the elastic and dissipative parameters of the particle material. We performed the stability analysis of the linearized set of equations and found that any inhomogeneities and vortices vanish after a long time and the system approaches the flow-free stage of homogeneous density. This behaviour is in contrast to that of a gas consisting of particles which interact via a (non-realistic) constant coefficient of restitution, for which inhomogeneities (clusters) and vortex patterns have been proven to arise and to continuously develop.