Search NASA⌕ Search

SEARCH · Search NASA

Results for “SATURATION”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗
DOE OSTI2026

Integration of the Biot–Gassmann Fluid Substitution Method and Machine Learning-Based Velocity–Stress Relationship for Estimating In Situ Stresses

Recent advancements have shown that in situ stresses can be reliably estimated through an integrated machine/deep learning (ML/DL)-based framework, which relies on models trained and validated using true triaxial ultrasonic velocity (TUV) experimental data that involve measurements of ultrasonic velocity in saturated rocks under varying stress configurations. However, when the goal is to interpret lower frequency measurements, it may be more appropriate to run experiments on dry rocks and then obtain Biot–Gassmann-derived equivalent saturated velocities (low-frequency approximation) and employ these quantities for training ML/DL models to predict in situ stress. Whether the dispersion effect of frequency on the velocity–stress relationship substantially impacts in situ stress prediction is an important and unresolved question. This work presents an enhancement of ML/DL-based workflow by training and implementing ML/DL models using equivalent saturated acoustic velocities (low-frequency) obtained by applying Biot–Gassmann fluid substitution on the ultrasonic velocities of dry cores. The models were trained on TUV data sets derived from three subsurface cores extracted from the geothermal well 16B(78)-32 at the Utah FORGE site. Each core was subjected to 75 unique stress configurations for velocity measurement in the dry state. The ML/DL trained on the TUV data set with equivalent saturated velocities demonstrated promising performance to predict in situ stress in subsurface geological rocks using velocity–stress relationships with R 2 of 0.86, 0.971, and 0.975 and root mean squared error (RMSE) of 2.59, 1.92, and 1.80 for validation/testing phases of vertical, minimum horizontal, and maximum horizontal stress models, respectively. Additionally, interpretation and explanation by Shapley additive explanations (SHAP) analysis further improved scientific validation and model reliability for estimating in situ stresses.

colloids↗
DOE OSTI2024

Microbial polyphenol metabolism is part of the thawing permafrost carbon cycle

Abstract With rising global temperatures, permafrost carbon stores are vulnerable to microbial degradation. The enzyme latch theory states that polyphenols should accumulate in saturated peatlands due to diminished phenol oxidase activity, inhibiting resident microbes and promoting carbon stabilization. Pairing microbiome and geochemical measurements along a permafrost thaw-induced saturation gradient in Stordalen Mire, a model Arctic peatland, we confirmed a negative relationship between phenol oxidase expression and saturation but failed to support other trends predicted by the enzyme latch. To inventory alternative polyphenol removal strategies, we built CAMPER, a gene annotation tool leveraging polyphenol enzyme knowledge gleaned across microbial ecosystems. Applying CAMPER to genome-resolved metatranscriptomes, we identified genes for diverse polyphenol-active enzymes expressed by various microbial lineages under a range of redox conditions. This shifts the paradigm that polyphenols stabilize carbon in saturated soils and highlights the need to consider both oxic and anoxic polyphenol metabolisms to understand carbon cycling in changing ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Benchmarking core turbulence and transport predictions for an inductive compact tokamak reactor plasma

Motivated by the need for accurate, timely, and efficient calculations of plasma transport, predictions of plasma turbulence properties made using different TGLF saturation rules are benchmarked against corresponding predictions from linear and nonlinear gyrokinetic CGYRO simulations. This benchmarking is carried out using parameters taken from an inductive burning plasma scenario in a hypothetical compact high-field (R maj = 4 m, B T = 8 T) tokamak, lying in a much different regime of parameter space than either the TGLF calibration regime or current-day experiments. The core turbulent transport in this scenario is predicted to be dominated by ion temperature gradient (ITG) turbulence. In general, the ITG critical gradients predicted by various TGLF saturation rules are quite close to the CGYRO predictions. Both codes predict similar linear ITG growth rates and frequency spectra, as well as their scaling with R/L T i = −Rd ln(T i )/dr. However, TGLF systematically predicts unstable trapped-electron modes (TEMs) above k y ρ s ≃ 0.5 not seen by CGYRO for the same parameters, due to TGLF predicting a lower threshold in R/L T e than CGYRO for TEM onset. It is shown that for this scenario, nonlinear CGYRO simulations predict stiffer ITG turbulence than the TGLF SAT0 and SAT1 saturation rules, with energy fluxes close in magnitude and scaling with R/L T i to what is predicted by the SAT2 saturation rule. Self-consistent core profiles calculated using nonlinear CGYRO flux predictions and the PORTALS transport solver are shown to agree fairly well with corresponding predictions made using the TGLF SAT2 model, including a similar level of density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Nonlinear magnetohydrodynamic modeling of ideal ballooning modes in high- β Wendelstein 7-X plasmas

We present nonlinear magnetohydrodynamic (MHD) simulations of high- β Wendelstein 7-X plasmas using the stellarator extension of the M3D- C 1 code, building on the recent work that shows benign saturation of ideal ballooning modes above the designed β limit in the standard configuration [Zhou et al., Phys. Rev. Lett. 133, 135102 (2024)]. First, we examine the results' sensitivity to the parallel thermal conductivity. It is found that while an increased parallel conductivity reduces the linear growth rate, the saturated pressure profile is barely affected. Second, we consider the dependence on the profile shape. It is shown that an equilibrium with a peaked pressure profile and lower β is subject to more significant change than a broad profile with higher β and a larger growth rate, suggesting that benign saturation, or nonlinear stability, is not guaranteed and not dictated by linear growth. Third, we study the influence of the magnetic configuration, with the equilibrium rotational transform varied by adjusting the planar coil current. With similar growth rates, similar magnitudes of profile change are found regardless of the presence of a low-order resonance, which implies that the saturation mechanism is not specific to a resonant or non-resonant mode. These results indicate that MHD stability should still be treated seriously in stellarator operation and design, for which nonlinear modeling using tools like M3D- C 1 can play an instrumental role.

Zhou, Yao [Shanghai Jiao Tong University (China)] ↗
DOE OSTI2025

Loading Capacity and Dilute Nitric Acid Rinse of Diglycolamide Resin for the Recovery of Transplutonium and Rare Earth Elements from Mark-18A Targets

N,N,N′,N′-tetraoctyldiglycolamide (TODGA) as a resin (“DGA resin”) produced by Eichrom Technologies will be used by Savannah River National Laboratory for the indiscriminate extraction of trivalent actinides and rare earth elements with the intent of recovering Cm and Am from dissolved irradiated 242 Pu Mark-18A targets in 7 to 9 M nitric acid. The extracted constituents will be recovered as an oxide by direct thermal decomposition of the loaded resin followed by calcination of the resultant residue. The characteristics of DGA resin with non-radiological feed simulant representative of the anticipated feed in the Mark-18A process, with Sm and Nd as surrogates for Cm and Am, respectively, including breakthrough point and saturation capacity were evaluated in this work. Additionally, this work examined the losses from the loaded resin by rinsing the resin bed with dilute acid to reduce the nitrate concentration in the resin bed prior to thermal decomposition operations to improve the safety posture of the process. A resin loading profile was developed, and the resin was determined to have a trivalent metal saturation capacity of 74 μmol/mL resin under the experimental conditions evaluated. Following a wash of the loaded resin bed with fresh 8 M HNO 3 , the trivalent metals retained was reduced to 68 μmol/mL resin, which represents the practical capacity of the resin for the Mark-18A process. Lighter lanthanides breakthrough the resin well before the saturation capacity is reached. Rinsing the saturated resin bed with 0.26 M HNO 3 was found to result in a rapid reduction in retention of rare earth elements by the resin. After 2.8 bed volumes of dilute acid rinse, the mean resin bed free acid concentration was reduced to 0.28 M and 3.1 bed volumes of dilute acid rinse resulted in a reduction of the cumulative rare earth element retention to 52 μmol/mL of resin.

Transplutonium separations↗
DOE OSTI2024

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

Enhanced Carbon Storage Forecasting via Cross-Geology Transfer Learning

Rapid simulation of the spatiotemporal evolution of pressure & saturation for SACROC 1. Neural operator was trained on only 153 simulation runs 2. Trained to account for heterogeneity and variations/uncertainties in engineering, fluids, and geology 3. Pressure forecast has less than 1% error 4. Saturation forecast has less than 2% error 5. Traditional simulator takes 1 hour for a single scenario, while neural operator takes less than 1 minute. Rapid simulation of the spatiotemporal evolution of pressure & saturation for IBDP 1. Transfer Learning was implemented on the SACROC-based Neural Operator that was trained on only 20 simulation runs for IBDP Site 2. SACROC and IBDP Sites have several significant differences in geology and engineering parameters. 3. Pressure forecast has less than 5 psi error 4. Saturation forecast has less than 7% error 5. Traditional simulator takes 1 hour for a single scenario, while neural operator takes less than 1 minute and only 20 simulations for training/validation.

Misra, Siddharth↗
DOE OSTI2025

Expedition UT-GOM2-2 Site H

Pressure and conventional cores were collected at Site H of the Walker Ridge Protracted Area Block 313 in the Terrebonne Basin, deepwater Gulf of America (Gulf of Mexico) during the University of Texas (UT) Deepwater Hydrate Coring Expedition (UT-GOM2-2). Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were recovered. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of about 200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. These ooze intervals also correspond to lighter sediment color, increased Ca content based on X-ray florescence (XRF) core scanning, and increased calcareous nannofossil abundance. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years), with a pronounced increase in sedimentation rate with depth. Below 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C1 /C2 ) and the methane to ethane plus propane (C1 /(C2 +C3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ13C isotopic signature of methane ranges between -69.9 and -78.5 ‰ relative to the Vienna Pee Dee Belemnite (VPDB) standard. Pressure core recovery in sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One pressure core was degassed and the average hydrate saturation in the core was determined to be 24%. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗
DOE OSTI2024

Coherent γ *-nucleus scattering and coherent nuclear states

In the context of a McLerran-Venugopalan (MV) model for a large nucleus, coherent scattering of a virtual photon on that nucleus is evaluated in the A − = 0 gauge, the gauge appropriate for the target nucleus. The evaluation of the scattering in A − = 0 gauge is very intricate compared to the usual A + = 0 gauge evaluation natural for the scattering process, but has the advantage of directly giving the scattering in terms of a partonic description of the nucleus. In the limit where a tagged forward jet puts the dipole-nucleus scattering in the saturation regime the coherent reactions are equal to the inelastic reactions. In terms of the nuclear wave function the coherent reactions come from color singlet and zero total transverse momentum quark-antiquark pairs in the wave function and in the saturation regime the nuclear wave function is a coherent state for these pairs. In the saturation region half of all quarks (or antiquarks) come from zero momentum and color charge pairs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Suppression of diffraction in deep-inelastic scattering on nuclei and dynamical mechanism of leading twist nuclear shadowing

Abstract Using the leading twist approach (LTA) to nuclear shadowing, we calculate the ratios of diffractive and usual parton distributions for a heavy nucleus (Pb) and the proton,$$ {R}_{A/p}=\left({f}_{i/A}^{D(3)}/{f}_{i/A}\right)/\left({f}_{i/p}^{D(3)}/{f}_{i/p}\right) $$ R A / p = f i / A D 3 / f i / A / f i / p D 3 / f i / p , for coherent and summed (coherent plus quasi-elastic) nuclear deep-inelastic scattering. We find thatR A/p ≈ 0.5 − 1 for quarks as well as for the ratio of the diffractive and total cross sections$$ {\left[\left({d\sigma}_{\textrm{diff}}/{d M}_X^2\right)/{\sigma}_{\textrm{tot}}\right]}_{eA}/{\left[\left({d\sigma}_{\textrm{diff}}/{d M}_X^2\right)/{\sigma}_{\textrm{tot}}\right]}_{ep} $$ dσ diff / dM X 2 / σ tot eA / dσ diff / dM X 2 / σ tot ep andR A/p ≈ 0.5 − 1.3 for gluons in a broad range ofx, including the kinematics of the Electron-Ion Collider, which reaffirms the difference from the nuclear enhancement ofR A/p predicted in the gluon saturation framework. We demonstrate that the magnitude ofR A/p is controlled by the cross section of the interaction of hadronic fluctuations of the virtual photon with target nucleons, which explains an enhancement ofR A/p in the color dipole model and its suppression in LTA. We argue that the black disk limit corresponds toR A/p = 1 and$$ {R}_{A/p}^{\textrm{coh}} $$ R A / p coh = 0.86 for the summed and coherent scattering, respectively. Relying on an intuitive definition of the saturation scale, we show that the ratio of the saturation scales of a heavy nucleus and proton$$ {Q}_{sA}^2(b)/{Q}_{sp}^2(b)\approx 1 $$ Q sA 2 b / Q sp 2 b ≈ 1 at small impact parametersbdue to the strong leading twist nuclear shadowing and diluteness of the nuclear density.

Physics↗
DOE OSTI2024

Not all that is β0 is β-function: the DGLAP resummation and the running coupling in NLO JIMWLK

Abstract We reanalyze the origin of the large transverse logarithms associated with the QCD one loopβfunction coefficient in the NLO JIMWLK Hamiltonian. We show that some of these terms are not associated with the running of the QCD coupling constant but rather with the DGLAP evolution. The DGLAP-like resummation of these logarithms is mandatory within the JIMWLK Hamiltonian, as long as the color correlation length in the projectile is larger than that in the target. This regime in fact covers the whole range of rapidities at which JIMWLK evolution is supposed to be applicable. We derive the RG equation that resums these logarithms to all orders inα s in the JIMWLK Hamiltonian. This is a nonlinear equation for the eikonal scattering matrixS(x). We solve this equation, and perform the DGLAP resummation in two simple cases: the dilute limit, where both the projectile and the target are far from saturation, and the saturated regime, where the target correlation length also determines its saturation momentum.

Physics↗
DOE OSTI2025

Most general EFTs from spurion analysis Hilbert series and minimal lepton flavor violation

We derive a saturation theorem for general Effective Field Theories (EFTs) constructed using spurion analysis. Let S be a set of spurion fields introduced to organize the breaking of a global symmetry G f , and H S be the subgroup of G f that remains unbroken under a generic vacuum expectation value ⟨S⟩; we show that the EFT Lagrangian constructed from the spurion analysis saturates the EFT Lagrangian without the spurions but restricted to H S invariance, provided that arbitrary powers of the spurion fields are allowed. As examples, we study several implementations of the Minimal Lepton Flavor Violation (MLFV) principle, corresponding to various origins of the neutrino masses. In each scenario, we compute the Hilbert series to obtain the numbers of independent lepton flavor covariants that appear in the corresponding EFT at mass dimension 6. These numbers agree with the number of H S invariants in the EFT without the spurions, demonstrating the saturation theorem.

Effective Field Theories↗
DOE OSTI2025

Exploring the impact of hydro-mechanical coupling on compaction localization

In deep subsurface environments, where porous rocks are saturated with fluids, hydromechanical coupling plays a critical role in the formation and propagation of localized deformation zones, such as compaction bands. In this study, we numerically investigate the influence of fluid-rock interactions on the development of compaction localization in fluid-saturated rocks. For this purpose, it uses a numerical framework that adapts the heat equation to model the diffusion of pore pressure caused by deformation. This is achieved with specialized subroutines for material behavior and heat transfer formulated within the ABAQUS platform and by allowing changes in rock permeability during strain localization. The results highlight that hydromechanical coupling significantly alters the compaction behavior compared to dry conditions. Under high confining stress, excess pore pressure inhibits the progression of compaction fronts, while in time-dependent processes like rock creep, pore pressure transients accelerate deformation by creating feedback loops between mechanical forces and fluid flow. These findings underscore the importance of saturation in controlling the dynamics of compaction localization and provide insights into the interplay between pore pressure, deformation, and permeability. The implications of this study extend to applications in geological storage, resource extraction, and subsurface geomechanics, offering a deeper understanding of coupled hydromechanical processes in porous rocks.

58 GEOSCIENCES↗
DOE OSTI2026

Forest aboveground biomass estimation through integration of sentinel-2 and PALSAR-2 time series: assessing models trained on GEDI and field inventory benchmarks

Accurate and spatially explicit forest Aboveground Biomass (AGB) mapping through remote sensing is critical for quantifying terrestrial carbon stocks and informing effective forest management strategies. However, AGB estimation in dense forests with complex terrain remains challenging due to satellite sensor signal saturation problem (saturation issue occurs in high biomass forests), structural complexity, and limited ground truth for calibration. This study presents a novel framework that integrates multi-temporal Sentinel-2 optical imagery, ALOS PALSAR-2 Synthetic Aperture Radar (SAR) data, and topographic variables with explainable Machine Learning to map AGB across mountainous forests within subtropical and temperate oceanic climate zones of Mexico. We evaluate the effects of temporal granularity and sensor synergy by comparing multiple temporal inputs and sensor configurations (Sentinel-2, PALSAR-2, and their fusion), and assess model performance using two reference datasets: NASA GEDI LiDAR-derived biomass and Mexico’s National Forest and Soil Inventory (INFyS). Our results showed that models trained on INFyS consistently outperformed those trained on GEDI, highlighting limitations in GEDI’s reliability in biomass estimates within this study region. Furthermore, the integration of Sentinel-2 and PALSAR-2 provided improved predictions compared to single-sensor models, particularly when combined with temporally explicit yearly statistics. The best-performing model, which was trained on INFyS data, and considered both Sentinel-2 and PALSAR-2 yearly statistics, as well as topographic variables, achieved an R2 of 0.64, RMSE of 51.10 Mg/ha, and relative RMSE (rRMSE) of 58.69%. Explainable ML analysis identified Sentinel-2 spectral indices and topographic features as key predictors, while PALSAR-2 metrics provided complementary information, partially mitigating saturation effects in high-biomass areas. Specifically, integrating both sensors substantially improved AGB estimation in high biomass forest (≥200 Mg/ha), yielding 98% gains over optical-only model, with resulting estimates exceeding GEDI L4B by 29% and ESA-CCI-BIOMASS by 174%. Terrain-stratified analysis indicated close agreement with GEDI in low-slope areas, with increasing divergence as slope steepness increased, while estimates remained consistently higher than ESA-CCI-BIOMASS across all slope classes. The proposed approach advances multi-sensor fusion and temporal feature engineering for AGB mapping using open-access satellite datasets, providing a scalable and reproducible framework for annual biomass monitoring in topographically complex mountainous forests. The resulting 25 m resolution biomass product has the potential to provide spatially detailed information for forest monitoring and may support applications in carbon accounting and forest management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A novel conditional generative model for efficient ensemble forecasts of state variables in large-scale geological carbon storage

Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.

58 GEOSCIENCES↗
DOE OSTI2024

Impact of supercritical carbon dioxide on the frictional strength of and the transport through thin cracks in shale

Understanding the mechanical and transport behavior of thin (i.e. small aperture) cracks slipping under supercritical carbon dioxide (sc-CO 2 ) conditions is essential to evaluate the integrity of sealing formations with buoyant sc-CO 2 below and the success of waterless fracturing. The two major items of interest in this work are frictional strength and permeability change of the crack. We used a triaxial cell that permits in situ visualization to conduct and monitor slippage along the faces of narrow cracks subjected to triaxial stresses. Such cracks are analogs to small geological faults. We tested carbonate-rich, 1-inch diameter Wolfcamp shale samples that are saw cut 30° to vertical to create a thin crack. Friction coefficients ranged from about 0.6 to 0.8 consistent with expectations for brittle rocks. The sc-CO 2 generally did not alter friction coefficient over the time scale of experiments. From a transport perspective, saturating cracks with sc-CO 2 substantially decreased permeability of the crack by 26%–52%, while slip resulted in a variety of permeability responses. Overall, the combined impact of sc-CO 2 saturation and slip reduced fault permeability for all tests. Our observations support the notion that the sealing capacity of some caprocks improves when saturated with sc-CO 2 and that some slip of small fractures is not necessarily detrimental to caprock integrity.

58 GEOSCIENCES↗
DOE OSTI2025

Enhanced pool boiling heat transfer using novel inverted trapezoid microgroove structures

Pool boiling heat transfer improvements with enhanced surfaces have been extensively studied for various hydrofluorocarbon (HFC) refrigerants. However, there is a limited number of studies in the open literature that focus on the pool boiling heat transfer characteristics of hydrofluoroolefin (HFO) refrigerants on enhanced surfaces. This study aims to address this gap by analyzing the pool boiling behavior of HFO refrigerants, specifically R1234yf and R1234ze(E). Experiments were conducted on both plain surfaces and novel inverted trapezoid microgroove structures. Further, the experimental results indicate that the microgroove structures can improve heat transfer coefficients by up to 80% compared to plain surfaces. Additionally, the performance of R1234yf and R1234ze(E) was compared to that of R134. It was observed that R1234yf exhibited comparable performance to R134a, while R1234ze(E) showed slightly lower performance. To assess the influence of saturation temperature on pool boiling heat transfer, tests were conducted at different saturation temperatures of 15℃, 25℃, and 35℃. The findings revealed that the pool boiling heat transfer coefficient increases with rising saturation temperature.

42 ENGINEERING↗