Symmetry-breaking instability in a prototypical driven granular gas

Reduced Thermodynamic Description of Phase Separation in a Quasi-One-Dimensional Granular Gas

We describe simulations of a quasi-one-dimensional, vibrated granular gas which exhibits an apparent phase separation into a liquidlike phase and a gaslike phase. In thermal equilibrium, such a phase separation in one dimension is prohibited by entropic considerations. We propose that the granular gas minimizes a function of the conserved mechanical variables alone: the particle number and volume. Simulations in small cells can be used to extract the equation of state and predict the coexisting pressure and densities, as confirmation of the minimization principle. Fluctuations in the system manifest themselves as persistent density waves but they do not destroy the phase-separated state.

Phase separation in driven granular gases: exploring the elusive character of nonequilibrium steady states

The emergence of patterns and phase separation in many-body systems far from thermal equilibrium is discussed using the example of driven granular gases. It is shown that phase separation follows a similar mechanism as in the systems of active Brownian particles. Depending on the quantities chosen for observation, it may or may not be easy to find functionals analogous to the free energy in equilibrium statistical physics. We argue that although such functionals can always be derived from the dynamics, it is of only limited value for predicting relevant aspects of the nonequilibrium steady state of the system. Consequently, although there is indeed a 'principle' governing the selection of collective nonequilibrium steady states (and the corresponding large deviation functional can be identified), it is not generally useful for predicting the behaviour of the system.

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.

The minimization of mechanical work in vibrated granular matter

Experiments and computer simulations are carried out to investigate phase separation in a granular gas under vibration. The densities of the dilute and the dense phase are found to follow a lever rule and obey an equation of state. Here we show that the Maxwell equal-areas construction predicts the coexisting pressure and binodal densities remarkably well, even though the system is far from thermal equilibrium. This construction can be linked to the minimization of mechanical work associated with density fluctuations without invoking any concept related to equilibrium-like free energies.

Inhomogeneous quasistationary state of dense fluids of inelastic hard spheres

We study closed dense collections of freely cooling hard spheres that collide inelastically with constant coefficient of normal restitution. We find inhomogeneous states (ISs) where the density profile is spatially nonuniform but constant in time. The states are exact solutions of nonlinear partial differential equations that describe the coupled distributions of density and temperature valid when inelastic losses of energy per collision are small. The derivation is performed without modeling the equations' coefficients that are unknown in the dense limit (such as the equation of state) using only their scaling form specific for hard spheres. Thus the IS is the exact state of this dense many-body system. It captures a fundamental property of inelastic collections of particles: the possibility of preserving nonuniform temperature via the interplay of inelastic cooling and heat conduction that generalizes previous results. We perform numerical simulations to demonstrate that arbitrary initial state evolves to the IS in the limit of long times where the container has the geometry of the channel. The evolution is like a gas-liquid transition. The liquid condenses in a vanishing part of the total volume but takes most of the mass of the system. However, the gaseous phase, which mass grows only logarithmically with the system size, is relevant because its fast particles carry most of the energy of the system. Remarkably, the system self-organizes to dissipate no energy: The inelastic decay of energy is a power law [1+t/t(c)](-2), where t(c) diverges in the thermodynamic limit. This is reinforced by observing that for supercritical systems the IS coincide in most of the space with the steady states of granular systems heated at one of the walls. We discuss the relation of our results to the recently proposed finite-time singularity in other container's geometries.

Hydrodynamics of a vibrated granular monolayer

We investigate the long-standing puzzle of phase separation in a granular monolayer vibrated from below. Although this system is three dimensional, an interesting dynamics occurs mostly in the horizontal plane, perpendicular to the direction of vibration. Experiments [Olafsen and Urbach, Phys. Rev. Lett. 81, 4369 (1998)] demonstrated that for a high amplitude of vibration the system is in the gaslike phase, but when the amplitude becomes smaller than a certain threshold, a phase separation occurs: A solidlike dense condensate of particles forms in the center of the system, surrounded by particles in the gaslike phase. We explain theoretically the experimentally observed coexistence of dilute and dense phases, employing Navier-Stokes granular hydrodynamics. We show that the phase separation is associated with a negative compressibility of granular gas.

Knudsen temperature jump and the Navier-Stokes hydrodynamics of granular gases driven by thermal walls

Thermal wall is a convenient idealization of a rapidly vibrating plate used for vibrofluidization of granular materials. The objective of this work is to incorporate the Knudsen temperature jump at thermal wall in the Navier-Stokes hydrodynamic modeling of dilute granular gases of monodisperse particles that collide nearly elastically. The Knudsen temperature jump manifests itself as an additional term, proportional to the temperature gradient, in the boundary condition for the temperature. Up to a numerical prefactor O(1) , this term is known from kinetic theory of elastic gases. We determine the previously unknown numerical prefactor by measuring, in a series of molecular dynamics (MD) simulations, steady-state temperature profiles of a gas of elastically colliding hard disks, confined between two thermal walls kept at different temperatures, and comparing the results with the predictions of a hydrodynamic calculation employing the modified boundary condition. The modified boundary condition is then applied, without any adjustable parameters, to a hydrodynamic calculation of the temperature profile of a gas of inelastic hard disks driven by a thermal wall. We find the hydrodynamic prediction to be in very good agreement with MD simulations of the same system. The results of this work pave the way to a more accurate hydrodynamic modeling of driven granular gases.

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.

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.

van der Waals normal form for a one-dimensional hydrodynamic model

Phase separation in a fluidized granular system is studied. We consider a one-dimensional hydrodynamic model that mimics a two-dimensional fluidized granular system with a vibrating wall and without gravity, which exhibits a phase separation. Close to the critical point, by means of an adiabatic elimination of the temperature, we deduce the van der Waals normal form, which is the equation that describes the slow dynamics of the system and predicts the qualitative behavior in different regions of parameters. This allows us to understand the origin of the effective viscosity and the spatial saturation at the onset of the bifurcation. The hydrodynamic model and van der Waals normal form exhibit a behavior similar to the one observed in molecular dynamics simulations.

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.

Collision statistics of driven granular materials

We present an experimental investigation of the statistical properties of spherical granular particles on an inclined plane that are excited by an oscillating side wall. The data is obtained by high-speed imaging and particle tracking techniques. We identify all particles in the system and link their positions to form trajectories over long times. Thus, we identify particle collisions to measure the effective coefficient of restitution and find a broad distribution of values for the same impact angles. We find that the energy inelasticity can take on values greater than one, which implies that the rotational degrees of freedom play an important role in energy transfer. We also measure the distance and the time between collision events in order to directly determine the distribution of path lengths and the free times. These distributions are shown to deviate from expected theoretical forms for elastic spheres, demonstrating the inherent clustering in this system. We describe the data with a two-parameter fitting function and use it to calculate the mean free path and collision time. We find that the ratio of these values is consistent with the average velocity. The velocity distributions are observed to be strongly non-Gaussian and do not demonstrate any apparent universal behavior. We report the scaling of the second moment, which corresponds to the granular temperature, and higher order moments as a function of distance from the driving wall. Additionally, we measure long-time correlation functions in both space and in the velocities to probe diffusion in a dissipative gas.

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.

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

Granular hydrodynamics predicts symmetry-breaking instability in a two-dimensional ensemble of nearly elastically colliding smooth hard disks driven, at zero gravity, by a rapidly vibrating sidewall. Supercritical and subcritical symmetry-breaking bifurcations of the stripe state are identified, and the supercritical bifurcation curve is computed. The cluster dynamics proceed as a coarsening process mediated by the gas phase. Well above the bifurcation point the final steady state, selected by coarsening, represents a single strongly localized densely packed "droplet."