Shear instabilities in granular mixtures

Dynamical instabilities in fluid mechanics are responsible for a variety of important common phenomena, such as waves on the sea surface or Taylor vortices in Couette flow. In granular media dynamical instabilities have just begun to be discovered. Here we show by means of molecular dynamics simulation the existence of a new dynamical instability of a granular mixture under oscillating horizontal shear, which leads to the formation of a striped pattern where the components are segregated. We investigate the properties of such a Kelvin-Helmholtz-like instability and show how it is connected to pattern formation in granular flow and segregation.     

Lattice-BGK approach to simulating granular flows

Many continuum theories for granular flow produce an equation of motion for the fluctuating kinetic energy density (granular temperature) that accounts for the energy lost in inelastic collisions. Apart from the presence of an extra dissipative term, this equation is very similar in form to the usual temperature equation in hydrodynamics. It is shown how a lattice-kinetic model based on the Bhatnagar-Gross-Krook (BGK) equation that was previously derived for a miscible two-component fluid may be modified to model the continuum equations for granular flow. This is done by noting that the variable corresponding to the concentration of one species follows an equation that is essentially analogous to the granular temperature equation. A simulation of an unforced granular fluid using the modified model reproduces the phenomenon of clustering instability, namely the spontaneous agglomeration of particles into dense clusters, which occurs generically in all granular flows. The success of the continuum theory in capturing the gross features of this basic phenomenon is discussed. Some shear flow simulations are also presented.     

Destabilization of a Saffman-Taylor fingerlike pattern in a granular suspension

We study the Saffman-Taylor instability in a granular suspension formed by micrometric beads immersed in a viscous liquid. When using an effective viscosity for the flow of the suspension in the Hele-Shaw cell to define the control parameter of the system, the results for the finger width of stable fingers are found to be close to the classical results of Saffman-Taylor. One observes, however, an early destabilization of the fingers that can be attributed to the discrete nature of the individual grains. Classically, the threshold of destabilization is linked to the noise in the cell and is thus difficult to quantify. We show that the grains represent a "controlled noise" and produce an initial perturbation of the interface with an amplitude proportional to the grain size. The finite amplitude instability mechanism proposed by Bensimon et al. allows us to link this perturbation to the value of the threshold observed.     

Granular dynamics in compaction and stress relaxation

Elastic and dissipative properties of granular assemblies under uniaxial compression are studied both experimentally and by numerical simulations. Following a novel compaction procedure at varying oscillatory pressures, the stress response to a step strain reveals an exponential relaxation followed by a slow logarithmic decay. Simulations indicate that the latter arises from the coupling between damping and collective grain motion predominantly through sliding. We characterize an analogous "glass transition" for packed grains, below which the system shows aging in time-dependent sliding correlation functions.     

The influence of thes-f interaction on the conduction band behaviour in ferromagnetic semiconductors. A theory of the red shift effect

The mechanism of the thermomechanical instability of an elastically slowly strained solid is explored, by a normal mode analysis of temperature and strain fluctuations. The strain dependence of the Grüneisen parameter is found to play a major role, since it controls the rates of damping of the lowest wavevector acoustic excitations pumped by the imposed strain. For -titanium the instability is found to occur at a critical strain comparable with that of the elastic limit. The instability can be envisaged in terms of a fast growing deformation, presumably associated with an acoustic emission.Gruppo Nazionale di Struttura della Materia del Consiglio Nazionale delle Ricerche, UR7     

Compactivity and transmission of stress in granular materials

We outline a statistical-mechanical theory of granular materials. Stress propagation and force fluctuations in static granular media are still poorly understood. We develop the statistical-mechanical theory that delivers the fundamental equations of stress equilibrium. The formalism is based on the assumptions that grains are rigid, cohesionless, and that friction is perfect. Since grains are assumed perfectly rigid, no strain or displacement field can enter the equations for static equilibrium of the stress field. The complete system of equations for the stress tensor is derived from the equations of intergranular force and torque balance, given the geometric specification of the material. These new constitutive equations are indeed fundamental and are based on relations between various components of the stress tensor within the material, and depend on the topology of the granular packing. The problem of incorporating into the formalism the "no tensile forces" constraint is considered. The compactivity concept is reviewed. We discuss the relation between the concept of compactivity and the problem of stress transmission. (c) 1999 American Institute of Physics.     

Nonlinear dynamics of density waves in granular flows through narrow vertical channels

We study granular flows through narrow channels driven by gravity in the framework of the kinetic theory for dissipative dense gases. We derive equations of motion for quasi-one-dimensional systems. In a certain range of flow density, the steady homogeneous regime is found to be unstable against the formation of density waves. We show moreover that near the onset of the instability, the governing equation for the flow density is a mixture of the Korteweg-de-Vries equation, which leads to soliton, and the Bürger equation which exhibits spatio-temporal chaos. The competition between chaos and solitons may lead either to regular spatially ordered density waves or to chaotic dynamics. We argue that these two types of dynamics can be encountered experimentally according to the channel width and the dissipative properties of the granular media. Received: 11 March 1998 / Revised and Accepted: 3 July 1998     

Jamming and fluctuations in granular drag

We investigate the dynamic evolution of jamming in granular media through fluctuations in the granular drag force. The successive collapse and formation of jammed states give a stick-slip nature to the fluctuations which is independent of the contact surface between the grains and the dragged object, thus implying that the stress-induced collapse is nucleated in the bulk of the granular sample. We also find that while the fluctuations are periodic at small depths, they become "stepped" at large depths, a transition which we interpret as a consequence of the long-range nature of the force chains.     

Convection, heaping, and cracking in vertically vibrated granular slurries

We address the situation in which vibration is applied to thin layers of granular material with water filling the pore space, but with air above. Beyond a critical drive parameter we observe heap-shaped convecting domains of nontrivial topology, which exhibit cracks as the forcing amplitude is further increased. We summarize these results in a phase diagram, understand the onset of the convecting heaps as a manifestation of the Rayleigh-Taylor instability for fluids, and measure the response of isolated convecting structures.     

Washboard road: the dynamics of granular ripples formed by rolling wheels

We report laboratory experiments on rippled granular surfaces formed under rolling wheels. Ripples appear above a critical speed and drift slowly in the driving direction. Ripples coarsen as they saturate and exhibit ripple creation and destruction events. All of these effects are captured qualitatively by 2D soft-particle simulations in which a disk rolls over smaller disks in a periodic box. The simulations show that compaction and segregation are inessential to the ripple phenomenon. We describe a simplified scaling model which gives some insight into the mechanism of the instability.     

Active thermal convection in vibrofluidized granular systems

We present a granular-hydrodynamic model that captures the essence of convection in a fully vibrofluidized granular system. The steady temperature distribution is solved analytically. Numerical simulation shows that the convection always develops through a supercritical bifurcation, with its energy about of the random (heat) one. A comparison calculation is performed for a normal fluid. The convection roll, or an active roll as we call it, has an angular velocity gradient from its interior to exterior. We conclude that active rolls are universal.Received: 25 March 2004, Published online: 9 September 2004PACS: 45.70.Mg Granular flow: mixing, segregation and stratification - 47.20.Bp Buoyancy-driven instability - 47.27.Te Convection and heat transfer     

Inhomogeneous Relaxation in Vibrated Granular Media: Consolidation Waves

We investigate the spatial dependence of the density of vibrated granular beds, using simulations based on a hybrid Monte Carlo algorithm. We find that the initial consolidation is typically inhomogeneous, both in the presence of a constant shaking intensity and when the granular bed is subjected to "annealed shaking". We also present a theoretical model which explains such inhomogeneous relaxation in terms of a "consolidation wave", in good qualitative agreement with our simulations. Our results are also in qualitative agreement with recent experiments.     

Convection in horizontally shaken granular material

In horizontally shaken granular material different types of pattern formation have been reported. We want to deal with the convection instability which has been observed in experiments and which recently has been investigated numerically. Using two dimensional molecular dynamics we show that the convection pattern depends crucially on the inelastic properties of the material. The concept of restitution coefficient provides arguments for the change of the behaviour with varying inelasticity. Received 3 March 1999     

Force fluctuations in a vertically pushed granular column

We present series of experiments on the resistance force encountered by a bottom piston pushing a vertical granular column confined in a two-dimensional cell. We show that, due to the presence of friction at the boundaries and between the grains, the signal shows many complex features. At slow driving velocities, we observe a transition to a stick-slip dynamic instability. Depending on the granular material used, the elementary stick-slip events may either be well characterized or largely distributed. We present a statistical study on the waiting time between events and the distribution of energy release as a function of the spring stiffness and the driving velocity. Received 5 August 1998 and Received in final form 22 October 1998     

Dense granular flow around a penetrating object: experiment and hydrodynamic model

We present in this Letter experimental results on the bidimensional flow field around a cylinder penetrating into dense granular matter, together with drag force measurements. A hydrodynamic model based on extended kinetic theory for dense granular flow reproduces well the flow localization close to the cylinder and the corresponding scalings of the drag force, which is found to not depend on velocity, but linearly on the pressure and on the cylinder diameter and weakly on the grain size. Such a regime is found to be valid at a low enough "granular" Reynolds number.     

A model for ripple instabilities in granular media

We extend the model of surface granular flow proposed in [#!bcre!#] to account for the effect of an external `wind', which acts as to dislodge particles from the static bed, such that a stationary state of flowing grains is reached. We discuss in detail how this mechanism can be described in a phenomenological way, and show that a flat bed is linearly unstable against ripple formation in a certain region of parameter space. We focus in particular on the (realistic) case where the migration velocity of the instability is much smaller than the grains' velocity. In this limit, the full dispersion relation can be established. We relate the critical wave vector to the mean hopping length and to the ratio of the flight time to the `stick' time. We provide an intuitive interpretation of the instability. Received: 30 January 1998 / Revised: 12 May 1998 / Accepted: 8 June 1998     

Protocol dependence of mechanical properties in granular systems

We study the protocol dependence of the mechanical properties of granular media by means of computer simulations. We control a protocol of realizing disk packings in a systematic manner. In 2D, by keeping material properties of the constituents identical, we carry out compaction with various strain rates. The disk packings exhibit the strain rate dependence of the critical packing fraction above which the pressure becomes non-zero. The observed behavior contrasts with the well-studied jamming transitions for frictionless disk packings. We also observe that the elastic moduli of the disk packings depend on the strain rate logarithmically. Our results suggest that there exists a time-dependent state variable to describe macroscopic material properties of disk packings, which depend on its protocol.     

Dilatancy and friction in sheared granular media

We introduce a simple model to describe the frictional properties of granular media under shear. We model the friction force in terms of the horizontal velocity and the vertical position z of the slider, interpreting z as a constitutive variable characterizing the contact. Dilatancy is shown to play an essential role in the dynamics, inducing a stick-slip instability at low velocity. We compute the phase diagram, analyze numerically the model for a wide range of parameters and compare our results with experiments on dry and wet granular media, obtaining a good agreement. In particular, we reproduce the hysteretic velocity dependence of the frictional force. Received 16 November 1999     

Granular size segregation in underwater sand ripples

We report an experimental study of a binary sand bed under an oscillating water flow. The formation and evolution of ripples is observed. The appearance of a granular segregation is shown to strongly depend on the sand bed preparation. The initial wavelength of the mixture is measured. In the final steady state, a segregation in volume is observed instead of a segregation at the surface as reported before. The correlation between this phenomenon and the fluid flow is emphasised. Finally, different exotic patterns and their geophysical implications are presented.Received: 25 July 2003, Published online: 25 March 2004PACS: 45.70.Qj Pattern formation - 45.70.Mg Granular flow: mixing, segregation and stratification - 47.20.Ma Interfacial instability     

Towards a Landau–Ginzburg-Type Theory for Granular Fluids

In this paper we show how, under certain restrictions, the hydrodynamic equations for the freely evolving granular fluid fit within the framework of the time dependent Landau–Ginzburg (LG) models for critical and unstable fluids. The granular fluid, which is usually modeled as a fluid of inelastic hard spheres (IHS), exhibits two instabilities: the spontaneous formation of vortices and of high density clusters. We suppress the clustering instability by imposing constraints on the system sizes, in order to illustrate how LG-equations can be derived for the order parameter, being the rate of deformation or shear rate tensor, which controls the formation of vortex patterns. From the shape of the energy functional we obtain the stationary patterns in the flow field. Quantitative predictions of this theory for the stationary states agree well with molecular dynamics simulations of a fluid of inelastic hard disks.