Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D Modeling”

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

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

At least 397 records · Page 22

SU-E-J-196: New Visualization Methods for Longitudinal MRI Registrations and Segmentations

Purpose: To develop visualization techniques to facilitate easy assessment of (a) registration and (b) tracking volumetric changes in structures during radiation therapy from MRI. Method: The frequently used method for visualizing registrations between scans is a multi-color overlay technique or deformation vector fields. However, the overlay technique is unintuitive and does not help to appreciate the quality of registration particularly when the registration mismatches are not very large. Similarly, the deformation fields give an indication of extent of deformation but do not help to assess the differences in registration. We present a mirroring and edge-augmented mirroring technique that places the fixed and moving image next to each other and allows the user to quickly assess the small differences in registration. Next, we present a volumetric intersection based 3D model to visualize the changes in diseased lymph node volumes in head and neck cancer. 3D model-based visualization provides more information about the location-specific changes in volume rather than the simplistic one dimensional information obtained from 2D plot of nodal volume changes. Result: We show results comparing our approach with the standard colorbased overlay method for comparing registrations of intra-patient registrations using T2-MRI. Upon comparing the mirroring technique with the color-overlay, one can more easily appreciate the differences in registration. Adding edge-based mirroring seems to further assist in evaluating the registration. Our approach for viewing registrations seems to be more intuitive and easy to use in order to help assess the quality of registration compared to color-based overlays. Similarly, the change volumetric model together with a 2D plot reveals more information including the locations undergoing changes and responding to treatment. Conclusions: Better approaches are necessary for assessing the quality of registrations and changes in diseased structures over time. Our approach presents some novel visualization solutions for registration and segmentation.

60 APPLIED LIFE SCIENCES↗

Skeletides: A Modular, Simplified Physical Model of Protein Secondary Structure

Three-dimensional (3D) models are essential for visualization and conceptual understanding of complex architectures such as protein structure. Although there is a plethora of software platforms that allow digital depictions of protein structure at an atomic level in silico, physical models are needed to convey an intuitive understanding of biomolecular architecture. However, it is a challenge to represent all the relevant features of proteins in a single physical model due to their sheer structural complexity. Here, we describe a modular protein model that focuses only on representation of the secondary structure—the underlying structural skeleton. The simplified model consists of amino acid units, which can be linked together to reproduce the two most fundamental structural features of protein secondary structure: the relative positions of the amino acid alpha carbon atoms, and the intra-main chain hydrogen bonding pattern. We use 3D printing and magnets to create a set of three modular amino acid building blocks, which when linked together into a chain, can faithfully represent alpha helices, beta sheets, and turns. These simple to make models can be used to quickly assemble the alpha carbon trace of an entire protein domain, which conveys a tactile experience of the complexity of protein skeletal architecture. These models highlight the modularity of protein structure: where a single structural unit, the amino acid, can be linked together to form a larger, regular secondary structure. These models also have a propensity to spontaneously organize into alpha helices and beta sheets, as demonstrated by their ability to autonomously assemble when placed in a circulating water tank. These models provide a missing educational tool to expand knowledge of protein structure, foster deeper insight into protein folding, and inspire greater interest in biomacromolecular architecture.

36 MATERIALS SCIENCE↗

Study of the muon component in the core-corona model using CONEX 3D

The discrepancy between models and data regarding the muon content in air showers generated by ultra-high energy cosmic rays still needs to be solved. The CONEX simulation framework provides a flexible tool to assess the impact of different interaction properties and thus address the muon puzzle. In this work, we present the multidimensional extension of CONEX and show its performance compared to CORSIKA by discussing muon-related air-shower features for three experiments: KASCADE, IceTop, and the Pierre Auger Observatory. We also implement an effective version of the core-corona model to demonstrate the impact of the core effect, as observed at the LHC, on the muon content in air showers produced by ultra-high energy cosmic rays. At a primary energy of E$_{0}$ = 10$^{19}$ eV, we obtain an increase of 15% to 20% in the muon content.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Utah FORGE: Seismic Velocity Models, February 2021

This dataset contains a map, showing the Utah FORGE seismic stations, and seismic velocity model data. There are 61 1-D velocity models which are in a compressed TAR file. A paper is referenced at the end of this description which discusses the use of these data in 3D modelling. The paper summary follows: We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2-5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ~22 degrees to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity. Reference: Zhang, H. and K. L. Pankow (2021). High-resolution Bayesian spatial auto-correlation (SPAC) pseudo-3D Vs model of Utah FORGE site with a dense geophone array, Geophys. Res. Int, https://doi.org/10.1093/gji/ggab049

15 GEOTHERMAL ENERGY↗

Supernova Explosions of the Lowest-mass Massive Star Progenitors

Abstract We here focus on the behavior of supernovae that technically explode in 1D (spherical symmetry). When simulated in 3D, however, the outcomes of representative progenitors of this class are quite different in almost all relevant quantities. In 3D, the explosion energies can be 2 to 10 times higher, and there are correspondingly large differences in the 56 Ni yields. These differences between the 3D and 1D simulations reflect in part the relative delay to explosion of the latter and in the former the presence of protoneutron star convection that boosts the driving neutrino luminosities by as much as ∼50% at later times. In addition, we find that the ejecta in 3D models are more neutron-rich, resulting in significant weak r -process and 48 Ca yields. Furthermore, we find that in 3D the core is an interesting, though subdominant, source of acoustic power. In summary, we find that though a model might be found theoretically to explode in 1D, one must perform supernova simulations in 3D to capture most of the associated observables. The differences between 1D and 3D models are just too large to ignore.

Astronomy & Astrophysics↗

batteryNET v0.0.1

Code to organize, prepare, and segment multi-scan, larger-than-memory, volumetric dataset. Deep learning segmentation uses pytorch-lightning for code organization, logging and metrics, and multi-GPU parallelization of training and prediction. Also uses moani for volumetric data augmentation and model implementation (currently U-Net++). Training can be done slice-by-slice with 2D models trained on 2D patches either aligned with or perpendicular to the electrode plane, or with 3D models on 3D patches. Sparse labels can be used to minimize required hand-labeling -- less than 0.1% of the entire dataset was labeled. Also includes code for data preparation (alignment/cropping/patching).

Ushizima, Daniela↗

Examination of Dendritic Growth During Solidification of Ternary Alloys via a Novel Quantitative 3D Cellular Automaton Model

A three-dimensional (3D) quantitative cellular automaton (CA) model was developed to simulate dendritic growth during solidification processing of ternary alloys. A detailed method was proposed to solve solute diffusion and calculate the solid fraction for ternary alloys during solidification. The present model has been shown to accurately predict dendrite morphologies and solute distributions of both single equiaxed dendrite and columnar dendrites during solidification. The model was also used to study the influence of the concentration of a third component (Mg in the Al-Si-Mg system) on dendritic growth. With increasing Mg concentration, the steady-state dendrite tip growth velocity was shown to decrease, resulting in a shorter primary dendrite length. The 3D CA simulation results agree well with the prediction of the LGK theoretical model. Multi-columnar dendrite growth was simulated with different cooling rates. Primary and secondary dendrite morphology was shown to have little variation at low cooling rates (1.0 to 2.0 K/s) and a constant undercooling. The average secondary dendrite arm spacing was shown to decrease with the increase of cooling rate within this range. The 3D CA simulation results are in good agreement with the experimental directional solidification experiments of a commercial A356 (Al-7 wt pct Si-0.3 wt pct Mg).

36 MATERIALS SCIENCE↗

A global three-dimensional model of the stratospheric sulfuric acid layer

A 3D model which encompasses SO2 production from OCS, followed by its oxidation to gaseous H2SO4, the condensation-evaporation equilibrium of gaseous and particulate H2SO4, and finally particle condensation and rainout, is presently used to study processes maintaining the nonvolcanically-perturbed stratosphere's sulfuric acid layer. A comparison of the results thus obtained with remotely sensed stratospheric aerosol extinction data shows the model to simulate the general behavior of stratospheric aerosol extinction.

Golombek, Amram↗

SURROGATE MODEL DEVELOPMENT AND VALIDATION FOR RELIABILITY ANALYSIS OF REACTOR PRESSURE VESSELS

In nuclear light water reactors (LWRs), the reactor coolant, core and shroud are contained within a massive, thick walled steel vessel known as a reactor pressure vessel (RPV). Given the tremendous size of these structures, RPVs typically contain a large population of pre-existing flaws introduced in the manufacturing process. After many years of operation, irradiation-induced embrittlement makes these vessels increasingly susceptible to fracture initiation at the locations of the pre-existing flaws. Because of the uncertainty in the loading conditions, flaw characteristics and material properties, probabilistic methods are widely accepted and used in assessing RPV integrity. The Fracture Analysis of Vessels – Oak Ridge (FAVOR) computer program developed by researchers at Oak Ridge National Laboratory is widely used for this purpose. This program can be used in order to perform deterministic and probabilistic risk-informed analyses of the structural integrity of an RPV subjected to a range of thermal-hydraulic events. FAVOR uses a one-dimensional representation of the global response of the RPV, which is appropriate for the beltline region, which experiences the most embrittlement, and employs an influence coefficient technique to rapidly compute stress intensity factors for axis-aligned surface-breaking flaws. The Grizzly code is currently under development at Idaho National Laboratory (INL) to be used as a general multiphysics simulation tool to study a variety of degradation mechanisms in nuclear power plant components. The first application of Grizzly has been to study fracture in embrittled RPVs. Grizzly can be used to model the thermo-mechanical response of an RPV under transient conditions observed in a pressurized thermal shock (PTS) scenario. The global response of the vessel provides boundary conditions for local 3D models of the material in the vicinity of a flaw. Fracture domain integrals are computed to obtain stress intensity factors, which can in turn be used to assess whether a fracture would initiate at a pre-existing flaw. To use Grizzly for probabilistic analysis, it is necessary to have a way to rapidly evaluate stress intensity factors. To accomplish this goal, a reduced order model (ROM) has been developed to efficiently represent the behavior of a detailed 3D Grizzly model used to calculate fracture parameters. This approach uses the stress intensity factor influence coefficient method that has been used with great success in FAVOR. Instead of interpolating between tabulated solutions, as FAVOR does, the ROM approach uses a response surface methodology to compute fracture solutions based on a sampled set of results used to train the ROM. The main advantages of this approach are that the process of generating the training data can be fully automated, and the procedure can be readily used to consider more general flaw configurations. This paper demonstrates the procedure used to generate a ROM to rapidly compute stress intensity factors for axis-aligned flaws. The results from this procedure are in good agreement with those produced using the traditional influence coefficient interpolation procedure, which gives confidence in this method. This paves the way for applying this procedure for more general flaw configurations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SU-E-J-01: 3D Fluoroscopic Image Estimation From Patient-Specific 4DCBCT-Based Motion Models

Purpose: 3D motion modeling derived from 4DCT images, taken days or weeks before treatment, cannot reliably represent patient anatomy on the day of treatment. We develop a method to generate motion models based on 4DCBCT acquired at the time of treatment, and apply the model to estimate 3D time-varying images (referred to as 3D fluoroscopic images). Methods: Motion models are derived through deformable registration between each 4DCBCT phase, and principal component analysis (PCA) on the resulting displacement vector fields. 3D fluoroscopic images are estimated based on cone-beam projections simulating kV treatment imaging. PCA coefficients are optimized iteratively through comparison of these cone-beam projections and projections estimated based on the motion model. Digital phantoms reproducing ten patient motion trajectories, and a physical phantom with regular and irregular motion derived from measured patient trajectories, are used to evaluate the method in terms of tumor localization, and the global voxel intensity difference compared to ground truth. Results: Experiments included: 1) assuming no anatomic or positioning changes between 4DCT and treatment time; and 2) simulating positioning and tumor baseline shifts at the time of treatment compared to 4DCT acquisition. 4DCBCT were reconstructed from the anatomy as seen at treatment time. In case 1) the tumor localization error and the intensity differences in ten patient were smaller using 4DCT-based motion model, possible due to superior image quality. In case 2) the tumor localization error and intensity differences were 2.85 and 0.15 respectively, using 4DCT-based motion models, and 1.17 and 0.10 using 4DCBCT-based models. 4DCBCT performed better due to its ability to reproduce daily anatomical changes. Conclusion: The study showed an advantage of 4DCBCT-based motion models in the context of 3D fluoroscopic images estimation. Positioning and tumor baseline shift uncertainties were mitigated by the 4DCBCT-based motion models, while they caused errors when using 4DCT-based motion models. Partially funded by Varian research grant.

60 APPLIED LIFE SCIENCES↗

Snake River Plain FORGE: Site Characterization Data

The site characterization data used to develop the conceptual geologic model for the Snake River Plain site in Idaho, as part of phase 1 of the Frontier Observatory for Research in Geothermal Energy (FORGE) initiative. This collection includes data on seismic events, groundwater, geomechanical models, gravity surveys, magnetics, resistivity, magnetotellurics (MT), rock physics, stress, the geologic setting, and supporting documentation, including several papers. Also included are 3D models (Petrel and Jewelsuite) of the proposed site. Data for wells INEL-1, WO-2, and USGS-142 have been included as links to separate data collections. These data have been assembled by the Snake River Geothermal Consortium (SRGC), a team of collaborators that includes members from national laboratories, universities, industry, and federal agencies, lead by the Idaho National Laboratory (INL). Other contributors include the National Renewable Energy Laboratory (NREL), Lawrence Livermore National Laboratory (LLNL), the Center for Advanced Energy Studies (CEAS), the University of Idaho, Idaho State University, Boise State University, University of Wyoming, University of Oklahoma, Energy and Geoscience Institute-University of Utah, US Geothermal, Baker Hughes Campbell Scientific Inc., Chena Power, US Geological Survey (USGS), Idaho Department of Water Resources, Idaho Geological Survey, and Mink GeoHydro.

15 GEOTHERMAL ENERGY↗

Neutron spectra from beam-target reactions in dense Z-pinches

The energy spectrum of neutrons emitted by a range of deuterium and deuterium-tritium Z-pinch devices is investigated computationally using a hybrid kinetic-MHD model. 3D MHD simulations are used to model the implosion, stagnation, and break-up of dense plasma focus devices at currents of 70 kA, 500 kA, and 2 MA and also a 15 MA gas puff. Instabilities in the MHD simulations generate large electric and magnetic fields, which accelerate ions during the stagnation and break-up phases. A kinetic model is used to calculate the trajectories of these ions and the neutron spectra produced due to the interaction of these ions with the background plasma. It is found that these beam-target neutron spectra are sensitive to the electric and magnetic fields at stagnation resulting in significant differences in the spectra emitted by each device. Most notably, magnetization of the accelerated ions causes the beam-target spectra to be isotropic for the gas puff simulations. It is also shown that beam-target spectra can have a peak intensity located at a lower energy than the peak intensity of a thermonuclear spectrum. A number of other differences in the shapes of beam-target and thermonuclear spectra are also observed for each device. Finally, significant differences between the shapes of beam-target DD and DT neutron spectra, due to differences in the reaction cross-sections, are illustrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of Electrochemical Recovery of Zirconium from Spent Nuclear Fuels

This project uses both modeling and experimental studies to design optimal electrochemical technology methods for recovery of zirconium from used nuclear fuel rods for more effective waste management. The objectives are to provide a means of efficiently separating zirconium into metallic high-level waste forms and to support development of a process for decontamination of zircaloy hulls to enable their disposal as low- and intermediate-level waste. Modeling work includes extension of a 3D model previously developed by Seoul National University for uranium electrorefining by adding the ability to predict zirconium behavior. Experimental validation activities include tests for recovery of zirconium from molten salt solutions and aqueous tests using surrogate materials. *This is a summary of the FY 2013 progress for I-NERI project # 2010-001-K provided to the I-NERI office.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Managing an Aging Pump-and-Treat Groundwater Remediation System at the Fernald Preserve - 19557

The Fernald Preserve, Ohio, Site underwent restoration and remediation under the Comprehensive Environmental Response, Compensation, and Liability Act (CERCLA). With the exception of the Great Miami Aquifer, the CERCLA remediation was declared physically complete on October 29, 2006, and the site was turned over to the US Department of Energy Office of Legacy Management for long-term care and maintenance. The Record of Decision under CERCLA specifies a concentration-based cleanup that requires the extraction and treatment of Great Miami Aquifer groundwater above Final Remediation Levels (FRLs) until the full, beneficial use potential of the aquifer is achieved, including use as a drinking water source. The area of the aquifer targeted for remediation is 126.5 hectares (312.70 acres). The groundwater pumping remediation was designed by modeling cleanup with the Variable Saturated Analysis Model in 3 Dimensions (VAM-3D) modeling code, developed by Hydrogeologic Inc. The original design called for 23 extraction wells operating at 379 liters per minute (lpm) (100 gallons per minute [gpm]) to 1135 lpm (300 gpm) each for a system total of 18,073 lpm (4775 gpm). In 2014, the system was optimized down to 20 extraction wells operating at 379 to 1893 lpm each for a system total of 19,180 lpm (5075 gpm). To date, pumping has been an effective remediation technology for the aquifer. At the end of 2017, the area of the aquifer remaining above FRLs was down to 38.2 hectares (94.4 acres). From August 1993 through December 2017 (24 years), 6235 kilograms (kg) (13,733 pounds) of uranium were removed from the aquifer by the pumping of 1.74 x 10{sup 11} liters (46.2 billion gallons) of water. Since 2010, treatment of groundwater to achieve regulatory discharge limits at the Great Miami River has not been needed on a routine basis. The ongoing remediation is now essentially a pumping operation, and this water does not require treatment. Metrics reported in the annual Site Environmental Report for the groundwater remediation indicate that, based on the model-predicted estimate of the mass of uranium that needs to be removed from the aquifer to achieve a concentration-based cleanup, the cleanup is approximately 86% - 87% complete. Current model predictions indicate that pumping will need to continue an additional 16 years (until 2035) in some of the wells. During the additional 16 years, it is estimated that approximately 1012 kg (2231 pounds) of uranium will be removed from the aquifer. System pumping rates are designed to decrease over the course of the remaining 16 years, dropping down to a system total of 11,638 lpm (3075 gpm) in 2022 and a system total of 4163 lpm (1100 gpm) in 2033. Even with the designed decreases in pumping rates, it is estimated that approximately 1 x 10{sup 11} liters (26.6 billion gallons) of water will be pumped from the aquifer to remove 1012 kg (2231 pounds) of uranium. This indicates that the remaining 16 years of operation will be much less efficient than the first 24 years. Decreasing efficiency is common for pumping operations. The worst-case scenario is that decreasing concentration trends stabilize above the cleanup FRLs and remediation goals are never achieved. Performance assessments of the pumping remedy will need to be routinely conducted to determine if pumping should be discontinued, if it should be continued but optimized, or if the pumping technology can be supplemented with a different technology that will improve the outcome of the remediation and make the operation more efficient and successful. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Mechanism-based Representative Volume Elements (RVEs) for Predicting Property Degradations in Multiphase Materials

Quantitative understanding of the evolving thermal-mechanical properties of a multi-phase material hinges upon the availability of quantitative statistically representative microstructure descriptions. Questions then arise as to whether a two-dimensional (2D) or a three-dimensional (3D) representative volume element (RVE) should be considered as the statistically representative microstructure. Although 3D models are more representative than 2D models in general, they are usually computationally expensive and difficult to be reconstructed. In this paper, we evaluate the accuracy of a 2D RVE in predicting the property degradations induced by different degradation mechanisms with the multiphase solid oxide fuel cell (SOFC) anode material as an example. Both 2D and 3D microstructure RVEs of the anodes are adopted to quantify the effects of two different degradation mechanisms: humidity-induced electrochemical degradation and phosphorus poisoning induced structural degradation. The predictions of the 2D model are then compared with the available experimental measurements and the results from the 3D model. It is found that the 2D model, limited by its inability of reproducing the realistic electrical percolation, is unable to accurately predict the degradation of thermo-electrical properties. On the other hand, for the phosphorus poisoning induced structural degradation, both 2D and 3D microstructures yield similar results, indicating that the 2D model is capable of providing computationally efficient yet accurate results for studying the structural degradation within the anodes.

degradation mechanisms↗

Three-Dimensional Geologic Characterization of a Great Basin Geothermal System: Astor Pass, Nevada

The Great Basin, western USA, exhibits anomalously high heat flow (~75±5 mWm-2) and active faulting and extension, resulting in ~430 known geothermal systems. Recent studies have shown that steeply dipping normal faults in transtensional pull-aparts are a common structural control of these Great Basin geothermal systems. The Astor Pass blind (no surface expression) geothermal system, Nevada, lies along the boundary between the Basin and Range to the east and the Walker Lane to the west. Across this boundary, strain is transferred from dextral shear in the Walker Lane to west-northwest directed extension in the Basin and Range, resulting in a transtensional setting consisting of both northwest-striking, left-stepping dextral faults and northerly striking normal faults. Previous studies indicate that Astor Pass was controlled by the intersection of a northwest-striking dextral normal fault and north-northwest striking normal-dextral fault bounding the western side of the Terraced Hills. Drilling (to ~1200 m) has revealed fluid temperatures of ~94°C, confirming a blind geothermal system. Expanding upon previous work and employing interpretation of 2D seismic reflection data, additional detailed geologic mapping, and well cuttings analysis, a 3-dimensional geologic model of the Astor Pass geothermal system was constructed. The 3D model indicates a complex interaction/intersection area of three discrete fault zones: a northwest-striking dextral-normal fault, a north-northwest-striking normal-dextral fault, and a north-striking west-dipping normal fault. These two discrete, critically-stressed intersection areas plunge moderately to steeply to the NW-NNW and probably act as conduits for upwelling geothermal fluids.

15 GEOTHERMAL ENERGY↗

Three-Dimensional Geothermal Fairway Mapping: Examples From the Western Great Basin, USA

Elevated permeability along fault systems provides pathways for circulation of geothermal fluids. Accurate location of such fluid flow pathways in the subsurface is crucial to future geothermal development in order to both accurately assess resource potential and mitigate drilling costs by increasing drilling success rates. Employing a variety of surface and subsurface data sets, we present detailed 3D geologic analyses of two Great Basin geothermal systems, the actively producing Brady’s geothermal system and a ‘greenfield’ geothermal prospect at Astor Pass, Nevada. 3D modeling provides the framework for quantitative structural analyses. We combine 3D slip and dilation tendency analysis along fault zones and calculations of fault intersection density in the two geothermal systems with the locations of lithologies capable of supporting dense, interconnected fracture networks. The collocation of these permeability promoting characteristics with elevated heat represent geothermal ‘fairways’, areas with ideal conditions for geothermal fluid flow. Location of geothermal fairways at high resolution in 3D space can help to mitigate the costs of geothermal exploration by providing discrete drilling targets and data-based evaluations of reservoir potential.

15 GEOTHERMAL ENERGY↗

Great Basin NV Play Fairway Analysis - Carson Sink

All datasets and products specific to the Carson Sink Basin. Includes a packed ArcMap (.mpk), individually zipped shapefiles, and a file geodatabase for the Carson Sink area; a GeoSoft Oasis montaj project containing GM-SYS 2D gravity profiles along the trace of our seismic reflection lines; a 3D model in EarthVision; spreadsheet of links to published maps; and spreadsheets of well data.

15 GEOTHERMAL ENERGY↗