SEARCH · Engineering Papers
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.
This dataset was used to support machine learning model development for wildfire impacts on utility-scale solar in the United States. The dataset includes information regarding energy generation, PM2.5, clearness index, temperature, wind speed, precipitation, and site size.
As the prospect of exceeding global temperature targets set forth in the Paris Agreement becomes more likely, methods of climate intervention are increasingly being explored. With this increased interest there is a need for an assessment process to understand the range of impacts across different scenarios against a set of performance goals in order to support policy decisions. The methodology and tools developed for Performance Assessment (PA) for nuclear waste repositories shares many similarities with the needs and requirements for a framework for climate intervention. Using PA, we outline and test an evaluation framework for climate intervention, called Performance Assessment for Climate Intervention (PACI) with a focus on Stratospheric Aerosol Injection (SAI). We define a set of key technical components for the example PACI framework which include identifying performance goals, the extent of the system, and identifying which features, events, and processes are relevant and impactful to calculating model output for the system given the performance goals. Having identified a set of performance goals, the performance of the system, including uncertainty, can then be evaluated against these goals. Using the Geoengineering Large Ensemble (GLENS) scenario, we develop a set of performance goals for monthly temperature, precipitation, drought index, soil water, solar flux, and surface runoff. The assessment assumes that targets may be framed in the context of risk-risk via a risk ratio, or the ratio of the risk of exceeding the performance goal for the SAI scenario against the risk of exceeding the performance goal for the emissions scenario. From regional responses, across multiple climate variables, it is then possible to assess which pathway carries lower risk relative to the goals. The assessment is not comprehensive but rather a demonstration of the evaluation of an SAI scenario. Future work is needed to develop a more complete assessment that would provide additional simulations to cover parametric and aleatory uncertainty and enable a deeper understanding of impacts, informed scenario selection, and allow further refinements to the approach.
Ridesharing has become an increasingly popular transportation method over the past decade. Transportation network companies such as Uber and Lyft generally provide two types of rideshare services: personal rideshare, in which users ride alone or with individuals they know, and pooled rideshare, in which users ride with passengers they do not know but share similar routes. Pooled rideshare is capable of reducing energy consumption and traffic in the transportation system in comparison to personal rideshare. Despite the growth in trip volume, ridesharing usage is still low compared to other popular transportation methods in the U.S., particularly traveling in one’s own personal vehicle. Furthermore, pooled rideshare usage is lower than personal rideshare. To understand riders’ preferences, a national survey (N = 2884) was conducted in the U.S. to investigate users’ choice behaviors in rideshare services examining personal versus pooled rideshare. Each survey respondent completed 20 stated-preference scenarios where participants choose between a personal or pooled rideshare option. Based on the responses, a mixed logit model was developed to capture the choice behavior preferences of the participants. The model unveiled the impact of demographic and trip attribute variables on users’ rideshare preferences. The discussion encompassed insights into demographic backgrounds and trip attributes, accompanied by a set of policy recommendations aimed at enhancing future pooled rideshare utilization.
Reliable simulation, prediction, and complete theoretical understanding of atmospheric blocking remain challenging despite its significant socio-economic impacts. Generations of climate models have notoriously underestimated blocking frequency, particularly over the Euro-Atlantic sector. Identifying factors controlling blocking frequency and dynamics is therefore essential for improving its simulation. Here, using a cloud-locking experiment, we show that cloud radiative effects (CREs) significantly increase the frequency of Euro-Atlantic blocking. CREs enhance upstream diabatic source of wave activity, both directly through longwave heating and indirectly through their feedback on latent heating, with the latter playing the dominant role. The resulting increase in the upstream diabatic source feeds into local wave activity downstream and promotes blocking formation. Qualitatively similar results are shown by multi-model experiments with radiatively inactive clouds to longwave radiation, albeit with a larger impact from mean-state changes. The results underscore the necessity of accurately representing cloud-radiation interactions in weather and climate models for improved prediction of blocking events.
The rate of sediment supply has significant impacts on river morphology, making it crucial to understand the geomorphic changes and grain size distribution dynamics in rivers. However, the effects of varying grain size sediment input on morphological changes in braided channels remain poorly understood. This study is the first to investigate the bar development and sediment sorting processes in braided channels with non-uniform sediment inputs using both numerical and experimental approaches. We applied a two-dimensional numerical model, Nays2DH to confirm and generalize experimental results. The model reproduced key experiment results, including 1) stream elevation changes, and 2) grain size distribution. Using this validated model, we explored the morphological changes and sorting process in a braided river with sediment inputs. The numerical experiments demonstrate that sediment input controls the elevation of the stream bed and the grain size distribution. Notably, both the elevation and grain-size distribution become relatively stable in downstream of the channel. Additionally, the simulation results suggest that an increased sediment supply leads to greater channel complexity, with bed surface armoring decreasing.
In 2018, the U.S. Department of Energy (DOE) held a workshop for the Atmospheric Radiation Measurement (ARM) (Mather and Voyles 2013) user facility to discuss critical climate challenges and locations where key ARM Mobile Facility (AMF) observational assets could impact Earth system modeling (ESM). As an outcome, the southeast United States (SE U.S.) was identified as a high-priority region to target climate-process studies that promote a deeper understanding of the climate system and bolster ARM interactions with the community to drive ESM advancement. The DOE ARM user facility is a globally recognized leader in deploying and operating strategically located observation sites around the world for studying the properties of aerosols and clouds and their interaction with radiation, precipitation, and the Earth’s surface. In partnering with the DOE Atmospheric System Research (ASR) program, ARM solicited a multi-agency Site Science Team approach to provide input and close interaction with ARM management towards a successful SE U.S. deployment of the ARM third Mobile Facility (AMF3) (Miller et al. 2016). These efforts included identifying key locations, science drivers and instruments, and measurement strategies to address the wider climate-process needs and ESM improvement. Community input served a vital role in establishing, refining, and informing the relevant drivers and decisions regarding this AMF3 deployment. The team has identified Northern Alabama (N. AL) as regionally representative to unlock the key opportunities that will improve our understanding and model representation of aerosol, cloud, and land surface processes and their couplings in the SE U.S. A defining aspect of the AMF3 deployment is its commitment to long-term (anticipated five-year) observations to mitigate potential seasonal-to-annual variability that often limits appropriate attribution of phenomena to local or larger-scale processes. The proposed location may leverage nearby surface networks and multi-agency and partner assets to enrich this multi-year deployment. One motivation is to understand the role of spatiotemporal variability (thermodynamic, land-surface) across aspects of the climate system, with our AMF3 team anticipating future demands on characterizing the relationships between local-to-regional cloud development and surface processes across a diverse patchwork of natural, managed, and urban landscapes as found throughout the N. AL regions. The main site targets an intact, representative, forested region – the Bankhead National Forest (BNF) – underscoring further team commitment to regionally important land atmosphere two-way interactive studies “from the canopy to the clouds”, with enhanced tower instrumentation augmenting traditional ARM capabilities adjacent to this site. Multiple supplemental sites will also be distributed across this region, prioritizing added needs for biodiversity. Anticipated high-priority cloud science themes will target N. AL as a regional SE U.S. hotbed for high-impact weather, convective cloud onset, and shallow to-deep cloud transitioning. Anticipated aerosol drivers will focus on chemical processes that control the evolution of organic aerosol, the seasonality and spatial distribution of water vapor and particle-phase water, and its role on aerosol optical properties. Anticipated land atmosphere drivers consider the two-way feedbacks between surface influence on aerosols, clouds, and precipitation properties and the associated radiative impacts on plant physiology and canopy-scale fluxes. Emphasis will include the study of the impact of surface processes on aerosols via precursor emission, and on clouds via moisture flux and thermal development.
In 2018, the U.S. Department of Energy (DOE) held a workshop for the Atmospheric Radiation Measurement (ARM) (Mather and Voyles 2013) user facility to discuss critical climate challenges and locations where key ARM Mobile Facility (AMF) observational assets could impact Earth system modeling (ESM). As an outcome, the southeast United States (SE U.S.) was identified as a high-priority region to target climate-process studies that promote a deeper understanding of the climate system and bolster ARM interactions with the community to drive ESM advancement. The DOE ARM user facility is a globally recognized leader in deploying and operating strategically located observation sites around the world for studying the properties of aerosols and clouds and their interaction with radiation, precipitation, and the Earth’s surface. In partnering with the DOE Atmospheric System Research (ASR) program, ARM solicited a multi-agency Site Science Team approach to provide input and close interaction with ARM management towards a successful SE U.S. deployment of the ARM third Mobile Facility (AMF3) (Miller et al. 2016). These efforts included identifying key locations, science drivers and instruments, and measurement strategies to address the wider climate-process needs and ESM improvement. Community input served a vital role in establishing, refining, and informing the relevant drivers and decisions regarding this AMF3 deployment. The team has identified Northern Alabama (N. AL) as regionally representative to unlock the key opportunities that will improve our understanding and model representation of aerosol, cloud, and land-surface processes and their couplings in the SE U.S. A defining aspect of the AMF3 deployment is its commitment to long-term (anticipated five-year) observations to mitigate potential seasonal-to-annual variability that often limits appropriate attribution of phenomena to local or larger-scale processes. The proposed location may leverage nearby surface networks and multi-agency and partner assets to enrich this multi-year deployment. One motivation is to understand the role of spatiotemporal variability (thermodynamic, land-surface) across aspects of the climate system, with our AMF3 team anticipating future demands on characterizing the relationships between local-to-regional cloud development and surface processes across a diverse patchwork of natural, managed, and urban landscapes as found throughout the N. AL regions. The main site targets an intact, representative, forested region – the Bankhead National Forest (BNF) – underscoring further team commitment to regionally important land-atmosphere two-way interactive studies “from the canopy to the clouds”, with enhanced tower instrumentation augmenting traditional ARM capabilities adjacent to this site. Multiple supplemental sites will also be distributed across this region, prioritizing added needs for biodiversity. Anticipated high-priority cloud science themes will target N. AL as a regional SE U.S. hotbed for high-impact weather, convective cloud onset, and shallow-to-deep cloud transitioning. Anticipated aerosol drivers will focus on chemical processes that control the evolution of organic aerosol, the seasonality and spatial distribution of water vapor and particle-phase water, and its role on aerosol optical properties. Anticipated land-atmosphere drivers consider the two-way feedbacks between surface influence on aerosols, clouds, and precipitation properties and the associated radiative impacts on plant physiology and canopy-scale fluxes. Emphasis will include the study of the impact of surface processes on aerosols via precursor emission, and on clouds via moisture flux and thermal development.
Fire emissions from the Maritime Continent (MC) over the western tropical Pacific are strongly influenced by El Niño–Southern Oscillation (ENSO), posing various climate effects to the Earth system. Here, in this study, we show that the historical biomass burning emissions of black carbon (BC bb ) aerosol in the dry season from the MC are strengthened in El Niño years due to the dry conditions. The eastern Pacific type of El Niño exerts a stronger modulation in BC bb emissions over the MC region than the central Pacific type of El Niño. Based on simulations using the fully coupled Community Earth System Model (CESM), the impacts of increased BC bb emissions on ENSO variability and frequency are also investigated in this study. With BC bb emissions from the MC scaled up by a factor of 10, which enables the identification of climate response from the internal variability, the increased BC bb heats the local atmosphere and changes land–sea thermal contrast, which suppresses the westward transport of the eastern Pacific surface water. It leads to an increase in sea surface temperature in the eastern tropical Pacific, which further enhances ENSO variability and increases the frequency of extreme El Niño and La Niña events. This study highlights the potential role of BC bb emissions on extreme ENSO frequency, and this role may be increasingly important in the warming future with higher wildfire risks.
Many nations responded to the COVID-19 pandemic by restricting travel and other activities during 2020, resulting in temporarily reduced emissions of CO 2 , other greenhouse gases and ozone and aerosol precursors. We perform a coordinated Intercomparison, CovidMIP, of Earth System model simulations to assess the impact on climate of these emissions reductions. Eleven models performed multiple initial-condition ensembles to produce over 280 simulations spanning both initial condition and model structural uncertainty. We find model consensus on reduced aerosol amounts (particularly over East Asia) and associated increases in surface shortwave radiation levels. However, any impact on near-surface temperature or rainfall during 2020-2024 is extremely small and is not detectable in this initial analysis. Regional analyses on a finer scale, and closer attention to extremes (especially linked to changes in atmospheric composition and air quality) are required to test the impact of COVID-19-related emission reductions on near-term climate.
The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.
During condensation of water on a superhydrophobic nanostructured surfaces, droplets coalesce and jump (droplet ejection) due to the release of excess surface energy. Meticulously designed nanostructured surfaces or coatings capable of droplet ejection can enhance heat transfer performance by easy removal of droplets during condensation. In the present work, the capabilities of the Nelumbo droplet ejection coatings technology currently used in commercial heat exchangers are explored and optimized for materials and conditions relevant to steam condensers. Specifically, this includes (1) adapting and optimizing the procedures and chemistries to deposit droplet ejecting coatings on materials used in the steam condenser, (2) testing of the heat transfer, durability, and fouling of the fabricated coatings under a variety of steam condenser relevant conditions, and (3) modeling of the impact of droplet dynamics and interfacial properties on heat transfer in steam condensers. We experimentally demonstrated that at low supersaturation conditions (S~1.02), developed superhydrophobic nanostructured surface gives a 40% higher peak heat transfer flux and a 17% higher peak overall heat transfer coefficient (U) with droplet ejection mechanism compared to the dropwise condensation mechanism on the same nanostructured surface. Though the experimentally observed droplet ejection mechanism was short-lived (~3 min) due to the limitation of the chiller to maintain a constant supersaturation condition in fabricated steam condenser. And large variations in water inlet temperature (+ to - 2 degrees C) caused flooding of nanostructure (Wenzel state) at higher supersaturation (S > or = 1.12). This is an important finding because typical supersaturation levels in field operation are over 1.8 and as such, either the flooding potential for these materials should be improved or operation at lower supersaturation may be required to maximize heat transfer efficiency. If supersaturation conditions can be maintained reliably in the steam condenser in addition to droplet ejection mechanism during condensation, this could provide a potential 2% reduction in water flow rate which equates to a savings of over 3900 gallons of water per minute for a 500 MW steam turbine and over $0.3MM savings annually given a 0.02 cent marginal water withdrawal cost. It is also worth noting that these savings are based on improvement of droplet ejection coatings over coatings capable of performing consistent dropwise mechanism in in-field surface condensers. Savings due to the droplet ejections coatings when compared to bare substrates could be much higher. Additionally, the present work provides the importance of steam condenser operating conditions and insights into the challenges in modeling the condensing droplet dynamics on surfaces.
Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.
With the surge in electric vehicle (EV) adoption and the need for extended driving ranges, optimizing energy efficiency, particularly through thermal management, is critical, especially in extreme weather. Managing the substantial energy needed for cabin climate control and battery temperature regulation can increase energy demands by over 50 %, severely limiting range. This study conducts a comparative analysis of thermal management systems (TMS) in three popular EV vehicles, 2020 Chevrolet Bolt, 2019 Nissan Leaf Plus, and 2020 Tesla Model 3, evaluating their distinct TMS configurations and performance under varied weather conditions. Using both numerical simulations and experimental data collected on a controlled test bench at Argonne National Laboratory, we assess how TMS architecture and operational modes influence energy consumption and range. A comprehensive TMS model was developed, integrating cabin and battery thermal sub-models in the Autonomie software platform, to simulate temperature fluctuations and range impacts. Cabin climate was modeled using a mono-zonal approach, while battery cell temperature distribution was estimated through a 2D nodal structure. Each vehicle's distinct TMS setup was evaluated: the Chevrolet Bolt and Tesla Model 3 use a dual evaporator vapor compression cycle with a PTC heater for the cabin and a coolant loop for battery thermal management; the Nissan Leaf Plus employs a heat pump with a PTC heater for the cabin and air-cooling for the battery. Tests conducted at ambient temperatures of 35°C, 22°C, -7°C, and -18°C reveal significant differences in energy use and range reduction across both configurations and conditions. At 35°C, the Tesla Model 3, Chevrolet Bolt, and Nissan Leaf Plus have a range reduction of 8%, 9%, and 13%, respectively, due to air conditioning. In winter, heating technology is paramount; at -7°C, the Nissan Leaf's heat pump configuration achieves a lower range reduction (19.3%) compared to the Tesla and Chevrolet Bolt PTC heaters, which reduce range by 28.3% and 31%, respectively. Further, this study provides valuable insights for automotive engineers, EV technology researchers, and thermal management system designers aiming to enhance electric vehicle performance by understanding how different weather conditions and TMS architectures impact energy consumption and driving range.
Agriculture yields are susceptible to changes in future temperature, precipitation, and other Earth system factors. Future changes to these physical Earth system attributes and their effects on agricultural yields are highly uncertain. United States agricultural producers will be affected by such changes whether they occur domestically or internationally (via international commodity markets). How important to U.S. agriculture are domestic changes due to climate compared to those occurring externally? Here we show that potential direct impacts to United States agriculture have financial impacts on U.S. producers that are roughly similar in magnitude but opposite in sign to impacts on U.S. producers due to changes in other parts of the world. This finding is robust across application of a wide range of potential future crop yield impacts applied to a multi-sector global integrated assessment model. Therefore, when examining country-specific impacts, it is necessary to model both the domestic and international impacts.
The use of blended fuel sources in land based gas turbine engines drives variations in the resulting operational profile (temperatures and pressures) which can impact engine reliability. Furthermore, variability in the manufacture of components affects the resulting microstructure which directly impacts material performance and reliability. Currently, data-driven models are typically used for maintaining and inspecting fleets of engines. Without explicitly capturing material and operational sources of variability conservatism must be used in developing component-level reliability models. Therefore, there exists an opportunity to use information from materials-scale physics models to better inform reliability modeling and reduce conservatism; the impact is more cost-efficient operation and maintenance of current and future fleets. Specifically, this work establishes a computational framework for evaluating the probabilistic high temperature creep performance of hot-section Ni-based superalloys where uncertainty comes from both microstructural and operational variability. A novel high-fidelity physics model which phenomenologically captures grain-boundary sensitive phenomena has been established. A probabilistic calibration procedure was used to calibrate the model and capture uncertainty in the parameterized model coefficients. A design of experiments methodology was established for identifying informative microstructural digital representations for suitable for forward model evaluation. Results show that training a machine-learning surrogate using this design criteria outperforms random selection of microstructural representations. Finally, two surrogate models were developed: (1) a deterministic surrogate model which predicts the local field response given microstructure, constitutive model parameters, and operating conditions (stress, temperature) and (2) a probabilistic model, where uncertainty comes from constitutive law uncertainty, built using denoising diffusion probabilistic models which samples responses given (1) microstructure and (2) operating conditions. These surrogate models enable partner Siemens Energy to rapidly perform UQ analysis specific to creep deformation across a range of microstructures and operating conditions. The impact is that these ML and physics codes can be used to establish more advanced reliability models for the inspection, servicing, and maintenance of land based gas turbine engines.
Measurements of the Higgs boson production times decay rates and differential cross-sections have recently been performed by the ATLAS experiment in several decay channels using up to 139 fb -1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV recorded at the Large Hadron Collider. This paper presents multiple interpretations of these Higgs boson measurements. Measurements of production-mode cross-sections, simplified template cross-sections and fiducial differential cross-sections in different decay channels are reparameterised in terms of the impact of Standard Model effective field theory operators, and constraints are reported on the corresponding Wilson coefficients. Production and decay rate measurements are interpreted in UV-complete extensions of the Standard Model, namely the two-Higgs-doublet model (2HDM) near the alignment limit and the Minimal Supersymmetric Standard Model (MSSM) for various MSSM benchmark scenarios. The constraints on the 2HDM parameters (cos(β - α), tanβ) and the MSSM parameters (m A , tanβ) are complementary to those obtained from direct searches for additional Higgs bosons.
The strength of materials is influenced by a range of external conditions, such as temperature and deformation rate. Consequently, materials that demonstrate substantial variations in their mechanical behavior due to fluctuations in temperature and strain rate require complex strength models to accurately predict material performance in real-world applications. To predict such complex behavior, a robust and flexible strength model is necessary. In this work, we utilize genetic programming-based symbolic regression (GPSR) to develop data-driven strength models that accurately represent the measured stress–strain responses of tin across a wide range of strain, strain rate and temperature regimes. The GPSR models are constrained by physically-informed conditions, which leads to significant improvement in extrapolation. The best model is integrated into a multi-physics code to perform Taylor impact simulations, validating the model’s accuracy and robustness. In conclusion, the model predictions showed excellent agreement with experimental results, particularly when compared to predictions using traditional strength models.
Understanding the neutron background is essential for determining the neutron yield from nuclear reactions. Here, the neutron backgrounds were determined for heavy-ion collision experiment using the shadow-bar method, where beams of 40,48 Ca at 56, 140 MeV/u impinged on targets of 58,64 Ni and 112,124 Sn. In the analysis presented here, brass shadow bars are placed in front of organic liquid scintillator neutron detectors to determine the energy-dependent neutron background fractions. The measurement of neutron spectra with and without shadow bars is important to determine the neutron background more accurately. The neutron background, along with its sources and systematic uncertainties, are explored with a focus on the impact of background models and their dependence on neutron energy.