Simulations of vibrated granular medium with impact-velocity-dependent restitution coefficient

Dynamics of a vibration-driven single disk

Granular particles exhibit rich collective behaviors on vibration beds, but the motion of an isolated particle is not well understood even for uniform particles with a simple shape such as disks or spheres. Here we measured the motion of a single disk confined to a quasi-two-dimensional horizontal box on a vertically vibrating stage. The translational displacements obey compressed exponential distributions whose exponent [Formula: see text] increases with the frequency, while the rotational displacements exhibit unimodal distributions at low frequencies and bimodal distributions at high frequencies. During short time intervals, the translational displacements are subdiffusive and negatively correlated, while the rotational displacements are superdiffusive and positively correlated. After prolonged periods, the rotational displacements become diffusive and their correlations decay to zero. Both the rotational and the translational displacements exhibit white noise at low frequencies, and blue noise for translational motions and Brownian noise for rotational motions at high frequencies. The translational kinetic energy obeys Boltzmann distribution while the rotational kinetic energy deviates from it. Most energy is distributed in translational motions at low frequencies and in rotational motions at high frequencies, which violates the equipartition theorem. Translational and rotational motions are not correlated. These experimental results show that the random diffusion of such driven particles is distinct from thermal motion in both the translational and rotational degrees of freedom, which poses new challenges to theory. The results cast new light on the motion of individual particles and the collective motion of driven granular particles.

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.

Controlled Viscosity in Dense Granular Materials

We experimentally investigate the fluidization of a granular material subject to mechanical vibrations by monitoring the angular velocity of a vane suspended in the medium and driven by an external motor. On increasing the frequency, we observe a reentrant transition, as a jammed system first enters a fluidized state, where the vane rotates with high constant velocity, and then returns to a frictional state, where the vane velocity is much lower. While the fluidization frequency is material independent, the viscosity recovery frequency shows a clear dependence on the material that we rationalize by relating this frequency to the balance between dissipative and inertial forces in the system. Molecular dynamics simulations well reproduce the experimental data, confirming the suggested theoretical picture.

Incorporation of velocity-dependent restitution coefficient and particle surface friction into kinetic theory for modeling granular flow cooling

Kinetic theory (KT) has been successfully used to model rapid granular flows in which particle interactions are frictionless and near elastic. However, it fails when particle interactions become frictional and inelastic. For example, the KT is not able to accurately predict the free cooling process of a vibrated granular medium that consists of inelastic frictional particles under microgravity. The main reason that the classical KT fails to model these flows is due to its inability to account for the particle surface friction and its inelastic behavior, which are the two most important factors that need be considered in modeling collisional granular flows. In this study, we have modified the KT model that is able to incorporate these two factors. The inelasticity of a particle is considered by establishing a velocity-dependent expression for the restitution coefficient based on many experimental studies found in the literature, and the particle friction effect is included by using a tangential restitution coefficient that is related to the particle friction coefficient. Theoretical predictions of the free cooling process by the classical KT and the improved KT are compared with the experimental results from a study conducted on an airplane undergoing parabolic flights without the influence of gravity [Y. Grasselli, G. Bossis, and G. Goutallier, Europhys. Lett. 86, 60007 (2009)10.1209/0295-5075/86/60007]. Our results show that both the velocity-dependent restitution coefficient and the particle surface friction are important in predicting the free cooling process of granular flows; the modified KT model that integrates these two factors is able to improve the simulation results and leads to better agreement with the experimental results.

Maximizing energy transfer in vibrofluidized granular systems

Using discrete particle simulations validated by experimental data acquired using the positron emission particle tracking technique, we study the efficiency of energy transfer from a vibrating wall to a system of discrete, macroscopic particles. We demonstrate that even for a fixed input energy from the wall, energy conveyed to the granular system under excitation may vary significantly dependent on the frequency and amplitude of the driving oscillations. We investigate the manner in which the efficiency with which energy is transferred to the system depends on the system variables and determine the key control parameters governing the optimization of this energy transfer. A mechanism capable of explaining our results is proposed, and the implications of our findings in the research field of granular dynamics as well as their possible utilization in industrial applications are discussed.

Low-frequency oscillations and convective phenomena in a density-inverted vibrofluidized granular system

Low-frequency oscillations (LFOs) are thought to play an important role in the transition between the Leidenfrost and convective states of a vibrated granular bed. This work details the experimental observation of LFOs, which are found to be consistently present for a range of driving frequencies and amplitudes, with particles of varying material and using containers of differing material properties. The experimentally acquired results show a close qualitative and quantitative agreement with both theory and simulations across the range of parameters tested. Strong agreement between experimental and simulation results was also observed when investigating the influence of sidewall dissipation on LFOs and vertical density profiles. This paper additionally provides evidence of two phenomena present in the Leidenfrost state: a circulatory motion over extended time periods in near-crystalline configurations, and a Leidenfrost-like state in which the dense upper region displays an unusual inverse thermal convection.

Simulations of granular gravitational collapse

A freely cooling granular gas in a gravitational field undergoes a collapse to a multicontact state in a finite time. Previous theoretical [D. Volfson et al., Phys. Rev. E 73, 061305 (2006)] and experimental work [R. Son et al., Phys. Rev. E 78, 041302 (2008)] have obtained contradictory results about the rate of energy loss before the gravitational collapse. Here we use a molecular dynamics simulation in an attempt to recreate the experimental and theoretical results to resolve the discrepancy. We are able to nearly match the experimental results, and find that to reproduce the power law predicted in the theory we need a nearly elastic system with a constant coefficient of restitution greater than 0.993. For the more realistic velocity-dependent coefficient of restitution, there does not appear to be a power-law decay and the transition from granular gas to granular solid is smooth, making it difficult to define a time of collapse.

Phase transitions in vibrated granular systems in microgravity

We numerically investigated various dynamical behaviors of a vibrated granular gas in microgravity. Using the parameters of an earlier Mini-Texus 5 experiment, three-dimensional simulations, based on molecular dynamics, efficiently reproduce experimental results. Using Kolmogorov-Smirnov tests, four dynamical regimes have been distinguished: gaseous state, partial clustering, complete clustering, and bouncing aggregates. Different grain radii and densities have been considered in order to describe a complete (r,η)-phase diagram. The latter exhibits rich features such as phase transitions and triple points. Our work emphasizes the complexity of diluted granular systems and opens fundamental perspectives.

Energy transport in a one-dimensional granular gas

We study heat conduction in one-dimensional granular gases. In particular, we consider two mechanisms of viscous dissipation during intergrain collisions. In one, the dissipative force is proportional to the grain's velocity and dissipates not only energy but also momentum. In the other, the dissipative force is proportional to the relative velocity of the grains and therefore conserves momentum even while dissipating energy. This allows us to explore the role of momentum conservation in the heat conduction properties of this one-dimensional nonlinear system. We find normal thermal conduction whether or not momentum is conserved.

Dissipative solitary waves in granular crystals

We provide a quantitative characterization of dissipative effects in one-dimensional granular crystals. We use the propagation of highly nonlinear solitary waves as a diagnostic tool and develop optimization schemes that allow one to compute the relevant exponents and prefactors of the dissipative terms in the equations of motion. We thereby propose a quantitatively accurate extension of the Hertzian model that encompasses dissipative effects via a discrete Laplacian of the velocities. Experiments and computations with steel, brass, and polytetrafluoroethylene reveal a common dissipation exponent with a material-dependent prefactor.

Short-pulse dynamics in strongly nonlinear dissipative granular chains

We study the energy decay properties of a pulse propagating in a strongly nonlinear granular chain with damping proportional to the relative velocity of the grains. We observe a wave disturbance that at low viscosities consists of two parts exhibiting two entirely different time scales of dissipation. One part is an attenuating solitary wave, dominated by discreteness and nonlinearity effects as in a dissipationless chain, and has the shorter lifetime. The other is a purely dissipative shocklike structure with a much longer lifetime and exists only in the presence of dissipation. The range of viscosities and initial configurations that lead to this complex wave disturbance are explored.

Energy dissipation and dispersion effects in granular media

The strong interactions between particles will make the energy within the granular materials propagate through the network of contacts and be partly dissipated. Establishing a model that can clearly classify the dissipation and dispersion effects is crucial for the understanding of the global behaviors in the granular materials. For particles with rate-independent material, the dissipation effects come from the local plastic deformation and can be constrained at the energy level by using energetic restitution coefficients. On the other hand, the dispersion effects should depend on the intrinsic nature of the interaction law between two particles. In terms of a bistiffness compliant contact model that obeys the energetical constraint defined by the energetic coefficients, our recent work related to the issue of multiple impacts indicates that the propagation of energy during collisions can be represented by a distributing law. In particular, this law shows that the dispersion effects are dominated by the relative contact stiffness and the relative potential energy stored at the contact points. In this paper, we will apply our theory to the investigation of the wave behavior in granular chain systems. The comparisons between our numerical results and the experimental ones by Falcon, [Eur. Phys. J. B 5, 111 (1998)] for a column of beads colliding against a wall show very good agreement and confirm some conclusions proposed by Falcon Other numerical results associated with the case of several particles impacting a chain, and the collisions between two so-called solitary waves in a Hertzian type chain are also presented.

Violation of the Porod law in a freely cooling granular gas in one dimension

We study a model of freely cooling inelastic granular gas in one dimension, with a restitution coefficient which approaches the elastic limit below a relative velocity scale delta. While at early times (t<<delta;{-1}) the gas behaves as a completely inelastic sticky gas conforming to predictions of earlier studies, at late times (t>>delta;{-1}) it exhibits a new fluctuation-dominated phase ordering state. We find distinct scaling behavior for the (i) density distribution function, (ii) occupied and empty gap distribution functions, (iii) the density structure function, and (iv) the velocity structure function, as compared to the completely inelastic sticky gas. The spatial structure functions (iii) and (iv) violate the Porod law. Within a mean-field approximation, the exponents describing the structure functions are related to those describing the spatial gap distribution functions.

Observation of two-wave structure in strongly nonlinear dissipative granular chains

In a strongly nonlinear viscous granular chain impacted by a single grain we observe a wave disturbance that consists of two parts exhibiting two time scales of dissipation. Above a critical viscosity there is no separation of the two pulses, and the dissipation and nonlinearity dominate the shocklike attenuating pulse.

Simulation study of granular compaction dynamics under vertical tapping

We study, by numerical simulation, the compaction dynamics of frictional hard disks in two dimensions, subjected to vertical shaking. Shaking is modeled by a series of vertical expansion of the disk packing, followed by dynamical recompression of the assembly under the action of gravity. The second phase of the shake cycle is based on an efficient event-driven molecular-dynamics algorithm. We analyze the compaction dynamics for various values of friction coefficient and coefficient of normal restitution. We find that the time evolution of the density is described by rho(t)=rho{infinity}-DeltarhoE{alpha}[-(ttau){alpha}], where E{alpha} denotes the Mittag-Leffler function of order 0<alpha<1 . The parameter tau is found to decay with tapping intensity Gamma according to a power law tau proportional, variantGamma{-gamma}, where parameter gamma is almost independent on the material properties of grains. Also, an expression for the grain mobility during the compaction process has been obtained. We characterize the local organization of disks in terms of contact connectivity and distribution of the Delaunay "free" volumes. Our analysis at microscopic scale provides evidence that compaction is mainly due to a decrease of the number of the largest pores. An interpretation of the memory effects observed for a discontinuous shift in tapping intensity Gamma is provided by the analysis of the time evolution of connectivity numbers and volume distribution of pores.

Flux through a hole from a shaken granular medium

We have measured the flux of grains from a hole in the bottom of a shaken container of grains. We find that the peak velocity of the vibration, v max, controls the flux, i.e., the flux is nearly independent of the frequency and acceleration amplitude for a given value of v max. The flux decreases with increasing peak velocity and then becomes almost constant for the largest values of v max. The data at low peak velocity can be quantitatively described by a simple model, but the crossover to nearly constant flux at larger peak velocity suggests a regime in which the granular density near the container bottom is independent of the energy input to the system.

Voronoi neighbor statistics of homogeneously sheared inelastic hard disks and hard spheres

In this work we extend our earlier Voronoi neighbor analysis [J. Chem. Phys. 123, 074502 (2005)] to homogeneously sheared inelastic hard-particle structures, the simplest model for rapid granular matter. The pair distribution function is partitioned into the nth neighbor coordination number (Cn), and the nth neighbor position distribution [fn(r)]. The distribution of the number of Voronoi faces (Pn) is also considered since C1 is its mean. We report the Cn and Pn for the homogeneously sheared inelastic hard-disk and hard-sphere structures. These statistics are sensitive to shear ordering transition, the nonequilibrium analogue of the freezing transition. In the near-elastic limit, the sheared fluid statistics approach that of the thermodynamic fluid. On shear ordering, due to the onset of order, the Cn for sheared structures approach that of the thermodynamic solid phase. The suppression of nucleation by the homogeneous shear is evident in these statistics. As inelasticity increases, the shear ordering packing fraction increases. Due to the topological instability of the isotropically perturbed face-centered cubic lattice, polyhedra with faces 12 to 18, with a mean at 14, coexist even in the regular close packed limit for the thermodynamic hard-sphere solid. In shear ordered inelastic hard-sphere structures there is a high incidence of 14-faceted polyhedra and a consequent depletion of polyhedra with faces 12, 13, 15-18, due to the formation of body-centered-tetragonal (bct) structures. These bct structures leave a body-centered-cubic-like signature in the Cn and Pn data. On shear ordering, close-packed layers slide past each other. However, with a velocity-dependent coefficient of restitution, at a critical shear rate these layers get disordered or amorphized. We find that the critical shear rate for amorphization is inversely proportional to the particle diameter, as compared to the inverse square scaling observed in dense colloidal suspensions.

Transport properties of dense dissipative hard-sphere fluids for arbitrary energy loss models

The revised Enskog approximation for a fluid of hard spheres which lose energy upon collision is discussed for the case that the energy is lost from the normal component of the velocity at collision but is otherwise arbitrary. Granular fluids with a velocity-dependent coefficient of restitution are an important special case covered by this model. A normal solution to the Enskog equation is developed using the Chapman-Enskog expansion. The lowest order solution describes the general homogeneous cooling state and a generating function formalism is introduced for the determination of the distribution function. The first order solution, evaluated in the lowest Sonine approximation, provides estimates for the transport coefficients for the Navier-Stokes hydrodynamic description. All calculations are performed in an arbitrary number of dimensions.