Energy decay in three-dimensional freely cooling granular gas

Lattice models for ballistic aggregation: Cluster-shape-dependent exponents

We study ballistic aggregation on a two-dimensional square lattice, where particles move ballistically in between momentum and mass conserving coalescing collisions. Three models are studied based on the shapes of the aggregates: In the first the aggregates remain point particles, in the second they retain the fractal shape at the time of collision, and in the third they assume a spherical shape. The exponents describing the power-law temporal decay of number of particles and energy as well as dependence of velocity correlations on mass are determined using large-scale Monte Carlo simulations. It is shown that the exponents are universal only for the point-particle model. In the other two cases, the exponents are dependent on the initial number density and correlations vanish at high number densities. The fractal dimension for the second model is close to 1.49.

High-energy velocity tails in uniformly heated granular materials

We experimentally investigate the velocity distributions of quasi-two-dimensional granular materials uniformly heated by an electromagnetic vibrator, where the translational velocity and the rotation of a single particle are Gaussian and independent. We observe the non-Gaussian distributions of particle velocity, with the density-independent high-energy tails characterized by an exponent of β=1.50±0.03 for volume fractions of 0.111≤ϕ≤0.832, covering a wide range of structures and dynamics. Surprisingly, our results are not only in excellent agreement with the prediction of the kinetic theories of granular gas but also hold for an extremely high-volume fraction of ϕ=0.832 where the granular material forms a crystalline solid and the kinetic theory of granular gas fails fantastically. Our experiment suggests that the density-independent high-energy velocity tails of β=1.50 are a characteristic of uniformly heated granular matter.

Dimension dependence of clustering dynamics in models of ballistic aggregation and freely cooling granular gas

Via event-driven molecular dynamics simulations we study kinetics of clustering in assemblies of inelastic particles in various space dimensions. We consider two models, viz., the ballistic aggregation model (BAM) and the freely cooling granular gas model (GGM), for each of which we quantify the time dependence of kinetic energy and average mass of clusters (that form due to inelastic collisions). These quantities, for both the models, exhibit power-law behavior, at least in the long time limit. For the BAM, corresponding exponents exhibit strong dimension dependence and follow a hyperscaling relation. In addition, in the high packing fraction limit the behavior of these quantities become consistent with a scaling theory that predicts an inverse relation between energy and mass. On the other hand, in the case of the GGM we do not find any evidence for such a picture. In this case, even though the energy decay, irrespective of packing fraction, matches quantitatively with that for the high packing fraction picture of the BAM, it is inversely proportional to the growth of mass only in one dimension, and the growth appears to be rather insensitive to the choice of the dimension, unlike the BAM.

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.

Early-stage aggregation in three-dimensional charged granular gas

Neutral grains made of the same dielectric material can attain considerable charges due to collisions and generate long-range interactions. We perform molecular dynamic simulations in three dimensions for a dilute, freely cooling granular gas of viscoelastic particles that exchange charges during collisions. As compared to the case of clustering of viscoelastic particles solely due to dissipation, we find that the electrostatic interactions due to collisional charging alter the characteristic size, morphology, and growth rate of the clusters. The average cluster size grows with time as a power law, whose exponent is relatively larger in the charged gas than the neutral case. The growth of the average cluster size is found to be independent of the ratio of characteristic Coulomb to kinetic energy, or equivalently, of the typical Bjerrum length. However, this ratio alters the crossover time of the growth. Both simulations and mean-field calculations based on Smoluchowski's equation suggest that a suppression of particle diffusion due to the electrostatic interactions helps in the aggregation process.

Ballistic aggregation in systems of inelastic particles: Cluster growth, structure, and aging

We study far-from-equilibrium dynamics in models of freely cooling granular gas and ballistically aggregating compact clusters. For both the cases, from event-driven molecular dynamics simulations, we have presented detailed results on structure and dynamics in space dimensions d=1 and 2. Via appropriate analyses it has been confirmed that the ballistic aggregation mechanism applies in d=1 granular gases as well. Aging phenomena for this mechanism, in both the dimensions, have been studied via the two-time density autocorrelation function. This quantity is demonstrated to exhibit scaling property similar to that in the standard phase transition kinetics. The corresponding functional forms have been quantified and the outcomes have been discussed in connection with the structural properties. Our results on aging establish a more complete equivalence between the granular gas and the ballistic aggregation models in d=1.

Rattler-induced aging dynamics in jammed granular systems

Granular materials jam when developing a network of contact forces able to resist the applied stresses. Through numerical simulations of the dynamics of the jamming process, we show that the jamming transition does not occur when the kinetic energy vanishes. Rather, as the system jams, the kinetic energy becomes dominated by rattler particles, which scatter within their cages. The relaxation of the kinetic energy in the jammed configuration exhibits a double power-law decay, which we interpret in terms of the interplay between backbone and rattler particles.

Shock propagation in locally driven granular systems

We study shock propagation in a system of initially stationary hard spheres that is driven by a continuous injection of particles at the origin. The disturbance created by the injection of energy spreads radially outward through collisions between particles. Using scaling arguments, we determine the exponent characterizing the power-law growth of this disturbance in all dimensions. The scaling functions describing the various physical quantities are determined using large-scale event-driven simulations in two and three dimensions for both elastic and inelastic systems. The results are shown to describe well the data from two different experiments on granular systems that are similarly driven.

Velocity distribution of a driven inelastic one-component Maxwell gas

The nature of the velocity distribution of a driven granular gas, though well studied, is unknown as to whether it is universal or not, and, if universal, what it is. We determine the tails of the steady state velocity distribution of a driven inelastic Maxwell gas, which is a simple model of a granular gas where the rate of collision between particles is independent of the separation as well as the relative velocity. We show that the steady state velocity distribution is nonuniversal and depends strongly on the nature of driving. The asymptotic behavior of the velocity distribution is shown to be identical to that of a noninteracting model where the collisions between particles are ignored. For diffusive driving, where collisions with the wall are modeled by an additive noise, the tails of the velocity distribution is universal only if the noise distribution decays faster than exponential.

Driven inelastic Maxwell gases

We consider the inelastic Maxwell model, which consists of a collection of particles that are characterized by only their velocities and evolving through binary collisions and external driving. At any instant, a particle is equally likely to collide with any of the remaining particles. The system evolves in continuous time with mutual collisions and driving taken to be point processes with rates τ(c)(-1) and τ(w)(-1), respectively. The mutual collisions conserve momentum and are inelastic, with a coefficient of restitution r. The velocity change of a particle with velocity v, due to driving, is taken to be Δv=-(1+r(w))v+η, where r(w)∈[-1,1] and η is Gaussian white noise. For r(w)∈(0,1], this driving mechanism mimics the collision with a randomly moving wall, where r(w) is the coefficient of restitution. Another special limit of this driving is the so-called Ornstein-Uhlenbeck process given by dv/dt=-Γv+η. We show that while the equations for the n-particle velocity distribution functions (n=1,2,...) do not close, the joint evolution equations of the variance and the two-particle velocity correlation functions close. With the exact formula for the variance we find that, for r(w)≠-1, the system goes to a steady state. Also we obtain the exact tail of the velocity distribution in the steady state. On the other hand, for r(w)=-1, the system does not have a steady state. Similarly, the system goes to a steady state for the Ornstein-Uhlenbeck driving with Γ≠0, whereas for the purely diffusive driving (Γ=0), the system does not have a steady state.

Instabilities in granular binary mixtures at moderate densities

A linear stability analysis of the Navier-Stokes (NS) granular hydrodynamic equations is performed to determine the critical length scale for the onset of vortices and clusters instabilities in granular dense binary mixtures. In contrast to previous attempts, our results (which are based on the solution to the inelastic Enskog equation to NS order) are not restricted to nearly elastic systems since they take into account the complete nonlinear dependence of the NS transport coefficients on the coefficients of restitution α(ij). The theoretical predictions for the critical length scales are compared to molecular dynamics (MD) simulations in flows of strong dissipation (α(ij) ≥ 0.7) and moderate solid volume fractions (ϕ ≤ 0.2). We find excellent agreement between MD and kinetic theory for the onset of velocity vortices, indicating the applicability of NS hydrodynamics to polydisperse flows even for strong inelasticity, finite density, and particle dissimilarity.