Search NASASearch

SEARCH · Search NASA

Results for “fluid structure interaction”

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.

89 records · Page 5

Pressure Spike in The LBNF Absorber Core s Gun Drilled Cooling Channel from an Accident Beam Pulse

The LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. The cooling water is supplied by the LBNF Absorber Radioactive Water (RAW) cooling system which is designed as per ASME B31.3 Normal Fluid Service [1]. An uninhibited beam accident pulse striking the water channels was identified as a credible accident scenario. In this study, it is assumed that the beam pulse hits the Absorber directly without interacting with any of the other upstream beamline components. The beam parameters used for the LBNF beam are 120 GeV, 2.4 MW with a 1.2 s cycle time. The accident pulse lasts for 10 μs. The maximum energy is deposited in the 3rd aluminum core block. For the sake of simplicity, it is assumed that the accident pulse strikes the 1 in. ID water channel directly. The analysis here simulates the pressure rise in the water during and after the beam pulse and its effects on the aluminum piping components that deliver water to the core blocks. The weld strengths as determined by the Load and Resistance Factor Design (LRDF) and the Allowable Strength Design (ASD) are compared to the forces generated in the weld owing to the pressure spike. A transient structural analysis was used to determine the equivalent membrane, peak, and bending stresses and they were com-pared to allowable limits.

43 PARTICLE ACCELERATORS

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Advanced System Thermal Fluids Solver Development for SAM

This work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key is the implementation of a high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. Leveraging the existing capabilities of the SAM code, significant code coverages were established in the finite volume method code. This in turn allows for a suite of test problems with different problem sizes and levels of complexity to be used to quantify the performance improvement of the finite volume method code. As evidently shown in this study, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE for the wide range of selected problems. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. In addition, for a complex reactor model, transient simulation was performed using the finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development. In this work, short-term priority development and testing items were identified, and long-term code adoption and integration plans were made for the eventual deployment of the finite volume method in the SAM code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE

Elucidating the Transition of 3D Morphological Evolution of Binary Alloys in Molten Salts with Metal Ion Additives

Molten salts serve as effective high-temperature heat transfer fluids and thermal storage media used in a wide range of energy generation and storage facilities, including concentrated solar power plants, molten salt reactors and high-temperature batteries. However, at the salt–metal interfaces, a complex interplay of charge-transfer reactions involving various metal ions, generated either as fission products or through corrosion of structural materials, takes place. Simultaneously, there is a mass transport of ions or atoms within the molten salt and the parent alloys. The precise physical and chemical mechanisms leading to the diverse morphological changes in these materials remain unclear. Here, to address this knowledge gap, this work employed a combination of synchrotron X-ray nanotomography and electron microscopy to study the morphological and chemical evolution of Ni-20Cr in molten KCl-MgCl 2 , while considering the influence of metal ions (Ni 2+ , Ce 3+ , and Eu 3+ ) and variations in salt composition. Our research suggests that the interplay between interfacial diffusivity and reactivity determines the morphological evolution. The summary of the associated mass transport and reaction processes presented in this work is a step forward toward achieving a fundamental comprehension of the interactions between molten salts and alloys. Overall, the findings offer valuable insights for predicting the diverse chemical and structural alterations experienced by alloys in molten salt environments, thus aiding in the development of protective strategies for future applications involving molten salts.

36 - MATERIALS SCIENCE

AMR-Wind: A Performance-Portable, High-Fidelity Flow Solver for Wind Farm Simulations

We present AMR-Wind, a verified and validated high-fidelity computational-fluid-dynamics code for wind farm flows. AMR-Wind is a block-structured, adaptive-mesh, incompressible-flow solver that enables predictive simulations of the atmospheric boundary layer and wind plants. It is a highly scalable code designed for parallel high-performance computing with a specific focus on performance portability for current and future computing architectures, including graphical processing units (GPUs). In this paper, we detail the governing equations, the numerical methods, and the turbine models. Establishing a foundation for the correctness of the code, we present the results of formal verification and validation. The verification studies, which include a novel actuator line test case, indicate that AMR-Wind is spatially and temporally second-order accurate. The validation studies demonstrate that the key physics capabilities implemented in the code, including actuator disk models, actuator line models, turbulence models, and large eddy simulation (LES) models for atmospheric boundary layers, perform well in comparison to reference data from established computational tools and theory. We conclude with a demonstration simulation of a 12-turbine wind farm operating in a turbulent atmospheric boundary layer, detailing computational performance and realistic wake interactions.

17 WIND ENERGY

Reactivity of Shale to Supercritical CO 2 : Insights from Microstructural Characterization and Mineral Phase Evolution in Caney Shales for CCUS Applications

Understanding mineral–fluid interactions in shale under supercritical CO 2 (scCO 2 ) conditions is relevant for assessing long-term geochemical containment. This study characterizes mineralogical transformations and elemental redistribution in five Caney Shale samples serving as proxies for reservoir (R1, R2, R3) and caprock (D1, D2) facies, subjected to 30-day static exposure to pure scCO 2 at 60 °C and 17.23 MPa (2500 psi), with no brine or impurities introduced. SEM-EDS analyses were conducted before and after exposure, with mineral phases classified into silicates, carbonates, sulfides, and organic matter. Initial compositions were dominated by quartz (38–47 wt.%), illite (16–23 wt.%), carbonates (12–18 wt.%), and organic matter (8–11 wt.%). Post-exposure, carbonate loss ranged from 15 to 40% in reservoir samples and up to 20% in caprock samples. Illite and K-feldspar showed depletion of Fe 2+ , Mg 2+ , and K + at grain edges and cleavages, while pyrite underwent oxidation with Fe redistribution. Organic matter exhibited scCO 2 -induced surface alteration and apparent sorption effects, most pronounced in R2 and R3. Elemental mapping revealed Ca 2+ , Mg 2+ , Fe 2+ , and Si 4+ mobilization near reactive interfaces, though no secondary mineral precipitates formed. Reservoir samples developed localized porosity, whereas caprock samples retained more structural clay integrity. The results advance understanding of mineral reactivity and elemental fluxes in shale-based CO 2 sequestration.

geochemical CO2 sequestration

Nanorod Pair Complexes Manipulated via Magnetic Casimir Forces

Here, we demonstrate a route to modulate Casmir–Lifshitz interactions between anisotropic nanoparticles and magnetic fluids. By semiclassical quantum electrodynamics, we study ground state dispersion forces for cylindrical dielectric nanorods made of polystyrene (PS) and zinc oxide (ZnO) embedded in toluene-based host media with gold-coated magnetite nanoparticles and also predict magnetic contributions to the fully retarded excited state interaction. The variation in magnetic permeability enables tuning between repulsive and attractive interactions, and measurable magnetic Casimir traps are predicted between a pair of ZnO–PS nanoparticles whose equilibrium position can be modulated over an order of magnitude with a small variation in the size of the magnetite nanoparticle. This provides an alternative magnetic Casimir-effect pathway to reversibly tune quantum electromagnetic forces at the nanoscale for the assembly and enhancement of colloidal stability.

Excited-state interaction

TCF - Sustainable Well Cement for Geothermal, Thermal Recovery, and Carbon Storage wells

Primary cementing is the most important operation performed on a subterranean well. Cement, placed in the annulus between the casing and the formations, serves as a hydraulic seal preventing fluids and gas migrations, protecting steel casing from corrosion, and supporting the well structure. Poor cementing jobs can be directly responsible for wells not reaching their full capacity, casing corrosion, compromised well integrity and, in the worst-case scenarios, well collapse. Geothermal, thermal recovery, and carbon storage wells offer environments that are especially difficult for cements to survive while the required lifespan of these wells can be years. In these wells currently used cements cannot provide durable well integrity, and new solutions are necessary. The necessity of stabilizing the electric grid, increasing its flexibility, and providing energy on demand will further expand the market of durable cement solutions for applications in high-temperature underground reservoir thermal energy storage wells (HT RTES). The design challenges of special cement systems for such wells originate from chemistry limitations of currently used well cements, aggressive environments, temperature and repeated shock conditions associated with them, and very weak formations that are not uncommon in these wells. During the last 70 years the most common systems for thermal-well construction have been Ordinary Portland Cement (OPC) (absolute majority of the wells), silica-lime system, and high aluminum cement-containing formulations. The major issues of calcium-silicate hydrates, that form during the hydration of OPC resulting in hardened material, are their poor chemical resistance to acids due to the calcium interactions with acid anions followed by its eventual dissolution, even in mild acids, such as carbonic acid from CO 2 dissolution, and inadequate bonding to the steel causing serious casing corrosion problems (Sugama & Pyatina, 2019). Performance of OPC-based cement may be significantly improved with organic additives; however, those are limited by their temperature stability. An alternative to calcium-silicate cement was developed by BNL and commercialized by Halliburton as ThermaLock TM cement. Calcium-aluminate-phosphate-based chemistry of the material allowed overcoming acid-resistance problems of OPC, especially in CO 2 -rich geothermal environments. In 2012-2015, building upon this experience, BNL developed Thermal Shock Resistant Cement (TSRC) with the support of Geothermal Technology Office (GTO) of DOE. TSRC has high-temperature stable chemistry, superior properties of cement-metal casing bond and corrosion protection (4-times better than currently used OPC-based formulation), significantly outperforms common well-cements in thermal-shock tests, has self-healing properties, and like ThermaLock TM is CO 2 resistant incorporating carbonate ions into stable hydrates (Gill et al., 2012).

15 GEOTHERMAL ENERGY

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Thermodynamic and Kinetic Modulation of Methylammonium Lead Bromide Crystallization Revealed by In Situ Monitoring

Hybrid organic–inorganic perovskite (HOIP) crystals are promising optoelectronic materials, but little is known about either the thermodynamic and kinetic controls on crystal growth or the underlying growth mechanism(s). Herein, we use fluid cell atomic force microscopy (AFM) and solution nuclear magnetic resonance (NMR) spectroscopy to investigate growth of the model HOIP crystal CH 3 NH 3 PbBr 3 (MAPbBr 3 ) and to determine how formic acid (HCOOH) modulates the thermodynamics and kinetics of growth. The results show that growth of MAPbBr 3 in dimethylformamide (DMF) proceeds through the classical pathway by the spreading of atomic crystal steps generated at screw dislocations on the {100} surface. Temperature dependent step velocity measurements demonstrate that with increasing concentration, HCOOH decreases both the solubility of MAPbBr 3 and the kinetic coefficient (b) of step movement. 1 H-NMR measurements indicate that HCOOH increases the lifetime of the methylammonium (MA + ) ions and promotes the association of MAPbBr 3 , thus tuning the solubility of the perovskite. HCOOH also alters the molecular tumbling motion and bulk diffusion of the MA + ions, possibly via H-bonding. Further, our findings establish a direct correlation between the mesoscale crystal growth kinetics and the molecular-scale interactions between organic additives and constituent ions, providing unprecedented insights for developing predictive syntheses of HOIP crystals with defined size, crystal habit and shape, and defect distribution.

36 MATERIALS SCIENCE

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS

Results of the second Ice Shelf–Ocean Model Intercomparison Project (ISOMIP+)

Ocean-driven basal melting of Antarctic ice shelves plays an important role in the mass loss of the Antarctic Ice Sheet. Ice shelf cavity-resolving ocean models are a valuable tool for understanding ice shelf-ocean interactions and for simulating projections of ice shelf and ocean states under future climate. Designed to assess the current state of ice shelf–ocean modelling, the second Ice Shelf–Ocean Model Intercomparison Project, ISOMIP+, consists of 12 ocean model configurations submitted with a common, idealised experimental setup. Here, we focus on the experiments Ocean0–2 (Asay-Davis et al., 2016), which are ocean models with idealised, static ice shelf geometries, but where the ocean reaches a balance with prescribed far-field ocean conditions. Different thermal transfer coefficient values (ranging from 0.011 to 0.2) are used for each model in the melting parameterisation to achieve a common, tuned melt rate since the models cover a range of types of vertical coordinates, ice–ocean boundary layer treatments, and numerical schemes. These model differences lead to spread in the resultant ocean properties, circulation, boundary-layer structure and spatial distribution of melting. We also highlight similarities between models, such as a shared linear relationship across most models between melt rate and overturning and barotropic streamfunctions during the spin-up and spin-down, demonstrating a robust relationship between melt and circulation across models and forcing conditions. The ISOMIP+ results provide a systematic comparison of ice shelf cavity-capable ocean models. However, we also demonstrate the need for realistic ice shelf–ocean model intercomparison projects (some already underway) to assess model biases and inter-model variation against sparse observations. Further research is needed to understand the differences between models and further improve our modelled representations of the ice–ocean boundary layer and ice shelf cavity circulation.

54 ENVIRONMENTAL SCIENCES

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING

Williston Basin Resource Study for Commercial-Scale Subsurface Hydrogen Storage

The Energy & Environmental Research Center (EERC), in partnership with the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL), the EERC’s State Energy Research Center (SERC), MPLX Operations LLC, and TC Energy Development Holdings Inc. (a subsidiary of TC Energy Corporation), studied the potential for subsurface hydrogen storage and recovery in the Williston Basin of western North Dakota. The project’s goal was to evaluate the feasibility of large-scale, secure geologic H 2 storage to support future hydrogen generation, storage, and use. This work included laboratory testing, H 2 –rock–fluid exposure experiments, literature reviews on H 2 embrittlement, and reservoir modeling and simulations. The study included an assessment of storage potential across three types of storage reservoirs using both reservoir simulation and DOE’s web-based tool SHASTA-HELP (Subsurface Hydrogen Assessment, Storage, and Technology Acceleration – Hydrogen Estimator for Logistical Planning), as well as investigation of potential H 2 production and markets for commercial-scale deployment. Building on prior EERC gas storage research, three storage options were selected for detailed evaluation: the Broom Creek Formation (a clastic saline reservoir), the Dickinson Lodgepole Mounds (DLM) complex (carbonate mud mound structures) of the Lodgepole Formation (an active oil and gas producing reservoir), and the Dunham Salt Interval of the Piper Formation (to be used for engineered salt cavern development). These targets were prioritized based on prior EERC research using datasets related to seal capacity, reservoir quality, mechanical integrity, and injectivity. Exposure tests on Broom Creek and DLM samples showed mineral dissolution and precipitation that increased brine salinity and altered reservoir rock surfaces. Although these results provide useful insight, they are limited by small sample sizes and short-term (30-day) exposure, requiring further study to assess long-term storage integrity. Salt formations were not tested because of their known nonreactivity and established mechanical stability. Results of reservoir simulations performed for a single site demonstrated that the Broom Creek Formation may be capable of receiving up to 42,000 tonnes of injected H 2 over 7 months via one well. H 2 recovery took place over 5 months, resulting in approximately 26,000 tonnes (~62% without cushion gas [CG]). This work suggests water production may be important and subsequent cycles of injection and production may perform more efficiently; however, significant site-specific work in the future is needed to assess actual reservoir performance of injection and withdrawal of H 2 storage. For oil reservoir potential, a multiple-well model was used to simulate injection of approximately 32,000 tonnes of H 2 into a single wellbore while simultaneously producing in place reservoir fluids from four offset wells to maintain reservoir pressure. The simulation results suggested a high recovery (~98%); in addition, cost advantages through existing infrastructure could be realized. Challenges in this reservoir include vi managing gas purity and leakage risks. In both scenarios, production of H 2 takes place in a single-well scenario with 10 cycles (7 months of injection and 5 months of production) over 10 years. Finally, the use of engineered caverns in the Dunhan Salt was evaluated, and the results suggest that while they have a smaller capacity (<1000 tonnes per cavern), they exhibit nearly complete gas recovery (>99%), fast response times, and low purity risk. While caverns in North Dakota may be smaller in capacity, fields can be developed in galleries to accommodate the volumetric needs and rapid turnaround times necessary to meet market demands. Geographic limitations and thin salt intervals in North Dakota may represent less total storage potential than salt domes elsewhere, but significant opportunities exist to expand this market for gas storage in North Dakota. A basinwide assessment was performed to estimate a first-of-its-kind value for H 2 storage on a large scale. DOE’s SHASTA-HELP, combined with EERC simulation work, was used to perform the assessment. Estimated H 2 storage potential varied widely for each formation type. The Broom Creek saline formation was estimated to have a storage potential of approximately 1.7–90.5 million tonnes (MMt). The DLM oil reservoirs were estimated to have 0.07–0.19 MMt of capacity. Notably, each of these estimates relies on significant assumptions regarding reservoir thickness, porosity, permeability, and CG needed for operation. Much research is needed to understand the true site-specific storage resource potential of each formation. Using the Dunham Salt Interval for cavern development may result in as much as 4.79 MMt (up to 2.87 MMt working gas) of H 2 storage potential. An important note for consideration is that the values presented here need significantly more geological characterization and engineering assessments prior to gaining confidence in performance. This will be a focal point for future research and development needs. The basinwide evaluation also indicated that North Dakota has significant H 2 generation potential, with estimates up to a possible 13 MMt annually, suggesting a substantial opportunity for H 2 market development and thus the need for commercial-scale H 2 storage to facilitate growth and resilience. Key Recommendations 1. Conduct detailed site characterization (3D seismic, well logs, core sampling) to reduce geologic uncertainty. 2. Perform techno-economic analyses incorporating market, regulatory, and incentive frameworks. 3. Investigate long-term interactions among H 2 , CGs, well materials, and formations to assess risks. 4. Develop pilot- and field-scale demonstrations to validate models and establish best practices.

03 NATURAL GAS