Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Network”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Extension of Generalized Fluid System Simulation Program's Fluid Property Database

This internship focused on the development of additional capabilities for the General Fluid Systems Simulation Program (GFSSP). GFSSP is a thermo-fluid code used to evaluate system performance by a finite volume-based network analysis method. The program was developed primarily to analyze the complex internal flow of propulsion systems and is capable of solving many problems related to thermodynamics and fluid mechanics. GFSSP is integrated with thermodynamic programs that provide fluid properties for sub-cooled, superheated, and saturation states. For fluids that are not included in the thermodynamic property program, look-up property tables can be provided. The look-up property tables of the current release version can only handle sub-cooled and superheated states. The primary purpose of the internship was to extend the look-up tables to handle saturated states. This involves a) generation of a property table using REFPROP, a thermodynamic property program that is widely used, and b) modifications of the Fortran source code to read in an additional property table containing saturation data for both saturated liquid and saturated vapor states. Also, a method was implemented to calculate the thermodynamic properties of user-fluids within the saturation region, given values of pressure and enthalpy. These additions required new code to be written, and older code had to be adjusted to accommodate the new capabilities. Ultimately, the changes will lead to the incorporation of this new capability in future versions of GFSSP. This paper describes the development and validation of the new capability.

Patel, Kishan↗
DOE OSTI2024

Improving the reliability of machine learned potentials for modeling inhomogeneous liquids

The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.

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

Relativistic analogue of the Newtonian fluid energy equation with nucleosynthesis

In Newtonian fluid dynamics simulations in which composition has been tracked by a nuclear reaction network, energy generation due to composition changes has generally been handled as a separate source term in the energy equation. Here, a relativistic equation in conservative form for total fluid energy, obtained from the spacetime divergence of the stress-energy tensor, in principle encompasses such energy generation; but it is not explicitly manifest. An alternative relativistic energy equation in conservative form—in which the nuclear energy generation appears explicitly, and that reduces directly to the Newtonian internal+kinetic energy in the appropriate limit—emerges naturally and self-consistently from the difference of the equation for total fluid energy and the equation for baryon number conservation multiplied by the average baryon mass m, when m is expressed in terms of contributions from the nuclear species in the fluid, and allowed to be mutable.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Flow Assisted Evaporative Cooling for Electric Motor

This paper examines a novel concept of flow assisted latent heat driven two-phase evaporative cooling (EC) confined in-between slot liner and active-winding of electric motor. Wicking micro-structure enhanced PDMS liner axially sucks coolant in the form of thin film between the PDMS liner and active-winding and eventually enables thin film evaporation on the outer surface of the active-winding. Therefore, EC based thermal management eliminates contact resistance between the winding and slot-liner and, enhances the heat extraction from the winding without compromising the electro-magnetic performance. Two-way coupled electro-magnetic (EM) – computational fluid dynamics/heat transfer (CFD/HT) and EM - lumped parameter thermal network (LPTN) models have been developed to assess the electro-thermal performance of the EC under steady and transient conditions. Taking a case study of a 125 kW jacket cooled BMW i3 motor and dielectric coolant FC-84, EC is shown to be capable of handling a maximum steady state rms current density of 26 A/mm2 at a evaporative heat transfer coefficient of 5,000 W/m2.K, which is about 78.7% higher compared to the traditional jacket cooling (JC). In case of EC, a maximum steady state and transient rms current density of 30 A/mm2 (106.2% higher compared to the JC), and 40.8 A/mm2 have been realized by using high thermal conductivity epoxy (1.9 W/m.K) impregnation material. Thermal performance of the EC is also assessed and compared with JC over a dynamic drive cycle. Lastly, a motorette testing has been performed to demonstrate the applicability of the proposed EC method and to validate the developed modeling framework.

42 ENGINEERING↗
DOE OSTI2016

A New Light Weight Structural Material for Nuclear Structures

Radiation shielding materials are commonly used in nuclear facilities to attenuate the background ionization radiations to a minimum level for creating a safer workplace, meeting regulatory requirements and maintaining high quality performance. The conventional radiation shielding materials have a number of drawbacks: heavy concrete contains a high amount of elements that are not desirable for an effective shielding such as oxygen, silicon, and calcium; a well known limitation of lead is its low machinability and toxicity, which is causing a major environmental concern. Therefore, an effective and environmentally friendly shielding material with increased attenuation and low mass density is desirable. Close-cell composite metal foams (CMFs) and open-cell Al foam with fillers are light-weight candidate materials that we have studied in this project. Close-cell CMFs possess several suitable properties that are unattainable by conventional radiation shielding materials such as low density and high strength for structural applications, high surface area to volume ratio for excellent thermal isolation with an extraordinary energy absorption capability. Open-cell foam is made up of a network of interconnected solid struts, which allows gas or fluid media to pass through it. This unique structure provided a further motive to investigate its application as radiation shields by infiltrating original empty pores with high hydrogen or boron compounds, which are well known for their excellent neutron shielding capability. The resulting open-cell foam with fillers will not only exhibit light weight and high specific surface area, but also possess excellent radiation shielding capability and good processability. In this study, all the foams were investigated for their radiation shielding efficiency in terms of X-ray, gamma ray and neutron. X-ray transmission measurements were carried out on a high-resolution microcomputed tomography (microCT) system. Gamma-emitting sources: 3.0mCi 60 Co, 1.8mCi 137 Cs, 13.5mCi 241 Am, and 5.0mCi 133 Ba were used for gamma-ray attenuation analysis. The evaluations of neutron transmission measurements were conducted at the Neutron Powder Diffractometer beam facility at North Carolina State University. The experimental results were verified theoretically through XCOM and Monte Carlo Z-particle Transport Code (MCNP). A mechanical investigation was performed by means of quasi-static compressive testing. Thermal characterizations were carried out through effective thermal conductivity and thermal expansion analyses in terms of high temperature guarded-comparative-longitudinal heat flow technique and thermomechanical analyzer (TMA), respectively. The experimental results were compared with analytical results obtained from, respectively, Brailsford and Major’s model and modified Turner’s model for verification. Flame test was performed in accordance with United States Nuclear Regulatory Commission (USNRC) standard. CMF sample and a 304L stainless steel control sample were subjected to a fully engulfing fire with an average flame temperature of 800°C for a period of 30 minutes. Finite Element Analysis was conducted to secure the credibility of the experimental results. This research indicates the potential of utilizing the light-weight close-cell CMFs and open-cell Al foam with fillers as shielding material replacing current heavy structures with additional advantage of high-energy absorption and excellent thermal characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Final Report: Mapping Fracture Network Creation with Microseismicity During EGS Demonstrations

Effective enhanced geothermal system (EGS) sites require optimal fracture networks for efficient heat transfer between hot rock and fluids. Microseismic mapping is a key tool used to infer the subsurface fracture geometry and is widely used in the petroleum industry to monitor fracture propagation during hydraulic stimulations. Traditional earthquake detection techniques (e.g. STA/LTA detectors) and location techniques are often employed to identify microearthquakes in geothermal regions. However, most commonly used algorithms may miss events if the seismic signal of an earthquake is small relative to the background noise level or if a microearthquake occurs within the coda of a larger event. Consequently, we have developed a set of algorithms that provide improved microearthquake detection even in difficult signal processing environments.

58 GEOSCIENCES↗
NASA NTRS1997

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

Applications are described of high-performance computing methods to the numerical simulation of complete jet engines. The methodology focuses on the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field elements. New partitioned analysis procedures to treat this coupled three-component problem were developed. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. The NASA-sponsored ENG10 program was used for the global steady state analysis of the whole engine. This program uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 was developed as well as the capability for the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames.

Felippa, C. A.↗
NASA NTRS2012

NASA Tech Briefs, March 2012

The topics include: 1) Spectral Profiler Probe for In Situ Snow Grain Size and Composition Stratigraphy; 2) Portable Fourier Transform Spectroscopy for Analysis of Surface Contamination and Quality Control; 3) In Situ Geochemical Analysis and Age Dating of Rocks Using Laser Ablation-Miniature Mass Spectrometer; 4) Physics Mining of Multi-Source Data Sets; 5) Photogrammetry Tool for Forensic Analysis; 6) Connect Global Positioning System RF Module; 7) Simple Cell Balance Circuit; 8) Miniature EVA Software Defined Radio; 9) Remotely Accessible Testbed for Software Defined Radio Development; 10) System-of-Systems Technology-Portfolio-Analysis Tool; 11) VESGEN Software for Mapping and Quantification of Vascular Regulators; 12) Constructing a Database From Multiple 2D Images for Camera Pose Estimation and Robot Localization; 13) Adaption of G-TAG Software for Validating Touch and Go Asteroid Sample Return Design Methodology; 14) 3D Visualization for Phoenix Mars Lander Science Operations; 15) RxGen General Optical Model Prescription Generator; 16) Carbon Nanotube Bonding Strength Enhancement Using Metal Wicking Process; 17) Multi-Layer Far-Infrared Component Technology; 18) Germanium Lift-Off Masks for Thin Metal Film Patterning; 19) Sealing Materials for Use in Vacuum at High Temperatures; 20) Radiation Shielding System Using a Composite of Carbon Nanotubes Loaded With Electropolymers; 21) Nano Sponges for Drug Delivery and Medicinal Applications; 22) Molecular Technique to Understand Deep Microbial Diversity; 23) Methods and Compositions Based on Culturing Microorganisms in Low Sedimental Fluid Shear Conditions; 24) Secure Peer-to-Peer Networks for Scientific Information Sharing; 25) Multiplexer/Demultiplexer Loading Tool (MDMLT); 26) High-Rate Data-Capture for an Airborne Lidar System; 27) Wavefront Sensing Analysis of Grazing Incidence Optical Systems; 28) Foam-on-Tile Damage Model; 29) Instrument Package Manipulation Through the Generation and Use of an Attenuated-Fluent Gas Fold; 30) Multicolor Detectors for Ultrasensitive Long-Wave Imaging Cameras; 31) Lunar Reconnaissance Orbiter (LRO) Command and Data Handling Flight Electronics Subsystem; and 32) Electro-Optic Segment-Segment Sensors for Radio and Optical Telescopes.

Source record↗
NASA NTRS2019

Multi-Node Modeling of Cryogenic Tank Pressurization System using Generalized Fluid System Simulation Program

This paper presents a multi-node model of autogenous pressurization of cryogenic propellant in a flight tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/Marshall Space Flight Center. Tests were conducted to measure the pressure and temperatures at the various axial locations of the stratified ullage at 75% and 45% fill level. Liquid nitrogen was pressurized by gaseous nitrogen from a supply tank while the drain valve from the tank remained closed during the pressurization process. The ullage was discretized into 25 uniformly distributed nodes: 5 in the radial direction and 5 in the axial direction assuming the flow to be axisymmetric. Heat and mass transfer between the liquid and vapor has been modeled at the liquid vapor interface. Heat transfer between wall and vapor at the ullage has been accounted for by assuming heat transfer occurs by natural convection. The model also accounts for heat leak to the tank through the insulation and metal wall by heat conduction. The predicted pressures and temperatures are compared with the measured data.

Bandyopadhyay, Alak↗
DOE OSTI2023

A Simulation Framework for Analyzing the Impact of Stochastic Occupant Behaviors on Demand Flexibility in Typical Commercial Buildings

As one of the primary users of the electric grid, buildings and building equipment, including heating, ventilation, and air conditioning (HVAC) systems, can be leveraged to provide the flexible demand needed to balance the grid. Typical strategies to achieve demand flexibility are to reduce electricity use during peak or critical periods by shutting down equipment or relaxing system setpoints, which will inevitably impact the occupants’ comfort. When occupants feel uncomfortable, they may take actions to regain their comfort, and some of those actions (such as turning on a personal fan) may have a negative impact on meeting the demand response goal. Therefore, it is important to incorporate occupant behaviors into the assessment ofthe building demand flexibility potential. In this study, a simulation framework that includes simulation of zone thermal loads, an HVAC system, and occupant behaviors, was developed to investigate the impact of occupant behaviors on demand flexibility. A case study was conducted using a small office model from the U.S. Department of Energy (DOE) Commercial Prototype Building Models to simulate the building envelope and zone loads. An agent-based occupant thermal behavior model was adapted to forecast occupants’ thermal comfort and their resulting thermal behaviors. An artificial neural network (ANN) based airflow model trained from a computational fluid dynamics (CFD) model of the zone was adopted to better predict the ambient environment of each occupant. An air-source heat pump simulation model that was calibrated from a real two-stage air-source heat pump system was used as the HVAC system. A typical load shedding event during peak hours was studied. Repeated simulations were conducted to capture the stochastic effects of occupant behaviors. The interplay between the demand flexibility, occupant comfort and behavior were analyzed by evaluating key performance indicators, including the energy use, occupant discomfort duration, and occupant behavior duration during the peak period. The results suggest that this framework can be used to analyze typical commercial buildings and their HVAC systems in terms ofdemand flexibility potential under the impact of occupant behaviors.

Chen, Zhelun↗
DOE OSTI2026

Recurrent convolutional neural networks for modeling nonadiabatic dynamics of quantum-classical systems

Recurrent neural networks (RNNs) have recently been extensively applied to model the time evolution in fluid dynamics, weather predictions, and even chaotic systems due to their ability to capture temporal dependencies and sequential patterns in data. Here we present an RNN model based on convolutional neural networks for modeling the nonlinear nonadiabatic dynamics of hybrid quantum-classical systems. The dynamical evolution of the hybrid systems is governed by equations of motion for classical degrees of freedom and von Neumann equation for electrons. The Physics-Aware Recurrent Convolution (PARC) neural network structure incorporates a differentiator-integrator architecture that inductively models the spatiotemporal dynamics of generic physical systems. Here, we apply our RNN approach to learn the space-time evolution of a one-dimensional semiclassical Holstein model after an interaction quench. For shallow quenches (small changes in electron-lattice coupling), the deterministic dynamics can be accurately captured using a single-CNN-based recurrent network. In contrast, deep quenches induce chaotic evolution, making long-term trajectory prediction significantly more challenging. Nonetheless, we demonstrate that the PARC-CNN architecture can effectively learn the statistical climate of the Holstein model under deep-quench conditions.

Holstein model↗
DOE OSTI2021

Multifidelity computing for coupling full and reduced order models

Hybrid physics-machine learning models are increasingly being used in simulations of transport processes. Many complex multiphysics systems relevant to scientific and engineering applications include multiple spatiotemporal scales and comprise a multifidelity problem sharing an interface between various formulations or heterogeneous computational entities. To this end, we present a robust hybrid analysis and modeling approach combining a physics-based full order model (FOM) and a data-driven reduced order model (ROM) to form the building blocks of an integrated approach among mixed fidelity descriptions toward predictive digital twin technologies. At the interface, we introduce a long short-term memory network to bridge these high and low-fidelity models in various forms of interfacial error correction or prolongation. The proposed interface learning approaches are tested as a new way to address ROM-FOM coupling problems solving nonlinear advection-diffusion flow situations with a bifidelity setup that captures the essence of a broad class of transport processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Data-driven recovery of hidden physics in reduced order modeling of fluid flows

In this article, we introduce a modular hybrid analysis and modeling (HAM) approach to account for hidden physics in reduced order modeling (ROM) of parameterized systems relevant to fluid dynamics. The hybrid ROM framework is based on using first principles to model the known physics in conjunction with utilizing the data-driven machine learning tools to model the remaining residual that is hidden in data. This framework employs proper orthogonal decomposition as a compression tool to construct orthonormal bases and Galerkin projection (GP) as a model to build the dynamical core of the system. Our proposed methodology hence compensates structural or epistemic uncertainties in models and utilizes the observed data snapshots to compute true modal coefficients spanned by these bases. The GP model is then corrected at every time step with a data-driven rectification using a long short-term memory (LSTM) neural network architecture to incorporate hidden physics. A Grassmannian manifold approach is also adopted for interpolating basis functions to unseen parametric conditions. The control parameter governing the system's behavior is thus implicitly considered through true modal coefficients as input features to the LSTM network. The effectiveness of the HAM approach is then discussed through illustrative examples that are generated synthetically to take hidden physics into account. Furthermore, our approach thus provides insights addressing a fundamental limitation of the physics-based models when the governing equations are incomplete to represent underlying physical processes.

42 ENGINEERING↗
DOE OSTI2022

Neural network–based pore flow field prediction in porous media using super resolution

Direct pore-scale simulations of fluid flow through porous media are computationally expensive to perform for realistic systems. Previous works have demonstrated using the geometry of the microstructure of porous media to predict the velocity fields therein based on neural networks. However, such trained neural networks do not perform well for unseen porous media with a large degree of heterogeneity. In this study we propose that incorporating a coarse velocity field in the input of neural networks is an effective way to improve the prediction performance. The coarse velocity field can be simulated with a low computational cost and provides global information to regularize the ill-posedness of the learning problem, which is usually caused by the use of local geometries due to the computational resource constraints. We show that incorporating the coarse-mesh velocity field significantly improves the prediction accuracy of the fine-mesh velocity field by comparison to the prediction that relies on geometric information alone, especially for the porous medium with a large interior vuggy pore space. We also show the flexibility of training the network in using coarse velocity fields with various resolutions. The results suggest that even using coarse velocity field with a very low resolution, the predictions are still enhanced and close to the ground truths. The feasibility of the method is further demonstrated by testing the trained network on real rocks. This study highlights the merits of incorporating a coarse-mesh velocity field into the input for neural networks, which provides global, physics-based information for the model, thereby improving the model's generalization capability.

42 ENGINEERING↗
DOE OSTI2023

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

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

Geomechanics of CO 2 Reservoir Seals

Chemically-assisted subcritical fracture growth within overlying shale caprock can threaten seal integrity of CO 2 sequestration reservoirs. Here, we integrate laboratory fracture mechanics testing conducted under a wide range of environmental conditions including dry and wet supercritical CO 2 , dry CO 2 gas, ambient air, deionized water, and aqueous solutions of varying pH and salinity, with numerical simulations of fracture network growth and caprock degradation calculations, to investigate the coupled chemo-mechanical processes associated with potential fracture-controlled leakage from subsurface CO 2 storage. Double torsion tests quantify in detail the influence of aqueous fluid chemistry and temperature on subcritical fracture growth of caprock at atmospheric pressure, while short rod tests derive fracture toughness of caprock under pressurized supercritical CO 2 environments. We find that lowered pH as expected after brine acidification enhances subcritical fracture growth properties in carbonate-rich shale. An increase in salinity that may result from scCO 2 -brine interaction would potentially suppress fracture growth in clay-rich shales. Dehydration by dry supercritical CO 2 significantly strengthens shale lithologies. Through boundary element fracture network models that account for subcritical fracture growth we demonstrate that changes in aqueous fluid composition can affect fracture network geometry and connectivity, and thus potentially degrading caprock efficiency at preventing upward migration of CO 2 . Calculations based on boundary element simulation suggest that caprock degradation rates may increase by 58% following acidification of formation water in contact with carbonate-rich caprock. With spatially and temporally variable CO 2 -brine interaction in storage reservoirs, chemo-mechanical interactions may control the location of potential leakage pathway evolution, with preexisting faults and reservoir heterogeneities providing the preferred locus for environmentally assisted fracture growth and subsequent caprock failure.

58 GEOSCIENCES↗
DOE OSTI2023

Utah FORGE 1-2551: Multi-Stage Fracturing System and Well Tractor to Enable Zonal Isolation During Stimulation and EGS Operations in Horizontal Wellbores - Workshop Presentation and Report

Included here are a presentation recording, slides, and report on the Multi-Stage Fracturing System and Well Tractor to Enable Zonal Isolation During Stimulation and EGS Operations in Horizontal Wellbores project by Colorado School of Mines, presented by Dr. William W. Fleckenstein. One goal of the project is to develop and demonstrate the use of low-cost, rapid multistage fracture stimulation technology with cemented casing frac sleeves that eliminates packers used in both conventional stimulation and injection flow control. The second goal is to develop and demonstrate the devices to effectively detect and control the flow of heat-carrier fluid solely from a long-reach injector through the network of induced and existing fractures and produced from an open-hole long-reach well to improve heat recovery. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 7, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗