Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “recent availability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 631 records · Page 35

Estimating scalar turbulent fluxes with slow-response sensors in the stable atmospheric boundary layer

Conventional and recently developed approaches for estimating turbulent scalar fluxes under stable atmospheric conditions are evaluated, with a focus on gases for which fast sensors are not readily available. First, the relaxed eddy accumulation (REA) classical approach and a recently proposed mixing length parameterization, labeled A22, are tested against eddy-covariance computations. Using high-frequency measurements collected from two contrasting sites (the frozen tundra near Utqiaġvik, Alaska, and a sparsely vegetated grassland in Wendell, Idaho, during winter), it is shown that the REA and A22 models outperform the conventional Monin–Obukhov similarity theory (MOST) utilized widely to infer fluxes from mean gradients. Second, scenarios where slow trace gas sensors are the only viable option in field measurements are investigated using digital filtering applied to fast-response sensors to simulate their slow-response counterparts. With a filtered scalar signal, the observed filtered eddy-covariance fluxes are referred to here as large-eddy-covariance (LEC) fluxes. A virtual eddy accumulation (VEA) approach, akin to the REA model but not requiring a mechanical apparatus to separate the gas flows, is also formulated and tested. A22 outperforms VEA and LEC in predicting the observed unfiltered (total) eddy-covariance (EC) fluxes; however, VEA can still capture the LEC fluxes well. This finding motivates the introduction of a sensor response time correction into the VEA formulation to offset the effect of sensor filtering on the underestimated net averaged fluxes. The only needed parameter for this correction is the mean velocity at the instrument height, a surrogate of the advective timescale. The VEA approach is very suitable and simple to use with gas sensors of intermediate speed (∼ 0.5 to 1 Hz) and with conventional open- or closed-path setups.

58 GEOSCIENCES↗

First-Principles Cost Analysis of Advanced High-Temperature Nuclear Plants

Due to the vast number of recent nuclear reactor innovations, particularly those pertaining to generation IV reactor types like high-temperature gas cooled reactors (HTGRs) and sodium-cooled fast reactors (SFRs), and the newer deployment strategies envisioned, such as use of small modular reactors (SMRs) or even microreactors, reliable, detailed, and complete costs of these nuclear innovations are needed in wide availability. The types of models that generally achieve these objectives are those incorporating the fundamental nature of the real-world systems they aspire to predict, such as first-principles models. Furthermore, first principles models typically offer predictiveness that is not attained by most other types of individual-models. However, detailed, first-principles cost modeling of nuclear reactors and entire nuclear plants is relatively limited. To address this limitation in the availability of detailed, predictive models based on fundamentals, we recently developed a range of cost models, mostly based on first-principles methodologies, to project full lifecycle costs (LCCs) of nuclear power plants (NPPs) based on multiple parallel SM-HTG-pebble bed reactors (PBRs) and SM-SFRs.

Prosser, Jacob H. [Strategic Analysis, Inc., Arlin↗

First-Principles Cost Analysis of Advanced High-Temperature Nuclear Plants

Due to the vast number of recent nuclear reactor innovations, particularly those pertaining to generation IV reactor types like high-temperature gas cooled reactors (HTGRs) and sodium-cooled fast reactors (SFRs), and the newer deployment strategies envisioned, such as use of small modular reactors (SMRs) or even microreactors, reliable, detailed, and complete costs of these nuclear innovations are needed in wide availability. The types of models that generally achieve these objectives are those incorporating the fundamental nature of the real-world systems they aspire to predict, such as first-principles models. Furthermore, first principles models typically offer predictiveness that is not attained by most other types of individual-models. However, detailed, first-principles cost modeling of nuclear reactors and entire nuclear plants is relatively limited. To address this limitation in the availability of detailed, predictive models based on fundamentals, we recently developed a range of cost models, mostly based on first-principles methodologies, to project full lifecycle costs (LCCs) of nuclear power plants (NPPs) based on multiple parallel SM-HTG-pebble bed reactors (PBRs) and SM-SFRs.

Prosser, Jacob H. [Strategic Analysis, Inc., Arlin↗

Enabling NLDAS-2 Anomaly Analysis Using Giovanni

A newly implemented feature in Giovanni (GES DISC Interactive Online Visualization and Analysis Interface) allows users to explore and visualize anomaly data from the NLDAS-2 Primary Forcing and Noah model data sets. For a given measurement and location, an anomaly describes how conditions for a particular time period compare to normal conditions, based on long-term averages. Analyzing anomalies is important for monitoring droughts, determining weather trends, and studying land surface processes relevant for meteorology, hydrology, and climate. Using Giovanni to analyze anomalies for NLDAS-2 data allows for these studies to be efficiently conducted for the central North American region. Phase 2 of NLDAS (NLDAS-2) currently runs at an 1/8th degree resolution, in near-real time, with data sets extending back to January 1979. NLDAS-2 provides data for soil moisture, precipitation, temperature, and other hydrology measurements. Hourly, monthly, and 30-year (1980-2009) monthly climatology data are available for several land surface models and forcing data sets. The Giovanni anomaly tool calculates monthly anomalies, for a given user-defined variable, as the difference between the NLDAS-2 monthly climatology data and the monthly data. The resulting anomaly describes how a chosen month compares to the 30-year monthly average. The presentation will demonstrate the capabilities and usefulness of Giovanni's anomaly tool, detail the recently added NLDAS-2 variables for which anomalies are available, and show how users can access the data.

drough↗

Thermodynamic re-modelling of the Cu–Nb–Sn system: Integrating the nausite phase

Currently available Cu–Nb–Sn phase diagrams lack the recently discovered nausite phase (Cu,Nb)Sn 2 , which is an important intermediate in the course of thermal processing of superconducting Nb 3 Sn wires. Processing decisively determines the resulting microstructure of Nb 3 Sn and, thus, its superconducting properties. Lack of suitable and complete phase diagrams, however, obstructs rational design of such thermal processing procedures. To close this gap and to obtain valid knowledge of homogeneity and stability range of nausite, various Cu–Nb–Sn samples, which are heat-treated between 300 °C and 500 °C, are investigated. By means of energy-dispersive X-ray spectroscopy (EDX), a temperature-dependent homogeneity range of nausite is observed, which covers average mole fractions of Cu between 0.09 and 0.15. This is correlated with a change in the mean atomic volume and can be seen in the lattice parameters determined by X-ray diffraction (XRD). Additionally performed first-principles calculations on different CuSn 2 and NbSn 2 model structures confirm this trend. Furthermore, the peritectic decomposition of nausite to NbSn 2 and liquid at 586 °C is determined by means of in situ XRD and differential scanning calorimetry (DSC). By using the CALPHAD (CALculation of PHase Diagrams) approach, all these findings are used to extend a previous thermodynamic description of the Cu–Nb–Sn system by including the nausite as an additional phase. Finally, with this noteworthy integration, the updated modelling of the Cu–Nb–Sn system can be used for optimizing the multistage heat-treatment steps during processing superconducting Nb 3 Sn wires.

36 MATERIALS SCIENCE↗

Neutron resonance transmission analysis prototype system for thorium fuel cycle safeguards

Emerging thorium-based reactor designs and fuel cycles present challenges to traditional non-destructive assay techniques used in international safeguards. Specifically, assaying the masses of 233 U and 235 U when they are present together in samples with high gamma ray backgrounds is difficult because of similar passive neutron signatures and relatively weak gamma-ray emissions of 233 U. The Pacific Northwest National Laboratory (PNNL) and the Massachusetts Institute of Technology (MIT) are developing a compact neutron resonance transmission analysis (NRTA) system as one potential solution to these challenges. The NRTA technique provides isotopic information for a sample via neutron time-of-flight (TOF) measurements that exploit a sample’s epithermal neutron resonance cross-sections. A recently developed portable NRTA system uses a commercially available, pulsed deuterium-tritium neutron generator with a ~2 m flight path and a GS20 lithium glass scintillator detector. Finally, this paper describes the prototype NRTA system design, a refined radiation transport model of the system, preliminary measurements with thorium and uranium sources, and demonstration of a quantitative isotopic estimation algorithm.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exploratory study reveals far reaching systemic and cellular effects of verapamil treatment in subjects with type 1 diabetes

Currently, no oral medications are available for type 1 diabetes (T1D). While our recent randomized placebo-controlled T1D trial revealed that oral verapamil had short-term beneficial effects, their duration and underlying mechanisms remained elusive. Now, our global T1D serum proteomics analysis identified chromogranin A (CHGA), a T1D-autoantigen, as the top protein altered by verapamil and as a potential therapeutic marker and revealed that verapamil normalizes serum CHGA levels and reverses T1D-induced elevations in circulating proinflammatory T-follicular-helper cell markers. RNA-sequencing further confirmed that verapamil regulates the thioredoxin system and promotes an anti-oxidative, anti-apoptotic and immunomodulatory gene expression profile in human islets. Moreover, continuous use of oral verapamil delayed T1D progression, promoted endogenous beta-cell function and lowered insulin requirements and serum CHGA levels for at least 2 years and these benefits were lost upon discontinuation. Thus, the current studies provide crucial mechanistic and clinical insight into the beneficial effects of verapamil in T1D.

60 APPLIED LIFE SCIENCES↗

Effect of injected flux and current temporal phasing on self-organization in the HIT-SI3 experiment

The HIT-SI3 device at the University of Washington uses three oscillating inductive helicity injectors to form and sustain spheromak plasma equilibria. By adjusting the temporal phase of the injector waveforms with respect to each other, the toroidal spectrum of the imposed perturbations can be controlled. Using a recently implemented GPU-based control system, the available mode spectra were explored experimentally by scanning the space of relative injector phasing. In this space, significant variation in the toroidal mode spectrum ($n$ = 1, 2, 3) of the perturbations was observed. Additionally, variation in characteristics of driven equilibria was also observed, including a ≈ $30$% range in toroidal current gain ($I$ $Φ$ / $I$ $Inj$ ). Experimental results are compared with both a composite-equilibria and nonlinear dynamic model, including extended MHD simulations using the NIMROD code and composite Taylor state equilibria computed using the PSI-Tet code. In conclusion, qualitative agreement is seen with the nonlinear models, but not with composite-equilibria models, suggesting the use of nonlinear models to better capture observed plasma dynamics and provide predictive use for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantification and assessment of the atmospheric boundary layer height measured during the AWAKEN experiment by a scanning LiDAR

The atmospheric boundary layer (ABL) height plays a key role in many atmospheric processes as one of the dominant flow length scales. However, a systematic quantification of the ABL height over the entire range of scales (i.e., with periods ranging from one minute to one year) is still lacking in literature. In this work, the ABL height is quantified based on high-resolution measurements collected by a scanning pulsed Doppler LiDAR during the recent American WAKE experimeNt (AWAKEN) campaign. The high availability of ABL height estimates (≈2200 collected over one year and each of them based on 10-min averaged statistics) allows to robustly assess five different ABL height models, i.e., one for convective thermal conditions and four for stable conditions. Thermal condition is quantified by a stability parameter spanning three orders of magnitude and probed by near-ground 3D sonic anemometry. The free-atmosphere stability, quantified by the Brunt–Väisälä frequency, is both calculated from simultaneous radiosonde measurements and obtained from the best fit of two of the chosen ABL height models. Good agreement is found between the data and three of the chosen models, quantified by mean absolute errors on the ABL height between 281 and 585 m. Furthermore, the seasonal variability of the convective ABL height model parameters (−15% to +23% with respect to the year baseline) agrees with the variability of buoyancy-generated turbulence caused by the variation in solar radiation throughout the year.

17 WIND ENERGY↗

Uncertainty Quantification for Data-Driven Machine Learning Models in Nuclear Engineering Applications: Where We Are and What Do We Need?

Machine learning (ML) has been leveraged to tackle a diverse range of tasks in almost all branches of nuclear engineering. Many of the successes in ML applications can be attributed to the recent performance breakthroughs in deep learning, the growing availability of computational power, data, and easy-to-use ML libraries. However, these empirical successes have often outpaced our formal understanding of the ML algorithms. An important but under-rated area is uncertainty quantification (UQ) of ML. ML-based models are subject to approximation uncertainty when they are used to make predictions, due to sources including but not limited to, data noise, data coverage, extrapolation, imperfect model architecture and the stochastic training process. The goal of this paper is to clearly explain and illustrate the importance of UQ of ML. We will elucidate the differences in the basic concepts of UQ of physics-based models and data-driven ML models. Various sources of uncertainties in physical modeling and data-driven modeling will be discussed, demonstrated, and compared. We will also present and demonstrate a few techniques to quantify the ML prediction uncertainties, including Monte Carlo dropout, deep ensemble, Bayesian neural networks, Gaussian Processes and conformal prediction. Lastly, we will discuss the need for building a verification, validation and UQ framework to establish ML credibility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

BLAST-Lite (Battery Lifetime Analysis and Simulation Tool - Lite) [SWR-22-69] Related to: BLAST aka: BLAST-Py

Battery Lifetime Analysis and Simulation Toolsuite (BLAST) provides a library of battery lifetime and degradation models for various commercial lithium-ion batteries from recent years. Degradation models are identified from publicly available lab-based aging data using NREL's battery life model identification toolkit. The battery life models predicted the expected lifetime of batteries used in mobile or stationary applications as functions of their temperature and use (state-of-charge, depth-of-discharge, and charge/discharge rates). Model implementation is in both Python and MATLAB programming languages. The MATLAB code also provides example applications (stationary storage and EV), climate data, and simple thermal management options. For more information on battery health diagnostics, prediction, and optimization, see NREL's Battery Lifespan webpage.

Smith, Kandler↗

Near-Real-Time Forecast of Satellite-Based Soil Moisture Using Long Short-Term Memory with an Adaptive Data Integration Kernel

Nowcasts, or near-real-time (NRT) forecasts, of soil moisture based on the Soil Moisture Active and Passive (SMAP) mission could provide substantial value for a range of applications including hazards monitoring and agricultural planning. To provide such a NRT forecast with high fidelity, we enhanced a time series deep learning architecture, long short-term memory (LSTM), with a novel data integration (DI) kernel to assimilate the most recent SMAP observations as soon as they become available. The kernel is adaptive in that it can accommodate irregular observational schedules. Testing over the CONUS, this NRT forecast product showcases predictions with unprecedented accuracy when evaluated against subsequent SMAP retrievals. It showed smaller error than NRT forecasts reported in the literature, especially at longer forecast latency. The comparative advantage was due to LSTM’s structural improvements, as well as its ability to utilize more input variables and more training data. The DI-LSTM was compared to the original LSTM model that runs without data integration, referred to as the projection model here. We found that the DI procedure removed the autocorrelated effects of forcing errors and errors due to processes not represented in the inputs, for example, irrigation and floodplain/lake inundation, as well as mismatches due to unseen forcing conditions. The effects of this purely data-driven DI kernel are discussed for the first time in the geosciences. Furthermore, this work presents an upper-bound estimate for the random component of the SMAP retrieval error.

54 ENVIRONMENTAL SCIENCES↗

Available Drawdowns for Each Oil Storage Cavern in the Strategic Petroleum Reserve (2021 Annual Report)

The Department of Energy maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns owned by the Strategic Petroleum Reserve (SPR). This information is important for assessing the SPRs ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. Sandia was directed to develop and implement a process to continuously assess and report the evolution of drawdown capacity, the subject of this report. A cavern has an available drawdown if after that drawdown, the long-term stability of the cavern, the cavern field, or the oil quality are not compromised. Thus, determining the number of available drawdowns requires the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect of enhanced creep on wellbore integrity, and the sympathetic stress effect of operations on neighboring caverns. A consensus has now been built regarding the assessment of drawdown capabilities and risks for the SPR caverns (Sobolik et al., 2014; Sobolik 2016). The process involves an initial assessment of the pillar - to - diameter (P/D) ratio for each cavern with respect to neighboring caverns. A large pillar thickness between adjacent caverns should be strong enough to withstand the stresses induced by closure of the caverns due to salt creep. The first evaluation of P/D includes a calculation of the evolution of P/D after a number of full cavern drawdowns. The most common storage industry standard is to keep this value greater than 1.0, which should ensure a pillar thick enough to prevent loss of fluids to the surrounding rock mass. However, many of the SPR caverns currently have a P/D less than 1.0 or will likely have a low P/D after one or two full drawdowns. For these caverns, it is important to examine the structural integrity with more detail using geomechanical models. Finite-element geomechanical models have been used to determine the stress states in the pillars following successive drawdowns. By computing the tensile and dilatant stresses in the salt, areas of potential structural instability can be identified that may represent red flags for additional drawdowns. These analyses have found that many caverns will maintain structural integrity even when grown via drawdowns to dimensions resulting in a P/D of less than 1.0. The analyses have also confirmed that certain caverns should only be completely drawn down one time. As the SPR caverns are utilized and partial drawdowns are performed to remove oil from the caverns (e.g., for occasional oil sales , purchases, or exchanges authorized by the Congress or the President), the changes to the cavern caused by these procedures must be tracked and accounted for so that an ongoing assessment of the caverns drawdown capacity may be continued. A proposed methodology for assessing and tracking the available drawdowns for each cavern was presented in Sobolik et al. (2018). This report is the latest in a series of annual reports, and it includes the baseline available drawdowns for each cavern, and the most recent assessment of the evolution of drawdown expenditure for several caverns .

02 PETROLEUM↗

Available Drawdowns for Each Oil Storage Cavern in the Strategic Petroleum Reserve (2022 Annual Report)

The Department of Energy maintains an up-to-date documentation of the number of available full drawdowns of each of the caverns owned by the Strategic Petroleum Reserve (SPR). This information is important for assessing the SPR's ability to deliver oil to domestic oil companies expeditiously if national or world events dictate a rapid sale and deployment of the oil reserves. Sandia was directed to develop and implement a process to continuously assess and report the evolution of drawdown capacity, the subject of this report. A cavern has an available drawdown if after that drawdown, the long-term stability of the cavern, the cavern field, or the oil quality are not compromised. Thus, determining the number of a vailable drawdowns requires the consideration of several factors regarding cavern and wellbore integrity and stability, including stress states caused by cavern geometry and operations, salt damage caused by dilatant and tensile stresses, the effect of enhanced creep on wellbore integrity, and the sympathetic stress effect of operations on neighboring caverns. A consensus has now been built regarding the assessment of drawdown capabilities and risks for the SPR caverns (Sobolik et al., 2014; Sobolik 2016). The process involves an initial assessment of the pillar-to-diameter (P/D) ratio for each cavern with respect to neighboring caverns. A large pillar thickness between adjacent caverns should be strong enough to withstand the stresses induced by closure of the caverns due to salt creep. The first evaluation of P/D includes a calculation of the evolution of P/D after a number of full cavern drawdowns. The most common storage industry standard is to keep this value greater than 1.0, which should ensure a pillar thick enough to prevent loss of fluids to the surrounding rock mass. However, many of the SPR caverns currently have a P/D less than 1.0 or will likely have a low P/D after one or two full drawdowns. For these caverns, it is important to examine the s tructural integrity with more detail using geomechanical models. Finite - element geomechanical models have been used to determine the stress states in the pillars following successive drawdowns. By computing the tensile and dilatant stresses in the salt, areas of potential structural instability can be identified that may represent "red flags" for additional drawdowns. These analyses have found that many caverns will maintain structural integrity even when grown via drawdowns to dimensions resulting in a P/D of less than 1.0. The analyses have also confirmed that certain caverns should only be completely drawn down one time. As the SPR caverns are utilized and partial drawdowns are performed to remove oil from the caverns (e.g., for occasional oil sales , purchases, or exchanges authorized by the Congress or the President), the changes to the cavern caused by these procedures must be tracked and accounted for so that an ongoing assessment of the cavern's drawdown capacity may be continued. A proposed methodology for assessing and tracking the available drawdowns for each cavern was presented in Sobolik et al. (2018). This report is the latest in a series of annual reports, and it includes the baseline available drawdowns for each cavern, and the most recent assessment of the evolution of drawdown expenditure for several caverns.

02 PETROLEUM↗

Aeroelastic Modeling for Distributed Wind Turbines: March 11, 2021 - November 10, 2021

Aeroelastic modeling (AM) is the primary methodology for structural and performance assessment of any wind turbine whereby providing an understanding of the impact of design parameters on its loading and power response before witnessing it in the field. Despite these advantages, the use of AM in the Distributed Wind Technology (DWT) sector is limited, especially within the less established manufacturers. This project represents an in-depth assessment of the status of AM and its role within the Standards for the DWT industry. The study gathered input and feedback from a large number of national and international stakeholders, reviewed technical strengths and weaknesses of the current edition of the design standards, analyzed recent industry workshops' and meetings' minutes, collected publicly available AM templates, and provided an evaluation of the existing AM codes. The study achieved several goals including providing strategies for the load assessment categorization of turbines based on rotor swept area and archetype, and guidance for AM verification and validation (V&V), which includes discussions of measurement requirements and a sample test-plan useful for future V&V campaigns and design standard development. This document summarizes the different tasks conducted in the course of the project and highlights the steps required to improve the AM adoption based on a multifaceted approach that encompasses: 1) augmenting AM software capabilities, 2) publishing AM best-practices and design-basis, 3) creating new model templates, 4) providing guidance for V&V of codes and specific turbine models leveraging field testing best-practice, and 5) addressing weaknesses in the current standards. Many of the future objectives identified in this study could leverage NREL's upcoming testing campaigns of three modern distributed wind turbines. Recommendations within this study will advance the value and the ease-of-use of AM, thereby allowing the industry to better capitalize this underutilized tool resulting in a more efficient design process, an easier path to certification, and overall better and more distributed reliable wind turbine products.

17 WIND ENERGY↗

The Hidden Flexibility of the Natural Gas Network for Electric Power Operations: A Case Study of a Near-Miss Winter Event

The U.S. power sector has become increasingly reliant on gas pipeline networks to deliver fuel to natural gas power plants. In addition to supplying relatively low-cost fuel, gas networks offer generators flexibility in their operations through the ability to deliver fuel when needed by using gas storage facilities or linepack if the gas network is at an operating point below its design capacity. However, disruptions or stress events on the gas network - like those occurring in the Northeast and Texas in recent years - can result in limitations on gas availability to generators at times when generation is in short supply. Here we examine a period of stress that occurred in the winter of 2022 in the Western United States. Using data on the region's natural gas pipeline network and electric generators, we build an integrated gas and electric model that closely replicates the actual dispatch of the period. We then evaluate the implications of removing flexibility employed by the gas network operator, which during that period curtailed scheduled gas deliveries to other parties to increase deliveries to natural gas power plants, which requested more gas than initially forecasted. We find that without the flexibility supplied by the gas network operator, there would have been curtailment of gas generation due to gas offtake constraints, requiring the power system operator to redispatch relying on more expensive generation or to potentially shed load. A sensitivity exploring a wind drought further exacerbates the strain, illustrating the potential challenge of managing gas and grid interactions as systems move to higher shares of variable renewable electricity. Based on this example, we discuss potential coordination strategies between the two system operators to ensure that the power system can successfully utilize and rely on the flexibility offered by natural gas networks.

03 NATURAL GAS↗

Cold Climate Air-Source Heat Pump Demonstration and Analysis: Experimental Study on Cold Temperature Performance

Space heating energy is the largest end use for U.S. residential buildings, accounting for nearly 45% of residential building energy consumption nationwide, and approximately 51% of space heating energy is consumed on-site by combusting natural gas (EIA 2020a). In Colorado, approximately 53% of the statewide residential building energy consumption goes to space heating, 76% of which is provided by natural gas (EIA 2020a). Heat pumps are an efficient, electric alternative space heating technology and have been proven viable for decades. However, historically air-source heat pumps (ASHPs) have been limited to temperate climates because of (1) subpar performance at extremely cold outdoor air temperatures, (2) the need for air conditioning in summer months, and (3) the availability of natural gas in colder climates. Relatively recent advances to cold climate ASHP technology, which typically relies on inverter-driven, variable-speed compressors and variable-speed fans, have significantly improved low-temperature heat pump performance, enabling the technology to potentially save energy for many homes in cold climates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗