Search NASA⌕ Search

SEARCH · Search NASA

Results for “rock damage”

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 OSTI2024

Blue Canyon Dome Emplacement Artifacts

The Blue Canyon Dome (BCD) experiments are a series of intermediate scale subsurface chemical explosions designed to explore the geophysical relationship between scaled yields and formation damage. Phase 2 of this experiment involved explosions designed to evaluate the effect of damage, demonstrating the effect of geologic features and previous damage on the near field seismic signatures, not just in amplitude, but also in spectral content. Based on the emplacement effect on the seismic signature, the moment-based yield estimates showed error between 10s-1000s of percent. This report details an alternative approach to emplacement-based filtering techniques of the seismic signatures to improve yield estimates for explosions in damaged rock. The results will be used to develop an improved inverse yield model that leverages the information in the seismic signature for each source based on emplacement type.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
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 OSTI2023

Thermo-Hydro-Mechanical Modeling of Brine Migration in a Heated Borehole Test in Bedded Salt

Abstract This research paper focuses on the thermo-hydro-mechanical (THM) modeling of brine migration in a heated borehole test conducted as part of the ongoing Brine Availability Test in Salt (BATS) at the Waste Isolation Pilot Plant (WIPP) in New Mexico. It is a component of the international collaboration project DECOVALEX-2023 (DEvelopment of COupled models and their VALidation against EXperiments), which aims to understand the THM processes governing brine flow in heated rock salt repositories through collaborative analysis by multiple research teams. Using the TOUGH–FLAC simulator, THM simulations were performed and compared with data from the BATS phase 1a. This experiment involved two identical horizontal-borehole arrays, one heated and one serving as a control, both equipped with sensor arrays. Analysis of measurements revealed water flow rate surges during heater power transitions, with the highest jump observed during cooling. Acoustic emission activity exhibited distinct patterns in response to heater power changes, suggesting that damage to rock salt is particularly pronounced during the cooling phase. The THM simulations successfully captured these phenomena, highlighting the significance of thermal effects, brine migration, and mechanical behavior in predicting brine availability in heated and damaged rock salt. Our modeling also revealed the critical interplay between heating and cooling-induced damage and its influence on flow properties, particularly affecting brine inflow estimation. Notably, we found that cooling-induced brine inflow spikes result from increased permeability due to tensile dilatancy. These findings have important implications for the development of robust containment strategies and enhance our understanding of the complex processes involved in repository performance.

58 GEOSCIENCES↗
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A micromechanical hydro-mechanical-damage coupled model for layered rocks considering multi-scale structures

An anisotropic Hydro-Mechanical-Damage coupled model was established in the framework of thermodynamics for saturated layered rocks, which contain arbitrarily-distributed cracks of smaller scales and orientated bedding planes of much larger scales. Here, the proposed model considers the anisotropic damage growth and frictional sliding of cracks and shear sliding and dilatancy of bedding planes, as well as the interaction between cracks and bedding planes. Numerical implementation of the nonlinear coupled model was developed based on TOUGHREACT, which is a well-established simulator for multiphase fluid flow and reactive transport analysis. The proposed model was first validated with laboratory water injection test results on rock blocks, and then used to study the effects of injection rates, in situ stresses and bedding planes on fluid injection processes at the field scale, demonstrating the importance to consider the variation of multi-scale structures for better understanding the coupled Hydro-Mechanical responses in layered rocks.

58 GEOSCIENCES↗
DOE OSTI2022

Lithologic controls on microfracturing from legacy underground nuclear explosions

Detection and verification of underground nuclear explosions (UNEs) can be improved with a better understanding of the nature and extent of explosion-induced damage in rock and the effect of this damage on radionuclide migration. Much of the previous work in this area has focused on centimeter- to meter-scale manifestations of damage, but to predict the effect of damage on permeability for radionuclide migration, observations at smaller scales are needed to determine deformation mechanisms. Based on studies of tectonic deformation in tuff, we expected that the heterogeneous tuff layers would manifest explosion-induced damage differently, with welded tuffs showing more fractures and nonwelded tuffs showing more deformation bands. In comparing post-UNE samples with lithologically matched pre-UNE equivalents, we observed damage in multiple lithologies of tuff through quantitative microfracture densities. We find that the texture (e.g., from deposition, welding, alteration, etc.) affects fracture densities, with stronger units fracturing more than weaker units. While we see no evidence of expected deformation bands in the nonwelded tuffs, we do observe, as expected, much larger microfracture densities at close range (<50 m) to the explosive source. We also observe a subtle increase in microfracture densities in post-UNE samples, relative to pre-UNE equivalents, in all lithologies and depths. The fractures that are interpreted to be UNE-induced are primarily transgranular and grain-boundary microfractures, with intragranular microfracture densities being largely similar to those of pre-UNE samples. This work has implications for models of explosion-induced damage and how that damage may affect flow pathways in the subsurface.

58 GEOSCIENCES↗
NASA NTRS1985

Stress assessment and spectral characterization of suspected acid deposition damage in red spruce (Picea Rubens) from Vermont

The effects of acid deposition on Picea rubens are studied. The Picea rubens located at Camels Hump Mt., Mt. Ascutney, and Ripton, VT were analyzed using stress level evaluations, in situ spectral data, pressure bomb analysis, and aircraft sensors. Spruce stress per circular plot and percent spruce mortality are calculated. The relation between stress levels and elevation and exposure and weather patterns is examined. It is observed that variations in the reflectance curves of the foliage and branches are related to cellular health, the type of cellular arrangement, and the degree of leaf tissue hydration; the leaf and twig specimens from high stress sites are more reflective in the red portion of the visible and less reflective in the NIR portion of the spectrum. The pressure bomb data reveal that the xylem water tension is higher in specimens from high stress sites. It is noted that remote sensing permits discrimination and mapping of suspected acid deposition damage.

Rock, B. N.↗
NASA NTRS1986

Use of TMS/TM data for mapping of forest decline damage in the northeastern United States

Remote sensing systems were used to monitor forest decline damage suspected of being due to air pollution. Field activities and aircraft overflights were centered on montane spruce/fir forest sites. Using aircraft data acquired with the Thematic Mapper Simulator (TMS) and LANDSAT Thematic Mapper (TM) during the growing season, extensive areas of forest decline damage were accurately mapped. Seven levels of decline damage are discrininated and mapped and the levels of discriminated damage agree well (rsq-0.94) with visual assessment conducted on the ground. New areas of high damage were discovered. A band ratio (TM5/TM4) is most useful in discriminating and quantifying the various levels of forest decline damage.

Rock, B. N.↗
DOE OSTI2025

Evolving physical and mechanical rock properties during exposure at Earth's surface

Rocks at Earth's surface experience progressive damage due to exposure to environmental conditions. Although environmental damage to rocks is recognized as being time-dependent, our quantitative understanding of its effects on physical and mechanical properties remains limited by field and laboratory studies constrained to short timescales. In this study, we aim to rigorously quantify how these properties evolve under longer natural environmental exposure spanning from the Holocene to the late Pleistocene. We conducted a ‘time series’ analysis of key physical and mechanical rock properties using multiple interior samples from 10 cogenetic granitoid boulders exposed for different amounts of time at two sites in Eastern California, USA. We find that over the late Pleistocene (~30 ka), rocks exhibit porosity increases of 40 ± 4 % and seismic velocity decreases of 7 ± 1 % at both sites. Similarly, over the same time period, rock tensile and compressive strengths declined by an average of 11 ± 4 %, and rock permeability increased on average by 104 ± 50 %. Microscopic analysis of thin sections of the same rocks indicates that microfracture density and intensity also progressively increased over the same timescales. These results suggest that time-dependent fracture propagation due to environmental exposure is likely a central contributor to the observed changes in rock properties. This study provides quantitative documentation that pronounced changes in the physical and mechanical properties of rock progressively occur within thousands of years of exposure.

Critical zone↗
DOE OSTI2020

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗
DOE OSTI2025

Experimental data for damage mechanics simulation challenge

While there are many computational approaches for simulating damage in rock and other materials, few have been ground truth tested with either known experimental data or with blind data sets. Here, in this work, we present a bench-mark laboratory data set for a damage mechanics challenge to compare computational approaches on damage evolution in brittle-ductile materials. The samples were fabricated through additive manufacturing to produce repeatable specimens designed to fail in controlled ways. The failure was induced in the samples using a 3-point bending test to produce different Modes such as Mode I and mixed Modes including I-II, I-III and I-II-III Modes to generate a calibration data set and a blind challenge data set. Data collected included spatial and temporal measurements from traditional digital load–displacement sensors, 2D digital image correlation measurement to map surface deformations, 3D X-ray microscopy to ground-truth the crack-failure geometry, and laser profilometry to capture surface roughness. The data sets are available, on a data repository, to the community to advance computational models to improve our ability to predict damage in brittle-ductile materials.

3-point bending↗
NASA NTRS2009

Cumulative Damage in Strength-Dominated Collisions of Rocky Asteroids: Rubble Piles and Brick Piles

Laboratory impact experiments were performed to investigate the conditions that produce large-scale damage in rock targets. Aluminum cylinders (6.3 mm diameter) impacted basalt cylinders (69 mm diameter) at speeds ranging from 0.7 to 2.0 km/s. Diagnostics included measurements of the largest fragment mass, velocities of the largest remnant and large fragments ejected from the periphery of the target, and X-ray computed tomography imaging to inspect some of the impacted targets for internal damage. Significant damage to the target occurred when the kinetic energy per unit target mass exceeded roughly 1/4 of the energy required for catastrophic shattering (where the target is reduced to one-half its original mass). Scaling laws based on a rate-dependent strength were developed that provide a basis for extrapolating the results to larger strength-dominated collisions. The threshold specific energy for widespread damage was found to scale with event size in the same manner as that for catastrophic shattering. Therefore, the factor of four difference between the two thresholds observed in the lab also applies to larger collisions. The scaling laws showed that for a sequence of collisions that are similar in that they produce the same ratio of largest fragment mass to original target mass, the fragment velocities decrease with increasing event size. As a result, rocky asteroids a couple hundred meters in diameter should retain their large ejecta fragments in a jumbled rubble-pile state. For somewhat larger bodies, the ejection velocities are sufficiently low that large fragments are essentially retained in place, possibly forming ordered "brick-pile" structures.

Housen, Kevin↗
DOE OSTI2022

Precursory Off-Fault Deformation in Restraining and Releasing Step Overs: Insights From Discrete Element Method Models

Accelerating geophysical activity is detected preceding some, but not all, large earthquakes. This observation may indicate that no precursors occur before some earthquakes, or that the instrumentation lacks the required sensitivity. Here, to aid crustal monitoring efforts, we use discrete element method models to identify the locations and styles of deformation that may provide useful information about approaching fault reactivation. We model the reactivation of two healed rough faults in a variety of step over configurations, embedded in a host rock with varying amounts of damage subject to shear velocity loading parallel to the faults. Both the fault geometry and ratio of fault to host rock strength control the amount of off-fault deformation. Consistent with field observations, models with larger steps and more preexisting host rock damage produce higher amounts of off-fault deformation. We assess the size of the continuous regions of high velocities and strains to compare the value of the precursory information of each velocity and strain component. Comparing the three components of the velocity vector suggests that the fault-parallel velocity produces the largest and most temporally continuous regions of elevated velocity. The size of these regions increases toward failure, indicating the usefulness of tracking this component. Comparing the volumetric and shear components of the three-dimensional strain tensor suggests that during most of the interseismic period, the shear strain provides more information about approaching fault slip than the volumetric strain. However, in the days and months preceding fault reactivation, both the shear and volumetric strains provide similarly valuable information.

58 GEOSCIENCES↗
DOE OSTI2019

JGR_Data_Reflection_Project

This research data is based upon the work supported by the U.S. Department of Energy, Office of Basic Energy Sciences, under Award Number DE‐SC0019117. Additional support provided by Colorado School of Mines is highly appreciated.

58 GEOSCIENCES↗
DOE OSTI2023

A constitutive framework for rocks undergoing solid dissolution

Here, we formulate a time-dependent damage theory for rocks subjected to mechanical deformation and solid dissolution. The constitutive description is inspired by the transition state theory, which states that the rate of dissolution is a function of the reactive surface area measured through the crack density in the volume. We use a gradient-enhanced damage framework in which damage depends on the deformation of the material as well as on the amount of solid mass dissolved over time. The gradient-enhanced formulation is characterized by a three-field variational formulation with the solid displacement, nonlocal equivalent strain, and nonlocal rate of solid dissolution as the basic state variables. Traditionally, time-independent damage theories have only allowed damage to increase with increasing external load. In the proposed framework, the degree of damage may increase due to solid dissolution even when the external load is held fixed. In this way, solid dissolution is viewed as a process that is responsible for bringing about rate-dependent effects such as creep and stress-relaxation, which are two common features of geomaterial behavior.

42 ENGINEERING↗