Search NASA⌕ Search

SEARCH · Search NASA

Results for “Landsat 10”

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,759 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Advanced Model Development for Large Eddy Simulation of Oxy-Combustion and Supercritical Carbon Dioxide Power Cycles

A joint experimental and numerical study is performed to observe the characteristics of a supercritical carbon dioxide turbulent mixing layer in the presence of strong nonlinearities in the thermodynamic and transport properties. A bespoke experimental setup is designed and employed for this purpose and provides insight into macroscopic mixing behavior. The mixing is experimentally observed using two techniques: shadowgraphy and spontaneous Raman scattering. Qualitative and quantitative intensity fields obtained via these techniques yield instantaneous and mean density data. Spanwise temperature data is also collected using analogue resistance temperature detectors. These measurements are used to quantify the level of mixed material within the field. The experimental data are supplemented by a companion high-fidelity numerical study. The numerical results are obtained through fully resolved, three-dimensional direct numerical simulation. The numerical dataset permits observation of the near-field mixing characteristics, which are difficult to measure experimentally due to the rapid dynamics and sharp thermophysical gradients in this area. Qualitative field visualizations are presented, followed by quantitative mixed material results and observations regarding thermodynamic property trends at select locations within the field. One-dimensional spectra of the turbulent kinetic energy and solenoidal dissipation are provided to observe the spectral characteristics of the flow. Reynolds stress anisotropy is analyzed graphically through anisotropy invariance maps (Lumley triangles). The mixing quantification, spectral data and anisotropy analysis of a flow at these thermodynamic conditions represent the main outcomes of the work.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Efficient Leaching of Metal Ions from Spent Li-Ion Battery Combined Electrode Coatings Using Hydroxy Acid Mixtures and Regeneration of Lithium Nickel Manganese Cobalt Oxide

Extensive use of Li-ion batteries in electric vehicles, electronics, and other energy storage applications has resulted in a need to recycle valuable metals Li, Mn, Ni, and Co in these devices. In this work, an aqueous mixture of glycolic and lactic acid is shown as an excellent leaching agent to recover these critical metals from spent Li-ion laptop batteries combined with cathode and anode coatings without adding hydrogen peroxide or other reducing agents. An aqueous acid mixture of 0.15 M in glycolic and 0.35 M in lactic acid showed the highest leaching efficiencies of 100, 100, 100, and 89% for Li, Ni, Mn, and Co, respectively, in an experiment at 120 °C for 6 h. Subsequently, the chelate solution was evaporated to give a mixed metal-hydroxy acid chelate gel. Pyrolysis of the dried chelate gel at 800 °C for 15 h could be used to burn off hydroxy acids, regenerating lithium nickel manganese cobalt oxide, and the novel method presented to avoid the precipitation of metals as hydroxide or carbonates. The Li, Ni, Mn, and Co ratio of regenerated lithium nickel manganese cobalt oxide is comparable to this metal ratio in pyrolyzed electrode coating and showed similar powder X-ray diffractograms, suggesting the suitability of α-hydroxy carboxylic acid mixtures as leaching agents and ligands in regeneration of mixed metal oxide via pyrolysis of the dried chelate gel.

Electrochemistry↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
DOE OSTI2026

Are Deep Energy Retrofits in Commercial Buildings Including Window Upgrades?

U.S. Commercial buildings account for about 20% of total U.S. energy consumption. Because the thermal performance of windows significantly affects building energy efficiency and HVAC system performance, best practice guidance often includes window and envelope improvements in conjunction with HVAC upgrades to optimize energy use and improve occupant comfort. It is an open question, however, regarding how often these best practices are implemented in the field. This paper aims to address that gap by conducting a literature review and a series of interviews with commercial building auditing and management professionals to explore the factors that drive window retrofits in commercial buildings. The paper explores a range of case studies from deep energy retrofits across the globe, comparing projects with and without window retrofits. The primary goals of this review are to: (1) provide data from real-world case studies illustrating the role of windows in deep energy renovations and HVAC upgrades, (2) conduct retrofit cost analyses for windows and high-performance HVAC systems and (3) offer insights into how window upgrade decisions are made and when they are implemented as part of deep energy retrofits. Most of the retrofit studies focused exclusively on high performance HVAC upgrades without considering how window upgrades might further enhance the overall energy efficiency of commercial buildings. Interviews with building industry experts shed light on the key factors influencing deep energy retrofit decisions and what factors tip the scales in favor of including window measures with more comprehensive retrofit projects.

Cort, Katherine↗
DOE OSTI2026

Large Language Model for Validation, Optical Calibration, and Learning (VOCAL) Distributed Temperature Sensing Interface

Distributed temperature sensing (DTS) using fiber optic sensors (FOS) offers a promising method for temperature measurements in advanced reactors, such as sodium fast reactors and molten salt cooled reactors. To support the calibration and validation of DTS measurements, Argonne National Laboratory developed the Validation, Optical Calibration, and Learning (VOCAL) software package. This report describes the integration of a local large language model (LLM) with a retrieval-augmented generation (RAG) system into the VOCAL interface to serve as an interactive user assistant. The LLM framework enhances the VOCAL platform’s accessibility to users by explaining interface components, clarifying inputs and outputs, and answering user queries dynamically in real-time. The accuracy of the LLM assistant performance was evaluated with 20 queries regarding the interface and its parameters using experimental data from the Thermal Hydraulic Experimental Test Article (THETA) facility. Results demonstrate that the LLM achieved a 95% accuracy rate, with a BERTScore of 0.8816 and SBERT value of 0.7417. Furthermore, validation of the RAG system within the LLM framework showed optimal accuracy with k-values between 1 and 2 using the k-refinement convergence test. The prompt perturbation analysis demonstrated good initial consistency for the RAG system, exhibiting the highest accuracy under punctuation variations and the greatest sensitivity under query reordering. Notably, the model’s errors were limited to data retrieval failures rather than factual hallucinations, reinforcing its baseline reliability. The integration of LLM provides a highly accurate, userfriendly enhancement to the VOCAL platform without disrupting its core computational capabilities for FOS calibration and validation.

Hong, Evan↗
DOE OSTI2026

Optimizing Porous Transport Layer Porosity for Proton Exchange Membrane Water Electrolysis

An empirical model is presented that describes anode-side losses related to porous transport layer (PTL) morphology in proton exchange membrane water electrolysis (PEMWE). The model is based on an advanced voltage breakdown analysis that links various overpotentials to PTL morphology. Custom Ti PTLs, spanning uncommonly low porosities (22 - 31%), were fabricated and analyzed with X-ray CT to obtain pore and particle size distributions. Particle size distributions were consistent across samples with an average particle diameter of 12.0?..mu..m, whereas average pore diameters ranged from 6.0 to 7.0?..mu..m. The PTLs were tested in standard PEMWE cell assemblies with anode catalyst loadings of 0.1 mgIr cm-2 to obtain polarization curves, electrochemical impedance spectra, and augmented Tafel analysis. The PTL-dependent anode side losses were deconvoluted and assigned to excess utilization, concentration, ion transport resistance, and electrical contact resistance overpotentials. The data and model reveal an optimal 20 - 28% PTL porosity region where utilization and contact resistance overpotentials are minimized without triggering concentration and ion transport losses related to water deprivation. The optimal PTL porosity depends on the operating current density and is demonstrated at realistic PEMWE water flow rates to establish PTL design guidance for operation at scale.

08 HYDROGEN↗
DOE OSTI2025

Interlayer Design for Halide Electrolytes in All‐Solid‐State Lithium Metal Batteries

Abstract All‐solid‐state lithium‐metal batteries (ASSLMBs) are promising for transportation electrification due to their superior safety and high energy density. Lithium halide electrolytes provide excellent processing flexibility, high ionic conductivity, and anodic stability (>4.1 V), making them highly compatible with high‐voltage cathodes, surpassing sulfide electrolytes (<2.1 V). Nevertheless, halide electrolytes suffer from low cathodic stability and form an electronically conductive interphase with lithium, resulting in a critical current density (CCD) of nearly zero. Herein, Li 3 YbCl 6 electrolytes are synthesized that are kinetically stable with lithium by forming an electronic insulating solid electrolyte interphase. Guided by critical overpotential criteria, a PI 3 interlayer is designed that transforms into Li 6 PI 3 upon contact with lithium, substantially reducing the interfacial resistance of Li 3 YbCl 6 against lithium to 34 Ω and achieving a high critical overpotential of 114 mV. By substituting Yb with Lu, Li 3 LuCl 6 electrolytes with Li 6 PI 3 interlayers reach a CCD of 1.0 mA cm −2 at a capacity of 1.0 mAh cm −2 , comparable to sulfide electrolytes but with higher oxidation stability. Additionally, Li 6 PI 3 enables stable cycling of Li//Li cells with Li 3 LuCl 6 electrolytes at 0.5 mA cm −2 for 400 cycles and maintains 86.5% capacity in Li//LiCoO 2 cells after 220 cycles at 30 °C, paving the way for high‐performance ASSLMBs.

Chemistry↗
DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗
DOE OSTI2024

Establishing accretion flares from supermassive black holes as a source of high-energy neutrinos

ABSTRACT The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

Astronomy & Astrophysics↗
DOE OSTI2023

Self‐Propelling Macroscale Sheets Powered by Enzyme Pumps

Nanoscale enzymes anchored to surfaces act as chemical pumps by converting chemical energy released from enzymatic reactions into spontaneous fluid flow that propels entrained nano‐ and microparticles. Enzymatic pumps are biocompatible, highly selective, and display unique substrate specificity. Utilizing these pumps to trigger self‐propelled motion on the macroscale has, however, constituted a significant challenge and thus prevented their adaptation in macroscopic fluidic devices and soft robotics. Using experiments and simulations, we herein show that enzymatic pumps can drive centimeter‐scale polymer sheets along directed linear paths and rotational trajectories. In these studies, the sheets are confined to the air/water interface. With the addition of appropriate substrate, the asymmetric enzymatic coating on the sheets induces chemically driven, buoyancy flows that controllably propel the sheet's motion on the air/water interface. The directionality and speed of the motion can be tailored by changing the pattern of the enzymatic coating, type of enzyme, and nature and concentration of the substrate. This work highlights the utility of biocompatible enzymes for generating motion in macroscale fluidic devices and robotics and indicates their potential utility for in vivo applications.

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

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Smart building HVAC control challenge: experience and solutions from the ADRENALIN project

A smart building HVAC control competition crowdsourced and compared algorithms on fair and equal ground using the standardized BOPTEST framework. The competition attracted 138 participants, but only 9% submitted valid solutions for the final stage, highlighting the complexity of advanced HVAC control design. The winning solutions showed significant potential to reduce energy use and cost by shifting demand, without compromising occupant comfort. Across scenarios, thermal energy cost reductions of 36–76% relative to a baseline, were achieved. In peak heat periods, the cost reduction leveraged limited energy use reduction (0–15%), but more significant energy price reduction (34–62%). This shows smart controls' ability to avoid as much as possible consumption during the morning peak hours, when spot prices are tendentially the highest. Hosting the competition has highlighted challenges in creating competitions that both are fair and promotes solutions that are transferable to real life implementation.

BOPTEST↗
DOE OSTI2026

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗
DOE OSTI2026

A physically interpretable precursor framework for sub-seasonal prediction of Northern Hemisphere flash flourishing

Flash flourishing describes rapid vegetation increases that can quickly reshape land–atmosphere exchanges and impacts on ecosystem, yet its large-scale precursors, circulation context, and sub-seasonal predictability remain poorly understood. Here, we identified onset-stage circulation regimes across northern extratropical latitudes (NEL; >30°N) using 200 and 1000 hPa geopotential height, and examined their regional expressions over eastern Asia, western North America, and Europe. Flash flourishing onset in East Asian was associated with a baroclinic circulation regime and was preceded by a North Atlantic sea surface temperature (SST) precursor at a four-pentad lead. In contrast, onset in western North American and European preferentially occurred under barotropic regimes, preconditioned by Great Plains soil moisture at three-pentad lead and North Atlantic SST at a four-pentad lead, respectively. Ridge regression forecasts revealed regime-dependent sub-seasonal predictability, with mean out-of-sample R 2 exceeding 0.3 up to lead times of two pentads in East Asia, three pentads in western North America, and four pentads in Europe. Together, these findings established a mechanistic and regionally specific framework for anticipating rapid vegetation greening at sub-seasonal timescales.

Kong, Xiangxu [Nanjing Univ. of Information Scienc↗
DOE OSTI2025

Invisible water in subducted crust: Lawsonite velocity anomalies under mantle conditions

Subduction zones transport significant amounts of water into Earth’s mantle, primarily through hydrous minerals such as lawsonite. However, the seismic detectability of lawsonite-bearing oceanic crust at mantle depths remains uncertain. To address this issue, we measured sound velocities of lawsonite up to 7.4 GPa and 600 °C. Both P- and S-wave velocities exhibited unexpected increases with temperature under high-pressure conditions. Our result suggest that hydrous oceanic crust exhibits higher seismic velocities than the surrounding mantle at depths of 150 to 250 km, resulting in high-velocity anomalies rather than the previously assumed low-velocity anomalies. Furthermore, the seismic velocity difference between hydrous and dry oceanic crust is less than 2%, making it challenging to distinguish between them using seismic velocities. This limitation may hinder the detection of the hydration state in subducted crust. In addition, lawsonite remains stable in 90% of subduction zones, and thus, such “seismically invisible water” may exist in most subducted slabs around the world.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Unraveling interphase-driven failure pathways in LiMn0.6Fe0.4PO4/graphite pouch cells

LiMnxFe1−xPO4 (LMFP) is a promising high-voltage, thermally stable, and earth-abundant cathode material, yet its practical application is limited by interphase instability and Mn dissolution. In this work, we systematically evaluate LiMn0.6Fe0.4PO4/graphite pouch cells using three electrolyte formulations including control carbonate electrolyte, control + 2 wt% vinylene carbonate (VC), and control + 2 wt% VC + 1 wt% 1,3,2-dioxathiolane 2,2-dioxide (DTD), to establish how electrolyte composition governs interphase chemistry and long-term degradation. Electrochemical testing shows that both additives are preferentially reduced prior to ethylene carbonate (EC) during cell formation, generating robust cathode-electrolyte interphase (CEI) and solid-electrolyte interphase (SEI) layers that suppress gas evolution and raise the first-cycle coulombic efficiency to 89.3%. Additionally, the dual-additive electrolyte delivers the most stable performance, retaining over 85% capacity after 600 cycles while minimizing impedance growth under long-term cycling at C/3 and 40 °C. Soft X-ray absorption spectroscopy confirms that VC + DTD effectively suppresses electrolyte oxidation at the cathode surface, and micro-X-ray fluorescence shows substantially reduced Mn dissolution and deposition on the graphite anode. Density functional theory simulations further provided insights into the structural and energetic influences of alkoxide species on the cathode surface, proposing a Mn2+ extraction mechanism. The combined experimental and computational findings establish a mechanistic link between electrolyte composition and interphase evolution, highlighting the effectiveness of electrolyte engineering for extending the operational lifetime of LMFP-based lithium-ion batteries.

Chak, Chanmonirath Michael↗
DOE OSTI2025

Tunable shear thickening, aging, and rejuvenation in suspensions of shape-memory-endowed liquid crystalline particles

The morphological features of particles, notably shape anisotropy, critically influence the rheological properties of dense suspensions, spanning both natural and engineered systems. This work explores the potential of using shape memory particles to dynamically regulate suspension fluid flow through controllable shape transformations. First, we synthesize shape-memory particles with programmable anisotropy from liquid crystal elastomers, such that the stiffness and shapes of the particles can be tuned by manipulating temperature. Our findings reveal that suspensions from such particles exhibit significant tunability in shear thickening behavior, transitioning from discontinuous shear thickening to a Newtonian-like response within a narrow temperature range of 60 ° C. This capability to modulate rheological responses in situ presents an approach for addressing processing challenges in many applications where control over flow behavior is paramount. Furthermore, we also show that suspensions composed of these anisotropic particles can undergo physical aging, and evolve into a glassy state. This state can be escaped upon activation of the shape memory effect. This reversibility underscores the potential for using such materials to engineer systems that can enter or come out of kinetic arrest by leveraging internal mechanical responses to external stimuli. The insights gained here not only broaden our understanding of the interplay between particle geometry and suspension dynamics but also pave the way for leveraging ensembles of stimuli-responsive objects to precisely control collective behaviors in many-body systems.

Science & Technology - Other Topics↗
DOE OSTI2024

Mitigating Sodium Ordering for Enhanced Solid Solution Behavior in Layered NaNiO 2 Cathodes

The O-type layered nickel oxides suffer from undesired cooperative Jahn–Teller distortion stemming from Ni 3+ ions and undergo multiple biphasic structural transformations during the insertion/extraction of large Na + ions, posing a significant challenge to stabilize the structural integrity. We present here a systematic investigation of the impact of substituting 5 % divalent (Mg 2+ ) or trivalent (Al 3+ or Co 3+ ) ions for Ni 3+ to alleviate Na + ion ordering and perturb the Jahn–Teller effect to enhance structural stability. Here, we gauge a fundamental understanding of the Mg−O and Na−O or Mg−O−Na bonding interactions, noting that the ionicity of the Mg−O bond deshields the electronic cloud of oxygen from Na + ions. Furthermore, calculations of the Van Vleck distortion modes reveal a relaxation of NiO 6 octahedra from Jahn–Teller distortion and a reduced electron density at the interlayer with Mg 2+ substitution. Long-range (operando X-ray diffraction) and short-range (magic angle spinning nuclear magnetic resonance) structural analyses provide insights into reduced ordering, allowing a stable continuous solid solution. Overall, Mg-substitution results in a high-capacity retention of ~96 % even after 100 cycles, showcasing the potential of this strategy for overcoming the structural instabilities and enhancing the performance of sodium-ion batteries.

25 ENERGY STORAGE↗