Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak estimation”

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

Laser absorption measurements of electron density in nanosecond-scale atmospheric pressure pulsed plasmas

We report on time-resolved measurements of electron number density by continuous-wave laser absorption in a low-energy nanosecond-scale laser-produced spark in atmospheric pressure air. Laser absorption is a result of free-free and bound-free electron excitation, with the absorption coefficient modeled and evaluated using estimates of the time-variation in electron temperature and probe laser absorption path length. Plasma electron number densities are determined to be as high as n e = 7 x 10 19 cm -3 , and decay to 1/e of their peak values over a period of about 50 ns following plasma formation using a 20 mJ, 10 ns pulse width frequency-doubled Nd:YAG laser. The measured plasma densities at later times are shown to be in reasonable agreement with Stark broadening measurements of the 3s[ 5 S o ]-3p[ 5 P] electronic transition in atomic oxygen at 777 nm. This study provides support for the use of such continuous wave laser absorption for time resolved electron density measurements in low energy spark discharges in air, provided that an estimate of the electron temperature and laser path length can be made by accompanying diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Computationally efficient CFD prediction of bubbly flow using physics-guided deep learning

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale framework was proposed in this paper by applying a physics-guided data-driven approach. Instrumental to this framework, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. In this work, by defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1979

Uranus - Disk structure within the 7300-A methane band

Orthogonal narrow-band (100 A) photoelectric slit scan photometry of Uranus has been used to infer the basic two-dimensional structure of the disk within the 7300-A methane band. Numerical image reconstruction and restoration techniques have been applied to quantitatively estimate the degrees of polar and limb brightening on the planet. Through partial removal of atmospheric smearing, an effective spatial resolution of approximately 0.9 arcsec has been achieved. Peak polar, limb, and central intensities on the disk are in the respective proportions 3:2:1. In addition, the bright polar feature is displaced from the geometric pole towards the equator of the planet.

Price, M. J.↗
NASA NTRS1985

Nonthermal Radio Emission and the HR Diagram

Perhaps the most reliable indicator of non-radiative heating/momentum in a stellar atmosphere is the presence of nonthermal radio emission. To date, 77 normal stellar objects have been detected and identified as nonthermal sources. These stellar objects are tabulated herein. It is apparent that non-thermal radio emission is not ubiquitous across the HR diagram. This is clearly the case for the single stars; it is not as clear for the binaries unless the radio emission is associated with their late-type components. Choosing to make this association, the single stars and the late-type components are plotted together. The following picture emerges: (1) there are four locations on the HR diagram where non-thermal radio stars are found; (2) the peak incoherent 5 GHz luminosities show a suprisingly small range for stars within each class; (3) the fraction of stellar energy that escapes as radio emission can be estimated by comparing the integrated maximum radio luminosity to the bolometric luminosity; (4) there are no apparent differences in L sub R between binaries with two cool components, binaries with one hot and one cool component, and single stars for classes C and D; and (5) The late-type stars (classes B, C, and D) are located in parts of the HR diagram where there is reason to suspect that the surfaces of the stars are being braked with respect to their interiors.

Gibson, D. M.↗
NASA NTRS1985

ROWS estimates of wave height, wind speed, and directional wave spectra for SIR-B underflights off Chile, October 11 and 12, 1984

An account is given of the successful operation of the Goddard Ku-band Radar Ocean Wave Spectrometer (ROWS) on four out of five of the NASA Wallops' P-3 underflights of the Shuttle Imaging Radar-B (SIR-B) off the coast of southern Chile during the period of October 7-12, 1984. ROWS provided surface wind speed, wave height, and directional wave height spectrum estimates from the instrument's nadir altimeter mode and off-nadir spectrometer mode. ROWS estimates of the directional spectra of 300-400-m wavelength swell-dominated 3-5 m seas on October 11 and 12 are presented and compared to SCR (Surface Contour Radar), SIR-B, and Navy GSOWM spectra. Good agreement is found with SCR height spectra for the 12th (when the SCR was operable). Also, surprisingly good agreement is found between ROWS slope spectra and SIR-B image spectra for both the 11th and 12th for both swell and local wind-sea components. ROWS altimeter mode estimates of maximum winds of about 10 m/s on the 11th and 12th are seen to be consistent with the ROWS spectrometer mode observed local wind-sea peak frequencies of about 0.12-0.14 Hz.

Jackson, F. C.↗
NASA NTRS1991

HF radar observations of Pc 5 field line resonances in the midnight/early morning MLT sector

The technique of F region sounding by HF coherent radar is applied to the study of field line resonances equatorward of a region of shear flow in the early morning sector. The motions were predominantly in the geomagnetic east-west direction, indicating north-south electric fields. These oscillations had pronounced peaks in their latitudinal power distribution. For the pulsation at 1.95 mHz, a latitudinal phase shift of 180 deg was observed across the peaks in all the look directions of the radar, and a longitudinal wavelength corresponding to an m value of 3 was obtained. For the 2.6-mHz pulsation, the phase shifts across the peaks had a variation with a look direction that indicated a significant longitudinal as well as latitudinal variation; for this activity, an m value of about 16 is estimated. These features are interpreted in terms of the field line resonance theory, and the possible sources of the pulsation energy are discussed.

Ruohoniemi, J. M.↗
NASA NTRS2011

Optimizing Muscle Parameters in Musculoskeletal Modeling Using Monte Carlo Simulations

Astronauts assigned to long-duration missions experience bone and muscle atrophy in the lower limbs. The use of musculoskeletal simulation software has become a useful tool for modeling joint and muscle forces during human activity in reduced gravity as access to direct experimentation is limited. Knowledge of muscle and joint loads can better inform the design of exercise protocols and exercise countermeasure equipment. In this study, the LifeModeler(TM) (San Clemente, CA) biomechanics simulation software was used to model a squat exercise. The initial model using default parameters yielded physiologically reasonable hip-joint forces. However, no activation was predicted in some large muscles such as rectus femoris, which have been shown to be active in 1-g performance of the activity. Parametric testing was conducted using Monte Carlo methods and combinatorial reduction to find a muscle parameter set that more closely matched physiologically observed activation patterns during the squat exercise. Peak hip joint force using the default parameters was 2.96 times body weight (BW) and increased to 3.21 BW in an optimized, feature-selected test case. The rectus femoris was predicted to peak at 60.1% activation following muscle recruitment optimization, compared to 19.2% activation with default parameters. These results indicate the critical role that muscle parameters play in joint force estimation and the need for exploration of the solution space to achieve physiologically realistic muscle activation.

Hanson, Andrea↗
NASA NTRS2017

Energy Release in the Solar Atmosphere from a Stream of Infalling Prominence Debris

Recent high-resolution and high-cadence extreme-ultraviolet (EUV) imaging has revealed a new phenomenon, impacting prominence debris, where prominence material from failed or partial eruptions can impact the lower atmosphere, releasing energy. We report a clear example of energy release and EUV brightening due to infalling prominence debris that occurred on 2011 September 7-8. The initial eruption of material was associated with an X1.8-class flare from AR 11283, occurring at 22:30 UT on 2011 September 7. Subsequently, a semicontinuous stream of this material returned to the solar surface with a velocity v greater than 150 km/s, impacting a region remote from the original active region between 00:20 and 00:40 UT on 2011 September 8. Using the Solar Dynamics Observatory/Atmospheric Imaging Assembly, the differential emission measure of the plasma was estimated throughout this brightening event. We found that the radiated energy of the impacted plasma was L(sub rad) approx. 10(exp 27) erg, while the thermal energy peaked at approx.10(exp 28) erg. From this we were able to determine the mass content of the debris to be in the range 2 x 10(exp 14) less than m less than 2 x 10(exp 15) g. Given typical prominence masses, the likely debris mass is toward the lower end of this range. This clear example of a prominence debris event shows that significant energy release takes place during these events and that such impacts may be used as a novel diagnostic tool for investigating prominence material properties.

Inglis, A. R.↗
NASA NTRS2019

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

A NASA robotic refueling mission experiment on board the International Space Station (ISS) was designed to repeatedly transfer simulated cryogenic propellant between two dewars. After each metered transfer, the cryogen was to be vented overboard via sublimation or evaporation. ISS payload providers must conduct analyses to demonstrate that any planned gaseous venting generate no more than a certain level of material that may interfere with optical measurements from other payloads that may be located nearby. This requirement is expressed in terms of a maximum column number density (CND). For the experiment under consideration, fluid droplets likely would accompany any such vapor. Earlier related studies led to development of analytical expressions for column density along general paths for a number of practical cases under steady conditions, including droplet evaporation. These expressions are not appropriate for rapid evaporation of small droplets however, since the droplets cannot sustain mass flow rates consistent with steady conditions. In this work, analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where it is found the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next, the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction, but at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed in this framework.

Woronowicz, M. S.↗
DOE OSTI2020

DEMONSTRATION OF A DATA-DRIVEN PHYSICS-BASED APPROACH FOR COMPUTATIONALLY EFFICIENT CFD PREDICTION OF TWO-PHASE BUBBLY FLOW

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale physics-guided data-driven approach, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. By defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

42 ENGINEERING↗
NASA NTRS1994

Gamma-ray burst constraints on the galactic frequency of extra-solar Oort clouds

With the strong CGRO/BATSE evidence that most gamma-ray bursts do not come from galactic neutron stars, models involving the accretion of a comet onto a neutron star (NS) no longer appear to be strong contenders for explaining the majority of bursts. If this is the case, then it is worth asking whether the lack of an observed galactic gamma-ray burst population provides a useful constraint on the number of comets and comet clouds in the galaxy. Owing to the previously unrecognized structural weakness of cometary nuclei, we find the capture cross sections for comet-NS events to be much higher than previously published estimates, with tidal breakup at distances R(sub b) approximately equals to 4 x 10(exp 10) cm from the NS. As a result, impacts of comets onto field NS's penetrating the Oort Clouds of other stars are found to dominate all other galactic NS-comet capture rates by a factor of 100. This in turn predicts that if comet clouds are common, there should be a significant population of repeater sources with (1) a galactic distribution, (2) space-correlated repetition, and (3) a wide range of peak luminosities and luminosity time histories. If all main sequences stars have Oort Clouds like our own, we predict approximately 4000 such repeater sources in the Milky Way at any time, each repeating on timescales of months to years. Based on estimates of the sensitivity of the CGRO/BATSE instrument and assuming isotropic gamma-ray beaming from such events, we estimate that a population of approximately 20-200 of these galactic NS-Oort Cloud gamma-ray repeater sources should be detectable by CGRO. In addition, if giant planet formation is common in the galaxy, we estimate that the accretion of isolated comets injected to the interstellar medium by giant planet formation should produce an additional source of galactic, nonrepeating events. Comparing these estimates to the three to four soft gamma-ray repeater sources detected by BATSE, one is forced to conclude that (1) comet impacts on NS's are inefficient at producing gamma-rays; or (2) the gamma-rays from such events are highly beamed; or (3) the fraction of stars in the galaxy with Oort Cloud like our own is not higher than a few percent.

Shull, J. Michael↗
DOE OSTI2025

Achieving a large transverse emittance ratio in the Relativistic Heavy Ion Collider: Beam experiments, analytical estimates, and implications for the Electron-Ion Collider

The Electron-Ion Collider (EIC), to be built at Brookhaven National Laboratory, will collide a polarized high-energy electron beam with polarized proton and light ion beams in the center-of-mass energy range of 20–140 GeV, achieving a maximum peak luminosity of up to 1 × 10 34 cm − 2 s − 1 for the proton and electron collisions. To achieve such a high luminosity, the EIC will employ small, flat beams at the interaction point, which requires a transverse emittance ratio of about 11:1. Since 2017, we have carried out a series of beam experiments to test the feasibility of achieving such a large emittance ratio in the Relativistic Heavy Ion Collider (RHIC). Thanks to the existing stochastic cooling and decoupling feedback system in RHIC, we successfully demonstrated a large transverse emittance ratio of 11:1 with gold ion beams in 2023. In this article, we will present more details and data analysis for these experiments, along with analytical estimates and numerical multiparticle simulations, and provide an outlook for future EIC coupling compensation. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2023

Estimating the Subsolar Magnetopause Position from Soft X-Ray Images Using a Low-Pass Image Filter

The Lunar Environment heliospheric X-ray Imager (LEXI) and Solar wind Magnetosphere Ionosphere Link Explorer (SMILE) missions will image the Earth’s dayside magnetopause and cusps in soft X-rays after their respective launches in the near future, to specify global magnetic reconnection modes for varying solar wind conditions. To support the success of these scientific missions, it is critical to develop techniques that extract the magnetopause locations from the observed soft X-ray images. In this research, we introduce a new geometric equation that calculates the subsolar magnetopause position ( R s ) from a satellite position, the look direction of the instrument, and the angle at which the X-ray emission is maximized. Two assumptions are used in this method: (1) The look direction where soft X-ray emissions are maximized lies tangent to the magnetopause, and (2) the magnetopause surface near the subsolar point is almost spherical and thus R s is nearly equal to the radius of the magnetopause curvature. We create synthetic soft X-ray images by using the Open Geospace General Circulation Model (OpenGGCM) global magnetohydrodynamic model, the galactic background, the instrument point spread function, and Poisson noise. We then apply the fast Fourier transform and Gaussian low-pass filters to the synthetic images to remove noise and obtain accurate look angles for the soft X-ray peaks. From the filtered images, we calculate R 2 and its accuracy for different LEXI locations, look directions, and solar wind densities by using the OpenGGCM subsolar magnetopause location as ground truth. Our method estimates R s with an accuracy of <0.3 R E when the solar wind density exceeds >10 cm -3 . The accuracy improves for greater solar wind densities and during southward interplanetary magnetic fields. The method captures the magnetopause motion during southward interplanetary magnetic field turnings. Consequently, the technique will enable quantitative analysis of the magnetopause motion and help reveal the dayside reconnection modes for dynamic solar wind conditions. This technique will support the LEXI and SMILE missions in achieving their scientific objectives.

Hyangpyo Kim↗
DOE OSTI2026

Machine-learning-based estimates of global natural vegetated wetland methane emissions (2000–2025)

Wetlands are the largest natural source of atmospheric methane (CH 4 ), yet comprehensive global budgets are typically delayed by years, preventing a timely understanding of CH 4 sources, sinks, and trends. To reduce this delay, we present a model emulator-driven framework and accompanying workflow that enable timely, continuous emission updates using a machine-learning emulator to reconstruct spatially explicit monthly emission fields at 1° × 1° resolution. We apply this framework to a global dataset of natural vegetated wetland CH 4 emissions to extend the most recent Global Methane Budget (GMB; Saunois et al., 2025) record that covers the 2000–2020 emissions through 2025. In the test data (∼ 30 % of the total dataset), the emulator achieved a global R 2 of 0.65 ± 0.003 (mean ± 95 % CI, hereafter) and an RMSE of 5.49 ± 0.12×10 -3 Tg CH 4 yr −1 . The emulator is trained on 35 GMB model estimates, including 22 process-based models and 13 atmospheric inversions, paired with 10 ensemble realizations of 11 gridded climate predictor variables from atmospheric reanalyses. Our results show that the global mean predicted wetland CH 4 emissions for 2021–2025 (157.8 ± 2.4 Tg CH 4 yr −1 ) are not significantly higher (∼ 0.05 Tg CH 4 yr −1 ) than the 2000–2020 baseline. However, this stability masks a significant hemispheric redistribution of emissions. We detect an increase in Northern Hemisphere (NH) emissions in 2021–2025, with mid- and high-latitudes increasing by 0.76 ± 0.07 and 0.35 ± 0.03 Tg CH 4 yr −1 , respectively, while the tropics and Southern Hemisphere (SH) extratropics show offsetting negative trends (−0.95 ± 0.19 and -0.11 ± 0.02 Tg CH 4 yr −1 , respectively). The predicted emissions are able to capture the low emissions in 2023 in South America linked to El Niño-related drought, as reported by recent studies (Ciais et al., 2026; Quinn et al., 2025). Furthermore, we identify a distinct seasonal amplification of global emission trends that peaks in late boreal summer. This new modeled dataset and operational framework bridge the gap between the latest updated budgets and low-latency monitoring, providing a scalable capacity to frequently update global emission estimates and critical early warnings of regional wetland feedback loops. The data are publicly available at https://doi.org/10.5281/zenodo.18870108 (Li et al., 2026).

Li, Mengze [National University of Singapore (Sing↗
NASA NTRS1995

Gamma-ray burst constraints on the galactic frequency of extrasolar Oort Clouds

With the strong Compton Gamma-Ray Observatory/Burst and Transient Source Experiment (CGRO/BATSE) evidence that most gamma-ray bursts do not come from galactic neutron stars, models involving the accretion of a comet onto a neutron star (NS) no longer appear to be strong contenders for explaining the majority of bursts. If this is the case, then it is worth asking whether the lack of an observed galactic gamma-ray burst population provides a useful constraint on the number of comets and comet clouds in the galaxy. Owing to the previously unrecognized structural weakness of cometary nuclei, we find the capture cross sections for comet-NS events to be much higher than previously published estimates, with tidal breakup at distances R(sub b) approx. equals 4 x 10(exp 10) cm from the NS. As a result, impacts of comets onto field NSs penetrating the Oort Clouds of other stars are found to dominate all other galactic NS-comet capture rates by a factor of 100. This in turn predicts that if comet clouds are common, there should be a significant population of repeater sources with (1) a galactic distribution, (2) space-correlated repetition, and (3) a wide range of peak luminosities and luminosity time histories. If all main sequence stars have Oort Clouds like our own, we predict approximately 4000 such repeater sources in the Milky Way at any time, each repeating on time scales of months to years. Based on estimates of the sensitivity of the CGRO/BATSE instrument and assuming isotropic gamma-ray beaming from such events, we estimate that a population of approximately 20-200 of these galactic NS-Oort Cloud gamma-ray repeater sources should be detectable by CGRO. In addition, if giant planet formation is common in the galaxy, we estimate that the accretion of isolated comets injected to the interstellar medium by giant planet formation should produce an additional source of galactic, nonrepeating, events. Comparing these estimates to the 3-4 soft gamma-ray repeater sources detected by BATSE, one is forced to conclude that (1) comet impacts on NSs are inefficient at producing gamma rays; or (2) the gamma rays from such events are highly beamed; or (3) the fraction of stars in the galaxy with Oort Clouds like our own is not higher than a few percent.

Shull, J. Michael↗
DOE OSTI2021

Modeling Within-Host Dynamics of SARS-CoV-2 Infection: A Case Study in Ferrets

The pre-clinical development of antiviral agents involves experimental trials in animals and ferrets as an animal model for the study of SARS-CoV-2. Here, we used mathematical models and experimental data to characterize the within-host infection dynamics of SARS-CoV-2 in ferrets. We also performed a global sensitivity analysis of model parameters impacting the characteristics of the viral infection. We provide estimates of the viral dynamic parameters in ferrets, such as the infection rate, the virus production rate, the infectious virus proportion, the infected cell death rate, the virus clearance rate, as well as other related characteristics, including the basic reproduction number, pre-peak infectious viral growth rate, post-peak infectious viral decay rate, pre-peak infectious viral doubling time, post-peak infectious virus half-life, and the target cell loss in the respiratory tract. These parameters and indices are not significantly different between animals infected with viral strains isolated from the environment and isolated from human hosts, indicating a potential for transmission from fomites. While the infection period in ferrets is relatively short, the similarity observed between our results and previous results in humans supports that ferrets can be an appropriate animal model for SARS-CoV-2 dynamics-related studies, and our estimates provide helpful information for such studies.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2003

Multiple Input Design for Real-Time Parameter Estimation in the Frequency Domain

A method for designing multiple inputs for real-time dynamic system identification in the frequency domain was developed and demonstrated. The designed inputs are mutually orthogonal in both the time and frequency domains, with reduced peak factors to provide good information content for relatively small amplitude excursions. The inputs are designed for selected frequency ranges, and therefore do not require a priori models. The experiment design approach was applied to identify linear dynamic models for the F-15 ACTIVE aircraft, which has multiple control effectors.

Morelli, Eugene↗
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↗