Triggering granular avalanches with ultrasound

Metastable state preceding shear zone instability: Implications for earthquake-accelerated landslides and dynamic triggering

Understanding the dynamic response of granular shear zones under cyclic loading is fundamental to elucidating the mechanisms triggering earthquake-induced landslides, with implications for broader fields such as seismology and granular physics. Existing prediction methods struggle to accurately predict many experimental and in situ landslide observations due to inadequate consideration of the underlying physical mechanisms. The mechanisms that influence landslide dynamic triggering, a transition from static (or extremely slow creeping) to rapid runout, remain elusive. Herein, we focus on the inherent physics of granular shear zones under dynamic loading using ring shear experiments. Except for coseismic slip caused by the dynamic load, varying magnitudes of postseismic creep with increasing cycles of dynamic loading are observed, highlighting the effects of coseismic weakening (shear zone fatigue) and subsequent postseismic healing. A metastable state, characterized by a significant increase in postseismic creep, typically precedes shear zone instability. The metastable state may arise as weakened shear resistance approaches the applied shear stress, demonstrating a phase transition from a solid-like state to a fluid state (plastic granular flow). The metastable state may potentially indicate the shear zone's stress state and serve as a precursor to impending instability. Furthermore, the proposed mechanisms offer a compelling explanation for the widespread postseismic landslide movement following earthquakes. Incorporating these mechanisms into the Newmark method has the potential to improve the prediction of earthquake-induced landslide displacement and enhance our understanding of dynamic triggering.

Numerical simulations of granular dam break: Comparison between discrete element, Navier-Stokes, and thin-layer models

Granular flows occur in various contexts, including laboratory experiments, industrial processes, and natural geophysical flows. To investigate their dynamics, different kinds of physically based models have been developed. These models can be characterized by the length scale at which dynamic processes are described. Discrete models use a microscopic scale to individually model each grain, Navier-Stokes models use a mesoscopic scale to consider elementary volumes of grains, and thin-layer models use a macroscopic scale to model the dynamics of elementary columns of fluids. In each case, the derivation of the associated equations is well-known. However, few studies focus on the extent to which these modeling solutions yield mutually coherent results. In this article, we compare the simulations of a granular dam break on a horizontal or inclined planes for the discrete model convex optimization contact dynamics (COCD), the Navier-Stokes model Basilisk, and the thin-layer depth-averaged model SHALTOP. We show that, although all three models allow reproducing the temporal evolution of the free surface in the horizontal case (except for SHALTOP at the initiation), the modeled flow dynamics are significantly different, and, in particular, during the stopping phase. The stresses measured at the flow's bottom, reflecting the flow dynamics, are in relatively good agreement, but significant variations are obtained with the COCD model due to complex and fast-varying granular lattices. Similar conclusions are drawn using the same rheological parameters to model a granular dam break on an inclined plane. This comparison exercise is essential for assessing the limits and uncertainties of granular flow modeling.