SEARCH · Engineering Papers

Results for “subsurface diffusion”

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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Jupiter System Data Analysis Program: Mechanisms, Manifestation, and Implications of Cryomagmatism on Europa

The objectives of the work completed under NASA Grant NAG5-8898 were (i) to document and characterize the low-albedo diffuse surfaces associated with triple bands and lenticulae, (ii) to determine their mechanisms of formation, and (iii) to assess the implications of these features for the resurfacing (in space and time) of the Europa and the nature of the Europan interior. The approach involved a combination of processing and analysis of Solid State Imaging data returned by the Galileo spacecraft during the primary and extended mission phases, together with numerical modeling of the physical processes interpreted to the observed features. We have modeled the formation of Europan triple explosive venting of cryoclastic material from bands and lenticulae halos by two processes: (i) a liquid layer in the Europan interior, and (ii) lag deposit formation by the thermal influence of subsurface cryomagmatic intrusions. We favor the latter hypothesis for explaining these features, and further suggest that a liquid water or brine intrusion is required to provide sufficient lateral heating of surface ice to explain the 25 km size of the largest features. (Solid ice diapirs, even under the most favorable conditions, become thermally exhausted before they heat significant lateral distances). We argue that water circulating in open fractures, or repeated cryomagmatic 'diking' events would provide sufficient thermal input to produce the observed features. Thus our work argues for the existence of a liquid beneath Europa's surface. Our results might most easily be explained by the presence of a continuous liquid layer (the putative Europan ocean); this would concur with the findings of the Galileo magnetometer team. However, we cannot rule out the possibility that discrete liquid pockets provide injections of fluid closer to the surface.

Fagents, Sarah A.↗

DOE OSTI2023

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗

DOE OSTI2022

Two-Phase Flow Mechanisms Controlling CO2 Intrusion into Shaly Caprock

Abstract Geologic carbon storage in deep saline aquifers has emerged as a promising technique to mitigate climate change. CO 2 is buoyant at the storage conditions and tends to float over the resident brine jeopardizing long-term containment goals. Therefore, the caprock sealing capacity is of great importance and requires detailed assessment. We perform supercritical CO 2 injection experiments on shaly caprock samples (intact caprock and fault zone) under representative subsurface conditions. We numerically simulate the experiments, satisfactorily reproducing the observed evolution trends. Simulation results highlight the dynamics of CO 2 flow through the specimens with implications to CO 2 leakage risk assessment in field practices. The large injection-induced overpressure drives CO 2 in free phase into the caprock specimens. However, the relative permeability increase following the drainage path is insufficient to provoke an effective advancement of the free-phase CO 2 . As a result, the bulk CO 2 front becomes almost immobile. This implies that the caprock sealing capacity is unlikely to be compromised by a rapid capillary breakthrough and the injected CO 2 does not penetrate deep into the caprock. In the long term, the intrinsically slow molecular diffusion appears to dominate the migration of CO 2 dissolved into brine. Nonetheless, the inherently tortuous nature of shaly caprock further holds back the diffusive flow, favoring safe underground storage of CO 2 over geological time scales.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

From the Interstellar Medium to Ocean Worlds: new challenges for Laboratory Astrophysics

From the Interstellar Medium to Ocean Worlds: new challenges for Laboratory Astrophysics For the past 35+ years, laboratory experiments in Astrophysics have been focusing on identifying and understanding the molecular species and chemistry occurring in interstellar environments, such as the diffuse interstellar medium (ISM) and dense molecular clouds. A significant portion of these laboratory efforts concerns investigations of carbon-based molecules ranging in size from a single carbon atom (e.g., CO, CO2, H2CO) to large polycyclic aromatic hydrocarbon (PAH) molecules potentially consisting of more than one hundred carbon atoms as well as their reactions with H2O and other interstellar species. The laboratory facilities built to conduct these investigations are capable of monitoring the chemistry and measuring the spectra under a wide range of interstellar conditions (e.g., low temperatures and pressures). The past 15 years of space exploration revealed the presence of at least nine other Ocean Worlds, besides Earth, in our Solar System. These worlds include moons and dwarf planets, ice-covered surfaces, subsurface oceans, as well as geysers venting H2O and other compounds, including organics in at least one environment, into space. Although these worlds represent environments that are significantly different from those traditionally considered the realm of Laboratory Astrophysics, their facilities are well situated to provide insight into the chemistry occurring in and around these worlds as well as experimental data for the interpretation of observations from missions visiting these Ocean World. This presentation will discuss the areas where Laboratory Astrophysics can provide experimental data to help understand the chemistry, and mission data of Ocean Worlds, and the challenges research into Ocean Worlds presents to our current facilities and how they can be overcome. Lastly, this presentation will introduce the new ICEE (In-situ Carbon Evolution Experiments) facility at NASA Ames Research Center as an example of how a lab dedicated to Laboratory Astrophysics can adapt itself to also meet the needs for Ocean Worlds research.

Andrew Lige Mattioda↗

DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗

DOE OSTI2021

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗

NASA NTRS1998

Sources of Hydrogen as Food for Deep Microbial Communities

To survive in deep subsurface environments autolithotrophic microbial communities require a sustainable food supply. One possible source is H2 which forms when H2O reacts with ferrous iron at rock surfaces or mineral grain boundaries to produce H2 plus ferric iron. The amount of H2 that can be supplied in this way, however, is relatively small and may not last for more than a few hundred or thousand years. A much larger reservoir of H2 exists in the rocks, inside mineral grains, arising from an as yet little-known redox conversion that affects OH- in nominally anhydrous minerals. These OH- represent small amounts of "water" that become incorporated during crystallization in H2O-laden environments. A corollary of the H2 formation from OH- is the formation of peroxy, an oxidized form of oxygen. While the peroxy become part of the mineral structure, the H2 molecules are diffusively mobile and can escape from within the mineral grains, entering the intergranular space. Assuming cautiously realistic number densities of OH- undergoing the in situ redox conversion to H2 plus peroxy, a 10 km deep rock column is expected to contain enough H2 to allow for a constant degassing rate of 50-100 nmole H2 per day per sq cm over 30 million years.

Freund, Friedemann↗

DOE OSTI2020

Nanoparticle Tracking to Probe Transport in Porous Media

From the granular and fractured subsurface environment to highly engineered polymer membranes used in pharmaceutical purification, porous materials are ubiquitous in nature and industrial applications. In particular, porous media are used extensively in processes including water treatment, pharmaceutical sterilization, food/beverage processing and heterogeneous catalysis, where hindered mass transport is either essential to the process or a necessary, but undesirable limitation. Unfortunately, there are currently no universal models capable of predicting mass transport based on a description of the porous material because real porous materials are complex, and because many coupled dynamic mechanisms (e.g. adsorption, steric effects, hydrodynamic effects, electrostatic interactions, etc.) give rise to the observed macroscopic transport phenomena. While classical techniques, like nuclear magnetic resonance and dynamic light scattering, provide useful information about mass transport in porous media at the ensemble level, they provide limited insight into the microscopic mechanisms that give rise to complex phenomena such as anomalous diffusion, hindered pore-space accessibility, and unexpected retention under flow, among many others. To address this issue, we have developed refractive index matching imaging systems, combined with single-particle tracking methods, allowing the direct visualization of single particle motion within a variety of porous materials. Here, we summarize our recent efforts to advance the understanding of nanoparticle transport in porous media using single-particle tracking methods in both fundamental and applied scenarios. First, we describe the basic principles for two-dimensional and three-dimensional single-particle tracking in porous materials. Then, we describe concrete examples of nanoparticle transport in porous materials from two perspectives: (1) understanding fundamental elementary particle transport processes in porous media, including pore accessibility and cavity escape, which limit transport in porous media; (2) facilitating applications in industrial processes, e.g. by understanding the mechanisms of particle fouling and remobilization in filtration membranes. Finally, we provide an outlook of opportunities associated with investigating other types of mass transport in confined environments using single-particle tracking methods, including electrophoretic and self-propelled motion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2002

Hydrogen in rocks: an energy source for deep microbial communities

To survive in deep subsurface environments, lithotrophic microbial communities require a sustainable energy source such as hydrogen. Though H2 can be produced when water reacts with fresh mineral surfaces and oxidizes ferrous iron, this reaction is unreliable since it depends upon the exposure of fresh rock surfaces via the episodic opening of cracks and fissures. A more reliable and potentially more voluminous H2 source exists in nominally anhydrous minerals of igneous and metamorphic rocks. Our experimental results indicate that H2 molecules can be derived from small amounts of H2O dissolved in minerals in the form of hydroxyl, OH- or O3Si-OH, whenever such minerals crystallized in an H2O-laden environment. Two types of experiments were conducted. Single crystal fracture experiments indicated that hydroxyl pairs undergo an in situ redox conversion to H2 molecules plus peroxy links, O3Si/OO\SiO3. While the peroxy links become part of the mineral structure, the H2 molecules diffused out of the freshly fractured mineral surfaces. If such a mechanism occurred in natural settings, the entire rock column would become a volume source of H2. Crushing experiments to facilitate the outdiffusion of H2 were conducted with common crustal igneous rocks such as granite, andesite, and labradorite. At least 70 nmol of H2/g diffused out of coarsely crushed andesite, equivalent at standard pressure and temperature to 5,000 cm3 of H2/m3 of rock. In the water-saturated, biologically relevant upper portion of the rock column, the diffusion of H2 out of the minerals will be buffered by H2 saturation of the intergranular water film.

Hydrogen/analysis↗

DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗

DOE OSTI2022

Chemical and isotopic evolution of flowback fluids from the Utica Gas Shale Play, Eastern Ohio USA

For this work, hydraulic fracturing flowback fluids were collected from two Utica/Point Pleasant well pads in eastern Ohio. One site was in the wet gas zone (UPPW), while the site ~50 km south consisted of four wells on the same pad in the dry gas zone (UPPS). Samples of input fluids also were collected before and during hydraulic fracturing. Flowback fluids are Na-Ca-Cl brines with total dissolved salt (TDS) concentrations that increase over several months from ~100 to 200 g/L. The slightly higher TDS of the dry gas fluids are in part due to recycled flowback used as input fluids, and in part due to lower volume of water used in hydraulic fracturing. Concentrations of most major ions (Ca 2+ , Mg 2+ , Na + , Sr2 + , Fe, Mn, Cl – , Br – ) are similar for the five wells sampled, although small but systematic changes occur in the major element ratios over time. Most notably, an increase in the Sr/Cl ratio corresponds to an increase in the 87 Sr/ 86 Sr ratio in the fluids, suggesting interactions with a more radiogenic Sr source in the subsurface. Dissolved Ba concentrations and Ra activities were different between the two sites, reflecting a high SO 4 2– fluid used at the UPPW site, and water-rock reactions occurring during hydraulic fracturing at the UPPS site. Water oxygen (δ 18 O) and hydrogen (δD) isotopes for input fluids used in the UPPW4 well and fresh water used for the UPPS wells fall on the Global Meteoric waterline (GMWL). Flowback fluids from both sites are relatively enriched in δ 18 O and δD compared to the input, but do not appear to follow a simple mixing trend, suggesting reaction and isotopic exchange with carbonates and fractionation due to imbibition in the rock. Cl isotopes, δ 37 Cl in the FP fluids varied from ~ –0.43 to +0.13‰, the largest variation was observed in the earlier stages of flowback, while in the later stages δ 37 Cl exhibited a small but systematic increase over time, suggesting diffusion control of isotopic composition. Some of the trace species measured (dissolved Fe, Mn, thiosulfate and organic acids) do not follow the same trends as the major ions, suggesting contributions from input fluids, or microbially mediated reactions are exerting control on their concentrations. Input water chemistry exerts an important control on the Sr and Ba concentrations in flowback water. High SO 4 2– in the input fluids used for hydraulically fracturing the wet gas well leads to precipitation of barite-celestite in flowback fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

The Role of Surface Hydrophobicity on the Structure and Dynamics of CO2 and CH4 Confined in Silica Nanopores

Advancing a portfolio of technologies that range from the storage of excess renewable natural gas for distributed use to the capture and storage of CO 2 in geological formation are essential for meeting our energy needs while responding to challenges associated with climate change. Delineating the surface interactions and the organization of these gases in nanoporous environments is one of the less explored approaches to ground advances in novel materials for gas storage or predict the fate of stored gases in subsurface environments. To this end, the molecular scale interactions underlying the organization and transport behavior of CO 2 and CH 4 molecules in silica nanopores need to be investigated. To probe the influence of hydrophobic surfaces, a series of classical molecular dynamics (MD) simulations are performed to investigate the structure and dynamics of CO 2 and CH 4 confined in OH-terminated and CH 3 -terminated silica pores with diameters of 2, 4, 6, 8, and 10 nm at 298 K and 10 MPa. Higher adsorption extents of CO 2 compared to CH 4 are noted on OH-terminated and CH 3 -terminated pores. The adsorbed extents increase with the pore diameter. Further, the interfacial CO 2 and CH 4 molecules reside closer to the surface of OH-terminated pores compared to CH 3 -terminated pores. The lower adsorption extents of CH 4 on OH-terminated and CH 3 -terminated pores result in higher diffusion coefficients compared to CO 2 molecules. The diffusivities of both gases in OH-terminated and CH 3 -terminated pores increase systematically with the pore diameter. The higher adsorption extents of CO 2 on OH-terminated and CH 3 -terminated pores are driven by higher van der Waals and electrostatic interactions with the pore surfaces, while CH 4 adsorption is mainly due to van der Waals interactions with the pore walls. These findings provide the interfacial chemical basis underlying the organization and transport behavior of pressurized CO 2 and CH 4 gases in confinement.

Mohammed, Sohaib↗

DOE OSTI2020

Interpretation of Hydraulic Fractures Based on Microseismic Response in the Marcellus Shale, Monongalia County, West Virginia, USA: Implications for Shale Gas Production

Summary Hydraulic fracturing is critical for extracting shale gas in the subsurface. The treatment technique of multistage hydraulic fracturing is widely used to maximize production. In multistage hydraulic fracturing, not all pumping stages make the same contribution for production, although the designed stimulation process is almost same in every pumping stage. In this study, we characterize the microseismic responses of hydraulic fracturing. Significant variations of microseismic characteristics are observed among the different pumping stages. Pre-existing natural fractures are examined along the horizontal well in selected stages. Combined with the Mohr’s circle analysis, results of the fracture study show that induced hydraulic fractures can be captured by pre-existing natural fractures. Induced hydraulic fractures are simulated by the unconventional fracture model (UFM), and the result reveals that the stress-shadow effect diverts the direction of hydraulic fractures. The diverted hydraulic fracture affects the development of the hydraulic fracture network, which has influence on production. The treatment of two-step pumping is investigated by application of hydraulic diffusivity. Hydraulic fracturing performance and production could be optimized by the treatment of two-step pumping in a single stage. The second pumping creates new fractures and fills the fractures with additional proppants to maintain production for a long duration.

Energy & Fuels↗

DOE OSTI2026

Effects of temperature and dose rate on ion-irradiated γ-LiAlO 2 pellets

Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by competing processes of defect production, driven by the dose rate, and defect recovery, controlled by diffusion, interaction, and annihilation. Here, this study investigates the effects of irradiation temperature and the dose rate on microstructural evolution, deuterium retention, and lithium volatilization in γ-LiAlO 2 pellets subjected to sequential He + and D + ion irradiation. Experiments were performed to a total fluence of 3 × 10 17 (He + + D + )/cm 2 at 623, 673, 723, and 773 K with an average He + dose rate of 7.7 × 10 −4 dpa/s, and to 2 × 10 17 (He + + D + )/cm 2 at 773 K with dose rates of 6.8 × 10 −5 , 2.9 × 10 −4 , and 7.3 × 10 −4 dpa/s. At 623 K, the microstructure was dominated by cavities and fractures with no observable precipitate formation, while small precipitates emerged at 673 K. Increasing the irradiation temperature to 723–773 K promoted the formation of larger, faceted LiAl 5 O 8 precipitates, and surface amorphization, accompanied by pronounced lithium depletion and H–D isotopic exchange. At 773 K, medium and high dose rates produced an amorphized surface layer over a crystalline subsurface containing LiAl 5 O 8 precipitates and blisters at the crystalline–amorphous interface, whereas low-dose-rate irradiation preserved surface crystallinity with cavities distributed in the matrix, around precipitates, and along grain boundaries. Precipitate morphology was anisotropic with limited size dependence on the dose rate. These results elucidate the coupled effects of temperature and the dose rate and demonstrate that sequential He + and D 2 + irradiation at 773 K reproduces key microstructural features and H isotope behavior observed in neutron-irradiated γ-LiAlO 2 at 573 K.

dose rate effects↗

DOE OSTI2021

Competitive formation of intercalated versus supported metal nanoclusters during deposition on layered materials with surface point defects

Intercalated metal nanoclusters (NCs) can be formed under the surface of graphite after sputtering to generate surface “portal” defects that allow deposited atoms to reach the subsurface gallery. Although, there is a competition between formation of supported NCs on top of the surface and intercalated NCs under the surface, the latter only dominating at sufficiently high temperature. A stochastic model incorporating appropriate system thermodynamics and kinetics is developed to capture this complex and competitive nucleation and growth process. Kinetic Monte Carlo simulation shows that the model captures experimental trends observed for Cu and other metals and reveals that higher temperatures are needed to facilitate detachment of atoms from supported NCs enabling them to reach the gallery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Targeted Dy intercalation under graphene/SiC for tuning its electronic band structure

Metal intercalation of graphene is a promising method to tune its electronic band structure and generate novel electronic and topological phases. The tuning depends critically on the ability to bond the intercalated atoms at predesigned, subsurface interlayer locations because the emerging band structure depends on metal location. In this work, we have studied Dy intercalation under single-layer graphene (SLG) on SiC using spot profile analysis–low-energy electron diffraction and scanning tunneling microscopy (STM). The experimental work is complemented with density-functional theory (DFT) analysis. Because different diffraction spots originate from different subsurface interlayer regions, it is possible to identify changes in the intercalation location by monitoring the spot intensity as a function of growth conditions. DFT calculations of the chemical potential as a function of intercalated Dy coverage support the variation of the stability of the intercalated phase at different intercalated locations. The preferred location is confirmed from STM studies showing the removal of the 6 × 6 moiré corrugation at the preferred location, observed at higher Dy coverage.

2-dimensional systems↗

DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗

NASA NTRS2009

Mid-Infrared Reflectance Imaging of Thermal-Barrier Coatings

An apparatus for mid-infrared reflectance imaging has been developed as means of inspecting for subsurface damage in thermal-barrier coatings (TBCs). The apparatus is designed, more specifically, for imaging the progression of buried delamination cracks in plasma-sprayed yttria-stabilized zirconia coatings on turbine-engine components. Progression of TBC delamination occurs by the formation of buried cracks that grow and then link together to produce eventual TBC spallation. The mid-infrared reflectance imaging system described here makes it possible to see delamination progression that is invisible to the unaided eye, and therefore give sufficiently advanced warning before delamination progression adversely affects engine performance and safety. The apparatus (see figure) includes a commercial mid-infrared camera that contains a liquid-nitrogen-cooled focal plane indium antimonide photodetector array, and imaging is restricted by a narrow bandpass centered at wavelength of 4 microns. This narrow wavelength range centered at 4 microns was chosen because (1) it enables avoidance of interfering absorptions by atmospheric OH and CO2 at 3 and 4.25 microns, respectively; and (2) the coating material exhibits maximum transparency in this wavelength range. Delamination contrast is produced in the midinfrared reflectance images because the introduction of cracks into the TBC creates an internal TBC/air-gap interface with a high diffuse reflectivity of 0.81, resulting in substantially higher reflectance of mid-infrared radiation in regions that contain buried delamination cracks. The camera is positioned a short distance (.12 cm) from the specimen. The mid-infrared illumination is generated by a 50-watt silicon carbide source positioned to the side of the mid-infrared camera, and the illumination is collimated and reflected onto the specimen by a 6.35-cm-diameter off-axis paraboloidal mirror. Because the collected images are of a steady-state reflected intensity (in contrast to the transient thermal response observed in infrared thermography), collection times can be lengthened to whatever extent needed to achieve desired signal-to-noise ratios. Each image is digitized, and the resulting data are processed in several steps to obtain a true mid-infrared reflectance image. The raw image includes thermal radiation emitted by the specimen in addition to the desired reflected radiation. The thermal-radiation contribution is eliminated by subtracting the image obtained with the illumination off from the image obtained with the illumination on. A flat-field correction is then made to remove the effects of non-uniformities in the illumination level and pixel-to-pixel variations in sensitivity. This correction is performed by normalizing to an image of a standard object that has a known reflectance at a wavelength of 4 microns. After correction, each pixel value is proportional to the reflectance (at a wavelength of 4-microns) at the corresponding location on the specimen. Mid-infrared reflectance imaging of specimens that were thermally cyc led for different numbers of cycles was performed and demonstrated t hat mid-infrared reflectance imaging was able to monitor the gradual delamination progression that occurs with continued thermal cycling . Reproducible values were obtained for the reflectance associated w ith an attached and fully delaminated TBC, so that intermediate refle ctance values could be interpreted to successfully predict the numbe r of thermal cycles to failure.

Edlridge, Jeffrey I.↗