SEARCH · Engineering Papers

Results for “Dynamic Reservoir Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Seawater Acidification and Bubble Plume Dispersion from Accidental Subsea CO 2 Pipeline Rupture: A Multiphase CFD Study

If a CO 2 reservoir or transmission pipeline were to leak, both the surrounding ecology and maritime traffic safety could be put at risk. To better understand and prepare for this risk, multiphase Computational Fluid Dynamics (CFD) models were built in ANSYS Fluent to capture the behavior of a leak once it enters the water. A 3D Eulerian–Eulerian model was used for validation, while a simplified 2D model was applied to simulate conditions at a 50-m depth. The models integrate bubble dynamics, gas holdup, CO 2 dissolution, dissolved species transport, and seawater acidification into a unified CFD framework. Mass transfer was calculated using the Hughmark correlation, and local seawater temperature and salinity were factored in to determine dissociation behavior and the relevant Henry’s Law constant. To confirm the 3D model’s accuracy, results were checked against two experimental datasets: the QICS field study and the Hauser Tank experiments. The team also modeled a hypothetical release scenario at the High Island 10L site and compared the results with earlier published work. The results show that at a depth of 50 m, the surrounding water column can completely absorb a CO 2 release at a rate of 35 kg/s, since the gas dissolves into the seawater as it rises toward the surface. Beyond confirming this mitigation capacity, the simulations shed light on how a leak would actually unfold in the environment, including the shape and movement of the rising bubble plume, how much CO 2 dissolves along the way, and the resulting shifts in seawater pH and pCO 2 . Together, this provides a practical framework for assessing how CO 2 leaks could affect marine environments in the Gulf of Mexico.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Revisiting the Role of Ocean Circulation Changes in Polar Ocean Heat Transport Anomalies under Global Warming

In response to greenhouse gas forcing, climate models predict that poleward ocean heat transport (OHT) weakens in the Southern Ocean but increases in the Arctic. The role of ocean circulation changes in this OHT response has been evaluated by decomposing OHT anomalies into a dynamic component (holding ocean temperature fixed while circulation evolves) and a thermodynamic component (holding ocean circulation fixed while temperature evolves). However, ocean temperature changes are themselves shaped by circulation changes through redistribution of the existing heat reservoir and subsequent effects on air–sea heat fluxes. The thermodynamic component can therefore be influenced by circulation changes, making the standard thermodynamic–dynamic decomposition incomplete for isolating the role of circulation changes in OHT anomalies. To address this issue, we use a passive–active decomposition to assess the relative contributions of ocean circulation and passive ocean temperature changes to polar OHT anomalies in a fully coupled climate model. Passive temperature changes are defined as those thermally forced by the atmosphere in the absence of circulation changes. In this passive–active decomposition, an advective term involving both circulation and passive temperature changes remains ambiguous—classifying it as active implies circulation changes dominate Southern Ocean OHT anomalies, whereas classifying it as passive implies temperature changes dominate. However, both interpretations imply that ocean circulation changes have a much weaker effect on polar OHT anomalies than inferred from the standard decomposition. In conclusion, these results help reconcile conclusions from studies using the standard decomposition with those using passive tracer methods to assess the role of circulation changes in polar OHT anomalies.

Arctic↗

DOE OSTI2025

Characterization of kerogen nanopores using 2D NMR relaxation and MD simulations

The characterization of kerogen nanopores is crucial for predicting the geostorage capacity and recoverability of natural gas in unconventional gas shale reservoirs. Towards this end, a powerful technique is presented which integrates 2D NMR T 1 -T 2 relaxation measurements with molecular dynamics (MD) simulations of hydrocarbons confined in the nanopores of kerogen. The integrated NMR-MD technique is demonstrated using T 1 -T 2 measurements of kerogen isolates and organic-rich chalks saturated with heptane, together with MD simulations of heptane completely dissolved in a realistic kerogen model. The NMR-MD results are used to extract the swelling ratio and nanopore size distribution of kerogen as a function of depth in the reservoir. The effects of organic nanoconfinement on the T 1 relaxation dispersion and T 2 residual dipolar coupling of heptane are investigated, as well as the effect of downhole effective stress on the kerogen nanopore size as a function of depth and compaction. Potential applications in partially depleted gas shale reservoirs are discussed, including CO 2 utilization/geostorage, geostorage of green H 2 , and integration of the NMR-MD technique with thermodynamic models for predicting the competitive sorption of gas mixtures in kerogen.

Compaction↗

DOE OSTI2026

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗

DOE OSTI2025

Disentangling climate and policy uncertainties for the Colorado River post-2026 operations

Abstract Lakes Mead and Powell in the Colorado River Basin underpin water and hydropower supply for the western United States. While the policies currently regulating the basin will expire by 2026, planning remains challenging due to intertwined climate variability and policy uncertainties. Based on streamflow projections from 10 dynamically downscaled CMIP6 global climate models and unique methods that add and remove internal variability, we evaluate future conditions at Powell and Mead under existing and alternative policies. Due to projected streamflow declines, under existing policy, both reservoirs will face substantial risks (>80% likelihood) of reaching dead pool before 2060. Adopting recently proposed alternative policies reduces but doesn’t eliminate such risks. All policies also exhibit tipping points where reservoir levels can change rapidly with a slight change in streamflow. A sustainable policy may require larger reductions to further reduce the reservoirs’ dead pool risks and provide better buffers from sudden changes.

Science & Technology - Other Topics↗

NASA NTRS1993

Chlorine catalyzed destruction of ozone - Implications for ozone variability in the upper stratosphere

The annual mean and the annual amplitude of ozone have been derived from ozone measurements from the SBUV and SBUV/2 spectrometers on board the Nimbus-7 and NOAA-11 satellites. These values differ significantly from values calculated using a two-dimensional model of stratospheric photochemistry and dynamics with standard chemistry. We have found that the differences between the calculated and data-derived values are considerably improved by changing the partitioning in the Cly family to create a larger reservoir of HCl and reducing ClO. This is accomplished by including a channel for the products HCl + O2 from the reaction ClO + OH in addition to the products Cl + HO2. This partitioning also improves the agreement between the calculated and measured values of ClO/HCl ratio.

Chandra, Sushil↗

DOE OSTI2024

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Xanthos-Lake Model Source Code

This repository contains the source code for Xanthos-Lake, a lake-modeling extension of the Xanthos framework that introduces a coupled lake component comprising the Xanthos-Lake Snow and Ice Model (xLSIM) and the Xanthos-Lake Water Balance Model (xLWBM). xLSIM is a basin-aware machine-learning model for lake snow, ice, and thermal conditions. It predicts monthly lake ice thickness, snow depth, snow-cover fraction, mixing-layer temperature, and lake ice fraction from meteorological forcing and lake surface-area information. It uses sequence-based deep-learning architectures, including Transformer and hybrid Long Short-Term Memory–Transformer (LSTM–Transformer) models, together with seasonal encoding, multi-lake learning, physical masking, and basin-level cryospheric and non-cryospheric classification. The training workflow uses Ray for scalable execution and includes optional Ray Tune hyperparameter optimization. Model predictions, observations, diagnostics, and feature-importance outputs are written in NetCDF. xLWBM is the water-balance component of the new lake framework. It simulates monthly lake storage, surface area, evaporation, inflow, outflow, and lake–groundwater exchange. It combines physical water-balance equations with calibrated bathymetric relationships, weir-based outlet flow, modified Penman open-water evaporation, groundwater head relaxation, Penman–Monteith snow and ice sublimation, and snow, ice, and thermal conditions supplied by xLSIM. The model calibrates lake parameters against satellite-derived surface-area data, using evaporation-based calibration where surface-area data are unavailable, and supports small, medium, and large lake classes. For large lakes, xLWBM is integrated with the managed-routing workflow so that lake storage and outflow interact directly with downstream river routing and reservoir operations. Together, xLSIM and xLWBM provide Xanthos with a coupled lake-modeling capability. xLSIM supplies the snow, ice, and thermal conditions that affect lake evaporation and snow- and ice-related water exchanges, while xLWBM translates those conditions into dynamic lake storage, surface area, evaporation, and discharge. In return, xLWBM supplies evolving lake surface area to xLSIM. This coupling enables Xanthos to represent lakes as active hydrologic components within basin-scale water-availability and routing simulations.

Machine Learning↗

DOE OSTI2025

COMPUTATION FLUID DYNAMICS ANALYSIS FOR GENERIC SMALL MODULAR REACTOR CONTAINMENT SEPARATE EFFECTS TEST

It is desirable for fourth-generation Small Modular Reactors to be passively cooled in standard and accident operations. Passive Containment Cooling Systems can reject heat from the containment structure, without using pumps or blowers. The targeted design containment structure is a large, domed, stainless steel, cylindrical vessel. In a postulated Design Basis Accident, steam will flash inside containment. Steam condensation occurs on the inner containment wall and transfers heat through the steel containment into a large body of water known as the annular reservoir (AR) surrounding the vessel serving as the ultimate heat sink. Natural circulation drives the flow in the AR and heat will be released to the environment by evaporation of water. Unique containment geometry requires a separate effects test (SET) facility for the verification and validation of the computer code and evaluation model development and assessment for reactor licensing efforts. In this study, STAR-CCM+, a computational fluid dynamics (CFD) code was used to inform the decision-making process on the design of the SET. The CFD simulation modeled, a two-phase turbulent flow with fluid film development and heat transfer for different containment geometries. The Reactor Excursion and Leak Analysis Program will also be used in a code-to-code verification against the CFD results.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

NASA NTRS1994

Hubble Space Telescope observations of planetary nebulae in the magellanic clouds. 2: SMP 85, a young planetary

We have obtained Hubble Space Telescope Planetary Camera images in the (O III) lambda 5007 A emission line, and Faint Object Spectrograph (FOS) UV spectrophotometry of the low-excitation planetary nebula SMP 85 in the Large Magellanic Cloud. By combining these results with existing optical spectrophotometry, absolute flux measurements, and dynamical and density information, we have been able to construct a fully self-consistent nebular model. This proves that SMP 85 is a dense, young, carbon-rich object which started to be ionized about 500-1000 years ago, and which contains a substantial inner reservoir of atomic or molecular gas, probably in the form of many small cloudlets. These cloudlets have been ejected at a velocity not exceeding 6 km/s, a result which, together with the morphology is an important clue to mass loss during late asymptotic giant branch (AGB) evolution. We have directly detected the central star through its UV continuum emission, and from both Zanstra techniques and nebular modeling derive a stellar temperature of 46000 +/- 2000 K, a stellar luminosity of 7300 +/- 700 solar luminosity, and a core mass of 0.63-0.67 solar mass. The nebular analysis also demonstrates that there is severe depletion of the nebular gases onto dust grains, most likely of the calcium magnesium silicate variety; a surprising result in view of the carbon-rich nature of the ionized nebula.

Dopita, Michael A.↗

NASA NTRS2004

High Silicate Crystalline-to-Amorphous Ratios in Comets C/2001 Q4 (NEAT) and Hale-Bopp

Crystalline silicates, by their apparent absence in the ISM, are dust grains that experienced high temperatures in the solar nebula. Mg-rich crystalline silicates formed either by condensation from hot nebular gases (1450 K) or by the annealing of Mg-rich amorphous silicates (approximately 1000 K) in shocks in the 5-10AU region or by radial transport into and out of the hot inner zones, e.g., T(sub d) greater than 1000K at r(sub h) less than 5AU, 10(exp -6) -10(exp -5) solar mass per year, alpha = 10(exp -4) of the early solar nebula. Mg-rich crystalline silicates are found in interplanetary dust particles (IDPs) and produce IR spectral features in many Oort cloud comets. In May 2004, we discovered strong crystalline silicate features in the dynamically new Oort cloud comet C/2001 Q4 (NEAT). Thermal emission modeling of comets Q4 and C/1995 O1 (Hale-Bopp) demonstrate that both these comets have similar, high silicate crystalline-toamorphous ratios of 2.4 and 2.1, respectively, indicating that these icy planetesimals aggregated from similar reservoirs of material or that crystalline silicates were widely distributed within the comet-forming zone. This argues for efficient annealing mechanisms and radial mixing.

Wooden, D. H.↗

NASA NTRS2004

High Silicate Crystalline-to-Amorphous Ratios in Comets C/2001 Q4 (NEAT) and Hale-Bopp

Crystalline silicates, by their apparent absence in the ISM, are dust grains that experienced high temperatures in the solar nebula. Mg-rich crystalline silicates formed either by condensation from hot nebular gases (1450 K) or by the annealing of Mg-rich amorphous silicates (approx. 1000 K) in shocks in the 5-10 AU region or by radial transport into and out of the hot inner zones, e.g., T(sub d) > 1000 K at r(sub h) < 5 AU, 10(exp -6) - 10(exp -5) M(sub O)/yr, alpha = 10(exp -4) of the early solar nebula. Mg-rich crystalline silicates are found in interplanetary dust particles (IDPs) and produce IR spectral features in many Oort cloud comets. In May 2004, we discovered strong crystalline silicate features in the dynamically new Oort cloud comet C/2001 Q4 (NEAT). Thermal emission modeling of comets Q4 and C/1995 O1 (Hale-Bopp) demonstrate that both these comets have similar, high silicate crystalline-to-amorphous ratios of 2.4 and 2.1, respectively, indicating that these icy planetesimals aggregated from similar reservoirs of material or that crystalline silicates were widely distributed within the comet-forming zone. This argues for efficient annealing mechanisms and radial mixing.

Wooden, D. H.↗

DOE OSTI2025

Exploring Water System Vulnerabilities in California's Central Valley Under the Late Renaissance Megadrought and Climate Change

Abstract California faces cycles of drought and flooding that are projected to intensify, but these extremes may impact water users across the state differently due to the region's natural hydroclimate variability and complex institutional framework governing water deliveries. To assess these risks, this study introduces a novel exploratory modeling framework informed by paleo and climate‐change based scenarios to better understand how impacts propagate through the Central Valley's complex water system. A stochastic weather generator, conditioned on tree‐ring data, produces a large ensemble of daily weather sequences conditioned on drought and flood conditions under the Late Renaissance Megadrought period (1550–1580 CE). Regional climate changes are applied to this weather data and drive hydrologic projections for the Sacramento, San Joaquin, and Tulare Basins. The resulting streamflow ensembles are used in an exploratory stress test using the California Food‐Energy‐Water System model, a highly resolved, daily model of water storage and conveyance throughout California's Central Valley. Results show that megadrought conditions lead to unprecedented reductions in inflows and storage at major California reservoirs. Both junior and senior water rights holders experience multi‐year periods of curtailed water deliveries and complete drawdowns of groundwater assets. When megadrought dynamics are combined with climate change, risks for unprecedented depletion of reservoir storage and sustained curtailment of water deliveries across multiple years increase. Asymmetries in risk emerge depending on water source, rights, and access to groundwater banks.

Gupta, Rohini S. [School of Civil and Environmenta↗

NASA NTRS2015

Pyrotechnically Actuated Gas Generator Utilizing Aqueous Methanol

A gas-generating device was developed to supplement the ram-air inflation of a supersonic ballute. The device is designed to initially pressurize the ballute following deployment, exposing and orienting its ram-air inlets to free-stream air for complete inflation. The supplemental pressurization decreases the total inflation time, and increases the likelihood of a successful inflation. The device contains a reservoir filled with an aqueous mixture of methanol that, when released in to the interior of the ballute, rapidly vaporizes due to the low ambient pressure. Upon activation of the device, a pair of redundant ring mechanisms initiate pyrotechnic charges that pressurize and rupture the reservoir, resulting in ejection of the methanol in to the ballute. In addition to its role in inflation, the device serves as the structural connection to the ballute. Analytical models were developed for the inflation capability of the device, which were verified using vacuum chamber testing of developmental hardware. Static, deployment, and environmental testing demonstrated the functionality of the ring mechanism and reservoir under several temperature and pressure conditions. Finally, the device was successfully operated during the first Supersonic Flight Dynamics Test (SFDT) of NASA's Low Density Supersonic Decelerator (LDSD) project. The design architecture is scalable to accommodate different quantities of gas generation, can be adjusted to operate in a variety of temperature and atmospheric pressure regimes, and provides a robust device that may be installed with minimal risk to personnel or hardware.

IA↗

DOE OSTI2026

Dynamically reconfigurable topological routing in nonlinear photonic systems

The propagation path of topologically protected states is bound to the interface between regions with different topology, and as such, the functionality of linear photonic devices leveraging these states is fixed during fabrication. Here, we propose a mechanism for dynamic control over a driven dissipative system’s local topology, yielding reconfigurable topological interfaces and thus tunable paths for protected routing. We illustrate our approach in non-resonantly pumped polariton lattices, where the nonlinear interaction between the polaritons and the exciton reservoir due to non-resonant pumping can yield picosecond-scale changes in the propagation paths of the chiral edge states. To analytically confirm the numerically observed topological dynamics, we generalize the spectral localizer framework to non-linear non-Hermitian Chern materials and apply this framework to a continuous model of the polariton system based on a driven-dissipative Gross-Pitaevskii equation. In doing so, we show that the local changes in the polariton lattice’s topology are captured by a local Chern marker. Looking forward, we anticipate such reconfigurable topological routing will enable the realization of novel classes of topological photonic devices.

Wong, Stephan [Sandia National Laboratories (SNL-N↗

DOE OSTI2026

Analytic rate theory of polariton relaxation that explains long polariton lifetime

Hybridization of a molecular exciton with a quantized photon creates a polariton. Despite extensive experimental investigations, the apparent lifetime of the exciton–polariton is not well-understood. For this work, we examined the steady-state population dynamics for a Holstein–Tavis–Cumming Hamiltonian to illuminate the long-term polaritonic dynamics and lifetime of the exciton–polariton in an optical cavity. For a realistic description of polariton relaxation, cavity loss and various exciton decay channels are included in the model. We found that in the presence of weak but finite exciton loss, the apparent lifetime of the lower polariton coincides with the out-of-cavity exciton lifetime and is independent of cavity-matter detuning. This is a simple explanation for the experimentally observed lifetimes for exciton polaritons and theoretically justifies the dark state reservoir hypothesis. Furthermore, if the upper polariton is initially populated, the system reaches the steady state very quickly, leading to single-exponential polariton relaxation. Starting from the lower polariton leads to a longer pre-steady-state time period, leading to double-exponential relaxation. Finally, we considered the effect of site orientational disorders and the exciton frequency disorderers. Under the collective limit, the effects of this disorder can be included in Fermi’s golden rule population dynamics without explicit sampling. For the exciton energy disorders, numerical calculations are needed. Our theoretical framework is applicable to interpret exciton–polariton experiments, especially related to the measured apparent lifetime of polaritons.

Chemical dynamics↗

DOE OSTI2025

Echo state network for coarsening dynamics of charge density waves

An echo state network (ESN) is a type of reservoir computer that uses a recurrent neural network with a sparsely connected hidden layer. Compared with other recurrent neural networks, one great advantage of ESN is the simplicity of its training process. Yet, despite the seemingly restricted learnable parameters, ESN has been shown to successfully capture the spatial-temporal dynamics of complex patterns. Here we build an ESN to model the coarsening dynamics of charge-density waves (CDWs) in a semiclassical Holstein model, which exhibits a checkerboard electron density modulation at half-filling stabilized by a commensurate lattice distortion. The inputs to the ESN are local CDW order parameters in a finite neighborhood centered around a given site, while the output is the predicted CDW order of the center site at the next time step. Special care is taken in the design of couplings between hidden layer and input nodes to ensure lattice symmetries are properly incorporated into the ESN model. Since the model predictions depend only on CDW configurations of a finite domain, the ESN is scalable and transferrable in the sense that a model trained on dataset from a small system can be directly applied to dynamical simulations on larger lattices. Furthermore, our work opens avenues for efficient dynamical modeling of pattern formations in functional electron materials.

2-dimensional systems↗

DOE OSTI2025

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↗