Symmetry breaking and coarsening of clusters in a prototypical driven granular gas

Noise induces rare events in granular media

The granular Leidenfrost effect [B. Meerson, et al., Phys. Rev. Lett. 91, 024301 (2003)PRLTAO0031-900710.1103/PhysRevLett.91.024301; P. Eshuis et al., Phys. Rev. Lett. 95, 258001 (2005)PRLTAO0031-900710.1103/PhysRevLett.95.258001] is the levitation of a mass of granular matter when a wall below the grains is vibrated, giving rise to a hot granular gas below the cluster. We find by simulation that for a range of parameters the system is bistable: the levitated cluster can occasionally break and give rise to two clusters and a hot granular gas above and below. We use techniques from the theory of rare events to compute the mean transition time for breaking to occur. This requires the introduction of a two-component reaction coordinate.

Liquid-gas phase separation in confined vibrated dry granular matter

A new phase transition is observed experimentally in a dry granular gas subject to vertical vibration between two horizontal plates. Molecular dynamics simulations of this system allow us to investigate the observed phase separation in detail. We find a high-density, low temperature liquid, coexisting with a low-density, high temperature gas moving coherently. The importance of the coherent motion for phase separation is investigated using frequency modulation.

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.

Volume fluctuations and compressibility of a vibrated granular gas

The volume fluctuations in the steady state reached by a vibrated granular gas of hard particles confined by a movable piston on the top are investigated by means of event-driven simulations. Also, a compressibility factor, measuring the response in volume of the system to a change in the mass of the piston, is introduced and measured. From the second moment of the volume fluctuations and the compressibility factor, an effective temperature is defined by using the same relation as obeyed by equilibrium molecular systems. The interpretation of this effective temperature and its relationship with the granular temperature of the gas, and also with the velocity fluctuations of the movable piston, is discussed. It is found that the ratio of the temperature based on the volume fluctuations to the temperature based on the piston kinetic energy obeys simple dependencies on the inelasticity and on the piston-particle mass ratio.

Exploring the limits of granular hydrodynamics: a horizontal array of inelastic particles

The limits of granular hydrodynamics are explored in the context of the one-dimensional granular system introduced by Du, Li, and Kadanoff [Phys. Rev. Lett. 74, 1268 (1995)]. The density profile of the characteristic steady state, in which a single particle commutes between the driving wall and a dense cluster, is well captured by a hydrodynamic description provided that the finite size of the particles is incorporated. The temperature, however, is not well described: since all energy exchange is located at the border of the cluster, it is precisely for this quantity that the continuum approach breaks down.

Vibrated granular gas confined by a piston

The steady state of a vibrated granular gas confined by a movable piston on top is discussed. Particular attention is given to the hydrodynamic boundary conditions to be used when solving the inelastic Navier-Stokes equations. The relevance of an exact general condition relating the grain fluxes approaching and moving away from each of the walls is emphasized. It is shown how it can be used to get a consistent hydrodynamic description of the boundaries. The obtained expressions for the fields do not contain any undetermined parameter. Comparison of the theoretical predictions with molecular-dynamics simulation results is carried out, and a good agreement is observed for low density and not too large inelasticity. A practical way of introducing small finite density corrections to the dilute limit theory is proposed to improve the accuracy of the theory.

van der Waals-like phase-separation instability of a driven granular gas in three dimensions

We show that the van der Waals-like phase-separation instability of a driven granular gas at zero gravity, previously investigated in two-dimensional settings, persists in three dimensions. We consider a monodisperse granular gas driven by a thermal wall of a three-dimensional rectangular container at zero gravity. The basic steady state of this system, as described by granular hydrodynamic equations, involves a denser and colder layer of granulate located at the wall opposite to the driving wall. When the inelastic energy loss is sufficiently high, the driven granular gas exhibits, in some range of average densities, negative compressibility in the directions parallel to the driving wall. When the lateral dimensions of the container are sufficiently large, the negative compressibility causes spontaneous symmetry breaking of the basic steady state and a phase separation instability. Event-driven molecular dynamics simulations confirm and complement our theoretical predictions.

Phase separation of a driven granular gas in annular geometry

This work investigates phase separation of a monodisperse gas of inelastically colliding hard disks confined in a two-dimensional annulus, the inner circle of which represents a "thermal wall." When described by granular hydrodynamic equations, the basic steady state of this system is an azimuthally symmetric state of increased particle density at the exterior circle of the annulus. When the inelastic energy loss is sufficiently large, hydrodynamics predicts spontaneous symmetry breaking of the annular state, analogous to the van der Waals-like phase separation phenomenon previously found in a driven granular gas in rectangular geometry. At a fixed aspect ratio of the annulus, the phase separation involves a "spinodal interval" of particle area fractions, where the gas has negative compressibility in the azimuthal direction. The heat conduction in the azimuthal direction tends to suppress the instability, as corroborated by a marginal stability analysis of the basic steady state with respect to small perturbations. To test and complement our theoretical predictions we performed event-driven molecular dynamics (MD) simulations of this system. We clearly identify the transition to phase separated states in the MD simulations, despite large fluctuations present, by measuring the probability distribution of the amplitude of the fundamental Fourier mode of the azimuthal spectrum of the particle density. We find that the instability region, predicted from hydrodynamics, is always located within the phase separation region observed in the MD simulations. This implies the presence of a binodal (coexistence) region, where the annular state is metastable. The phase separation persists when the driving and elastic walls are interchanged, and also when the elastic wall is replaced by weakly inelastic one.

Convection in vibrated annular granular beds

The response to vibration of a granular bed, consisting of a standard cylindrical geometry but with the addition of a dissipative cylindrical inner wall, has been investigated both experimentally (using positron emission particle tracking) and numerically (using hard sphere molecular dynamics simulation). The packing fraction profiles and granular temperature distributions (in both vertical and horizontal directions) were determined as a function of height and distance from the axis. The two sets of results were in reasonable agreement. The molecular dynamics simulations were used to explore the behavior of the granular bed in the inner wall-outer wall coefficient of restitution phase space. It was observed that one could control the direction of the toroidal convection rolls by manipulating the relative dissipation at the inner and outer walls via the coefficients of restitution, and with several layers of grains it was seen that double convection rolls could also be formed, a result that was subsequently confirmed experimentally.

Local equation of state and velocity distributions of a driven granular gas

We present event-driven simulations of a granular gas of inelastic hard disks with incomplete normal restitution in two dimensions between vibrating walls (without gravity). We measure hydrodynamic quantities such as the stress tensor, density and temperature profiles, as well as velocity distributions. Relating the local pressure to the local temperature and local density, we construct a local constitutive equation. For strong inelasticities the local constitutive relation depends on global system parameters, like the volume fraction and the aspect ratio. For moderate inelasticities the constitutive relation is approximately independent of the system parameters and can hence be regarded as a local equation of state, even though the system is highly inhomogeneous with heterogeneous temperature and density profiles arising as a consequence of energy injection. With respect to local velocity distributions we find that they do not scale with the square root of the local granular temperature. Moreover the high-velocity tails are different for the distribution of the x and the y components of the velocity, and even depend on the position in the sample, the global volume fraction, and the coefficient of restitution.

Phase diagram of van der Waals-like phase separation in a driven granular gas

Equations of granular hydrostatics are used to compute the phase diagram of the recently discovered van der Waals-like phase separation in a driven granular gas. The model two-dimensional system consists of smooth hard disks in a rectangular box, colliding inelastically with each other and driven by a "thermal" wall at zero gravity. The spinodal line and the critical point of the phase separation are determined. Close to the critical point, the spinodal and binodal (coexistence) lines are determined analytically. Effects of the finite size of the confining box in the direction parallel to the thermal wall are investigated. These include suppression of the phase separation by heat conduction in the lateral direction and a change from supercritical to subcritical bifurcation.

Giant fluctuations at a granular phase separation threshold

We investigate a phase separation instability that occurs in a system of nearly elastically colliding hard spheres driven by a thermal wall. If the aspect ratio of the confining box exceeds a threshold value, granular hydrostatics predict phase separation: the formation of a high-density region coexisting with a low-density region along the wall that is opposite to the thermal wall. Event-driven molecular dynamics simulations confirm this prediction. The theoretical bifurcation curve agrees with the simulations quantitatively well below and well above the threshold. However, in a wide region of aspect ratios around the threshold, the system is dominated by fluctuations, and the hydrostatic theory breaks down. Two possible scenarios of the origin of the giant fluctuations are discussed.

Onset of thermal convection in a horizontal layer of granular gas

The Navier-Stokes granular hydrodynamics is employed for determining the threshold of thermal convection in an infinite horizontal layer of granular gas. The dependence of the convection threshold, in terms of the inelasticity of particle collisions, on the Froude and Knudsen numbers is found. A simple necessary condition for convection is formulated in terms of the Schwarzschild's criterion, well known in thermal convection of (compressible) classical fluids. The morphology of convection cells at the onset is determined. At large Froude numbers, the Froude number drops out of the problem. As the Froude number goes to zero, the convection instability turns into a recently discovered phase-separation instability.