skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Deformation and Frictional Failure of Granular Media in 3D Analog and Numerical Experiments
Abstract Frictional sliding along grain boundaries in brittle shear zones can result in the fragmentation of individual grains, which ultimately can impact slip dynamics. During deformation at small scales, stick–slip motion can occur between grains when existing force chains break due to grain rearrangement or failure, resulting in frictional sliding of granular material. The rearrangement of the grains leads to dilation of the granular package, reducing the shear stress and subsequently leading to slip. Here, we conduct physical experiments employing HydroOrbs, an elasto-plastic material, to investigate grain comminution in granular media under simple shear conditions. Our findings demonstrate that the degree of grain comminution is dependent on both the normal force and the size of the grains. Using the experimental setup, we benchmark Discrete Element Method (DEM) numerical models, which are capable of simulating the movement, rotation, and fracturing of elasto-plastic grains subjected to simple shear. The DEM models successfully replicate both grain comminution patterns and horizontal force fluctuations observed in our physical experiments. They show that increasing normal forces correlate with higher horizontal forces and more fractured grains. The ability of our DEM models to accurately reproduce experimental results opens up new avenues for investigating various parameter spaces that may not be accessible through traditional laboratory experiments, for example, in assessing how internal friction or cohesion affect deformation in granular systems.  more » « less
Award ID(s):
1843676
PAR ID:
10502454
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Pure and Applied Geophysics
Volume:
181
Issue:
7
ISSN:
0033-4553
Format(s):
Medium: X Size: p. 2083-2105
Size(s):
p. 2083-2105
Sponsoring Org:
National Science Foundation
More Like this
  1. Deformation experiments on hematite characterize its slip‐rate dependent frictional properties and deformation mechanisms. These data inform interpretations of slip behavior from exhumed hematite‐coated faults and present‐day deformation at depth. We used a rotary‐shear apparatus to conduct single‐velocity and velocity‐step experiments on polycrystalline specular hematite rock (∼17 μm average plate thickness) at slip rates of 0.85 μm/s to 320 mm/s, displacements of primarily 1–3 cm and up to 45 cm, and normal stresses of 5 and 8.5 MPa. The average coefficient of friction is 0.70; velocity‐step experiments indicate velocity‐strengthening to velocity‐neutral behavior at rates <1 mm/s. Scanning electron microscopy showed experimentally generated faults develop in a semi‐continuous, thin layer of red hematite gouge. Angular gouge particles have an average diameter of ∼0.7 μm, and grain size reduction during slip yields a factor of 10–100 increase in surface area. Hematite is amenable to (U‐Th)/He thermochronometry, which can quantify fault‐related thermal and mechanical processes. Comparison of hematite (U‐Th)/He dates from the undeformed material and experimentally produced gouge indicates He loss occurs during comminution at slow deformation rates without an associated temperature rise required for diffusive loss. Our results imply that, in natural fault rocks, deformation localizes within coarse‐grained hematite by stable sliding, and that hematite (U‐Th)/He dates acquired from ultracataclasite or highly comminuted gouge reflect minor He loss unrelated to thermal processes. Consequently, the magnitude of temperature rise and associated thermal resetting in hematite‐bearing fault rocks based on (U‐Th)/He thermochronometry may be overestimated if only diffusive loss of He is considered. 
    more » « less
  2. Abstract Localized frictional sliding on faults in the continental crust transitions at depth to distributed deformation in viscous shear zones. This brittle‐ductile transition (BDT), and/or the transition from velocity‐weakening (VW) to velocity‐strengthening (VS) friction, are controlled by the lithospheric thermal structure and composition. Here, we investigate these transitions, and their effect on the depth extent of earthquakes, using 2D antiplane shear simulations of a strike‐slip fault with rate‐and‐state friction. The off‐fault material is viscoelastic, with temperature‐dependent dislocation creep. We solve the heat equation for temperature, accounting for frictional and viscous shear heating that creates a thermal anomaly relative to the ambient geotherm which reduces viscosity and facilitates viscous flow. We explore several geotherms and effective normal stress distributions (by changing pore pressure), quantifying the thermal anomaly, seismic and aseismic slip, and the transition from frictional sliding to viscous flow. The thermal anomaly can reach several hundred degrees below the seismogenic zone in models with hydrostatic pressure but is smaller for higher pressure (and these high‐pressure models are most consistent with San Andreas Fault heat flow constraints). Shear heating raises the BDT, sometimes to where it limits rupture depth rather than the frictional VW‐to‐VS transition. Our thermomechanical modeling framework can be used to evaluate lithospheric rheology and thermal models through predictions of earthquake ruptures, postseismic and interseismic crustal deformation, heat flow, and the geological structures that reflect the complex deformation beneath faults. 
    more » « less
  3. We report on an open-source, publicly available C++ software module called Chrono::GPU, which uses the Discrete Element Method (DEM) to simulate large granular systems on Graphics Processing Unit (GPU) cards. The solver supports the integration of granular material with geometries defined by triangle meshes, as well as co-simulation with the multi-physics simulation engine Chrono. Chrono::GPU adopts a smooth contact formulation and implements various common contact force models, such as the Hertzian model for normal force and the Mindlin friction force model, which takes into account the history of tangential displacement, rolling frictional torques, and cohesion. We report on the code structure and highlight its use of mixed data types for reducing the memory footprint and increasing simulation speed. We discuss several validation tests (wave propagation, rotating drum, direct shear test, crater test) that compare the simulation results against experimental data or results reported in the literature. In another benchmark test, we demonstrate linear scaling with a problem size up to the GPU memory capacity; specifically, for systems with 130 million DEM elements. The simulation infrastructure is demonstrated in conjunction with simulations of the NASA Curiosity rover, which is currently active on Mars. 
    more » « less
  4. Abstract The operation of fracture, diffusion, and intracrystalline‐plastic micromechanisms during semibrittle deformation of rock is directly relevant to understanding mechanical behavior across the brittle‐plastic transition in the crust. An outstanding question is whether (1) the micromechanisms of semibrittle flow can be considered to operate independently, as represented in typical crustal strength profiles across the brittle to plastic transition, or (2) the micromechanisms are coupled such that the transition is represented by a distinct rheology with dependency on effective pressure, temperature, and strain rate. We employ triaxial stress‐cycling experiments to investigate elastic‐plastic and viscoelastic behaviors during semibrittle flow in two distinctly different monomineralic, polycrystalline, synthetic salt‐rocks. During semibrittle flow at high differential stress, granular, low‐porosity, work‐hardened salt‐rocks deform predominantly by grain‐boundary sliding and wing‐crack opening accompanied by minor intragranular dislocation glide. In contrast, fully annealed, near‐zero porosity salt‐rocks flow at lower differential stress by intragranular dislocation glide accompanied by grain‐boundary sliding and opening. Grain‐boundary sliding is frictional during semibrittle flow at higher strain rates, but the associated dispersal of water from fluid inclusions along boundaries can activate fluid‐assisted diffusional sliding at lower strain rates. Changes in elastic properties with semibrittle flow largely reflect activation of sliding along closed grain boundaries. Observed microstructures, pronounced hysteresis and anelasticity during cyclic stressing after semibrittle flow, and stress relaxation behaviors indicate coupled operation of micromechanisms leading to a distinct rheology (hypothesis 2 above). 
    more » « less
  5. Abstract Faults can release energy via a variety of different slip mechanisms ranging from steady creep to fast and destructive earthquakes. Tying the rheology of the crust to various slip dynamics is important for our understanding of plate tectonics and earthquake generation. Here we propose that the interplay of fractures and viscous flow leads to a spectrum between stick‐slip and creep. We use an elasto‐visco‐plastic rock analog (Carbopol U‐21) where we vary the yield stress to investigate its impact on slip dynamics in shear experiments. The experiments are performed using a simple shear apparatus, which provides distributed shear across the entire width of the experiment and allows in situ observations of deformation. We record force and displacement during deformation and use time lapse photography to document fracture development. A low yield stress (25 Pa) leads to creep dynamics in the absence of fractures. An intermediate yield stress (144 Pa) leads to the development and interaction of opening (mode I) and shear (mode II) fractures. This interaction leads to a spectrum in slip dynamics ranging from creep to stick‐slip. A high yield stress (357 Pa) results in the development of many mode I fractures and a deformation signal dominated by stick‐slip. These results show that bulk yield stress, fracture formation, and slip dynamics are closely linked and can lead to a continuum between creep and stick‐slip. We suggest that rheology should be considered as an additional mechanism to explain the broad range of slip dynamics in natural faults. 
    more » « less