Cluster-Growth in Freely Cooling Granular Media

When dissipative particles are left alone, their fluctuation energy decays due to collisional interactions, clusters build up and grow with time until the system size is reached. When the effective dissipation is strong enough, this may lead to the `inelastic collapse', i.e. the divergence of the collision frequency of some particles. The cluster growth is an interesting physical phenomenon, whereas the inelastic collapse is an intrinsic effect of the inelastic hard sphere (IHS) model used to study the cluster growth - involving only a negligible number of particles in the system. Here, we extend the IHS model by introducing an elastic contact energy and the related contact duration t_c. This avoids the inelastic collapse and allows to examine the long-time behavior of the system. For a quantitative description of the cluster growth, we propose a burning - like algorithm in continuous space, that readily identifies all particles that belong to the same cluster. The criterion for this is here chosen to be only the particle distance. With this method we identify three regimes of behavior. First, for short times a homogeneous cooling state (HCS) exists, where a mean-field theory works nicely, and the clusters are tiny and grow very slowly. Second, at a certain time which depends on the system's properties, cluster growth starts and the clusters increase in size and mass until, in the third regime, the system size is reached and most of the particles are collected in one huge cluster.Comment: 16 pages, 21 figures. Chaos 9(3) (in press, 1999

The Simplest Piston Problem II: Inelastic Collisions

We study the dynamics of three particles in a finite interval, in which two light particles are separated by a heavy ``piston'', with elastic collisions between particles but inelastic collisions between the light particles and the interval ends. A symmetry breaking occurs in which the piston migrates near one end of the interval and performs small-amplitude periodic oscillations on a logarithmic time scale. The properties of this dissipative limit cycle can be understood simply in terms of an effective restitution coefficient picture. Many dynamical features of the three-particle system closely resemble those of the many-body inelastic piston problem.Comment: 8 pages, 7 figures, 2-column revtex4 forma

The energy flux into a fluidized granular medium at a vibrating wall

We study the power input of a vibrating wall into a fluidized granular medium, using event driven simulations of a model granular system. The system consists of inelastic hard disks contained between a stationary and a vibrating elastic wall, in the absence of gravity. Two scaling relations for the power input are found, both involving the pressure. The transition between the two occurs when waves generated at the moving wall can propagate across the system. Choosing an appropriate waveform for the vibrating wall removes one of these scalings and renders the second very simple.Comment: 5 pages, revtex, 7 postscript figure

Phase transition in inelastic disks

This letter investigates the molecular dynamics of inelastic disks without external forcing. By introducing a new observation frame with a rescaled time, we observe the virtual steady states converted from asymptotic energy dissipation processes. System behavior in the thermodynamic limit is carefully investigated. It is found that a phase transition with symmetry breaking occurs when the magnitude of dissipation is greater than a critical value.Comment: 9 pages, 6 figure

Towards a continuum theory of clustering in a freely cooling inelastic gas

We performed molecular dynamics simulations to investigate the clustering instability of a freely cooling dilute gas of inelastically colliding disks in a quasi-one-dimensional setting. We observe that, as the gas cools, the shear stress becomes negligibly small, and the gas flows by inertia only. Finite-time singularities, intrinsic in such a flow, are arrested only when close-packed clusters are formed. We observe that the late-time dynamics of this system are describable by the Burgers equation with vanishing viscosity, and predict the long-time coarsening behavior.Comment: 7 pages, 5 eps figures, to appear in Europhys. Let

Energy flows in vibrated granular media

We study vibrated granular media, investigating each of the three components of the energy flow: particle-particle dissipation, energy input at the vibrating wall, and particle-wall dissipation. Energy dissipated by interparticle collisions is well estimated by existing theories when the granular material is dilute, and these theories are extended to include rotational kinetic energy. When the granular material is dense, the observed particle-particle dissipation rate decreases to as little as 2/5 of the theoretical prediction. We observe that the rate of energy input is the weight of the granular material times an average vibration velocity times a function of the ratio of particle to vibration velocity. `Particle-wall' dissipation has been neglected in all theories up to now, but can play an important role when the granular material is dilute. The ratio between gravitational potential energy and kinetic energy can vary by as much as a factor of 3. Previous simulations and experiments have shown that E ~ V^delta, with delta=2 for dilute granular material, and delta ~ 1.5 for dense granular material. We relate this change in exponent to the departure of particle-particle dissipation from its theoretical value.Comment: 19 pages revtex, 10 embedded eps figures, accepted by PR

Influence of Canal Geometry and Dynamics on Controllability

This paper presents the results of the Task Committee on Canal Automation Algorithms with regard to the influence of canal properties on the controllability of irrigation canals. While the control provided by individual algorithms was not evaluated, studies were performed to illustrate inherent hydraulic limitations—the inability of canal pools to recover rapidly from disturbances or flow perturbations. Studies were performed in nondimensional form to develop a better understanding of how pool properties influence pool response. Three such studies were performed. First, nondimensional backwater curves were developed for a range of canal conditions. The second study involved the propagation of waves initiated at the upstream end of a canal pool, as this is influenced by downstream boundary conditions. Finally, the response of pools to downstream withdrawals was examined in terms of their sluggish recovery even when the correct flow change is applied upstream. These results will help in understanding how canal properties influence the ability of operators to effectively control a canal either manually or automatically, and should influence future design practices

Dynamics of Freely Cooling Granular Gases

We study dynamics of freely cooling granular gases in two-dimensions using large-scale molecular dynamics simulations. We find that for dilute systems the typical kinetic energy decays algebraically with time, E(t) ~ t^{-1}, in the long time limit. Asymptotically, velocity statistics are characterized by a universal Gaussian distribution, in contrast with the exponential high-energy tails characterizing the early homogeneous regime. We show that in the late clustering regime particles move coherently as typical local velocity fluctuations, Delta v, are small compared with the typical velocity, Delta v/v ~ t^{-1/4}. Furthermore, locally averaged shear modes dominate over acoustic modes. The small thermal velocity fluctuations suggest that the system can be heuristically described by Burgers-like equations.Comment: 4 pages, 5 figure

Spatial Correlations in Compressible Granular Flows

For a freely evolving granular fluid, the buildup of spatial correlations in density and flow field is described using fluctuating hydrodynamics. The theory for incompressible flows is extended to the general, compressible case, including longitudinal velocity and density fluctuations, and yields qualitatively different results for long range correlations. The structure factor of density fluctuations shows a maximum at finite wavenumber, shifting in time to smaller wavenumbers and corresponding to a growing correlation length. It agrees well with two-dimensional molecular dynamics simulations.Comment: 12 pages, Latex, 3 figure