Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability analysis”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients in Consequence Management: Propagated Dose Uncertainty in ICRP 66 Human Respiratory Tract Model

Reference inhalation dose models rely on deterministic biokinetics and reference computational phantoms, limiting their applicability to the variability present in population-specific exposures encountered in emergency response scenarios. Here, this study introduces REDCAL, a Python-based computational framework developed to propagate uncertainty in inhalation dose coefficients using the International Commission on Radiological Protection (ICRP) Publication 66 Human Respiratory Tract Model. REDCAL integrates ICRP deposition and clearance models, systemic biokinetics, and governing physics principles, and leverages Sandia National Laboratories’ Dakota toolkit for uncertainty quantification via Latin Hypercube Sampling. REDCAL was validated against DCAL, with biokinetic retention results differing by less than 1% and effective dose coefficients by less than 2% across all tested radionuclides. Stochastic sampling introduced variability in dose coefficients, with geometric standard deviations (GSD) in committed effective dose coefficients (CEDC) ranging from 1.0 to 1.5, based on lognormal distribution fits. Analysis demonstrated that variations in the activity median aerodynamic diameter (AMAD) notably influenced the computed CEDC values. Smaller particles (<1 µm) increased doses by 20–30% due to deeper lung deposition and prolonged retention for alpha emitting radionuclides, such as 241 Am and 239 Pu. Radionuclides with fast clearance, such as 133 I, demonstrated a dose reduction exceeding 50%, as AMAD increased beyond 5 µm due to upper airway deposition and rapid mucociliary clearance. The greatest GSD among the radionuclides reported in this study was for 241 Am. In most cases, the largest GSDs in the CEDC were associated with larger particle sizes, an expected outcome, as ICRP Publication 66 defines GSD in particle size as a function of AMAD, resulting in an extended tail of the lognormal distribution. The findings support improved inhalation dose assessments and enhance consequence management strategies for the U.S. Federal Radiological Monitoring and Assessment Center by quantifying uncertainty in dose coefficients and strengthening decision-making for emergency response scenarios.

Biokinetic Modeling↗
DOE OSTI2025

The varying temperature- and pressure-induced phase transition pathways in hybrid improper ferroelectric Sr 3 Sn 2 O 7

A variable-temperature and pressure single-crystal diffraction study of hybrid improper ferroelectric Sr 3 Sn 2 O 7 is reported. In combination with symmetry analysis, we reveal that the application of pressure and temperature induce distinct phase transition pathways, driven by a differing response of the octahedral rotations to these stimuli. Contrary to what has been previously predicted, we observe the ferroelectric to paraelectric phase transition between 10.17(18) and 12.13(14) GPa, meaning the hybrid improper ferroelectric phase remains stable to significantly higher pressures than expected.

Ferroelectrics↗
DOE OSTI2024

Ni( II )-binding affinity of Cc NikZ-II and its homologs: the role of the HH -prong and variable loop revealed by structural and mutational studies

Extracytoplasmic Ni(II)-binding proteins (NiBPs) are molecular shuttles involved in cellular nickel uptake. Here, we determined the crystal structure of apo CcNikZ-II at 2.38 Å, which revealed a Ni(II)-binding site comprised of the double His (HH-)prong (His511, His512) and a short variable (v-)loop nearby (Thr59-Thr64, TEDKYT). Mutagenesis of the site identified Glu60 and His511 as critical for high affinity Ni(II)-binding. Phylogenetic analysis showed 15 protein clusters with two groups containing the HH-prong. Metal-binding assays with 11 purified NiBPs containing this feature yielded higher Ni(II)-binding affinities. Replacement of the wild type v-loop with those from other NiBPs improved the affinity by up to an order of magnitude. This work provides molecular insights into the determinants for Ni(II) affinity and paves way for NiBP engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Reduced Erosion Augments Soil Carbon Storage Under Cover Crops

ABSTRACT Cover crops, a promising strategy to increase soil organic carbon (SOC) storage in croplands and mitigate climate change, have typically been shown to benefit soil carbon (C) storage from increased plant C inputs. However, input‐driven C benefits may be augmented by the reduction of C outputs induced by cover crops, a process that has been tested by individual studies but has not yet been synthesized. Here we quantified the impact of cover crops on organic C loss via soil erosion (SOC erosion) and revealed the geographical variability at the global scale. We analyzed the field data from 152 paired control and cover crop treatments from 57 published studies worldwide using meta‐analysis and machine learning. The meta‐analysis results showed that cover crops widely reduced SOC erosion by an average of 68% on an annual basis, while they increased SOC stock by 14% (0–15 cm). The absolute SOC erosion reduction ranged from 0 to 18.0 Mg C −1 ha −1 year −1 and showed no correlation with the SOC stock change that varied from −8.07 to 22.6 Mg C −1 ha −1 year −1 at 0–15 cm depth, indicating the latter more likely related to plant C inputs. The magnitude of SOC erosion reduction was dominantly determined by topographic slope. The global map generated by machine learning showed the relative effectiveness of SOC erosion reduction mainly occurred in temperate regions, including central Europe, central‐east China, and Southern South America. Our results highlight that cover crop‐induced erosion reduction can augment SOC stock to provide additive C benefits, especially in sloping and temperate croplands, for mitigating climate change.

Huang, Wenjuan [Department of Ecology, Evolution, ↗
DOE OSTI2023

Dayflow: CONUS Daily Streamflow Reanalysis, Version 2 (DayflowV2)

The DayflowV2 dataset provides multiple meteorologic forcings driven hourly streamflow information for approximately 2.7 million NHDPlusV2 stream reaches in the conterminous US (CONUS). DaymetV4, Stage-IV, and Analysis of Period of Record for Calibration (AORC) forcings and their corresponding hybrids drive a nationally scalable modeling framework integrating the simulated runoff from the Variable Infiltration Capacity (VIC) model with the Routing Application for Parallel computatIon of Discharge (RAPID) routing model. Streamflow with (Assimilated) and without (Naturalized) streamflow assimilation at US Geological Survey (USGS) streamflow monitoring sites are included in DayflowV2. A comprehensive evaluation of streamflow at 7,526 USGS gauges is performed for both streamflow types. The resulting key evaluation metrics are also included in the Dayflow dataset. The reanalysis data are available for variable periods; 36 years (1980-2015) for DaymetV4 (DayflowV1), 18 years (2002-2019) for Stage-IV and its hybrids, and 40 years (1980-2019) for AORC and its hybrids.

13 HYDRO ENERGY↗
DOE OSTI2025

Myco-CORPSE simulations assessing mycorrhizal carbon allocation across U.S. forests and global change scenarios

Plants allocate a substantial portion of their fixed carbon belowground to mycorrhizal fungi in exchange for nutrients and other benefits. However, most current ecosystem models omit mycorrhizal processes, limiting our ability to predict plant–soil carbon dynamics under environmental change. To address this gap, we used a mycorrhiza-explicit soil biogeochemical model, Myco-CORPSE (Mycorrhizal Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment), to simulate tree carbon allocation to arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) fungi in temperate forests.The dataset includes outputs from two sets of model simulations:1. Perturbation experiments: Simulations across gradients of ECM dominance (0–100%), nitrogen deposition, soil temperature, and net primary productivity (NPP) to test how these factors affect mycorrhizal C allocation and nutrient cycling.2. FIA-based simulations: Model applications to over 1,800 U.S. forest sites using site-specific data from the U.S. Forest Inventory and Analysis (FIA) program, including vegetation composition, mycorrhizal type, climate, litter traits, soil properties, and N deposition.Model outputs include simulated mycorrhizal carbon allocation and related biogeochemical variables, such as soil and microbial carbon and nitrogen stocks. Data are provided in CSV format and organized by experiment type (in separate ZIP files). Python scripts for running simulations, plotting, and spatial mapping are also included and organized similarly. No proprietary software is required. These outputs support a peer-reviewed study and were used to generate figures and tables in the associated publication.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Nuclear Energy in Long-Term System Models: A Multi-Model Perspective

Long-term energy system models–including electric sector capacity expansion models–are widely used tools for informing planning, technology assessment, and policy analysis. Recent decarbonization goals and rapid technological change have increased the need to appropriately represent economic characteristics and technical details of energy system resources, including variable renewable energy, energy storage technologies, carbon-capture-equipped capacity, and nuclear energy. Nuclear power represents about 20% of electricity generation and 50% of carbon-free electricity in the United States as of 2021. However, there are many perspectives on the role of existing and new nuclear in the future U.S. energy system, which is reflected in the broad range of potential contributions reported in the literature. This project aims to understand how issues central to nuclear energy are represented in long-term energy models. Building on earlier collaborations that focused on variable renewable energy and energy storage, this project convenes four modeling teams that use national-scale long-term energy system models from the Electric Power Research Institute, the National Renewable Energy Laboratory, the U.S. Energy Information Administration, and the U.S. Environmental Protection Agency to share methods and data, update models, run coordinated scenarios, and identify research needs. Improving tools can provide more insightful analyses and ensure that methods are more transparent.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Wholesale Electricity Analysis via Simulation & Learning Experiments (WEASLE): Platform Development and Pilot Competition

This document reports the development of the Wholesale Electricity Analysis via Simulation and Learning Experiments (WEASLE) platform and the pilot competition that was conducted to test the platform. Due to the increasing reliance on variable renewable energy resources for bulk power, the pilot competition, called the Energy Storage Participation Algorithm Competition (ESPA-Comp), was used to test the effect of various market designs on storage utilization and market efficiency. Basic details of the platform are provided, including an overview of the market clearing engine, the battery dispatch and degradation models, electric grid topology and resource mix, and software architecture. Two market designs were tested: a two-settlement market analogous to typical ISO design today, and a multi-settlement market that allows additional forward-trading periods during the real-time market. Results from the pilot competition show that the storage bidding problem is nontrivial and is well suited for future challenges. We find that: 1) all four teams utilized different approaches to the bidding problem, 2) different methodological approaches led to substantially different offer behaviors, 3) resource profits are clustered by team and methodological approach, 4) simulated offers reduced market surplus by about 0.5%, 5) substantially different prices between two-settlement and multi-settlement markets albeit minimal difference in overall market surplus.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Genetic Subtypes and Natural Resistance Mutations in HCV Genotype 4 Infected Saudi Arabian Patients

This study aimed to characterize the HCV genetic subtypes variability and the presence of natural occurring resistance-associated substitutions (RASs) in Saudi Arabia patients. A total of 17 GT patients were analyzed. Sequence analysis of NS3, NS5A, and NS5B regions was performed by direct sequencing, and phylogenetic analyses were used to determine genetic subtypes, RAS, and polymorphisms. Nine patients were infected by GT 4a, two with GT 4o and three with GT 4d. Two patients were infected with apparent recombinant virus (4a/4o/4a in NS3/NS5A/NS5B), and one patient was infected with a previously unknown, unclassifiable, virus of GT 4. Natural RASs were found in six patients (35%), including three infected by GT 4a, two by GT 4a/GT 4o/GT 4a, and one patient infected by an unknown, unclassifiable, virus of GT 4. In particular, NS3-RAS V170I was demonstrated in three patients, while NS5A-RASs (L28M, L30R, L28M + M31L) were detected in the remaining three patients. All patients were treated with sofosbuvir plus daclatasvir; three patients were lost to follow-up, whereas 14 patients completed the treatment. A sustained virological response (SVR) was obtained in all but one patient carrying NS3-RAS V170I who later relapsed. GT 4a is the most common subtype in this small cohort of Saudi Arabia patients infected with hepatitis C infection. Natural RASs were observed in about one-third of patients, but only one of them showed a treatment failure.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Hourly Electricity Demand Projections for Eight Combined Climate and Socioeconomic Scenarios

This dataset contains 40 years (1980-2019) of simulated historical hourly electricity demand (i.e., loads) and 80 years (2020-2099) of projected hourly loads for 54 Balancing Authorities (BAs) and 48 states plus the District of Columbia. Details about the scenarios and variables included in this dataset are in the readme.pdf file. The two primary models that created the dataset are a version of the Global Change Analysis Model with detailed sectoral resolution over the United States (GCAM-USA) and the Total ELectricity Loads (TELL) model. Links to the model source code and workflow for deriving the dataset are provided in an accompanying meta-repository: https://github.com/IMMM-SFA/burleyson-etal_2023_applied_energy. Projections are for four future climate scenarios that represent combinations of Representative Concentration Pathways (RCPs) 4.5 and 8.5 combined with two levels of climate model sensitivities: rcp45cooler, rcp45hotter, rcp85cooler, and rcp85hotter. The four climate scenarios are crossed with Shared Socioeconomic Pathways (SSPs) 3 and 5 to yield eight different future load projections: rcp45cooler_ssp3, rcp45cooler_ssp5, rcp45hotter_ssp3, rcp45hotter_ssp5, rcp85cooler_ssp3, rcp85cooler_ssp5, rcp85hotter_ssp3, and rcp85hotter_ssp5. The climate scenarios are from the IM3 Thermodynamic Global Warming (TGW) dataset which is linked below in the related metadata. The related metadata also contains links to a repository containing the raw GCAM-USA output files.

Climate Change↗
DOE OSTI2024

Hourly Electricity Demand Projections for Eight Combined Climate and Socioeconomic Scenarios

This dataset contains 40 years (1980-2019) of simulated historical hourly electricity demand (i.e., loads) and 80 years (2020-2099) of projected hourly loads for 54 Balancing Authorities (BAs) and 48 states plus the District of Columbia. Details about the scenarios and variables included in this dataset are in the readme.pdf file. The two primary models that created the dataset are a version of the Global Change Analysis Model with detailed sectoral resolution over the United States (GCAM-USA) and the Total ELectricity Loads (TELL) model. Links to the model source code and workflow for deriving the dataset are provided in an accompanying meta-repository: https://github.com/IMMM-SFA/burleyson-etal_2024_applied_energy. Projections are for four future climate scenarios that represent combinations of Representative Concentration Pathways (RCPs) 4.5 and 8.5 combined with two levels of climate model sensitivities: rcp45cooler, rcp45hotter, rcp85cooler, and rcp85hotter. The four climate scenarios are crossed with Shared Socioeconomic Pathways (SSPs) 3 and 5 to yield eight different future load projections: rcp45cooler_ssp3, rcp45cooler_ssp5, rcp45hotter_ssp3, rcp45hotter_ssp5, rcp85cooler_ssp3, rcp85cooler_ssp5, rcp85hotter_ssp3, and rcp85hotter_ssp5. The climate scenarios are from the IM3 Thermodynamic Global Warming (TGW) dataset which is linked below in the related metadata. The related metadata also contains links to a repository containing the raw GCAM-USA output files.

Climate Change↗
DOE OSTI2020

Incorporating the Impacts of Climate Change on Hydrology in a Performance Assessment Model - 20403

The evidence of climate change is increasingly well-documented and impacts should be incorporated in performance assessment studies. The current climate literature provides both observational evidence and climate model projections of climate trends and/or climate change in the late 20. and early 21. centuries for North America and the northeast United States. Probabilistic modeling is a core requirement for quantifying uncertainty and evaluating its impacts. Not evaluating future climate states in a performance assessment because of the existence of uncertainty is contradictory to good modeling practices - the most uncertain issues and parameters require the most attention in effective probabilistic modeling. Excluding climate change limits development of modeling information that could aid in effective decision making. In this work we develop methods to use the output from hydrologic models and analysis of historical aerial imagery to quantify and implement the impacts of climate change on hydrologic processes at a nuclear waste site in West Valley, New York. Specifically, we used the HELP (Hydraulic Performance of Landfill Performance) model to characterize key hydrologic processes under both current and future climate conditions to assess the impacts of changing climate on hydrology. A suite of previous climatic models were reviewed and synthesized to produce a cohesive representation of the current state of knowledge of the impact of climate change on important model inputs such as precipitation. Output from the HELP simulations was coupled to the GoldSim model that was used to develop the Probabilistic Performance Assessment (PPA) approach through the application of a novel 'nearest neighbor' technique. First, several thousand realizations were generated from the HELP model using a Latin Hypercube experimental design to ensure adequate coverage of the parameter space of explanatory variables used to drive HELP. For example, porosity is a physical parameter that is used as an input to both HELP and the GoldSim PA model. We then ran sensitivity analysis (SA) algorithms on the output of HELP for each of the responses of interest. For each predictor, each time we build an SA model we get a different value for the sensitivity index (SI). From the collection of all the SA models, the average was computed among all of the SIs to represent the predictor within the context of the nearest neighbor approach. That is, we conduct SA on each HELP outcome for each scenario. This gives us parameter sensitivity indices for the outcomes. We average the parameter sensitivity indices across the outcomes to get the average SI for a scenario. For each realization that is generated from the Goldsim PA model, Goldsim generates random values for physical/empirical parameters that HELP uses as well. For each vector of physical/empirical parameters that Goldsim generates, the vector from the 5,000 HELP runs that is most 'similar' to the Goldsim vector is computed using the nearest neighbor approach. In this context 'similar' means minimization of the SA-weighted sum of the absolute differences among the 5,000 values computed for this statistic, where each value corresponds to a different HELP realization. In order to account for the impacts of climate change, this process was repeated using the spatially downscaled future climate projections. For each of the key parameters of interest, it was assumed that a linear change depicted the relationship between the values for the present day and those for 2100. In this way, the climatically-driven changes in key parameters used to inform the GoldSim model are quantified and incorporated into the PA model output for the future. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Geothermal Deep Direct Use for Turbine Inlet Cooling in East Texas

The National Renewable Energy Laboratory (NREL), the Southern Methodist University Geothermal Laboratory (SMU), Eastman Chemical (Longview, TX), and TAS (Houston, TX) evaluated the feasibility of using geothermal heat to improve the performance of a natural-gas power plant in East Texas. The area of interest is the Eastman Chemical plant in Longview, Texas, which is on the northwestern margin of a geologic region known as the Sabine Uplift. The feasibility study focused on determining the potential for accessing a subsurface hot-water geothermal resource within a 10-km radius of the site to provide thermal energy for absorption chillers. Wells within a 20-km radius are included for broader geological comparison to determine the heat flow, temperature-at-depth, field porosity and permeability. The lithologies of most interest are the Lower Cretaceous Trinity Group and Upper Jurassic Cotton Valley Group. The deeper Cotton Valley formations are hotter (averaging 117 to 130°C), yet permeability and porosity are low. The shallower Trinity Group contains more variability in permeability and porosity and lower temperatures averaging about 98 to 117°C. The shallower formations are considered despite the lower temperature because of increased ability to produce larger volumes of water and extract enough heat before reinjection. The complete SMU analysis is available in the National Geothermal Data System (NGDS). Tapping such deep geothermal sources for direct heating (as opposed to power generation) is known as geothermal deep direct use (DDU). Geothermal DDU has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for what are typically small-scale, variable-demand projects. This project examines the feasibility of geothermal energy integration in a natural-gas combined cycle power station in East Texas. The DDU resource is tapped to drive absorption chillers (24/7) for production of chilled water at 5-10°C (41-50°F). This chilled water is stored until needed, which allows for continuous operation with a relatively small-capacity geothermal/absorption chiller system. When conditions are favorable, the chilled water is dispatched to cool the air entering the compressor stage of a gas combustion turbine. This process, known as turbine inlet cooling (TIC), boosts power production during periods of high temperature and high-power demand. Such systems can enhance grid reliability and reduce the cost for peak-demand power. A simulation model of the power plant was developed in IPSEpro software and validated against operational data from the plant. This model allowed the team to estimate the additional power that could be produced by applying TIC under different operating and ambient conditions. Absorption chiller performance was estimated from vendor sources to determine the production rate of chilled water from the geothermal resource. Geothermal drilling and development costs were estimated using NREL's GEOPHIRES 2.0. The expected lower drilling costs in this region led to an estimated cost of geothermal heat of about $4/MMBtu (1.4 cents/kWh t ). The estimated cost for the absorption chillers and TIC hardware were obtained from literature sources and project partners. Hourly data were obtained for weather, natural gas and electricity prices, and plant operating state for 2017, which served as a representative year. NREL estimated the capital cost, operating cost, and additional electricity production and revenue for different combinations of geothermal capacity, chiller capacity, and water storage-tank size. The analysis drove toward smaller geothermal and chiller systems to reduce equipment cost. A relatively low-cost water storage tank accumulated the near-continuous chilled water output for later use when TIC was most valued.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Co-expression networks in Chlamydomonas reveal significant rhythmicity in batch cultures and empower gene function discovery

The unicellular green alga Chlamydomonas reinhardtii is a choice reference system for the study of photosynthesis and chloroplast metabolism, cilium assembly and function, lipid and starch metabolism, and metal homeostasis. Despite decades of research, the functions of thousands of genes remain largely unknown, and new approaches are needed to categorically assign genes to cellular pathways. Growing collections of transcriptome and proteome data now allow a systematic approach based on integrative co-expression analysis. We used a dataset comprising 518 deep transcriptome samples derived from 58 independent experiments to identify potential co-expression relationships between genes. We visualized co-expression potential with the R package corrplot, to easily assess co-expression and anti-correlation between genes. We extracted several hundred high-confidence genes at the intersection of multiple curated lists involved in cilia, cell division, and photosynthesis, illustrating the power of our method. Surprisingly, Chlamydomonas experiments retained a significant rhythmic component across the transcriptome, suggesting an underappreciated variable during sample collection, even in samples collected in constant light. Our results therefore document substantial residual synchronization in batch cultures, contrary to assumptions of asynchrony. We provide step-by-step protocols for the analysis of co-expression across transcriptome data sets from Chlamydomonas and other species to help foster gene function discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Design and performance of a variable gap system for thermal conductivity measurements of high temperature, corrosive, and reactive fluids

We report high-temperature fluids such as molten salts, liquid metals, and gasses are being proposed for many advanced energy systems including thermal energy storage devices, concentrating solar plants, and advanced nuclear reactor designs. However, the chemical behavior and thermophysical properties of many of these fluids have not been well characterized, which hinders the design, modeling, safety analysis, and deployment of these systems. Thermal conductivity is a property that is especially limited by existing measurement capabilities, which are subject to errors caused by convection, material interaction, radiative heat transfer, and instrument degradation. Therefore, there is a lack of standard, systematic measurement techniques for high-temperature, reactive, and corrosive fluids. In this work, the development of a variable gap thermal conductivity measurement system is detailed. The system is designed to measure the thermal conductivity of highly corrosive and reactive fluids, and survive operation between 100 °C and 800 °C. The effects of convection are minimized by limiting the thickness of the specimen to thin sizes (<0.3 mm). Corrections for radiative heat transfer were included in the working equations to consider specimens with varying optical properties. The design, construction, instrumentation, operating principles, and data analysis techniques are discussed in detail. The system was tested up to 500 °C using helium gas and molten KNO 3 -NaNO 3 to verify the measurement technique and determine the sources of error. At 300 and 400 °C KNO 3 -NaNO 3 , results showed maximum relative error of 6% when compared to results in the literature. The helium results were within 13% of those in the literature at 300 and 400 °C. Higher errors were observed at 500 °C for both fluids, and the sources of these errors are discussed.

42 ENGINEERING↗
DOE OSTI2023

Bulk Transfer Coefficients Estimated From Eddy-Covariance Measurements Over Lakes and Reservoirs

The drag coefficient, Stanton number and Dalton number are of particular importance for estimating the surface turbulent fluxes of momentum, heat and water vapor using bulk parameterization. Although these bulk transfer coefficients have been extensively studied over the past several decades in marine and large-lake environments, there are no studies analyzing their variability for smaller lakes. Here, we evaluated these coefficients through directly measured surface fluxes using the eddy-covariance technique over more than 30 lakes and reservoirs of different sizes and depths. Our analysis showed that the transfer coefficients (adjusted to neutral atmospheric stability) were generally within the range reported in previous studies for large lakes and oceans. All transfer coefficients exhibit a substantial increase at low wind speeds (<3 m s -1 ), which was found to be associated with the presence of gusts and capillary waves (except Dalton number). Stanton number was found to be on average a factor of 1.3 higher than Dalton number, likely affecting the Bowen ratio method. At high wind speeds, the transfer coefficients remained relatively constant at values of 1.6·10 -3 , 1.4·10 -3 , 1.0·10 -3 , respectively. We found that the variability of the transfer coefficients among the lakes could be associated with lake surface area. In flux parameterizations at lake surfaces, it is recommended to consider variations in the drag coefficient and Stanton number due to wind gustiness and capillary wave roughness while Dalton number could be considered as constant at all wind speeds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Equatorial Origin of the Observed Tropical Pacific Quasi-Decadal Variability From ENSO Nonlinearity

Quasi-decadal (QD) climate variability is detected in the tropical Pacific based on the recent 70 years of observations. This QD variability is identified in equatorial sea surface temperatures (SSTs), the pattern of which resembles the Central Pacific (CP) El Niño-Southern Oscillation (ENSO) but extends further meridionally to the northeastern subtropical Pacific. Whereas equatorial upper-ocean heat content and SSTs are in quadrature for ENSO, these two quantities are almost in phase on the QD timescale. Further analysis shows that nonlinear dynamical heating, primarily originating from strong El Niño events, tends to lead QD SSTs by a quarter of its dominant period (approximately 30 months) and shapes the dominant QD periodicity in observations. Here, our results suggest that the observed QD variability largely originates from ENSO nonlinearity and thus is fundamentally different from ENSO's oscillatory nature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Techno-Economic Analysis for a Potential Geothermal District Heating System in Tuttle, Oklahoma: Preprint

Geothermal deep direct use (DDU) has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for a large-scale and variable heat demand. The National Renewable Energy Laboratory (NREL) and University of Oklahoma evaluated the feasibility of a geothermal district heating (GDH) and cooling system in two schools and 250 houses by utilizing existing oil and gas (O&G) wells in Tuttle, Oklahoma. Heating and cooling demand in the two schools and a typical single-family house were modeled using EnergyPlus building energy simulation software. The modeling results indicated that annual heating demand in two schools and 250 houses is approximately 2.61 GWhth, and cooling demand in the two schools is approximately 2.65 GWhth. In this scope, the techno-economic analysis (TEA) was conducted using the GEOPHIRES tool combined with the TOUGH2 reservoir simulator. The reservoir performance, including geothermal heat production capacity, was modeled by the reservoir simulator TOUGH2. Then, levelized cost of heat (LCOH) was calculated using GEOPHIRES version 3.0, which includes new features such as hourly heat load optimization and peak performance evaluation. Geothermal reservoir temperature was estimated as 90.5 degrees C at a total depth of 3.3 km by the regional average temperature gradient of 22.8 degrees C/km and validated by cation geothermometer calculations. Four production scenarios with two different well configurations and two different heat load profiles have been developed for well flow rates ranging between 3.1 kg/s and 9.3 kg/s. The LCOH of the district heating and cooling system was calculated between $95 and $210/MWh ($28/MMBtu to $62/MMBtu) for four different production scenarios. Typical natural gas prices for residential customers in Oklahoma have ranged from 9 to 19 $/MMBtu over the past decade, which indicates a challenge for deployment of such a GDH system.

deep direct-use↗