Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal system 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 91 records · Page 5

Supo Thermal Model Development IV

This report describes the continuation of the Computational Fluid Dynamics (CFD) model of the Supo cooling system described in the reports, Supo Thermal Model Development, II, and III. In the previous reports, the solution natural convection of Supo was simulated using a 2D, axisymmetric model that altered the actual geometry. A 3D geometry of Supo will more accurately estimate the natural convection heat transfer coefficient (HTC) of the system since Supo contains several non-symmetric internal components that can only be modeled in 3D. This, in turn, will improve the natural convection HTC results for comparison purposes with past operational data and previous 2D CFD results. Unfortunately, the complexities of performing multiphase 3D CFD simulations has proven too large of a task for our current capabilities. Instead, this report will focus on the 3D model development for the 2.39 kW/L (31.3 kW) case and limitations and challenges associated with the 3D simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗

Development of a Fast Spectrum Test Loop at the Advanced Test Reactor

The United States currently has no infrastructure to provide a fast neutron flux for research purposes. With the increased interest in Sodium Fast Reactor (SFR) technology, this type of infrastructure is necessary for further development of key components such as Versatile Test Reactor. The purpose of this research is to characterize the thermal hydraulics of a test loop capable of providing SFR like operating conditions. This includes showcasing feasibility in mechanical design, reaching a high fast to thermal neutron flux ratio ideally near 15, an inlet temperature of 600 K, a temperature increase of 50 K from inlet to outlet, and test chamber radial temperature profiles similar to SFR fuel rods. Existing infrastructure must be used to provide the neutron flux for the test loop, and it was decided that the Advanced Test Reactor (ATR) was the optimal choice due to its availability and previous work conducted similar to this project. To achieve the goals of this project, several mechanical designs were put forward with the most feasible design being chosen. The neutronics and thermal hydraulics analysis were conducted in tandem, with an iterative process used to maximize efficiency of both systems. Analysis of the neutronics system was done using Attila and MCNP while the thermal hydraulic analysis was done using RELAP5. This document outlines the thermal hydraulic analysis, though it includes the results from the mechanical design and neutronics analysis. The loop design includes a primary system with molten sodium, a set of “booster fuel plates” to increase the fast flux within the test section, a water system to cool the booster fuel plates, a helium system to separate the sodium and water systems, and a thermal neutron shield to decrease the thermal flux. Each of these systems is modeled in the neutronics and thermal hydraulic analyses. The RELAP5 model showcased the feasibility of the design proposed by this report. For all designs tested, the desired 50 K temperature increase from inlet to outlet was reached. The fuel rod radial temperature profile was lower than was expected, likely due to a lower rod thermal power than expected. A parametric study was also conducted to characterize the system with three fuel pins. This study showed the effect of thermal power, fluid pressure, fluid velocity, and inlet temperature on the system. ©Copyright by Anton Higgins December 11, 2019 All Rights Reserved

99 GENERAL AND MISCELLANEOUS↗

Thermal Stress Modeling and Analysis of Packed-bed Thermocline Energy Storage Tank for INL Thermal Energy Distribution System (TEDS)

The Thermal Energy Distribution System (TEDS) at Idaho National Laboratory (INL) is a thermal-hydraulic flow loop to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit for TEDS. The packed-bed TES is adopted in TEDS because of its benefit as a low-cost single-tank storage option compared to the traditional two-tank storage. However, thermal ratcheting is one potential design concern which is caused by the rearrangement of granular filler inside a packed-bed tank during continuous thermal cycling operation of the packed-bed TES tank. If the thermally induced stress exceeds yield strength of the tank wall, it may cause catastrophic consequences like rupture of the thermal storage tank. Thus, it is crucial to understand the phenomenon to ensure the robust operation. Based on the temperature boundary conditions given by transient thermal analyses with computational fluid dynamics (CFD) simulations, the thermal ratcheting analysis is then conducted to evaluate the hoop stress and resultant thermal ratcheting potential of the TES tanks with two different modeling approaches: (1) infinite rigidity model and (2) Drucker-Prager (DP) model. The validity of each modeling method was examined by comparing the numerical simulation with the experimental data obtained from the packed-bed TES tank for Solar One Plant and evaluate the thermal ratcheting potential of the TEDS TES tank.

25 ENERGY STORAGE↗

A new dynamic zOnal model with air-diffuser (DOMA) - Application to thermal comfort prediction

A new Dynamic zOnal Model with Air-diffuser (DOMA) was developed. Several case studies were investigated and tested to evaluate and validate this program using measurement data. This new model was integrated into a TRaNsient SYstems Simulation program library and coupled with the multi-zone thermal model. The DOMA/TRNSYS coupled model was then used to predict room temperature distribution over an entire day of a single-zone building. The results show that increasing the heating outputs of the electric floor system, for example, from 75 to 200 W/m 2 , would not effectively improve the indoor thermal comfort, since the thermostat will reach the set point first and then turn off the system before the room gets enough heat and reach a comfortable level. This indicates the importance of selecting an appropriate location and set point for the thermostat when using a floor heating system. This potential thermal comfort issue can only be identified through the two-node model with a dynamic zonal model rather than the conventional PMV model, which thus suggests that for optimizing indoor thermal comfort of a building equipped with a time-sensitive control strategy and/or HVAC system, the TSENS results obtained from the two-node model integrated with DOMA are more appropriate than PMVs.

Construction & Building Technology↗

Thermal Analysis of a Solid Particle Light-Trapping Planar Cavity Receiver Using Computational Fluid Dynamics

Concentrated solar power (CSP) is one of the most effective ways of harnessing solar power to create efficient, durable, and resilient energy systems. This study entails thermal modeling and analysis of a novel central tower receiver configuration. This receiver uses solid particles as the heat transfer fluid (HTF), a promising option for third-generation CSP systems. The configuration considered here is the light-trapping planar cavity receiver (LTPCR) introduced by the National Renewable Energy Laboratory. While heat transfer studies of various LTPCR subsystems have been done, system-level thermal analysis of the LTPCR receiver has not been attempted. This study also presents important sensitivity analyses of the operating parameters of the CSP system, which can help guide the design of future central tower receivers. This study employs Ansys Fluent as a computational fluid dynamics (CFD) tool to model fluid dynamics and heat transfer in the receiver, intending to quantify its thermal performance. The model seamlessly integrates Monte Carlo ray tracing data, which generates absorbed solar flux profiles from the heliostat field design, with the heat transfer characteristics of the fluidized particle bed. This unified model is designed to accurately predict the thermal behavior of the LTPCR. Analysis of preliminary results reveals that the primary loss mechanisms are radiative and natural convective losses, in that order. Based on observations from a baseline case, several strategies are suggested and numerically tested. These solutions include selective cooling of high-temperature regions and manipulation of particle bed parameters. Selective cooling of high-temperature regions reduced the peak temperature by 151 degrees C and decreased thermal losses by 0.9%. Improving the particle-wall heat transfer coefficient (P-W HTC) of the particle bed decreased the thermal losses by 1.7% and decreased the peak temperatures by 57 degrees C. Decreasing the particle inlet temperature (PIT) also reduced thermal losses by 3.5% and decreased peak temperatures by 29 degrees C. Compounding these strategies improved the thermal losses of the receiver from 13.5% in the baseline case to 7.5%. Additionally, the study explores the variation in thermal performance across different locations of the receiver, where a variation of thermal losses from 12.9% to 17.3% is found. This allows a comprehensive evaluation of potential improvements in efficiency and temperature management.

computational fluid dynamics↗

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Improving AHU Performance by Minimizing Approach Temperature, Reducing Air Maldistribution, and Efficiently Handling Sensible and Latent Loads (Phase 1 Interim Final Technical Report)

Residential air handling units (AHUs) have stayed the same in form and efficiency for the past 30+ years, with incremental improvements made to address safety, functionality, and energy-efficiency. The purpose of this research in Phase I, Topic 9a: Next Generation Residential Air Handlers, was to improve AHU performance by minimizing heat exchanger (HX) approach temperature, reducing air maldistribution, and developing alternative system configurations which more efficiently handle sensible and latent loads. In this research Optimized Thermal Systems (OTS) developed, modeled, and evaluated multiple alternative system concepts. A dual vapor compression system separate sensible and latent cooling (SSLC) concept was studied to inform work on alternative concepts and to show best-case performance benefit. System concepts included ejector enhanced vapor compression cycles, desiccant assisted dehumidification, dual evaporator SSLC, and alternative AHU HX configurations. A dual vapor compression system showed COP improvement of 20%, however, required additional components, increased unit size, and increased cost. Two types of ejector enhanced vapor compression cycles with dual evaporators improved system COP by 9 to 11%, and reduced AHU losses by as much as 18%, with design changes limited to the AHU, no unit physical size increase, and a moderate increase to system first cost. Desiccant assisted air-conditioning required increased air flow rate resulting in higher fan power and the desiccant wheel increased sensible heat load leading to increased compressor power and reduced system COP. Dual evaporator cycles were found to degrade performance due to increased expansion losses. Optimized single slab HX designs used in place of the traditional A-coil HX led to 44–49% reduction in aluminum, 47–60% less refrigerant charge, and improved HX velocity distribution.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Machine Learning in Nuclear Thermal Hydraulics, Part 2: Opportunities and Perspectives

This paper summarizes the issues in nuclear thermal hydraulics (NTH), mainly from the three named aspects, and illustrates potential opportunities of machine learning (ML) applications in dealing with these issues. These issues are problems in simultaneously realizing computational efficiency and accuracy of system thermal hydraulic (STH) modeling and simulation, uncertainty in the scaling analysis for the best-estimate calculations of transient and accident analysis, and difficulties in optimal data production and usage for the development and assessment of evaluation models. Due to ML uncertainty, data insufficiency, and lack of cognition about how to efficiently incorporate knowledge and data, challenges of adapting ML techniques still exist. Accordingly, several perspectives are proposed to provide insights for ML applications in solving these issues.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal Modeling of Advanced Test Reactor Fuel in a Generalized Dry Storage System Under Hypothetical Accident Conditions

Work was performed to investigate the thermal behavior of a sealed aluminum-clad spent nuclear fuel (ASNF) dry storage configuration under hypothetical accident conditions during long-term storage. The considered system consists of a concrete dry storage overpack; a welded over canister, backfilled with argon; and nine Department of Energy (DOE) standardized canisters (DOESCs) loaded with ASNF, backfilled with helium gas. One technical concern associated with the long-term storage of ASNF includes radiolytic hydrogen generation. The yield of large quantities of hydrogen could lead to DOESC over pressurization or a flammable internal atmosphere. While flammability concerns can be resolved by preventing the ingress of oxygen, the canister pressure is partially controlled by the atmosphere temperature. This memorandum summarizes the results of modeling and simulation (M&S) work completed with Star CCM+ (a computational fluid dynamics [CFD] software) considering the thermal effects of a fully engulfing flame and another scenario with the vent ports of the concrete dry storage overpack completely blocked. The goal was to determine the bounding maximum temperature of the DOESC internal atmosphere to evaluate a worst-case pressure scenario, as well as the critical vent port blockage durations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Research Highlight - Dispatch Optimization, System Design and Cost Benefit Analysis of a Nuclear Reactor with Molten Salt Thermal Storage

We highlight our work from our previous publication titled "Dispatch Optimization, System Design, and Cost Benefit Analysis of a Nuclear Reactor with Molten Salt Thermal Storage" published in MDPI Energies 2022 (doi.org/10.3390/en15103599). Variable renewable energy availability has increased the volatility in energy prices in most markets. Nuclear power plants, with a large ratio of capital to variable costs, have historically operated as base load energy suppliers but the need for more flexible operation is increasing. We simulate the techno-economic performance of a 950 MWt nuclear power plant, based on the Westinghouse lead-cooled fast reactor, coupled with molten salt thermal storage as a method for flexible energy dispatch. We use the System Advisor Model to model the nuclear reactor thermal power input and power cycle operating modes. We combine this robust engineering model with a mixed-integer linear program model for optimized dispatch scheduling. We then simulate the coupled nuclear and thermal storage system under different market scenarios with varying price volatility. We find that the coupled plant outperforms the base plant under markets where energy price peaks fluctuate by a factor of two or more about the mean price. We show that a calculated power purchase agreement price for the plant improves by up to 10% when operating under California energy market conditions. Sensitivity analysis on the thermal storage cost shows that the optimal design remains unchanged even when doubling costs.

97 MATHEMATICS AND COMPUTING↗

Dynamic modeling of heat pipe integrated thermal battery latent heat storage system experiment validation

A heat pipe integrated thermal battery system has been constructed to investigate a high-temperature latent heat thermal energy storage technology that takes advantage of near isothermal operation of latent heat storage and heat pipes to potentially enable high-energy isothermal heat storage. A dynamic model constructed in Modelica has been validated, showing errors between 2.5 °C–39.7 °C across 10-h to 47-h simulations against experiment results, showing good prediction capability of experiment output, especially against phase change time. Model calibrations showing vessel heat-up capability of 3 kW and heat pipes combining to provide 600 W each during experiment operation validate experiment circumstances including reduced material loading and reduced power capability. The experiment configuration uses an Al-Mg-Zn eutectic metal as the storage material, heated via heat tape wrapped around the vessel and guide tubes to bring the system to operation range (>400 °C) and to simulate charging heat exchange, respectively, with heat rejection occurring through the surfaces of the material and facilitated via guide tubes with less insulation wrapping. The model is available in the open-source repository HYBRID on Github.

25 - ENERGY STORAGE↗

Forecasting Solar-Thermal Systems Performance under Transient Operation Using a Data-Driven Machine Learning Approach Based on the Deep Operator Network Architecture

Modeling and prediction of the dynamic behavior of thermal systems operating under intermittent energy input and variable load requirements represent one of the greatest challenges in the development of efficient and reliable renewable-based power generation technologies. In this work, a data-driven machine learning modeling framework was developed based on a modified version of the Deep Operator Network architecture where the time coordinate in the trunk net is replaced with historical data of the predicting quantity. The modeling framework can be used to accurately predict the performance of renewable-based energy conversion technologies including wind- and solar-based power plants. This novel framework was applied on a solar-thermal system that consists of a solar collection loop using a flat plate collector, a power generation loop comprising an Organic Rankine Cycle, and a thermal energy storage tank connecting both loops. Variable solar irradiance, air temperature, and power load profiles were used by the Deep Operator Network to predict the State-of-Charge and the efficiency of the thermal system for several days. The results were compared with the State-of-Charge and efficiency functions calculated using a physics-based model. For a simple operation scenario, characterized by a clear sky solar irradiance profile and constant load, the standard deviation in the State-of-Charge prediction by Deep Operator Network is below 0.9% during a seven-day prediction time horizon. For the most realistic operation scenario that considers real solar irradiance and a rough load profile, the maximum standard deviation in the predictions for the State-of-Charge and efficiency are below 6.8% and 2.5%, respectively. A comparison between Deep Operator Network and Long Short Term Memory network was also performed. In general, both networks predict very well the State-of-Charge for different data density conditions; however, a higher accuracy, with a standard deviation below 2.0%, is obtained by the Deep Operator Network during three and half days using sparser training data of 20-minute points. The same accuracy for the State-of-Charge prediction with the Long Short Term Memory network is achieved only for 14 h. Average standard deviations for the State-of-Charge prediction of 1.1% with the Deep Operator Network and 1.5% with the Long Short Term Memory network are obtained for a four-day prediction time using a denser training data of 5-minute points.

DeepONet↗

Liouvillian dynamics of the open Schwinger model: String breaking and kinetic dissipation in a thermal medium

Understanding the dynamics of bound state formation is one of the fundamental questions in confining quantum field theories such as Quantum Chromodynamics (QCD). One hadronization mechanism that has garnered significant attention is the breaking of a string initially connecting a fermion and an antifermion. Deepening our understanding of real-time string-breaking dynamics with simpler, lower dimensional models like the Schwinger model can improve our understanding of the hadronization process in QCD and other confining systems found in condensed matter and statistical systems. Here in this paper, we consider the string-breaking dynamics within the Schwinger model and investigate its modification inside a thermal medium, treating the Schwinger model as an open quantum system coupled to a thermal environment. Within the regime of weak coupling between the system and environment, the real-time evolution of the system can be described by a Lindblad evolution equation. We analyze the Liouvillian gaps of this Lindblad equation and the time dependence of the system’s von Neumann entropy. We observe that the late-time relaxation rate decreases as the environment correlation length increases. Moreover, when the environment correlation length is infinite, the system exhibits two steady states, one in each of the sectors with definite charge-conjugation-parity (CP) quantum numbers. For parameter regimes where an initial string breaks in vacuum, we observe a delay of the string breaking in the medium, due to kinetic dissipation effects. Conversely, in regimes where an initial string remains intact in vacuum time evolution, we observe string breaking (melting) in the thermal medium. We further discuss how the Liouvillian dynamics of the open Schwinger model can be simulated on quantum computers and provide an estimate of the associated Trotter errors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Performance Evaluation and Optimal Design Analysis of Continuous-Operation Solar-Driven Cooling Absorption Systems With Thermal Energy Storage

A solar absorption cooling system consisting of a flat plate collector, thermal energy storage tank, and absorption chiller is analyzed in this work. A dimensionless model is developed from the energy balance on each component and the chiller's characteristic performance curves. The model is used to determine the interaction and influence of different parameters such as tank size, solar collector area, chiller size, cooling load, cooling temperature, heat loss, and mass flow rates on the performance. From the analysis, smaller solar collector areas are required for lower cooling loads and smaller tank volumes. A specific cooling load of 1 x 10^-5 will require a specific solar collector area between two and six times larger, depending on the initial tank temperature, than the area required for a baseline system that considers typical commercial design and operation parameters. A similar behavior was observed for the specific tank volume. For the baseline system, the minimum specific area of the collector of 9.57 is achieved for an initial tank temperature of 1.19. For a cooling load of 1 x 10^-5, the optimum initial tank temperature will be 1.11 that results in a minimum specific solar collector area of 25.26. A specific tank volume of 4 x 10^-4 will also have an optimum initial tank temperature of 1.11 that minimizes the specific solar collector area to a value of 28.18. The approach and analysis in this work can be used to determine design parameters for solar absorption cooling systems based on a proper relation among system's dimensions to achieve optimum operation.

cooling load↗

HYBRID Modeling of DETAIL Experimental Facility

Idaho National Laboratory continues to be at the forefront of advanced reactor systems and integrated energy systems research. Through the Thermal Energy Distribution System (TEDS) and Microreactor Agile Non-nuclear Test Bed (MAGNET) experiments in the Dynamic Energy Transport and Innovation Laboratory (DETAIL) within the Engineering Systems Laboratory, a mix of digital and physical testing systems can explore the operation of thermally and electrically integrated systems. Modeling of these systems has been done using the HYBRID modeling repository of the Integrated Energy Systems program. DETAIL also contains high temperature steam electrolysis (HTSE) units to produce hydrogen and generate the load imposed on MAGNET and TEDS. The dynamic modeling capabilities within the IES program describe the DETAIL components as-built and can be used to evaluate controls, process flows, and overall system conditions.

08 HYDROGEN↗

A hierarchical framework for aggregating grid-interactive buildings with thermal and battery energy storage

The behind-the-meter (BTM) thermal and battery energy storage can help improve energy efficiency, reduce energy costs, and enhance energy resilience, particularly in rural areas and for disadvantaged communities. Aggregating numerous BTM energy storage systems can act as a price influencer with a significant source of load shifting and peak demand reduction. An integrated and scalable control mechanism is required to effectively utilize energy storage systems and flexible building loads to maximize the economic benefits, considering various distribution system constraints. Here, this paper presents an innovative hierarchical coordination framework for energy storage and flexible load in buildings, considering various factors such as electricity prices, thermal comfort, and distribution system modeling and constraints. At the upper level, a distribution system operator optimizes the power flow to minimize its power procurement costs from the electricity wholesale market, while at the lower level, aggregators determine the optimal dispatch of battery and thermal energy storage systems in multiple buildings on behalf of end-users to minimize operating costs according to the power prices. These problems are solved using a game-theoretic approach through negotiations between the distribution system operator and aggregators as a bi-level decision model. Simulation case studies have been performed for a test distribution network with a number of building end-users using energy storage systems to quantify the performance of aggregators. The results demonstrate that the proposed strategy can reduce peak load for a reliable electricity distribution network while saving electricity bills for customers.

25 ENERGY STORAGE↗