Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Condition Monitoring”

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

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

WIRE-21 Sensor Irradiation Experiment Ready for HFIR Insertion

The ability to deploy new nuclear fuels for current or future reactor concepts requires a wealth of data regarding fuel performance during normal operation, anticipated operational occurrences, and design-basis accidents. Most of these data have historically been collected during experiments in materials test reactors, ideally with online instrumentation to collect as much data as possible. However, advanced instrumentation could also allow for in situ monitoring of fuel operating conditions during commercial reactor operation to maximize fuel utilization, reduce unnecessary conservativism in design margins, and improve operator understanding of limiting peaking factors. The latter approach would complicate fuel handling, particularly during refueling, unless the instrumentation could be placed inside the fuel rods and transmitted wirelessly to a receiver located outside the fuel’s primary pressure boundary. To this end, Westinghouse Electric Company (WEC) developed wireless sensors based on inductive coupling that can transmit information regarding fuel centerline temperatures and rod internal pressures wirelessly from within a fuel rod to a nearby instrument thimble. After testing these sensors in lower-power university research reactors, the next step is to perform high neutron fluence testing to characterize the performance of these wireless sensors under conditions that are more representative of the intended application—in this case, light-water reactors (LWRs). The removable Be (RB) positions of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL) provide the neutron flux, experiment volume, and access to instrument leads required to achieve these sensor testing goals. This report summarizes the design, analysis, and assembly of the Wireless Instrumented RB Experiment 2021 (WIRE-21). This is the most highly instrumented irradiation experiment ever performed in HFIR. The experiment will use seven different sensing techniques to measure temperature, pressure, neutron flux, and neutron fluence during reactor operation. In addition to WEC’s wireless temperature and pressure sensors, WIRE-21 includes an array of thermocouples, self-powered neutron detectors, spatially distributed fiber optic temperature sensors, passive SiC temperature monitors, and flux wires. The design of WIRE-21 and the cabling that was installed in HFIR also provide the infrastructure to enable accelerated, economical testing of advanced sensor technologies while leveraging the extremely high neutron flux that is available in HFIR. The containment for WIRE-21 is similar to previous RB irradiation vehicles but includes a few modifications, most notably the use of integrated compression seals to pass a larger number of sensor leads through the experiment’s pressure boundary. In addition to the sensor leads, inert gas lines are passed into the experiment to enable active temperature control and the ability to pneumatically actuate a bellows-driven pressure sensor. WIRE-21 is targeting component temperatures (300–350°C) and neutron fluence levels (~10 22 n/cm 2 ) that would be expected in the plenum region of LWR fuels, except for the active sensing region of the wireless temperature sensor, which is targeting LWR fuel centerline temperatures (~800–1,100°C). WIRE-21 was successfully assembled, passed all nondestructive examination, and was delivered to HFIR for insertion during upcoming cycle 498 (April 2022).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Innovations in Sensors and Controls for Building Energy Management: Research and Development Opportunities Report for Emerging Technologies

Sensors, actuators, and controllers, which collectively serve as the backbone of cyberphysical systems for building energy management, are one of the core technical areas of investment for achieving the U.S. Department of Energy (DOE) Building Technologies Office's (BTO's) goals for energy affordability in the national building stock - both commercial and residential. In fact, an aggregated annual energy savings of 29% is estimated in the commercial sector alone through the implementation of efficiency measures using current state-of-the-art sensors and controls to retune buildings by optimizing programmable settings based on occupant schedules and comfort requirements, as well as detecting and diagnosing equipment operation and installation problems (Fernandez et al. 2017). Monitoring and control of building conditions and operations has advanced significantly, from the invention of the modern thermostat just before the start of the 20th century to the midcentury incorporation of direct digital control into devices, the introduction of open protocols and network communications at the end of the century, and finally the invention of cloud-based computing and additional advancements that have enabled remote operation and a proliferation of connected and intelligent devices in building automation. Despite this potential, however, two main challenges hinder widespread adoption of sensors and controls in building operations that can ensure savings for high-efficiency components and equipment (e.g., heat pumps, windows, and lighting devices), as well as additional savings from more sophisticated control architectures and algorithms. energy savings of 29% is estimated in the commercial sector alone through the implementation of efficiency measures using current state-of-the-art sensors and controls to retune buildings by optimizing programmable settings based on occupant schedules and comfort requirements, as well as detecting and diagnosing equipment operation and installation problems (Fernandez et al. 2017). Monitoring and control of building conditions and operations has advanced significantly, from the invention of the modern thermostat just before the start of the 20th century to the midcentury incorporation of direct digital control into devices, the introduction of open protocols and network communications at the end of the century, and finally the invention of cloud-based computing and additional advancements that have enabled remote operation and a proliferation of connected and intelligent devices in building automation. Despite this potential, however, two main challenges hinder widespread adoption of sensors and controls in building operations that can ensure savings for high-efficiency components and equipment (e.g., heat pumps, windows, and lighting devices), as well as additional savings from more sophisticated control architectures and algorithms.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Annual Site Environmental Report (2019)

Los Alamos National Laboratory’s (the Laboratory’s) annual site environmental reports are prepared by the Laboratory’s environmental organizations, as required by U.S. Department of Energy Order 231.1B, Administrative Change 1, Environment, Safety, and Health Reporting, and Order 458.1, Administrative Change 3, Radiation Protection of the Public and the Environment. The following chapters in this report discuss our success in complying with environmental laws, regulations, and orders (Chapter 2, Compliance Summary); how we manage the Laboratory’s environmental performance (Chapter 3, Environmental Programs); how we monitor for air emissions of radioactive materials and climate conditions (Chapter 4, Air Quality); how we monitor for effects of Laboratory operations on groundwater quality (Chapter 5, Groundwater Protection); how we monitor the movement of chemicals and radionuclides by storm water runoff and the levels of chemicals and radionuclides in deposited sediment (Chapter 6, Watershed Quality); how we monitor for the presence, levels, and effects of chemicals and radionuclides in plants, animals, and soil (Chapter 7, Ecosystem Health); and finally, what radionuclide dose or risk from chemical exposure members of the public may experience as a result of Laboratory operations (Chapter 8, Public Dose and Risk Assessment).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

ASER Annual Site Environmental Report 2020

Los Alamos National Laboratory’s (the Laboratory’s) annual site environmental reports are prepared each year by the Laboratory’s environmental organizations, as required by U.S. Department of Energy Order 231.1B, Administrative Change 1, Environment, Safety, and Health Reporting, and Order 458.1, Administrative Change 4, Radiation Protection of the Public and the Environment. The following chapters in this report discuss our success in complying with environmental laws, regulations, and orders (Chapter 2, Compliance Summary); how we manage the Laboratory’s environmental performance (Chapter 3, Environmental Programs and Analytical Data Quality); how we monitor for air emissions of radioactive materials and climate conditions (Chapter 4, Air Quality); how we monitor for effects of Laboratory operations on groundwater quality (Chapter 5, Groundwater Protection); how we monitor the movement of chemicals and radionuclides by storm water runoff and the levels of chemicals and radionuclides in deposited sediment (Chapter 6, Watershed Quality); how we monitor for the presence, levels, and effects of chemicals and radionuclides in plants, animals, soil, and vegetation (Chapter 7, Ecosystem Health); and finally, what radionuclide dose or risk from chemical exposure members of the public may experience as a result of Laboratory operations (Chapter 8, Public Dose and Risk Assessment).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Deep Vadose Zone Monitoring Test Bed (FY23 Status Report)

A subsurface air injection at the Hanford Site’s Deep Vadose Zone Monitoring Test Bed was completed to realize a change in subsurface hydrologic conditions in accordance with a soil desiccation remedy. The injection mimicked a previous injection at the site that was performed in accordance with a vadose zone treatability study. Unlike the previous test which relied on electrical methods only, the change in hydrologic conditions during the recent test was also monitored using cross-hole seismic sensing methods to assess the ability of the seismic methods to evaluate changes in moisture conditions of desiccated sediment at the field scale. Data from in situ neutron probes indicates a reduction of soil moisture in the vicinity of the injection well due to the air injection. Similar changes were observed in the time-lapse electrical resistivity and seismic data, which indicates a loss of soil moisture over time. Tomographic inversions of the time-lapse geophysical data illustrate the 2D and 3D features of the soil moisture distribution over time. Time-lapse electrical resistivity tomography (ERT) results show a reduction in the electrical conductivity of the subsurface in the vicinity of the injection well, with most changes occurring within the screened interval. Similar patterns are observed in the seismic tomography results, with both methods illustrating two lobe-shaped features of reduced soil moisture. Use of seismic and ERT technologies in tandem takes advantage of two complementary geophysical monitoring technologies, providing increased sensitivity to specific hydrologic conditions. The multiphysics approach, therefore, has the potential to improve the ability to estimate subsurface moisture conditions from sensor-based and remotely sensed geophysical data that will ultimately improve the ability of remediation contractors to evaluate remedy performance.

58 GEOSCIENCES↗
DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 2. Quantitative Moisture Dynamics Estimation Model

Abstract The long‐term containment of high‐level radioactive waste in geological disposal repositories relies on Engineered Barrier Systems (EBS), with bentonite clay emerging as a candidate material due to its unique properties. Understanding moisture dynamics within bentonite buffers is crucial for EBS performance, as it directly influences the material's swelling capacity, thermal and hydraulic conductivity, mechanical properties, and long‐term evolution under complex thermal‐hydrological‐mechanical (THM) processes. This study develops an advanced Electrical Resistivity Tomography (ERT)‐based framework to quantitatively monitor moisture dynamics under THM conditions. Our framework extends the Waxman‐Smits model to incorporate the coupled effects of temperature, water content, fluid chemistry, and mechanical changes on bentonite's electrical properties. Utilizing HotBENT‐Lab data from our companion paper, which includes electrical conductivity, CT density, and thermocouple measurements, this study offers a novel methodological framework bridging different scales of the model. Our results show that the extended model can estimate water content from ERT data, capturing spatial and temporal variations in moisture distribution within bentonite columns. However, the model tends to overestimate water content compared to CT density‐derived measurements. We address this discrepancy by incorporating a simplified swelling effect model, which improves agreement between ERT and CT density‐based water content estimates. We also discuss model limitations, including simplified treatment of swelling and micropore effects, and propose a conceptual framework for transitioning from laboratory to field applications, addressing challenges such as parameter scalability, field validation methods, and integration of diverse data sources. This ERT‐based framework can potentially advance real‐world moisture monitoring of bentonite‐based EBS in nuclear waste repositories. Plain Language Summary Safely containing high‐level radioactive waste depends on barriers made from materials like bentonite clay, which is effective because it swells and seals in the waste. To ensure these barriers work well over time, it's important to understand how moisture moves through the clay. Our study developed a new method using ERT to monitor moisture levels in bentonite under conditions that mimic those in actual storage sites, including changes in temperature, water content, and mechanical stress. This study improved an existing model to better account for how these factors affect the clay, allowing us to create more accurate moisture maps. Initially, the proposed model overestimated the amount of water in the clay, but its accuracy was improved by factoring in how the clay swells when wet. This study also identified some limitations of the model and suggested ways to adapt it for use in real‐world waste storage sites. This new approach could lead to better monitoring and safety checks for nuclear waste storage systems, helping to ensure long‐term containment. Key Points This work develops an ERT‐based framework extending the Waxman‐Smits model to monitor bentonite moisture dynamics during coupled THM processes The extended model accurately estimates water content from Electrical Resistivity Tomography data, incorporating swelling effects to improve precision This work proposes a conceptual framework for transitioning from laboratory to field applications, advancing EBS monitoring in nuclear waste repositories

Chen, Hang↗
DOE OSTI2020

Spectroradiometric detection of competitor diatoms and the grazer Poteriochromonas in algal cultures

To address challenges in early detection of pond pests, we have extended a spectroradiometric monitoring method, initially demonstrated for measurement of pigment optical activity and biomass, to the detection of algal competitors and grazers. The method relies upon measurement and interpretation of pond reflectance spectra spanning from the visible into the near-infrared. Reflectance spectra are acquired every 5 min with a multi-channel, fiber-coupled spectroradiometer, providing monitoring of algal pond conditions with high temporal frequency. The spectra are interpreted via numerical inversion of a reflectance model, in which the above-water reflectance is expressed in terms of the absorption and backscatter coefficients of the cultured species, with additional terms accounting for the pigment fluorescence features and for the water-surface reflection of sunlight and skylight. With this method we demonstrate detection of diatoms and the predator Poteriochromonas in outdoor cultures of Nannochloropsis oceanica and Chlorella vulgaris, respectively. Further, the relative strength of these signatures is compared to microscopy and sequencing analysis. Spectroradiometric detection of diatoms is then further assessed on beaker-contained mixtures of Microchloropsis salina with Phaeodactylum tricornutum, Thalassiosira weissflogii, and Thalassiosira pseudonana, respectively, providing an initial evaluation of the sensitivity and specificity of detecting pond competitors.

09 BIOMASS FUELS↗
DOE OSTI2024

Adaptive Narrowband Damping for Improving Harmonic Stability of Modular Multilevel Converter

Harmonic instability events between modular multilevel converter (MMC) and ac systems have been widely reported in recent years. To damp harmonic resonance, this paper proposes an adaptive narrowband damping control that automatically programs, adds, and adjusts the damping around the oscillation frequency when an oscillation is detected. First, the paper presents a low-pass filter design for MMC control loops that pushes all negative damping of the MMC impedance down to the medium-frequency range (< ~ 1000 Hz). Then, an adaptive damping control that uses online oscillation detection is proposed, which can automatically configure the narrowband damper to provide positive damping to the MMC around the detected oscillation frequency. In contrast to existing narrowband damping methods, the proposed adaptive narrowband damper dynamically adjusts the damping gain and the width of the damping range based on continuous monitoring of system resonance conditions (e.g., adjust damping gain to zero when the system resonance disappears). Electromagnetic transient simulation results validate the efficacy of the proposed method in two typical MMC-based power systems.

active damping↗
DOE OSTI2025

WEST L-mode record long pulses guided by predictions using Integrated Modeling

A new record was set on the WEST Tokamak, designed to operate long duration plasmas in a tungsten (W) environment, with an injected energy of 1.15 GJ and a plasma duration 364s. Scenario development was supported by integrated modeling using the High Fidelity Plasma Simulator (HFPS), the European IMAS-coupled version of JETTO/JINTRAC, which integrates physics-driven modules into a unified framework. In particular, a reduced model for Lower-Hybrid heating and Current-Drive (LHCD) and the quasi-linear turbulent transport model TGLF are crucial for long pulses predictions up to the Last Closed Flux Surface (LCFS). Using this workflow, a 100 s reference discharge was modeled and plasma kinetic profiles and loop voltage were quantitatively well matched. In preparation for the recent long duration experiments, non-inductive current-drive actuators (I P , n e , P LHCD ) were varied to determine the operational domain going towards fully non-inductive discharges. In particular, decreasing the plasma current is shown to ease the access to such conditions, with a careful monitoring of (n e , P LHCD ) to avoid machine limitations. In addition, post-prediction experiments conducted within the investigated parameter range validated the predicted dependencies and were shown to be in quantitative agreement. Exploratory work on the use of ECCD for MHD stability purpose is also introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Evaluating factors influencing infrasonic signal detection and automatic processing performance utilizing a regional network

Physical and deployment factors that influence infrasound signal detection and assess automatic detection performance for a regional infrasound network of arrays in the Western U.S. are explored using signatures of ground truth (GT) explosions (yields). Despite these repeated known sources, published infrasound event bulletins contain few GT events. Arrays are primarily distributed toward the south-southeast and south-southwest at distances between 84 and 458 km of the source with one array offering azimuthal resolution toward the northeast. Events occurred throughout the spring, summer, and fall of 2012 with the majority occurring during the summer months. Depending upon the array, automatic detection, which utilizes the adaptive F-detector successfully, identifies between 14% and 80% of the GT events, whereas a subsequent analyst review increases successful detection to 24%–90%. Combined background noise quantification, atmospheric propagation analyses, and comparison of spectral amplitudes determine the mechanisms that contribute to missed detections across the network. This analysis provides an estimate of detector performance across the network, as well as a qualitative assessment of conditions that impact infrasound monitoring capabilities. Finally, the mechanisms that lead to missed detections at individual arrays contribute to network-level estimates of detection capabilities and provide a basis for deployment decisions for regional infrasound arrays in areas of interest.

58 GEOSCIENCES↗
DOE OSTI2021

Collective Analysis of Alpha Particle Losses Due to Self-Absorption by Mass Loading on Radioactive Particulate Glass Fiber Filters

In this study, we derived a relationship between filter mass loading and the percent loss during analysis using the mass loading data collected from six previous studies of self-absorption. Components of mass loading include particulate dust, radioactive particulates, and filter material. In a research report published in 1984, Higby calculated a minimum burial depth for an alpha particle to be lost due to absorption (100% loss) of about 3.7 mg/cm 2 based on calculations for the range of 239-Pu alpha particles in glass fiber filters. From there, Higby concluded that a correction factor of 0.85 assumes approximately 15% losses in the count rate of both alpha and beta particles. In 2000, Luetzelschwab et al. recommended assuming a 40% loss at a loading of 3.3 mg/cm 2 and a 28% loss for a loading of 2.3 mg/cm 2 which included the frontal face mass of the filter. More recently, the 100% losses due to absorption were reported to be in the 10 mg/cm 2 range. Presented here is a trinomial relationship method of relating percent loss due to self-absorption to filter mass loading, based on data reported by Higby, Luetzelschwab et al., Huang et al., Barnett et al., Smith et al., and Hogue et al. Under normal operating conditions at the stacks monitored by Effluent Management, the mass loading of sample filters averages 0.09 ± 0.12 (2s) mg/cm 2 (excluding negative values and outliers) and ranges from 0 mg/cm 2 to 0.24 mg/cm 2 . Based on current mass loading results for Effluent Management stack sample filters, the forced-zero trinomial relationship method estimated self-absorption losses of less than 5%. Because American National Standards Institute/Health Physics Society N13.1-2011 guidelines indicate a correction factor should be used when the penetration of radioactive material into the collection media or self-absorption of radiation by the material collected would reduce the count rate by more than 5%, it is possible continued application of a correction factor to the Effluent Management stack samples is no longer necessary. Nevertheless, continuing to assign a correction factor at the 5% threshold (i.e., 0.95) would be a conservative approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

RCT Continuing Training (Q3, 2021) [Slides]

Radiological Control Technicians (RCTs) serve an important part in the mission of Los Alamos National Laboratory by monitoring the diverse radiological conditions at the Lab. The role of RCTs is critical in maintaining personnel safety from radiological hazards and in keeping occupational radiation exposure As low As Reasonably Achievable (ALARA). This training content consists of two parts: viewing an online lecture and completing the associated exercise guide. The online lecture, presented here, will review performing job coverage and all of the necessary steps associated with it.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

RCT Continuing Training Presentation Transcript (Q3, 2021)

Radiological Control Technicians (RCTs) serve an important part in the mission of Los Alamos National Laboratory by monitoring the diverse radiological conditions at the Lab. This role is critical in maintaining personnel safety from radiological hazards, and in keeping occupational radiation exposure As Low As Reasonably Achievable (ALARA). This training presentation will review performing job coverage and all of the necessary steps associated with it.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Gaussian process for calibration and control of GlueX Central Drift Chamber

We have developed and implemented a machine learning based system to calibrate and control the GlueX Central Drift Chamber at Jefferson Lab, VA, in near real-time. The system monitors environmental and experimental conditions during data taking and uses those as inputs to a Gaussian process (GP) with learned prior. The GP predicts calibration constants in order to recommend a high voltage (HV) setting for the detector that maintains consistent detector performance (gain and resolution) throughout data taking. This approach is in stark contrast to traditional detector operations in which the detector operates at fixed HV and its calibration parameters vary quite considerably with time. Additionally, the ML based system utilizes uncertainty quantification to correct the recommended control parameters when appropriate. We will present results from the ML system autonomously during the Charged Pion Polarizability (CPP) experiment conducted in Hall D at Jefferson Lab.

McSpadden, Helen↗
DOE OSTI2024

Ensemble Manufacturing Techniques for Steam Turbine Components Across Length Scales

Faster design to manufacturing requirements were sought for steam turbine components to meet the changing demands of today’s power grid of improved efficiency through operating temperature increases and enhanced operational flexibility from baseload to cyclic operations. Advances in multiple advanced manufacturing (AM) process enabled redesign of turbine components for extreme environments with the potential to reduce cost. AM is of particular interest to improve component functionality, higher temperature capability, and superior durability in large scale steam turbine applications. AM methods have an immense potential to open-up the design space by working directly with the 3D model to produce near-net shapes, enable fast design-manufacturing iterations, thereby significantly reducing product cost and lead-time up to 25 % from current baseline. However, these benefits cannot be fully realized due to the potential for unknown AM processing defects and their resultant effect upon component performance in service. Siemens is partnering with Oak Ridge National Laboratory (ORNL), Electric Power Research Institute (EPRI), and Connecticut Center for Advanced Technology (CCAT) to advance the knowledge of complex process-material interactions for desired microstructures and properties that are closely interlinked to component geometries across different length scales. The proposed program utilized an ensemble of multidisciplinary technologies to accelerate the development of materials, high-throughput experiments for their qualification, and design flexibility/topology optimization for repair/redesign of components to address critical failure mechanisms for improved performance and increased reliability of existing power plant components. The proposed activities, if successfully demonstrated for identified components, will enable paradigm shift in customized manufacturing and accelerated qualification/certification towards increased steam turbine component durability, increased turbine efficiency, and reduced CO 2 emissions in load-following environments compared to today’s technology. Technology maturation is built into the project as successful research will include customized process-component down-selection enabling AM methodologies to be incorporated directly into the existing supply chain. The specific activities of the proposed effort are: 1. Topology optimization of down-selected steam turbine parts that are amenable to additive and hybrid manufacturing for cost/performance improvement. 2. Process-structure-property relationships for five AM processes for steam turbine materials of interest to compare with conventional materials. 3. Perform part/assembly build process using advanced additive/hybrid machine tools followed by quality inspection of the built components for insight into qualification for production scale-up. Steam turbine rig testing of printed components under targeted, well monitored and characterized environmental conditions of for performance comparison of baseline and redesigned components.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analysis of Rig Parameter Data Using Drilling Process Modeling Constraints, Volume 5: Utah FORGE Well 16B(78)-32

Drill rig parameter measurements are routinely used during deep well construction to monitor and guide drilling conditions for improved performance and reduced costs. While insightful into the drilling process, these measurements are of reduced value without a standard to aid in data evaluation and decision making. In the main body of this work (Volume 1), a method is demonstrated whereby rock reduction model constraints are used to interpret drilling response parameters; the method could be applied in real-time to improve decision-making in the field and to further discern technology performance during post-drilling evaluations. Drilling parameters are evaluated using laboratory-validated rock reduction models for predicting the phenomenological response of drag bits (Detournay and Defourny, 1992) in computational algorithms. The method presented has applicability to development of advanced analytics on future geothermal wells using real-time electronic data recording for improved performance and reduced drilling costs. A drilling cost model is also used to show the tradeoff between rate of penetration and bit life and the influence on interval drilling costs. Details of the bit specifications and performance are cataloged in an independent volume, documented under separate cover, for each of the four wells, and include Volume 2: Utah FORGE 16A(78)-32; Volume 3: Utah FORGE 56-32; Volume 4: Utah FORGE 78B-32 and Volume 5: Utah FORGE 16B(78)-32.

15 GEOTHERMAL ENERGY↗