Breakdown of hydrodynamics in a one-dimensional system of inelastic particles

Rescaling invariance and anomalous energy transport in a small vertical column of grains

It is well known that energy dissipation and finite size can deeply affect the dynamics of granular matter, often making usual hydrodynamic approaches problematic. Here we report on the experimental investigation of a small model system, made of ten beads constrained into a 1D geometry by a narrow vertical pipe and shaken at the base by a piston excited by a periodic wave. Recording the beads motion with a high frame rate camera allows to investigate in detail the microscopic dynamics and test hydrodynamic and kinetic models. Varying the energy, we explore different regimes from fully fluidized to the edge of condensation, observing good hydrodynamic behavior down to the edge of fluidization, despite the small system size. Density and temperature fields for different system energies can be collapsed by suitable space and time rescaling, and the expected constitutive equation holds very well when the particle diameter is considered. At the same time, the balance between dissipated and fed energy is not well described by commonly adopted dependence due to the up-down symmetry breaking. Our observations, supported by the measured particle velocity distributions, show a different phenomenological temperature dependence, which yields equation solutions in agreement with experimental results.

Kapitza resistance at a domain boundary in linear and nonlinear chains

We explore Kapitza thermal resistance on the boundary between two homogeneous chain fragments with different characteristics. For a linear model, an exact expression for the resistance is derived. In the generic case of frequency mismatch between the domains, the Kapitza resistance is well defined in the thermodynamic limit. At the same time, in the linear chain, the resistance depends on the thermostat properties and therefore is not a local property of the considered domain boundary. Moreover, if the temperature difference at the ends of the chain is fixed, then neither the temperature drop at the domain boundary nor the heat flux depend on the system size; for the normal transport, one expects the scaling N^{-1} for both. For specific assessment, we consider the case of an isotopic boundary-only the masses in different domains are different. If the domains are nonlinear, but integrable (Toda lattice, elastically colliding particles), the anomalies are similar to the case of linear chain, with the addition of well-articulated thermal dependence of the resistance. For the case of elastically colliding particles, this dependence follows a simple scaling law R_{k}∼T^{-1/2}. For Fermi-Pasta-Ulam domains, both the temperature drop and the heat flux decrease with the chain length, but with different exponents, so the resistance vanishes in the thermodynamic limit. For the domains comprised of rotators, the thermal resistance exhibits the expected normal behavior.

Blast in a One-Dimensional Cold Gas: From Newtonian Dynamics to Hydrodynamics

A gas composed of a large number of atoms evolving according to Newtonian dynamics is often described by continuum hydrodynamics. Proving this rigorously is an outstanding open problem, and precise numerical demonstrations of the equivalence of the hydrodynamic and microscopic descriptions are rare. We test this equivalence in the context of the evolution of a blast wave, a problem that is expected to be at the limit where hydrodynamics could work. We study a one-dimensional gas at rest with instantaneous localized release of energy for which the hydrodynamic Euler equations admit a self-similar scaling solution. Our microscopic model consists of hard point particles with alternating masses, which is a nonintegrable system with strong mixing dynamics. Our extensive microscopic simulations find a remarkable agreement with Euler hydrodynamics, with deviations in a small core region that are understood as arising due to heat conduction.

Speed-dispersion-induced alignment: A one-dimensional model inspired by swimming droplets experiments

We investigate the collective dynamics of self-propelled droplets, confined in a one-dimensional microfluidic channel. On the one hand, neighboring droplets align and form large trains of droplets moving in the same direction. On the other hand, the droplets condensate, leaving large regions with very low density. A careful examination of the interactions between two "colliding" droplets demonstrates that local alignment takes place as a result of the interplay between the dispersion of their speeds and the absence of Galilean invariance. Inspired by these observations, we propose a minimalistic 1D model of active particles reproducing such dynamical rules and, combining analytical arguments and numerical evidences, we show that the model exhibits a transition to collective motion in 1D for a large range of values of the control parameters. Condensation takes place as a transient phenomena, which tremendously slows down the dynamics, before the system eventually settles into a homogeneous aligned phase.

Thermal conductivity at the high-density limit and the levitating granular cluster

The granular Leidenfrost state consists of a dense granular cluster levitating above a hot granular gas. The density of particles inside the cluster can be very high and even close to the density of crystalline packing. To describe this state theoretically, one needs to know the density dependence of constitutive relations (pressure, heat losses, thermal conductivity) at these very high densities. However, the accurate expression for the coefficient of thermal conductivity is lacking. In this work, the constitutive relations were measured at high densities in molecular dynamics simulations in three different settings: a uniform freely cooling dense granulate (to measure heat losses), a uniform ensemble of elastically colliding particles (to measure pressure), and a dense granular medium between two thermal walls under gravity (to measure thermal conductivity). Next, the hydrodynamic equations with the resulting expressions were solved to describe the levitating cluster state in various parameter regimes. Separate molecular dynamics simulations were performed to test the theoretical predictions and measure the density and temperature profiles of the granular Leidenfrost state, and a good agreement with theoretical results was observed.

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.

Experiments and characterization of low-frequency oscillations in a granular column

The behavior of a vertically vibrated granular bed is reminiscent of a liquid in that it exhibits many phenomena such as convection and Faraday-like surface waves. However, when the lateral dimensions of the bed are confined such that a quasi-one-dimensional geometry is formed, the only phenomena that remain are bouncing bed and the granular Leidenfrost effect. This permits the observation of the granular Leidenfrost state for a wide range of energy injection parameters and more specifically allows for a thorough characterization of the low-frequency oscillation (LFO) that is present in this state. In both experiments and particle simulations we determine the LFO frequency from the power spectral density of the center-of-mass signal of the grains, varying the amplitude and frequency of the driving, the particle diameter, and the number of layers in the system. We thus find that the LFO frequency (i) is inversely proportional to the fast inertial timescale and (ii) decorrelates with a typical decay time proportional to the slow dissipative timescale in the system. The latter is consistent with the view that the LFO is driven by the inherent noise that is present in the granular Leidenfrost state with a low number of particles.

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.

Burrowing behaviour of robotic bivalves with synthetic morphologies

Several bivalve species burrow into sandy sediments to reach their living position. There are many hypotheses concerning the functional morphology of the bivalve shell for burrowing. Observational studies are limited and often qualitative and should be complemented by a synthetic approach mimicking the burrowing process using a robotic emulation. In this paper we present a simple mechatronic set-up to mimic the burrowing behaviour of bivalves. As environment we used water and quartz sand contained in a glass tank. Bivalve shells were mathematically modelled on the computer and then materialized using a 3D printer. The burrowing motion of the shells was induced by two external linear motors. Preliminary experiments did not expose any artefacts introduced to the burrowing process by the set-up. We tested effects of shell size, shape and surface sculpturing on the burrowing performance. Neither the typical bivalve shape nor surface sculpture did have a clear positive effect on burrowing depth in the performed experiments. We argue that the presented method is a valid and promising approach to investigate the functional morphology of bivalve shells and should be improved and extended in future studies. In contrast to the observation of living bivalves, our approach offers complete control over the parameters defining shell morphology and motion pattern. The technical set-up allows the systematic variation of all parameters to quantify their effects. The major drawback of the built set-up was that the reliability and significance of the results was limited by the lack of an optimal technique to standardize the sediment state before experiments.

Periodic orbits of inelastic particles on a ring

We consider the dynamics of N rigid particles of arbitrary mass that are constrained to move on a frictionless ring. Collisions between particles are inelastic with a constant coefficient of restitution e, and between collisions the particles move with constant velocity. We study sequences of collisions that are self-similar in the sense that the relative positions return to their original relative positions after the collision sequence while the relative velocities are reduced by a constant factor. For a given collision sequence, we develop the analytic machinery to determine the particle velocities and the locations of collisions, and we show that the problem of determining self-similar orbits reduces to solving an eigenvalue problem to obtain the particle velocities and solving a linear system to obtain the locations of interparticle collisions. For inelastic systems, we show that the collision locations can always be uniquely determined. We also show that this is in sharp contrast to the case of elastic systems in which infinite families of self-similar orbits can coexist.

Effective particles and classification of the dynamics of homogeneous granular chains with no precompression

We develop a systematic methodology for classifying the periodic orbits of homogeneous ordered granular chains with no dissipation, under the assumption that all granules oscillate with the same frequency. The analysis is based on the idea of balancing linear momentum for sets of auxiliary models consisting of "effective particles." The auxiliary models may be defined for any given finite, ordered granular chain composed of n identical granules (beads) that interact with each other through strongly nonlinear Hertzian interaction law. In turn, the auxiliary models may be effectively used for theoretically predicting the total number of periodic orbits and the corresponding amplitude ratios of the granules. Good correspondence between the theoretical models and results of direct numerical simulations is reported. The results presented herein can be used to understand the complex intrinsic dynamics of ordered granular media, and to systematically study the generation of mode localization in these strongly nonlinear systems. The derived analytical models can be utilized to predict the response of the effective particles, and based on that, to predict primary pulse transmission in periodic layered media with granular interfaces. Moreover, our analysis can be extended to the general class of nonlinear chains of particles with smooth interacting potentials and possible separation between particles during the motion.

Sudden chain energy transfer events in vibrated granular media

In a mixture of two species of grains of equal size but different mass, placed in a vertically vibrated shallow box, there is spontaneous segregation. Once the system is at least partly segregated and clusters of the heavy particles have formed, there are sudden peaks of the horizontal kinetic energy of the heavy particles, that is otherwise small. Together with the energy peaks the clusters rapidly expand and the segregation is partially lost. The process repeats once segregation has taken place again, either randomly or with some regularity in time depending on the experimental or numerical parameters. An explanation for these events is provided based on the existence of a fixed point for an isolated particle bouncing with only vertical motion. The horizontal energy peaks occur when the energy stored in the vertical motion is partly transferred into horizontal energy through a chain reaction of collisions between heavy particles.

Degenerate orbit transitions in a one-dimensional inelastic particle system

Continuous transitions between different periodic orbits in a one-dimensional inelastic particle system are investigated. We show that continuous transitions that occur when adding or subtracting a single collision are, generically, of co-dimension 2. We give a full mechanical description of the system and explain why this is the case. Surprisingly, we also show that there are an infinite set of degenerate transitions of co-dimension 1. We provide a theoretical analysis that gives a simple criteria to classify which transitions are degenerate purely using the discrete set of collisions that occur in the orbits. Our analysis allows us to understand the nature of the degeneracy. We also show that higher degrees of degeneracy can occur, and provide an explanation.

Onset of convection in strongly shaken granular matter

Strongly vertically shaken granular matter can display a density inversion: A high-density cluster of beads is elevated by a dilute gaslike layer of fast beads underneath ("granular Leidenfrost effect"). For even stronger shaking the granular Leidenfrost state becomes unstable and granular convection rolls emerge. This transition resembles the classical onset of convection in fluid heated from below at some critical Rayleigh number. The same transition is seen in molecular dynamics (MD) simulations of the shaken granular material. The critical shaking strength for the onset of granular convection can be calculated from a linear stability analysis of a hydrodynamiclike model of the granular flow. Experiment, MD simulations, and theory quantitatively agree.

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.

Bioconvection and front formation of Paramecium tetraurelia

We have investigated the bioconvection of Paramecium tetraurelia in high-density suspensions made by centrifugal concentration. When a suspension is kept at rest in a Hele-Shaw cell, a crowded front of paramecia is formed in the vicinity of the bottom and it propagates gradually toward the water-air interface. Fluid convection occurs under this front, and it is driven persistently by the upward swimming of paramecia. The roll structures of the bioconvection become turbulent with an increase in the depth of the suspension; they also change rapidly as the density of paramecia increases. Our experimental results suggest that lack of oxygen in the suspension causes the active individual motions of paramecia to induce the formation of this front.

Coefficient of restitution for viscoelastic disks

The dissipative collision of two identical viscoelastic disks is studied. By using a known law for the elastic part of the interaction force and the viscoelastic damping model an analytical solution for the coefficient of restitution is given. The coefficient of restitution depends significantly on the impact velocity. It approaches 1 for small velocities and decreases for increasing velocities.

Theory of thermostatted inhomogeneous granular fluids: A self-consistent density functional description

The authors present a study of the nonequilibrium statistical properties of a one dimensional hard-rod fluid dissipating energy via inelastic collisions and subject to the action of a Gaussian heat bath, simulating an external driving mechanism. They show that the description of the fluid based on the one-particle phase-space reduced distribution function, in principle necessary because of the presence of velocity dependent collisional dissipation, can be contracted to a simpler description in configurational space. Indeed, by means of a multiple-time-scale method the authors derive a self-consistent governing equation for the particle density distribution function. This equation is similar to the dynamic density functional equation employed in the study of colloids, but contains additional terms taking into account the inelastic nature of the fluid. Such terms cannot be derived from a Liapunov generating functional and contribute not only to the relaxational properties, but also to the nonequilibrium steady state properties. A validation of the theory against molecular dynamics simulations is presented in a series of cases, and good agreement is found.

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.

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.

Transport of a heated granular gas in a washboard potential

We study numerically the motion of a one dimensional array of Brownian particles in a washboard potential, driven by an external stochastic force and interacting via short range repulsive forces. In particular, we investigate the role of instantaneous elastic and inelastic collisions on the system dynamics and transport. The system displays a locked regime, where particles may move only via activated processes and a running regime where particles drift along the direction of the applied field. By tuning the value of the friction parameter controlling the Brownian motion we explore both the overdamped dynamics and the underdamped dynamics. In the two regimes we considered the mobility and the diffusivity of the system as functions of the tilt and other relevant control parameters such as coefficient of restitution, particle size, and total number of particles. We find that while in the overdamped regime the results for the interacting systems present similarities with the known noninteracting case, in the underdamped regime the inelastic collisions determine a rich variety of behaviors among which is an unexpected enhancement of the inelastic diffusion.

Collapse of periodic orbits in a driven inelastic particle system

The dynamical behavior of a one-dimensional inelastic particle system with particles of unequal mass traveling between two walls is investigated. The system is driven by adding energy at one of the walls while the other wall is stationary and does not add energy. By deriving analytic solutions for the periodic orbits of this system, we show that there are a countable infinity of critical mass ratios at which the particle dynamics become highly degenerate in the following sense. As the mass ratio passes through these critical points, large numbers of stable periodic orbits can collapse onto a single trivial orbit. We show that the widely studied equal-mass systems represent one of these critical points and are therefore such a degenerate case. We also show that in the elastic limit the number of orbits that collapse onto the single trivial orbit can become arbitrarily large.

Granular fountains: convection cascade in a compartmentalized granular gas

This paper extends the two-compartment granular fountain [D. van der Meer, P. Reimann, K. van der Weele, and D. Lohse, Phys. Rev. Lett. 92, 184301 (2004)] to an arbitrary number of compartments: the tendency of a granular gas to form clusters is exploited to generate spontaneous convective currents, with particles going down in the well-filled compartments and going up in the diluted ones. We focus upon the bifurcation diagram of the general -compartment system, which is constructed using a dynamical flux model and which proves to agree quantitatively with results from molecular dynamics simulations.

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.

Inelastic Takahashi hard-rod gas

We study a one-dimensional fluid of hard rods interacting with each other via binary inelastic collisions and a short-ranged square-well potential. Upon tuning the depth and the sign of the well, we investigate the interplay between dissipation and cohesive or repulsive forces. Molecular-dynamics simulations of the cooling regime indicate that the presence of this simple interparticle interaction is sufficient to significantly modify the energy dissipation rates expected by Haff's law for the free cooling. The simplicity of the model makes it amenable to an analytical approach based on the Boltzmann-Enskog transport equation which allows deriving the behavior of the granular temperature. Furthermore, in the elastic limit, the model can be solved exactly to provide a full thermodynamic description. A meaningful theoretical approximation explaining the properties of the inelastic system in interaction with a thermal bath can be directly extrapolated from the properties of the corresponding elastic system, upon a proper redefinition of the relevant observables. Simulation results both in the cooling and driven regimes can be fairly interpreted according to our theoretical approach and compare rather well to our predictions.

Granular Leidenfrost effect: experiment and theory of floating particle clusters

Granular material is vertically vibrated in a 2D container: above a critical shaking strength, and for a sufficient number of beads, a crystalline cluster is elevated and supported by a dilute gaseous layer of fast beads underneath. We call this phenomenon the granular Leidenfrost effect. The experimental observations are explained by a hydrodynamic model featuring three dimensionless control parameters: the energy input S, the number of particle layers F, and the inelasticity of the particle collisions epsilon. The (S,F) phase diagram, in which the Leidenfrost state lies between the purely solid and gas phases, shows accurate agreement between experiment and theory.

Velocity distributions in dilute granular systems

We investigate the idea that velocity distributions in granular gases are determined mainly by eta, the coefficient of restitution and q, which measures the relative importance of heating (or energy input) to collisions. To this end, we study by numerical simulation the properties of inelastic gases as functions of eta, concentration phi, and particle number N with various heating mechanisms. For a wide range of parameters, we find Gaussian velocity distributions for uniform heating and non-Gaussian velocity distributions for boundary heating. Comparison between these results and velocity distributions obtained by other heating mechanisms and for a simple model of a granular gas without spatial degrees of freedom, shows that uniform and boundary heating can be understood as different limits of q, with q>>1 and q < or approximately 1 respectively. We review the literature for evidence of the role of q in the recent experiments.

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.

Structural transitions in two-dimensional hard-sphere systems

We spread randomly noncharged steel particles (diameter, 1.59 mm) on a silicon wafer to form a two-dimensional hard-sphere system. The particle structure versus the particle coverage was monitored. We observed the particle structural transition from liquidlike to triangular-lattice crystal-like with increasing particle coverage by analyzing the particle structure factor. The particle coverage at which the structural transition occurs was quantified by the curves of S(max) (A) and G6 (A); S(max) is the amplitude of the first peak of the structure factor (depicting the particle positional order), and G6 is the bond orientation order parameter. We also conducted a Monte Carlo simulation study. The Monte Carlo simulation results show good agreement with the experimental results at low particle area fractions. However, at high area fractions, the experimentally observed particle structure is less organized than that generated by simulations.

Are Brazil nuts attractive?

We present event-driven simulation results for single and multiple intruders in a vertically vibrated granular bed. Under our vibratory conditions, the mean vertical position of a single intruder is governed primarily by a buoyancylike effect. Multiple intruders also exhibit buoyancy governed behavior; however, multiple neutrally buoyant intruders cluster spontaneously and undergo horizontal segregation. These effects can be understood by considering the dynamics of two neutrally buoyant intruders. We have measured an attractive force between such intruders which has a range of five intruder diameters, and we provide a mechanistic explanation for the origins of this force.

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.

Crucial role of sidewalls in velocity distributions in quasi-two-dimensional granular gases

Our experiments and three-dimensional molecular dynamics simulations of particles confined to a vertical monolayer by closely spaced frictional walls (sidewalls) yield velocity distributions with non-Gaussian tails and a peak near zero velocity. Simulations with frictionless sidewalls are not peaked. Thus interactions between particles and their containers are an important determinant of the shape of the distribution and should be considered when evaluating experiments on a constrained monolayer of particles.

Inelastic hard rods in a periodic potential

A simple model of inelastic hard rods subject to a one-dimensional array of identical wells is introduced. The energy loss due to inelastic collisions is balanced by the work supplied by an external stochastic heat bath. We explore the effect of the spatial nonuniformity on the steady states of the system. The spatial variations of the density, granular temperature, and pressure induced by the gradient of the external potential are investigated and compared with the analogous variations in an elastic system. Finally, we study the clustering process by considering the relaxation of the system starting from a uniform homogeneous state.

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.

Pressure measurement in two-dimensional horizontal granular gases

A two-dimensional granular gas is produced by vibrating vertically a partial layer of beads on a horizontal plate. Measurements of the force applied by the granular gas to the sidewalls of the container, or granular pressure, are used to study the effect of the shaking strength, density, bead-plate restitution coefficient, and particle size on the steady properties of the gas.

Velocity distributions in dissipative granular gases

Motivated by recent experiments reporting non-Gaussian velocity distributions in driven dilute granular materials, we study by numerical simulation the properties of 2D inelastic gases. We find theoretically that the form of the observed velocity distribution is governed primarily by the coefficient of restitution eta and q=N(H)/N(C), the ratio between the average number of heatings and the average number of collisions in the gas. The differences in distributions we find between uniform and boundary heating can then be understood as different limits of q, for q>>1 and q less, similar 1, respectively.

Fluid-like behavior of a one-dimensional granular gas

We study the properties of a one-dimensional (1D) granular gas consisting of N hard rods on a line of length L (with periodic boundary conditions). The particles collide inelastically and are fluidized by a heat bath at temperature Tb and viscosity gamma. The analysis is supported by molecular dynamics simulations. The average properties of the system are first discussed, focusing on the relations between granular temperature Tg=mv2, kinetic pressure, and density rho=N/L. Thereafter, we consider the fluctuations around the average behavior obtaining a slightly non-Gaussian behavior of the velocity distributions and a spatially correlated velocity field; the density field displays clustering: this is reflected in the structure factor which has a peak in the k approximately 0 region suggesting an analogy between inelastic hard core interactions and an effective attractive potential. Finally, we study the transport properties, showing the typical subdiffusive behavior of 1D stochastically driven systems, i.e., </x(t)-x(0)/2> approximately Dt(1/2), where D for the inelastic fluid is larger than the elastic case. This is directly related to the peak of the structure factor at small wave vectors.

Speckle visibility spectroscopy and variable granular fluidization

We introduce a dynamic light scattering technique capable of resolving motion that changes systematically, and rapidly, with time. It is based on the visibility of a speckle pattern for a given exposure duration. Applying this to a vibrated layer of glass beads, we measure the granular temperature and its variation with phase in the oscillation cycle. We observe several transitions involving jammed states, where the grains are at rest during some portion of the cycle. We also observe a two-step decay of the temperature on approach to jamming.

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."

Symmetry-breaking instability in a prototypical driven granular gas

Symmetry-breaking instability of a laterally uniform granular cluster (strip state) in a prototypical driven granular gas is investigated. The system consists of smooth hard disks in a two-dimensional box, colliding inelastically with each other and driven, at zero gravity, by a "thermal" wall. The limit of nearly elastic particle collisions is considered, and granular hydrodynamics with the Jenkins-Richman constitutive relations is employed. The hydrodynamic problem is completely described by two scaled parameters and the aspect ratio of the box. Marginal stability analysis predicts a spontaneous symmetry-breaking instability of the strip state, similar to that predicted recently for a different set of constitutive relations. If the system is big enough, the marginal stability curve becomes independent of the details of the boundary condition at the driving wall. In this regime, the density perturbation is exponentially localized at the elastic wall opposite the thermal wall. The short- and long-wavelength asymptotics of the marginal stability curves are obtained analytically in the dilute limit. The physics of the symmetry-breaking instability is discussed.

Scaling, multiscaling, and nontrivial exponents in inelastic collision processes

We investigate velocity statistics of homogeneous inelastic gases using the Boltzmann equation. Employing an approximate uniform collision rate, we obtain analytic results valid in arbitrary dimension. In the freely evolving case, the velocity distribution is characterized by an algebraic large-velocity tail, P(v,t) approximately v(-sigma). The exponent sigma(d,epsilon), a nontrivial root of an integral equation, varies continuously with the spatial dimension d and the dissipation coefficient epsilon. Although the velocity distribution follows a scaling form, its moments exhibit multiscaling asymptotic behavior. Furthermore, the velocity autocorrelation function decays algebraically with time, A(t)=<v(0).v(t)> approximately t(-alpha), with a nonuniversal dissipation-dependent exponent alpha=1/epsilon. In the forced case, the steady state Fourier transform is obtained via a cumulant expansion. Even in this case, velocity correlations develop and the velocity distribution is non-Maxwellian.

Sudden collapse of a granular cluster

Single clusters in a vibro-fluidized granular gas in N connected compartments become unstable at strong shaking. They are experimentally shown to collapse very abruptly. The observed cluster lifetime (as a function of the driving intensity) is analytically calculated within a flux model, making use of the self-similarity of the process. After collapse, the cluster diffuses out into the uniform distribution in a self-similar way, with an anomalous diffusion exponent 1/3.

Breaking of equipartition in one-dimensional heat-conducting systems

Using information-theoretical methods, we studied how energy equipartition is broken in one-dimensional systems under a heat flow composed of alternating particles of two different masses. The average energy stored in particles of different masses is seen to be different in both ideal gases and harmonic lattices.

Characterization of the stationary states of a dilute vibrofluidized granular bed

This paper reports two phenomena in an event driven simulation of a dilute vibrofluidized granular material in two dimensions. Both phenomena show inhomogeneity in the horizontal direction. They are convection rolls similar to the Rayleigh-Bénard thermal convection in fluids, and a clustering instability, where the bed spontaneously phase separates into coexisting dense and dilute regions. Detailed investigations show that these are different from the known instabilities in a vibrated granular medium. Characterization of these instabilities is carried out with a phase diagram using suitable parameters from the kinetic theory of vibrofluidized beds.

Bifurcations of a driven granular system under gravity

The molecular dynamics study on the granular bifurcation in a simple model is presented. The model consists of hard disks, which undergo inelastic collisions; the system is under the uniform external gravity and is driven by the heat bath. The competition between the two effects, namely, the gravitational force and the heat bath, is carefully studied. We found that the system shows three phases, namely, the condensed phase, the locally fluidized phase, and the granular turbulent phase, upon increasing the external control parameter. We conclude that the transition from the condensed phase to the locally fluidized phase is distinguished by the existence of fluidized holes, and the transition from the locally fluidized phase to the granular turbulent phase is understood by the destabilization transition of the fluidized holes due to mutual interference.