Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Initial Fermi orbital descriptors for FLOSIC calculations: The quick-FOD method

Fermi orbital descriptors (FODs) play a key role in Fermi-Löwdin orbital self-interaction correction (FLOSIC) calculations used to remove self-interaction from approximate density functionals. Optimal FODs are obtained by minimizing the self-interaction-corrected total energy, and, in this process, identifying initial sets of FODs becomes crucial. Here we propose, implement and test a novel method for automatically initializing FODs, quick-FOD, based on the minimization of an empirical energy expression that involves a Coulomb-like FOD-electron density attraction, an FOD-FOD short-range repulsion, and an exchange-like FOD-density repulsion. Quick-FOD successfully reproduces FOD arrangements in qualitatively good agreement with Lewis theory and with full-fledged FLOSIC calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Operation of helium sub-atmospheric multistage cryogenic centrifugal compressor trains: Part 1 – Steady state modeling and speed ratio selection

Helium cryogenic systems which can provide cooling below the normal boiling point of helium (approximately 4.2 K) are often required by superconducting radio-frequency niobium resonators utilized in modern high-energy particle accelerators. Achieving temperatures below 4.2 K generally involves operating a cryogenic vessel with liquid helium under sub-atmospheric conditions, thereby lowering the saturation pressure and corresponding saturation temperature. Over the last several decades, multi-stage cryogenic centrifugal compressor trains (CC’s) have been operated efficiently and reliably within large-scale cryogenic systems to continuously evacuate helium vapor generated by a device within the vessel, maintaining sub-atmospheric conditions in the vessel while pressurizing the return vapor to above atmospheric conditions. Traditionally, these CC systems have been operated using empirically derived control philosophies and insight gathered from previous operational experience. Recent efforts at the Facility for Rare Isotope Beams (FRIB) have been aimed at the development of a theoretical basis to characterize the operation of multi-stage cryogenic centrifugal compressor train and utilizing predictive model results to generate control parameters. The objective of this research was identifying operational points which adequately balance cryogenic system efficiency, stability, and overall ease of operation. Furthermore, this manuscript provides an overview of the predictive model development, characterization of the FRIB cryogenic centrifugal compressors and implementation of the predicted performance results during steady-state system operation.

Compressor train control↗
DOE OSTI2023

Dark matter is the new BBN

Measurements of the primordial element abundances provide us with an important probe of our universe’s early thermal history, allowing us to constrain the expansion rate and composition of our universe as early as ~ 1 s after the Big Bang. Prior to this time, we have essentially no empirical information on which to base any such claims. In this paper, we imagine a future time in which we have not only detected the particles that make up the dark matter, but have measured their mass and annihilation cross section with reasonable precision. In analogy to the light element abundances, the dark matter abundance in this scenario could be used to study and constrain the expansion rate and composition of our universe at the time of dark matter freeze out, which for a standard thermal relic occurs at T f ~ m χ /20, corresponding to t ~ 4 x 10 -10 s x (TeV/m χ ) 2 , many orders of magnitude prior to the onset of Big Bang nucleosynthesis. Here we consider how such measurements could be used to constrain scenarios which feature exotic forms of radiation or matter, a ultralight scalar, or modifications to gravity, each of which have the potential to be much more powerfully probed with dark matter than with the light element abundances.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

High salinity seawater boiling point elevation: Experimental verification

Boiling point elevation (BPE) for high salinity solutions is a crucial parameter that significantly affects the overall thermal efficiency of desalination, water treatment, and many other industrial processing systems, especially when minimum or zero liquid discharge is a key requirement. This paper addresses the discrepancies among the empirical correlations reported in the literature that were commonly used to determine BPE for seawater solutions. An evaluation of the experimental data with low solution concentration reveals that the data cannot be extrapolated to determine the BPE at high concentration (X > 140 ppt or g/kg). Besides, disparities are established within different equations used to determine the BPE for high salinity seawater. As a result of the minimal availability of experimental data of BPE for X > 140 ppt and the need to eliminate the discrepancies, an experimental setup was built, which can produce reliable and repeatable BPE data for any solution concentration. Boiling points were obtained experimentally for seawater solutions with a wide range of salinity (X between 20 and 280 ppt) and brine solution temperature (t b between 20 and 80°C). The experimental results of boiling point elevations for seawater are further compared with the experimental data and analytical results reported in the literature. The newly obtained experimental results show that the minimum and maximum BPE values for the brine solutions tested are 0.175°C at t b = 80°C and X = 20 ppt, and 11.6°C at t b = 80°C and X = 280 ppt, respectively. Furthermore, the measured BPE data for seawater concentration over 80 ppt are fitted and presented using an appropriate fitting model. The new correlation can accurately predict the BPE for high salinity seawater brines and thus should be very helpful in the design of high-energy-efficient evaporation systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Computational and experimental study of different brines in temperature swing solvent extraction desalination with amine solvents

Rapid global urbanization and high-salinity wastewater disposal from industrial activities have exerted significant pressure on water resources. Over the past few years, temperature swing solvent extraction (TSSE) has been identified as a promising technique to desalinate hypersaline brines. Despite its potential, the TSSE desalination literature has been mainly based on empirical insights, and the limited molecular simulation studies have primarily focused on NaCl brines. Herein, we use molecular dynamics (MD) simulations to study the TSSE desalination of four different brines, namely, KCl, KBr, NaCl, and NaBr using diisopropylamine as the solvent. Based on both bulk and interfacial brine-diisopropylamine MD results, here we investigate the qualitative and quantitative performance of the simulations by benchmarking these results against our experimental evaluations of these same systems. MD results provide satisfactory qualitative agreement with the experimental data of water solubilization in the organic phase and amine solubilization in the aqueous phase for the KBr, KCl, and NaBr brines. Also, the molecular mechanism of solvation of ionic species by water molecules over diisopropylamine suggested by the MD simulations is in agreement with our experimental data. However, larger qualitative and quantitative deviations were observed for the NaCl brines, and this is likely due to polarization and charge transfer effects, as quantified by our quantum chemical calculations.

42 ENGINEERING↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2025

Bipolar transport in diamond under photo-excitation: Evidence of free charge scattering by excitons

The fundamental characteristics of the Coulomb e-h interactions in ultra-pure diamond are investigated through time-resolved photoconductivity measurements in homogeneously photoexcited samples with ohmic-like low resistive B + sandwich contacts. Here, we demonstrate that, in the electric field, the bipolar charge mobility exhibits drastic reduction for photo-injected densities of 10 14 ‐10 16 cm ‐3 that cannot be explained only by mutual scattering between free e-h particles. This phenomenon is clarified by adding scattering caused by excitons whose density grows in the photo-exited system. The comparison of drift mobilities at 300 K in diamond obtained in this study and in pure Si and Ge obtained from literature, where exciton scattering is insignificant, is presented. We derive the empirical expression for the density dependences of bipolar mobility for the electric fields up to 1 kV/cm that could be useful for the development of diamond electronic devices.

Bipolar drift↗
DOE OSTI2021

Linking habitat suitability with a longleaf pine-hardwood model: Building a species-predictive fire-land management framework

Active management of fire-dependent ecosystems for specific species leads to complex tradeoffs, which affect conservation outcomes to other species. Therefore a multi-species evaluation of management actions is required. Habitat Suitability Models (HSMs) can help in predicting the likelihood of species occurrence using corresponding environmental variables and empirical relationships that link occurrence with specific environmental conditions. Incorporating multiple species into HSMs and relating them to habitat dynamics is crucial for ecosystems that require active management with prescribed fire. To address this issue, we developed multi-species HSM driven within an existing population model of the longleaf pine-hardwood ecosystem to assess the suitability of an ecosystem given different fire management strategies and environmental conditions. The population model used in this study provides spatial and temporal changes of longleaf pine-hardwood habitat structure in response to fire. These habitat values are used by the HSM to calculate habitat suitability for three threatened and endangered faunal species of this ecosystem, which all thrive with frequent fire, but have unique habitat requirements. Transient habitat conditions are traced to predict longleaf pine ecosystem trajectories under various management strategies, thereby evaluating current land management actions, such as thinning or prescribed fire frequencies. We tested a suite of environmental conditions to emphasize the sensitivity of the species to different fire management actions. The results of our modeling suggest that maximum suitable habitat for all three species can be achieved with fire frequency occurring at approximately once every three years. The modeling results support current management actions and provide a new habitat assessment tool that incorporates ecological factors for multiple species, thus providing for habitat optimization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling the effect of habitat availability and quality on endangered winter-run Chinook salmon ( Oncorhynchus tshawytscha ) production in the Sacramento Valley

Reconciliation between water uses and sustaining aquatic species populations requires an effort to identify and quantify essential habitat characteristics for ecosystem health and valuation methods to predict ecosystem response to restoration actions. This process is particularly challenging for anadromous fish species such as California's Sacramento River winter-run Chinook salmon, due to their limited geographic range and diverse life history habitat requirements. Tools, such as life-cycle models, are needed to manage population dynamics and quantify the composite effects of processes across space and time. Nevertheless, complex institutions can hinder result interpretation and communication, and limit model use in decision-making. This paper focuses on the federally endangered and endemic Sacramento River winter-run Chinook (Oncorhynchus tshawytscha) by developing a Winter-Run Habitat-based Population Model (WRHAP). Additionally, WRHAP is a conceptual, freshwater rearing stage model that includes alternative rearing habitats reported in the literature (e.g., floodplains, off-channel and tributaries), defines rules of habitat use based on instream conditions, and incorporates a juvenile growth module that combines bioenergetics modeling with empirical growth rates. Model outputs reasonably follow observed out-migration patterns and provide a realistic smolt size distribution arriving at the San Francisco Bay. This effort demonstrates the importance of currently “non-critical” habitats (as defined by the Endangered Species Act) for juvenile development, with floodplain habitat contributing to a quarter of out-migrating biomass (despite < 18% availability and <10% of total rearing days), and off-channel growth being one of the most sensitive parameters (explaining ~13% of average juvenile weight variance). The model shows the utility of a simple population model to explore relationships between habitat quality/quantity and juvenile development, and to assist water/environmental management and decision-making processes focused on species recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Test of model of equivalence of tree height growth and transpiration rates in percolation-based phenomenology for root-soil interaction

We recently introduced a phenomenological model to represent root-soil interactions that accounts for variability in tree growth on local scales and within species (Hunt et al., 2020). Here, using publicly available datasets, we examine the developed model to provide additional empirical evidence, which also gives novel insight into two particular tree species, Eucalyptus regnans and Eucalyptus sieberri. In particular, we discovered two novel aspects of the model: a power-law time dependence with a different power-law exponent than usually heretofore applied, and apparently a direct equivalence of results for transpiration and growth rates. Also, although our original expectation was that growth and transpiration were merely proportional, comparison of the entire BAAD data set (Falster et al., 2015) for tree height with our model (Hunt et al., 2020; Fig. 1) already suggested approximate equality of these disparate measurements (volume/area, height). In this present manuscript we find further evidence for equivalence of the two measurements and provide additional theoretical basis for a model assumption regarding constancy of sap flow rates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Development of physics-consistent conditional diffusion model to overcome data scarcity in critical heat flux

Deep generative modeling provides a powerful pathway to overcome data scarcity in energy-related applications where experimental data are often limited. By learning the underlying probability distribution of the training dataset, deep generative models, such as the diffusion model, can generate high-fidelity synthetic samples that statistically resemble the training data. Such synthetic data generation can significantly enrich the size and diversity of the available training data, and more importantly, improve the robustness of downstream machine learning models in predictive tasks. The objective of this paper is to investigate the effectiveness of diffusion models for overcoming data scarcity in nuclear energy applications. By leveraging a public dataset on critical heat flux which covers a wide range of commercial nuclear reactor operational conditions, we developed a diffusion model that can generate an arbitrary amount of synthetic samples. Since a vanilla diffusion model can only generate samples randomly, we also developed a conditional diffusion model capable of generating targeted critical heat flux data under user-specified thermal-hydraulic conditions. The performance of the diffusion model was evaluated based on its ability to capture empirical feature distributions and pair-wise correlations, as well as to maintain physical consistency. The results showed that both the diffusion model and conditional diffusion model can successfully generate realistic and physics-consistent critical heat flux data. Furthermore, uncertainty quantification results demonstrate that the conditional diffusion model is highly effective in augmenting critical heat flux data while maintaining acceptable levels of uncertainty.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

2p x-ray absorption spectroscopy of 3d transition metal systems

This review provides an overview of the different methods and computer codes that are used to interpret 2p x-ray absorption spectra of 3d transition metal ions. We first introduce the basic parameters and give an overview of the methods used. We start with the semi-empirical multiplet codes and compare the different codes that are available. A special chapter is devoted to the user friendly interfaces that have been written on the basis of these codes. Next we discuss the first principle codes based on band structure, including a chapter on Density Functional theory based approaches. We also give an overview of the first-principle multiplet codes that start from a cluster calculation and we discuss the wavefunction based methods, including multi-reference methods. We end the review with a discussion of the link between theory and experiment and discuss the open issues in the spectral analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Core-to-core X-ray emission spectra from Wannier based multiplet ligand field theory

Recent advances using Density Functional Theory (DFT) to augment Multiplet Ligand Field Theory (MLFT) have led to ab-initio calculations of many formerly empirical parameters. This development makes MLFT more predictive instead of interpretive, thus improving its value for studies of highly correlated 3 d , 4 d , and ƒ-electron systems. Synchrotron time is always at a premium, and tools that provide predictive capabilities have clear value when it comes to planning studies. Here, in this work, we develop a DFT + MLFT based approach for core-to-core Kα x-ray emission spectra (XES) and evaluate its performance for a range of transition metal systems. We find good agreement between theory and experiment, as well as the ability to capture key spectral trends related to spin and oxidation state. We also discuss limitations of the model in the context of the remaining free parameters and suggest directions forward.

Ab initio theoretical treatment↗
DOE OSTI2020

Sequence-to-sequence neural networks for short-term electrical load forecasting in commercial office buildings

The U.S. power grid is transforming to become smarter, cleaner, and more effi- cient. This is leading to the addition of significant distributed variable renew- able generation. Due to the variable nature of renewable generation, the short- and long-term supply-demand imbalances are less predictable, and conventional approaches to mitigating the imbalance will not be efficient or cost-effective. To address this challenge, transactive control technologies have been proposed which balance energy generation and consumption with market activity and in- frastructural limitations. Transactive control requires the ability of individual end-use loads to express flexibility as a function of a transactive signal (e.g., price). Empirical gray- and black-box models have been widely used to express flexibility, and although these approaches are generally easy to construct and simple to use, they do not capture the non-linear behavior that some end-use loads represent . Machine learning approaches have been proposed to address this limitation. Although deep learning approaches for forecasting end-use loads have been explored, certain aspects of the application of deep models to load forecasting are not well understood. These aspects include how much training data is required, and how models should be structured and trained. To that end, this work explores how to approach applying deep recurrent neural networks to short-term electrical load forecasting with a case study of four commercial office buildings. We identify data requirements for training accurate models of whole building electricity use conditioned on outdoor temperature, provide insight into model hyperparameter sensitivity, and demonstrate how readily models can be generalized to unseen buildings.

Skomski, Elliott↗
DOE OSTI2021

Evaluation of Phase Change Plaster/Paste Composites for Building Envelopes

Thermal energy storage (TES) materials increase a building’s thermal mass and provide the flexibility to shave and time-shift energy use. These materials are needed to enable grid-interactive efficient buildings, allowing building occupants to lower costs by leveraging free ambient heating/cooling as well as shifting consumption to cheaper electricity periods. Phase change materials (PCMs) that store heat through a solid-liquid/solid-solid transformation are of particular interest for TES because of their ability to store large amounts of energy across a constant temperature or very small temperature range. Furthermore, previous studies on PCM building envelopes evaluated products that may only be incorporated as part of new construction and are not easily applicable to existing, energy-intensive buildings. In this work, we leverage laboratory PCM characterization and whole-building energy simulation to study the energy and cost savings of PCM plaster/paste coatings, which can be applied post construction to both new builds and retrofits of existing buildings (about 40% of residential buildings in the United States fall into the retrofit category). We evaluated the energy impact of PCM plaster/paste composites for a single-family residential building in the United States, including layer design (location and thickness), composite properties (phase change temperature, energy density, and shape of the enthalpy curve), and U.S. climate zone. The enthalpy curve shape of PCM shows that phase change temperature and energy density are not the only two parameters that must be considered for PCM selection. Maximizing the energy density amid the effective TES zone is key to boosting energy savings. We also studied PCM price targets with the aim of achieving reasonable payback periods based on both block and time-of-use electricity rates. Moreover, we developed an empirical equation to provide a rough estimate of the energy savings potential of incorporating PCM without the need to perform detailed simulations. These results will help guide the selection and design of PCM plaster/paste composite elements for residential building envelopes.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗
DOE OSTI2024

Long-Term field testing of the accuracy and HVAC energy savings potential of occupancy presence sensors in A Single-Family home

The energy-saving potential of occupancy-centric smart thermostats has been extensively explored in simulations but lacked field testing for energy savings quantification and sensor performance assessment in real buildings. This paper presents a long-term field study conducted in a single-family home in Texas, U.S. to evaluate the performance of occupancy-centric controls (OCC) of HVAC (heating, ventilation, and air-conditioning) system in terms of energy savings, sensor accuracy, and impact on electric peak demand. The test site was equipped with a commercial off-the-shelf (COTS) smart thermostat and multiple occupancy presence sensors for OCC implementation. Additionally, a sub-metering system was installed to monitor electricity consumption of various end-use equipment, including the HVAC system. A supplementary device was installed to track the ground-truth occupancy for the accuracy evaluation of the occupancy presence sensor. Scenarios of baseline and OCC controls were alternated weekly over the 20-month testing period. The results indicated an effective OCC execution, as evidenced by indoor temperature profiles. During the 2023 cooling season, OCC achieved total energy savings of 1,958 kWh, corresponding to a 17.6% energy savings ratio. Under certain conditions, daily HVAC energy savings reached as high as 17 kWh, with a savings ratio of 35%. Sensor performance showed an overall accuracy of 83.8%, a False Positive Rate (FPR) of 12.8%, and a False Negative Rate (FNR) of 47.4%. A key limitation was the sensor’s inability to detect stationary occupants during sleep, leading to a midnight FNR of nearly 100% and significantly compromising thermal comfort. Additionally, the implementation of OCC resulted in extended periods of high electricity demand on summer afternoons, affecting occupant’s thermal comfort and posing potential challenges to community-level grid operations if OCC were widely adopted. Furthermore, this study addresses a critical research gap by empirically investigating energy-saving potential and occupancy sensor performance in residential buildings. Through a comprehensive field-testing study, the research examines the interrelationship between sensor accuracy, energy savings, and thermal comfort, an area that has received limited attention in the current literature.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Analytical-based simulation approach for an anion exchange membrane fuel cell

An analytical and empirical-based 1-D, non-isothermal, steady-state model for anion exchange membrane fuel cell capable of capturing two-phase phenomena is presented in this study. Coupled multi-physics including mass and charge transport, electrochemical reactions, heat transfer, and two-phase water transport are considered in the model and the simulated results are compared to experimental data. To better represent actual material properties and localized conditions, the model applies multilayer discretization in the gas diffusion electrode to enhance prediction accuracy. The model successfully predicts the baseline performance at 70 °C, 131 kPa abs., 92% RH with pure H 2 /O 2 gas as well as the limiting current at 10% H 2 . The robust simulation approach allows for simplistic and accurate estimation of cell performance without the complications of applying two-phase parameters and expensive computational need for numerical models. In addition, the results from the sensitivity studies of material properties and operating conditions provide valuable insights on water management strategies and optimal component design for advancing anion exchange membrane fuel cell technology.

1-D model↗
DOE OSTI2024

Optimal environmental and economic performance trade-offs for fifth generation district heating and cooling network topologies with waste heat recovery

Network topology greatly influences both the economic and environmental performance of fifth generation district heating and cooling (5GDHC) systems. In this study the optimal trade-offs between the environmental and economic performance of 5GDHC network topologies for a five-building district with waste heat recovery were explored. A life cycle assessment method was used to calculate the total life cycle CO 2 emissions (LCCO2) associated with the installation and operation of various network topologies. Twelve months of empirical data from a data center cooling system were analyzed to assess its suitability for integration into a 5GDHC system. The most suitable method for utilizing this waste heat was selected based on the ambient loop warm pipe setpoint, waste heat temperature, and district energy system configuration. A multi-objective optimization algorithm was used to select the 5GDHC network topology that provided the optimal trade-off between LCCO2 and life cycle cost (LCC). A trade-off parameter was employed to weigh the importance of each objective in the selection process. The results showed waste heat from the data center was suitable for integration into the 5GDHC system due to its availability and consistent temperature profile. When return temperatures of 25 °C or higher were available from the liquid-cooled system, direct pre-heating of the ambient loop warm pipe was found to be the most effective waste heat integration method. The selection of the network topology that provided the optimal trade-off between LCCO2 and LCC (optimal trade-off topology) was highly dependent on factors such as fuel prices, CO 2 prices, electricity CO 2 emissions factors, availability of waste heat, embodied CO 2 emissions associated with network installation and network infrastructure costs. Optimal trade-off topologies produced substantial LCCO2 reductions relative to corresponding LCC increases. LCCO2 reduction to LCC increase ratios from 5.78 to 117.79 were identified with CO 2 offset costs ranging from 4.77 to 60.08 ($/tCO 2 e).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗