Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum correlations in quantum information”

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.

1,525 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Improving Marine Energy Production Through Commercialization of a Low Cost, Drag-Reducing Slippery Coating (CRADA 679 Abstract)

Marine energy capture systems offer great promise for providing clean energy, but they operate in a challenging and dynamic environment and must be optimized to the highest extent possible. Computational studies predict that drag reduction on marine energy and blue economy systems will result in meaningful improvements in energy efficiency. Based on extensive coating experience, PNNL is proposing to bring to market a new drag-reducing coating called Superhydrophobic Lubricant-Infused Drag-Efficient Coating – SLIDE-Coat. PNNL has a deep knowledge base regarding this class of slippery coatings, and we have a group of enthusiastic industry partners who have committed to partnering, conducting field testing, and providing well over 50% cost share. In addition to validating the technology, the team will create a commercialization roadmap to ensure the commercial success of the technology after the government-sponsored two-year program is complete. The main technical goal is to design, manufacture, and experimentally validate a new coating that can reduce hydrodynamic skin friction drag by 20%. The main commercialization goal is to develop a roadmap that identifies activities needed to complete SLIDE-Coat’s development after the conclusion of this project. Year 1 will focus on adaptation and modification of the existing SLIC coating system to optimize drag reduction and assessment and quantification of drag reduction on relevant materials in a laboratory setting. Year 2 will focus on optimizing and demonstrating drag reduction with preferred coatings. The manufacturability, ease of application, adhesion to relevant surfaces, and consistency will be assessed. Quantification of drag reduction on prototype materials in a relevant marine field setting will be executed.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Dynamic Interfacial Design in Adaptive Hybrid Materials Enables Reversible and Tunable Mechano-Optic Smart Responses

Next-generation polymeric materials are shifting toward adaptive and interactive behaviors of living systems; however, designing materials that can reversibly modulate optical properties under mechanical deformation while maintaining mechanical robustness remains a key challenge. Here, we report a mechanically robust vitrimer-based adaptive hybrid material (AHM) that exhibits a stretch-induced reversible transparency-to-opacity transition, enabled by the integration of dynamic interactions at the polymer–silica nanoparticle interface and controlled nanoparticle self-assembly. The AHM combines boronic ester–functionalized polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene (S-Bpin) with diol-functionalized silica nanoparticles (diol-SiNPs) to form a hybrid network hosting both dynamic boronic ester and hydrogen-bonding interactions. These reversible linkages facilitate controlled nanoparticle self-assembly and enable strain-induced nanoparticle alignment/aggregation. Upon stretching, SiNP-rich domains align and aggregate within the polymer matrix, while local modulus mismatch between stiff aggregated SiNP/borylated-styrene-rich regions and the softer elastomeric midblock induces surface microwrinkle formation. These internal aggregates and surface wrinkles cooperatively enhance light scattering, producing the opaque state under strain. Furthermore, the tailored AHM exhibits high toughness, thermomechanical stability, reprocessability, and programmable shape-memory behavior. This work presents a dynamic interfacial design strategy for mechanically robust, optically reconfigurable, and reusable soft materials for adaptive optics, smart windows, sensing, soft robotics, and circular smart-material platforms.

adaptive hybrid materials↗
DOE OSTI2024

An Overview of Electric Vehicle Load Modeling Strategies for Grid Integration Studies

The adoption of electric vehicles (EVs) has emerged as a solution to reduce greenhouse gas emissions in the transportation sector, which has motivated the implementation of public policies to promote their use in several countries. However, the high adoption of EVs poses challenges for the electricity sector, as it would imply an increase in energy demand and possible impacts on the power quality (PQ) of the power grid. Therefore, it is important to conduct EV integration studies in the power grid to determine the amount that can be incorporated without causing problems and identify the areas of the power sector that will require reinforcements. Accurate EV load patterns are required for this type of study that, through mathematical modeling, reflect both the dynamic behavior and the factors that influence the decision to recharge EVs. This article aims to present an overview of EVs, examine the different factors considered in the literature for modeling EV load patterns, and review modeling methods. EV load modeling methods are classified into deterministic, statistical, and machine learning. The article shows that each modeling method has its advantages, disadvantages, and data requirements, ranging from simple load modeling to more accurate models requiring large datasets.

Computer Science↗
DOE OSTI2025

Mild-Annealed Molecular Layer Deposition (MLD) Tincone Thin Film as Photoelectrochemically Stable and Efficient Electron Transport Layer for Si Photocathodes

Metalcone thin films, composed of inorganic–organic hybrids, are synthesized using molecular layer deposition (MLD) through reactions between organometallic precursors (e.g., Sn, Al, and Ti) and organic reactants (e.g., ethylene glycol and glycerol). Despite their unique properties, metalcones exhibit significant vulnerability to water due to their organic components, limiting their potential in electrochemical applications. This study focuses on enhancing the photoelectrochemical stability of tincone thin films in aqueous electrolyte while preserving their hybrid characteristics through mild annealing in air at 250 °C. As-deposited and vacuum-annealed tincone thin films exhibited significant degradation under these conditions, while high-temperature-annealed (500 °C) tincone thin films offered improved stability with a significant decline in charge transfer efficiency. In contrast, mild annealing in air maintained the C–O bond at half level and improved the stability and charge transport without compromising the unique characteristics of tincone. This was confirmed by ellipsometry, X-ray photoelectron spectroscopy (XPS), and Fourier transform infrared spectroscopy (FTIR). Mild-annealed tincone deposited on a lightly doped p-type silicon (p-Si) photocathode produced a 20-fold increase in CO volume compared to high-temperature annealed tincone in a CO 2 -saturated potassium bicarbonate (KHCO 3 ) electrolyte with dispersed graphene oxide–cobalt phthalocyanine (GO-CoPc) under 1 sun illumination at 0.9 V vs reversible hydrogen electrode (RHE), while maintaining the faradaic efficiency for CO and H 2 . These results suggest that mild-annealed tincone thin films hold significant potential as protective charge transport layers on silicon photocathodes for the aqueous CO 2 reduction reaction (CO 2 RR).

Annealing (metallurgy)↗
DOE OSTI2023

Polymer grafted aramid nanofiber reinforces immiscible waste polypropylene/poly(ethylene terephthalate)

Polypropylene (PP) and poly(ethylene terephthalate) (PET) are plastics commonly used for packaging because of their excellent barrier and mechanical properties. The properties of these plastics are often diminished after mechanical recycling, inevitably causing down-cycling. Furthermore, this problem is exacerbated when different kinds of polymers mix. Aramid nanofibers have the potential to improve the mechanical properties of polymers due to their excellent mechanical properties but their poor dispersion in polymers is a challenge. Grafting polymers onto nanofibers can help address this challenge. In this work, different loading levels (1%, 2%, and 5%) of polymer grafted aramid nanofibers (ANFs) are blended with waste PP/PET (90/10), simulating a PP waste stream containing traces of PET contaminants. Scanning electronic microscopy, rheology, and differential scanning calorimetry results show the affinity of PP functionalized aramid nanofibers (PP_ANF) towards the PP matrix. At 1 wt% of the nanofiber, the size of the PET droplets in the PP matrix of the PP_ANF blend range from 0.2 to 2.0 μm while that of unmodified ANF and PET_ANF blends are in the range of 0.1–6.2 and 0.5–7.4 μm, respectively. In summary, polymer grafted ANFs have the tendency of improving properties of its like polymers due to similarity in the grafting polymer and the polymer matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

A Sensitivity-driven Wide Area Protection (SWAP) Coordination Tool for High Penetration of Inverter-based Resources (IBR)

Traditionally, power system generation sources have been composed of synchronous generators, of which the fault current behavior is understood with minimal differences between generation size and types due to the physics of their construction. Present protection schemes and modeling methods are based upon these understood characteristics. Most renewable generation is composed of inverter-based resources (IBR), in which fault current is determined by switching control software and hardware limitations, each of which can vary between manufacturers and even between models of the same manufacturer. The resulting fault current is low in magnitude, low in negative-sequence current, unpredictable phase angles, and is a challenge to model. These characteristics also result in a challenge to traditional protection schemes and fault simulation software. To address several of these concerns, the project has the following goals: 1. Improve IBR models: Improve IBR models used in short circuit (SC) programs to accurately capture the response of IBRs at the bulk power system (BPS) level for fault and protection studies. 2. Develop automation tool: Develop an automation tool that allows engineers to identify protection coordination and sensitivity issues by performing SC and protection coordination studies in a high IBR-penetrated grid by applying variations to the IBR models, faults, contingencies, etc. 3. Develop schemes: Develop new protection mitigation solution schemes that complement the existing protection systems to ensure safe operation of the BPS with higher IBR penetration levels. The project team did not achieve this final goal, as the Department of Energy (DOE) stopped the project early due to changes in DOE funding priorities. The termination notice came at the beginning of the final project phase, while the team was identifying and beginning to investigate protection issues. It should be noted that the team discussed a 100% penetration scenario. However, this scenario would require the use of grid-forming IBR models that are not presently available. Since developing these models requires additional effort, the 100% penetration scenario was not pursued during this project. In the future, developing the methodology and models for the 100% scenario could benefit the industry.

14 SOLAR ENERGY↗
DOE OSTI2025

Decoupling thermal and irradiation effects on grain boundary segregation

Radiation-induced segregation (RIS) is most often measured by peak solute concentration at a boundary. However, this may give an incomplete picture of segregation quantity and phenomena. Radiation-induced and thermal segregation at grain boundaries was investigated in Fe-9.6 at.% Cr after 9 MeV Fe 3+ ion irradiation at 400 °C. The experimental results were compared to kinetic Monte Carlo (KMC) simulations. The study revealed that Cr enrichment (peak segregation) at the grain boundaries was comparable in both the irradiated and non-irradiated conditions, although irradiation resulted in broader segregation profiles in both experiments and simulations, indicating an overall increase in grain boundary segregation due to irradiation. This broadening is attributed to back diffusion into the grain interior. While it is an established phenomenon, this study offers a quantitative evaluation using experimental data and KMC modeling. Further, these results emphasize the importance of analyzing the entire segregation profile and decoupling the thermal and irradiation contributions to solute segregation.

Grain boundary segregation↗
DOE OSTI2026

Divergent trait controls on soluble sugars and starch underlie global strategies of tree carbohydrate storage

Nonstructural carbohydrate (NSC) stores buffer tree metabolism, osmotic regulation, and defense, thereby mediating tolerance and survival under climate extremes. Yet, the functional and evolutionary determinants of interspecific variation in NSC remain elusive, limiting understanding and prediction of forest carbon allocation and mortality under global change. Here, we present a cross-species synthesis of NSC concentrations across multiple organs for 281 woody species from 102 mixed forest communities worldwide, where we quantified species-specific deviations from community means to disentangle intrinsic trait effects from environmental and methodological variation. We found phylogenetic signals in NSC deviations, with coniferous gymnosperms and evergreen species consistently maintaining lower stem soluble sugars and starch concentrations than co-occurring angiosperms and deciduous species, respectively. A global pattern emerged where greater stomatal sensitivity to leaf water potential was associated with declines in the relative concentrations of both sugars and starch. In contrast, xylem hydraulic safety traits showed weak and organ-dependent relationships with NSC concentrations. Sugars increased with photosynthetic capacity and declined with wood density, whereas starch showed the reverse pattern, which aligned with the distinct functional-metabolic roles of sugars and starch. By integrating trait-based ecology with a community-centered framework, our study provides global evidence that stomatal regulation, photosynthetic capacity, specific leaf area, and wood density jointly govern interspecific NSC variation, through contrasting effects on sugars and starch. These are among the most broadly measured traits globally, thus the emergent carbohydrate–trait relationships can have broad applications toward understanding and predicting forest growth and survival under climate change.

tropic system↗
DOE OSTI2026

Elimination of detrimental grain boundary segregation in garnets

Garnet Li 7 La 3 Zr 2 O 12 electrolyte is considered a key enabler of solid-state batteries with Li metal electrodes, but the grain boundaries impair its performance. To date, the understanding of grain boundary structures and its impact on performance remains elusive. Here, we show that element segregation at Li 7 La 3 Zr 2 O 12 grain boundaries critically governs Li transport and nucleation. During conventional sintering, Al, Ta, and La segregate at grain boundaries, locally depleting Li and creating space-charge layers that lower total ionic conductivity. Simultaneously, this segregation leads to higher electronic conductivity along grain boundaries, which promotes Li nucleation at grain boundary edges with increased risk of dendrite formation. The underlying mechanism of segregation is governed by both thermodynamic driving forces and diffusion kinetics. Building on this understanding, we develop a strategy to achieve segregation-free grain boundaries through a rapid sintering protocol that utilizes the onset of solid-state softening. This approach yields transparent, polycrystalline Li 7 La 3 Zr 2 O 12 with negligible grain boundary impedance and enhanced dendrite tolerance. By elucidating the structural origins and electrochemical consequences of grain boundary segregation, this work provides a guidance for the rational optimization of solid electrolytes.

Energy - Storage↗
DOE OSTI2024

Insights into the thermo-hydraulic properties of compacted MX80 bentonite during hydration under elevated temperature

In high-level radioactive waste geological repositories, compacted bentonite undergoes coupled thermo-hydraulic processes due to heat released from a central waste canister and groundwater imbibition from the surrounding host rock. An understanding of these processes is essential for long-term simulations of radionuclide migration and canister corrosion, which requires an understanding of temperature effects on the coupled thermo-hydraulic properties governing these processes. In this study, a tank-scale radial infiltration test was used to investigate water imbibition processes in compacted bentonite under a central heater temperature of 200 °C that simulates high thermal gradients in a repository. Interpretation of this test focuses on evaluation of the liquid water wetting front during hydration and the interpretation of the temperature-dependent transient soil water retention curve (SWRC), thermal conductivity function (TCF), and hydraulic conductivity function (HCF). The SWRC during imbibition follows a temperature-dependent wetting path. Temperature effects on the hydraulic conductivity of bentonite in saturated conditions had the greatest effect on the shape of the HCF, with minimal temperature effects at higher suctions. The transient thermal conductivity data matched well with a new TCF linked with the SWRC shape. Although the bentonite layer was restrained, local deformations during hydration may have affected the shapes of the TCF and HCF.

Engineering↗
DOE OSTI2026

Simple device modification simultaneously enhances indoor passive evaporative cooling power and duration

Passive evaporative cooling from a damp cloth or vessel can be used for localized cooling of perishables such as food and medicine, drinking water, and even people. Such applications have two key objectives, which are normally in tension with each other: minimizing the water consumption rate and maximizing the cooling (steady-state temperature swing between the cool cloth and the warm surroundings). Here, we present a simple device that delivers higher performance in both metrics simultaneously by adding a perforated aluminum foil sheet suspended over a stagnant air layer to thermally shield the evaporating surface while also facilitating mass transfer. Experimental results demonstrate a simultaneous 10%–25% increase in temperature swing and 2–3× reduction in water consumption rate as compared to conventional evaporative cooling. This improvement is achieved through careful modeling to optimize the coupled heat and mass transfer through the airgap (thermally insulating and vapor permeable) and perforated foil (infrared reflective and vapor permeable).

Chen, Sarah Mizuno↗
DOE OSTI2026

Metrology for femtosecond pulsed x-ray heating in diamond anvil cell experiments at the European XFEL: Revisiting the iron phase diagram up to 150 GPa

The development of pulsed intense x-ray sources, such as free electron laser, offers new avenues for high pressure experiments. Here, we study the feasibility and metrology of x-ray heating in diamond anvil cells at the European x-ray free electron laser. This method enables one to volumetrically heat the sample while inhibiting chemical migration and probing the crystallographic structure of the sample throughout the heating with a high repetition rate. We focus our study on iron, whose phase diagram is well established up to 100 GPa, to explore the possibilities and limitations of this technique. We volumetrically heat iron samples at starting pressures ranging from 10 to 138 GPa, using the x-ray beam pulsed at 4.5 MHz in a serial pump-and-probe experimental design. Experimental challenges arise from temperature gradients within the sample, changes in temperature at the 100 ns timescale, the difficulty of direct temperature estimates, the effect of thermal pressure, and the presence of metastable crystallites due to rapid cycles of heating and cooling. Hence, we develop a multi-crystal-like data processing method that allows us to account for sample heterogeneity in probed conditions. We then calibrate our measurements using known physical properties of iron under pressure. Thermal pressure in our experiments increases from 4% of the isochoric prediction at 10 GPa to 23% at 138 GPa, and we show that our data are in agreement with most previous observations of iron in this pressure range. The method can now be implemented at higher pressures and temperatures and on materials with unknown phase diagrams.

Materials science↗
DOE OSTI2025

Active Fluids Form System-Spanning Filamentary Networks

Recent experimental realizations of liquid-liquid phase separation of active liquid crystals have offered an insight into the interaction between phase separation, ubiquitous in soft matter and biology, and chaotic active flows. Here, in this Letter, we use continuum theory to examine phase separation of an active liquid crystal and a passive fluid and report two new results. First, we provide an analytical derivation of the activity-induced suppression of the phase boundary of the coexistence region—a result first reported in simulations and experiments. We show that the shift in the critical point is a result of the balance between self-stirring active flows and phase-separating diffusive fluxes. Second, we show that this same balance is responsible for dramatically changing the morphology of the phase separated state, resulting in the emergence of a new mixed active phase consisting of a dynamical filamentous active network that invades the entire system area, trapping droplets of passive material. This structure exists even for very low volume fractions of active material. Our work provides an important step towards the goal of understanding how to use activity as a new handle for sculpting interfaces.

active nematics↗
DOE OSTI2024

Prediction of BiS2-type pnictogen dichalcogenide monolayers for optoelectronics

Abstract In this work, we introduce a 2D materials family with chemical formula MX 2 (M={As, Sb, Bi} and X={S, Se, Te}) having a rectangular 2D lattice. This materials family has been predicted by systematic ab-initio structure search calculations in two dimensions. Using density-functional theory and many-body perturbation theory, we study the structural, vibrational, electronic, optical, and excitonic properties of the predicted MX 2 family. Our calculations reveal that the predicted SbX 2 and BiX 2 monolayers are stable while the AsX 2 layers exhibit an in-plane ferroelectric instability. All materials display strong excitonic effects and good optical absorption within the infrared-to-visible range. Hence, these monolayers can harvest solar energy and serve in optoelectronics applications. Furthermore, our results indicate that exfoliation of the predicted MX 2 monolayers from their bulk counterparts is experimentally viable.

Materials Science↗
DOE OSTI2026

Solar Heat for Industrial Processes: Integration with Chemical Reactors

The integration of solar thermal systems with chemical reactors has been proposed as part of a larger effort to develop and deploy solar heat for industrial processes (SHIP) technologies. A strong motivation for SHIP processes and technologies is the potential for high thermal efficiency coupled with low-cost thermal energy storage (TES) which can enable commercial deployment of such systems. While there are different ways to categorize SHIP technologies, one important such distinction is between directly irradiated systems and indirect off-sun process driven by a SHIP system. While directly irradiated systems can provide high thermal efficiencies and high fluxes, they usually require complex engineering solutions due to the need for redesigning the established processes and unit operations. In most cases, it is also more challenging to couple such a process to a TES system, losing some of the benefits of SHIP. On the other hand, using a SHIP system to drive an industrial process off-sun can allow better integration with existing process chains, easier TES capabilities, and potential for more applications fitting a specific SHIP technology. However, the integration of SHIP systems with the industrial processes is not fully explored in detail, especially in the case of high-temperature processes such as reforming, cracking, cement manufacturing, and iron/steelmaking. Many of these systems require heating fluxes of >50 kW/m^2, supplied via combustion of hydrocarbons in a fire box and benefitting from radiative heat transfer between the flue gases and the reaction zones. As such, using SHIP systems for such processes is more complex than providing the same thermal input in the form of a heat transfer medium (HTM) entering the reactor, kiln, or furnace. Moreover, in case convective heat transfer using SHIP is envisioned, for example using supercritical CO2 as the HTM from a particle receiver, the thermal integration might be more challenging than initially envisioned: lower heat transfer coefficients and limited approach temperature might require large flow rates, causing a mismatch between the process thermal requirements and the thermal capacity of the SHIP system. In addition, even if the heat exchange between SHIP and reactor is effective, there is still a cold leg HTM at the reaction temperature or slightly below it. Chemical plants usually include a set of heat exchangers, heat recovery steam generators, and even power generation units - in a tightly integrated design - to recover the flue gases which are eventually vented. With SHIP systems mostly operating on a closed HTM loop, bottoming the cold leg is crucial. In this talk, we will present different modeling results for a variety of syngas production reactions, using catalytic and chemical looping processes, and discuss some of the challenges and design considerations for off-sun chemical reactors using SHIP systems.

14 SOLAR ENERGY↗
DOE OSTI2026

Virtual Growth of SRF Materials

Niobium's native surface oxide affects SRF cavity and superconducting qubit performance, motivating interest in controlling its crystalline structure. We combine a literature-derived machine-learning analysis with temperature-dependent XRD to study crystalline ordering in Nb2O5. Random Forest models, trained on 74 processing conditions from 17 papers and validated by leave-one-group-out cross-validation, predicted broad crystallinity outcomes well (balanced accuracy 0.809), but struggled with specific polymorph identity (0.577). Annealing temperature was the dominant predictor across all targets; oxygen partial pressure showed negligible importance, reflecting narrow literature coverage rather than physical irrelevance. Temperature-dependent XRD on anodized and H2O2-treated Niobium showed structural evolution consistent with the machine learning predictions. Our model and overall approach provide a data-driven framework for identifying and optimizing conditions that promote crystallization in initially amorphous oxides. This framework can guide the selection of growth and post-annealing conditions for Nb surfaces by narrowing the experimental parameter space, thereby reducing trial-and-error efforts in developing oxide structures relevant to SRF applications.

Tilkin, Anthony [Fermilab]↗
DOE OSTI2026

Value of Geothermal Energy Storage for Supply-Side and Demand-Side Applications

This report presents the results of a study examining the value potential for geothermal energy storage (GES), a long-duration energy storage resource that stores thermal and/or geomechanical energy in the subsurface. GES could benefit the overall U.S. power system by temporally shifting electricity generation (supply-side) or meeting building heating and cooling load (demand-side). This report analyzes supply-side and demand-side opportunities independently because of differences in applications and models. Currently there is significant uncertainty about the development costs for GES, with only a limited number of demonstration plants for electric energy storage and building heating and cooling storage developments. In this report, we estimate the value of supply-side and demand-side GES to the bulk power system in the contiguous United States. Because of the significant uncertainty about GES development costs, this analysis does not consider GES deployment costs but instead focuses on the value of GES to the U.S. electricity system. The estimated values of GES provide reference points for economically competitive commercial cost targets. Supply-side GES is modeled as part of an enhanced geothermal system (EGS) generation plant in NREL's Regional Energy Deployment System (ReEDS) capacity expansion model (Ho et al. 2021). In contrast to conventional geothermal plants, which generate constant power, EGS plants have unique features that may allow for in-reservoir energy storage for flexible generation. Demand-side GES for heating and cooling, including seasonal hot and cold storage and short-duration heat pump storage, is incorporated into a price-taker model using Cambium electricity marginal cost projections. To establish an upper bound for the value of GES, analysis focused on favorable scenarios for storage with high generation from zero marginal cost, variable renewable energy resources. High penetrations of variable renewable energy generation can increase hourly electricity price variability, which increases the value of temporal energy arbitrage for storage technologies like GES.

15 GEOTHERMAL ENERGY↗