Sedimentation instabilities: impact of the fluid compressibility and viscosity

Drainage explains soil liquefaction beyond the earthquake near-field

Earthquake-induced soil-liquefaction is a devastating phenomenon associated with loss of soil rigidity due to seismic shaking, resulting in catastrophic liquid-like soil deformation. Traditionally, liquefaction is viewed as an effectively undrained process. However, since undrained liquefaction only initiates under high energy density, most earthquake liquefaction events remain unexplained, since they initiate far from the earthquake epicenter, under low energy density. Here we show that liquefaction can occur under drained conditions at remarkably low seismic-energy density, offering a general explanation for earthquake far-field liquefaction. Drained conditions promote interstitial fluid flow across the soil during earthquakes, leading to excess pore pressure gradients and loss of soil strength. Drained liquefaction is triggered rapidly and controlled by a propagating compaction front, whose velocity depends on the seismic-energy injection rate. Our findings highlight the importance of considering soil liquefaction under a spectrum of drainage conditions, with critical implications for liquefaction potential assessments and hazards.

Sink versus tilt penetration into shaken dry granular matter: The role of the foundation

We study the behavior of cylindrical objects as they sink into a dry granular bed fluidized due to lateral oscillations. Somewhat unexpectedly, we have found that, within a large range of lateral shaking powers, cylinders with flat bottoms sink vertically, while those with a "foundation" consisting of a shallow ring attached to their bottom, tilt besides sinking. The latter scenario seems to dominate independently from the nature of the foundation when strong enough lateral vibrations are applied. We are able to explain the observed behavior by quasi-2D numerical simulations, which also demonstrate the influence of the intruder's aspect ratio. The vertical sink dynamics is explained with the help of a Newtonian equation of motion for the intruder. Our findings may shed light on the behavior of buildings and other manmade structures during earthquakes.

Sinking during earthquakes: Critical acceleration criteria control drained soil liquefaction

This article focuses on liquefaction of saturated granular soils, triggered by earthquakes. Liquefaction is defined here as the transition from a rigid state, in which the granular soil layer supports structures placed on its surface, to a fluidlike state, in which structures placed initially on the surface sink to their isostatic depth within the granular layer. We suggest a simple theoretical model for soil liquefaction and show that buoyancy caused by the presence of water inside a granular medium has a dramatic influence on the stability of an intruder resting at the surface of the medium. We confirm this hypothesis by comparison with laboratory experiments and discrete-element numerical simulations. The external excitation representing ground motion during earthquakes is simulated via horizontal sinusoidal oscillations of controlled frequency and amplitude. In the experiments, we use particles only slightly denser than water, which as predicted theoretically increases the effect of liquefaction and allows clear depth-of-sinking measurements. In the simulations, a micromechanical model simulates grains using molecular dynamics with friction between neighbors. The effect of the fluid is captured by taking into account buoyancy effects on the grains when they are immersed. We show that the motion of an intruder inside a granular medium is mainly dependent on the peak acceleration of the ground motion and establish a phase diagram for the conditions under which liquefaction happens, depending on the soil bulk density, friction properties, presence of water, and peak acceleration of the imposed large-scale soil vibrations. We establish that in liquefaction conditions, most cases relax toward an equilibrium position following an exponential in time. We also show that the equilibrium position itself, for most liquefaction regimes, corresponds to the isostatic equilibrium of the intruder inside a medium of effective density. The characteristic time to relaxation is shown to be essentially a function of the peak ground velocity.

Pneumatic fractures in confined granular media

We perform experiments where air is injected at a constant overpressure P_{in}, ranging from 5 to 250 kPa, into a dry granular medium confined within a horizontal linear Hele-Shaw cell. The setup allows us to explore compacted configurations by preventing decompaction at the outer boundary, i.e., the cell outlet has a semipermeable filter such that beads are stopped while air can pass. We study the emerging patterns and dynamic growth of channels in the granular media due to fluid flow, by analyzing images captured with a high speed camera (1000 images/s). We identify four qualitatively different flow regimes, depending on the imposed overpressure, ranging from no channel formation for P_{in} below 10 kPa, to large thick channels formed by erosion and fingers merging for high P_{in} around 200 kPa. The flow regimes where channels form are characterized by typical finger thickness, final depth into the medium, and growth dynamics. The shape of the finger tips during growth is studied by looking at the finger width w as function of distance d from the tip. The tip profile is found to follow w(d)∝d^{β}, where β=0.68 is a typical value for all experiments, also over time. This indicates a singularity in the curvature d^{2}d/dw^{2}∼κ∼d^{1-2β}, but not of the slope dw/dd∼d^{β-1}, i.e., more rounded tips rather than pointy cusps, as they would be for the case β>1. For increasing P_{in}, the channels generally grow faster and deeper into the medium. We show that the channel length along the flow direction has a linear growth with time initially, followed by a power-law decay of growth velocity with time as the channel approaches its final length. A closer look reveals that the initial growth velocity v_{0} is found to scale with injection pressure as v_{0}∝P_{in}^{3/2}, while at a critical time t_{c} there is a cross-over to the behavior v(t)∝t^{-α}, where α is close to 2.5 for all experiments. Finally, we explore the fractal dimension of the fully developed patterns. For example, for patterns resulting from intermediate P_{in} around 100-150 kPa, we find that the box-counting dimensions lie within the range D_{B}∈[1.53,1.62], similar to viscous fingering fractals in porous media.

Note: "Lock-in accelerometry" to follow sink dynamics in shaken granular matter

Understanding the penetration dynamics of intruders in granular beds is relevant not only for fundamental physics, but also for geophysical processes and construction on sediments or granular soils in areas potentially affected by earthquakes. While the penetration of intruders in two dimensional (2D) laboratory granular beds can be followed using video recording, this is useless in three dimensional (3D) beds of non-transparent materials such as common sand. Here, we propose a method to quantify the sink dynamics of an intruder into laterally shaken granular beds based on the temporal correlations between the signals from a reference accelerometer fixed to the shaken granular bed, and a probe accelerometer deployed inside the intruder. Due to its analogy with the working principle of a lock-in amplifier, we call this technique lock-in accelerometry.

Dynamic aerofracture of dense granular packings

A transition in hydraulically induced granular displacement patterns is studied by means of discrete numerical molecular dynamics simulations. During this transition the patterns change from fractures and fingers to finely dispersed bubbles. The dynamics of the displacement patterns are studied in a rectangular Hele-Shaw cell filled with a dense but permeable two-dimensional granular layer. At one side of the cell the pressure of the compressible interstitial gas is increased. At the opposite side from the inlet of the cell a semipermeable boundary is located. This boundary is only permeable towards the gas phase while preventing grains from leaving the cell. The imposed pressure gradient compacts the grains. In the process we can identify and describe a mechanism that controls the transition of the emerging displacement patterns from fractures and fingers to finely dispersed bubbles as a function of the interstitial gas's properties and the characteristics of the granular phase.