Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape Decomposition”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mycorrhiza better predict soil fungal community composition and function than aboveground traits in temperate forest ecosystems

We sampled soil in the organic and mineral horizons beneath two AM-associated (Fraxinus americana, Thuja occidentalis) and two ECM-associated tree species (Betula alleghaniensis, and Tsuga canadensis), with an evergreen and deciduous species in each mycorrhizal group. To characterize fungal communities and organic matter decomposition beneath each tree species, we sequenced the ITS1 region of fungal DNA and measured the potential activity of carbon and nitrogen-targeting extracellular enzymes. Each tree species harbored distinct fungal communities, supporting the need to consider both mycorrhizal type and leaf habit. However, between tree characteristics, mycorrhizal type better predicted fungal communities. Across fungal guilds, saprotrophic fungi were the most important group in shaping fungal community differences in soils beneath all tree species. The effect of leaf habit on carbon and nitrogen-targeting hydrolytic enzymes depended on tree mycorrhizal association in the organic horizon, while oxidative enzyme activities were higher beneath EcM-associated trees across both soil horizons and leaf habits. These data include extracellular hydrolytic and oxidative enzyme activities, ITS sequencing fungal community data, soil carbon to nitrogen, and soil pH. Site level data include climate (mean annual temperature, precipitation), elevation, and soil series and order information.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effect of transport mechanism on the coarsening of bicontinuous structures: A comparison between bulk and surface diffusion

Here, coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to t 1/4 for surface diffusion and t 1/3 for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamic mode decomposition for gyrokinetic eigenmode analysis

Dynamic mode decomposition (DMD) is a post-processing approach to decompose a complex time series into a set of modes via spectral analysis. DMD provides a new and powerful method to recover gyrokinetic drift-wave eigenfrequencies and eigenfunctions based only on the solution of the gyrokinetic-Maxwell initial value problem with almost no added cost to the initial value solver. In the present paper, DMD is applied to the CGYRO gyrokinetic code using a newly-developed CGYRO-DMD post-processor. CGYRO-DMD is numerically efficient, even on a single CPU. It does not set any restrictions on the plasma shape, beta (ratio of the plasma pressure to the magnetic field pressure), collisionality or number of species, and allows one to resolve numerous eigenmodes, even of comparable growth rates. In addition, DMD is not limited to unstable modes, but rather can capture stable and unstable branches simultaneously. In this work, we illustrate the accuracy of DMD through gyrokinetic analysis of mode transition for electromagnetic drift wave instabilities.

drift-wave eigenmodes↗
DOE OSTI2023

Spatially quasi-periodic bifurcations from periodic traveling water waves and a method for detecting bifurcations using signed singular values

We present a method of detecting bifurcations by locating zeros of a signed version of the smallest singular value of the Jacobian. This enables the use of quadratically convergent root-bracketing techniques or Chebyshev interpolation to locate bifurcation points. Only positive singular values have to be computed, though the method relies on the existence of an analytic or smooth singular value decomposition (SVD). The sign of the determinant of the Jacobian, computed as part of the bidiagonal reduction in the SVD algorithm, eliminates slope discontinuities at the zeros of the smallest singular value. We use the method to search for spatially quasi-periodic traveling water waves that bifurcate from large-amplitude periodic waves. The water wave equations are formulated in a conformal mapping framework to facilitate the computation of the quasi-periodic Dirichlet-Neumann operator. We find examples of pure gravity waves with zero surface tension and overhanging gravity-capillary waves. In both cases, the waves have two spatial quasi-periods whose ratio is irrational. We follow the secondary branches via numerical continuation beyond the realm of linearization about solutions on the primary branch to obtain traveling water waves that extend over the real line with no two crests or troughs of exactly the same shape. The pure gravity wave problem is of relevance to ocean waves, where capillary effects can be neglected. Such waves can only exist through secondary bifurcation as they do not persist to zero amplitude. The gravity-capillary wave problem demonstrates the effectiveness of using the signed smallest singular value as a test function for multi-parameter bifurcation problems. This test function becomes mesh independent once the mesh is fine enough.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Global distribution, formation and fate of mineral‐associated soil organic matter under a changing climate: A trait‐based perspective

Abstract Soil organic matter (SOM) is the largest actively cycling reservoir of terrestrial carbon (C), and the majority of SOM in Earth's mineral soils (~65%) is mineral‐associated organic matter (MAOM). Thus, the formation and fate of MAOM can exert substantial influence on the global C cycle. To predict future changes to Earth's climate, it is critical to mechanistically understand the processes by which MAOM is formed and decomposed, and to accurately represent this process‐based understanding in biogeochemical and Earth system models. In this review, we use a trait‐based framework to synthesize the interacting roles of plants, soil micro‐organisms, and the mineral matrix in regulating MAOM formation and decomposition. Our proposed framework differentiates between plant and microbial traits that influence total OM inputs to the soil (‘feedstock traits’) versus traits that influence the proportion of OM inputs that are ultimately incorporated into MAOM (‘MAOM formation traits’). We discuss how these feedstock and MAOM formation traits may be altered by warming, altered precipitation and elevated carbon dioxide. At a planetary scale, these feedstock and MAOM formation traits help shape the distribution of MAOM across Earth's biomes, and modulate biome‐specific responses of MAOM to climate change. We leverage a global synthesis of MAOM measurements to provide estimates of the total amount of MAOM‐C globally (~840–1540 Pg C; 34%–51% of total terrestrial organic C), and its distribution across Earth's biomes. We show that MAOM‐C concentration is highest in temperate forests and grasslands, and lowest in shrublands and savannas. Grasslands and croplands have the highest proportion of soil organic carbon (SOC) in the MAOM fraction (i.e. the MAOM‐C:SOC ratio), while boreal forests and tundra have the lowest MAOM‐C:SOC ratio. Drawing on our trait framework, we then review experimental data and posit the effects of climate change on MAOM pools in different biomes. We conclude by discussing how MAOM is integrated into soil C models, and how feedstock and MAOM formation traits may be included in these models. We also summarize the projected fate of MAOM under climate change scenarios (Representative Concentration Pathways 4.5 and 8.5) and discuss key model uncertainties. Read the free Plain Language Summary for this article on the Journal blog.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessment of four strain energy decomposition methods for phase field fracture models using quasi-static and dynamic benchmark cases

Abstract Strain energy decomposition methods in phase field fracture models separate strain energy that contributes to fracture from that which does not. However, various decomposition methods have been proposed in the literature, and it can be difficult to determine an appropriate method for a given problem. The goal of this work is to facilitate the choice of strain decomposition method by assessing the performance of three existing methods (spectral decomposition of the stress or the strain and deviatoric decomposition of the strain) and one new method (deviatoric decomposition of the stress) with several benchmark problems. In each benchmark problem, we compare the performance of the four methods using both qualitative and quantitative metrics. In the first benchmark, we compare the predicted mechanical behavior of cracked material. We then use four quasi-static benchmark cases: a single edge notched tension test, a single edge notched shear test, a three-point bending test, and a L-shaped panel test. Finally, we use two dynamic benchmark cases: a dynamic tensile fracture test and a dynamic shear fracture test. All four methods perform well in tension, the two spectral methods perform better in compression and with mixed mode (though the stress spectral method performs the best), and all the methods show minor issues in at least one of the shear cases. In general, whether the strain or the stress is decomposed does not have a significant impact on the predicted behavior.

Zhang, Shuaifang↗
DOE OSTI2023

Manipulation of Scattering Spectra with Topology of Light and Matter

Structured lights, including beams carrying spin and orbital angular momenta, radially and azimuthally polarized vector beams, as well as spatiotemporal optical vortices, have attracted significant interest due to their unique amplitude, phase front, polarization, and temporal structures, enabling a variety of applications in optical and quantum communications, micromanipulation, and super-resolution imaging. In parallel, structured optical materials, metamaterials, and metasurfaces consisting of engineered unit cells—meta-atoms, opened new avenues for manipulating the flow of light and optical sensing. While several studies explored structured light effects on the individual meta-atoms, their shapes are largely limited to simple spherical geometries. However, the synergy of the structured light and complex-shaped meta-atoms has not been fully explored. Here, in this paper, the role of the helical wavefront of Laguerre–Gaussian beams in the excitation and suppression of higher-order resonant modes inside all-dielectric meta-atoms of various shapes, aspect ratios, and orientations, is demonstrated and the excitation of various multipolar moments that are not accessible via unstructured light illumination is predicted. The presented study elucidates the role of the complex phase distribution of the incident light in shape-dependent resonant scattering, which is of utmost importance in a wide spectrum of applications ranging from remote sensing to spectroscopy.

42 ENGINEERING↗
DOE OSTI2025

Modeling Plant Nutrient Acquisition Strategies Alters Projections of Carbon and Nitrogen Dynamics in Bioenergy Agroecosystems

Plant strategies to acquire nutrients from limited environments help shape ecosystem carbon (C) and nitrogen (N) cycling and response to environmental change. The effects of plant strategies on ecosystem dynamics are largely uncharacterized in bioenergy agroecosystems, where the impacts could determine bioenergy's ability to meet its sustainability goals of storing C and reducing N loss. We used FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms and Protection), a plant–microbe interaction model of coupled plant nutrient uptake and soil organic matter decomposition, to simulate the effects of nutrient acquisition strategies on soil microbial activity and ecosystem nutrient cycling in bioenergy feedstocks miscanthus (Miscanthus × giganteus) and sorghum (Sorghum bicolor (L.) Moench). We examined the model's ability to reproduce the relative effects of belowground nutrient uptake on microbial activity using a reanalysis of empirical data showing that miscanthus root exudation provoked a larger soil microbial response than sorghum. From baseline model simulations, we found that the ability of miscanthus to retranslocate N resulted in higher N uptake at a lower C cost than the sorghum/soybean rotation and that soil C and N pools increased under perennial (miscanthus) and decreased under annual (sorghum/soybean) cultivation. The model also predicted that greater root exudation increased soil C accumulation, highlighting the role of roots in forming stable soil C. Overall, the baseline model was unable to reproduce field observations of miscanthus root exudation stimulating microbial activity more than sorghum. To improve the model, we updated the soil microbial parameters in miscanthus to have faster decomposition, a higher C/N ratio, and greater carbon use efficiency. These changes improved the simulated soil microbial response to miscanthus root exudation, supporting the hypothesis that miscanthus soils foster a microbial community that is more responsive to root exudation than that of sorghum.

09 BIOMASS FUELS↗
DOE OSTI2025

Metasurfaces with Multipolar Resonances and Enhanced Light–Matter Interaction

Metasurfaces, composed of engineered nanoantennas, enable unprecedented control over electromagnetic waves by leveraging multipolar resonances to tailor light–matter interactions. This review explores key physical mechanisms that govern their optical properties, including the role of multipolar resonances in shaping metasurface responses, the emergence of bound states in the continuum (BICs) that support high-quality factor modes, and the Purcell effect, which enhances spontaneous emission rates at the nanoscale. These effects collectively underpin the design of advanced photonic devices with tailored spectral, angular, and polarization-dependent properties. This review discusses recent advances in metasurfaces and applications based on them, highlighting research that employs full-wave numerical simulations, analytical and semi-analytic techniques, multipolar decomposition, nanofabrication, and experimental characterization to explore the interplay of multipolar resonances, bound and quasi-bound states, and enhanced light–matter interactions. A particular focus is given to metasurface-enhanced photodetectors, where structured nanoantennas improve light absorption, spectral selectivity, and quantum efficiency. By integrating metasurfaces with conventional photodetector architectures, it is possible to enhance responsivity, engineer photocarrier generation rates, and even enable functionalities such as polarization-sensitive detection. The interplay between multipolar resonances, BICs, and emission control mechanisms provides a unified framework for designing next-generation optoelectronic devices. This review consolidates recent progress in these areas, emphasizing the potential of metasurface-based approaches for high-performance sensing, imaging, and energy-harvesting applications.

Kerker effect↗
DOE OSTI2023

Roots selectively decompose litter to mine nitrogen and build new soil carbon

Abstract Plant–microbe interactions in the rhizosphere shape carbon and nitrogen cycling in soil organic matter (SOM). However, there is conflicting evidence on whether these interactions lead to a net loss or increase of SOM. In part, this conflict is driven by uncertainty in how living roots and microbes alter SOM formation or loss in the field. To address these uncertainties, we traced the fate of isotopically labelled litter into SOM using root and fungal ingrowth cores incubated in a Miscanthus x giganteus field. Roots stimulated litter decomposition, but balanced this loss by transferring carbon into aggregate associated SOM. Further, roots selectively mobilized nitrogen from litter without additional carbon release. Overall, our findings suggest that roots mine litter nitrogen and protect soil carbon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sea ice and atmospheric circulation shape the high-latitude lapse rate feedback

Abstract Arctic amplification of anthropogenic climate change is widely attributed to the sea-ice albedo feedback, with its attendant increase in absorbed solar radiation, and to the effect of the vertical structure of atmospheric warming on Earth’s outgoing longwave radiation. The latter lapse rate feedback is subject, at high latitudes, to a myriad of local and remote influences whose relative contributions remain unquantified. The distinct controls on the high-latitude lapse rate feedback are here partitioned into “upper” and “lower” contributions originating above and below a characteristic climatological isentropic surface that separates the high-latitude lower troposphere from the rest of the atmosphere. This decomposition clarifies how the positive high-latitude lapse rate feedback over polar oceans arises primarily as an atmospheric response to local sea ice loss and is reduced in subpolar latitudes by an increase in poleward atmospheric energy transport. The separation of the locally driven component of the high-latitude lapse rate feedback further reveals how it and the sea-ice albedo feedback together dominate Arctic amplification as a coupled mechanism operating across the seasonal cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A fast multi-fidelity method with uncertainty quantification for complex data correlations: Application to vortex-induced vibrations of marine risers

Here we develop a fast multi-fidelity modeling method for very complex correlations between high- and low-fidelity data by working in modal space to extract the proper correlation function. We apply this method to infer the amplitude of motion of a flexible marine riser in cross-flow, subject to vortex-induced vibrations (VIV). VIV are driven by an absolute instability in the flow, which imposes a frequency (Strouhal) law that requires a matching with the impedance of the structure; this matching is easily achieved because of the rapid parametric variation of the added mass force. As a result, the wavenumber of the riser spatial response is within narrow bands of uncertainty. Hence, an error in wavenumber prediction can cause significant phase-related errors in the shape of the amplitude of response along the riser, rendering correlation between low- and high-fidelity data very complex. Working in modal space as outlined herein, dense data from low-fidelity data, provided by the semi-empirical computer code VIVA, can correlate in modal space with few high-fidelity data, obtained from experiments or fully-resolved CFD simulations, to correct both phase and amplitude and provide predictions that agree very well overall with the correct shape of the amplitude response. We also quantify the uncertainty in the prediction using Bayesian modeling and exploit this uncertainty to formulate an active learning strategy for the best possible location of the sensors providing the high fidelity measurements.

42 ENGINEERING↗
DOE OSTI2026

Hydrological Control on Soil Redox Condition and Carbon Loss of Coastal Wetland Under Sea-Level Rise

Coastal wetlands are critical carbon sinks with their biogeochemical and ecological functioning shaped by dynamic hydrological conditions that are increasingly influenced by climate change. A key unresolved question is how hydrologic flow, vegetation response, and rising sea levels interact to regulate soil redox condition and carbon loss in coastal wetlands. Using a field-tested hydrological–biogeochemical–ecological modeling framework, we reveal how the interplay between terrestrial groundwater discharge and tidal fluctuations generates complex groundwater flow patterns at the terrestrial–aquatic interface, and how these patterns modulate soil redox conditions, in turn influencing soil organic matter decomposition and carbon loss. Notably, rising sea levels suppress soil CO2 emissions while reducing lateral dissolved carbon losses, thereby enhancing litter carbon sequestration under anoxic conditions. As vegetation responds to sea-level rise and carbon inputs diminish, litter carbon subsequently declines. These findings underscore a critical hydrological control on carbon cycling, advancing our understanding of coastal ecosystem resilience in a warming world.

Chen, Kewei [ORNL] (ORCID:000000032580514X)↗
DOE OSTI2021

Gravitational effective field theory islands, low-spin dominance, and the four-graviton amplitude

Abstract We analyze constraints from perturbative unitarity and crossing on the leading contributions of higher-dimension operators to the four-graviton amplitude in four spacetime dimensions, including constraints that follow from distinct helicity configurations. We focus on the leading-order effect due to exchange by massive degrees of freedom which makes the amplitudes of interest infrared finite. In particular, we place a bound on the coefficient of the R 3 operator that corrects the graviton three-point amplitude in terms of the R 4 coefficient. To test the constraints we obtain nontrivial effective field-theory data by computing and taking the large-mass expansion of the one-loop minimally-coupled four-graviton amplitude with massive particles up to spin 2 circulating in the loop. Remarkably, we observe that the leading EFT coefficients obtained from both string and one-loop field-theory amplitudes lie in small islands. The shape and location of the islands can be derived from the dispersive representation for the Wilson coefficients using crossing and assuming that the lowest-spin spectral densities are the largest. Our analysis suggests that the Wilson coefficients of weakly-coupled gravitational physical theories are much more constrained than indicated by bounds arising from dispersive considerations of 2 → 2 scattering. The one-loop four-graviton amplitudes used to obtain the EFT data are computed using modern amplitude methods, including generalized unitarity, supersymmetric decompositions and the double copy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The Effect of Topology on Phase Behavior under Confinement

This work studies how morphology (i.e., the shape of a structure) and topology (i.e., how different structures are connected) influence wall adsorption and capillary condensation under tight confinement. Numerical simulations based on classical density functional theory (cDFT) are run for a wide variety of geometries using both hard-sphere and Lennard-Jones fluids. These cDFT computations are compared to results obtained using the Minkowski functionals. It is found that the Minkowski functionals can provide a good description of the behavior of Lennard-Jones fluids down to small system sizes. In addition, through decomposition of the free energy, the Minkowski functionals provide a good framework to better understand what are the dominant contributions to the phase behavior of a system. Lastly, while studying the phase envelope shift as a function of the Minkowski functionals it is found that topology has a different effect depending on whether the phase transition under consideration is a continuous or a discrete (first-order) transition.

42 ENGINEERING↗
DOE OSTI2022

B/PS bulges in DESI Legacy edge-on galaxies – I. Sample building

ABSTRACT We present the biggest up-to-date sample of edge-on galaxies with boxy/peanut-shaped (B/PS) bulges and X-structures. The sample was prepared using images from the Dark Energy Spectroscopic Instrument (DESI) Legacy catalogue and contains about 2000 galaxies. To find suitable candidates in catalogue, we made the assumption that the residues (original images minus model) of galaxies with B/PS bulges should exhibit a characteristic X-shape. Galaxies with such features were selected by eye and then used as input data for a neural network training, which was applied to a bigger sample of edge-on galaxies. Using the available data and the photometric models from the literature, we investigated the observational and statistical properties of the sample created. Comparing the B/D ratios for galaxies with and without B/PS bulges, we found that the B/D ratio for galaxies from our sample is statistically higher, with typical values in the range ≈0.2–0.5 depending on the decomposition procedure. We studied how the opening angles φ of the X-structure and the length of its rays are distributed in the formed sample and found them to be consistent with previous measurements and predictions from N-body models, e.g. $\varphi \gtrsim 25{^\circ }$, but measured here for a much larger number of galaxies. We found a sharp increase in the B/PS bulge fraction for stellar masses log M⋆ ≳ 10.4, but for edge-on galaxies, which complements the results of previous works. The sample can be used in future work to test various bar models and their relationship with B/PS bulges, as well as to study their stability and evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Effects of tensile loading during annealing of alnico melt spun ribbons

Conventional magnetic annealing (MA) of the permanent magnet alloy alnico involves application of an external magnetic field at temperatures within the spinodal decomposition range. This field biases the growth of the Fe-Co rich, ferromagnetic α 1 -phase in an energetically favorable 〈001〉 direction in alignment with the applied field within an Al-Ni rich, paramagnetic α 2 -phase. Utilizing a magnetic field to bias the α 1 -phase may limit alnico from reaching theoretical coercivity due to (1) the field having maximum biasing ability at temperatures near the Curie temperature where large α 1 -phase nanorods form and (2) connectivity of the α 1 -phase occurs unavoidably during MA. Both decrease the effective shape anisotropy of the α 1 -phase, thereby reducing coercivity. Herein, we explore tensile-loading as a biasing mechanism to control and optimize the final alnico nanostructure beyond that achieved by MA. Two samples of melt-spun alnico were heat-treated at 860 °C for 5 minutes: one sample was subjected to 10 MPa tensile stress for comparison with a stress-free control sample. Structural and magnetic characterization revealed that the stress-annealed ribbon sample possessed expected phase assemblages, but was distinguished by a ∼2× larger grain diameter and an elongated anisotropic α 1 -phase within grains that were oriented to a shear stress along 〈001〉 directions at an angle of ∼45° relative to the loading direction. Both types of annealing produced a similar increase in the coercivity and remanence, but a decrease in saturation magnetization.

Rinko, E. A.↗
DOE OSTI2021

Degradation Diagnostics from the Subsurface of Lithium-Ion Battery Electrodes

Despite the long-established rocking-chair theory of lithium-ion batteries (LIBs), developing novel characterization methodology with higher spatiotemporal resolution facilitates a better understanding of the solid electrolyte interphase studies to shape the reaction mechanisms. In this work, we develop a Xenon ion plasma focused ion beam (Xe + PFIB)-based characterization technique to probe the cross-sectional interface of both ternary cathode and graphite anode electrodes, with the focus on revealing the chemical composition and distribution underneath the electrode surface by in-depth analysis of secondary ions. Particularly, the lithium fluoride is detected in the pristine cathode prior to contact with the electrolyte, reflecting that the electrode degradation is in the form of the loss of lithium inventory during electrode preparation. This degradation is related to the hydrolysis of the cathode material and the decomposition of the PVDF binder. Here, through the quantitative analysis of the transition-metal degradation products, manganese is found to be the dominant element in the newly formed inactive fluoride deposition on the cathode, while no transition metal signal can be found inside the anode electrode. These insights at high resolution implemented via a PFIB-based characterization technique not only enrich the understanding of the degradation mechanism in the LIBs but also identify and enable a high-sensitivity methodology to obtain the chemical survey at the subsurface, which will help remove the capacity-fade observed in most LIBs.

25 ENERGY STORAGE↗