SEARCH · Engineering Papers

Results for “heat map”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Spatial mapping of the localized corrosion behavior of a magnesium alloy AZ31B tungsten inert gas weld

Sections of a magnesium alloy, AZ31B, joined with tungsten inert gas (TIG) welding, were examined with scanning electrochemical microscopy (SECM) and scanning Kelvin probe force microscopy (SKPFM) to investigate corrosion mechanisms by correlating observed corrosion behavior with weld-affected microstructural variations. Insight into the changing nature of the galvanic couples between weld zones and at localized microgalvanic sites were investigated using SECM and SKPFM to map both electrochemically active regions and Volta potential differences across the weld-affected zones. The formation of an Al-Zn solidification network in the fusion zone (FZ) at and near the TIG weld epicenter differs from the outer heat-affected zone (HAZ), where intermetallic particles (IMPs) are the notable secondary phase from the magnesium matrix. These microstructures were mapped with SKPFM before and after brief exposure to a salt solution, revealing micro-galvanic couples as the main driving force to corrosion initiation and propagation within each zone. The IMPs and Al-Zn solidification network act as strong cathodes and govern the corrosion processes. The galvanic coupling and evolution of the intrinsic corrosion behavior between the weld zones is explained by monitoring the hydrogen evolution reaction (HER) with SECM over time. Anodically induced cathodic activation is confirmed for this welded material, as micro-galvanic couples between microstructural features are found to transition over time to broad electrochemically active areas within the weld-affected zones, resulting in polarity reversal as time of exposure proceeds.

Anodically induced cathodic activation↗

DOE OSTI2025

A segmented approach to modeling building height: Delineating high-rise and low-rise buildings for enhanced height estimation

Understanding building height is imperative to the overall study of energy efficiency, population distribution, urban morphologies, emergency response, among others. Currently, existing approaches for modeling building height at scale are hindered by two pervasive issues. First, there is no consistent approach to quantify what a high-rise building is at a macro scale, leaving researchers unable to accurately compare results across geographies and domains. Second, high-rise buildings represent a small fraction of the built environment, implying data imbalance challenges that negatively affect current approaches. This is a problem of practical relevance since information on high-rise buildings is important for studies on urban heat islands, population dynamics, and pollution dispersion. Here, we introduce a novel approach to map building height which first identifies two distinct distributions within the built environment, with one being composed of low-rise buildings and one composed of high-rise buildings. We then develop an ensemble scheme where discrete specialist models are trained for each subset of low-rise buildings and high-rise buildings to infer building height from morphology features. For experiments mapping heights of 4.85 million buildings in Japan, we show an increase of 34 % in accuracy within 3m error when compared to the current state-of-the-art when modeling high-rise buildings, which based on KNN experimentation we define as any building > 12m . Our findings show that such an ensemble framework outperforms the current state-of-the-art approaches, which is especially relevant in relation to inferring height for high-rise buildings, a prominent issue of existing approaches for mapping the built environment.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Mechanically induced thermal runaway severity analysis of Li-ion batteries and continuous energy release monitoring

The large-scale deployment of Li-ion batteries in stationary energy storage and electrical vehicle applications demands a strong focus on safety, particularly on the thermal runaway risk and severity evaluation. A standardized single-side mechanical indentation test protocol was developed to induce an internal short-circuit (ISC) and evaluate cells' thermal runaway severity at different state of charge (SOC). The observed hazard severity (OHS in five categories) and evaluated scores in this work have a comprehensive consideration of each cell's capacity, initial voltage, SOC, temperature and voltage change, allowing a better evaluation of the cells' thermal runaway potential. This method was applied to about 200 Li-ion batteries in order to build an extensive thermal runaway database covering various SOCs, capacities and chemistries. In this study, we monitored the transitions of stored electrochemical energy and applied mechanical energy into both thermal energy and acoustic emissions (AE). The surface temperature and mechanical failures were monitored by infrared imaging and AE to capture critical events within battery cells throughout the mechanical indentation tests. Furthermore, the initial temperature maps can predict two types of follow-up events: thermal runaway or gradual heat release via conduction. Analyzing each cell's severity, AEs, and leveraging the evolving database offer insights into predicting occurrences of thermal runaway. The test method, thermal runaway severity evaluation and prediction, and the corresponding database provide battery designers, manufacturers, and end-users a clear overview of Li-ion batteries' thermal runaway potential under mechanical abuse, advancing the safety design of Li-ion batteries.

Acoustic emission↗

DOE OSTI2024

Global scale assessment of urban precipitation anomalies

Urbanization has accelerated dramatically across the world over the past decades. Urban influence on surface temperatures is now being considered as a correction term in climatological datasets. Although prior research has investigated urban influences on precipitation for specific cities or selected thunderstorm cases, a comprehensive examination of urban precipitation anomalies on a global scale remains limited. This research is a global analysis of urban precipitation anomalies for over one thousand cities worldwide. We find that more than 60% of the global cities and their downwind regions are receiving more precipitation than the surrounding rural areas. Moreover, the magnitude of these urban wet islands has nearly doubled in the past 20 y. Urban precipitation anomalies exhibit variations across different continents and climates, with cities in Africa, for example, exhibiting the largest urban annual and extreme precipitation anomalies. Cities are more prone to substantial urban precipitation anomalies under warm and humid climates compared to cold and dry climates. Cities with larger populations, pronounced urban heat island effects, and higher aerosol loads also show noticeable precipitation enhancements. This research maps global urban rainfall hotspots, establishing a foundation for the consideration of urban rainfall corrections in climatology datasets. This advancement holds promise for projecting extreme precipitation and fostering the development of more resilient cities in the future.

Science & Technology - Other Topics↗

DOE OSTI2022

Interplay of strain and phase evolution of laser powder bed fusion Ti–6Al–4V

While additive manufacturing (AM) provides a method of producing geometrically complex and highly detailed structures, the generation of residual strain in AM processes like laser powder bed fusion (L-PBF) can negatively impact performance-enabling properties. In applications such as orthopedic implants, specific performance windows require optimized microstructures in order to obtain desirable properties from multi-phase alloys like Ti–6Al–4V. This research aims to quantify the microscale origins of strain in L-PBF manufactured Ti–6Al–4V by understanding how strain is distributed at the grain and sub-grain scale, the interplay between phase evolution and strain, and examining post-processing strain relief strategies to control these features. Model spinal cage implants were manufactured from Ti–6Al–4V powder via L-PBF and then subjected to strain relieving heat treatment cycles above and below the Ti–6Al–4V β transus as a function of time and cooling rate. Residual strain was then studied via high resolution electron backscatter diffraction (HR-EBSD), and 2D strain maps with sub-micron resolution were generated for each post-processing state. It was found that macroscale thermal strains decreased with heat treatment time, but additional contributions from phase stabilizing residual strains retained primarily in the α' grains as lattice distortive strain remained. Additionally, the retention of β phase significantly changed the strain and dislocation distribution while reducing overall residual strain. In conclusion, these results were validated and reinforced with 3D mesoscopic micromechanical modeling of strain behavior across simulated microstructures, confirming that the local lattice dilation of α’ martensite is a primary contributor of microscale strain generation and retention in L-PBF Ti–6Al–4V.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Solar Thermal Energy Planner (STEP 1): A New Decision Support Tool for Solar Industrial Process Heat Applications

Solar thermal technologies are a promising technology to supply low-cost thermal energy to industrial processes, but there are often significant barriers to entry to industrial owners considering these technologies for their energy demands. To overcome this barrier and convey economic value to customers, NREL and Sandia National Laboratories developed Solar Thermal Energy Planner (STEP 1), a new web-based decision support tool for solar industrial process heat systems. At SolarPACES 2024, the STEP 1 tool was still under development; progress, methodologies, and a preliminary case study was presented. With the STEP 1 tool launch in May 2025, in this work, the initial version of the full public tool will be presented with demonstrations of its capabilities using a few case studies. First, the user's process heat needs such as location, process media (e.g., steam, air), process temperature, land availability, electricity and fuel costs, among other parameters. STEP 1 features a mapping interface that allows users to draw land and roof boundaries. The process media and temperature inform technology selection criteria modules that determine the appropriate solar thermal collection technologies, as well as congruent heat transfer media (e.g., hot water, oil, salt). Once the solar thermal technology selected, its nominal thermal production for the given site is characterized using NREL's System Advisor Model (SAM). Then, a modified version of NREL's REopt optimal sizing and dispatch optimization tool determines cost-optimal sizing. Within minutes, the user receives the results of the technoeconomics analysis, including the size and performance of the cost-optimal solar-plus-storage system. The cost of the system is compared to business-as-usual (e.g., an existing, standalone natural gas boiler). Users can download key results to store for sensitivity analyses. Examples of flat plate collector, parabolic trough, and molten salt tower applications with and without PV hybridization for different industrial facility types are presented in this work. The STEP 1 tool aims to reduce barriers to the adoption of solar heating solutions stemming from a lack of familiarity and technical background with solar system design options and costs among industry stakeholders.

14 SOLAR ENERGY↗

DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Raman spectroscopic characterization of SiO 2 phase transformation and Si substrate stress relevant to EBC performance

To accurately model the long-term durability of environmental barrier coatings (EBCs), a more complete understanding of the phase composition and transformations of the thermally grown oxide SiO 2 (TGO) is desired. For the TGO formed during thermal cycling in steam, cristobalite formation and the subsequent β- to α-cristobalite transformation has been identified as a potentially life-limiting mechanism. In this study, Raman micro-spectroscopy was used to quantify the cristobalite transformation on a polycrystalline Si coupon that was exposed to steam at 1350°C for 100 h. The phase transformation was mapped at 200–260°C on the TGO surface at different ramp rates using a heating stage and a micro-positioning stage. The stress in the Si substrate was also determined using Raman spectroscopy by measuring the stress induced peak shift. The α→β phase transformation produced a 300–500 MPa tensile stress in the Si substrate, which compared well to the stress predicted from the volumetric expansion of the cristobalite. In conclusion, quantifying the phase transformation and residual stress are critical tools in developing the next generation of high performance EBCs.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Develop Accurate Techniques for Passive SiC Temperature Monitoring of Miniature Samples for Cross-Cutting Applications

Passive thermometry is critically important because most fuels and materials irradiation experiments are not instrumented, and it is necessary to understand the irradiation temperature to properly interpret any post-irradiation examination data, including evolving properties and/or microstructures. The standard passive thermometry approach uses continuous dilatometry to evaluate changes in the instantaneous coefficient of thermal expansion during post-irradiation thermal annealing. This approach has limitations in terms of sample size (minimum length requirements) and the maximum irradiation temperature that can be accurately determined, which is limited by the reduced swelling (and therefore recovery) following higher temperature irradiation and limitations on the furnaces used with push-rod dilatometers. This work evaluates two new proposed techniques for post-irradiation evaluation of passive SiC temperature monitors: differential scanning calorimetry (DSC) and Raman spectroscopy. DSC is an extremely sensitive technique that can be used for any specimen geometry and is capable of higher temperature operation. Raman spectroscopy is similar in that it is a surface technique capable of examining extremely small samples (submillimeter), can be used with a heated stage up to 1,500°C (planned for future work), and is capable of mapping local irradiation temperatures throughout a sample. Existing SiC samples that were previously irradiated over a wide range of temperatures were cut into multiple pieces to allow for annealing studies using multiple different techniques: dilatometry, DSC, and Raman spectroscopy. This approach mitigates the concern that samples analyzed using one technique may have a slightly different irradiation history than those analyzed using a different technique. Recovery was clearly observed during annealing using both DSC and dilatometry. In some cases, a direct comparison could not be made due to some of the DSC runs accidentally including material from multiple specimens and issues with using an alternative DSC sample holder for the highest temperature annealing studies. Nevertheless, one trend was clear: the DSC runs resulted in higher irradiation temperatures compared to those of the dilatometry runs. Part of this could be attributed to the higher temperature ramp rates used during the DSC runs, which are often preferred to reduce noise in the measurements. By comparison, dilatometry has previously been shown to produce better data at lower ramp rates. Future work should further investigate the ideal ramp rate for both techniques to produce consistent results. Additional work should evaluate the best holder material to use for DSC runs exceeding 1,000°C to provide reliable data while preventing interactions between SiC and the holder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Develop Accurate Techniques for Passive SiC Temperature Monitoring of Miniature Samples for Cross-Cutting Applications

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Modeling and Simulation of Air-Source CO2 Heat Pump Water Heater

Carbon dioxide (CO2) has been widely used as working fluid for the vapor-compression refrigeration systems in large marine device. Due to the potential energy efficiency and the favorable environmental properties of CO2 as a working fluid, CO2 heat pump water heater (HPWH) systems are regarded a promising technology for centralized domestic hot water (DHW) heating in residential and commercial buildings. However, there is still at the early stage of appropriately optimizing and improving the energy performance of CO2 HPWH. This requires CO2 HPWH simulation tools capable of capturing the accurate impact of the emerging compressor, throttle device, and heat exchanger technology on CO2 heat transfer and energy efficiency. In this study, high efficiency components (compressors, pumps, fans, heat exchangers) were identified and applied to the state-of-art CO2 HPWH designs and analyzed their performance by using numerical simulation. This was done by simulating the performance of CO2 HPWH using ACMODEL design model combined with the component models developed at Oak Ridge National Laboratory (ORNL) for orifice tube, map-based compressor, and tube-in-tube gas cooler. ACMODEL is an equipment design model for CO2-based air conditioners and heat pumps developed by Purdue University to account for the details of each component. The simulated CO2 HPWH performance was then compared with the heat pump water heater using conventional refrigerants.

Gao, Zhiming↗

DOE OSTI2022

Wireless sensing and control of temperature using magnetic fields

A method and system for sensing and controlling temperature with magnetic fields are provided. The method comprises placing a compound in thermal communication with a number of temperature or heat sources and placing a number of magnets in thermal communication with the compound. A number of magnetic sensors are placed in electromagnetic communication with the number of magnets. Changes in the magnetic field of the magnets are detected by the sensors and used to determine the temperature of the compound according to a model that maps magnetic field characteristics to temperature. The amount of cure of the compound can then be estimated from the temperature. The temperature or heat sources are controlled in response to the temperature measurement and the estimated amount of cure of the compound.

Chen Mazumdar, Yi↗

DOE OSTI2022

Replicating HCCI-like autoignition behavior: What gasoline surrogate fidelity is needed?

This work seeks to characterize the fidelity needed in a gasoline surrogate with the intent to replicate the complex autoignition behavior exhibited within advanced combustion engines, and specifically Homogeneous Charge Compression Ignition (HCCI). A low-temperature gasoline combustion (LGTC) engine operating in HCCI mode and a rapid compression machine (RCM) are utilized to experimentally quantify fuel reactivity, through autoignition and preliminary heat release characteristics. Fuels considered include a research grade E10 U.S. gasoline (RD5-87), three multi-component surrogates (PACE-1, PACE-8, PACE-20), and a binary surrogate (PRF88.4). Each fuel was studied at lean/HCCI-like conditions covering a wide range of temperatures and pressures that are representative of naturally aspirated to high boost engine operation. Detailed chemical kinetic modeling is also undertaken using a recently updated gasoline surrogate kinetic model to simulate the RCM experiments and to provide chemical insight into surrogate-to-surrogate differences. The LGTC engine experiments demonstrate nearly identical reactivity between PACE-20 and RD5-87 across studied conditions, while faster phasing is seen for both PACE-1 and PACE-8 due to their stronger intermediate- and low-temperature heat release (ITHR/LTHR) at naturally aspirated and boosted conditions, respectively. The RCM experiments reveal typical low-temperature, negative temperature coefficient (NTC) and intermediate-temperature autoignition behaviors at all pressure conditions for RD5-87, which are qualitatively reproduced by all surrogates. Quantitative discrepancies in both autoignition and preliminary heat release are observed for all surrogates, while their ability to replicate RD5-87 autoignition behavior follows the order of PACE-20 > PACE-1 > PACE-8 > PRF88.4. Excellent mapping is obtained between the LGTC engine and the RCM, where the engine pressure-time trajectories can be characterized by the regimes represented by the RCM autoignition isopleths. The kinetic model performs commendably when simulating both autoignition and preliminary heat release of PACE-20, while typically overpredicting ignition delay times for PACE-1, PACE-8 and PRF88.4 at high-pressure and low-temperature/NTC conditions. Sensitivity and rate of production (ROP) analyses highlight surrogate-to-surrogate differences in the governing chemical kinetics where n-pentane initiates rapid OH branching at a faster rate and an earlier timing for PACE-20 than iso-pentane does for PACE-1 and PACE-8, making it computationally more reactive than the other surrogates. The current study highlights the need to include non-standardized properties, such as the lean/HCCI-like autoignition characteristics, in addition to ASTM properties (e.g., RON, MON) as metrics of fuel reactivity and targets to be matched when formulating high-fidelity surrogates that fully capture gasoline advanced combustion behavior such as HCCI-like autoignition.

42 ENGINEERING↗

DOE OSTI2024

Learning the boundary-to-domain mapping using Lifting Product Fourier Neural Operators for partial differential equations

Neural operators such as the Fourier Neural Operator (FNO) have been shown to provide resolution-independent deep learning models that can learn mappings between function spaces. For example, an initial condition can be mapped to the solution of a partial differential equation (PDE) at a future time-step using a neural operator. Despite the popularity of neural operators, their use to predict solution functions over a domain given only data over the boundary (such as a spatially varying Dirichlet boundary condition) remains unexplored. In this paper, we refer to such problems as boundary-to-domain problems; they have a wide range of applications in areas such as fluid mechanics, solid mechanics, heat transfer etc. We present a novel FNO-based architecture, named Lifting Product FNO (or LP-FNO) which can map arbitrary boundary functions defined on the lower-dimensional boundary to a solution in the entire domain. Specifically, two FNOs defined on the lower-dimensional boundary are lifted into the higher dimensional domain using our proposed lifting product layer. We demonstrate the efficacy and resolution independence of the proposed LP-FNO for the 2D Poisson equation.

Kashi, Aditya↗

DOE OSTI2022

Thermal Storage Integrated into Air-Source Heat Pumps to Leverage Building Electrification: A Systematic Literature Review

Air-source heat pumps (ASHPs) can support a decarbonized economy by replacing combustion appliances in homes and electrifying heating systems in buildings. However, ASHPs have not seen significant adoption primarily due to deteriorated performance under cold conditions - at very low temperatures they require auxiliary resistance heating to meet the heating demand and defrost the evaporator. The additional heat lowers the system efficiency. Thermal energy storage (TES) is a candidate technology to help overcome some of these issues. This paper presents a systematic literature review to map the existing research on the integration of TES into ASHPs. Our review of 59 publications indicates that thermal storage increases the ASHP coefficient of performance by 27% on average, albeit with higher initial cost compared to conventional fossil-fueled heating systems. Phase change materials may be ideal to be integrated with ASHPs due to their high energy density and compact design, but only a few publications address TES sizing and design. First and Second Laws of Thermodynamics are widely used to create metrics to assess ASHP-TES integration, and only recently have cost and environmental impact been explored. This literature review suggests that more comprehensive metrics are needed to evaluate the potential benefits of ASHP-TES systems.

air source heat pump↗

DOE OSTI2025

Dual-season thermal energy storage-integrated heat pump system for load shifting: experiments, modeling and nationwide analysis

Here, this paper proposes and validates a dual-season thermal energy storage-integrated heat pump (TES-HP) system that shifts building electric load in both cooling and heating modes using a single phase-change TES. The design uses the TES as a heat sink for the condenser in summer and a heat source for the evaporator in winter, thereby reducing on-peak electric demand while ensuring compatible with existing air-distribution systems. A 14-kW TES-HP prototype with six operating modes and novel refrigerant charge management strategies was experimentally tested under varied conditions. Subsequently, data-driven polynomial performance curves were developed and validated against experiments, then coupled with a phase change material (PCM) model and a calibrated building model within a rule-based controller. Single-day case studies demonstrated effective on-peak demand reduction while maintaining thermal comfort. During cooling, hourly power was lowered by up to 1.5 kW; in heating, average hourly power decreased from over 3.5 kW to 1.8 kW, compared to an HP-only system, achieving 64.5 % load shifting, and the need for resistance heating was eliminated. Seasonal simulations showed typical on-peak electricity savings of 25–35 % in cooling and 40–65 % in heating, with the largest benefits on the hottest and coldest days. Extended response-surface analysis and nationwide mapping quantified load shifting as functions of ambient and TES temperatures, revealing a seasonal trade-off in TES phase-change temperature. These results demonstrate the TES-HP system as a practical and scalable solution for grid-interactive buildings that reduces on-peak demand and electricity use while maintaining thermal comfort and enhancing grid flexibility and reliability.

Grid-interactive buildings↗

DOE OSTI2023

Spatially Distributed Ramp Reversal Memory in VO 2

Ramp‐reversal memory has recently been discovered in several insulator‐to‐metal transition materials where a non‐volatile resistance change can be set by repeatedly driving the material partway through the transition. This study uses optical microscopy to track the location and internal structure of accumulated memory as a thin film of VO 2 is temperature cycled through multiple training subloops. These measurements reveal that the gain of insulator phase fraction between consecutive subloops occurs primarily through front propagation at the insulator‐metal boundaries. By analyzing transition temperature maps, it is found, surprisingly, that the memory is also stored deep inside both insulating and metallic clusters throughout the entire sample, making the metal‐insulator coexistence landscape more rugged. This non‐volatile memory is reset after heating the sample to higher temperatures, as expected. Diffusion of point defects is proposed to account for the observed memory writing and subsequent erasing over the entire sample surface. By spatially mapping the location and character of non‐volatile memory encoding in VO 2 , this study results enable the targeting of specific local regions in the film where the full insulator‐to‐metal resistivity change can be harnessed in order to maximize the working range of memory elements for conventional and neuromorphic computing applications.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Updated U.S. Low-Temperature Heating and Cooling Demand by County and Sector

This dataset includes U.S. low-temperature heating and cooling demand at the county level in major end-use sectors: residential, commercial, manufacturing, agricultural, and data centers. Census division-level end-use energy consumption, expenditure, and commissioned power database were dis-aggregated to the county level. The county-level database was incorporated with climate zone, numbers of housing units and farms, farm size, and coefficient of performance (COP) for heating and cooling demand analysis. This dataset also includes a paper containing a full explanation of the methodologies used and maps. Residential data were updated from the latest Residential Energy Consumption Survey (RECS) dataset (2015) using 2020 census data. Commercial data were baselined off the latest Commercial Building Energy Consumption Survey (CBECS) dataset (2012). Manufacturing data were baselined off the latest Manufacturing Energy Consumption Survey (MECS) dataset (2021).

15 GEOTHERMAL ENERGY↗