Nonlinear waves in dry and wet Hertzian granular chains

A one-dimensional dry granular medium, a chain of beads which interact via the nonlinear Hertz potential, exhibits strongly nonlinear behaviors. When such an alignment further contains some fluid in the interstices between grains, it may exhibit new interesting features. We report some recent experiments, analysis and numerical simulations concerning nonlinear wave propagation in dry and wet chains of spheres. We consider first a monodisperse chain as a reference case. We then analyze how the pulse characteristics are modified in the presence of an interstitial viscous fluid. The fluid not only induces dissipation but also strongly affect the intergrain stiffness: in a wet chain, wave speed is enhanced and pulses are shorter. Simple experiments performed with a single sphere colliding a wall covered by a thin film of fluid confirm these observations. We demonstrate that even a very small amount of fluid can overcome the Hertzian potential and is responsible for a large increase of contact stiffness. Possible mechanisms for wet contact hardening are related to large fluid shear rate during fast elastohydrodynamic collision between grains.     

A reflectivity equation for acoustic media with absorption

The wave equation which governs the propagation of dilatational stress (pressure) in viscoacoustic media can be reduced to the lossless acoustic wave equation by the introduction of a complex traveltime coordinate. An interesting property of this equation is that it contains a term which is proportional to the reflectivity function in its single scattering approximation.     

Jamming, yielding, and rheology of weakly vibrated granular media

We establish that the rheological curve of dry granular media is nonmonotonic, both in the presence and absence of external mechanical agitations. In the presence of weak vibrations, the nonmonotonic flow curves govern a hysteretic transition between slow but steady and fast, inertial flows. In the absence of vibrations, the nonmonotonic flow curve governs the yielding behavior of granular media. Finally, we show that nonmonotonic flow curves can be seen in at least two different flow geometries and for several granular materials.     

A new empirical model for the acoustic properties of loose granular media

This paper examines physical parameters of loose granular mixes and their empirical relations to the acoustic performance of these mixes. In this work a new classification of granular media has been proposed which is related to the characteristic particle dimension and the specific density of the grain base. It has been shown that this classification is a useful characteristic for rapid evaluation of the acoustic performance of loose granular mixes. The characteristic impedance and propagation constant have been measured for a representative selection of grain mixes and used to develop a new empirical model. This model relates the above acoustic characteristics to the characteristic particle dimension, porosity, tortuosity and specific density of the grain base, which are routinely measurable parameters. A very good agreement with the experimental data is illustrated in the frequency range of 250-4000 Hz for materials with the grain base of 0.4-3.5 mm and specific densities between 200 and 1200 kg/m³. Unlike many theoretical models for the prediction of the acoustic properties of porous media, the proposed expressions do not involve any special functions of complex argument, empirical shape factors or sophisticated characteristics of porous structure. These are practical enough to be of interest to acoustic and noise control engineers and material manufacturers.     

Segregation transitions in wet granular matter

We report the effect of interstitial fluid on the extent of segregation by imaging the pile that results after bidisperse color-coded particles are poured into a silo. Segregation is sharply reduced and preferential clumping of small particles is observed when a small volume fraction of fluid V(f) is added. We find that viscous forces in addition to capillary forces have an important effect on the extent of segregation s and the angle of repose straight theta. We show that the sharp initial change and the subsequent saturation in s and straight theta occurs over similar V(f). We also find that a transition back to segregation can occur when the particles are completely immersed in a fluid at low viscosities.     

Avalanche dynamics in wet granular materials

We have studied the dynamics of avalanching wet granular media in a rotating drum apparatus. Quantitative measurements of the flow velocity and the granular flux during avalanches allow us to characterize novel avalanche types unique to wet media. We also explore the details of viscoplastic flow (observed at the highest liquid contents) in which there are lasting contacts during flow, leading to coherence across the entire sample. This coherence leads to a velocity-independent flow depth at high rotation rates and novel robust pattern formation in the granular surface.     

Fluctuations in granular media

Dense slowly evolving or static granular materials exhibit strong force fluctuations even though the spatial disorder of the grains is relatively weak. Typically, forces are carried preferentially along a network of "force chains." These consist of linearly aligned grains with larger-than-average force. A growing body of work has explored the nature of these fluctuations. We first briefly review recent work concerning stress fluctuations. We then focus on a series of experiments in both two- and three-dimension [(2D) and (3D)] to characterize force fluctuations in slowly sheared systems. Both sets of experiments show strong temporal fluctuations in the local stress/force; the length scales of these fluctuations extend up to 10(2) grains. In 2D, we use photoelastic disks that permit visualization of the internal force structure. From this we can make comparisons to recent models and calculations that predict the distributions of forces. Typically, these models indicate that the distributions should fall off exponentially at large force. We find in the experiments that the force distributions change systematically as we change the mean packing fraction, gamma. For gamma's typical of dense packings of nondeformable grains, we see distributions that are consistent with an exponential decrease at large forces. For both lower and higher gamma, the observed force distributions appear to differ from this prediction, with a more Gaussian distribution at larger gamma and perhaps a power law at lower gamma. For high gamma, the distributions differ from this prediction because the grains begin to deform, allowing more grains to carry the applied force, and causing the distributions to have a local maximum at nonzero force. It is less clear why the distributions differ from the models at lower gamma. An exploration in gamma has led to the discovery of an interesting continuous or "critical" transition (the strengthening/softening transition) in which the mean stress is the order parameter, and the mean packing fraction, gamma, must be adjusted to a value gamma(c) to reach the "critical point." We also follow the motion of individual disks and obtain detailed statistical information on the kinematics, including velocities and particle rotations or spin. Distributions for the azimuthal velocity, V(theta), and spin, S, of the particles are nearly rate invariant, which is consistent with conventional wisdom. Near gamma(c), the grain motion becomes intermittent causing the mean velocity of grains to slow down. Also, the length of stress chains grows as gamma-->gamma(c). The 3D experiments show statistical rate invariance for the stress in the sense that when the power spectra and spectral frequencies of the stress time series are appropriately scaled by the shear rate, Omega, all spectra collapse onto a single curve for given particle and sample sizes. The frequency dependence of the spectra can be characterized by two different power laws, P proportional, variant omega(-alpha), in the high and low frequency regimes: alpha approximately 2 at high omega; alpha<2 at low omega. The force distributions computed from the 3D stress time series are at least qualitatively consistent with exponential fall-off at large stresses. (c) 1999 American Institute of Physics.     

Mixing and condensation in a wet granular medium

We have studied the effect of small amounts of added liquid on the dynamic behavior of a granular system consisting of a mixture of glass beads of two different sizes. Segregation of the large beads to the top of the sample is found to depend in a nontrivial way on the liquid content. A transition to viscoplastic behavior occurs at a critical liquid content, which depends upon the bead size. We show that this transition can be interpreted as a condensation due to the hysteretic liquid bridge forces connecting the beads, and we provide the corresponding phase diagram.     

Sound wave propagation in weakly polydisperse granular materials

Dynamic simulations of wave propagation are performed in dense granular media with a narrow polydisperse size-distribution and a linear contact-force law. A small perturbation is created on one side of a static packing and its propagation, for both P- and S-waves, is examined. A size variation comparable to the typical contact deformation already changes sound propagation considerably. The transmission spectrum becomes discontinuous, i.e., a lower frequency band is transmitted well, while higher frequencies are not, possibly due to attenuation and scattering. This behaviour is qualitatively reproduced for (i) Hertz non-linear contacts, for (ii) frictional contacts, (iii) for a range of smaller amplitudes, or (iv) for larger systems. This proves that the observed wave propagation and dispersion behaviour is intrinsic and not just an artifact of (i) a linear model, (ii) a frictionless packing, (iii) a large amplitude non-linear wave, or (iv) a finite size effect.     

Computation of acoustic absorption in media composed of packed microtubes exhibiting surface irregularity

A multi-scale homogenization technique and a finite element-based solution procedure are employed to compute acoustic absorption in smooth and rough packed microtubes. The absorption considered arises from thermo-viscous interactions between the fluid media and the microtube walls. The homogenization technique requires geometric periodicity, which for smooth tubes is invoked using the periodicity of the finite element mesh; for rough microtubes, the periodicity invoked is that associated with the roughness. Analysis of the packed configurations, for the specific microtube radii considered, demonstrates that surface roughness does not appreciably increase the overall absorption, but instead shifts the peaks and values of the absorption curve. Additionally, the effect of the fluid media temperature on acoustic absorption is also explored. The results of the investigation are used to make conclusions about tailored design of acoustically absorbing microtube-based materials.     

Conditions for forming weakly diverging acoustic bundles in ocean waveguides

Necessary and sufficient conditions for the dependence of the cycle length of the Brilloin waves or rays of geometrical acoustics on the ray parameter that is inversely proportional to their phase velocity are formulated. The formulated conditions determine the existence of weakly diverging acoustic bundles in vertically stratified oceanic waveguides, even with a point sound source.     

Strain localisation in granular media

This paper discusses strain localisation in granular media by presenting experimental, full-field analysis of mechanical tests on sand, both at a continuum level, as well as at the grain scale. At the continuum level, the development of structures of localised strain can be studied. Even at this scale, the characteristic size of the phenomena observed is in the order of a few grains. In the second part of this paper, therefore, the development of shear bands within specimen of different sands is studied at the level of the individual grains, measuring grains kinematics with x-ray tomography. The link between grain angularity and grain rotation within shear bands is shown, allowing a grain-scale explanation of the difference in macroscopic residual stresses for materials with different grain shapes. Finally, rarely described precursors of localisation, emerging well before the stress peak are observed and commented.     

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     

Structure and cluster formation in granular media

The two most important phenomena at the basis of granular media are excluded volume and dissipation. The former is captured by the hard sphere model and is responsible for, e.g., crystallization, the latter leads to interesting structures like clusters in non-equilibrium dynamical, freely cooling states. The freely cooling system is examined concerning the energy decay and the cluster evolution in time. Corrections for crystallization and multi-particle contacts are provided, which become more and more important with increasing density.     

Unclustering transition in freely cooling wet granular matter

The free cooling of one-dimensional wet granular matter is presented in the framework of the minimal capillary model. We demonstrate analytically and by extensive simulations that above a critical density, the clustering of wet granular matter is not monotonic in time, but undergoes a sharp unclustering transition. This precipitation of granular droplets out of the gas takes place when the granular temperature comes close to the energy scale set by the capillary interaction.     

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     

Structures and non-equilibrium dynamics in granular media

In granular media, dissipation leads to interesting phenomena like cluster formation in non-equilibrium dynamical states. As an example, the freely cooling system is examined concerning the energy decay and the cluster evolution with time. Furthermore, the probability distribution of the collision frequency is discussed. Uncorrelated events lead to a Poisson distribution for the collision frequencies in the homogeneous system, whereas cooperative phenomena can be related to a power-law decay of the collision probability per unit time. To cite this article: S. Luding, C. R. Physique 3 (2002) 153–161.     

A k-space Green's function solution for acoustic initial value problems in homogeneous media with power law absorption

An efficient Green's function solution for acoustic initial value problems in homogeneous media with power law absorption is derived. The solution is based on the homogeneous wave equation for lossless media with two additional terms. These terms are dependent on the fractional Laplacian and separately account for power law absorption and dispersion. Given initial conditions for the pressure and its temporal derivative, the solution allows the pressure field for any time t>0 to be calculated in a single step using the Fourier transform and an exact k-space time propagator. For regularly spaced Cartesian grids, the former can be computed efficiently using the fast Fourier transform. Because no time stepping is required, the solution facilitates the efficient computation of the pressure field in one, two, or three dimensions without stability constraints. Several computational aspects of the solution are discussed, including the effect of using a truncated Fourier series to represent discrete initial conditions, the use of smoothing, and the properties of the encapsulated absorption and dispersion.