Dynamics of freely cooling granular gases

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.

Granular cooling of ellipsoidal particles in microgravity

A three-dimensional granular gas of ellipsoids is established by exposing the system to the microgravity environment of the International Space Station. We use two methods to measure the dynamics of the constituent particles and report the long-time development of the granular temperature until no further particle movement is detectable. The resulting cooling behavior can be well described by Haff’s cooling law with time scale τ. Different analysis methods show evidence of particle clustering towards the end of the experiment. By using the kinetic theory for ellipsoids we compare the translational energy dissipation of individual collision events with the overall cooling time scale τ. The difference from this comparison indicates how energy is distributed in different degrees of freedom including both translation and rotation during the cooling.

Hydrodynamics of granular particles on a line

We investigate a lattice model representing a granular gas in a thin channel. We deduce the hydrodynamic description for the model from the microscopic dynamics in the large-system limit, including the lowest finite-size corrections. The main prediction from hydrodynamics, when finite-size corrections are neglected, is the existence of a steady "uniform longitudinal flow" (ULF), with the granular temperature and the velocity gradient both uniform and directly related. Extensive numerical simulations of the system show that such a state can be observed in the bulk of a finite-size system by attaching two thermostats with the same temperature at its boundaries. The relation between the ULF state and the shocks appearing in the late stage of a cooling gas of inelastic hard rods is discussed.

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.

Energy decay in a tapped granular column: Can a one-dimensional toy model provide insight into fully three-dimensional systems?

The decay of energy within particulate media subjected to an impulse is an issue of significant scientific interest, but also one with numerous important practical applications. In this paper, we study the dynamics of a granular system exposed to energetic impulses in the form of discrete taps from a solid surface. By considering a one-dimensional toy system, we develop a simple theory, which successfully describes the energy decay within the system following exposure to an impulse. We then extend this theory so as to make it applicable also to more realistic, three-dimensional granular systems, assessing the validity of the model through direct comparison with discrete particle method simulations. The theoretical form presented possesses several notable consequences; in particular, it is demonstrated that for suitably large systems, effects due to the bounding walls may be entirely neglected. We also establish the existence of a threshold system size above which a granular bed may be considered fully three dimensional.

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.

Clustering and velocity distributions in granular gases cooling by solid friction

We present large-scale molecular dynamics simulations to study the free evolution of granular gases. Initially, the density of particles is homogeneous and the velocity follows a Maxwell-Boltzmann (MB) distribution. The system cools down due to solid friction between the granular particles. The density remains homogeneous, and the velocity distribution remains MB at early times, while the kinetic energy of the system decays with time. However, fluctuations in the density and velocity fields grow, and the system evolves via formation of clusters in the density field and the local ordering of velocity field, consistent with the onset of plug flow. This is accompanied by a transition of the velocity distribution function from MB to non-MB behavior. We used equal-time correlation functions and structure factors of the density and velocity fields to study the morphology of clustering. From the correlation functions, we obtain the cluster size, L, as a function of time, t. We show that it exhibits power law growth with L(t)∼t^{1/3}.

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.

Blast Dynamics in a Dissipative Gas

The blast caused by an intense explosion has been extensively studied in conservative fluids, where the Taylor-von Neumann-Sedov hydrodynamic solution is a prototypical example of self-similarity driven by conservation laws. In dissipative media, however, energy conservation is violated, yet a distinctive self-similar solution appears. It hinges on the decoupling of random and coherent motion permitted by a broad class of dissipative mechanisms. This enforces a peculiar layered structure in the shock, for which we derive the full hydrodynamic solution, validated by a microscopic approach based on molecular dynamics simulations. We predict and evidence a succession of temporal regimes, as well as a long-time corrugation instability, also self-similar, which disrupts the blast boundary. These generic results may apply from astrophysical systems to granular gases, and invite further cross-fertilization between microscopic and hydrodynamic approaches of shock waves.

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.

Higher-order effects on orientational correlation and relaxation dynamics in homogeneous cooling of a rough granular gas

The orientational or angular correlation between the directions of the translational and rotational motions is analyzed theoretically for the homogeneous cooling state of a rough granular gas. The dynamical equations are derived using an approximate form of the single-particle distribution function that incorporates angular correlations. The goal is to assess the effects of higher-order angular corrections for which both quadratic- and quartic-order terms (in translational and rotational velocities of particles) are retained in the perturbation expansion of the distribution function. We show that higher-order corrections can markedly affect the steady-state orientational correlation when the normal restitution coefficient is moderate or small, and this effect is more prominent for nearly smooth particles. The transient evolution of orientational correlation is found to be significantly affected by higher-order terms. In particular the higher-order orientational correlations can dominate over the leading-order contribution during short times even in the quasi-elastic limit, although the steady correlation remains unaffected by such corrections in the same limit. The buildup of correlations during the transient stage seems to be closely tied to the evolution of the ratio between the rotational and translational temperatures. It is demonstrated that the transient dynamics of the temperature ratio and its steady state remain insensitive to higher-order angular correlation.

Energy decay in three-dimensional freely cooling granular gas

The kinetic energy of a freely cooling granular gas decreases as a power law t(-θ) at large times t. Two theoretical conjectures exist for the exponent θ. One based on ballistic aggregation of compact spherical aggregates predicts θ=2d/(d+2) in d dimensions. The other based on Burgers equation describing anisotropic, extended clusters predicts θ=d/2 when 2≤d≤4. We do extensive simulations in three dimensions to find that while θ is as predicted by ballistic aggregation, the cluster statistics and velocity distribution differ from it. Thus, the freely cooling granular gas fits to neither the ballistic aggregation or a Burgers equation description.

Shock propagation in granular flow subjected to an external impact

We analyze a recent experiment [Boudet, Cassagne, and Kellay, Phys. Rev. Lett. 103, 224501 (2009)] in which the shock created by the impact of a steel ball on a flowing monolayer of glass beads is studied quantitatively. We argue that radial momentum is conserved in the process and hence show that in two dimensions the shock radius increases in time t as a power law t{1/3}. This is confirmed in event driven simulations of an inelastic hard sphere system. The experimental data are compared with the theoretical prediction and are shown to compare well at intermediate times. At long times the experimental data exhibit a crossover to a different scaling behavior. We attribute this to the problem becoming effectively three dimensional due to the accumulation of particles at the shock front and propose a simple hard sphere model that incorporates this effect. Simulations of this model capture the crossover seen in the experimental data.

Coarse-grained dynamics of the freely cooling granular gas in one dimension

We study the dynamics and structure of clusters in the inhomogeneous clustered regime of a freely cooling granular gas of point particles in one dimension. The coefficient of restitution is modeled as r(0)<1 or 1, depending on whether the relative speed is greater or smaller than a velocity scale δ. The effective fragmentation rate of a cluster is shown to rise sharply beyond a δ-dependent time scale. This crossover is coincident with the velocity fluctuations within a cluster becoming order δ. Beyond this crossover time, the cluster-size distribution develops a nontrivial power-law distribution, whose scaling properties are related to those of the velocity fluctuations. We argue that these underlying features are responsible for the recently observed nontrivial coarsening behavior in the one-dimensional freely cooling granular gas.

Navier-Stokes hydrodynamics of thermal collapse in a freely cooling granular gas

We show that, in dimension higher than one, heat diffusion and viscosity cannot arrest thermal collapse in a freely evolving dilute granular gas, even in the absence of gravity. Thermal collapse involves a finite-time blowup of the gas density. It was predicted earlier in ideal, Euler hydrodynamics of dilute granular gases in the absence of gravity, and in nonideal, Navier-Stokes granular hydrodynamics in the presence of gravity. We determine, analytically and numerically, the dynamic scaling laws that characterize the gas flow close to collapse. We also investigate bifurcations of a freely evolving dilute granular gas in circular and wedge-shaped containers. Our results imply that, in general, thermal collapse can only be arrested when the gas density becomes comparable with the close-packing density of grains. This provides a natural explanation to the formation of densely packed clusters of particles in a variety of initially dilute granular flows.

Dilute wet granular particles: nonequilibrium dynamics and structure formation

We investigate a gas of wet granular particles covered by a thin liquid film. The dynamic evolution is governed by two-particle interactions, which are mainly due to interfacial forces in contrast to dry granular gases. When two wet grains collide, a capillary bridge is formed and stays intact up to a certain distance of withdrawal when the bridge ruptures, dissipating a fixed amount of energy. A freely cooling system is shown to undergo a nonequilibrium dynamic phase transition from a state with mainly single particles and fast cooling to a state with growing aggregates such that bridge rupture becomes a rare event and cooling is slow. In the early stage of cluster growth, aggregation is a self-similar process with a fractal dimension of the aggregates approximately equal to Df approximately 2 . At later times, a percolating cluster is observed which ultimately absorbs all the particles. The final cluster is compact on large length scales, but fractal with Df approximately 2 on small length scales.

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.

Experimental investigation of the freely cooling granular gas

We investigate the dynamics of the freely cooling granular gas. For this purpose we diamagnetically levitate the grains providing a terrestrial milligravity environment. At early times we find good agreement with Haff's law, where the time scale for particle collisions can be determined from independent measurements. At late times, clustering of particles occurs. This can be included in a Haff-like description taking into account the decreasing number of free particles. At very late times, only a single particle determines the dynamics, which is again described by a version of Haff's law. With this a good description of the data is possible over the whole time range.

Attempted density blowup in a freely cooling dilute granular gas: hydrodynamics versus molecular dynamics

It has been recently shown [I. Fouxon, Phys. Rev. E 75, 050301(R) (2007); I. Fouxon, Phys. Fluids 19, 093303 (2007)] that, in the framework of ideal granular hydrodynamics (IGHD), an initially smooth hydrodynamic flow of a granular gas can produce an infinite gas density in a finite time. Exact solutions that exhibit this property have been derived. Close to the singularity, the granular gas pressure is finite and almost constant. We report molecular dynamics (MD) simulations of a freely cooling gas of nearly elastically colliding hard disks, aimed at identifying the "attempted" density blowup regime. The initial conditions of the simulated flow mimic those of one particular solution of the IGHD equations that exhibits the density blowup. We measure the hydrodynamic fields in the MD simulations and compare them with predictions from the ideal theory. We find a remarkable quantitative agreement between the two over an extended time interval, proving the existence of the attempted blowup regime. As the attempted singularity is approached, the hydrodynamic fields, as observed in the MD simulations, deviate from the predictions of the ideal solution. To investigate the mechanism of breakdown of the ideal theory near the singularity, we extend the hydrodynamic theory by accounting separately for the gradient-dependent transport and for finite density corrections.

Nonlinear theory of nonstationary low Mach number channel flows of freely cooling nearly elastic granular gases

We employ hydrodynamic equations to investigate nonstationary channel flows of freely cooling dilute gases of hard and smooth spheres with nearly elastic particle collisions. This work focuses on the regime where the sound travel time through the channel is much shorter than the characteristic cooling time of the gas. As a result, the gas pressure rapidly becomes almost homogeneous, while the typical Mach number of the flow drops well below unity. Eliminating the acoustic modes and employing Lagrangian coordinates, we reduce the hydrodynamic equations to a single nonlinear and nonlocal equation of a reaction-diffusion type. This equation describes a broad class of channel flows and, in particular, can follow the development of the clustering instability from a weakly perturbed homogeneous cooling state to strongly nonlinear states. If the heat diffusion is neglected, the reduced equation becomes exactly soluble, and the solution develops a finite-time density blowup. The blowup has the same local features at singularity as those exhibited by the recently found family of exact solutions of the full set of ideal hydrodynamic equations [I. Fouxon, Phys. Rev. E 75, 050301(R) (2007); I. Fouxon,Phys. Fluids 19, 093303 (2007)]. The heat diffusion, however, always becomes important near the attempted singularity. It arrests the density blowup and brings about previously unknown inhomogeneous cooling states (ICSs) of the gas, where the pressure continues to decay with time, while the density profile becomes time-independent. The ICSs represent exact solutions of the full set of granular hydrodynamic equations. Both the density profile of an ICS and the characteristic relaxation time toward it are determined by a single dimensionless parameter L that describes the relative role of the inelastic energy loss and heat diffusion. At L>>1 the intermediate cooling dynamics proceeds as a competition between "holes": low-density regions of the gas. This competition resembles Ostwald ripening (only one hole survives at the end), and we report a particular regime where the "hole ripening" statistics exhibits a simple dynamic scaling behavior.

Violation of the Porod law in a freely cooling granular gas in one dimension

We study a model of freely cooling inelastic granular gas in one dimension, with a restitution coefficient which approaches the elastic limit below a relative velocity scale delta. While at early times (t<<delta;{-1}) the gas behaves as a completely inelastic sticky gas conforming to predictions of earlier studies, at late times (t>>delta;{-1}) it exhibits a new fluctuation-dominated phase ordering state. We find distinct scaling behavior for the (i) density distribution function, (ii) occupied and empty gap distribution functions, (iii) the density structure function, and (iv) the velocity structure function, as compared to the completely inelastic sticky gas. The spatial structure functions (iii) and (iv) violate the Porod law. Within a mean-field approximation, the exponents describing the structure functions are related to those describing the spatial gap distribution functions.

Dynamical collision network in granular gases

We address the problem of recollisions in cooling granular gases. To this aim, we dynamically construct the interaction network in a granular gas, using the sequence of collisions collected in an event driven simulation of inelastic hard disks from time 0 until time t . The network is decomposed into its k -core structure: particles in a core of index k have collided at least k times with other particles in the same core. The difference between cores k+1 and k is the so-called k -shell, and the set of all shells is a complete and nonoverlapping decomposition of the system. Because of energy dissipation, the gas cools down: its initial spatially homogeneous dynamics, characterized by the Haff law, i.e., a t{-2} energy decay, is unstable toward a strongly inhomogeneous phase with clusters and vortices, where energy decays as t{-1} . The clear transition between those two phases appears in the evolution of the k -shells structure in the collision network. In the homogeneous state the k -shell structure evolves as in a growing network with a fixed number of vertices and randomly added links: the shell distribution is strongly peaked around the most populated shell, which has an index k{max} approximately 0.9(d) with (d) the average number of collisions experienced by a particle. During the final nonhomogeneous state a growing fraction of collisions is concentrated in small, almost closed, communities of particles: k{max} is no more linear in (d) and the distribution of shells becomes extremely large developing a power-law tail approximately k{-3} for high shell indexes. We conclude proposing a simple algorithm to build a correlated random network that reproduces, with few essential ingredients, the whole observed phenomenology, including the t{-1} energy decay. It consists of two kinds of collisions (links): single random collisions with any other particle and long chains of recollisions with only previously encountered particles. The algorithm disregards the exact spatial arrangement of clusters, suggesting that the observed stringlike structures are not essential to determine the statistics of recollisions and the energy decay.

Formation and evolution of density singularities in hydrodynamics of inelastic gases

We use hydrodynamics to investigate clustering in a gas of inelastically colliding spheres. The hydrodynamic equations exhibit a finite-time density blowup, where the gas pressure remains finite. The density blowup signals the formation of close-packed clusters. The blowup dynamics is universal and describable by exact analytic solutions continuable beyond the blowup time. These solutions show that dilute hydrodynamic equations yield a powerful effective description of a granular gas flow with close-packed clusters, described as finite-mass pointlike singularities of the density. This description is similar in spirit to the description of shocks in ordinary ideal gas dynamics.

Dry and wet granular shock waves

The formation of a shock wave in one-dimensional granular gases is considered, for both the dry and the wet cases, and the results are compared with the analytical shock wave solution in a sticky gas. Numerical simulations show that the behavior of the shock wave in both cases tends asymptotically to the sticky limit. In the inelastic gas (dry case) there is a very close correspondence to the sticky gas, with one big cluster growing in the center of the shock wave, and a step-like stationary velocity profile. In the wet case, the shock wave has a nonzero width which is marked by two symmetric heavy clusters performing breathing oscillations with slowly increasing amplitude. All three models have the same asymptotic energy dissipation law, which is important in the context of the free cooling scenario. For the early stage of the shock formation and asymptotic oscillations we provide analytical results as well.

Velocity distributions and aging in a cooling granular gas

We use large-scale molecular dynamics simulations to study freely evolving granular gases with dimensionality d=2,3 . The system dissipates kinetic energy (or cools) due to inelastic collisions between granular particles. The density and velocity fields are approximately homogeneous at early times, and the system is said to be in a homogeneous cooling state (HCS). However, fluctuations in the density and velocity fields grow, and the system evolves into an inhomogeneous cooling state (ICS). We study the nature of velocity distributions in both the HCS and ICS. We also investigate the aging property of the velocity autocorrelation function.

Free cooling of the one-dimensional wet granular gas

The free cooling behavior of a wet granular gas is studied in one dimension. We employ a particularly simple model system in which the interaction of wet grains is characterized by a fixed energy loss assigned to each collision. Macroscopic laws of energy dissipation and cluster formation are studied on the basis of numerical simulations and mean-field analytical calculations. We find a number of remarkable scaling properties which may shed light on earlier unexplained results for related systems.

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.

Power-law velocity distributions in granular gases

The kinetic theory of granular gases is studied for spatially homogeneous systems. At large velocities, the equation governing the velocity distribution becomes linear, and it admits stationary solutions with a power-law tail, f (v) approximately v(-sigma) . This behavior holds in arbitrary dimension for arbitrary collision rates including both hard spheres and Maxwell molecules. Numerical simulations show that driven steady states with the same power-law tail can be realized by injecting energy into the system at very high energies. In one dimension, we also obtain self-similar time-dependent solutions where the velocities collapse to zero. At small velocities there is a steady state and a power-law tail but at large velocities, the behavior is time dependent with a stretched exponential decay.

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.

Hydrodynamic singularities and clustering in a freely cooling inelastic gas

We employ hydrodynamic equations to follow the clustering instability of a freely cooling dilute gas of inelastically colliding spheres into a well-developed nonlinear regime. We simplify the problem by dealing with a one-dimensional coarse-grained flow. We observe that at a late stage of the instability the shear stress becomes negligibly small, and the gas flows solely by inertia. As a result the flow formally develops a finite-time singularity, as the velocity gradient and the gas density diverge at some location. We argue that flow by inertia represents a generic intermediate asymptotic of unstable free cooling of dilute inelastic gases.

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.

Percolation and Burgers' dynamics in a model of capillary formation

Capillary networks are essential in vertebrates to supply tissues with nutrients. Experiments of in vitro capillary formation show that cells randomly spread on a gel matrix autonomously organize to form vascular networks. Cells form disconnected networks at low densities and connected ones above a critical density. Above the critical density the network is characterized by a typical mesh size approximately 200 microm, which is approximately constant on a wide range of density values. In this paper we present a full characterization of a recently proposed model which reproduces the main features of the biological system, focusing on its dynamical properties, on the fractal properties of patterns, and on the percolative phase transition. We discuss the relevance of the model in relation with some experiments in living beings and proposed diagnostic methods based on the measurement of the fractal dimension of vascular networks.

Role of friction in pattern formation in oscillated granular layers

Particles in granular flows are often modeled as frictionless (smooth) inelastic spheres; however, there exist no frictionless grains, just as there are no elastic grains. Our molecular dynamics simulations reveal that friction is essential for realistic modeling of vertically oscillated granular layers: simulations of frictionless particles yield patterns with an onset at a container acceleration about 30% smaller than that observed in experiments and simulations with friction. More importantly, even though square and hexagonal patterns form for a wide range of the oscillation parameters in experiments and in our simulations of frictional inelastic particles, only stripe patterns form in the simulations without friction, even if the inelasticity is increased to obtain as much dissipation as in frictional particles. We also consider the effect of particle friction on the shock wave that forms each time the granular layer strikes the container. While a shock wave still forms for frictionless particles, the spatial and temporal dependence of the hydrodynamic fields differ for the cases with and without friction.

Cluster growth in two- and three-dimensional granular gases

Dissipation in granular media leads to interesting phenomena such as cluster formation and crystallization in nonequilibrium dynamical states. The freely cooling system is examined concerning the energy decay and the cluster evolution in time, both in two and three dimensions. We also suggest an interpretation of the three-dimensional cluster growth in terms of percolation theory, but this point deserves further study.

Steady-state velocity distributions of an oscillated granular gas

We use a three-dimensional molecular dynamics simulation to study the single particle distribution function of a dilute granular gas driven by a vertically oscillating plate at high accelerations (15g-90g). We find that the density and the temperature fields are essentially time-invariant above a height of about 40 particle diameters, where typically 20% of the grains are contained. These grains form the nonequilibrium steady-state granular gas with a Knudsen number unity or greater. In the steady-state region, the probability distribution function of the horizontal velocity c(x) (scaled by the local horizontal temperature) is found to be nearly independent of height, even though the hydrodynamic fields vary with height. We find that the high energy tails of the distribution function are described by a stretched exponential approximately exp(-Bcalphax), where alpha depends on the restitution coefficient e and falls in the range 1.2<alpha<1.6. However, alpha does not vary significantly for a wide range of friction coefficient values. We find that the distribution function of a frictionless inelastic hard sphere model can be made similar to that of a frictional model by adjusting e. However, there is no single value of e that mimics the frictional model over a range of heights.

Granular clustering in a hydrodynamic simulation

We examine the hydrodynamics of a granular gas using numerical simulation. We demonstrate the appearance of shearing and clustering instabilities predicted by linear stability analysis, and show that their appearance is directly related to the inelasticity of collisions in the material. We discuss the rate at which these instabilities arise and the manner in which clusters grow and merge.

Coarsening of granular clusters: Two types of scaling behaviors

We report on an experimental study of small cluster dynamics during the coarsening process in driven granular submonolayers of 120-microm bronze particles. The techniques of electrostatic and vertical mechanical vibration were employed to excite the granular gas. We measure the scaling exponent for the evaporation of small clusters during coarsening. It was found that the surface area of small clusters S vs time t behaves as S to (t(0)-t)(2/3) for lower frequencies and S to (t(0)-t) for higher frequencies. We argue that the change in the scaling exponent is related to the transition from three-dimensional (3D) to 2D character of motion in the granular gas.