Sand stirred by chaotic advection

We study the spatial structure of a granular material, N particles subject to inelastic mutual collisions, when it is stirred by a bidimensional smooth chaotic flow. A simple dynamical model is introduced where four different time scales are explicitly considered: i) the Stokes time, accounting for the inertia of the particles, ii) the mean collision time among the grains, iii) the typical time scale of the flow, and iv) the inverse of the Lyapunov exponent of the chaotic flow, which gives a typical time for the separation of two initially close parcels of fluid. Depending on the relative values of these different times a complex scenario appears for the long-time steady spatial distribution of particles, where clusters of particles may or not appear.Comment: 4 pages, 3 figure

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 to 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.Comment: 11 pages, 14 figure

Effect of boundaries on the force distributions in granular media

The effect of boundaries on the force distributions in granular media is illustrated by simulations of 2D packings of frictionless, Hertzian spheres. To elucidate discrepancies between experimental observations and theoretical predictions, we distinguish between the weight distribution {\cal P} (w) measured in experiments and analyzed in the q-model, and the distribution of interparticle forces P(f). The latter one is robust, while {\cal P}(w) can be obtained once the local packing geometry and P(f) are known. By manipulating the (boundary) geometry, we show that {\cal P}(w) can be varied drastically.Comment: 4 pages, 4 figure

Granular discharge and clogging for tilted hoppers

We measure the flux of spherical glass beads through a hole as a systematic function of both tilt angle and hole diameter, for two different size beads. The discharge increases with hole diameter in accord with the Beverloo relation for both horizontal and vertical holes, but in the latter case with a larger small-hole cutoff. For large holes the flux decreases linearly in cosine of the tilt angle, vanishing smoothly somewhat below the angle of repose. For small holes it vanishes abruptly at a smaller angle. The conditions for zero flux are discussed in the context of a {\it clogging phase diagram} of flow state vs tilt angle and ratio of hole to grain size

Hydrodynamic theory for granular gases

A granular gas subjected to a permanent injection of energy is described by means of hydrodynamic equations derived from a moment expansion method. The method uses as reference function not a Maxwellian distribution $f_{\sf M}$ but a distribution $f_0 = \Phi f_{\sf M}$, such that $\Phi$ adds a fourth cumulant $\kappa$ to the velocity distribution. The formalism is applied to a stationary conductive case showing that the theory fits extraordinarily well the results coming from our molecular dynamic simulations once we determine $\kappa$ as a function of the inelasticity of the particle-particle collisions. The shape of $\kappa$ is independent of the size $N$ of the system.Comment: 10 pages, 9 figures, more about our research in http://www.cec.uchile.cl/cinetica

Generation of Porous Particle Structures using the Void Expansion Method

The newly developed "void expansion method" allows for an efficient generation of porous packings of spherical particles over a wide range of volume fractions using the discrete element method. Particles are randomly placed under addition of much smaller "void-particles". Then, the void-particle radius is increased repeatedly, thereby rearranging the structural particles until formation of a dense particle packing. The structural particles' mean coordination number was used to characterize the evolving microstructures. At some void radius, a transition from an initially low to a higher mean coordination number is found, which was used to characterize the influence of the various simulation parameters. For structural and void-particle stiffnesses of the same order of magnitude, the transition is found at constant total volume fraction slightly below the random close packing limit. For decreasing void-particle stiffness the transition is shifted towards a smaller void-particle radius and becomes smoother.Comment: 9 pages, 8 figure

Phenomenological glass model for vibratory granular compaction

A model for weakly excited granular media is derived by combining the free volume argument of Nowak et al. [Phys. Rev. E 57, 1971 (1998)] and the phenomenological model for supercooled liquids of Adam and Gibbs [J. Chem. Phys. 43, 139 (1965)]. This is made possible by relating the granular excitation parameter \Gamma, defined as the peak acceleration of the driving pulse scaled by gravity, to a temperature-like parameter \eta(\Gamma). The resulting master equation is formally identical to that of Bouchaud's trap model for glasses [J. Phys. I 2, 1705 (1992)]. Analytic and simulation results are shown to compare favourably with a range of known experimental behaviour. This includes the logarithmic densification and power spectrum of fluctuations under constant \eta, the annealing curve when \eta is varied cyclically in time, and memory effects observed for a discontinuous shift in \eta. Finally, we discuss the physical interpretation of the model parameters and suggest further experiments for this class of systems.Comment: 2 references added; some figure labels tweaked. To appear in PR

On the validity of the Boltzmann equation to describe low density granular systems

The departure of a granular gas in the instable region of parameters from the initial homogeneous cooling state is studied. Results from Molecular Dynamics and from Direct Monte Carlo simulation of the Boltzmann equation are compared. It is shown that the Boltzmann equation accurately predicts the low density limit of the system. The relevant role played by the parallelization of the velocities as time proceeds and the dependence of this effect on the density is analyzed in detail

A nonlinear hydrodynamical approach to granular materials

We propose a nonlinear hydrodynamical model of granular materials. We show how this model describes the formation of a sand pile from a homogeneous distribution of material under gravity, and then discuss a simulation of a rotating sandpile which shows, in qualitative agreement with experiment, a static and dynamic angle of repose.Comment: 17 pages, 14 figures, RevTeX4; minor changes to wording and some additional discussion. Accepted by Phys. Rev.

Confined granular packings: structure, stress, and forces

The structure and stresses of static granular packs in cylindrical containers are studied using large-scale discrete element molecular dynamics simulations in three dimensions. We generate packings by both pouring and sedimentation and examine how the final state depends on the method of construction. The vertical stress becomes depth-independent for deep piles and we compare these stress depth-profiles to the classical Janssen theory. The majority of the tangential forces for particle-wall contacts are found to be close to the Coulomb failure criterion, in agreement with the theory of Janssen, while particle-particle contacts in the bulk are far from the Coulomb criterion. In addition, we show that a linear hydrostatic-like region at the top of the packings unexplained by the Janssen theory arises because most of the particle-wall tangential forces in this region are far from the Coulomb yield criterion. The distributions of particle-particle and particle-wall contact forces $P(f)$ exhibit exponential-like decay at large forces in agreement with previous studies.Comment: 11 pages, 11 figures, submitted to PRE (v2) added new references, fixed typo