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

The Bernoulli effect in horizontal granular flows

The Bernoulli effect that commonly occurs in continuous fluids is simultaneously increasing a fluid's velocity and decreasing static pressure or the fluid's gravitational potential energy. Although, the Bernoulli effect has already been extensively explored, there is a lack of research on the relationship between flow velocity and pressure in a discrete medium. In the present study, this relationship in horizontal granular flows excited by vertical vibration is experimentally studied. It was found that the random motion and horizontal directed motion of the granules restrict each other so that the total pressure remains almost constant with respect to time and height. In fact, it implies that the Bernoulli effect occurs in the granular flows. It was also found that the pressure constant of the Bernoulli effect depends on the vibrating intensity and frequency, which reflects the energy transfer in the granular flows. Our results show a dynamic property of the granular flows, which is different from continuous fluids, even though it is similar to some extent.

Reticulate collisional structure in boundary-driven granular gases

We report a peculiar head-on collision network between two vibrating boundaries in experiments performed during a parabolic flight and in a laboratory using horizontal vibration. This structure is a new ordering, which is due to an orientation correlation between the relative position and velocity of any particle pair. It weakens the collision frequency and produces a long-range boundary effect. Moreover, we find the molecular chaos assumption is violated in a larger portion of the phase space. Using an anisotropic distribution model, we modify angular integration results and compare them to the results of the kinetic theory.

Hydrodynamics of granular gases with a two-peak distribution

Vibrating walls, frequently employed to maintain the temperature (i.e., average velocity) in a granular gas, modify the system strongly, rendering it dissimilar to a molecular one in various aspects. As evidenced by microgravity experiments employing a quasi-two-dimensional (quasi-2D) rectangular box and by 2D simulations, the one-peak velocity distribution is split into two, rendering the stress both nonuniform and anisotropic-without a shear flow and in the absence of gravitation. To account for this, granular hydrodynamics (as first proposed by Haff and later derived employing the kinetic theory) is generalized by introducing two additional variables, with one accounting for the distance between the two peaks and a second for the difference between the average velocities along different directions. The hydrodynamic theory thus generalized relates the velocity distribution to the stress, yielding results that agree with experiments and simulations.

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.

Hydrodynamic correlation functions of a driven granular fluid in steady state

We study a homogeneously driven granular fluid of hard spheres at intermediate volume fractions and focus on time-delayed correlation functions in the stationary state. Inelastic collisions are modeled by incomplete normal restitution, allowing for efficient simulations with an event-driven algorithm. The incoherent scattering function F(incoh)(q,t) is seen to follow time-density superposition with a relaxation time that increases significantly as the volume fraction increases. The statistics of particle displacements is approximately Gaussian. For the coherent scattering function S(q,ω), we compare our results to the predictions of generalized fluctuating hydrodynamics, which takes into account that temperature fluctuations decay either diffusively or with a finite relaxation rate, depending on wave number and inelasticity. For sufficiently small wave number q we observe sound waves in the coherent scattering function S(q,ω) and the longitudinal current correlation function C(l)(q,ω). We determine the speed of sound and the transport coefficients and compare them to the results of kinetic theory.

Model for density waves in gravity-driven granular flow in narrow pipes

A gravity-driven flow of grains through a narrow pipe in vacuum is studied by means of a one-dimensional model with two coefficients of restitution. Numerical simulations show clearly how density waves form when a strikingly simple criterion is fulfilled: that dissipation due to collisions between the grains and the walls of the pipe is greater per collision than that which stems from collisions between particles. Counterintuitively, the highest flow rate is observed when the number of grains per density wave grows large. We find strong indication that the number of grains per density wave always approaches a constant as the particle number tends to infinity, and that collapse to a single wave, which was often observed also in previous simulations, occurs because the number of grains is insufficient for multiple wave formation.

Traveling shock front in quasi-two-dimensional granular flows

While the density profile of a granular shock front can be obtained by the conventional treatment of supersonic fluids, its temperature profile is very different from that in ordinary shocks. We study the density and temperature profiles of a traveling granular shock generated by piling up metal spheres in a closed bottom quasi-two-dimensional channel. We successfully account for the temperature profile in the granular shock using a simple kinetic theory in terms of energy transfer from the mean flow direction to the transverse direction. Contrary to ordinary fluids and previous granular shock experiments, the granular shock width is found to increase with the inflow rate.

Particle kinematics in a dilute, three-dimensional, vibration-fluidized granular medium

We report an experimental study of particle kinematics in a three-dimensional system of inelastic spheres fluidized by intense vibration. The motion of particles in the interior of the medium is tracked by high-speed video imaging, yielding a spatially resolved measurement of the velocity distribution. The distribution is wider than a Gaussian and broadens continuously with increasing volume fraction. The deviations from a Gaussian distribution for this boundary-driven system are different in sign and larger in magnitude than predictions for homogeneously driven systems. We also find correlations between velocity components which grow with increasing volume fraction.

Collision statistics in sheared inelastic hard spheres

The dynamics of sheared inelastic-hard-sphere systems is studied using nonequilibrium molecular-dynamics simulations and direct simulation Monte Carlo. In the molecular-dynamics simulations Lees-Edwards boundary conditions are used to impose the shear. The dimensions of the simulation box are chosen to ensure that the systems are homogeneous and that the shear is applied uniformly. Various system properties are monitored, including the one-particle velocity distribution, granular temperature, stress tensor, collision rates, and time between collisions. The one-particle velocity distribution is found to agree reasonably well with an anisotropic Gaussian distribution, with only a slight overpopulation of the high-velocity tails. The velocity distribution is strongly anisotropic, especially at lower densities and lower values of the coefficient of restitution, with the largest variance in the direction of shear. The density dependence of the compressibility factor of the sheared inelastic-hard-sphere system is quite similar to that of elastic-hard-sphere fluids. As the systems become more inelastic, the glancing collisions begin to dominate over more direct, head-on collisions. Examination of the distribution of the times between collisions indicates that the collisions experienced by the particles are strongly correlated in the highly inelastic systems. A comparison of the simulation data is made with direct Monte Carlo simulation of the Enskog equation. Results of the kinetic model of Montanero [J. Fluid Mech. 389, 391 (1999)] based on the Enskog equation are also included. In general, good agreement is found for high-density, weakly inelastic systems.

Shock front width and structure in supersonic granular flows

The full structure of a shock front around a blunt body in a quasi-two-dimensional granular flow is studied. Two features, a large density gradient and a very small thickness of the front, characterize this shock and make it different from shocks in molecular gases. Both of these features can be understood using a modified version of the granular kinetic theory. Our model separates the particles into two subpopulations: fast particles having experienced no collisions and randomly moving particles. This separation is motivated by direct measurements of the particle velocities which show a bimodal distribution. Our results not only shed new light on the use of the granular kinetic theory under extreme conditions (shock formation) but bring new insight into the physics of shocks in general.

Boltzmann equation for dissipative gases in homogeneous states with nonlinear friction

Combining analytical and numerical methods, we study within the framework of the homogeneous nonlinear Boltzmann equation a broad class of models relevant for the dynamics of dissipative fluids, including granular gases. We use the method presented in a previous paper [J. Stat. Phys. 124, 549 (2006)] and extend our results to a different heating mechanism: namely, a deterministic nonlinear friction force. We derive analytically the high-energy tail of the velocity distribution and compare the theoretical predictions with high-precision numerical simulations. Stretched exponential forms are obtained when the nonequilibrium steady state is stable. We derive subleading corrections and emphasize their relevance. In marginal stability cases, power-law behaviors arise, with exponents obtained as the roots of transcendental equations. We also consider some simple Bhatnagar-Gross-Krook models, driven by similar heating devices, to test the robustness of our predictions.

Chaoticity of the wet granular gas

In this work we derive an analytic expression for the Kolmogorov-Sinai entropy of dilute wet granular matter, valid for any spatial dimension. The grains are modeled as hard spheres and the influence of the wetting liquid is described according to the capillary model, in which dissipation is due to the hysteretic cohesion force of capillary bridges. The Kolmogorov-Sinai entropy is expanded in a series with respect to density. We find a rapid increase of the leading term when liquid is added. This demonstrates the sensitivity of the granular dynamics to humidity, and shows that the liquid significantly increases the chaoticity of the granular gas.

Shear-induced crystallization in jammed systems

Simulations are used to address the effects of oscillating shear strain on jammed systems composed of spherical particles. The simulations show that shear oscillations with amplitudes of more than a few percent lead to substantial crystallization of the system. To ensure that the conclusions are independent of the simulation methodology, a range of simulations are carried out that use both molecular dynamics and athermal dynamics methods, soft and hard potentials, potentials with and without attractive forces, and systems with and without surrounding walls. The extent of crystallization is monitored primarily by the Q(6) order parameter, but also in some simulations by the potential energy and the radial distribution function, and by direct visual inspection. A mechanism is proposed for shear-induced crystallization of jammed systems, based on fold catastrophes of the free-energy landscape.

Correlation and response in a driven dissipative model

We consider a simple dissipative system with spatial structure in contact with a heat bath. The system always exhibits correlations except in the cases of zero and maximal dissipation. We explicitly calculate the correlation function and the nonlocal response function of the system and show that they have the same spatial dependence.

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.

Local heat flux and energy loss in a two-dimensional vibrated granular gas

We performed event-driven simulations of a two-dimensional granular gas between two vibrating walls and directly measured the local heat flux and local energy dissipation in the stationary state. Describing the local heat flux as a function of the coordinate in the direction perpendicular to the driving walls, we test a generalization of Fourier's law, q(x)=-kappa inverted delta T(x)+mu inverted delta rho(x), by relating the local heat flux to the local gradients of the temperature and density. This ansatz accounts for the fact that heat flux can also be generated by density gradients, not only by temperature gradients. Assuming the transport coefficients kappa and mu to be independent of x, we check the validity of this assumption and test the generalized Fourier law in the simulations. Both kappa and mu are determined for different system parameters, in particular, for a wide range of coefficients of restitution. We also compare our numerical results to existing hydrodynamic theories. Agreement is found for kappa for very small inelasticities only, i.e., when the gradients are small. Beyond this region, kappa and mu exhibit a striking nonmonotonic behavior. This may hint that hydrodynamics to Navier-Stokes order cannot be applied to moderately inelastic vibrated systems.