Search NASA⌕ Search

SEARCH · Search NASA

Results for “high mass loading”

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,404 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1958

High Mach Number, Low-Cowl-Drag, External-Compression Inlet With Subsonic Dump Diffuser

A zero-cowl-drag, all-external-compression inlet with an abrupt area discontinuity or subsonic dump diffuser has been proposed and demonstrated for high Mach number application. Isentropic compression was focused at the lip of a cylindrical cowl and limited by the requirements for shock attachment. Theoretical calculations of the performance of such an inlet showed that it may be suited for application at Mach numbers of 4 and above. In an experimental investigation at Mach 3.85, a 4.75-inch-diameter inlet attained a total-pressure recovery of 0.41 at a corresponding mass-flow ratio of 0.92. Subsonic diffusion over an area ratio of about 3.5 to 1.0 was accomplished in 1.25 inlet diameters. Essentially no difference in performance was obtained in a comparison of the dump and a conventional subsonic diffuser. On a range basis, the over-all performance of this inlet was at least comparable with that of the best on-design isentropic external-compression inlet previously reported. A variable Mach number configuration utilizing a telescoping spike with a subsonic dump diffuser is also included.

James F Connors↗
DOE OSTI2024

Massive red spiral galaxies in SDSS-IV MaNGA survey

ABSTRACT Massive red spiral galaxies (MRSGs) are supposed to be the possible progenitors of lenticular galaxies (S0s). We select a large sample of MRSGs ($M_*\gt 10^{10.5}\rm {\rm M}_{\odot }$) from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) DR17 using the g − r colour versus stellar mass diagram, along with control samples of blue spirals and S0s. Our main results are as follows: (1) After comparing the Sérsic index, concentration parameter, asymmetry parameter distribution, size–mass relation, and Σ1 (stellar mass surface density within the central 1 kpc)−mass relation, we find MRSGs are similar to S0s and have more compact and symmetric structures than blue spirals. MRSGs also resemble S0s in Dn4000, metallicity, Mgb/$\rm \left\langle Fe \right\rangle$, and V/σ radial profile. (2) By using MaNGA 2D spectra data, we separate the spatial regions into inner (R < 0.8Re) and outer (0.8 < R < 1.5Re) regions, and detect residual star formation in the outer regions of MRSGs. (3) When we select a sub-sample of MRSGs with NUV − r > 5, we find that they are completely star formation quenched in both inner and outer regions. Compared to optically selected MRSGs, NUV − r selected MRSGs appear to be more concentrated and have more massive dark matter haloes. The similarities between S0s and MRSGs suggest the possible evolutionary trend between MRSGs and S0s.

Astronomy & Astrophysics↗
DOE OSTI2026

Hazard and risk analysis framework for nuclear power plant–based integrated energy systems

Employing integrated energy systems (IESs) with nuclear power plants (NPPs) can improve NPP utilization by leveraging dedicated thermal and electric power delivery, but it may also increase operational safety risks. This paper presents a framework to identify and quantify hazards and risks for such IESs. The framework combines accidentology to review past industrial accidents with failure modes and effects analysis (FMEA) to identify potential future incidents. Hydrogen explosion and toxic chemical release hazards are of particular concern. Explosion consequences are quantified using the Bauwens-Dorofeev (Bauwens) and trinitrotoluene equivalent mass (TNT-EM) methods, while chemical release consequences are computed using the Gaussian atmospheric dispersion method. Operational disturbances from direct electrical and thermal integration that may affect NPP safety are modeled using probabilistic risk analysis (PRA). Hazards and risks are then evaluated for regulatory compliance. The framework is applied to IESs comprising pressurized or boiling water reactors supplying three levels of thermal and electrical power to industrial customers. Case studies include high-temperature steam electrolysis hydrogen plants of varying capacities and a synthetic fuel production plant. Sensitivity analysis examines piping component failures in the PRA model as a precursor to cost estimation for thermal extraction line design. Additionally, Fussel-Vessely (FV) and risk increase importance (RII) measures identify risk-informed design improvements for the thermal extraction system. FMEA highlights hazards such as loss of offsite power, prompt loss of electrical load, loss of thermal output, and immediate steam diversion, in addition to hydrogen explosions and toxic chemical releases. Both Bauwens and TNT-EM methods suggest maintaining several hundred meters of separation between the NPP and hydrogen facility to mitigate explosion risks. PRA results show a maximum initiating event frequency increase of 1.15% and an overall risk increase of 0.28%. Importance measure analysis identifies upstream pipe leak isolation components as critical. Evaluating the results against safety regulations, it is concluded that hazards and risks can be managed to comply with regulations through risk-informed thermal and electrical connection designs, component selection, maintenance programs, and safe separation distances between NPPs and integrated industrial facilities.

08 - HYDROGEN↗
NASA NTRS2026

Petermann Glacier on the Brink: Progress, Challenges and Insights

Petermann Glacier, the largest marine-terminating glacier in northern Greenland based on catchment area and ice discharge, plays a key role in regulating ice discharge from the Greenland Ice Sheet into the Arctic Ocean. With an upstream catchment connected to the ice sheet interior via a deep subglacial canyon, its future stability has major implications for sea level rise. In this review, we synthesize recent advances in understanding Petermann’s dynamics across three critical interfaces: the ice–ocean, ice–atmosphere, and ice–bed boundaries. At the surface, observations show that reanalysis products underestimate air temperatures and melt, while regional climate models diverge significantly in their estimates of surface mass balance, underscoring the need for improved in situ data and models. At the ocean boundary, enhanced basal melting driven by both subglacial runoff and Atlantic water intrusions is identified as the dominant driver of recent mass loss of Petermann Glacier, with continued warming posing a serious threat to the stability of the floating tongue. At the bed, new geophysical synthesis reveals complex geology, likely spatial variability in geothermal heat flux, and the influence of the megacanyon on seasonal hydrology and velocity fluctuations. Petermann’s mass balance has been negative in recent decades without corresponding flow acceleration. However, the glacier has undergone significant calving events, and the current rifting that began in September 2025 highlights its vulnerability. This upcoming calving event underscores the timeliness of this review, as it will put Petermann Glacier’s terminus at its most retreated position since records began in 1923. The anticipated retreat of the ice tongue also reduces buttressing and brings the terminus closer to the grounding zone, and modeling studies suggest that calving within 12 km of the grounding zone could potentially trigger dynamic retreat, accelerating ice discharge, and a doubling of flow speeds. We conclude that improved observations, sustained monitoring of oceanographic and atmospheric properties, and high-resolution modeling are critical to constraining projections of Petermann Glacier’s future and its role in the stability of the Greenland Ice Sheet.

Dominik Fahrner↗
DOE OSTI2026

Molecular origin of anisotropic shear elastoplasticity in chitin

Chitin nano- and mesoscale structures present in the exoskeleton of crustaceans exhibit exceptional longitudinal stiffness and toughness, rivaling or even exceeding that of many synthetic polymer architectures. Here, we reveal the origin of the asymmetric shear response in chitin multiscale architectures, marked by pronounced anisotropy in deformation. Under shear aligned with the molecular axis, chitin accommodates strain through coherent atomic rearrangements that enable elastic recovery. In contrast, shear applied perpendicular to this axis induces liquid-like plasticity via localized sliding. These results demonstrate the intrinsic mechanical anisotropy of chitin, underpinning a dual function in resisting repetitive loading while dissipating internal stress during high-strain events. Our findings establish the molecular basis of shear elastoplasticity in multiscale chitin structures, wherein axial elasticity supports energy storage in load-bearing regions, whereas transverse plasticity enables controlled energy dissipation. These atomic-scale insights lay a foundation for the predictive design of chitin-based materials with tunable strength-toughness profiles.

Wan, Zhangmin [University of British Columbia, Van↗
NASA NTRS2022

Optical properties of morphologically complex black carbon aerosols: Effects of coatings

Optical properties are computed for fractal-like aggregate black carbon (BC) aerosols coated with different substances. Two models are used for these aerosols: (i) the coated aggregate model (model I), where coating is added to the voids and surroundings of BC fractal-like aggregates; and (ii) the closed-cell aggregate model (model II), where coating is added concentrically to each monomer of BC fractal-like aggregates. Our results favor choosing the coated aggregate model (model I) to simulate scattering and absorption by coated BC aerosols because this model is morphologically more realistic, and because this model yields mass absorption cross section (MAC) and backscattering linear depolarization ratio (LDR) values that are consistent with field measurements. Moreover the corresponding computed degree of linear polarization (DoLP) and LDR values are very sensitive to changes in the coating volume fraction (fvol) and the coating refractive index (m). With the same absorbing BC core, the MAC value increases steadily with increasing fvol and m values. For example, using a density of 1.8 g/cm 3 and BC refractive index 1.95 + i0.79, the calculated MAC values for uncoated BC aggregates range from 6.2 to 6.8 m 2 /g at 0.55 µm. When coating material is applied to the BC aggregates, the calculated MAC values for model I particles increase to between 9.8 and 13.2 m 2 /g (depending on m) when fvol = 87.5%. The backscattering LDR values also tend to increase with the increasing m values for the shapes and sizes considered in our study. For model I particles, the backscattering LDR values span a wide range of 4.2–27.8% at a wavelength of 0.35 µm at fvol = 87.5% when m increases from 1.33 to 1.55. Our results are relevant to analyses of polarimetric and lidar observations of smoke particles, especially when these particles undergo hygroscopic growth.

Black carbon aerosols↗
NASA NTRSDate not supplied

Artemis II Mission Performance of the Orion Active Thermal Control System

The Orion spacecraft has recently completed Artemis II, the first crewed mission around the moon in 54 years. This test flight demonstrated vehicle capability as the foundation for all future missions under the Artemis program. As a major system of Orion, the Active Thermal Control System (ATCS) successfully managed crew and vehicle heat loads to provide adequate cooling and thermal comfort for the duration of the Artemis II mission. Thermal control on the Orion vehicle is managed through two redundant cooling loops in the Crew Module (CM) that absorb heat loads from the cabin and Air Revitalization System (ARS) heat exchangers, the Liquid Cooling Garment (LCG) heat exchanger, and the avionics cold plates, and then transfer those heat loads to the Service Module (SM) radiator loops and/or the Phase Change Material (PCM) and ammonia boilers when required through thermal configuration or mission phase. During Artemis II, Orion did not necessitate any supplemental cooling from ammonia boilers until service module separation, allowing the maximum availability of ammonia cooling prior to entry. This presentation will provide an overview of the Orion ATCS major subsystems, highlight the performance during key mission phases, and comment on the applicability of the observed thermal response to future Artemis missions.

Ellie Thurston↗
NASA NTRS1979

Dynamic analysis of a photovoltaic power system with battery storage capability

A photovolataic power system with a battery storage capability is analyzed. A dual battery current control concept is proposed, which enables the battery to either supply or accept power depending upon system environment and load conditions. A simulation of the power system, including the battery current control, is developed and evaluated. The evaulation demonstrate the visbility of the battery control concept of switch the battery from a charge to discharge mode and back as required by load and environmental conditions. An acceptable system operation is demonstrated over the entire insolation range. Additionally, system sensitivity, bandwidth, and damping characteristics of the battery control are shown to be acceptable for a projected hardware implementation.

Merrill, W. C.↗
NASA NTRS1982

Contaminant characterization of five satellite materials

An extensive laboratory test program was performed to characterize outgassing of five satellite materials. The materials were Chemglaze Z-306 over 9922 primer, M-773 adhesive, multilayer insulation, polyurethane foam, and silverized Teflon. The souces were prepared to the specifications of a typical satellite program. Dynamic thermogravimetric mass loss characteristics of these five materials were obtained in vacuum with a beam microbalance. The temperature of the material was linearly raised from 25 C to over 650 C while monitoring the mass loss, rate of mass loss, temperature, and the composition of the outgassed material by residual gas analysis. Isothermal source emission/capture coefficients/reemission parameters were obtained with an array of quartz crystal microbalances (QCM). A detailed test matrix in which the temperatures of the QCMs (-160 C, - 100 C, -40 C, and +10 C) and the source material (125 C, and 90 C, and 50 C) were varied was performed. The array of QCMs was also used to measure the spatial distribution of the source emission. The effects of vacuum ultraviolet radiation on the deposition and reemission parameters was determined.

J A Muscari↗
DOE OSTI2024

Multiple Peaks and a Long Precursor in the Type IIn Supernova 2021qqp: An Energetic Explosion in a Complex Circumstellar Environment

Abstract We present optical photometry and spectroscopy of the Type IIn supernova (SN) 2021qqp. Its unusual light curve is marked by a long precursor for ≈300 days, a rapid increase in brightness for ≈60 days, and then a sharp increase of ≈1.6 mag in only a few days to a first peak ofM r ≈ −19.5 mag. The light curve then declines rapidly until it rebrightens to a second distinct peak ofM r ≈ −17.3 mag centered at ≈335 days after the first peak. The spectra are dominated by Balmer lines with a complex morphology, including a narrow component with a width of ≈1300 km s −1 (first peak) and ≈2500 km s −1 (second peak) that we associate with the circumstellar medium (CSM) and a P Cygni component with an absorption velocity of ≈8500 km s −1 (first peak) and ≈5600 km s −1 (second peak) that we associate with the SN–CSM interaction shell. Using the luminosity and velocity evolution, we construct a flexible analytical model, finding two significant mass-loss episodes with peak mass loss rates of ≈10 and ≈5M ⊙ yr −1 about 0.8 and 2 yr before explosion, respectively, with a total CSM mass of ≈2–4M ⊙ . We show that the most recent mass-loss episode could explain the precursor for the year preceding the explosion. The SN ejecta mass is constrained to be ≈5–30M ⊙ for an explosion energy of ≈(3–10) × 10 51 erg. We discuss eruptive massive stars (luminous blue variable, pulsational pair instability) and an extreme stellar merger with a compact object as possible progenitor channels.

Astronomy & Astrophysics↗
NASA NTRS2006

Effect of Thermomechanical Processing on the Microstructure, Properties, and Work Behavior of a Ti50.5 Ni29.5 Pt20 High-Temperature Shape Memory Alloy

TiNiPt shape memory alloys are particularly promising for use as solid state actuators in environments up to 300 °C, due to a reasonable balance of proper ties, including acceptable work output. However, one of the challenges to commercializing a viable high-temperature shape memory alloy (HTSMA) is to establish the appropriate primary and secondary processing techniques for fabrication of the material in a required product form such as rod and wire. Consequently, a Ti 50.5 Ni 29.5 Pt 20 alloy was processed using several techniques including single- pass high-temperature extrusion, multiple-pass high-temperature extrusion, and cold drawing to produce bar stock, thin rod, and fine wire, respectively. The effects of heat treatment on the hardness, grain size, room temperature tensile properties, and transformation temperatures of hot- and cold-worked material were examined. Basic tensile properties as a function of temperature and the strain-temperature response of the alloy under constant load, for the determination of work output, were also investigated for various forms of the Ti 50.5 Ni 29.5 Pt 20 alloy, including fine wire.

Grain Size↗
NASA NTRSDate not supplied

Li-ion Battery Aging with Hybrid Physics-Informed Neural Networks and Fleet-wide Data

In this work, we propose a hybrid model for Li-ion battery discharge and aging prediction that leverages fleet-wide data to predict future capacity drops.The model is built upon an hybrid approach merging physics-based and empirical equations, as well as neural network models in a recurrent neural network cell. The hybrid physics-informed neural network can predict voltage discharge cycles given the loading profile, and estimate the used capacity of the battery under random-loading conditions by tracking aging parameters connected to the residual capacity of the battery. By merging information on the battery aging parameters with existing fleet-wide aging data, the model can predict the future residual capacity of the battery that is being monitored, and therefore enable predictions of voltage discharge curves far ahead in the battery life cycle. We validated the approach using the NASA Prognostics Data Repository Battery data-set, which contains experimental data on Li-ion batteries discharged at random loading conditions in a controlled environment. The approach also allows the identification of discrepancies between the battery aging trend and the trend observed at the fleet level, so that batteries behaving differently from the rest of the fleet can be subject to closer monitoring and further testing to refine predictions.

PINN↗
NASA NTRS2020

An Ultra-long Life, High-performance, Flexible Li-CO2Battery Based on Multifunctional Carbon Electrocatalysts

Integrating CO2 utilization and renewable energy delivery/storage, the rechargeable Li–CO2 battery has been considered as a promising candidate for next-generation secondary batteries. However, high-performance catalyst(s) for efficient formation and decomposition of the discharge product, Li2CO3, are an imperative part of a Li–CO2 battery. The development of flexible Li–CO2 batteries extends their applications into compliant and wearable devices/systems, but at the same time imposes a big challenge for battery fabrication and lifetime enhancement. In this study, a rechargeable quasi-solidus flexible Li–CO2 battery was designed and fabricated using highly active N,S-doped carbon nanotubes (N,S-doped CNTs) as the cathode catalyst, and a smart polymer gel as the flexible electrolyte. This newly-developed flexible Li–CO2 battery exhibited a capacity as high as 23560 mAh g−1 based on the catalyst mass and an ultra-long lifetime of up to 538 cycles with excellent mechanical flexibility. This work provides a platform for the design and development of high-performance flexible Li–CO2 batteries from low-cost, earth-abundant, carbon-based multifunctional cathode catalysts.

Superior stability↗
NASA NTRS2012

Returning an Entire Near-Earth Asteroid in Support of Human Exploration Beyond Low-Earth Orbit

This paper describes the results of a study into the feasibility of identifying, robotically capturing, and returning an entire Near-Earth Asteroid (NEA) to the vicinity of the Earth by the middle of the next decade. The feasibility of such an asteroid retrieval mission hinges on finding an overlap between the smallest NEAs that could be reasonably discovered and characterized and the largest NEAs that could be captured and transported in a reasonable flight time. This overlap appears to be centered on NEAs roughly 7 m in diameter corresponding to masses in the range of 250,000 kg to 1,000,000 kg. The study concluded that it would be possible to return a approx.500,000-kg NEA to high lunar orbit by around 2025. The feasibility is enabled by three key developments: the ability to discover and characterize an adequate number of sufficiently small near-Earth asteroids for capture and return; the ability to implement sufficiently powerful solar electric propulsion systems to enable transportation of the captured NEA; and the proposed human presence in cislunar space in the 2020s enabling exploration and exploitation of the returned NEA. Placing a 500-t asteroid in high lunar orbit would provide a unique, meaningful, and affordable destination for astronaut crews in the next decade. This disruptive capability would have a positive impact on a wide range of the nation's human space exploration interests. It would provide a high-value target in cislunar space that would require a human presence to take full advantage of this new resource. It would offer an affordable path to providing operational experience with astronauts working around and with a NEA that could feed forward to much longer duration human missions to larger NEAs in deep space. It represents a new synergy between robotic and human missions in which robotic spacecraft would retrieve significant quantities of valuable resources for exploitation by astronaut crews to enable human exploration farther out into the solar system. The capture, transportation, examination, and dissection of an entire NEA would provide valuable information for planetary defense activities that may someday have to deflect a much larger near-Earth object. Transportation of the NEA to lunar orbit with a total flight time of 6 to 10 years would be enabled by a ~40-kW solar electric propulsion system with a specific impulse of 3,000 s. The flight system could be launched to low-Earth orbit (LEO) on a single Atlas V-class launch vehicle, and return to lunar orbit a NEA with at least 28 times the mass launched to LEO. Longer flight times, higher power SEP systems, or a target asteroid in a particularly favorable orbit could increase the mass amplification factor from 28-to-1 to 70-to-1 or greater. The NASA GRC COMPASS team estimated the full life-cycle cost of an asteroid capture and return mission at approx.$2.6B.

asteroid mining↗
DOE OSTI2026

Development of a Framework and Methodology for an Advanced Reactor Materials Environmental Effects Design Guide

Advanced non-light-water reactor components may operate at elevated temperature while experiencing cyclic loading, significant neutron irradiation, and exposure to reactor coolant. ASME Boiler and Pressure Vessel Code, Section III, Division 5, provides design rules for elevated-temperature service but does not include specific procedures to account for environmental effects on material properties. This report develops an initial framework and methodology for an Environmental Effects Design Guide (EEDG) focused on neutron irradiation; coolant-environment effects are reserved for future work. The proposed approach treats irradiation as a property-based overlay on the existing Division 5 design process, with two routes: a sparse-data route applying two reduction factors — FCR on creep-rupture strength and FF on fatigue life — for the creep-fatigue evaluations that typically control the design of advanced high-temperature reactor components, and a fuller framework developing the property-to-rule chain across the four Division 5 checks (primary load, strain limits and ratcheting, creep-fatigue, and buckling), together with swelling and weldments as scope items. Both routes are scoped by an in-pile qualification that restricts the use of post-irradiation-examination-derived properties in regimes where an in-pile mechanism could control the design outcome. Illustrative outputs derived on a compiled annealed Type 316 database — FCR ≈ 0.78–0.86 and FF ≈ 0.4 — demonstrate the calculation method within that specific dataset. The framework is an initial, testable design-rule concept; it identifies a practical path for preliminary design evaluations under sparse data and the material data and testing needed to develop the framework further.

Barua, Bipul (ORCID:0000000247184113)↗
NASA NTRS1992

The Utilization of Advanced Composites in Military Aircraft

The purpose of this workshop was to identify the current state of the art in key issues related to compression loading and fluid effects in composite materials. In the area of compression loading, there was considerable concern over the different results obtained from various test methods. It was agreed that failure modes produced by the various test methods along with a better fundamental understanding of compression failure were key issues in the development of compression test methods. In the area of fluid effects, a lack of a comprehensive data base hampers identification of key mechanisms leading to fluid degradation. This is further complicated by the fact that interactions depend on the fluid and composite under consideration.

Advisory Group for Aerospace Research and Developm↗
NASA NTRS1994

Elastic-plastic models for multi-site damage

This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model.

Actis, Ricardo L.↗
DOE OSTI2026

Storing Affordability: Battery Storage as an Asset to Reduce Data Center Cost Shifts

This report examines how battery energy storage systems (BESS) can help utilities accommodate large load growth while protecting affordability for existing ratepayers. Rapid growth in electricity demand from artificial intelligence (AI) data centers is straining the U.S. grid. Furthermore, many new data centers are entering rural markets, which could offer economic benefits but may also pose implementation challenges for smaller utilities. At the same time, retail electricity prices are increasing faster than inflation, elevating customer affordability as a key challenge. While data centers have not been the primary driver of increases in residential prices to date, they have pushed wholesale energy and capacity prices higher in several markets. Fundamental utility cost-allocation principles show that data center growth can be rate-positive for existing customers only if new peak demand grows faster than the costs a utility must incur to serve it. Several factors, including a utility’s degree of wholesale market exposure, forecast uncertainty and stranded-asset risk, and tariff design can determine the outcome of load growth on retail rates. Energy storage can make several affordability contributions in the face of this landscape of uncertainty and market volatility, including deferral of higher-cost grid investments through improved utilization of existing assets and flexibility of new large loads, insulation from volatile wholesale prices through peak shaving, and reliability support to address grid risks stemming from the behavior of AI data center loads. Different potential BESS deployment pathways—utility-scale front-of-the-meter systems, aggregated small-scale storage installations, and data center-sited behind-the-meter storage—are compared against each other and against conventional capacity alternatives. This framework is intended as a conceptual resource to utilities, particularly smaller public utilities with rural service territories, who may be considering the role that energy storage can play in insulating existing ratepayers from data center cost shifts.

25 ENERGY STORAGE↗