Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Geothermal play fairway analysis, part 1: Example from the Snake River Plain, Idaho

The Snake River Plain (SRP) volcanic province overlies the track of the Yellowstone hotspot, a thermal anomaly that extends deep into the mantle. Most of the area is underlain by a basaltic volcanic province that overlies a mid-crustal intrusive complex, which in turn provides the long-term heat flux needed to sustain geothermal systems. Previous studies have identified several known geothermal resource areas within the SRP. For the geothermal study presented herein, our goals were to: (1) adapt the methodology of Play Fairway Analysis (PFA) for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal PFA model for the SRP and identify the most promising plays, using GIS-based software tools that are standard in the petroleum industry. The study focused on identifying three critical resource parameters for exploitable hydrothermal systems in the SRP: heat source, reservoir and recharge permeability, and cap or seal. Data included in the compilation for heat source were heat flow, distribution and ages of volcanic vents, groundwater temperatures, thermal springs and wells, helium isotope anomalies, and reservoir temperatures estimated using geothermometry. Reservoir and recharge permeability was inferred from the analysis of stress orientations and magnitudes, post-Miocene faults, and subsurface structural lineaments based on magnetics and gravity data. Data for cap or seal included the distribution of impermeable lake sediments and clay-seal associated with hydrothermal alteration below the regional aquifer. These data were used to compile Common Risk Segment maps for heat, permeability, and seal, which were combined to create a Composite Common Risk Segment map for all southern Idaho that reflects the risk associated with geothermal resource exploration and identifies favorable resource tracks. Our regional assessment indicated that undiscovered geothermal resources may be located in several areas of the SRP. Two of these areas, the western SRP and Camas Prairie, were selected for more detailed assessment, during which heat, permeability, and seal were evaluated using newly collected field data and smaller grid parameters to refine the location of potential resources. These higher resolution assessments illustrate the flexibility of our approach over a range of scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Photothermal Atomic Force Microscopy Coupled with Infrared Spectroscopy (AFM-IR) Analysis of High Extinction Coefficient Materials: A Case Study with Silica and Silicate Glasses

Photothermal atomic force microscopy coupled with infrared spectroscopy (AFM-IR) brings significant value as a spatially resolved surface analysis technique for disordered oxide materials such as glasses, but additional development and fundamental understanding of governing principles is needed to interpret AFM-IR spectra, since the existing theory described for organic materials does not work for materials with high extinction coefficients for infrared (IR) absorption. This paper describes theoretical calculation of a transient temperature profile inside the IR-absorbing material considering IR refraction at the interface as well as IR adsorption and heat transfer inside the sample. Furthermore, this calculation explains the differences in peak positions and amplitudes of AFM-IR spectra from those of specular reflectance and extinction coefficient spectra. It also addresses the information depth of the AFM-IR characterization of bulk materials. AFM-IR applied to silica and silicate glass surfaces has demonstrated novel capability of characterizing subsurface structural changes and surface heterogeneity due to mechanical stresses from physical contacts, as well as chemical alterations manifested in surface layers through aqueous corrosion.

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

Creation of Laboratory-scale Testbed for Autonomous Monitoring and Control of Subsurface Systems

This report describes the results of an LDRD project aiming to develop a laboratory-scale testbed for autonomous subsurface monitoring and control systems. Modifications were made to an existing polyaxial loading frame to incorporate electrical resistivity tomography (ERT) and acoustic emission (AE) monitoring systems. This allows engineered subsurface processes to be studied under realistic temperatures and stresses in a scaled-down laboratory environment with similar monitoring systems as are deployed in field-scale systems. The paper outlines some challenges that were faced in implementing the ERT system at the laboratory scale and the solutions that were devised. Results are shown from two demonstration tests—one focusing on the ERT system, and one focusing on the AE system.

42 ENGINEERING↗
DOE OSTI2024

Electrical Resistivity Tomography based monitoring of stress perturbations to optimize placement of high-precision strain meters

The Center for Understanding Subsurface Signals and Permeability is a new U.S. Department of Energy Earthshot Center focused on understanding and predicting the long-term evolution of permeability in enhanced geothermal systems. The center will use a highly instrumented testbed within the Sanford Underground Research Facility to conduct field scale experiments that elucidate and test capabilities to simulate geochemical-geomechanical interactions and permeability evolution. Here we demonstrate initial developments using previously collected electrical resistivity tomography (ERT) monitoring data with high-performance multi-physics modelling advancements to inform the optimal location of two new monitoring boreholes. Specifically, ERT monitoring data collected during shear stimulation testing shows marked responses to changes in stress during borehole pressurization. We demonstrate how the same response is being simulated, ultimately to train a machine-learning algorithm to estimate rock properties and enable enhanced prediction of stress and strain responses anticipated during future testing campaigns.

Stress, EGS, CUSSP, 3D Electrical Imaging↗
DOE OSTI2021

Subsurface Characterization of Hydraulic Fracture Test Site-2 (HFTS-2), Delaware Basin

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp Formation of the Permian (Delaware) Basin. This paper focuses on integration, advanced geological characterization, and 3D subsurface modeling of the comprehensive HFTS-2 dataset. The study showcases a multidisciplinary reservoir characterization approach that incorporates geology, petrophysics, geochemistry, geomechanics, microseismic, and subsurface engineering analysis. Subsurface characterization of organic-rich mudstone formations requires understanding complex hydraulic fracture network growth in relation to inherent lithology, geomechanical properties, and interaction with pre-existing natural fractures. This paper presents a characterization workflow incorporating pre- and post-stimulation subsurface data, unique to the HFTS-2 dataset. The study integrated: (1) rock properties from logs, cores, and thin sections; (2) natural and hydraulic fracture descriptions from cores and image logs; (3) local and regional stresses; (4) geomechanics; (5) microseismic; (6) fiber optic (FO) and bottomhole pressure gauge (BHPG) response; and (7) produced fluids analysis. During a stimulation treatment, creation of the stimulated rock volume (SRV) is influenced by several subsurface factors. Key contributing factors include structural context, stress conditions, lithology, facies architecture, pre-existing natural fractures, and geomechanical properties. The HFTS-2 subsurface data integration indicates that the SRV is comprised of a complex juxtaposition of hydraulic fracture swarms, as evidenced by image logs analysis, core description, and microseismic monitoring. The HFTS-2 microseismic event density was used to generate 3D heat maps that serve as a representative SRV footprint, corroborated by secondary datasets. These maps were further integrated with petrophysical and geomechanical characteristics, as well as responses from FO and BHPG, to estimate the lateral and vertical dimensions of the effective fractures. The geological characterization for the HFTS-2 dataset combined with 3D modeling for petrophysical and geomechanical properties provides a strong foundation for subsurface simulation and optimization studies. Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S005R001/2477501/urtec-2021-5243-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5243 The workflow improved our understanding of HFTS-2 hydraulic fracture propagation and characteristics in relation to offset pressure depletion and interaction with pre-existing natural fractures. Analysis showed that fracture geometry varies by stage and by well, and a complex fracture network is generated with varying fracture density. The multidisciplinary workflow presented herein for integration and characterization serves as a foundation to evaluate completion efficiency and estimate areal and vertical stimulation and depletion extent for the project. Furthermore, the workflow and learnings can also be transferred to other unconventional plays.

58 GEOSCIENCES↗
DOE OSTI2021

The subsurface structure of abraded Al–Zn–Mg–Cu alloy

Surface abrasion has been shown to introduce heavy shear stress and produce a few hundred nanometers thick altered surface layer (ASL) on Al alloys. Such ASLs exhibit microstructural changes at room temperature and different corrosion resistance than the underlying substrate. Here, we report how sample dimensions and temperature affect the deformed subsurface microstructure evolution on an abraded Al–Zn–Mg–Cu alloy using transmission electron microscopy (TEM) and atom probe tomography. The ASL microstructure evolves differently in the bulk sample than in the thin TEM specimen. In the freshly abraded ASL, the original η'/η precipitates and grain boundaries were deformed toward the abrasion direction. After 7-day natural aging of the freshly abraded bulk sample, well-defined subgrains formed locally in the ASL, and unusual Al 2 Cu (θ) and Mg x Zn y phases precipitated at the subgrain boundaries. However, during natural aging of this TEM specimen taken from the abraded and aged bulk samples, much faster changes occurred, including the continued formation of subgrains and θ particles in the entire ASL, and the diffusion of Zn atoms out of Mg x Zn y particles to form a Zn particle. On the other hand, no obvious microstructural changes occurred in the ASL during 3-month natural aging of the TEM specimen immediately taken from the freshly abraded sample. Finally, this work advances the understanding of how Al alloy microstructures respond during and following any shear deformation, and complex and unexpected effects of sample dimensions on the behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A new perspective on density and strength loss profiles at the surface of thermally oxidized nuclear graphite

Oxidation of graphite components could influence their designed life in a high-temperature nuclear reactor. The oxidized regions could potentially lower the allowed stress capacity. The American Society of Mechanical Engineers rules for the design and construction of graphite-moderated reactors recommend that subsurface regions that might become excessively damaged by oxidation during reactor operation be identified and excluded from geometry and stress calculations. Identification of oxidation-affected regions is possible, in principle, through complex modeling exercises of reactor behavior during hypothetical accident scenarios coupled with graphite oxidation models, but this procedure may not have the precision needed for informed decisions. Here, this paper proposes an alternate method, based on interpretation of a series of well-designed oxidation experiments, which could augment the designer's tools. The procedure is illustrated by data on oxidation by air of several graphite grades (NBG-18, PCEA, IG-110, R4-650) that are corroborated with independent literature information, when available. The Wichner model for graphite oxidation used for this analysis provides conservative results that could be quickly implemented in the design process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulating Isotropic Rayleigh Wave Horizontal to Vertical Ratios at Rock Valley

Elastic anisotropy in the subsurface can manifest from pre-existing structure or stress-induced cracking – both phenomena are of interest at the Rock Valley site in Nevada. The azimuthal dependence of Rayleigh wave ellipticity, also known as H/V, can reveal elastic anisotropy, but the effect of isotropic variability in H/V must first be accounted for. To that end, we simulate synthetic seismograms from a ring of sources in a fully 3-D isotropic velocity model that captures the geology at Rock Valley. We measure H/V from these waveforms and compare them with H/V predictions obtained using 1D models extracted at specific locations from the 3D model. We find that an average of the azimuthal H/V measurements from the 3-D simulations and 1-D predictions agree, but that the azimuthal variability of H/V is a large fraction of the mean value. These results estimate the noise character in azimuthal H/V so it can be factored out of future anisotropic simulations.

58 GEOSCIENCES↗
DOE OSTI2026

Fracture Intersections under Stress: Laboratory Data and Code [Data set]

The connectivity of natural and induced fractures governs the injection and withdrawal of fluids from subsurface reservoirs. Connectivity depends on intersections that control how fluids mix and move through the entire system. Here, we present data sets from 3D X-ray microscopy measurements of simple fracture networks under stress. 3D printing was used to create prismatic blocks that formed fracture networks composed of 2 orthogonal fractures. The network orientation was either "x" or "+" relative to an applied vertical stress. 3D data sets were collected for normal loads of 25, 100 and 200 Newtons for samples with fracture surfaces with either correlated or uncorrelated asperity distributions. The file contains data from the 12 samples analyzed along with an example code used to extract the intersection geometry. Additional experimental details can be found in the manuscript "Geologic Stress Modulates Fluid Mixing at Fracture Intersections" (10.1038/s43247-026-03525-9)and supplemental information to appear in Communications Earth & Environment in 2026.

02 PETROLEUM↗
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↗
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↗