Search NASA⌕ Search

SEARCH · Search NASA

Results for “Enthalpy of formation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Desorption lifetimes and activation energies influencing gas–surface interactions and multiphase chemical kinetics

Abstract. Adsorption and desorption of gases on liquid or solid substrates are involved in multiphase processes and heterogeneous chemical reactions. The desorption energy (Edes0), which depends on the intermolecular forces between adsorbate and substrate, determines the residence time of chemical species at interfaces. We show how Edes0 and temperature influence the net uptake or release of gas species, the rates of surface–bulk exchange and surface or bulk reactions, and the equilibration timescales of gas–particle partitioning. Using literature data, we derive a parameterization to estimate Edes0 for a wide range of chemical species based on the molecular mass, polarizability, and oxygen-to-carbon ratio of the desorbing species independent of substrate-specific properties, which is possible because of the dominant role of the desorbing species' properties. Correlations between Edes0 and the enthalpies of vaporization and solvation are rooted in molecular interactions. The relation between Edes0 and desorption kinetics reflects the key role of interfacial exchange in multiphase processes. For small molecules and semi-volatile organics (VOC, IVOC, SVOC), Edes0 values around 10–100 kJ mol−1 correspond to desorption lifetimes around nanoseconds to days at room temperature. Even higher values up to years are obtained at low temperatures and for low volatile organic compounds (LVOC, ELVOC/ULVOC) relevant for secondary organic aerosols (SOA). Implications are discussed for SOA formation, gas–particle partitioning, organic phase changes, and indoor surface chemistry. We expect these insights to advance the mechanistic and kinetic understanding of multiphase processes in atmospheric and environmental physical chemistry, aerosol science, materials science, and chemical engineering.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Applications of Low Density Flow Techniques and Catalytic Recombination at the Johnson Space Center

The talk presents a brief background on defInitions of catalysis and effects associated with chemically nonequilibrium and low-density flows of aerospace interest. Applications of catalytic recombination on surfaces in dissociated flow are given, including aero heating on reentry spacecraft thermal protection surfaces and reflection of plume flow on pressure distributions associated with the space station. Examples include aero heating predictions for the X-38 test vehicle, the inlet of a proposed gas-sampling probe used in high enthalpy test facilities, and a parabolic body at angle of attack. The effect of accommodation coefficients on thruster induced pressure distributions is also included. Examples of tools used include simple aero heating formulas based on boundary layer solutions, an engineering approximation that uses axisymmetric viscous shock layer flow to simulate full three dimensional flow, full computational fluid dynamics, and direct simulation Monte-Carlo calculations. Methods of determining catalytic recombination rates in arc jet flow are discus ed. An area of catalysis not fully understood is the formation of single-wall carbon nanotubes (SWNT) with gas phase or nano-size metal particles. The Johnson Space Center is making SWNTs using both a laser ablation technique and an electric arc vaporization technique.

Scott, Carl D.↗
DOE OSTI2024

Unifying thermochemistry concepts in computational heterogeneous catalysis

Thermophysical properties of adsorbates and gas-phase species define the free energy landscape of heterogeneously catalyzed processes and are pivotal for an atomistic understanding of the catalyst performance. These thermophysical properties, such as the free energy or the enthalpy, are typically derived from density functional theory (DFT) calculations. Enthalpies are species-interdependent properties that are only meaningful when referenced to other species. The widespread use of DFT has led to a proliferation of new energetic data in the literature and databases. However, there is a lack of consistency in how DFT data is referenced and how the associated enthalpies or free energies are stored and reported, leading to challenges in reproducing or utilizing the results of prior work. Additionally, DFT suffers from exchange–correlation errors that often require corrections to align the data with other global thermochemical networks, which are not always clearly documented or explained. In this review, we introduce a set of consistent terminology and definitions, review existing approaches, and unify the techniques using the framework of linear algebra. This set of terminology and tools facilitates the correction and alignment of energies between different data formats and sources, promoting the sharing and reuse of ab initio data. Standardization of thermochemistry concepts in computational heterogeneous catalysis reduces computational cost and enhances fundamental understanding of catalytic processes, which will accelerate the computational design of optimally performing catalysts.

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

Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

Thermodynamics of Silicon-Hydroxide Formation in H2O Containing Atmospheres

The formation of volatile silicon-hydroxide species from SiO2 in water containing atmospheres has been identified as a potentially important mode of degradation of Si-based ceramics. Availability of thermodynamic data for these species is a major problem. This study is part of an ongoing effort to obtain reliable, experimentally determined thermodynamic data for these species. The transpiration method was used to measure the pressure of Si-containing vapor in equilibrium with SiO2 (cristobalite) and Ar + H2O(g) with various mole fractions of water vapor, X(sub H2O), at temperatures ranging from 1000 to 1780 K. Enthalpies and entropies for the reaction, SiO2(s) + 2H2O(g) = Si(OH)4(g), were obtained, at X(sub H2O) = 0.15 and 0.37, from the variation of lnK with 1/T according to the 'second law method'. The following data were obtained: delta(H)deg = 52.9 +/- 3.7 kJ/mole and delta(S)deg = -68.6 +/- 2.5 J/mole K at an average temperature of 1550 K, and delta(H)deg = 52.5+/-2.0 kJ/mole and delta(S)deg= -69.7 +/- 1.5 J/moleK at an average temperature of 1384 K, for X(sub H2O)= 0.15 and 0.37, respectively. These data agree with results from the literature obtained at an average temperature of 1600 K, and strongly suggest Si(OH)4(g) is the dominant vapor species. Contradictory results were obtained with the determination of the dependence of Si-containing vapor pressure on the partial pressure of water vapor at 1187 and 1722 K. These results suggested the Si-containing vapor could be a mixture of Si(OH)4 + SiO(OH)2. Further pressure dependent studies are in progress to resolve these issues.

Copland, Evan↗
DOE OSTI2020

Glass transition temperature studies of planetary ball milled glasses: Accessing the rapidly cooled glassy state in Na 4 P 2 S 7-x O x , 0 ≤ x ≤ 7, Oxy-thio phosphate glasses

This paper details the very first in-depth study of the glass transition temperature (T g ) for planetary balled milled (PBM) Na 4 P 2 S 7-x O x (NaPSO), where 0 ≤ x ≤ 7, glasses. For 0 ≤ x ≤ 5, fully amorphous, homogeneous, and chemically reacted compositions can be prepared. For samples with 5 < x ≤ 7, partially amorphous but fully chemically reacted glass-ceramics can be produced. The reproducible, thermally cyclable onset T g and the onset crystallization temperature, T c , were investigated as a function of the composition and of the short-range order (SRO) structures of the glasses. Additionally, we examined the extreme sub-T g exothermic relaxations present due to the high energy milling processing and correspondingly high effective quenching rate. Using a differential scanning calorimeter (DSC) and scanning from well below the onset T g , the glasses are observed to relax exothermally below T g . The integrated relaxation enthalpy, ΔH rel , is ~ 0 for the x = 2 glass, but is significantly exothermic for lower and higher values of x. Once the high fictive state of these glasses has been erased by scanning to above the T g , but below the T c of the glass, the DSC scans of the glasses recover normal behavior for a normally cooled glass without any relaxation exotherm. The T g s of these glasses remain relatively constant for 0 ≤ x ≤ 3, but for glasses with 3 < x ≤ 5, T g increases sharply following the combined effects of the phosphate chain length increasing caused by the formation of bridging oxygens and the increasing concentration of more charge dense non-bridging oxygens (NBOs). For x > 5, both T g and Tc decrease sharply, becoming nearly identical at x = 7 where glass formation essentially ceases. For smaller values of x, T c increases with both chain length and the structural complexity of the composition but decreases sharply for glass-ceramics of x > 5.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Low-temperature formation of Ti2AlN during post-deposition annealing of reactive multilayer systems

M n+1 AXn-phase Ti 2 AlN thin-films were synthesized using reactive sputtering-based methods involving the deposition of single-layer TiAlN, and Ti/AlN and TiN/TiAl multilayers of various modulation periods at ambient temperature and subsequent annealing at elevated temperatures. Ex situ and in situ x-ray diffraction measurements were used to characterize the Ti 2 AlN formation temperature and phase fraction. During annealing, Ti/AlN multilayers yielded Ti 2 AlN at a significantly lower in situ temperature of 650 °C compared to TiN/TiAl multilayers or single-layer TiAlN (750 °C). The results suggest a reactive multilayer mechanism whereby distinct Ti and AlN layers react readily to release exothermic energy resulting in lower phase transition temperatures compared to TiN and TiAl layers or mixed TiAlN. With a modulation period of 5 nm, however, Ti/AlN multilayers yielded Ti 2 AlN at a higher temperature of 750 °C, indicating a disruption of the reactive multilayer mechanism due to a higher fraction of low-enthalpy interfacial TiAlN within the film.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Entropy is an important design principle in the photosystem II supercomplex

Photosystem II (PSII) can achieve near-unity quantum efficiency of light harvesting in ideal conditions and can dissipate excess light energy as heat to prevent the formation of reactive oxygen species (ROS) under light stress. Understanding how this pigment–protein complex accomplishes these opposing goals is a topic of great interest that has so far been explored primarily through the lens of the system energetics. Despite PSII’s known flat energy landscape, a thorough consideration of the entropic effects on energy transfer in PSII is lacking. In this work, we aim to discern the free energetic design principles underlying the PSII energy transfer network. To accomplish this goal, we employ a structure-based rate matrix and compute the free energy terms in time following a specific initial excitation to discern how entropy and enthalpy drive ensemble system dynamics. We find that the interplay between the entropy and enthalpy components differ among each protein subunit, which allows each subunit to fulfill a unique role in the energy transfer network. This individuality ensures that PSII can accomplish efficient energy trapping in the reaction center (RC), effective nonphotochemical quenching (NPQ) in the periphery, and robust energy trapping in the other-monomer RC if the same-monomer RC is closed. We also show that entropy, in particular, is a dynamically tunable feature of the PSII free energy landscape accomplished through regulation of LHCII binding. These findings help rationalize natural photosynthesis and provide design principles for more efficient solar energy harvesting technologies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Influence of Carboxymethyl Cellulose as a Thickening Agent for Glauber’s Salt-Based Low Temperature PCM

This work is focused on a novel, promising low temperature phase change material (PCM), based on the eutectic Glauber’s salt composition. To allow phase transition within the refrigeration range of temperatures of +5 °C to +12 °C, combined with a high repeatability of melting–freezing processes, and minimized subcooling, the application of three variants of sodium carboxymethyl cellulose (Na-CMC) with distinct molecular weights (700,000, 250,000, and 90,000) is considered. The primary objective is to optimize the stabilization of this eutectic PCM formulation, while maintaining the desired enthalpy level. Preparation methods are refined to ensure repeatability in mixing components, thereby optimizing performance and stability. Additionally, the influence of Na-CMC molecular weight on stabilization is examined through differential scanning calorimetry (DSC), T-history, and rheology tests. The PCM formulation of interest builds upon prior research in which borax, ammonium chloride, and potassium chloride were used as additives to sodium sulfate decahydrate (Glauber’s salt), prioritizing environmentally responsible materials. The results reveal that CMC with molecular weights of 250 kg/mol and 90 kg/mol effectively stabilize the PCM without phase separation issues, slowing crystallization kinetics. Conversely, CMC of 700 kg/mol proved ineffective due to the disruption of gel formation at its low gel point, hindering higher concentrations. Calculations of ionic concentration indicate higher Na ion content in PCM stabilized with 90 kg/mol CMC, suggesting increased ionic interactions and gel strength. A tradeoff is discovered between the faster crystallization in lower molecular weight CMC and the higher concentration required, which increases the amount of inert material that does not participate in the phase transition. After thermal cycling, the best formulation had a latent heat of 130 J/g with no supercooling, demonstrating excellent performance. This work advances PCM’s reliability as a thermal energy storage solution for diverse applications and highlights the complex relationship between Na-CMC molecular weight and PCM stabilization.

25 ENERGY STORAGE↗
DOE OSTI2020

Thermal Activation of a Copper-Loaded Covalent Organic Framework for Near-Ambient Temperature Hydrogen Storage and Delivery

Copper(II) formate is efficiently incorporated into the pores of a 2D imine-based covalent organic framework (COF) via coordination with the phenol and imine groups. The coordinated metal ion is then reduced to Cu(I) with a thermal treatment that evolves CO 2 . After loading with hydrogen gas, the majority of H 2 desorbs from the coordinatively saturated Cu(II) COF at temperatures < -100 degrees C. However, the activated Cu(I) COF retains adsorbed H 2 above room temperature. Adsorption/desorption of H 2 was highly reversible. Diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) strongly supports a molecular hydrogen interaction with Cu(I). A Kissinger analysis of variable ramp rate desorption experiments estimates the enthalpy of H 2 desorption from Cu(I) at 15 kJ mol -1 . The results represent an advance toward practical H 2 storage and delivery in a lightweight, stable, and highly versatile material.

08 HYDROGEN↗
DOE OSTI2026

Technoeconomic analysis of hydrogen storage using 1,4-butanediol (BDO)/γ-butyrolactone (GBL) as a stationary backup power system

Liquid organic hydrogen carriers (LOHCs) are compounds that store and release hydrogen in stable forms at high density. While one-way carriers such as methanol and ammonia have gained attention, their economic advantages are often realized from their use as an export product and direct use as a fuel. Liquid carrier materials that can instead be cycled for energy storage have promise for stationary power applications. In particular, the reversible LOHC system 1,4-butanediol (BDO, H 2 -rich) and gamma-butyrolactone (GBL, H 2 -lean) has a lower enthalpy of dehydrogenation compared to conventional cyclic hydrocarbons and can utilize non-precious metal copper-based catalysts. Here, in this study, BDO/GBL system capital and operating expenses are characterized for vapor phase versus liquid phase hydrogenation and dehydrogenation in a 10 MW backup power application corresponding to sizing of Tier 2 datacenters as well as other critical infrastructure such as hospitals. Costs are benchmarked against two incumbent technologies: a well-established methylcyclohexane/toluene carrier system and compressed gas storage. BDO/GBL storage costs are found to differ substantially between operating modes, with liquid phase hydrogenation coupled with liquid phase dehydrogenation leading to the lowest LCOS of $\$$4.58/kg H 2 in the absence of byproduct formation. In this bounding case, LCOS for the BDO/GBL system is lower than for MCH/TOL ($\$$6.97/kg H 2 ) and compressed gas ($\$$8.48/kg H 2 at 170 bara and $\$$12.05/kg H 2 at 350 bara). However, escalating costs of carrier replacement due to byproduct formation (ranging from an added $\$$6–11/kg H 2 ) illustrate the need for highly selective catalysts to ensure BDO/GBL carrier viability.

BDO/GBL↗
DOE OSTI2024

Defect-Limited Mobility and Defect Thermochemistry in Mixed A-Cation Tin Perovskites: (CH3NH3)1-xCsxSnBr3

Hybrid organic-inorganic semiconductors crystallizing in the perovskite structure present a significant opportunity for realizing defect-tolerant semiconductors. In this work, we examine the solid solution, (CH3NH3)1-xCsxSnBr3, and identify the thermochemistry dictating the intrinsic carrier concentrations and how local structural distortions influences this electronic behavior. This family of compounds exhibits the expected systematic trend in decreasing optical gap with the cesium to methylammonium A-site mixing ratio, x, in the visible region. However, the carrier mobility, as determined from time-resolved microwave conductivity measurements, trends opposite to that expected from first-principles calculations combined with Boltzmann scattering theory calculations of the carrier mobility. We propose that this is a result of increasing carrier scattering with x in (CH3NH3)1-xCsxSnBr3 due to a significant increase in the carrier density with x. By examining the dependence of the carrier density as a function of x, we infer the compositional-dependence of the average enthalpy and nonconfigurational entropy per defect. While diffraction reveals a cubic aristotypic perovskite structure as a function of x at room temperature, the pair distribution functions obtained from synchrotron X-ray total scattering from these materials are better described by symmetry-adapted displacement modes of the Pm3m crystal structure, which we attribute to a large degree of anharmonic dynamics of the atom positions. This analysis shows that CH3NH3+-rich compositions retain mostly linear Sn-Br-Sn bonding environments through the displacements, while Cs+-rich compositions lead to more significantly bent Sn-Br-Sn environments. We propose that these bent bonding arrangements in Cs-rich compositions yield a higher propensity for defect formation. This also provides a rationale for carrier trapping that gives rise to anomalous microwave transients. Together, these results provide insight into the structure-dynamics-properties relationships in this highly anharmonic system with high amplitude atomic motions and low defect formation energies.

defect thermochemistry↗
DOE OSTI2025

A critical analysis of U-Pu-Zr phase transitions using calorimetric, microstructural, and phase equilibria data

Metallic fuels consisting primarily of uranium, plutonium, and zirconium (U-Pu-Zr) are a leading material candidate for fast-spectrum nuclear reactors. Early demonstration programs proved the principle of safe and efficient fast reactor operation, however there is still considerable uncertainty regarding the phase equilibria and microstructural evolution across the ternary composition space. Quantitative phase formation and identification measurements are scarce and often incomplete, with studies reporting either phase transition temperatures or phase identification data, but not both from the same specimens. In this study, we critically compared experimental and calculated phase transition data and correlated with the microstructure and phase characterization data of as-cast and annealed U-Pu-Zr alloys. Differential scanning calorimetry (DSC) was used to measure phase transitions in the subsolidus regions (723−948 K) of three ternary U-Pu-Zr alloys with similar plutonium concentrations but various U/Zr ratios. Due to sluggish kinetics and narrow ranges of phase stability, complex peaks required the use of a Frazier-Suzuki peak fitting algorithm to deconvolute and calculate transition peak temperatures and enthalpies. We also identified trends of phase transition behavior by critically comparing our DSC data with previous phase transition measurements as well as historical and calculated phase equilibrium diagrams. In conclusion, this provides a critical approach for benchmarking and assessing the quality of new U-Pu-Zr phase equilibria data prior to its incorporation into nuclear material databases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Entropic Control of HD Exchange Rates over Dilute Pd-in-Au Alloy Nanoparticle Catalysts

Dilute Pd-in-Au alloy catalysts are promising materials for selective hydrogenation catalysis. Extensive surface science studies have contributed mechanistic insight on the energetic aspect of hydrogen dissociation, migration and recombination on dilute alloy systems. Yet, translating these fundamental concepts to the kinetics and free energy of hydrogen dissociation on nanoparticle catalysts operating at ambient pressures and temperatures remains challenging. Here, the effect of the Pd concentration and Pd ensemble size on the catalytic activity, apparent activation energy and rate limiting process is addressed by combining experiment and theory. Experiments in a flow reactor show that a compositional change from 4 to 8 atm% Pd of the Pd-in-Au alloy catalyst leads to strong increase in activity, even exceeding the activity per Pd atom of monometallic Pd under the same conditions, albeit with an increase in apparent activation energy. First-principles calculations show that the rate and apparent activation enthalpy for HD exchange increase when increasing the Pd ensemble size from single Pd atoms to Pd trimers in a Au surface, suggesting that the ensemble size distribution shifts from mainly single Pd atoms on the 4 atm% Pd alloy to larger Pd ensembles of at least three atoms for the 8 atm% Pd/Au catalyst. The DFT studies also indicated that the rate-controlling process is different: H 2 (D 2 ) dissociation determines the rate for single atoms whereas recombination of adsorbed H and D determines the rate on Pd trimers. Both experiment and theory suggest that the increased reaction rate with increasing Pd content and ensemble size stems from an entropic driving force. Lastly, our results support hydrogen migration between Pd sites via Au and indicate that the dilute alloy design prevents the formation of subsurface hydrogen, which is crucial in achieving high selectivity in hydrogenation catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

Effects of Wall Cooling on Supersonic Modes in High-Enthalpy Hypersonic Boundary Layers over a Cone

High-enthalpy hypersonic boundary-layer transition plays an important role in many entry/descent vehicles. Aerothermodynamic performance of these vehicles strongly depends on the transition location on the surface. However, detailed transition flow physics in these chemically reacting boundary layers are poorly understood and transition estimates during the design phase rely heavily on empirically-derived transition criteria such as 𝑹𝒆(𝜽)/𝑴(𝒆). One of the most intriguing characteristics in hypersonic boundary layers is the presence of unstable supersonic modes, first identified in the 1990s. Due to a recent surge in hypersonic applications, there has been renewed interest in studying the flow physics of supersonic modes using either theory or direct numerical simulations. This paper investigates the rise of supersonic modes in a high-enthalpy hypersonic flow over a 𝟓° half-angle cone at various wall temperatures using both quasiparallel linear stability theory and linear PSE. It was found that supersonic modes exist in all wall temperature conditions including the adiabatic wall case. Cooler wall temperature causes the second Mack mode to become an unstable supersonic mode naturally downstream of the upper-branch neutral location when the wall is sufficiently cooled. In terms of the integrated growth, the second mode is still the dominant mode. Nonetheless, supersonic modes can cause additional series of relatively weaker growth than that of the second mode beyond the peak amplitude location. According to the present linear PSE results, contrary to what was speculated in the literature, supersonic-mode pressure disturbance structure radiated into the freestream is nonacoustic in nature and the formation of unstable supersonic modes is mainly associated with the synchronization of phase speed between the instability and acoustic waves in a nonp ayer, not due to nonlinear modal interaction as suggested in the literature.

Chang, Chau-Lyan↗
DOE OSTI2020

Bond Synergy Model for Bond Energies in Alloy Oxides

In this work we introduce a metal-oxide bond-energy model for alloy oxides based on pure-phase bond energies and bond synergy factors that describe the effect of alloying on the bond energy between cations and oxygen, an important quantity to understand the formation of alloy oxides and their composition. This model is parameterized for binary cation-alloy oxides using density-functional theory energies and is shown to be directly transferable to multi-component alloy oxides. We parameterized the model for alloy oxide energies with metal cations that form the basis of corrosion resistant alloys, including Fe, Ni, Cr, Mo, Mn, W, Co, and Ru. We find that isoelectronic solutes allow quantification of pure-phase bond energies in oxides and the calculated bond energy values give sensible results compared to common experience, including the role of Cr as the passive-layer former in Fe–Ni–Cr alloys for corrosion applications. Additionally, the bond synergy factors give insights into the mutual strengthening and weakening effects of alloying on cation-oxygen bonds and can be related to enthalpy of mixing and charge neutrality constraints. We demonstrate how charge neutrality can be identified and achieved by the oxidation states that the different cations assume depending on alloy composition and the presence of defects.

Electrochemistry↗
DOE OSTI2023

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗