Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock failure”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Time-Dependent Directional Damage Theory for Brittle Rocks Considering the Kinetics of Microcrack Growth

In this work, a novel microcrack damage theory for describing the time-dependent behavior of brittle rocks is proposed. Instead of using scalar or tensorial variables to approximate the effect of microcracks, the concept of directional damage is introduced by upscaling the directional distribution density of microcracks through the noninteracting homogenization scheme. The contribution of an individual crack on global strain is first derived from the crack-opening displacements including sliding, closure, and dilation. The model is then cast in the free energy form, and the range of parameters ensuring compliance with the second law of thermodynamics is derived rigorously. Considering that the subcritical propagation of microcracks is the dominant mechanism behind the time-dependent deformation and failure of brittle rocks, the driving force for directional damage is defined as the corresponding energy release rate, and the effect of time is introduced through a damage evolution law inspired from the kinetics of subcritical crack growth. The model is applied to simulate the behavior of basalt rock under various loading conditions. In addition to predicting the stress–strain–time response, the model offers enhanced resolution in representing the anisotropic characteristics of microcrack-induced damage in brittle rocks and their time-dependent evolutions.

58 GEOSCIENCES↗
DOE OSTI2022

Evolution of tensile and shear cracking in crystalline rocks under compression

Here we performed laboratory-scale unconfined compression experiments on Barre granite specimen with a single pre-existing flaw to study the evolution of tensile and shear cracking that occurs during different stress levels. Moment tensor inversion of acoustic emission (AE) events in combination with the nonelastic strain component obtained through the 2D digital image correlation (DIC) technique are used to track the changes in the source mechanisms of the stress-induced cracks. The mode of deformation computed from the image based strain profiles, enables visual comparison of the nucleation, growth and interaction of the tensile and shear cracks with the microseismic source mechanism observed by moment tensor inversion of the AE events. Comparing the results obtained from the two techniques on the same dataset demonstrates a quantitative correlation between the acoustic and visual observations in terms of the cracking mechanisms obtained at different stages of cracking from crack initiation at the flaw tips to the failure of the rock specimen.

58 GEOSCIENCES↗
DOE OSTI2020

Mjölnir: a miniature triaxial rock deformation apparatus for 4D synchrotron X-ray microtomography

An X-ray transparent experimental triaxial rock deformation apparatus, here named 'Mjölnir', enables investigations of brittle-style rock deformation and failure, as well as coupled thermal, chemical and mechanical processes relevant to a range of Earth subsurface environments. Designed to operate with cylindrical samples up to 3.2 mm outside-diameter and up to 10 mm length, Mjölnir can attain up to 50 MPa confining pressure and in excess of 600 MPa axial load. The addition of heaters extends the experimental range to temperatures up to 140°C. Deployment of Mjolnir on synchrotron beamlines indicates that full 3D datasets may be acquired in a few seconds to a few minutes, meaning full 4D investigations of deformation processes can be undertaken. Mjölnir is constructed from readily available materials and components and complete technical drawings are included in the supporting information.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Predicting the proximity to macroscopic failure using local strain populations from dynamic in situ X-ray tomography triaxial compression experiments on rocks

Predicting the proximity of large-scale dynamic failure is a critical concern in the engineering and geophysical sciences. Here we use evolving contractive, dilatational, and shear strain deformation preceding failure in dynamic X-ray tomography experiments to examine which strain components best predict the proximity to failure. We develop machine learning models to predict the proximity to failure using time series of three-dimensional local incremental strain tensor fields acquired in rock deformation experiments under stress conditions of the upper crust. Three-dimensional scans acquired in situ throughout triaxial compression experiments provide a distribution of density contrasts from which we estimate the three-dimensional incremental strain that accumulates between each scan acquisition. Training machine learning models on multiple experiments of six rock types provides suites of feature importance that indicate the predictive power of each feature. Comparing the average importance of groups of features that include information about each strain component quantifies the ability of the contractive, dilatational and shear strain to predict the proximity of macroscopic failure. A total of 24 models of four machine learning algorithms with six rock types indicate that 1) the dilatational strain provides the best predictive power of the strain components, and 2) the intermediate values (25th-75th percentile) of the strain population provide the best predictive power of the statistics of the strain populations. In addition, the success of the predictions of models trained on one rock type and tested on other rock types quantifies the similarities and differences of the precursory strain accumulation process in the six rock types. These similarities suggest the potential existence of a unified theory of brittle rock deformation for a range of rock types.

58 GEOSCIENCES↗
DOE OSTI2020

Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method

Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

58 GEOSCIENCES↗
DOE OSTI2021

Attention Network Forecasts Time-to-Failure in Laboratory Shear Experiments

Rocks under stress deform by creep mechanisms that include formation and slip on small-scale internal cracks. Intragranular cracks and slip along grain contacts release energy as elastic waves termed acoustic emissions (AE). AEs are thought to contain predictive information that can be used for fault failure forecasting. Here, we present a method using unsupervised classification and an attention network to forecast labquakes using AE waveform features. Our data were generated in a laboratory setting using a biaxial shearing device with granular fault gouge intended to mimic the conditions of tectonic faults. Here, we analyzed the temporal evolution of AEs generated throughout several hundred laboratory earthquake cycles. We used a Conscience Self-Organizing Map (CSOM) to perform topologically ordered vector quantization based on waveform properties. The resulting map was used to interactively cluster AEs. We examined the clusters over time to identify those with predictive ability. Finally, we used a variety of LSTM and attention-based networks to test the predictive power of the AE clusters. By tracking cumulative waveform features over the seismic cycle, the network is able to forecast the time-to-failure (TTF) of lab earthquakes. Our results show that analyzing the data to isolate predictive signals and using a more sophisticated network architecture are key to robustly forecasting labquakes. In the future, this method could be applied on tectonic faults to monitor earthquakes and augment early warning systems.

58 GEOSCIENCES↗
DOE OSTI2024

Intermediate Principal Stress Effects on the 3D Cracking Behavior of Flawed Rocks Under True Triaxial Compression

Crack initiation, growth, and coalescence in flawed rocks have been extensively studied for 2D (planar, penetrating) flaws under uniaxial/biaxial compression. However, little is known as to the mechanisms and processes of cracking from 3D flaws under true triaxial compression, where the intermediate principal stress (σ 2 ) is distinguished from the major and minor principal stresses. In this work, we systematically investigate the effects of σ 2 on the 3D cracking behavior of rock specimens with preexisting flaws, through the use of mechanistic simulations of mixed-mode fracture in rocks. We explore how two characteristics of σ 2 , namely, (i) its orientation with respect to the flaw and (ii) its magnitude, affect two aspects of the cracking behavior, namely, (i) the cracking pattern and (ii) the peak stress. Results show that the orientation of σ 2 exerts more control over the cracking pattern than the flaw inclination angle. The peak stress becomes highest when σ 1 is parallel to the flaw, whereas it becomes lowest when σ 2 is parallel to the flaw. Also, the effects of σ 2 magnitude are more significant when σ 2 becomes more oblique to the flaw plane. On the basis of our observations, we propose mechanisms underlying the cracking behavior of 3D flawed rocks under true triaxial compression.

58 GEOSCIENCES↗
DOE OSTI2023

Material Property Determinations for Rock Valley Legacy Samples in Support of Source Physics Experiment

The Source Physics Experiment (SPE) is a long-term NNSA research and development effort designed to improve nonproliferation verification and monitoring capabilities. The overarching goals of the SPE program are to improve understanding of prompt signals and physical signatures that develop from underground chemical explosions and associated modeling capabilities. Our work focuses on a primary factor controlling chemical explosion induced signals and signatures: the material properties of the rocks in which the chemical explosion takes place. This document reports on material property determinations of legacy core USGS Test Well F and outcrop analogs for the subsurface stratigraphy for the third phase of SPE in the Rock Valley (RV) area of the Nevada National Security Site (NNSS). The objective of this work is to establish a baseline set of lithologic descriptions and material properties expected prior to observatory borehole drilling in support of the SPE-RVDC (Rock Valley Direct Comparison) experiment. We determine for each rock type the compressional failure envelope, elastic properties as a function of stress (bulk modulus versus mean stress, shear modulus versus shear stress, Young’s modulus versus axial stress and Poisson’s ratio versus axial stress), indirect tensile strength, and porosity. Geologic characterization, both at the core-scale and microscale, provides context for using the data in modeling efforts and to inform interpretations for the material properties testing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A scale-invariant bonded particle model for simulating large deformation and failure of continua

In this paper, a scale-invariant bond model for simulating two-dimensional large deformation of continua is presented. An accurate torque form caused by relative rotation between bonded particles is derived based on beam theory. Validations of the proposed bond model are carried out by large deformation simulations of a simply supported beam subjected to central loading and a cantilever beam with end loading. Comparison of simulations with analytical solutions show a very good agreement. Later, Brazilian tests are simulated to demonstrate the ability of the framework to model failure and fracture process of rock.

58 GEOSCIENCES↗
DOE OSTI2023

Coupled Hydromechanical Modeling and Assessment of Induced Seismicity at FWU: Utilizing Time-Lapse VSP and Microseismic Data

The objective of this work is to utilize integrated geomechanics, field vertical seismic profile (VSP) and microseismic data to characterize the complex subsurface stress conditions at the Farnsworth Unit (FWU). The model is based on a five-spot sector model extracted from a primary geomechanical model. The five-spot well injection pattern is characterized by extensive reservoir characterization data, such well logs, extracted cores and borehole geophone data, to facilitate the detailed examination of stress changes and microseismic event occurrences. The study utilizes field vertical seismic volumes acquired from the injection well 13-10A. The seismic volumes successfully provided snapshots of the behavior of the reservoir at distinct times. The use of VSP and microseismic data provided direct and indirect estimates of the dynamic stress changes occurring in the overburden, reservoir and underburden rock formations. In order to illuminate the stress regions and identify rocks that have undergone inelastic failure, microseismic event occurrences were utilized. Microseismic activity has been detected at the FWU; further study of its locations, timing, and magnitude was needed to deduce the nature of the changing stress state. The results of the study revealed that microseismic events were successfully modeled within the Morrow B formation. Moment magnitudes of seismic events were within the same magnitudes for events in the reservoir, suggesting the suitability of the model. The results of the study showed that the computed moment magnitudes for seismic events were insignificant to warrant safety concerns. The study findings showed the usefulness of coupled hydromechanical models in predicting the subsurface stress changes associated with CO 2 injection. The knowledge gained from this study will serve as a guideline for industries planning to undertake underground CO 2 storage, and characterize the subsurface stress changes.

02 PETROLEUM↗
DOE OSTI2021

Volumetric and shear strain localization throughout triaxial compression experiments on rocks

Deformation localization is a widely observed, but rarely quantified process in the crust. Recent observations suggest that the localization of seismicity and fracture networks can help identify the approach to catastrophic failure. Here, we quantify the localization processes of the volumetric and deviatoric strain components in twelve triaxial compression experiments imaged with X-ray tomography. We capture three-dimensional images of the rock cores during triaxial compressing toward failure, and then calculate the local strain components using digital volume correlation. The divergence and curl of the incremental displacement vector field provide the volumetric and deviatoric components of the strain field. We quantify localization using the proportion of the rock occupied by high magnitudes of the volumetric and deviatoric strains, and the Gini coefficient of these high magnitude strains, which measures the deviation from a uniform process. We find that the vast majority, but not all, of the experiments experience strain localization toward failure. The rocks typically experience their maximum degree of strain localization not immediately preceding failure, but on average at 90% of the failure stress. The volumetric strain tends to localize more than the deviatoric strain. Finally, these observations support using the localization of the volumetric strain, along with the deviatoric strain, to identify the evolution of the precursory phase preceding earthquakes.

58 GEOSCIENCES↗
DOE OSTI2021

Failure in Confined Brazilian Tests on Sandstone

Strength of rocks in the confined tension region, where the minimum principal stress is tensile, has only infrequently been measured and is not well understood. Quasi-static confined Brazilian tests under a range of confining stresses (2.76 to 27.58 MPa) where used to determine the strength of sandstone in the confined tension region. The test results indicate that the strength in the confined tension region was a strong function of the intermediate principal stress: increasing the intermediate principal stress significantly increased the strength of the sandstone. The strength data were well fit by the Mogi–Coulomb criterion, which accounts for the intermediate principal stress. Unconfined Brazilian strength data were not well fit to the Mogi–Coulomb criterion derived from the confined Brazilian test data, consistent with a transition from tensile to shear processes dominating failure with increasing confining pressure. Observations of post-failure fracture surfaces reveal more indication of shear processes with increasing confining pressure. Numerical simulations from combined finite-discrete element method are compared to the experimental results and reflect similar conditions for failure compared to the experimental tests in the confined tension region.

42 ENGINEERING↗
DOE OSTI2023

Towards commercialization of fluorinated cation-disordered rock-salt Li-ion cathodes

Cation-disordered rock-salt cathodes (DRX) are promising materials that could deliver high capacities (>250 mAh g −1 ) with Earth abundant elements and materials. However, their electrochemical performances, other than the capacity, should be improved to be competitive cathodes, and many strategies have been introduced to enhance DRXs. Fluorination has been shown to inhibit oxygen loss and increase power density. Nevertheless, fluorinated cation-disordered rock-salts still suffer from rapid material deterioration and low scalability which limit their practical applications. This mini-review highlights the key challenges for the commercialization of fluorinated cation-disordered rock-salts, discusses the underlying reasons behind material failure and proposes future development directions.

25 ENERGY STORAGE↗
DOE OSTI2026

Stress and Strain Heterogeneity and Persistence in Uniaxially‐ and Triaxially‐Loaded Sandstone

Two critical questions in brittle rock mechanics are how rocks developed localized strains and to what extent internal stress heterogeneity controls this localization and subsequent macroscopic failure. Definitive answers have not yet been found, but would provide insight into rock fracture mechanics as relevant to hydrocarbon extraction and sequestration. Here, we use synchrotron X‐ray tomography (XRT) and 3D X‐ray diffraction (3DXRD) during uniaxial and triaxial tests on Nugget and Bentheimer sandstones to examine strain and stress localization prior to mechanical failure. 3DXRD was used to measure intra‐granular lattice strains which were used to compute elastic stress tensors of each grain. Digital volume correlation (DVC) was applied to XRT images to determine the strain field in the sample. Both samples featured marked spatial heterogeneity, localization, and temporal persistence of elevated stresses and strains during their mechanical deformation toward failure. Both samples featured a majority of grains with at least one principal stress component that was tensile, a signature of the influence of heterogeneity on stress transmission. Measurements further revealed that compressive stress orientations and statistics evolved in a similar manner to those of inter‐particle forces in loose granular materials, with triaxially‐compressed rock exhibiting enhanced grain stress heterogeneity compared to uniaxially‐compressed rock. Our results complement recent work by others who employed XRT and scanning 3DXRD to study triaxially‐compressed sandstone, but extend those results to uniaxial compression, sandstones of varied porosity, and grain stress measurements throughout the 3D full extent of the samples rather than in a single layer examined with scanning 3DXRD.

58 GEOSCIENCES↗
DOE OSTI2022

Directional Cooling-Induced Fracturing: Directional Cooling Data Under True-Triaxial Stress State

Data includes Directional Cooling-Induced Fracturing (DCIF) testing data using westerly granite blocks. This submission includes data from two samples of westerly granite, lab sample 7 and 8. Files contain stress, temperature and acoustic emission data acquired during polyaxial, laboratory testing of westerly granite blocks for each sample. FILES: .tradb -- files containing acoustic emission waveforms; sqlite3 database .pridb -- files containing basic acoustic emission information (no waveforms); sqlite3 database .geom -- geometry of the AE sensor network (ASCII)

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Seismic Diffusivity and the Influence of Heterogeneity on Injection‐Induced Seismicity

Abstract The spatiotemporal patterns of injection‐induced seismicity (IIS) are commonly interpreted with the concept of a triggering front, which propagates in a diffusion‐like manner with an associated diffusivity parameter. Here, we refer to this diffusivity as the “seismic diffusivity.” Several previous studies implicitly assume that seismic diffusivity is equivalent to the effective hydraulic diffusivity of the subsurface, which describes the behavior of the mean pressure field in heterogeneous porous media. Seismicity‐based approaches for hydraulic characterization or simulations of IIS using domains of homogeneous equivalent porous media are implicitly based on this assumed equivalence. However, seismicity is expected to propagate with the threshold triggering pressure, and thus not be controlled by the evolution of the mean pressure field. We present numerical simulations of fluid injection to compare the seismic and effective hydraulic diffusivities in heterogeneous formations (including fractured rock). The numerical model combines uncoupled, linear pressure diffusion with the Mohr‐Coulomb failure criterion to simulate IIS. We demonstrate that connected pathways of relatively high hydraulic diffusivity in heterogeneous media (particularly in fractured rock domains) allow the threshold triggering pressure to propagate more rapidly than predicted by the effective hydraulic diffusivity. As a result, the seismic diffusivity is greater than the effective hydraulic diffusivity in heterogeneous porous media, possibly by an order of magnitude or more. Additionally, we present a case study of IIS near Soultz‐sous‐Forêts where seismic diffusivity is found to be at least one order of magnitude larger than the effective hydraulic diffusivity.

Haagenson, Ryan↗