Search NASA⌕ Search

SEARCH · Search NASA

Results for “Subsurface Stress”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗
DOE OSTI2024

Nitrogen Vacancy Centers in Diamond for Stress Sensing Applications: Theory and Experiments

Sensing of stress under high elevated environmental conditions with high resolution has critical importance for range of applications including earth’s subsurface scanning and geological CO2 storage monitoring. Using first principles density functional theory (DFT) approach combined with the theoretical modelling of low energy Hamiltonian, here we investigate a novel approach to detect unprecedented level of pressure by taking advantage of solid-state electron’s spin of Nitrogen vacancy (NV) centers in nanodiamond. We computationally explore the effect of strain on the defect band edges and band gaps by varying the lattice parameters of diamond supercell hosting a single NV center. A low energy Hamiltonian is developed that includes effect of stress on the energy level of ±1 spin manifold at the ground state. By quantifying the energy level shift and split, we predict the pressure sensing of up to 0.3 MPa/√Hz. We show superiority of the quantum sensing over traditional optical sensing techniques by discussing our results from DFT and theoretical modelling for the frequency shift per unit pressure. The proposed quantum stress sensing could be useful to measure earth’s subsurface vibrations. Our results open avenues for development of sensing technology with a high sensitivity and resolution under extreme pressure conditions that potentially has a wider applicability than existing pressure sensing technologies.

Paudel, Hari P.↗
DOE OSTI2024

Aromatic amino acid metabolism and active transport regulation are implicated in microbial persistence in fractured shale reservoirs

Abstract Hydraulic fracturing has unlocked vast amounts of hydrocarbons trapped within unconventional shale formations. This large-scale engineering approach inadvertently introduces microorganisms into the hydrocarbon reservoir, allowing them to inhabit a new physical space and thrive in the unique biogeochemical resources present in the environment. Advancing our fundamental understanding of microbial growth and physiology in this extreme subsurface environment is critical to improving biofouling control efficacy and maximizing opportunities for beneficial natural resource exploitation. Here, we used metaproteomics and exometabolomics to investigate the biochemical mechanisms underpinning the adaptation of model bacterium Halanaerobium congolense WG10 and mixed microbial consortia enriched from shale-produced fluids to hypersalinity and very low reservoir flow rates (metabolic stress). We also queried the metabolic foundation for biofilm formation in this system, a major impediment to subsurface energy exploration. For the first time, we report that H. congolense WG10 accumulates tyrosine for osmoprotection, an indication of the flexible robustness of stress tolerance that enables its long-term persistence in fractured shale environments. We also identified aromatic amino acid synthesis and cell wall maintenance as critical to biofilm formation. Finally, regulation of transmembrane transport is key to metabolic stress adaptation in shale bacteria under very low well flow rates. These results provide unique insights that enable better management of hydraulically fractured shale systems, for more efficient and sustainable energy extraction.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Corrosive Influence of Carbon Dioxide on Crack Initiation in Quartz: Comparison With Liquid Water and Vacuum Environments

The stimulation of crack growth in quartz and siliceous materials by injecting carbon dioxide (CO 2 ) represents a key technology in long-term carbon storage and in the development of natural gas wells. While this technology is widely used, the molecular impact of CO 2 interactions on the solid matrix is only incompletely understood. In this work, we employ reactive molecular dynamics simulations to study how the CO 2 fluid environment affects the mechanical properties of pre-cracked single-crystal quartz. The thermodynamic conditions of interest are those relevant to subsurface reservoirs. We report how structural properties of quartz—bond length distribution and crack tip shape—evolve upon introduction of a fluid. These properties are directly related to macroscopic quantities of the global stress-strain curves, thus reaffirming the inherent coupling across multiple scales for fluid-solid interactions in the subsurface. Here, we find that CO 2 reduces the fracture toughness of quartz by 12.1% compared to that of quartz in vacuum, thereby promoting crack growth and enhancing fluid transport in the subsurface.

58 GEOSCIENCES↗
DOE OSTI2021

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗
NASA NTRS2006

Crack Initiation and Growth in Rigid Polymeric Closed-Cell Foam Cryogenic Applications

Cryogenic vessels, such as the Space Shuttle External Tank, are often insulated with closed-cell foam because of its low thermal conductivity. The coefficient of thermal expansion mismatch between the foam and metallic substrate places the foam under a biaxial tension gradient through the foam thickness. The total foam thickness affects the slope of the stress gradient and is considered a significant contributor to the initiation of subsurface cracks. Rigid polymeric foams are brittle in nature and any subsurface cracks tend to propagate a finite distance toward the surface. This presentation investigates the relationship between foam thickness and crack initiation and subsequent crack growth, using linear elastic fracture mechanics, in a rigid polymeric closed-cell foam through analysis and comparison with experimental results.

Sayyah, Tarek↗
NASA NTRS1992

Nucleation and growth of rolling contact failure of 440C bearing steel

A 'two-body' elasto-plastic finite element model of 2-dimensional rolling and rolling-plus-sliding was developed to treat the effect of surface irregularities. The model consists of a smooth cylinder in contact with a semi-infinite half-space that is either smooth or fitted with one of 0.4 microns deep or 7 microns deep groove, or a 0.4 microns high ridge-like asperity. The model incorporates elastic-linear-kinematic hardening-plastic (ELKP) and non-linear-kinematic hardening-plastic (NLKP) material constitutive relations appropriate for hardened bearing steel and the 440C grade. The calculated contact pressure distribution is Hertzian for smooth body contact, and it displays intense, stationary, pressure spikes superposed on the Hertzian pressure for contact with the grooved and ridged surface. The results obtained for the 0.4 microns deep groove compare well with those reported by Elsharkawy and Hamrock for an EHD lubricated contact. The effect of translating the counterface on the half space as opposed to indenting the half space with the counter face with no translation is studied. The stress and strain values near the surface are found to be similar for the two cases, whereas they are significantly different in the subsurface. It is seen that when tiny shoulders are introduced at the edge of the groove in the finite element model, the incremental plasticity and residual stresses are significantly higher in the vicinity of the right shoulder (rolling direction is from left to right) than at the left shoulder. This may explain the experimental observation that the spall nucleation occurs at the exit end of the artificially planted indents. Pure rolling calculations are compared with rolling + sliding calculations. For a coefficient of friction, mu = 0.1, the effect of friction is found to be small. Efforts were made to identify the material constitutive relations which best describe the deformation characteristics of the bearing steels in the initial few cycles. Elastic-linear-kinematic hardening-plastic (ELKP) material constitutive relations produce less net plastic deformation in the initial stages for a given stress, than seen in experiments. A new set of constitutive relations: non-linear-kinematic hardening-plastic (NLKP) was used. This material model produces more plasticity than the ELKP model and shows promise for treating the net distortions in the early stages. Techniques for performing experimental measurements that can be compared with the finite element calculations were devised. The measurements are being performed on 9mm-diameter, 440C steel cylindrical rolling elements in contact with 12.5 mm-diameter, 52100 steel balls in a 3-ball-rod fatigue test machine operating at 3600 RPM. Artificial, 7 microns deep, indents were inserted on the running track of the cylindrical rolling elements and profilometer measurements of these indents made, before and after the rolling. These preliminary measurements show that the indents are substantially deformed plastically in the process of rolling. The deformations of the groove calculated with the finite element model are comparable to those measured experimentally.

Gupta, V.↗
DOE OSTI2024

Molecular mechanisms and environmental adaptations of flagellar loss and biofilm growth of Rhodanobacter under environmental stress

Biofilms aid bacterial adhesion to surfaces via direct and indirect mechanisms, and formation of biofilms is considered as an important strategy for adaptation and survival in suboptimal environmental conditions. However, the molecular underpinnings of biofilm formation in subsurface sediment/groundwater ecosystems where microorganisms often experience fluctuations in nutrient input, pH, and nitrate or metal concentrations are underexplored. Here, we examined biofilm formation under different nutrient, pH, metal, and nitrate regimens of 16 Rhodanobacter strains isolated from subsurface groundwater wells spanning diverse levels of pH (3.5 to 5) and nitrates (13.7 to 146 mM). Eight Rhodanobacter strains demonstrated significant biofilm growth under low pH, suggesting adaptations for survival and growth at low pH. Biofilms were intensified under aluminum stress, particularly in strains possessing fewer genetic traits associated with biofilm formation, findings warranting further investigation. Through random barcode transposon-site sequencing (RB-TnSeq), proteomics, use of specific mutants, and transmission electron microscopy analysis, we discovered flagellar loss under aluminum stress, indicating a potential relationship between motility, metal tolerance, and biofilm growth. Comparative genomic analyses revealed the absence of flagella and chemotaxis genes and the presence of a putative type VI secretion system in the highly biofilm-forming strain FW021-MT20. In this study we identified genetic determinants associated with biofilm growth under metal stress in a predominant environmental genus, Rhodanobacter, and identified traits aiding survival and adaptation to contaminated subsurface environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Natural nanoparticle complexes at water-water interfaces

Aqueous two-phase systems (ATPSs) with multifunctional attributes have significant promise as biomimetic materials, but current approaches do not harness their full potential. Here, we show that ATPSs can be stabilized to form microcapsules by the interfacial assembly of chitin nanofibers (ChNF) and rod-like cellulose nanocrystals (CNC). The high structural stability of the ChNF/CNC complexes integrates permeability and transport across the membrane at the liquid/liquid interface. Driven by density and osmotic stress gradients, the microcapsules show switchable motility, including cyclic meniscus-climbing and subsurface transport. These observations demonstrate a self-regulating system with potential for cargo transfer, cell biomimicry, as well as ATPSs-based microreactors and microrobots.

Bioinspired materials↗
DOE OSTI2026

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior (Final Technical Report)

Enhanced Geothermal Systems (EGS) produce clean energy by circulating fluid through hot rock deep underground and bringing that heat to the surface to generate electricity. For this process to work reliably, fluids must be able to move efficiently through networks of natural or engineered fractures in the rock. Enhancing and maintaining subsurface permeability over time is essential for sustainable energy production. However, fluid injection changes the underground temperature, pressure, rock stress, and chemistry, which can alter permeability and sometimes trigger earthquakes. Predicting these interconnected processes remains a key challenge. To address this, we combined high-temperature laboratory experiments with high-fidelity simulations to better understand how fractures in geothermal reservoirs evolve over time. Our experiments measured how fractures respond to stress, slip, slip rate, and chemical reactions under geothermal conditions. These data were integrated into coupled thermal-hydrological-mechanical-chemical and earthquake (THMC+E) models tailored to the Utah FORGE site. The validated modeling framework improves predictions of reservoir performance and seismic response and helps guide operational decisions. This work reduces technical risk and strengthens the scientific foundation needed to make geothermal energy a reliable and scalable clean energy resource.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗
DOE OSTI2023

Skillful Multi‐Month Predictions of Ecosystem Stressors in the Surface and Subsurface Ocean

Abstract Anthropogenic carbon emissions and associated climate change are driving rapid warming, acidification, and deoxygenation in the ocean, which increasingly stress marine ecosystems. On top of long‐term trends, short term variability of marine stressors can have major implications for marine ecosystems and their management. As such, there is a growing need for predictions of marine ecosystem stressors on monthly, seasonal, and multi‐month timescales. Previous studies have demonstrated the ability to make reliable predictions of the surface ocean physical and biogeochemical state months to years in advance, but few studies have investigated forecast skill of multiple stressors simultaneously or assessed the forecast skill below the surface. Here, we use the Community Earth System Model (CESM) Seasonal to Multiyear Large Ensemble (SMYLE) along with novel observation‐based biogeochemical and physical products to quantify the predictive skill of dissolved inorganic carbon (DIC), dissolved oxygen, and temperature in the surface and subsurface ocean. CESM SMYLE demonstrates high physical and biogeochemical predictive skill multiple months in advance in key oceanic regions and frequently outperforms persistence forecasts. We find up to 10 months of skillful forecasts, with particularly high skill in the Northeast Pacific (Gulf of Alaska and California Current Large Marine Ecosystems) for temperature, surface DIC, and subsurface oxygen. Our findings suggest that dynamical marine ecosystem prediction could support actionable advice for decision making.

Environmental Sciences & Ecology↗
NASA NTRS1994

The global distribution of giant radiating dike swarms on Venus: Implications for the global stress state

Magellan radar data of Venus reveal 163 large radial lineament systems composed of graben, fissure, and fracture elements. On the basis of their structure, plan view geometry, and volcanic associations, at least 72% are interpreted to have formed primarily through subsurface dike swarm emplacement, the remainder through uplift or a combination of these two mechanisms. The population of swarms is used to determine regional and global stress orientation. The stress configuration recorded from 330-210 deg E (Aphrodite Terra) is best explained by isostatic compensation of existing long wavelength topography or coupling between mantle flow and the lithosphere. The rest are correlated with concentrations of rifting and volcanism in the Beta-Atla-Themis region. The global stress field on Venus is different than that of Earth, where plate boundary forces dominate.

Grosfils, Eric B.↗
NASA NTRS2012

Magnetic Signature of the Lunar South Pole-Aitken Basin: Character, Origin, and Age

A new magnetic map of the Moon, based on Lunar Prospector (LP) magnetometer observations, sheds light on the origin of the South Pole-Aitken Basin (SPA), the largest and oldest of the recognized lunar basins. A set of WNW-trending linear to arcuate magnetic features, evident in both the radial and scalar observations, covers much of a 1000 km wide region centered on the NW portion of SPA. The source bodies are not at the surface because the magnetic features show no first-order correspondence to any surface topographic or structural feature. Patchy mare basalts of possible late Imbrianage are emplaced within SPA and are inferred to have been emplaced through dikes, directly from mantle sources. We infer that the magnetic features represent dike swarms that served as feeders for these mare basalts, as evident from the location of the Thomson/ Mare Ingenii, Van de Graaff, and Leeuwenhoek mare basalts on the two largest magnetic features in the region. Modeling suggests that the dike zone is between 25 and 50 km wide at the surface, and dike magnetization contrasts are in the range of 0.2 A/m. We theorize that the basaltic dikes were emplaced in the lunar crust when a long-lived dynamo was active. Based on pressure, temperature, and stress conditions prevalent in the lunar crust, dikes are expected to be a dominantly subsurface phenomenon, consistent with the observations reported here.

South Pole-Aitken↗
NASA NTRS1957

Stress-life relation of the rolling-contact fatigue spin rig

The rolling-contact fatigue spin rig was used to test groups of SAE 52100 9.16-inch-diameter balls lubricated with a mineral oil at 600,000-, 675,000-, and 750,000-psi maximum Hertz stress. Cylinders of AISI M-1 vacuum and commercial melts and MV-1 (AISI M-50) were used as race specimens. Stress-life exponents produced agree closely with values accepted in industry. The type of failure obtained in the spin rig was similar to the subsurface fatigue spells found in bearings.

SURFACES, CONTACT - ROLLING↗
DOE OSTI2021

Effect of Shear Displacement and Stress Changes on Fracture Hydraulic Aperture and Flow Anisotropy

Fluid flow in fractures has the potential to drastically change the economic, environmental, and safety risks associated with a subsurface operation. This work focuses on the analysis of experimental shear fracturing and subsequent shear displacement of a shale specimen that contains multiple pre-existing natural fractures. A triaxial direct shear apparatus with concurrent µCT imaging allowed continuous monitoring of sample permeability and changes in fracture geometry during fracture creation, displacement, and changes in stress. Steady-state, low Reynolds number fluid flow simulations were used to determine the changes of each individual fracture’s hydraulic aperture due to displacement and effective stress. The development of anisotropy in the hydraulic aperture of a generated shear fracture was determined for shear displacements of 0.38, 0.74, and 2.04 mm. At high confining stress (30 MPa), shear displacement significantly reduced the hydraulic aperture of the shear fracture to below flow detectable limits and X-ray imaging resolution (23 × 23 × 23 µm 3 ). The shear generated fracture was found to be more sensitive to changes in effective stress compared to the pre-existing natural fractures. Existing models derived from fluid flow equations accounting for stochastic fracture roughness were able to partially fit the simulated flow properties of the pre-existing fractures, but no model was found to fit the flow properties of the newly formed shear fracture. These results are the first report of simultaneous changes in the geometry and flow properties of surrounding, parallel natural fractures as a shear fracture is generated with varying stress and shear displacement conditions.

58 GEOSCIENCES↗
DOE OSTI2021

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗
DOE OSTI2024

Sensing at the Nanoscale Using Nitrogen-Vacancy Centers in Diamond: A Model for a Quantum Pressure Sensor

The sensing of stress under harsh environmental conditions with high resolution has critical importance for a range of applications including earth’s subsurface scanning, geological CO2 storage monitoring, and mineral and resource recovery. Using a first-principles density functional theory (DFT) approach combined with the theoretical modelling of the low-energy Hamiltonian, here, we investigate a novel approach to detect unprecedented levels of pressure by taking advantage of the solid-state electronic spin of nitrogen-vacancy (NV) centers in diamond. We computationally explore the effect of strain on the defect band edges and band gaps by varying the lattice parameters of a diamond supercell hosting a single NV center. A low-energy Hamiltonian is developed that includes the effect of stress on the energy level of a ±1 spin manifold at the ground state. By quantifying the energy level shift and split, we predict pressure sensing of up to 0.3 MPa/Hz using the experimentally measured spin dephasing time. We show the superiority of the quantum sensing approach over traditional optical sensing techniques by discussing our results from DFT and theoretical modelling for the frequency shift per unit pressure. Importantly, we propose a quantum manometer that could be useful to measure earth’s subsurface vibrations as well as for pressure detection and monitoring in high-temperature superconductivity studies and in material sciences. Our results open avenues for the development of a sensing technology with high sensitivity and resolution under extreme pressure limits that potentially has a wider applicability than the existing pressure sensing technologies.

36 MATERIALS SCIENCE↗