Model for the propagation of sound in granular materials

Softening of stressed granular packings with resonant sound waves

We perform numerical simulations of a two-dimensional bidisperse granular packing subjected to both a static confining pressure and a sinusoidal dynamic forcing applied by a wall on one edge of the packing. We measure the response experienced by a wall on the opposite edge of the packing and obtain the resonant frequency of the packing as the static or dynamic pressures are varied. Under increasing static pressure, the resonant frequency increases, indicating a velocity increase of elastic waves propagating through the packing. In contrast, when the dynamic amplitude is increased for fixed static pressure, the resonant frequency decreases, indicating a decrease in the wave velocity. This occurs both for compressional and for shear dynamic forcing and is in agreement with experimental results. We find that the average contact number Zc at the resonant frequency decreases with increasing dynamic amplitude, indicating that the elastic softening of the packing is associated with a reduced number of grain-grain contacts through which the elastic waves can travel. We image the excitations created in the packing and show that there are localized disturbances or soft spots that become more prevalent with increasing dynamic amplitude. Our results are in agreement with experiments on glass bead packings and earth materials such as sandstone and granite and may be relevant to the decrease in elastic wave velocities that has been observed to occur near fault zones after strong earthquakes, in surficial sediments during strong ground motion, and in structures during earthquake excitation.

Elastic wave propagation in confined granular systems

We present numerical simulations of acoustic wave propagation in confined granular systems consisting of particles interacting with the three-dimensional Hertz-Mindlin force law. The response to a short mechanical excitation on one side of the system is found to be a propagating coherent wave front followed by random oscillations made of multiply scattered waves. We find that the coherent wave front is insensitive to details of the packing: force chains do not play an important role in determining this wave front. The coherent wave propagates linearly in time, and its amplitude and width depend as a power law on distance, while its velocity is roughly compatible with the predictions of macroscopic elasticity. As there is at present no theory for the broadening and decay of the coherent wave, we numerically and analytically study pulse propagation in a one-dimensional chain of identical elastic balls. The results for the broadening and decay exponents of this system differ significantly from those of the random packings. In all our simulations, the speed of the coherent wave front scales with pressure as p1/6; we compare this result with experimental data on various granular systems where deviations from the p1/6 behavior are seen. We briefly discuss the eigenmodes of the system and effects of damping are investigated as well.

Characterization of soliton damping in the granular chain under gravity

A soliton created in the horizontal granular chain damps due to gravity in the vertical chain. We show that there are two types of propagating modes, quasisolitary and oscillatory, in the vertical chain, depending on the strength of impulse. We find that the type of damping is a power law in depth or time. We also find that the absolute value of the exponent of the power law decreases as the strength of the initial impulse increases in the quasisolitary regime. In the oscillatory regime, however, in which the initial impulse is weak, the power-law exponent is independent of the strength of the initial impulse. We show that the power-law damping is caused by the gravitation which results in the change of the force constant at each contact.