Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPECTATION VALUE”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Quantum phase transitions in long-range interacting hyperuniform spin chains in a transverse field

Hyperuniform states of matter are characterized by anomalous suppression of long-wavelength density fluctuations. While most of the interesting cases of disordered hyperuniformity are provided by complex many-body systems such as liquids or amorphous solids, classical spin chains with certain long-range interactions have been shown to demonstrate the same phenomenon. Such spin-chain systems are ideal models for exploring the effects of quantum mechanics on hyperuniformity. It is well-known that the transverse field Ising model shows a quantum phase transition (QPT) at zero temperature. Under the quantum effects of a transverse magnetic field, classical hyperuniform spin chains are expected to lose their hyperuniformity. High-precision simulations of these cases are complicated because of the presence of highly nontrivial long-range interactions. We perform an extensive analysis of these systems using density matrix renormalization group simulations to study the possibilities of phase transitions and the mechanism by which they lose hyperuniformity. Even for a spin chain of length 30, we see discontinuous changes in properties like the “τ order metric” of the ground state, the measure of hyperuniformity, and the second cumulant of the total magnetization along the x-direction, all suggestive of first-order QPTs. An interesting feature of the phase transitions in these disordered hyperuniform spin chains is that, depending on the parameter values, the presence of a transverse magnetic field may lead remarkably to an increase in the order of the ground state as measured by the “τ order metric,” even if hyperuniformity is lost. Therefore, it would be possible to design materials to target specific novel quantum behaviors in the presence of a transverse magnetic field. Our numerical investigations suggest that these spin chains can show no more than two QPTs. We further analyze the long-range interacting spin chains via the Jordan-Wigner mapping onto a system of spinless fermions, showing that under the pairwise-interaction approximation and a mean-field treatment, there can be at most two quantum phase transitions. Here, based on these numerical and theoretical explorations, we conjecture that for spin chains with long-range pair interactions that have convergent cosine transforms, there can be a maximum of two quantum phase transitions at zero temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Cost details from front-end engineering design of piperazine with the advanced stripper

This Front-End Engineering Design (FEED) was funded by the U.S. Department of Energy (DOE) to estimate the cost to capture and compress 90% of the CO 2 from an existing natural gas combined cycle (NGCC) in Denver City, Texas, USA. The FEED used the PZAS (Piperazine with the Advanced Stripper) 2G amine scrubbing technology developed and modeled by The University of Texas at Austin. This FEED is unique in providing more public cost details than other FEEDs funded by DOE. The primary objective of the FEED was to provide a comprehensive estimate for the total installed cost of the capture plant. The estimated capital cost of NGCC at the Mustang Station is $\$727$ million for a capacity of 460 MW and 1.6 million tonnes CO 2 /yr. This includes a contingency of $\$104.6$ million and a contractor’s profit of $\$60.1$ million. The total direct field cost is $\$384.1$ million. With an optimistic fuel value of $\$3$/MMBtu, the estimated cost of capture varies from $\$85$/t at 4% IRR/85% load to $\$170$/t at 10% IRR/52% load. Air cooling is technically feasible but expensive. The air cooling systems for the water wash and pump-around intercooling account for 23.4% of the direct cost. The gas-fired boilers represent only 4% of the direct costs, but steam extraction would reduce energy cost, free up cooling water, and reduce the direct costs of processing additional flue gas and CO 2 from the boilers. The absorbers represent 9.6% of the direct cost in this FEED with no direct contact cooler and only 7.6 m of packing. Furthermore, the solvent cross exchangers are less expensive than expected (2.5% of direct costs). Doubling the number of these exchangers could reduce the heat duty from 3.0 to 2.5 GJ/t CO 2 .

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

Experimental Study on the Transfer of Cs, Sr, Pu and Water in Concrete for the Estimation of Contamination - 20085

As part of the 'Nuclear Energy Science and Technology and Human Resources Development project- through concentrating wisdom-by MEXT' in Japan, the major results from the second year of a project on the analysis of radionuclide contamination mechanisms of concrete and the estimation of contamination distribution at the Fukushima Daiichi nuclear power station active during 2018-2020 are presented. The transfer of Cs, Sr, Pu, and water through concrete was studied experimentally. In a saturated condition, the clay did not affect the transfer and Cs penetrated much faster than Sr. In the dried condition, clay retarded Cs transfer. Carbonation of the concrete retarded the transfer of Cs and Sr. Pu was expected to be precipitated immediately in a high pH condition of concrete. Water suction depends on the temperature. Further, the detailed mechanism was analyzed by {sup 1}H-NMR. (authors)

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

(Ir)Relevance of Disorder for Superconductivity in Cuprates

In cuprates that exhibit high-temperature superconductivity (HTS), as the doping level 𝑝 is increased beyond the optimal value, the superfluid density 𝑛 𝑠⁢0 decreases and eventually vanishes, closely tracking the critical temperature 𝑇 𝑐 . This has been interpreted using a dirty-𝑑-wave extension of the Bardeen-Cooper-Schrieffer theory, assuming that 𝑇 𝑐 and 𝑛 𝑠⁢0 decrease with increasing disorder and pair breaking. Here, to test this hypothesis, we tuned the doping level in overdoped La 2-x ⁢Sr x ⁢CuO 4−𝛿 films by electrolyte gating, measured several parameters used to quantify the level of disorder (the residual resistivity ratio RRR, the mean free path 𝑙 0 , and the electron mobility 𝜇), and studied how 𝑇 𝑐 and 𝑛 𝑠⁢0 scale with these. The experimentally measured dependence turned out to be the opposite of the expected—𝑇 𝑐 and 𝑛 𝑠⁢0 increased with increasing disorder. This brings the question of the true origin of the demise of 𝑇 𝑐 and 𝑛 𝑠⁢0 with overdoping back to the center stage of the quest to decipher the HTS enigma.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Intermediate scattering potential strength in electron-irradiated YBa 2 Cu 3 O 7-δ from London penetration depth measurements

Temperature-dependent London penetration depth, λ(T), of a high quality optimally-doped YBa 2 Cu 3 O 7-δ single crystal was measured using tunnel-diode-resonator technique. Controlled artificial disorder was induced by low-temperature 20 K irradiation by 2.5 MeV electrons at two large doses of 3.8 x 10 19 and 5.3 x 10 19 electrons per cm 2 . The irradiation caused significant suppression of the superconductor's critical temperature, T c , from 94.6 K to 90.0 K, and to 78.7 K, respectively. The low-temperature behavior of λ(T) evolves from a T-linear in pristine state to a T 2 behavior after irradiation, expected for a line-nodal d-wave superconductor. However, the original theory that explained such behavior assumed a unitary limit of the scattering potential, whereas usually in normal metals and semiconductors, Born scattering is sufficient to describe the experiment. To estimate the scattering potential strength, we calculated the superfluid density, ρ s (t = T/T c ) = λ 2 (0)/λ 2 (T), varying the amount and strength of non-magnetic scattering using a self-consistent t-matrix theory. Fitting the obtained curves to a power law, ρ s = 1 - Rt n , and to a polynomial, ρ s = 1 - At - Bt 2 , and comparing the coefficients n in one set and A and B in another with the experimental values, we estimate the phase shift to be around 70 and 65°, respectively. We correlate this result with the evolution of the density of states with nonmagnetic disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A Component-Level Bottom-Up Cost Model for Pumped Storage Hydropower

Pumped Storage Hydropower (PSH) is currently the largest source of utility-scale electricity storage in the U.S. and worldwide. As the accelerating deployment of variable renewable technologies creates opportunity and value for energy storage, it has become increasingly important to characterize PSH costs to understand how it competes. Site-specific considerations and limited cost data in the public domain make it difficult to estimate capital costs for potential new PSH sites. This report documents a spreadsheet-based tool that addresses this challenge and creates a component-level bottom-up cost model for PSH that can be made publicly available for widespread use. It uses detailed site-level physical characteristics and design specifications to calculate key performance and cost parameters for individual components and the project as a whole. The model was developed in consultation with HDR, Inc. and Small Hydro Consulting to ensure it aligns with industry expectations. It enables PSH cost exploration across a wide range of system assumptions and could be customized or extended for the needs of a variety of users.

cost model↗
DOE OSTI2023

Dielectron production at midrapidity at low transverse momentum in peripheral and semi-peripheral Pb–Pb collisions at ${\sqrt{s}}_{\textrm{NN}}$ = 5.02 TeV

The first measurement of the e + e – pair production at low lepton pair transverse momentum (p T,ee ) and low invariant mass (m ee ) in non-central Pb–Pb collisions at ${\sqrt{s}}_{\textrm{NN}}$ = 5.02 TeV at the LHC is presented. The dielectron production is studied with the ALICE detector at midrapidity (|η e | < 0.8) as a function of invariant mass (0.4 ≤ m ee < 2.7 GeV/c 2 ) in the 50–70% and 70–90% centrality classes for p T,ee < 0.1 GeV/c, and as a function of p T,ee in three m ee intervals in the most peripheral Pb–Pb collisions. Below a p T,ee of 0.1 GeV/c, a clear excess of e + e – pairs is found compared to the expectations from known hadronic sources and predictions of thermal radiation from the medium. The m ee excess spectra are reproduced, within uncertainties, by different predictions of the photon–photon production of dielectrons, where the photons originate from the extremely strong electromagnetic fields generated by the highly Lorentz-contracted Pb nuclei. Lowest-order quantum electrodynamic (QED) calculations, as well as a model that takes into account the impact-parameter dependence of the average transverse momentum of the photons, also provide a good description of the p T,ee spectra. The measured $\sqrt{\left\langle {p}_{\textrm{T},\textrm{ee}}^2\right\rangle }$ of the excess p T,ee spectrum in peripheral Pb–Pb collisions is found to be comparable to the values observed previously at RHIC in a similar phase-space region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Accounting for uncertainty in complex alluvial aquifer modeling by Bayesian multi-model approach

Alluvial aquifers by nature are complex caused by varied depositional environments. Developing a reliable groundwater model to represent an alluvial aquifer is non-trivial. Also, relying on a single best calibrated model may not be sufficient because of an inadequate choice of model parameter values. To better understand groundwater dynamics and improve model prediction reliability, this study presents a Bayesian multi-model uncertainty quantification (BMMUQ) framework to account for model parameter uncertainty in complex alluvial groundwater modeling. The methodology was applied to the agriculturally intensive Mississippi River alluvial aquifer (MRAA), Northeast Louisiana. An aquifer architecture was first constructed using 7,259 well logs in the MRAA area which covers three fluvial deposits (alluvium, braided-stream terrace, and braided-stream terrace-loess). A 12-layer MODFLOW model was then developed to address the alluvial aquifer complexity and well calibrated through a genetic algorithm. This study quantified model parameter uncertainty in hydraulic conductivity and specific storage of sand facies. Bayesian model averaging (BMA) with the Expectation Maximization (EM) algorithm was adopted to derive posterior model weights and head variances of 50 alternative conceptual groundwater flow models, and thereby obtains BMA ensemble model predictions instead of only relying on the best calibrated conceptual model. Overall, results show that an estimated around 950 million m3 of groundwater storage loss occurs in 2015 with respect to the beginning of 2004, due to high groundwater demand for irrigation in the MRAA area. Explicitly quantifying model uncertainty can produce more reliable groundwater level predictions from BMA ensemble model. The presented groundwater modeling framework improves our understanding of the MRAA and provides a valuable tool to assist agricultural water management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Anion photoelectron spectroscopy and chemical bonding of ThS 2 – and ThSO –

Anion photoelectron spectra of ThSO – and ThS 2 – were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO – , ThS 2 , and ThS 2 – , as well as the oxygen congeners ThO 2 and ThO 2 – . The ab initio calculations, which included complete basis set extrapolations, spin–orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO 2 – that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO 2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck–Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Altogether the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1–2 kcal/mol.

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

Measurement of the background in the CMS muon detector in ${p}{p}$-collisions at $\sqrt{s} = 13$ TeV

The CMS detector, including its muon system, has been operating at the CERN LHC in increasingly challenging conditions for about 15 years. The muon detector was designed to provide excellent triggering and track reconstruction for muons produced in proton–proton collisons at an instantaneous luminosity ($\mathcal{L}$) of 1 x 10 34 cm –2 s –1 . During the Run 2 data-taking period (2015–2018), the LHC achieved an instantaneous luminosity of twice its design value, resulting in larger background rates and making the efficient detection of muons more difficult. While some backgrounds result from natural radioactivity, cosmic rays, and interactions of the circulating protons with residual gas in the beam pipe, the dominant source of background hits in the muon system arises from proton–proton interactions themselves. Charged hadrons leaving the calorimeters produce energy deposits in the muon chambers. In addition, high-energy particles interacting in the hadron calorimeter and forward shielding elements generate thermal neutrons, which leak out of the calorimeter and shielding structures, filling the CMS cavern. We describe the method used to measure the background rates in the various muon subsystems. These rates, in conjunction with simulations, can be used to estimate the expected backgrounds in the High-Luminosity LHC. This machine will run for at least 10 years starting in 2029 reaching an instantaneous luminosity of $\mathcal{L}$ = 5 x 10 34 cm –2 s –1 and increasing ultimately to $\mathcal{L}$ = 7.5 x 10 34 cm –2 s –1 . These background estimates have been a key ingredient for the planning and design of the muon detector upgrade.

Tytgat, M. (ORCID:0000000239902074)↗
DOE OSTI2020

Changes in Convective Available Potential Energy and Convective Inhibition under Global Warming

Atmospheric convective available potential energy (CAPE) is expected to increase under greenhouse gas–induced global warming, but a recent regional study also suggests enhanced convective inhibition (CIN) over land although its cause is not well understood. In this study, a global climate model is first evaluated by comparing its CAPE and CIN with reanalysis data, and then their future changes and the underlying causes are examined. The climate model reasonably captures the present-day CAPE and CIN patterns seen in the reanalysis, and projects increased CAPE almost everywhere and stronger CIN over most land under global warming. Over land, the cases or times with medium to strong CAPE or CIN would increase while cases with weak CAPE or CIN would decrease, leading to an overall strengthening in their mean values. These projected changes are confirmed by convection-permitting 4-km model simulations over the United States. The CAPE increase results mainly from increased low-level specific humidity, which leads to more latent heating and buoyancy for a lifted parcel above the level of free convection (LFC) and also a higher level of neutral buoyancy. The enhanced CIN over most land results mainly from reduced low-level relative humidity (RH), which leads to a higher lifting condensation level and a higher LFC and thus more negative buoyancy. Over tropical oceans, the near-surface RH increases slightly, leading to slight weakening of CIN. Over the subtropical eastern Pacific and Atlantic Ocean, the impact of reduced low-level atmospheric lapse rates overshadows the effect of increased specific humidity, leading to decreased CAPE.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗
DOE OSTI2020

Thermodynamic signatures of an antiferromagnetic quantum critical point inside a superconducting dome

Recent experiments in unconventional superconductors, and in particular iron-based materials, have reported evidence of an antiferromagnetic quantum critical point (AFM-QCP) emerging inside the superconducting dome of the phase diagram. Fluctuations associated with such an AFM-QCP are expected to promote unusual temperature dependencies of thermodynamic quantities. Here, we compute the T dependence of the specific heat C(T) deep inside a fully gapped s +– superconducting state as the AFM-QCP is approached. We find that at the AFM-QCP, the specific heat C(T) vanishes quadratically with temperature, as opposed to the typical exponential suppression seen in fully gapped BCS superconductors. This robust result is due to a nonanalytic contribution to the free energy arising from the general form of the bosonic (AFM) propagator in the SC state. Away from the AFM-QCP, as temperature is lowered, C(T) shows a crossover from a T2 behavior to an exponential behavior, with the crossover temperature scale set by the value of the superconducting gap and the distance to the QCP. As a result, we argue that these features in the specific heat can be used to unambiguously determine the existence of AFM-QCPs inside the superconducting domes of iron-based and other fully gapped unconventional superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Innovating the next generation of commercial smart building software

Nearly 30% of commercial building energy use is wasted due to equipment faults and HVAC controls problems. The result is increased emissions, compromised comfort and productivity, and less reliable coordination of building power needs with a clean grid. The energy impact alone represents $17 billion in potential savings. Today’s smart building software provides a robust solution to address these operational deficiencies. Energy management and information systems (EMIS) are saving up to 9% on average, with two-year paybacks. They are being incorporated into energy management processes, commissioning services, and utility programs. As effective as they are, two barriers prevent even deeper benefits; limited personnel to fix problems once they are identified, and the expense and time to manually implement changes in control systems. In partnership with the research community, the EMIS industry is developing new capabilities to overcome these barriers. Moving beyond siloed products for either fault detection and diagnostics, or optimal control, these new capabilities empower users to not only automatically identify faults, but also to push corrective action, and control improvements to their buildings. In this paper, several areas for enhancements are documented: ‘one-time’ correction of faults such as setpoints, schedules, and economizer lockouts; short-term active testing for automated proportional integral derivative (PID) loop tuning and functional testing; and continuous supervisory control for demand flexibility and year-round efficiency. Results are presented from a pair of partner implementations out of a dozen providers integrating these enhancements into their products, including field tests from across the country, and insights into operator acceptance and integration into operations and maintenance practices.

Casillas, Armando↗
DOE OSTI2021

Strain hardening magnesium-silicate-hydrate composites (SHMSHC) reinforced with short and randomly oriented polyvinyl alcohol microfibers

Magnesium-silicate-hydrate (M-S-H) cement has emerged as an alternative binder with potentially lower energy requirements and emissions. Due to the intrinsically low pH values of its matrix, however, normal steel reinforcement is not appropriate for M-S-H system. To toughen the matrix and to overcome the brittle nature of the material, a new strain hardening magnesium-silicate-hydrate composite (SHMSHC) is developed for the first time by incorporating 2 vol% short and randomly oriented polyvinyl alcohol (PVA) microfibers. The resulting SHMSHC exhibits significant strain hardening with a tensile strain capacity of more than 3%, a compressive strength beyond 50 MPa and a tensile strength of around 3 MPa. Remarkably, saturated multiple cracking with a tight crack width less than 10 μm invisible to the human naked eyes was observed. The newly developed SHMSHC addresses the fundamental challenge of M-S-H binder and is expected to widen possible application areas of M-S-H system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Radiation-induced bowing of SiC/SiC composites under neutron flux gradients—integral experimental data for model validation

Here, the radiation-induced swelling of SiC and its composites, including strong dependencies on temperature and dose, can drive significant lateral bowing in the presence of temperature and/or dose gradients. In recent years, simulations have been performed to assess the extent of bowing in SiC composite light-water reactor (LWR) fuel cladding and boiling water reactor (BWR) channel boxes. However, to date, no integral experimental data exist to validate these models. This work provides the first experimental bowing evaluation of three ∼380 mm long SiC composite specimens irradiated under varying neutron dose gradients (∼50°C–60°C, 0.03–0.06 dpa): two tubes (∼9.8 mm diameter) and a miniature BWR channel box (∼30 mm square). The measured radiation-induced length swelling (∼0.3%–0.7% linear) was consistently 10%–21% higher than values obtained from 3D finite element structural analyses with inputs from 3D radiation transport calculations. This discrepancy could be at least partially explained by differences in dose rate (∼10 -8 dpa/s) compared to the literature data (∼10-6 dpa/s) used to establish the dose-to-swelling correlations in the model. Nevertheless, the modeled bowing magnitudes (<2 mm) obtained from finite element analyses and simple analytical equations were within the bounds of the experimental measurements for all specimens. With improved confidence in the ability to predict the structural response and measure the macroscopic deformations, future experiments will target transient bowing under neutron flux gradients at representative LWR temperatures and assess whether grid spacers can mitigate the tens of millimeters of bowing that would otherwise be expected in ∼4 m long LWR components.

bowing↗
DOE OSTI2020

Seasonal and CO 2 -Induced Shifts of the ITCZ: Testing Energetic Controls in Idealized Simulations with Comprehensive Models

The Tropical Rain Belts with an Annual Cycle and Continent Model Intercomparison Project (TRACMIP) ensemble—a multimodel ensemble of slab-ocean simulations in idealized configurations—provides a test of the relationship between the zonal mean ITCZ and the cross-equatorial atmospheric energy transports (AHT eq ). In a gross sense, the ITCZ position is linearly related to AHT eq , as expected from the energetic framework. Yet, in many aspects, the TRACMIP model simulations do not conform to the framework. Throughout the annual cycle there are large excursions in the ITCZ position unrelated to changes in the AHT eq and, conversely, substantial variations in the magnitude of the AHT eq while the ITCZ is stationary at its northernmost position. Variations both in the net vertical energy input at the ITCZ location and in the vertical profile of ascent play a role in setting the model behavior apart from the conceptual framework. Nevertheless, a linear fit to the ITCZ–AHT eq relationship captures a substantial fraction of the seasonal variations in these quantities as well as the intermodel or across-climate variations in their annual mean values. The slope of the ITCZ–AHT eq linear fit for annual mean changes across simulations with different forcings and configurations varies in magnitude and even sign from model to model and we identify variations in the vertical profile of ascent as a key factor. A simple sea surface temperature–based index avoids the complication of changes in the vertical structure of the atmospheric circulation and provides a more reliable diagnostic for the ITCZ position.

54 ENVIRONMENTAL SCIENCES↗
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↗