Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow forming”

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 649 records · Page 36

Jet energy scale and resolution measured in proton–proton collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

Jet energy scale and resolution measurements with their associated uncertainties are reported for jets using 36–81 fb -1 of proton–proton collision data with a centre-of-mass energy of √s=13 TeV collected by the ATLAS detector at the LHC. Jets are reconstructed using two different input types: topo-clusters formed from energy deposits in calorimeter cells, as well as an algorithmic combination of charged-particle tracks with those topo-clusters, referred to as the ATLAS particle-flow reconstruction method. The anti-k t jet algorithm with radius parameter R=0.4 is the primary jet definition used for both jet types. This result presents new jet energy scale and resolution measurements in the high pile-up conditions of late LHC Run 2 as well as a full calibration of particle-flow jets in ATLAS. Jets are initially calibrated using a sequence of simulation-based corrections. Next, several in situ techniques are employed to correct for differences between data and simulation and to measure the resolution of jets. The systematic uncertainties in the jet energy scale for central jets (|η|<1.2) vary from 1% for a wide range of high-p T jets (250 T <2000 GeV), to 5% at very low p T (20 GeV) and 3.5% at very high p T (>2.5 TeV). The relative jet energy resolution is measured and ranges from (24±1.5)% at 20 GeV to (6±0.5)% at 300 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multifidelity computing for coupling full and reduced order models

Hybrid physics-machine learning models are increasingly being used in simulations of transport processes. Many complex multiphysics systems relevant to scientific and engineering applications include multiple spatiotemporal scales and comprise a multifidelity problem sharing an interface between various formulations or heterogeneous computational entities. To this end, we present a robust hybrid analysis and modeling approach combining a physics-based full order model (FOM) and a data-driven reduced order model (ROM) to form the building blocks of an integrated approach among mixed fidelity descriptions toward predictive digital twin technologies. At the interface, we introduce a long short-term memory network to bridge these high and low-fidelity models in various forms of interfacial error correction or prolongation. The proposed interface learning approaches are tested as a new way to address ROM-FOM coupling problems solving nonlinear advection-diffusion flow situations with a bifidelity setup that captures the essence of a broad class of transport processes.

59 BASIC BIOLOGICAL SCIENCES↗

Modular Hybrid Solid State Transformer (H-SST) for Next Generation Flexible and Adaptable Large Power Transformer (LPT) (Final Technical Report)

A 500 kVA HSST, the world first of such system, has been successfully developed. The HSST has the capability to regulate the power flow of the transmission or distribution power lines, or provide voltage regulation. The HSST design is modular therefore it has the ability to be connected in series of parallel to form a large power transformer with smart controllability. The project has enabled a high TRL level prototype be demonstrated at University of Texas at Austin which can lead to additional commercialization and demonstration opportunities. A large number of students were trained under this program that will benefit the power grid industry. A number of publications are published at well known international conferences and journals, further accomplish the mission of knowledge development and dissemination.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Turbulence Modeling with Nek5000/RS, SOD2D and Alya

We present Validation and Verification (V & V) study of two high-order spectral element method (SEM) based Computational Fluid Dyanimcs (CFD) codes that we will utilize for work on turbulence modeling: Nek5000 and SOD2D. While the former solves the incompressible form of Navier-Stokes equation, the latter works with compressible set of equations and uses an entropy-viscosity formulation to account for the discontinuities for high Mach number flows. We demonstrate the accuracy of these codes for two benchmark sub-sonic turbulent flows: periodic channel and pipe flow, by carrying out first and second order statistical analysis, grid convergence and turbulent structure analysis using wall-resolved large eddy simulations (WRLES). Later, we report on the implementation and testing of various wall-modeling strategies for large eddy simulation of turbulent flows in Nek5000. These include both classical log-law based and decision-tree based machine-learning models. Accuracy of these closure strategies are analyzed and necessary future work is outlined.

97 MATHEMATICS AND COMPUTING↗

Pushing the Envelope-Moving Dynamic Building Envelope Thermal Energy Storage Systems Mainstream: Preprint

Buildings contribute to nearly 40% of the U.S. national energy consumption and a significant portion of CO2 emissions. More importantly, disadvantaged communities are disproportionately affected by energy burden and thermal discomfort in their homes. This paper will discuss two novel DOE's BTO supported thermal energy storage (TES) integrated dynamic building envelope technologies, their ability to harvest ambient energy, and improve energy efficiency by reducing HVAC loads and peak electricity demand while enhancing energy and thermal resilience in buildings. The paper will also discuss the recent advancements that have made these systems more affordable and easier to integrate into new and existing buildings. The first solution is a thermally anisotropic building envelope (TABE) system that can redirect ambient thermal energy (heat or coolness) from diurnal outdoor conditions, solar irradiance, and night sky radiation from the envelope to a hydronic loop. The redirected thermal energy can be stored in a TABE-integrated thermal energy storage system and use the stored energy to offset HVAC energy use and peak demand. The second solution is an innovative plug-and-play thermal switch in the form of insertable plugs integrated with a phase change material (PCM). The plug can vary its thermal resistance based on the indoor and outdoor conditions, thus allowing preferential directional heat flow, and enhancing utilization of free ambient cooling and heating to charge/discharge the PCM, much like a solid-state economizer. While the fist solution can be actively controlled, the second solution is passive, requiring no external power, and work solely based on the ambient temperature.

anisotropic envelope↗

Pushing the Envelope—Moving Dynamic Building Envelope Thermal Energy Storage Systems Mainstream

Buildings contribute to nearly 40% of the U.S. national energy consumption and a significant portion of CO2 emissions. More importantly, disadvantaged communities are disproportionately affected by energy burden and thermal discomfort in their homes. This paper will discuss two novel DOE's BTO supported thermal energy storage (TES) integrated dynamic building envelope technologies, their ability to harvest ambient energy, and improve energy efficiency by reducing HVAC loads and peak electricity demand while enhancing energy and thermal resilience in buildings. The paper will also discuss the recent advancements that have made these systems more affordable and easier to integrate into new and existing buildings. The first solution is a thermally anisotropic building envelope (TABE) system that can redirect ambient thermal energy (heat or coolness) from diurnal outdoor conditions, solar irradiance, and night sky radiation from the envelope to a hydronic loop. The redirected thermal energy can be stored in a TABE-integrated thermal energy storage system and use the stored energy to offset HVAC energy use and peak demand. The second solution is an innovative plug-and-play thermal switch in the form of insertable plugs integrated with a phase change material (PCM). The plug can vary its thermal resistance based on the indoor and outdoor conditions, thus allowing preferential directional heat flow, and enhancing utilization of free ambient cooling and heating to charge/discharge the PCM, much like a solid-state economizer. While the fist solution can be actively controlled, the second solution is passive, requiring no external power, and work solely based on the ambient temperature.

Shrestha, Som↗

On the Flow of a Cement Suspension: The Effects of Nano-Silica and Fly Ash Particles

Additives such as nano-silica and fly ash are widely used in cement and concrete materials to improve the rheology of fresh cement and concrete and the performance of hardened materials and increase the sustainability of the cement and concrete industry by reducing the usage of Portland cement. Therefore, it is important to study the effect of these additives on the rheological behavior of fresh cement. In this paper, we study the pulsating Poiseuille flow of fresh cement in a horizontal pipe by considering two different additives and when they are combined (nano-silica, fly ash, combined nano-silica, and fly ash). To model the fresh cement suspension, we used a modified form of the power-law model to demonstrate the dependency of the cement viscosity on the shear rate and volume fraction of cement and the additive particles. The convection–diffusion equation was used to solve for the volume fraction. After solving the equations in the dimensionless forms, we conducted a parametric study to analyze the effects of nano-silica, fly ash, and combined nano-silica and fly ash additives on the velocity and volume fraction profiles of the cement suspension. According to the parametric study presented here, larger nano-silica content results in lower centerline velocity of the cement suspension and larger non-uniformity of the volume fraction. Compared to nano-silica, fly ash exhibits an opposite effect on the velocity. Larger fly ash content results in higher centerline velocity, while the effect of the fly ash on the volume fraction is not obvious. For cement suspension containing combined nano-silica and fly ash additives, nano-silica plays a dominant role in the flow behavior of the suspension. The findings of the study can help the design and operation of the pulsating flow of fresh cement mortars and concrete in the 3D printing industry.

36 MATERIALS SCIENCE↗

Convectively Forced Diurnal Gravity Waves in the Maritime Continent

Long-lived, zonally propagating diurnal rainfall disturbances are a highly pronounced and common feature in the Maritime Continent (MC). A recent study argues that these disturbances can be explained as diurnally phase-locked gravity waves. Here we explore the origins of these waves through regional cloud-permitting numerical model experiments. The gravity waves are reproduced and isolated in the model framework through the combined use of realistic geography and diurnally cyclic lateral boundary conditions representative of both characteristic easterly and westerly background zonal flow regimes. These flow regimes are characteristic of the Madden–Julian oscillation (MJO) suppressed and active phase in the MC, respectively. Tests are conducted wherein Borneo, Sumatra, or both islands and/or their orography are removed. These tests imply that the diurnal gravity waves are excited and maintained directly by latent heating from the vigorous mesoscale convective systems (MCSs) that form nocturnally in both Borneo and Sumatra. Removing orography has only a secondary impact on both the MCSs and the gravity waves, implying that it is not critical to these waves. We therefore hypothesize that diurnal gravity waves are fundamentally driven by mesoscale organized deep convection, and are only sensitive to orography to the measure that the convection is affected by the orography and its mesoscale flows. Factor separation further reveals that the nonlinear interaction of synchronized diurnal cycles in Sumatra and Borneo slightly amplifies this gravity wave mode compared to if either island existed in isolation. This nonlinear feedback appears most prominently at longitudes directly between the two islands.

Ruppert, Jr., James H.↗

Nanopore networks in colloidal silica assemblies characterized by XCT for confined fluid flow modeling

Fluid flow through nanopore structures has exhibited different behavior than that described by Darcy's law, derived from pore networks of micrometer scale. Here, pressure assemblies of colloidal silica spheres of specified diameters from 94 to 620 nanometers were prepared as model nano-porous media. The smaller silica spheres tended to form random packing, but some substantial ordering was found in the packing of the larger silica spheres. X-ray computed tomography (XCT), with a voxel resolution of 16 nm, was used to characterize the nanopore networks. The porosity of the Nano-XCT pore network (34%) was less than the total porosity of the assemblies as determined by X-ray attenuation (52%), possibly due to the limits of voxel resolution during segmentation. Fluid flow in the nanopore networks was simulated using the single-phase Lattice Boltzmann Method, and the simulated permeability was compared with the empirical and experimental values. Based on our characterization, although a well-ordered packing of silica spheres was not achieved, it was found that the nanopore networks in the colloidal silica assemblies had pore size distributions corresponding to the particle sizes. The simulated permeability was less than the experimental measurement for water flow, but the complex packing of silica spheres and surface chemistry issues need to be considered in future research.

02 PETROLEUM↗

Bayesian Physics Informed Spatio-Temporal Network for Streamflow Data Imputation

Reliable reconstruction of incomplete streamflow records is critical for improving hydrological forecasting, flood preparedness, and water resource management. However, large observational gaps and uncertainties in governing physical parameters limit the accuracy of traditional statistical and machinelearning imputation frameworks. To address these challenges, we develop a Bayesian Physics-Informed Spatio-Temporal Network (BPI-STNet) that jointly captures spatial and temporal dependencies while enforcing hydrologic consistency through embedded physical constraints. The framework integrates a GraphSAGE-LSTM architecture to model spatial connectivity across gauges and temporal flow dynamics, coupled with a Bayesian update mechanism to estimate uncertain parameters in a simplified water-balance framework. Unlike conventional physics-informed networks that rely on sampling-based posterior estimation, BPI-STNet derives an analytic solution to the inverse problem, allowing closed-form Bayesian updates of uncertain parameters Λ={α,β,k} using Gaussian priors and likelihoods. Applied to daily observations from the Susquehanna River Basin (1980-2022), BPI-STNet achieves substantial improvements over a purely data-driven RGNN baseline, which reduced RMSE by 23 % and MAE by 9 %, and achieving an average NSE values up to 0.96. The results demonstrate that coupling Bayesian inference with physics-informed learning yields physically consistent, uncertainty-aware reconstructions that preserve the temporal persistence and statistical distribution of observed flows. The proposed framework establishes a generalizable paradigm for data-sparse hydrologic systems where both data fidelity and physical interpretability are essential.

Krishnan Kutty Ambika, Anukesh [ORNL] (ORCID:00000↗

Probing triaxial deformation of atomic nuclei in high-energy heavy ion collisions

Most atomic nuclei are deformed with a quadrupole shape described by its overall strength β 2 and triaxiality γ. The deformation can be accessed in high-energy heavy ion collisions by measuring the collective flow response of the produced quark-gluon plasma to the eccentricity ε 2 and the density gradient d ⊥ in the initial state. Using an analytical estimate and a Glauber model, I show that the variances $\langleε^{2}_{2}\rangle$ or $\langle$(δd ⊥ /d ⊥ ) 2 $\rangle$ and skewnesses $\langleε^{2}_{2}δd_{⊥}/d_{⊥}\rangle$ or $\langle(δd_{⊥}/d_{⊥})^{3}\rangle$ have a simple analytical form of a'+b'β$^{2}_{2}$ and a'+[b'+c' cos(3γ)]β$^{3}_{2}$, respectively. From these, I constructed several normalized skewnesses to isolate the γ dependence from that of β 2 and show that the correlations between a normalized skewness and a variance can constrain simultaneously β 2 and γ. Assuming a linear relation with elliptic flow v 2 and mean-transverse momentum [p T ] of final-state particles, v 2 ∝ ε 2 and δ [p T ] / [p T ] ∝ δd ⊥ /d ⊥ , similar conclusions are also expected for the variances and skewnesses of v 2 and [p T ], i.e., a + bβ$^{2}_{2}$ for $\langle v^{2}_{2}\rangle$ and $\langle(δ[p_{T}] / [p_{T}])^{2}\rangle$ and a + [b + c cos (3γ)] β$^{3}_{2}$ for $\langle v^{2}_{2}δ [p_{T}] / [p_{T}]\rangle$ or $\langle(δ[p_{T}] / [p_{T}])^{3}\rangle$. My findings motivate a dedicated system scan of high-energy heavy ion collisions at RHIC and LHC to measure triaxiality of atomic nuclei: one first determines the coefficients b and c by collisions of isobaric near prolate nuclei, cos (3γ) ≈ 1, and near oblate nuclei, cos (3γ) ≈ –1 , with known β 2 values, followed by collisions of other species of interest with similar mass number. The (β 2 , γ) values for this species can be inferred directly from the measured variance and skewness observables from these collisions. Finally, the results demonstrate the unique opportunities offered by high-energy collisions as a tool to perform interdisciplinary nuclear physics studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Alteration of phycobilisome excitation energy transfer properties in response to attenuations in peripheral electron flow

In Synechocystis sp. PCC 6803 ( S . 6803), two types of phycobilisome (PBS) complexes, CpcG-PBS and CpcL-PBS, function to harvest light energy for photosynthetic reaction centers (RCs), photosystem I (PSI) and photosystem II (PSII). The compositional differences between these two forms of PBS and their specificity for RCs have led to suggestions that they may differ in function. To address this question, we examined how PBS-RC interactions, and the transfer of excitation energy from PBS to RCs, might be adjusted under conditions where electron demand and photon availability are modulated. The CpcG-PBS, CpcL-PBS, and RC complexes were isolated from a S . 6803 strain defective in expression of flavodiiron 1 (oxygen reduction reaction 1, ORR1) grown under varied light regimes. The energy transfer preference from CpcL-PBS to either PSI or PSII was investigated by in vitro crosslinking and 77 K fluorescence emission spectroscopy to assess energy transfer efficiency under photoexcitation. While the results demonstrate that the transfer of excitation energy from CpcL-PBS favors PSI over PSII in WT strains as previously shown, the preference of CpcL-PBS switches from PSI to PSII in ORR1 strains. Surprisingly, this change in preference was reproduced when ORR1 CpcL-PBS was crosslinked with WT RCs, or when WT CpcL-PBS was cross-crosslinked with ORR1 RCs, indicating there are physical modifications to both PBS and RCs that mediate the preference switch. In contrast, the analysis with ORR1 CpcG-PBS shows similar preferences to WT. Additionally, PBS populations in ORR1 shifted to a greater proportion of CpcL-PBS relative to CpcG-PBS. These results demonstrate that under conditions where electron utilization changes, there is a tuning of the excitation energy allocation from CpcL-PBS to RCs to manage the energy distribution for photosynthesis under dynamic flux conditions.

59 BASIC BIOLOGICAL SCIENCES↗

Protecting Customer Privacy Through Distributed Energy Resource Anonymization

Due to their stochastic nature, the increase of Renewable Energy Resources (RERs) as a primary source of energy for power grids creates challenges regarding the reliability and resilience of the system. In order to combat these obstacles, expansion of Distributed Energy Resources (DERs) and their participation in Demand Response (DR) programs is necessary. Widespread participation requires prioritizing customer privacy and addressing concerns that may arise regarding communication between DERs and the Grid Service Provider (GSP). This paper discusses the use of flow reservation resources to split the operating cycles of DER load profiles into unique phases. The splitting of phases increases anonymization of the DERs by making it more difficult to determine the individual characteristics of the device. We discuss an example of this using simulated DER load profile data and examine the resulting effectiveness by using a machine learning algorithm for classification, called Support Vector Machine (SVM).

Distributed Energy Resource, Anonymization, Renewa↗

Relational bulk reconstruction from modular flow

Abstract The entanglement wedge reconstruction paradigm in AdS/CFT states that for a bulk qudit within the entanglement wedge of a boundary subregion$$ \overline{A} $$ A ¯ , operators acting on the bulk qudit can be reconstructed as CFT operators on$$ \overline{A} $$ A ¯ . This naturally fits within the framework of quantum error correction, with the CFT states containing the bulk qudit forming a code protected against the erasure of the boundary subregionA. In this paper, we set up and study a framework for relational bulk reconstruction in holography: given two code subspaces both protected against erasure of the boundary regionA, the goal is to relate the operator reconstructions between the two spaces. To accomplish this, we assume that the two code subspaces are smoothly connected by a one-parameter family of codes all protected against the erasure ofA, and that the maximally-entangled states on these codes are all full-rank. We argue that such code subspaces can naturally be constructed in holography in a “measurement-based” setting. In this setting, we derive a flow equation for the operator reconstruction of a fixed code subspace operator using modular theory which can, in principle, be integrated to relate the reconstructed operators all along the flow. We observe a striking resemblance between our formulas for relational bulk reconstruction and the infinite-time limit of Connes cocycle flow, and take some steps towards making this connection more rigorous. We also provide alternative derivations of our reconstruction formulas in terms of a canonical reconstruction map we call the modular reflection operator.

Physics↗

Modifying the resolving cysteine affects the structure and hydrogen peroxide reactivity of peroxiredoxin 2

Peroxiredoxin 2 (Prdx2) is a thiol peroxidase with an active site Cys (C52) that reacts rapidly with H 2 O 2 and other peroxides. The sulfenic acid product condenses with the resolving Cys (C172) to form a disulfide which is recycled by thioredoxin or GSH via mixed disulfide intermediates or undergoes hyperoxidation to the sulfinic acid. C172 lies near the C terminus, outside the active site. It is not established whether structural changes in this region, such as mixed disulfide formation, affect H 2 O 2 reactivity. To investigate, we designed mutants to cause minimal (C172S) or substantial (C172D and C172W) structural disruption. Stopped flow kinetics and mass spectrometry showed that mutation to Ser had minimal effect on rates of oxidation and hyperoxidation, whereas Asp and Trp decreased both by 100-fold. To relate to structural changes, we solved the crystal structures of reduced WT and C172S Prdx2. The WT structure is highly similar to that of the published hyperoxidized form. C172S is closely related but more flexible and as demonstrated by size exclusion chromatography and analytical ultracentrifugation, a weaker decamer. Size exclusion chromatography and analytical ultracentrifugation showed that the C172D and C172W mutants are also weaker decamers than WT, and small-angle X-ray scattering analysis indicated greater flexibility with partially unstructured regions consistent with C-terminal unfolding. We propose that these structural changes around C172 negatively impact the active site geometry to decrease reactivity with H 2 O 2 . This is relevant for Prdx turnover as intermediate mixed disulfides with C172 would also be disruptive and could potentially react with peroxides before resolution is complete.

59 BASIC BIOLOGICAL SCIENCES↗

Convex Programs for Minimal-Area Problems

The minimal-area problem that defines string diagrams in closed string field theory asks for the metric of least area on a Riemann surface with the condition that all non-contractible closed curves have length at least 2π. This is an extremal length problem in conformal geometry as well as a problem in systolic geometry. Here, we consider the analogous minimal-area problem for homology classes of curves and, with the aid of calibrations and the max flow-min cut theorem, formulate it as a local convex program. We derive an equivalent dual program involving maximization of a concave functional. These two programs give new insights into the form of the minimal-area metric and are amenable to numerical solution. We explain how the homology problem can be modified to provide the solution to the original homotopy problem.

97 MATHEMATICS AND COMPUTING↗

Co-existence of atomically dispersed Ru and Ce 3+ sites is responsible for excellent low temperature N 2 O reduction activity of Ru/CeO 2

Nitrous oxide N 2 O reduction is a big challenge due to high global warming potential of N 2 O. (~300 times higher compared with CO 2 ). The best known catalysts, such as Rh/ceria, require relatively high temperatures for N 2 O decomposition. Herein, we report that Ru/ceria catalysts with low Ru loading of ~0.25 wt% efficiently catalyze low temperature N 2 O reduction by CO starting at 100 °C (full N 2 O conversion below 200 °C) under industrially relevant flow rates and gas concentrations. Further, this remarkable performance stems from maintaining isolated Ru cations even on reduced ceria surface and, simultaneously, the propensity of Ru to affect ceria surface to form labile surface oxygen thereby creating large number of oxygen vacancies (Ce +3 cations) in the presence of CO. In contrast, for Rh/CeO 2 catalysts with equivalent metal loading, the activity is much lower because atomically dispersed Rh sinters into metallic clusters at the onset reaction temperature (~200 °C): these clusters are much less effective than isolated single Ru ions, with lower Ce +3 concentration maintained on reduced Rh/CeO 2 catalyst. Our study highlights the benefits of gaining molecular-level insight into the dynamic nature of catalytically active sites under reaction conditions for preparing catalysts containing low loading of precious metals with unsurpassed low temperature activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The value of in-reservoir energy storage for flexible dispatch of geothermal power

Geothermal systems making use of advanced drilling and well stimulation techniques have the potential to provide tens to hundreds of gigawatts of clean electricity generation in the United States by 2050. With near-zero variable costs, geothermal plants have traditionally been envisioned as providing “baseload” power, generating at their maximum rated output at all times. However, as variable renewable energy sources (VREs) see greater deployment in energy markets, baseload power is becoming increasingly less competitive relative to flexible, dispatchable generation and energy storage. Herein we conduct an analysis of the potential for future geothermal plants to provide both of these services, taking advantage of the natural properties of confined, engineered geothermal reservoirs to store energy in the form of accumulated, pressurized geofluid and provide flexible load-following generation. We develop a linear optimization model based on multi-physics reservoir simulations that captures the transient pressure and flow behaviors within a confined, engineered geothermal reservoir. We then optimize the investment decisions and hourly operations of a power plant exploiting such a reservoir against a set of historical and modeled future electricity price series. Further, we find that operational flexibility and in-reservoir energy storage can significantly enhance the value of geothermal plants in markets with high VRE penetration, with energy value improvements of up to 60% relative to conventional baseload plants operating under identical conditions. Across a range of realistic subsurface and operational conditions, our modeling demonstrates that confined, engineered geothermal reservoirs can provide large and effectively free energy storage capacity, with round-trip storage efficiencies comparable to those of leading grid-scale energy storage technologies. Optimized operational strategies indicate that flexible geothermal plants can provide both short- and long-duration energy storage, prioritizing output during periods of high electricity prices. Sensitivity analysis assesses the variation in outcomes across a range of subsurface conditions and cost scenarios.

15 GEOTHERMAL ENERGY↗