SEARCH · Engineering Papers

Results for “energy simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Impact of Increased Latent Generations on Sensitivity Calculations with SCALE

Analyses of cross section sensitivity data from systems with fissile material allow analysts to associate an importance for each material, nuclide, reaction, and neutron energy by simulating real world criticality scenarios. Although criticality safety validation efforts can be guided by the cross-section sensitivity and uncertainty data generated for a particular system, these calculations can often be computationally expensive and sometimes cumbersome without proper guidance. The TSUNAMI suite within the SCALE code package has several methods for generating sensitivity data, including multigroup and continuous energy (CE) capabilities. The release of SCALE 6.3 has three different CE methods for generating cross section sensitivity data: (1) the Iterated Fission Probability (IFP) method with the KENO Monte Carlo transport solver, (2) the IFP method with the Shift Monte Carlo transport solver, and (3) the Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method with the KENO Monte Carlo transport solver. Although the CLUTCH method has additional parameters for generating sensitivity data files relative to the IFP method, all three methods use latent generations, which are the generations between an event (i.e., fission) and the assessment of importance based on the asymptotic population of progeny neutrons. Increasing the number of latent generations in a calculation leads to increased discrimination of the sensitivity coefficients but at the cost of the increased uncertainty associated with those generated values. Analysts must balance the accuracy of the sensitivity calculations and its uncertainty with the associated computational cost involved in generating the values. This paper discusses the impact of adjusting the latent generation parameter for a range of sensitivity values and how these changes compare with the direct perturbation values obtained from a change of ±0.5% Δ k in both benchmark and safety application models. Two benchmarks from the International Handbook of Evaluated Criticality Safety Benchmark Experiments and the MPC-32 dual purpose canister for spent nuclear fuel are used for analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Neutrino Transport with Monte Carlo Method. I. Toward Fully Consistent Implementation of Nucleon Recoils in Core-collapse Supernova Simulations

The small energy exchange via nucleon recoils in neutrino–nucleon scattering is now supposed to be one of the important factors for successful explosions of core-collapse supernovae (CCSNe), as they can change neutrino spectra through accumulation of a large number of scatterings. In deterministic methods employed for neutrino transport in CCSN simulations, we normally cannot afford to deploy a large enough number of energy bins needed to resolve this small energy exchange, and subgrid techniques are employed one way or another. In this paper, we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by ~15% for heavy-lepton neutrinos and much smaller amounts for other types of neutrinos in a typical postbounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the deterministic methods. In order to mimic the latter calculation, we redistribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected, depending on the energy resolution. We also discuss the possible implications of our results for the deterministic methods.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Molecular dynamics study of silicon atomic layer etching by chlorine gas and argon ions

Classical molecular dynamics (MD) is used to simulate atomic layer etching processes of silicon by alternating exposure to chlorine gas and argon ions. In order to validate our model, a rigorous comparison is done with ion beam experiments found in the literature [Park et al., Jpn. J. Appl. Phys. 44, 389 (2005)]. Here it is shown that the etch per cycle (EPC) as a function of argon ion energy from simulations is in quantitative agreement with experimental results if the correct argon ion fluence is used in the simulations. The EPC as a function of ion irradiation time and amount of chlorine exposure also show good agreement with the experiment. The MD simulations also show the formation of an amorphous silicon region with chlorine atoms mixed uniformly throughout following ion bombardment. Finally, the etch products during the ion irradiation step are analyzed and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Generation and validation of comprehensive synthetic weather histories using auto-regressive moving-average models

As energy system design moves to more complex methods of optimization including machine learning there is a significant need for more weather data than is available. One method to solve this is using synthetic data models such as the auto-regressive moving-average (ARMA) model which has been frequently utilized to create such data. This paper looks at extending the ARMA algorithm to generate solar components through the use of clearsky detrending, maintaining vector relationships and by leveraging physical relationships. The method for the creation of entirely synthetic weather data files including key weather variables for energy system analysis is presented. Furthermore, a detailed comparison of energy system simulations utilizing both real and synthetic data is made using NREL’s System Advisor Model. Whilst good agreement is made for the solar variables, and other weather variables, ARMA methods often fail to capture the standard deviation and skew of annual weather distributions. Vector-ARMA is shown to maintain correlations between variables and thus generate data sets that perform similarly in energy system design. Here, it is finally shown that the ARMA method fails to preserve day-today correlations in weather variables and thus over-predicts optimal energy storage by 21% for a residential solar application.

42 ENGINEERING↗

DOE OSTI2021

A Parameterization of Sub-grid Topographical Effects on Solar Radiation in the E3SM Land Model (Version 1.0): Implementation and Evaluation Over the Tibetan Plateau

Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth System Models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to evaluate the performance of ELM. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover and surface temperature especially in the high-elevation and snow-cover regions over the TP. Topography has non-negligible effects on surface energy budget, snow cover, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20 W/m 2 , 10 W/m 2 , and 5 W/m 2 , respectively. The differences in land surface albedo, snow cover fraction, and surface temperature between TOP and PP exceed 0.1, 20%, and 1K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations, and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are more significant with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

58 GEOSCIENCES↗

DOE OSTI2023

HydroChrono: An Open-Source Hydrodynamics Package for Project Chrono

In this paper we present the development and verification of HydroChrono, a hydrodynamics package for the Project Chrono physics engine. This package includes the implementation of hydrodynamics equations, the added mass for multibody systems, the development of I/O functions as well as a Python API, and comparison against standard reference cases and other existing tools. HydroChrono provides a flexible, fully open-source solution for simulating wave energy converters (WECs), floating offshore wind turbines (FOWTs) platforms, and other hydrodynamic systems. Here we show, via comparisons with existing tools for benchmark verification cases, that HydroChrono accurately models hydrodynamic forces - making it a useful tool for the design and optimization of these systems. Additionally, the integration of HydroChrono with Project Chrono offers access to finite element modeling capabilities and high-fidelity modelling - with Chrono's existing coupling to CFD and SPH codes. This provides numerical modelers with a multifidelity simulation framework for designing and validating these systems. The development of HydroChrono provides a new, open-source solution for simulating hydrodynamic systems. Its compatibility with other simulation tools enables a more streamlined and efficient design process, advancing the field and providing new opportunities for innovation in this area.

BEM↗

DOE OSTI2021

ComStock/COVID-19: U.S. Commercial Building Stock Analysis of COVID-19 Mitigation Strategies

ComStock is a DOE platform, with NREL as the lead developer, that models the U.S. commercial building stock using OpenStudio/EnergyPlus to perform physics-based building energy modeling simulations. The ASHRAE Epidemic Taskforce (ASHRAE-ETF) has released several COVID-19 mitigation strategies for commercial buildings. ComStock was used to analyze the impact of four strategies on the existing commercial building stock. These strategies are upgrading to MERV-13 filters, disabling demand control ventilation (DCV), increasing minimum outdoor air to 100%, and HVAC flushing mode operation. Results show the impact of each measure on energy consumption, electric consumption, gas consumption and peak electric.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Reassessing energy deposition for the ITER 5 MA vertical displacement event with an improved DINA model

The beryllium (Be) main chamber wall interaction during a 5 MA/1.8 T upward, unmitigated VDE scenario, first analysed in [J. Coburn et al., Phys. Scr. T171 (2020) 014076] for ITER, has been re-evaluated using the latest energy deposition analysis software. Updates to the DINA disruption model are summarized, including an improved numerical convergence for the 0D power balance, limitations on the safety factor within the plasma core, and the choice to maintain a constant plasma +halo poloidal cross-section. Such updates result in a broad halo region and higher radiated power fractions compared to previous models. The new scenario lasts for ~75 ms and deposits ~29 MJ of energy, with the radial distribution of parallel heat flux q‖(r)resembling an exponential falloff with an effective λ_q=75-198 mm. A maximum halo width w_h of 0.52 m at the outboard midplane is observed. SMITER field line tracing and energy deposition simulations calculate a q_{⊥,max} of ~83 MW/m^2 on the upper first wall panels (FWP). Heat transfer calculations with the MEMOS-U code show that the FWP surface temperature reaches ~1000 K, well below the Be melt threshold. Variations of this 5 MA scenario with Be im-purity densities from 0 to 3∙10^19 m^-3 also remain below the melt threshold despite differences in energy deposition and duration. These results are in contrast to the early study which predicted melt damage to the first wall [J. Coburn et al., Phys. Scr. T171 (2020) 014076], and emphasize the importance of accurate models for the halo width w_h and the heat flux distribution q‖(r)within that halo width. The 2020 halo model in DINA has been compared with halo current experiments on COMPASS, JET, and Alcator C-Mod, and the preliminary results build confidence in the broad halo width predictions. Results for the 5 MA VDE are compared with those for a 15 MA equivalent, generated using the new DINA model. At the higher current, significant melting of the upper FWP is to be expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Optical and electrochemical effects of H 2 and O 2 bubbles at upward-facing Si photoelectrodes

The effects of the size, contact-angle, and coverage of gas bubbles on solar fuels devices were characterized at cm-scale, upward-facing planar and microwire-array Si photoelectrodes in stagnant electrolytes. Experimental measurements were supported by ray-tracing simulations of surface attached gas bubble films. A dilute, redox-active tracer allowed for the quantification of the mass-transport effects of bubble coverage during photoanodic O 2 (g) evolution at upward-facing photoanodes in 1.0 M KOH(aq.). Measurements of the gas coverage at upward-facing p-Si photocathodes in 0.50 M H 2 SO 4 (aq.) allowed for the nucleation rate and contact angle of H 2 (g) bubbles to be evaluated for systems having various surface free energies. Under simulated solar illumination, the rapid departure of small O 2 (g) bubbles produced stable photocurrents at upward-facing oxygen-evolving Si photoanodes and yielded increased mass-transport velocities relative to a stagnant electrolyte, indicating that bubbles can provide a net benefit to the photoelectrochemical performance of an upward-facing photoanode in solar fuels devices.

14 SOLAR ENERGY↗

DOE OSTI2022

Thermal Instability in the Manufacturing of Wind Turbine Blade Spar Caps – Phase 2

As the wind turbine industry continues to develop, the technology surrounding the design, development, and manufacturing of blades requires advancement as well. Wind turbine blades continue to increase in size, increasing the cost and time associated with manufacturing. Simulation driven manufacturing cycle design is required for the industry to maintain manufacturing schedules and optimize the manufacturing process. A recent effort to reduce the cycle time of wind blade spar caps encountered a recurring manufacturing defect during the exotherm and curing of the resin system after filling in the vacuum-assisted resin transfer molding (VARTM) process. Thermal waves developed in sections of the spar cap resulting in unacceptable part tolerances, requiring the parts to be thrown out. The first phase of this IACMI project focused on characterizing the fabric and epoxy system used in manufacturing and simulating the manufacturing process to determine the root cause of the thermal waves. The second phase of this project refined and validated the manufacturing simulations. The simulation workflow was then used to analyze multiple cure cycle options to determine if the thermal waves would appear. The estimated filling time from the filling simulation was shown to match well to experiment, and the temperature history from the curing simulation lined up well with experiment. This indicates that these simulations could be used to evaluate future curing cycles to be used with manufacturing. It is recommended that the distortion simulation work be completed for this project, as it would allow for further manufacturing cycle evaluation through the prediction of stresses and strains in the part. This could provide indicators of defects like the thermal waves seen in manufacturing earlier.

17 WIND ENERGY↗

DOE OSTI2021

Nonlinear dynamics and transport driven by energetic particle instabilities using a gyro-Landau closure model *

Abstract Energetic particle (EP) destabilized Alfvén eigenmode (AE) instabilities are simulated for a DIII-D experimental case with a pulsed neutral beam using a gyro-Landau moments model which introduces EP phase-mixing effects through closure relations. This provides a computationally efficient reduced model which is applied here in the nonlinear regime over timescales that would be difficult to address with more complete models. The long timescale nonlinear evolution and related collective transport losses are examined including the effects of zonal flow/current generation, nonlinear energy cascades, and EP profile flattening. The model predicts frequencies and mode structures that are consistent with experimental observations. These calculations address issues that have not been considered in previous modelling: the EP critical gradient profile evolution in the presence of zonal flows/currents, and the dynamical nature of the saturated state. A strong level of intermittency is present in the predicted instability-driven transport; this is connected to the zonal flow growth and decay cycles and nonlinear energy transfers. Simulation of intermittent AE-enhanced EP transport will be an important issue for the protection of plasma facing components in the next generation of fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Framework to select robust energy retrofit measures for residential communities

Residential building energy retrofits are essential for enhancing environmental sustainability and reducing energy costs. The selection of retrofit measures is influenced by factors such as building systems, occupant behavior, government policy, weather variability, and climate change, all of which can significantly impact energy performance. Compared to retrofitting individual homes, evaluating and selecting optimal retrofit solutions for an entire community is challenging due to diverse residential compositions and variability present. Therefore, engineering robustness is crucial for ensuring consistent energy performance and resilience across different conditions. In this context, robustness refers to the ability of a retrofit measure to maintain its functionality and remain an optimal choice despite external disturbances or changes in inputs and conditions. This study presents a framework for evaluating the robustness of multiple retrofit measures across various building systems, occupant behaviors, and environmental scenarios at the community level. The framework comprises five key steps: scenario model development, integration of the National Residential Efficiency Measures database, energy performance simulation, cost-benefit aggregation, and retrofit solution selection. Each step enhances the framework’s robustness by incorporating the diversity of building characteristics, occupant behaviors, environmental conditions, retrofit options, and evaluation criteria. The framework’s effectiveness is demonstrated through a case study in southern Michigan in the United States, which includes 63 one-story single-family houses, 121 two-story single-family houses, and 8 townhouses. The study identifies furnace retrofits as the most robust solution for the entire community, consistently achieving source energy reductions of 4.7 %–8.0 % and payback period of 10–20 years across various scenarios. These findings are consistent with previous research, indicating the framework’s potential for broader applications in optimizing community-scale residential energy retrofits.

Shu, Lei↗

DOE OSTI2020

MBMS1.0: An Open-Source Code for Modeling and Simulation of Membrane-Based Dehumidification and Energy Recovery

Membrane-based dehumidification is currently being considered as a promising solution for the building application due to its low cost and very limited energy consumption. Developing a simple and efficient open-source code simulation tool is important for boosting the optimization and evaluation of such device in HVAC community. This paper reports a first-order physics based model which accounts for the fundamental heat and mass transfer of humid-air vapor at feed side to flow stream at permeate side. The current model comprises two membrane mass transfer submodels (i.e. microstructure model and performance map model); and it adopts a segment-by-segment methodology for discretizing heat and mass transfer governing equations. The model is capable of simulating both dehumidifiers and energy recovery ventilators with parallel-flow cross-flow, and counter-flow configurations. The model was validated with the measurements at appropriate device. The practices in dehumidification and energy recovery exchangers are also discussed. The model and open-source codes are expected to become a solid fundament for developing a more comprehensive and accurate membrane-based dehumidification in the future.

Gao, Zhiming↗

DOE OSTI2020

A nonparametric approach to modeling dependence among production parameters for energy cost estimation

Dependence existing among production parameters can significantly affect energy costs. For material removal processes, processing times and material amounts to be cut have positive dependence, and the dependence can affect energy costs by varying a peak electricity load. Existing methods to model the dependence of production parameters are based on parametric approaches that cannot accurately represent the dependence due to the rigid nature of parametric models. To address this issue, this study proposes a method to quantify dependence among manufacturing parameters through the application of the empirical copula to the dependence between the milling processing time and the amount of volume to be cut. The proposed method is illustrated by a case study of a manufacturing facility consisting of milling machines; a total of 27 scenarios are simulated for energy cost estimations. The case study clearly shows that the proposed method can capture the dependence between the milling parameters more accurately than a conventional dependence measure from parametric correlation models. Furthermore, the findings from this study would be useful to estimate energy costs in that the proposed method provides a better fit to real production data than conventional parametric approaches can.

42 ENGINEERING↗

DOE OSTI2022

How close are urban scale building simulations to measured data? Examining bias derived from building metadata in urban building energy modeling

Residential and commercial buildings in the United States accounted for 40% of total energy in 2020. Building energy modeling (BEM) is a useful tool that allows individuals, researchers, companies, or utilities to save energy by optimizing buildings through estimation of building technology savings and performance projection of building energy under various environmental conditions. Urban building energy modeling (UBEM) expands the scope beyond individual buildings to the buildings in a neighborhood, city, utility and more. Yet there is a knowledge gap in the literature as to how these models compare to measured data on an individual and aggregated basis. As UBEM data and methods continue to develop, it is important to consider the accuracy, bias, and limitations of the models. Here, nation-scale data and UBEM software suite named Automatic Building Energy Modeling (AutoBEM) was used to model 50,843 buildings in Chattanooga, Tennessee. The uncalibrated simulation results were compared to aggregated 15-minute electricity data for the year 2019 with visualizations highlighting sources of bias in building data and the AutoBEM framework while considering how they relate to other UBEM methods. Estimation of building type and year of constructions are found to be the major sources of bias. Accounting for the amount of conditioned area per building significantly improves the overall fit of the simulated energy use intensity. it was found that inherent variation in building energy use contributes to R 2 values between 0.008 and 0.095 across building types but slope values near 1 for the total number of buildings. This indicates the need for building aggregation for representative building energy modeling with data sources available at an urban scale while illustrating the need for additional individual building data and model improvement beyond the originally produced UBEM models for individual building analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Basin-Size Mapping: Prediction of Metastable Polymorph Synthesizability Across TaC–TaN Alloys

The sizes of the basins of attraction on the potential energy surface are helpful indicators in determining the experimental synthesizability of metastable phases. In principle, these basins can be controlled with changes in thermodynamic conditions such as composition, pressure, and surface energy. Herein, we use random structure sampling to computationally study how alloying smoothly perturbs basin of attraction sizes. The TaC 1-x N x pseudobinary is an ideal test system given the structural and polymorphic contrast of its parent compounds and their technological relevance as epitaxial substrates for Al 1-x Ga x N. While we find limited thermodynamic stability across all computationally observed phases, random structure sampling shows a significant composition region where the rocksalt basin dominates. As such, we predict the potential for the nonequilibrium synthesis of metastable rocksalt TaC 1-x N x alloys as substrates for Al 1-x Ga x N. At higher nitrogen concentrations, other low-energy metastable polymorphs emerge that continue to retain the hexagonal close packing suitable for III-N growth. Confidence in these trends was established through uncertainty quantification of the basin sizes and energy distributions; such analysis utilized the Beta and Dirichlet distributions. In conclusion, we also find (a) polymorph basin sizes can be rationalized in terms of energetic preferences for different coordination environments; and (b) basin sizes universally shrink with increasing nitrogen content, making the system more prone to amorphous growth.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Improving energy efficiency while reducing anthropogenic heat from buildings: how retrofits influence the building stock and urban microclimate in Los Angeles

Anthropogenic heat (AH) from buildings contributes to urban overheating, especially during heat waves, yet building retrofit studies usually evaluate energy savings without assessing impacts on AH. This study quantifies how common building retrofit measures affect both building energy use and AH emissions across the City of Los Angeles. Using a bottom-up urban building energy modeling framework coupled with high-resolution local weather from the Weather Research and Forecasting model with Building Effect Parameterization (WRF-BEP), we evaluate eleven retrofit measures and two multi-measure retrofit packages. HVAC and LED lighting retrofits provide the largest city-wide annual site energy savings, while roof coating is most effective for reducing AH. A package optimized for energy savings reduces summer site energy use by about 32% (2.3 TWh), while a package incorporating AH-focused measures reduces the total AH by over 50% (137 PJ) with minimal difference in energy savings. The AH-aware package produces substantially greater urban cooling, reducing mean near-surface air temperature by up to 0.62 ℃ and peak temperature by up to 3.79 ℃. These results show that retrofit strategies selected only for energy savings may overlook major opportunities for urban heat mitigation. The study provides a framework for integrating AH into building retrofit planning and urban heat resilience policy.

Anthropogenic heat↗