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

Early calendar life and health prediction of silicon batteries via machine learning with uncertainty quantification

Lithium-ion batteries with silicon anodes promise high energy density but are limited by calendar lifetime. Reducing the long iteration time to obtain experimental results requires predicting calendar lifetime early in a cell's life. In this study, we demonstrate that lightweight machine learning models with feature engineering can provide calendar lifetime estimates from early electrochemical signals. After 1 month of electrochemical aging, the best models achieve 10% error in calendar-life prediction and can separate "bad" from "good" lifetime cells with a mean F1 score of 0.857. As battery systems exhibit inherent variability, four methods for uncertainty quantification are compared, and confidence intervals are demonstrated with an uncertainty of +-3.6 months in lifetime prediction. A feature importance analysis indicates that early patterns in voltage decay are the strongest indicators of calendar lifetime. Finally, this modeling approach has high error when generalizing to new electrode chemistries or testing conditions but with appropriately low confidence.

25 ENERGY STORAGE↗
DOE OSTI2024

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Validation of a Global Geospace Model With a Systems Science Approach Based on Canonical Correlation Analysis

A systems science approach based on canonical correlation analysis (CCA) is applied as a new, behavioral way to validate global geospace models. The biggest novelty of the technique is that it validates models at a system level, whereby a side‐by‐side comparison is performed of CCA applied to a 30‐day observational and the corresponding simulation data sets comprising quiet, moderate and active times. The simulation used the Multiscale Atmosphere‐Geospace Environment (MAGE) model. It is shown that (a) CCA must be combined with sensitivity analysis to be effective, (b) the MAGE model generally reproduces the observed behavior (more so for quieter time intervals), quantified by the intercorrelations between different variables and (c) the technique identifies the SuperMAG SML index as a quantity for which refinements of the model are needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Enzalutamide response in a panel of prostate cancer cell lines reveals a role for glucocorticoid receptor in enzalutamide resistant disease

Representative in vitro model systems that accurately model response to therapy and allow the identification of new targets are important for improving our treatment of prostate cancer. Here we describe molecular characterization and drug testing in a panel of 20 prostate cancer cell lines. The cell lines cluster into distinct subsets based on RNA expression, which is largely driven by functional Androgen Receptor (AR) expression. KLK3 , the AR-responsive gene that encodes prostate specific antigen, shows the greatest variability in expression across the cell line panel. Other common prostate cancer associated genes such as TMPRSS2 and ERG show similar expression patterns. Copy number analysis demonstrates that many of the most commonly gained (including regions containing TERC and MYC ) and lost regions (including regions containing TP53 and PTEN) that were identified in patient samples by the TCGA are mirrored in the prostate cancer cell lines. Assessment of response to the anti-androgen enzalutamide shows a distinct separation of responders and non-responders, predominantly related to status of wild-type AR. Surprisingly, several AR-null lines responded to enzalutamide. These AR-null, enzalutamide-responsive cells were characterized by high levels of expression of glucocorticoid receptor (GR) encoded by NR3C1 . Treatment of these cells with the anti-GR agent mifepristone showed that they were more sensitive to this drug than enzalutamide, as were several of the enzalutamide non-responsive lines. This is consistent with several recent reports that suggest that GR expression is an alternative signaling mechanism that can bypass AR blockade. This study reinforces the utility of large cell line panels for the study of cancer and identifies several cell lines that represent ideal models to study AR-null cells that have upregulated GR to sustain growth.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2019

Preoperative Radiomic Approach to Evaluate Tumor-Infiltrating CD8{sup +} T Cells in Hepatocellular Carcinoma Patients Using Contrast-Enhanced Computed Tomography

Background: To help identify potential hepatocellular carcinoma (HCC) candidates for immunotherapies, we aimed to develop and validate a radiomics-based biomarker (Rad score) to predict the infiltration of tumor-infiltrating CD8{sup +} T cells in HCC patients, and to evaluate the correlation of Rad score with tumor immune characteristics. Methods: Overall, 142 HCC patients (n = 100 and n = 42 in the training and validation sets, respectively) were subjected to radiomic feature extraction. Imaging features and immunochemistry data of patients in the training set were subjected to elastic-net regularized regression analysis to predict the level of CD8{sup +} T cell infiltration. Results: A Rad score for CD8{sup +} T-cell infiltration, which contained seven variables, was developed and was validated in the validation set (area under the curve [AUC]: training set 0.751, 95% confidence interval [CI] 0.656–0.846; validation set 0.705, 95% CI 0.547–0.863). The decision curve indicated the clinical usefulness of the Rad score. A higher Rad score correlated with superior overall and disease-free survival outcomes (p = 0.012 and 0.0088, respectively). Using the pathological slides, we found that the Rad score positively correlated with the percentage of tumor-infiltrating lymphocytes (TILs; Spearman rho = 0.51, p < 0.0001). Moreover, the Rad score could also discriminate inflamed tumors from immune-desert and immune-excluded tumors (Kruskal–Wallis, p < 0.0001), and higher Rad scores could be found in patients with positive programmed cell death ligand 1 expression in tumor/immune cells, as well as those with positive programmed cell death protein 1 expression. Conclusion: The newly developed Rad score was a powerful predictor of CD8{sup +} T-cell infiltration, which could be useful in identifying potential HCC patients who can benefit from immunotherapies when validated in large-scale prospective cohorts.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
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↗