Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pore pressure”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Microsized Pore Structure Determination in EPDM Rubbers Using High-Pressure 129 Xe NMR Techniques

Microsized pore parameters, such as pore size and distance between pores in a series of model EPDM rubbers, were determined in situ under the pressure of 500 psi using 129 Xe nuclear magnetic resonance (NMR) techniques: spin–lattice (T 1 ) and spin–spin (T 2 ) relaxation measurements, pulsed-field gradient (PFG) NMR, and two-dimensional exchange spectroscopy (2D EXSY). The T 1 /T 2 (>>1) ratio for the xenon confined in the pores is larger than that for nonconfined free xenon. This suggests that almost the entire pore surface interacts with xenon atoms like a closed pore. While these pores still connect each other through very narrow diffusion/exchange channels, it is possible to observe the echo decay in PFG-NMR and cross-peaks in 2D EXSY. The results show that both diffusion (D pore ≈ 2.1 × 10 –10 m 2 /s) and exchange (exchange rate, τ exch = a few tens of milliseconds) of xenon between a pore within the material and outer surface are prolonged. The exchange distances (l), which correspond to the xenon gas penetration depth, were estimated to be 70–100 μm based on the measured diffusion coefficients and exchange rate (1/τ exch ). NMR diffraction analysis reveals that pore size (a) and pore distance (b) are on the order of magnitude of micrometers and tens of micrometers, while the diffusion coefficients of xenon gas in the diffusion channels (D eff ) are about 10 –8 m 2 /s. Overall, this study suggests that the pores with a few micrometers connected through very narrow flowing channels with the length of several tens of micrometers are developed 70 to 100 μm below the rubber surface. Furthermore, the overall steady-state diffusion of xenon is slower, approximately 2 orders of magnitudes, than the diffusion in the channel between the pores. Finally, the pore and exchange distances correlated with the composition of rubbers showed that the properties of EPDM rubber as a high-pressure gas barrier could be improved by reducing the size of cracks and the depth of gas penetration by the addition of both carbon black and silica fillers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A Dynamic Pore Network Model for Imbibition Simulation Considering Corner Film Flow

Wetting films can develop in the corners of angular pores under strong wetting conditions. Modeling the dynamics of corner film remains elusive using direct numerical simulations because of the significant scale difference between main meniscus and corner film flow. In this paper, the modified interacting capillary bundle model (ICB), developed in our previous work to describe accurately corner film dynamics in a single square tube, is incorporated into a single-pressure dynamic pore network model (DPNM) to simulate imbibition in strongly wetting porous media with corner film flow. The traditional pore network is decomposed into several layers of interacting subpore networks where the 0th layer of subpore network simulates the main meniscus flow and higher layers the corner film flow. The fluid flow between different layers is captured by interlayer throats. In addition, the snap-off mechanism caused by the thickening of wetting corner film is considered. The accuracy of the developed model is validated for four cases: spontaneous imbibition in a single square tube, wetting fluid redistribution through corner films under a capillary pressure difference, snap off in a narrow throat connecting two large pores, and imbibition dynamics in a real microfluidic porous geometry. The validated model is then used to simulate both spontaneous and controlled imbibition in a pore network with random pore size distribution. Finally, the interaction between corner film and main meniscus flow in porous media is analyzed from a pore-scale perspective.

58 GEOSCIENCES↗
DOE OSTI2022

Frictional Behavior Downdip Along the Subduction Megathrust: Insights From Laboratory Experiments on Exhumed Samples at In Situ Conditions

Subduction megathrusts exhibit a spectrum of slip modes, including catastrophic earthquakes. Although the mechanical and frictional properties of materials sampled from subduction zones have been studied extensively, few datasets have been collected for compositions and at pressure and temperature conditions representative of those in situ. The Nankai subduction zone in southwest Japan is a well-studied margin, and abundant data provide an opportunity to advance our understanding of fault and earthquake physics. Here, we use samples exhumed in the Shimanto and Sanbagawa Belts on Shikoku Island of southwest Japan that represent analogs for materials along the present-day megathrust at depths of ~5–>25 km, and we shear these at their peak in situ pressure-temperature (P-T) conditions. Effective normal stresses range from 28 to 192 MPa, and temperatures from 105°C to 470°C. We used pore fluid pressures of 45–240 MPa, corresponding to fluid overpressure ratios λ of 0.65 and 0.90. Slip velocities of 0.1–100 μm/s were used, in order to focus on the nucleation of instability and earthquakes. We found predominantly velocity-strengthening (inherently stable) behavior under all conditions for λ = 0.65. For λ = 0.90, velocity-weakening behavior was observed at 350°C, with velocity-strengthening behavior at lower and higher temperatures. The rate/state frictional stability parameter (a-b) increases with slip velocity at temperatures up to ~200°C and remains constant or decreases with slip velocity at higher temperatures. Altogether, our results demonstrate the potentially important roles of both temperature and slip velocity in controlling the distribution of stress and frictional rheology along subduction thrusts.

58 GEOSCIENCES↗
DOE OSTI2021

Sub-nanometer confinement enables facile condensation of gas electrolyte for low-temperature batteries

Abstract Confining molecules in the nanoscale environment can lead to dramatic changes of their physical and chemical properties, which opens possibilities for new applications. There is a growing interest in liquefied gas electrolytes for electrochemical devices operating at low temperatures due to their low melting point. However, their high vapor pressure still poses potential safety concerns for practical usages. Herein, we report facile capillary condensation of gas electrolyte by strong confinement in sub-nanometer pores of metal-organic framework (MOF). By designing MOF-polymer membranes (MPMs) that present dense and continuous micropore (~0.8 nm) networks, we show significant uptake of hydrofluorocarbon molecules in MOF pores at pressure lower than the bulk counterpart. This unique property enables lithium/fluorinated graphite batteries with MPM-based electrolytes to deliver a significantly higher capacity than those with commercial separator membranes (~500 mAh g −1 vs. <0.03 mAh g −1 ) at −40 °C under reduced pressure of the electrolyte.

25 ENERGY STORAGE↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The evolution of seafloor venting from hydrate-sealed gas reservoirs

Here we use 3D seismic data to show that three rows of seafloor gas mounds can be traced downward to leak points that lie at the hydrate-gas contact within three individual dipping coarse-grained sand bodies in the Terrebonne Basin, Walker Ridge block 313, northern Gulf of Mexico. We predict the overpressure within the sand bodies by assuming that the gas pressure at the vent points equal the least principal stress. We interpret that free gas accumulates at the base of the hydrate stability zone, trapped by the overlying sand which has a high methane hydrate saturation. The free gas accumulates until the gas pressure at the base of the hydrate stability zone reaches the least principal stress in the overlying mudrock, whereupon hydraulic fractures form and fluids are vented to the surface. The warm rising fluids and perhaps localized exothermic formation of hydrate raise the local salinity and temperature. This process progressively shifts the base of the hydrate stability zone to shallower depths and dissociates the hydrate seal within the sand, which creates new leak points and results in the observed migration of the seafloor vents. Within the southwest corner of the Terrebonne Basin, this process has repeated multiple times within the Blue, Orange, and Green sands. Here this study shows how the hydrate stability zone can influence the location of fluid expulsion and in turn be affected by the warm, rising, saline fluids sealed by hydrate.

58 GEOSCIENCES↗
DOE OSTI2024

Enforcing global constraints for the dispersion closure problem: τ 2 -SIMPLE algorithm

Permeability and effective dispersion tensors are critical parameters to characterize flow and transport in porous media at the continuum scale. Homogenization theory defines a framework in which such effective properties are first computed from solving a closure problem in a repeating unit cell of the periodic microstructure and then used in a macroscopic formulation for efficient computation. The closure problem is formulated as a local boundary value problem subjected to global constraints, which guarantee the uniqueness of the solution and can be difficult to satisfy for complex geometries and at high flow conditions. These constraints also ensure that pore-scale pressure, velocity, and concentration fields can be accurately reconstructed from the closure variable. Building on a previous work, here we present a framework that allows to satisfy global constraints associated to both the permeability and the dispersion closure problems by introducing two artificial time scales. The algorithm, called τ 2 -SIMPLE, computes both permeability and effective dispersion given an arbitrarily complex geometry and flow condition. Furthermore, this algorithm is demonstrated to be accurate for both 2D and 3D geometries across varying flow conditions, and thus it can be used to quickly characterize effective properties from porous media images in many applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Innovations in underground hydrogen storage with multiphysics simulations, optimization, and monitoring: A review

Underground Hydrogen Storage (UHS) is a promising solution for large-scale energy storage and a critical component in advancing low-carbon energy system. Ensuring the safety and efficiency of UHS necessitates a comprehensive understanding of multiphysical interactions driven by cyclic pore fluid pressure fluctuations and coupled physicochemical processes. Here, this review examines the key geomechanical responses in UHS, including rock property variations under cyclic loading, fracture evolution and propagation, reservoir stress sensitivity, and fault stability. It also explores the impact of geochemical and microbial reactions on geomechanical characteristics. We provide an in-depth analysis of Thermal-Hydraulic-Mechanical-Chemical (THMC) coupled numerical simulations, highlighting their potential for future multi-scale modeling. Limitations of current machine learning (ML) approaches in addressing UHS challenges are highlighted, emphasizing the need for innovative ML-based methodologies. Operational strategies for hydrogen injection and production are reviewed, focusing on safety, efficiency, and economic viability. The necessity for multi-objective optimization (MOO) to balance storage efficiency, risk mitigation, and cost-effectiveness is also discussed. Current monitoring technologies are evaluated to ensure safe and efficient UHS operations. Finally, this review identifies critical knowledge gaps and underscores the importance of advancing geomechanical understanding under multiphysics-coupling. We highlight the need for ML-driven multiphysics theories, enhanced modeling techniques, and robust optimization strategies to improve UHS performance. This study serves as a comprehensive reference for future research and the large-scale implementation of UHS systems.

25 ENERGY STORAGE↗
DOE OSTI2021

Seasonal Seismicity in the Lake Biwa Region of Central Japan Moderately Modulated by Lake Water Storage Changes

Seasonal seismicity in Japan has been observed and attributed to variations induced by hydrological loading, as well as atmospheric and pore-fluid pressure. We examine the seasonal variations of crustal (<15 km depth) and low magnitude (M < 4.5) earthquakes in central Japan near the Biwako-Seigan Fault Zone (BSFZ), and analyze their correlation with the annual hydrological cycle of nearby Lake Biwa, the largest freshwater lake in Japan. The 2002–2018 lake water storage is estimated using water-level gauge and satellite altimetry data. Regional surface mass loading from soil moisture and snow predicted by the Global Land Data Assimilation System is used to augment the lake mass variation. We compute the stress and stressing-rate changes induced by the surface mass loading from lake storage, snow, and soil moisture at seismogenic depth (10 km) along the BSFZ using a finite element model incorporating a 3D Earth structure. The results show that seasonal surface mass loading generates 1–5 kPa peak-to-peak Coulomb stress change and 2.5–15 kPa/year stressing-rate variations on the fault segments along the BSFZ. Two of the four fault segments present significant correlation between local seismicity rate and loading-derived stress and/or stressing rate changes. Among the hydrological load components, lake water, snow, and soil moisture exhibit different annual phases, and the lake water storage dominates the total annual stress changes. Our analysis indicates that seasonal surface hydrological loads, especially lake water variations, modulate the seismicity rate along the BSFZ in central Japan.

58 GEOSCIENCES↗
DOE OSTI2024

Insights on seasonal solifluction processes in warm permafrost Arctic landscape using a dense monitoring approach across adjacent hillslopes

Solifluction processes in the Arctic are highly complex, introducing uncertainties in estimating current and future soil carbon storage and fluxes, and assessment of hillslope and infrastructure stability. This study aims to enhance our understanding of triggers and drivers of soil movement of permafrost-affected hillslopes in the Arctic. To achieve this, we established an extensive soil deformation and temperature sensor network, covering 48 locations across multiple hillslopes within a 1 km² watershed on the Seward Peninsula, AK. We report depth-resolved measurements down to 1.8 m depth for May to September 2022, a period conducive to soil movement due to deepening thaw layers and frequent rain events. Over this period, surface movements of up to 334 mm were recorded. In general, these movements occur close to the thawing front, and are initiated as thawing reaches depths of 0.4 to 0.75 m. The largest movements were observed at the top of the south-east facing slope, where soil temperatures are cold (mean annual soil temperatures averaging -1.13°C) and slopes are steeper than 15°. Our analysis highlights three primary factors influencing movements: slope angle, soil thermal conditions, and thaw depth. The latter two significantly impact the generation of pore water pressures at the thaw–freeze interface. Specifically, soil thermal conditions govern the liquid water content, while thaw depth influences both the height of the water column and, consequently, the pressure at the thawing front. These factors affect soil properties, such as cohesion and internal friction angle, which are crucial determinants of slope stability. This underscores the significance of a precise understanding of subsurface thermal conditions, including spatial and temporal variability in soil temperature and thaw depth, when assessing and predicting slope instabilities. Based on our observations, we developed a Factor of Safety proxy that consistently falls below the triggering threshold for all probes exhibiting displacements exceeding 50 mm. This study offers novel insights into patterns and triggers of hillslope movements in the Arctic and provides a venue to evaluate their impact on soil redistribution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Opening-mode fracturing and cementation during hydrocarbon generation in shale: An example from the Barnett Shale, Delaware Basin, West Texas

Relative timing of fracturing is a key input for predictive fracture models, but timing information for fractures is commonly difficult to obtain. In this study, we used crosscutting relations and fluid inclusion assemblage temperatures from fracture cements from a few well-documented sampled fractures, combined with a one-dimensional burial history model, to establish timing for three generations of opening-mode fractures in a Barnett Shale core from the southern part of the Delaware Basin, Pecos County, West Texas. A burial history model is presented for the cored well and matched to measured vitrinite reflectance in samples from the core, and bottomhole temperature in the well. Here, the earliest fractures (group 1) likely formed due to early fluid-expulsion events (ca. 300 Ma) and were folded during host-rock compaction. Later group 2 fractures are sealed with fibrous barite containing primary, liquid hydrocarbon inclusions (mean homogenization temperature [T h ] = –9°C) and aqueous fluid inclusions (mean T h = 108.1°C). Group 2 fractures likely formed in response to fluid overpressure associated with cracking of type II kerogen to oil. Group 3 vertical fractures are up to 2 m in height with kinematic apertures ranging from less than 0.05 to 1.4 mm, partly open, and strike dominantly 010°–020°. Sequentially trapped aqueous fluid inclusions in fracture-spanning quartz cement bridges (mean T h = 110°C in crack-seal texture and 128°C in post-crack-seal fracture cement) record fracture opening under increasing temperature, inferred to reflect increasing burial, with continued overpressuring during the Triassic to Late Cretaceous. Some group 3 fractures may have continued to fill during Cenozoic uplift.

02 PETROLEUM↗
DOE OSTI2020

Utah FORGE: Discrete Fracture Network (DFN) Data

The FORGE team is making these fracture models available to researchers wanting a set of natural fractures in the FORGE reservoir for use in their own modeling work. They have been used to predict stimulation distances during hydraulic stimulation at the open toe section of well 16A(78)-32. These fracture sets are fully stochastic and do not contain the deterministic set that matches the pilot well 58-32 FMI data. Well 58-32 has been completed and 16A(78)-32 is to be drilled as part of Phase 3. The original .fab files are not included due to redundancy. The *.fabgz data for the 800m and 1200m depth areas are in the native FracMan format and have been compressed using Gzip. Filtered data for the 800m depth area includes .csv spreadsheets, native FracMan (.fab), and GOCAD (.ts) files that are in a compressed zip format. The file titled "SGW 2020 Finnila and Podgorney DFN fracture files on GDR.pdf" is a description of the data and should be reviewed prior to data use.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Utah FORGE Phase 3 Native State Model: 2022 Update

This is the Phase 3 native state model update. The Phase 3 numerical model represents a significant subsurface volume below the FORGE site footprint. The model domain of 4.0 km x 4.0 km x 4.2 km is located approximately between depths of 4000 to 4200 meters below land surface. This data archive consists of 10 files, 4 of which are simulation input files and the remaining 6 are simulation output files. There is an included readme.txt file that contains details on each of the data files. The input files include meshes, FALCON code inputs, tabulated data of water properties, temperature values, and model boundaries. The output files include simulation outfiles and point data of modeled material properties.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Utah FORGE 5-2557: Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes - Workshop Presentation

This is a presentation on the Role of Fluid and Temperature in Fracture Mechanics and Coupled Thermo-Hydro-Mechanical-Chemical (THMC) Processes for Enhanced Geothermal Systems project by Purdue University, presented by Distinguished Professor of Physics & Astronomy, Laura J. Pyrak-Nolte. The project's objective was to develop and validate a macroscopic model that accounts for local deformation/frictional behavior, seismic/aseismic behavior, chemical reactions, and determine the adequacy of classic Coulomb failure vs. rate-and-state friction. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Investigation of UV, ns-laser damage resistance of hafnia films produced by electron beam evaporation and ion beam sputtering deposition methods

Laser-induced damage in coating materials with a high index of refraction, such as hafnia, limits the performance of high power and high energy laser systems. Understanding the underlying physics responsible for laser damage holds the key for developing damage-resistant optical films. Previous studies have reported a substantial difference in laser damage onset for hafnia films produced by different deposition methods, yet the underlying mechanisms for the observed difference remain elusive. We combined laser damage testing with analytical characterizations and theoretical simulations to investigate the response of hafnia films produced by electron (e-) beam evaporation vs ion beam sputtering (IBS) methods upon UV ns-laser exposure. We found that e-beam produced hafnia films were overall more damage resistant; in addition, we observed a polarization anisotropy associated with the onset of damage in the e-beam films, while this effect was absent in the latter films. The observed differences can be attributed to the stark contrast in the pressure inside the pores inherent in both films. The high pressure inside the IBS-induced nanobubbles has been shown to reduce the threshold for laser-induced plasma breakdown leading to film damage. The polarization effects in the e-beam coatings can be related to the asymmetric electric field intensification induced by the columnar void structure. Our findings provide a fundamental basis for developing strategies to produce laser damage-resistant coatings for UV pulsed laser applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Pore-scale observations of natural hydrate-bearing sediments via pressure core sub-coring and micro-CT scanning

Abstract Both intra-pore hydrate morphology and inter-pore hydrate distribution influence the physical properties of hydrate-bearing sediments, yet there has been no pore-scale observations of hydrate habit under pressure in preserved pressure core samples so far. We present for the first time a pore-scale micro-CT study of natural hydrate-bearing cores that were acquired from Green Canyon Block 955 in UT-GOM2-1 Expedition and preserved within hydrate pressure–temperature stability conditions throughout sub-sampling and imaging processes. Measured hydrate saturation in the sub-samples, taken from units expected to have in-situ saturation of 80% or more, ranges from 3 ± 1% to 56 ± 11% as interpreted from micro-CT images. Pore-scale observations of gas hydrate in the sub-samples suggest that hydrate in silty sediments at the Gulf of Mexico is pore-invasive rather than particle displacive, and hydrate particles in these natural water-saturated samples are pore-filling with no evidence of grain-coating. Hydrate can form a connected 3D network and provide mechanical support for the sediments even without cementation. The technical breakthrough to directly visualize particle-level hydrate pore habits in natural sediments reported here sheds light on future investigations of pressure- and temperature-sensitive processes including hydrate-bearing sediments, dissolved gases, and other biochemical processes in the deep-sea environment.

58 GEOSCIENCES↗
DOE OSTI2021

Tools for pressure core sub-coring and pore-scale micro-CT (computed tomography) scans

Abstract. The pore habits of gas hydrate in natural sediment matrices provide essential clues for understanding physical (mechanical, thermal, hydraulic, and electrical) properties of hydrate-bearing sediments, yet there are no tools that can directly visualize the pore habits of natural gas hydrate other than indirect interpretation based on core-scale or field-scale observations. A significant challenge is to obtain a mini-core from pressure cores retrieved from natural reservoirs for high-resolution micro-CT (computed tomography) scans while maintaining pressure and temperature conditions required for stability of gas hydrate during all operational steps including manipulation, cutting, transferring, sub-coring and CT scanning. We present a new set of tools for pore-scale micro-CT imaging of natural hydrate-bearing sediments while maintaining pressure and temperature control. The tests with laboratory-prepared cores and pressure cores successfully demonstrate the capability of this set of tools to subsample a mini-core from pressure cores, transfer the mini-core to an X-ray transparent core holder, and conduct micro-CT scans. Successfully obtained CT images prove the functionality of this set of tools.

Seol, Yongkoo↗