Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir simulation”

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 721 records · Page 40

Possible impacts of molten salt reactors on the International Monitoring System

Molten salt reactors (MSRs) are gaining support as many countries look for ways to increase power generation and replace aging nuclear energy production facilities. MSRs have inherently safe designs, are scalable in size, can burn transuranic wastes from traditional solid fuel nuclear reactors, can store excess heat in thermal reservoirs for water desalination, and can be used to produce medical isotopes as part of the real-time liquid-fuel recycling process. The ability to remove 135 Xe in real time from the fuel improves the power production in an MSR because 135 Xe is the most significant neutron-absorbing isotope generated by nuclear fission. Xenon-135, and other radioactive gases, are removed by sparging the fuel with an inert gas while the liquid fuel is recirculated from the reactor inner core through the heat exchangers. Without effective abatement technologies, large amounts of radioactive gas could be released during the sparging process. This work examines the potential impact of radioxenon releases on samplers used by the International Monitoring System (IMS) to detect nuclear explosions. Here, atmospheric transport simulations from seven hypothetical MSRs on different continents were used to evaluate the holdup time needed before release of radioxenon so IMS samplers would register few detections. Abatement technologies that retain radioxenon isotopes for at least 120 d before their release will be needed to mitigate the impacts from a molten salt breeder reactor used to replace a nuclear power plant. A holdup time of about 150 d is needed to reduce emissions to the average level of current nuclear power plants.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Efficient and generalizable nested Fourier-DeepONet for three-dimensional geological carbon sequestration

Geological carbon sequestration (GCS) involves injecting CO2 into subsurface geological formationsfor permanent storage. Numerical simulations could guide decisions in GCS projects by predictingCO 2 migration pathways and the pressure distribution in storage formation. However, these simula-tions are often computationally expensive due to highly coupled physics and large spatial-temporalsimulation domains. Surrogate modelling with data-driven machine learning has become a promis-ing alternative to accelerate physics-based simulations. Among these, the Fourier neural operator(FNO) has been applied to three-dimensional synthetic subsurface models. Despite its good accuracyin simulating CO 2 plume migration, it requires large computational resources in training and alsolacks generalizability. Here, to further improve performance, we have developed a nested Fourier-DeepONet by combining the expressiveness of the FNO with the modularity of a deep operatornetwork (DeepONet). This new framework is twice as efficient as a nested FNO for training and has atleast 80% lower GPU memory requirement due to its flexibility to treat temporal coordinates sepa-rately. These performance improvements are achieved without compromising prediction accuracy.In addition, the generalization and extrapolation ability of nested Fourier-DeepONet beyond thetraining range has been thoroughly evaluated. Nested Fourier-DeepONet outperformed the nestedFNO for extrapolation in time with more than 50% reduced error. It also exhibited good extrapolationaccuracy beyond the training range in terms of reservoir properties, number of wells, and injectionrate.

Lee, Jonathan E. [Department of Chemical and Envir↗

Mobilization of trace metals from caprock and formation rocks at the Illinois Basin – Decatur Project demonstration site under geological carbon dioxide sequestration conditions

One concern for geologic CO 2 sequestration is the potential leakage of CO 2 or CO 2 -saturated brine containing trace metals into overlying aquifers, which poses the risk of adversely affecting underground sources of drinking water. In this work, rock and brine samples collected in the Illinois Basin – Decatur Project (IBDP), a large-scale geological CO 2 sequestration demonstration project in Decatur, Illinois, USA, as well as synthetic brine samples were used to understand the mobilization of trace metals as a result of CO 2 –rock–brine interactions under IBDP-specific conditions (50 °C and 20.6 MPa). Rock sample characterization indicated that, at the IBDP site, trace metal concentrations were greatest in cap rock samples from the Eau Claire Formation, compared with those from formations above and below this Formation. The natural brine samples collected from the Mt. Simon Sandstone reservoir were highly reducing and saline, with trace metal concentrations up to 680 times greater than the U.S. Environmental Protection Agency-prescribed drinking water standards. Batch leaching experiments indicated that both trace metal mobilization from rock and immobilization from the brine occurred when high-pressure CO 2 was introduced into the rock-brine system. The amounts of metals mobilized from the rock generally accounted for <5% of the total metals in the rock, but for some metals, including Ni, Pb, and Tl, up to 63% of the metals in the rock were mobilized to the brine. Leaching of trace metals into synthetic brines was different from that into the natural brine. The results of this study provided current information on the trace element sources and will help in risk simulations and experimental design to further evaluate potential impact of leakage on groundwater quality, which is important for future GCS projects to consider for monitoring and containment assurance purposes.

58 GEOSCIENCES↗

A Novel Relative Permeability Model for Gas and Water Flow in Hydrate-Bearing Sediments With Laboratory and Field-Scale Application

In a producing gas hydrate reservoir the effective porosity available for fluid flow constantly changes with dissociation of gas hydrate. Therefore, accurate prediction of relative permeability using legacy models (e.g. Brooks-Corey (B-C), van Genuchten, etc.) that were developed for conventional oil and gas reservoirs would require empirical parameters to be calibrated at various S h over its range of variation, but such calibrations are precluded because of lack of experimental relative permeability data. This study proposes a new relative permeability model for gas hydrate-bearing media that is a function of maximum capillary pressure, capillary entry pressure, pore size distribution index, residual saturations, hydrate saturation, and four other constants. The three novel features of the proposed model are: (i) requires fitting its six empirical parameters only once using experimental data from any single S h , and the same set of empirical parameters predict relative permeability at all S h , (ii) includes the effect of capillarity, and (iii) includes the effect of pore-size distribution. From practical standpoint, the model can be used to simulate multiphase flow in gas hydrate-bearing sediments where the proposed relative permeability can account for the evolving hydrate saturation. The proposed model is implemented in a numerical simulator and the wall time required to perform simulations using the proposed model is shown to be similar to the time it takes to run same simulations with the B-C model. The proposed model is a step forward towards achieving the goal of physically accurate modeling of multiphase flow for gas hydrate-bearing sediments that accounts for the effect of gas hydrate saturation change on relative permeability.

54 ENVIRONMENTAL SCIENCES↗

Design and Build of Reactor Simulator for Fission Surface Power Technology Demonstrator Unit

The Nuclear Systems Team at NASA Marshall Space Flight Center (MSFC) focuses on technology development for state of the art capability in non-nuclear testing of nuclear system and Space Nuclear Power for fission reactor systems for lunar and Mars surface power generation as well as radioisotope power systems for both spacecraft and surface applications. Currently being designed and developed is a reactor simulator (RxSim) for incorporation into the Technology Demonstrator Unit (TDU) for the Fission Surface Power System (FSPS) Program, which is supported by multiple national laboratories and NASA centers. The ultimate purpose of the RxSim is to provide heated NaK to a pair of Stirling engines in the TDU. The RxSim includes many different systems, components, and instrumentation that have been developed at MSFC while working with pumped NaK systems and in partnership with the national laboratories and NASA centers. The main components of the RxSim are a core, a pump, a heat exchanger (to mimic the thermal load of the Stirling engines), and a flow meter for tests at MSFC. When tested at NASA Glenn Research Center (GRC) the heat exchanger will be replaced with a Stirling power conversion engine. Additional components include storage reservoirs, expansion volumes, overflow catch tanks, safety and support hardware, instrumentation (temperature, pressure, flow) for data collection, and power supplies. This paper will discuss the design and current build status of the RxSim for delivery to GRC in early 2012.

Godfroy, Thomas↗

Seeding Mechanism for High-Pressure Nozzles

NASA tests new launch and entry vehicle configurations in wind tunnels, where flow visualizations and quantitative flow field measurements are often desired. Some of these vehicles have rocket motors for propulsion, retro-propulsion, or reaction control. Unseeded,high-pressure air is often used to simulate these rocket motor plumes. However, it is difficult to make off-body measurements in these regions, for several reasons. First, the plumes themselves are usually unseeded and at such low pressures that the flow velocity cannot currently be measured. Second, the plumes displace seeded tunnel air, preventing measurements such as laser sheet visualizations. Third, they force shock waves ahead of the vehicle, which melts the ice crystal fog commonly used for many visualization and measurement techniques. The Virginia Tech Advanced Propulsion and Power Lab has designed and tested a novel method for seeding the flow in these small-scale, high-pressure nozzles to allow quantitative measurement with a particle-based measurement such as particle image velocimetry (PIV). The method involves a Venturi contraction to draw the seed liquid out of a reservoir and into the nozzle channel, wherein shearing forces atomize the seed into particles. The concept was tested with a laser sheet visualization, which showed the nozzle spray pattern and the controllability of the seeding flow rate. This relatively inexpensive and simple technique may prove useful in wind tunnel experiments involving particle-based laser diagnostics and small, high-pressure nozzles.

particle image velocimetry↗

Continuous integration data-driven platform of industrial-scale subsurface storage for real-time analytics

This project helped address the growing need for efficient and scalable models to support geological carbon and energy storage, which are crucial for achieving net-zero emissions. Traditionally accurate high-fidelity numerical models have been used to simulate relevant storage processes under a handful of processes, however such models are computationally demanding, making uncertainty quantification impractical. Consequently, we first developed a machine learning framework, based on Graph Neural Operators (GNOs), to improving the accuracy of model predictions for a fixed computational budget. We then developed an Ensemble of Improved Neural Operators (ENO), which uses bagging and Monte Carlo dropout techniques, to further improve prediction accuracy. Lastly, we developed the way to explain progressive transfer learning methods to reduce the amount of training data and computational cost of training (i.e., reduce trainable parameters) when using our models for multiple storage sites. Our numerical investigation, which used real-world case studies, demonstrated that our framework can significantly improve the safety and efficiency of geological storage operations, with potential applications in other domains such as geothermal reservoirs and climate modeling.

54 ENVIRONMENTAL SCIENCES↗

A Peculiar Class of Debris Disks from Herschel/DUNES: A Steep Fall Off in the Far Infrared

Context. The existence of debris disks around old main sequence stars is usually explained by continuous replenishment of small dust grains through collisions from a reservoir of larger objects. Aims. We present photometric data of debris disks around HIP 103389 (HD199260), HIP 100350 (HN Peg, HD206860), and HIP 114948 (HD 219482), obtained in the context of our Herschel Open TIme Key Program DUNES (DUst around NEarby Stars). Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 30 sigma sensitivity of a few mJy at l00 micron and 160 micron. In addition, we obtained Herschel/PACS photometric data at 70 micron for HIP 103389. These observations are complemented by a large variety of optical to far-infrared photometric data. Two different approaches are applied to reduce the Herschel data to investigate the impact of data reduction on the photometry. We fit analytical models to the available spectral energy distribution (SED) data using the fitting method of simulated therma1 annealing as well as a classical grid search method. Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths >= 70 micron. We investigate the significance of the peculiar shape of these SEDs and the impact on models of the disks provided it is real. Using grain compositions that have been applied successfully for modeling of many other debris disks, our modeling reveals that such a steep decrease of the SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a standard equilibrium collisional cascade. In contrast, a steep grain size distribution or, alternatively an upper grain size in the range of few tens of micrometers are implied. This suggests that a very distinct range of grain sizes would dominate the thermal. emission of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete and that the significance of our results depends on the version of the data reduction pipeline used. Conclusions. A new mechanism to produce the dust in the presented debris disks, deviations from the conditions required for a standard equilibrium collisional cascade (grain size exponent of -3.5), and/or significantly different dust properties would be necessary to explain the potentially steep SED shape of the three debris disks presented.

Ertel, S.↗

Performing a Launch Depressurization Test on an Inflatable Space Habitat

In July, 2014 JPL's Environmental Test Laboratory successfully performed a launch depressurization test on an inflatable space habitat proposed to be installed on the International Space Station. The inflatable habitat is to be launched in the SpaceX Dragon Trunk. During the launch, the unpressurized Dragon Trunk will rapidly change from ground level atmospheric pressure to the vacuum of space. Since the inflatable habitat is tightly folded during launch with multiple layers of bladder, Kevlar fabric sections, and micro-meteoroid shielding, it was not possible to analyze or simulate how the residual air pockets would behave during the launch. If the inflatable habitat does not vent adequately and expands, it could rupture the payload bay of the launch vehicle. A launch depressurization test was chosen as the best way to qualify the inflatable habitat. When stowed, the inflatable habitat measured approximately 241 cm (95 inches) in diameter by 152 cm (60 inches) high and weighed close to 1361 kg (3,000 pounds). Two vacuum chambers connected by a large vacuum line were used to perform this test. The inflatable habitat was mounted in the smaller chamber, which was 396 cm (13 feet) in diameter and 1128 cm (37 feet) high. The larger chamber, which was 823 cm (27 feet) in diameter and 2,591 cm (85 feet) high, was rough pumped and used as a vacuum reservoir. A two stage axial type compressor and ten Stokes vacuum pumps were also used during the depressurization. Opening a butterfly valve on the vacuum line, at the smaller chamber, was manually controlled so that the smaller chamber's depressurization rate matched the launch pressure profile.

ETL↗

Impacts of urbanization, antecedent rainfall event, and cyclone tracks on extreme floods at Houston reservoirs during Hurricane Harvey

The objective of this study is to evaluate the effects of urbanization, antecedent rainfall event (ARE), and varying cyclone tracks on the streamflow—and thus the subsequent reservoir status—during the floods caused by Hurricane Harvey in August-September 2017. Through a hydrological modeling approach, we examined how these factors influenced the inflows, peak pool elevations and outflows of the two most important detention reservoirs in the Houston region, the Addicks and Barker Reservoirs. A high-resolution rainfall reanalysis dataset for extreme storm events, along with a suite of synthetic rainfall values from a variety of storm tracks, were adopted to represent both the truth and the maximum possible rainfall during the Hurricane Harvey period. Results showed the following: Urbanization only led to slight increases in peak inflows, not necessarily leading to an increase in peak pool elevations. ARE contributed to the peak inflow and pool elevation slightly. In contrast, if the cyclone had followed the most adverse track consistent with earlier forecasts (all else being equal), the total volumetric flow into the two reservoirs could have been significantly larger (37.00% and 48.70% respectively), thus increasing the peak pool elevations by 1.06 and 1.37 meters respectively. These results suggest that large uncertainties exist for flood management at a watershed scale during hurricanes, because of the uncertainties with regard to the cyclone track. This is true even if storm-relative precipitation rates could be predicted perfectly.

54 ENVIRONMENTAL SCIENCES↗

Comparing Exciton-Polariton Couplings in Two-Dimensional MoS2 Cavities

Exciton-polaritons (EPs) are an emerging approach to achieve strong light-matter coupling and robust quantum phenomena, with application to areas such as photonic neuromorphic computing for instance reservoir computing. 2D transition metal dichalcogenides (TMDs) are promising candidates to host EPs due to their potential for strong coupling as evidenced by large Rabi splitting energies, which in the case of few-layer MoS2 has been demonstrated to reach 293 meV when coupled to the C exciton, placing it in the ultrastrong coupling regime, at room temperature. Here we examine how the Rabi splitting of CVD-grown monolayer MoS2 in a cavity can be optimized by tuning two variables: cavity geometry (thickness of silver and dielectric spacers) and the targeted exciton for coupling (A at about 1.9 eV, B at about 2.06 eV, or C at about 2.9 eV). With the cavity geometry, there is a tradeoff between optical field confinement and transmission. With the targeted exciton, there may be a tradeoff between Rabi splitting and lifetime. Transfer matrix simulations are used to inform how thick the cavity layers need to be, and at what angle measurements should be performed, to achieve Rabi splitting at a certain exciton energy. Experimentally, the Rabi splitting is measured by placing the cavity at a specified angle and measuring the optical reflectance and transmittance. By fitting the data, we find Rabi splittings of 124 meV, 194 meV, and 288 meV for the A, B, and C excitons, respectively, which places them in the strong, and on the edge of the ultrastrong, coupling regime.

exciton-polaritons↗

Mapping of deciduous forest cover using simulated Landsat-D TM data

An evaluation is presented of the supervised vegetation classification images of heavily forested terrain in the eastern US produced from simulated Landsat-D Thematic Mapper data as part of the Joint NASA/Geosat Test Case study of the Lost River, West Virginia, gas field. This investigation utilized data supplied by the NS-001 aircraft multispectral scanner (15 m IFOV) of NASA. The instruments recognized a total of 9 vegetation classes and one soil class. A hybrid parallelepiped-Bayesian maximum likelihood classifier was employed to produce the supervised vegetation classifications. The data were obtained at the height of the fall foliage display for this portion of West Virginia. Results show that both classifications are highly accurate, based on a comparison of assignment of vegetation classes on each of the images with actual ground conditions. It is noted that the anomalous distribution of certain woody species found in the resulting supervised classifications may be related to the microseepage of methane from the gas reservoir.

Rock, B. N.↗

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗

Contrast Between 1992 and 1997 High Latitude Spring Haloe Observations of Lower Stratospheric HCl

HCl measurements from HALOE in the northern hemisphere during mid-May 1997 revealed vortex fragments in which the chlorine reservoir partitioning was strongly pushed toward HCl (approx. 90% HCl, approx. 10% ClONO2), similar to partitioning previously observed in the Antarctic vortex region. In contrast, observations of ClONO2 and HCl in the northern polar spring, 1992, and in other years, show these species established the balance typical for gas phase photochemical reactions in this region (approx. 60% HCl, approx. 40% ClONO2). Annually, chlorine reservoirs in the winter lower stratosphere polar vortex are converted to chlorine radicals via heterogeneous reactions on particle surfaces at very cold temperatures (less than about 200 K). As temperatures warm in spring, the heterogeneous processes become insignificant compared with gas phase reactions, and the chlorine reservoirs are reformed. Measurements through the northern winter/spring in 1992 show rapid formation of ClONO2, followed by steady loss of ClONO2 and increasing HCl. Although ClONO2 measurements are not available for 1997, the HCl increase in 1997 is observed to be much more rapid and the eventual HCl mixing ratio is about 50% greater than that of 1992. The observations are examined through comparison with the Goddard three-dimensional chemistry and transport model. This model utilizes winds and temperatures from the Goddard Earth Observing System Data Assimilation System and a complete integration scheme for stratospheric photochemistry. Analysis of the evolution of HCl and ClONO2 shows that the observed difference in the overall rate of HCl formation is explained by the sensitivity of the gas-phase chemistry to the ozone mixing ratio and the temperature. The results show that the model accurately simulates HCl and ClONO2 evolution during these two winters. Model validity is further supported by comparisons with O3 and reactive nitrogen species NO and NO2. This analysis provides a sensitive test of the lower stratospheric chlorine photochemistry, particularly because the analysis considers constituent evolution at a time when the HCl and ClONO2 are far from a photochemical stationary state.

Douglass, A. R.↗

Observations of stratospheric hydrogen fluoride by halogen occultation experiment (HALOE)

The Halogen Occultation Experiment (HALOE) Hydrogen Fluoride (HF) channel on the Upper Atmospheric Research Satellite (UARS) is providing the first global measurements of stratospheric HF, the dominant flourine reservoir in the atmosphere. This paper describes the latitudinal and seasonal variations of HALOE-observed HF in terms of vertical profiles, altitude/latitude cross sections, and column abundances. The HF global distribution shows a 'tracerlike' structure and its column amount increases with latitude, in agreement with previous aircraft measurements of the HF column amount. A comparison between the HALOE HF column above 20 km and the ATMOS 1985 measurements is used to estimate the annual rate of increase of stratospheric HF. Exponential rates of 4.9-6.6%/yr and linear growth rates of 6-8.6%/yr in 1985 and 4.3-5.5%/yr in 1992-1993 are found. HALOE HF measurements during the 1993 Antarctic spring are briefly described. This species behaves like a conserved tracer and its distribution shows an area of enhanced mixing ratios correlated with the polar vortex that has a clear latitude boundary. Finally, simulated HF distributions by the National Center for Atmospheric Research (NCAR) two-dimensional model are used to compare with HALOE observations of HF. Reasonable agreements in the global structure and the absolute amount of HF are found. The differences between the model and the observed results indicate the need for improving treatment of atmospheric dynamics and fluorine-related chemical parameters in the model simulations.

Luo, M.↗

Jackson, L., Johnson, M.B., Latrach, A., Grimes, D., Martinez, C., and Mclaughlin, J.F., 2024, Multidisciplinary geotechnical data collection, curation, and analysis for conformity with the regulatory framework for geologic carbon storage in Wyoming, USA: Geological Society of America Abstracts with Programs. Vol. 56, No. 5, 2024, doi: 10.1130/abs/2024AM-405024

Title: Multidisciplinary Geotechnical Data Collection, Curation, and Analysis for Conformity with the Regulatory Framework for Geologic Carbon Storage in Wyoming, USA. Text: Construction and operation of wells for geologic sequestration of carbon dioxide necessitate that they are permitted under the Environmental Protection Agency’s Underground Injection Control Class VI requirements. Class VI wells conform to stringent requirements to ensure long-term safety and integrity of the storage site and the protection of Underground Sources of Drinking Water. Entities pursuing Class VI permitting must provide comprehensive geologic site characterization, including regional geologic structure and stratigraphy, aquifer information, reservoir and confining unit geomechanical properties, geochemical analyses, assessment of trapping capacity and mechanisms, and a variety of other of multidisciplinary geotechnical data. The Wyoming Class VI Site Characterization Database Project is focused on developing a geologic site characterization database of geotechnical information, which has been compiled and verified from established, public databases/entities and scientific literature to expedite Class VI permitting in Sweetwater County within the Greater Green River Basin of southern Wyoming. The preliminary suite of compiled data from 14,000 wells includes 8,000 wells with logs and 7,250 wells with formation tops, ~70 wells with core data (e.g., X-Ray diffraction, petrographic, and petrophysical data), ~2,500 water analyses, ~740 seismic events data, and ~520 bottom-hole temperature measurements. Future work on—and stemming from—this project will include new core analyses, calculation and interpolation of subsurface temperature gradients, mechanical earth models, geochemical simulations, storage capacity estimation, stratigraphic column generation and correlation, and construction of subsurface maps. Finally, this work will help to inspire and facilitate subsurface data compilation and curation beyond Sweetwater County, Wyoming.

42 ENGINEERING↗

Reservoir Thermal Energy Storage Benchmarking (Rev. 3)

A benchmarking analysis of RTES research funded by GTO through the Beyond Batteries projects was conducted against the ESGC to see where they fit within the identified ESGC Use Cases. The projects were found to advance knowledge in multiple ESGC use cases, either directly or in some cases, indirectly as enabling technologies. This analysis is helpful to understand where RTES and associated research fits into the larger discussion around energy storage technologies. Also, a retrospective analysis of the Beyond Batteries projects was conducted to evaluate what the projects learned and how the results can be applied to advance the value of RTES. Major results of each of the studies are summarized in Table 2. Additionally, a comparative metrics analysis for RTES was completed to understand where RTES lies within the energy storage industry. Metrics for evaluation of RTES and its comparison to other storage technologies were selected and ranges of their values compiled. The selected metrics – LCOE (levelized cost of energy), capital costs, roundtrip efficiency, energy storage capacity, and storage time – were chosen based on data availability and have a particularly strong influence on the potential deployment of a storage technology. Charts which compare the metrics are presented in section 4.3 and show ranges for each of the 10 selected technologies. However, due to a lack of domestic operational facilities, values for RTES and for portions of the remaining technologies are based on theoretical modeling and studies of best-case scenarios. LCOE estimates for RTES fall within the lower reaches of Figure 15, but nevertheless amount to 2 – 5 times the ESGC Roadmap goal for LCOE, for example in the Facilitating and Evolving Grid Use Case. Capital costs for RTES sit on the higher end (Figure 16) but are expected to decrease as new projects are developed and the technology is refined. The theoretical roundtrip efficiency reported for RTES varies from mid to high percentages (Figure 17) with efficiencies upwards of 93% in modeled scenarios in the Portland Basin (Bershaw et al.,2020). RTES is also expected to have the largest energy storage capacities and longest storage times, likely matched only by lower efficiency hydrogen storage. To better assess the role that RTES could play in energy storage we examined it’s potential in the U.S. The potential depends on many factors. Recently, many researchers have started looking at deep sedimentary basins, depleted oil and gas fields, and basalt formations as potential targets for RTES development. The United States Geological Survey (USGS) has analyzed various cities and shown substantial RTES potential in the cooling sector (Pepin et al., 2021). By modeling RTES in low-quality groundwater (e.g., brackish), it is shown to be favorable across the U.S. with particular suitability in the Illinois Basin, Coastal Plains, and Basin and Range regions. Seasonal RTES operations have also been modeled in the Portland Basin by those at the USGS and Portland State University to simulate an RTES system supplying heating loads needed for the Oregon Health and Science University. Simulations suggest that high conductive heat loss in the initial years exists but tends to decrease with increasing time and development of the resource due to self-insulating nature of the basalts (Burns et al., 2020). Other national laboratory efforts are taking a close look at many of the technical issues involved with RTES (McLing et al., 2019, McLing et al., 2022). These include difficulties in understanding geochemical, hydrogeological, mechanical, and microbiological changes at such elevated temperatures and operational scenarios. Major gaps in research are identified and suggested for future work. With this increased focus to understand how to make RTES successful in the U.S., this technology could be a potential solution to many of the nation’s energy storage problems. For the energy independence of this country, the DOE should prioritize de-risking this technology by making future investments in pilot-scale demonstrations to attract potential investors.

15 GEOTHERMAL ENERGY↗

Control development and sizing analysis for a 5th generation district heating and cooling network using Modelica

5th generation district heating and cooling system (5GDHC) are a relatively new concept. They use a single district loop near ambient temperature to provide heating and cooling. This paper improves on the modelling and control of a 5GDHC system called the reservoir network. The study updates the sewage heat exchanger plant model to more realistically represent seasonal changes, uses refined pump models with variable efficiency, introduces a ground coupled district pipe model to consider the inertia of the district network and implements a new control strategy for geothermal storage and sewage heat exchanger. The new approach reduced operating costs, mainly due to pumping cost for storage, sewage heat exchanger plant and distribution pump, while increasing the overall robustness of the approach in different sizing conditions. Thanks to the new controller, the pumping consumption was reduced by 21% with respect to the original baseline. Furthermore, the new control makes the system take better advantage of design changes, when reducing borehole field size and increasing the sewage heat exchanger size, the pumping energy savings become 29% with respect to the original baseline. Lastly, borehole field temperature stability was analyzed through 40 years of simulation.

Zanetti, Ettore↗