Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shear Zone”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Enhancing NDE Reliability for Grade 91 Steel Welds: Ultrasonic Imaging and Microstructural Correlations

Ensuring the integrity of Grade 91 (9Cr-1Mo-V) steel welds is vital for the safe and reliable operation of fossil fuel–fired and nuclear power plants. Here, this study applies an imaging technique for the ultrasonic characterization of two Grade 91 steel welds created with cold metal transfer and flux-cored arc welding processes. Ultrasonic immersion testing in the through-transmission configuration was employed to generate shear waves, which helped identify the weld metal, heat-affected zone, and base metal regions. These weld microstructures were also correlated to their ultrasonic images using metallography, ultrasonic amplitude, hardness measurements, and grain size. The findings from this study can assist practitioners in developing new nondestructive evaluation technologies, improving the inspection reliability of creep strength–enhanced ferritic steel welds by potentially identifying weld microstructure regions susceptible to creep-type failures.

Grade 91 steel↗

Analysis of Fault Zone Resonance Modes Recorded by a Dense Seismic Array Across the San Jacinto Fault Zone at Blackburn Saddle

We present observations and modeling of spatial eigen-functions of resonating waves within fault zone waveguide, using data recorded on a dense seismic array across the San Jacinto Fault Zone (SJFZ) in southern California. The array consists of 5-Hz geophones that cross the SJFZ with ~10-30 m spacing at the Blackburn Saddle near the Hemet Stepover. Wavefield snapshots after the S wave arrival are consistent for more than 50 near-fault events, suggesting that this pattern is controlled by the fault zone structure rather than source properties. Data from example event with high signal to noise ratio shows three main frequency peaks at ~1.3 Hz, ~2.0 Hz, and ~2.8 Hz in the amplitude spectra of resonance waves averaged over stations near the fault. The data are modeled with analytical expressions for eigen-functions of resonance waves in a low-velocity layer (fault zone) between two quarter-spaces. Using a grid search-based method, we investigate the possible width of the waveguide, location within the array, and shear wave velocities of the media that fit well the resonance signal at ~1.3 Hz. The results indicate a ~300 m wide damaged fault zone layer with ~65% S wave velocity reduction compared to the host rock. The SW edge of the low-velocity zone is near the mapped fault surface trace, indicating that the damage zone is asymmetrically located at the regionally faster NE crustal block. The imaging resolution of the fault zone structure can be improved by modeling fault zone resonance modes and trapped waves together.

58 GEOSCIENCES↗

Graph-learning approach to combine multiresolution seismic velocity models

SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.

Geochemistry & Geophysics↗

Material Characterization of Grade 91 Steel Welds using, Micro-resolution Ultrasonic Imaging System

Grade 91 (9Cr-1Mo-V) steel is a common structural material for fossil-fuel fired and nuclear power plants because of its creep strength and oxidation resistance at high temperatures. In this research, we studied the feasibility of correlating the high-resolution images obtained from the immersion ultrasonic imaging technique to microstructure regions in Grade 91 steel plate weldments. The high-resolution ultrasonic images presented in this paper were generated using a 20-MHz focused transducer as the transmitter, and a laser vibrometer as the receiver in the through transmission configuration to produce longitudinal and shear waves. This nondestructive evaluation (NDE) method was used to identify the Fusion Zone (FZ), Heat-Affected Zone (HAZ), and Base Metal (BM) regions of Grade 91 cold metal transfer (CMT) and flux core arc welds (FCAW). The ultrasonic imaging results were analyzed using beam propagation amplitudes and homogeneity, which in turn were used to relate the different microstructure regions of the FCAW and CMT welded samples. With this study, we were able to identify highly susceptible microstructures in the weld that could potentially promote creep type failures.

20 FOSSIL-FUELED POWER PLANTS↗

Coupled Thermo–Hydro–Mechanical Modeling of Hydro-Shearing Stimulation in an Enhanced Geothermal System in the Raft River Geothermal Field, USA

The Raft River Enhanced Geothermal System (EGS) demonstration project aims to improve the geothermal production by enhancing the reservoir permeability via shear stimulation in well RRG-9. In this paper, we performed a series of 3D thermo– hydro–mechanical (THM) simulations to investigate the influence of hydraulic and thermal effects on the development of this EGS. The model includes synthetic fracture populations based on borehole televiewer images and in-situ stress measurements from well RRG-9. Fracture permeability evolution is determined using an empirical permeability law developed from laboratory experiments. The model was calibrated by comparing the hydraulic response to field observations, including wellhead pressure, injection rate, and well injectivity. Particularly, we analyzed the enhancement of reservoir permeability and the spatial extent of the stimulation zone for the given injection schedule. Our results indicate that the permeability enhancement of fractured geothermal reservoir is caused by the combined effects of injection-induced cooling and fluid pressure increase. The decrease of temperature plays a dominant role in reactivation of natural fractures under the hydro-shearing mechanism, while the higher injection pressure promotes shear failure and enlarges the stimulation zone. Finally, for the specified extensional stress state, the model favors greatest permeability enhancement along the maximum horizontal principal stress, a moderate vertical enhancement, and a smaller gain along the minimum horizontal principal stress.

58 GEOSCIENCES↗

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗

The High-Frequency Signature of Slow and Fast Laboratory Earthquakes

Tectonic faults fail through a spectrum of slip modes, ranging from slow aseismic creep to rapid slip during earthquakes. Understanding the seismic radiation emitted during these slip modes is key for advancing earthquake science and earthquake hazard assessment. In this work, we use laboratory friction experiments instrumented with ultrasonic sensors to document the seismic radiation properties of slow and fast laboratory earthquakes. Stick-slip experiments were conducted at a constant loading rate of 8 μm/s and the normal stress was systematically increased from 7 to 15 MPa. We produced a full spectrum of slip modes by modulating the loading stiffness in tandem with the fault zone normal stress. Acoustic emission data were recorded continuously at 5 MHz. We demonstrate that the full continuum of slip modes radiate measurable high-frequency energy between 100 and 500 kHz, including the slowest events that have peak fault slip rates <100 μm/s. The peak amplitude of the high-frequency time-domain signals scales systematically with fault slip velocity. Stable sliding experiments further support the connection between fault slip rate and high-frequency radiation. Experiments demonstrate that the origin of the high-frequency energy is fundamentally linked to changes in fault slip rate, shear strain, and breaking of contact junctions within the fault gouge. Our results suggest that having measurements close to the fault zone may be key for documenting seismic radiation properties and fully understanding the connection between different slip modes.

58 GEOSCIENCES↗

Strain-rate dependent deformation mechanisms in single-layered Cu, Mo, and multilayer Cu/Mo thin films

Here, strain-rate sensitivity and rate-dependent hardness, over a range of 10 -2 to 10 2 s -1 , of sputter-deposited single-layered Cu, Mo, and 5 nm Cu/ 5 nm Mo, and 100 nm Cu/ 100 nm Mo multilayer films with a total film thickness of 5 μm were measured using nanoindentation. The plastic zone underneath the nanoindents was characterized via cross-sectional transmission electron microscopy (XTEM). The multilayer films exhibited enhanced hardness but slightly reduced strain-rate sensitivity with decreasing layer thickness from 100 nm to 5 nm. Only the 5 nm Cu/ 5 nm Mo multilayer film exhibited shear bands underneath the nanoindents, and the size of the shear bands increased with increasing strain rate. In contrast, the 100 nm Cu/ 100 nm Mo multilayer film exhibited material pile-up around the indents and significant nanotwinning within Cu grains. The effect of strain rate and layer thickness on the hardness and strain rate sensitivity of the multilayer thin films is interpreted using a modified confined layer slip (CLS) model. The reduced rate sensitivity at 5 nm as compared to 100 nm correlates with abundant growth nanotwins in the Cu grains in 100 nm and formation of shear bands in 5 nm multilayers. In single layer films, a substructure with a high density of dislocations was observed consistent with the plastic strain gradient in the indent plastic zone. No evidence of deformation twins was noted in any of the samples.

36 MATERIALS SCIENCE↗

A new mechanism for filament disconnection at the X-point: poloidal shear in radial E × B velocity

Plasma filaments generated by turbulence in the scrape-off layer are sometimes observed in experiments and simulations to disconnect from the target plate in the vicinity of the X-point, resulting in a ‘quiescent zone’ with reduced fluctuations, possibly affecting the fluctuating part of the divertor heat loads. This phenomena is usually explained by the flux tubes squeezing, induced by flux expansion, down to scales where collisions dominate. In this work, we consider an additional mechanism spontaneously arising at the X-point: the poloidal shear in radial E × B velocity. Through a simple model, tested against 3D global turbulence simulations, herein we quantify the effectiveness of this disconnection mechanism and identify ways to control it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control Mechanisms for Self‐Sealing in Activated Clay‐Rich Faults Through Controlled Hydraulic Injection Experiment

Abstract In a high‐pressure injection fault activation experiment conducted at the Mont Terri underground research laboratory in Switzerland, the transmissivity of the Opalinus Clay fault significantly increased due to opening and shearing. The fluid injection, spanning a few hours, generated a 10 m radius fault activation patch. Subsequent pressure pulse tests conducted bi‐weekly for a year revealed the gradual return of fault transmissivity to its initial state. The study utilized fluid pressure decay analysis, optical fiber monitoring, continuous active source seismic measurements and borehole displacement sensors for measuring fault displacements. The fault zone exhibited a dilation of approximately 1.4 mm, associated with both normal and tangential movements during activation, resulting in a sudden transmissivity increase from 1 × 10 −12 to 3.2 × 10 −7 m 2 /s. Early post‐activation, transient compaction and the subsequent slow compaction were observed, transitioning to an extension regime. The pressure pulse tests demonstrated a rapid transmissivity drop by more than two orders of magnitude within the first 10 days, followed by a gradual and less pronounced decrease. Plastic shear and compaction dominated the transmissivity evolution until 70 days after injection ended, followed by a period where additional factors, such as clay mineral swelling, influenced the behavior. Extrapolation suggested a sealing process taking at least 50 years after the initial activation. Plain Language Summary A field‐scale fault activation experiment offers valuable insights into the elasto‐plastic processes governing the sealing of shale faults. The experiment reveals a rapid increase in the fault's transmissivity by approximately five orders of magnitude during activation. Subsequent observations show a gradual transmissivity decrease by about three orders of magnitude post‐activation, with slow long‐term plastic shear and compaction of the fault competing against secondary processes, notably clay mineral swelling. All conceptual models employed to interpret these field data converge on the estimation that the fault's return to its initial low transmissivity state would require a minimum of 50 years. Key Points High‐pressure injection fault activation experiment at the Mont Terri underground research laboratory Continuous transmissivity measurements record self‐sealing inside a clay‐rich fault zone Transmissivity undergoes a phase of domination by slow plastic compaction and shearing during the initial post‐activation period, with mineral swelling exerting its influence over the long term

Guglielmi, Yves↗

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗

Machine Learning for Joint Quality Control

The use of lightweight material combinations has been highly demanded in manufacturing automotive structures. However, making robust dissimilar material joints of such lightweight materials is still challenging. A significant barrier to achieving high-quality and repeatable joint performance is a deficient understanding of the relationship between the welding process, joint attributes, and joint performance. In this context, welding factors refer to material, equipment, environment, and process parameters, while joint features comprise specific microstructural attributes of the weld such as nugget size, heat affected zone (HAZ) topology, intermetallic layer thickness, and sheet thickness reduction. Joint performance is quantified in terms of strength (e.g., tensile shear, coach peel, cross-tension), weld size, and hardness, among other factors. While there have been many attempts to establish this process-structure-property relationship by developing a model derived from the associated physics and first principles, the complexity of the joining processes compounded by the complex interactions with different materials in an automotive assembly line environment, has hindered the usefulness of such attempts. The complexity is further exacerbated using different stacking materials, especially comprising dissimilar material combinations. In practice, the common approach has been the laborious process of creating welds, characterizing them, and then physically testing them through experimentation. With the emergence of artificial intelligence (AI) methods, an alternative pathway to eliciting the desired process-structure-property relationship at an accelerated pace is to use a data-driven approach by employing machine-learning (ML) techniques. This approach is benefitted by the availability of large streams of data, generated through years of research and testing by original equipment manufacturers, in the form of material, process, environmental, equipment, microstructural, and bulk-scale performance information from multimodal, multiscale sensors making measurements from laboratory-scale to production-scale processes. During Phase I efforts, which ended in fiscal year (FY) 2021, the Oak Ridge National Laboratory and Pacific Northwest National Laboratory (ORNL/PNNL) team demonstrated the effectiveness of different ML/AI frameworks in modeling complex relationships between resistance spot welding (RSW) process parameters, weld attributes, and joint properties using a subset of data from General Motors (GM). In FY 2022, the project team further refined and expanded their respective ML models to analyze additional welds with new weld stack-ups and materials to enhance the ML model predictive capability. ORNL extended its unified deep neural networks (DNN) ML training and prediction framework with new data streams of process parameters, and PNNL extended its model describing RSW process parameters’ associations with weld attributes. In FY 2023, the project team completed the development of the AI/ML architecture for analyzing aluminum/steel joints manufactured by GM via RSW and transitioned into the inline welding quality monitoring task for steel/steel RSW joints provided by GM.

36 MATERIALS SCIENCE↗

Elasticity of davemaoite as a primary contributor to lower-mantle heterogeneities

Geophysical detection of subducted mid–ocean ridge basalt (MORB) in the lower mantle is hindered by uncertainties in the elasticity of Fe,Al,Mg,Ti–bearing davemaoite, a key MORB component. Using Brillouin spectroscopy and x-ray diffraction, we determined the elasticity of a Ca 0.906(1) Fe 2+ 0.027(1) Fe 3+ 0.042(1) Mg 0.033(1) Al 0.072(1) Ti 0.020(1) Si 0.912(1) O 3 davemaoite up to 113 gigapascals and 2294 K. Here, we found that it exhibited a shear wave velocity 10 to 20% slower than end-member davemaoite, making it the slowest phase among major lower-mantle minerals. Our models show that MORB, containing 20 to 25 volume percent davemaoite, potentially contributes to large low-shear-velocity provinces (LLSVPs), whereas a cumulate layer enriched in davemaoite crystallized from basal magma ocean may comprise ultralow-velocity zones (ULVZs). Davemaoite’s ability to host incompatible and heat-producing elements possibly links LLSVPs and ULVZs to mantle plume initiation and geochemical signatures of ocean island basalts.

Zhou, Wen-Yi [Texas A & M Univ., College Station, ↗

Numerical and Experimental Investigations on the Ignition Behavior of OME

On the path towards climate-neutral future mobility, the usage of synthetic fuels derived from renewable power sources, so-called e-fuels, will be necessary. Oxygenated e-fuels, which contain oxygen in their chemical structure, not only have the potential to realize a climate-neutral powertrain, but also to burn more cleanly in terms of soot formation. Polyoxymethylene dimethyl ethers (PODE or OMEs) are a frequently discussed representative of such combustibles. However, to operate compression ignition engines with these fuels achieving maximum efficiency and minimum emissions, the physical-chemical behavior of OMEs needs to be understood and quantified. Especially the detailed characterization of physical and chemical properties of the spray is of utmost importance for the optimization of the injection and the mixture formation process. The presented work aimed to develop a comprehensive CFD model to specify the differences between OMEs and dodecane, which served as a reference diesel-like fuel, with regards to spray atomization, mixing and auto-ignition for single- and multi-injection patterns. The simulation results were validated against experimental data from a high-temperature and high-pressure combustion vessel. The sprays’ liquid and vapor phase penetration were measured with Mie-scattering and schlieren-imaging as well as diffuse back illumination and Rayleigh-scattering for both fuels. To characterize the ignition process and the flame propagation, measurements of the OH* chemiluminescence of the flame were carried out. Significant differences in the ignition behavior between OMEs and dodecane could be identified in both experiments and CFD simulations. Liquid penetration as well as flame lift-off length are shown to be consistently longer for OMEs. Zones of high reaction activity differ substantially for the two fuels: Along the spray center axis for OMEs and at the shear boundary layers of fuel and ambient air for dodecane. Additionally, the transient behavior of high temperature reactions for OME is predicted to be much faster.

33 ADVANCED PROPULSION SYSTEMS↗

Environmental sensitivities of shallow-cumulus dilution – Part 2: Vertical wind profile

This second part of a numerical study on shallow-cumulus dilution focuses on the sensitivity of cloud dilution to changes in the vertical wind profile. Insights are obtained through large-eddy simulations of maritime and continental cloud fields. In these simulations, the speed of the initially uniform geostrophic wind and the strength of geostrophic vertical wind shear in the cloud and subcloud layer are varied. Increases in the cloud-layer vertical wind shear (up to 9 ms -1 km -1 ) lead to 40%–50% larger cloud-core dilution rates compared to their respective unsheared counterparts. When the background wind speed, on the other hand, is enhanced by up to 10 ms -1 and subcloud-layer vertical wind shear develops or is initially prescribed, the dilution rate decreases by up to 25%. The sensitivities of the dilution rate are linked to the updraft strength and the properties of the entrained air. Increases in the wind speed or vertical wind shear result in lower vertical velocities across all sets of experiments with stronger reductions in the cloud-layer wind shear simulation (27%–47%). Weaker updrafts are exposed to mixing with the drier surrounding air for a longer time period, allowing more entrainment to occur (i.e., the “core-exposure effect”). However, reduced vertical velocities, in concert with increased cloud-layer turbulence, also assist in widening the humid shell surrounding the cloud cores, leading to entrainment of more humid air (i.e., the “core–shell dilution effect”). In the experiments with cloud-layer vertical wind shear, the core-exposure effect dominates and the cloud-core dilution increases with increasing shear. Conversely, when the wind speed is increased and subcloud-layer vertical wind shear develops or is imposed, the core–shell dilution effect dominates to induce a buffering effect. The sensitivities are generally stronger in the maritime simulations, where weaker sensible heat fluxes lead to narrower, more tilted, and, therefore, more suppressed cumuli when cloud-layer shear is imposed. Moreover, in the experiments with subcloud wind shear, the weaker baseline turbulence in the maritime case allows for a larger turbulence enhancement, resulting in a widening of the transition zones between the cores and their environment, leading to the entrainment of more humid air.

54 ENVIRONMENTAL SCIENCES↗

Gouge formation and dilation impacts to flow during fracture shearing

We report flow in fractures is of critical importance to fluid flow in low-permeability shales, particularly when they restrict vertical migration from reservoirs. The Mt. Simon Sandstone is a target for geologic CO 2 sequestration and the overlying Eau Claire Formation is the sealing formation capping this reservoir. The lower Eau Claire, a primary barrier to upward migration of CO 2 from the Mt. Simon, is lithologically heterogenous with a mixture of sandstone and shale. We use a series of flow-through experiments on fractured samples from the Eau Claire to document the morphological changes taking place in the fracture and surrounding rock as it is sheared in a stepwise fashion. The samples were imaged using computed tomography after each shearing episode and flow resistance was calculated based on the pressure differential across each core. We found that gouge accumulations in the fracture plane led to specific instances of self-sealing behavior despite the fact that most samples increased in average fracture aperture and transmissivity. Gouge limited the increase or, in some cases, decreased flow through these samples. Fractures adjacent to shale-rich zones and concordant with bedding were more likely to show evidence of self-sealing, while increased sandstone content and fracture orientation perpendicular to bedding planes was more likely to result in fracture dilation and fracture permeability increases. Reactivation and shearing of vertical fractures in lithologically heterogenous caprocks is more likely to pose risks to long term CO 2 storage scenarios.

58 GEOSCIENCES↗

Acoustic Energy Release During the Laboratory Seismic Cycle: Insights on Laboratory Earthquake Precursors and Prediction

Abstract Machine learning can predict the timing and magnitude of laboratory earthquakes using statistics of acoustic emissions. The evolution of acoustic energy is critical for lab earthquake prediction; however, the connections between acoustic energy and fault zone processes leading to failure are poorly understood. Here, we document in detail the temporal evolution of acoustic energy during the laboratory seismic cycle. We report on friction experiments for a range of shearing velocities, normal stresses, and granular particle sizes. Acoustic emission data are recorded continuously throughout shear using broadband piezo‐ceramic sensors. The coseismic acoustic energy release scales directly with stress drop and is consistent with concepts of frictional contact mechanics and time‐dependent fault healing. Experiments conducted with larger grains (10.5 μm) show that the temporal evolution of acoustic energy scales directly with fault slip rate. In particular, the acoustic energy is low when the fault is locked and increases to a maximum during coseismic failure. Data from traditional slide‐hold‐slide friction tests confirm that acoustic energy release is closely linked to fault slip rate. Furthermore, variations in the true contact area of fault zone particles play a key role in the generation of acoustic energy. Our data show that acoustic radiation is related primarily to breaking/sliding of frictional contact junctions, which suggests that machine learning‐based laboratory earthquake prediction derives from frictional weakening processes that begin very early in the seismic cycle and well before macroscopic failure.

58 GEOSCIENCES↗

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗