Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic simulation”

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 73 records · Page 4

Enabling Nuclear and Solar Thermal Propulsion with the Computational Laboratory: Faster and Cheaper Materials Characterization

Nuclear and solar thermal propulsion offer enhanced efficiencies for in-space travel which may open up the outer Solar System to both crewed and automated craft. However, these advanced propulsion systems require cutting-edge materials which are often not easily tested in the laboratory. Ab initio thermodynamics and other computational methods offer a way to test materials at high pressures and temperatures more quickly, and with reduced cost. Here, several case studies are presented which highlight the utility of modern simulation techniques. First, an extensive thermodynamic study was performed at the ab initio level to elucidate the role of hot hydrogen flow on the erosion of four refractory carbides to examine their appropriateness for use as channel coatings. A detailed analysis of reaction products informs estimates of erosion which are validated in comparison with heritage NERVA data. Specific material suggestions are made for both nuclear and solar designs. Secondly, we examined the mechanical behavior and chemical stability of various components of an idealized nuclear cermet core design using techniques spanning from the atomic to microscale. Next, the mechanical and chemical behavior of tungsten (a potential matrix material) and the stability of uranium mononitride (a potential ceramic fuel element), are characterized in the presence of hydrogen and at elevated temperatures with ab initio thermodynamics. Finally, the mechanical response and failure of tungsten under tensile strain is simulated at the micron scale with a combined finite element and dislocation dynamics approach. The resultant stress-strain data agrees well with experiments and leads to a workflow which can incorporate ab initio thermodynamic data into micron-scale simulations and thus provide real-world relevant estimates of erosion and stress-strain behavior. The power of this framework and our plans to use it in the future will be discussed.

William C Tucker↗

Probabilistic Analysis of Gas Turbine Field Performance

A gas turbine thermodynamic cycle was computationally simulated and probabilistically evaluated in view of the several uncertainties in the performance parameters, which are indices of gas turbine health. Cumulative distribution functions and sensitivity factors were computed for the overall thermal efficiency and net specific power output due to the thermodynamic random variables. These results can be used to quickly identify the most critical design variables in order to optimize the design, enhance performance, increase system availability and make it cost effective. The analysis leads to the selection of the appropriate measurements to be used in the gas turbine health determination and to the identification of both the most critical measurements and parameters. Probabilistic analysis aims at unifying and improving the control and health monitoring of gas turbine aero-engines by increasing the quality and quantity of information available about the engine's health and performance.

Gorla, Rama S. R.↗

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between liquid and vapor phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures, temperatures, and pressurant requirements agree with experimental data, with pressurant mass predictions within 16\% of observed values. This study highlights the importance of selecting appropriate turbulence models to accurately simulate complex flow and heat transfer phenomena during tank pressurization and expulsion. By enhancing the accuracy and reliability of CFD models for liquid hydrogen under cryogenic conditions, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between vapor and liquid phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved autogenous pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures and temperatures agree well with experimental data, whereas the pressurant mass predictions are within 16% of observed values. This study emphasizes the importance of selecting appropriate turbulence models to accurately simulate the complex flow and heat transfer phenomena during tank autogenous pressurization and expulsion. By improving the accuracy and reliability of CFD models for these processes, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗

Modeling the Deterioration of Engine and Low Pressure Compressor Performance During a Roll Back Event Due to Ice Accretion

The main focus of this study is to apply a computational tool for the flow analysis of the engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) of NASA Glenn Research Center. A data point was selected for analysis during which the engine experienced a full roll back event due to the ice accretion on the blades and flow path of the low pressure compressor. The computational tool consists of the Numerical Propulsion System Simulation (NPSS) engine system thermodynamic cycle code, and an Euler-based compressor flow analysis code, that has an ice particle melt estimation code with the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Decreasing the performance characteristics of the low pressure compressor (LPC) within the NPSS cycle analysis resulted in matching the overall engine performance parameters measured during testing at data points in short time intervals through the progression of the roll back event. Detailed analysis of the fan-core and LPC with the compressor flow analysis code simulated the effects of ice accretion by increasing the aerodynamic blockage and pressure losses through the low pressure compressor until achieving a match with the NPSS cycle analysis results, at each scan. With the additional blockages and losses in the LPC, the compressor flow analysis code results were able to numerically reproduce the performance that was determined by the NPSS cycle analysis, which was in agreement with the PSL engine test data. The compressor flow analysis indicated that the blockage due to ice accretion in the LPC exit guide vane stators caused the exit guide vane (EGV) to be nearly choked, significantly reducing the air flow rate into the core. This caused the LPC to eventually be in stall due to increasing levels of diffusion in the rotors and high incidence angles in the inlet guide vane (IGV) and EGV stators. The flow analysis indicating compressor stall is substantiated by the video images of the IGV taken during the PSL test, which showed water on the surface of the IGV flowing upstream out of the engine, indicating flow reversal, which is characteristic of a stalled compressor.

Veres, Joseph P.↗

Modeling the Deterioration of Engine and Low Pressure Compressor Performance During a Roll Back Event Due to Ice Accretion

The main focus of this study is to apply a computational tool for the flow analysis of the engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) of NASA Glenn Research Center. A data point was selected for analysis during which the engine experienced a full roll back event due to the ice accretion on the blades and flow path of the low pressure compressor. The computational tool consists of the Numerical Propulsion System Simulation (NPSS) engine system thermodynamic cycle code, and an Euler-based compressor flow analysis code, that has an ice particle melt estimation code with the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Decreasing the performance characteristics of the low pressure compressor (LPC) within the NPSS cycle analysis resulted in matching the overall engine performance parameters measured during testing at data points in short time intervals through the progression of the roll back event. Detailed analysis of the fan-core and LPC with the compressor flow analysis code simulated the effects of ice accretion by increasing the aerodynamic blockage and pressure losses through the low pressure compressor until achieving a match with the NPSS cycle analysis results, at each scan. With the additional blockages and losses in the LPC, the compressor flow analysis code results were able to numerically reproduce the performance that was determined by the NPSS cycle analysis, which was in agreement with the PSL engine test data. The compressor flow analysis indicated that the blockage due to ice accretion in the LPC exit guide vane stators caused the exit guide vane (EGV) to be nearly choked, significantly reducing the air flow rate into the core. This caused the LPC to eventually be in stall due to increasing levels of diffusion in the rotors and high incidence angles in the inlet guide vane (IGV) and EGV stators. The flow analysis indicating compressor stall is substantiated by the video images of the IGV taken during the PSL test, which showed water on the surface of the IGV flowing upstream out of the engine, indicating flow reversal, which is characteristic of a stalled compressor.

Jorgenson, Philip C. E.↗

Performance Analysis of Optimized STARC-ABL Designs Across the Entire Mission Profile

Boundary layer ingestion (BLI) offers the potential for significant fuel burn reduction by exploiting strong aeropropulsive interactions. NASA’s STARC–ABL concept uses an electri- cally powered BLI tail cone thruster on what is otherwise a traditional airframe. Despite the traditional airframe of this configuration, aeropropulsive integration is critical to the perfor- mance of the BLI propulsor. Furthermore, due to being electrically powered, the fan pressure ratio and efficiency of the BLI tail cone thruster vary widely across the flight envelope, and this variation in fan performance must be accounted for with the aeropropulsive integration of the BLI system. Thus, accurate performance prediction for this novel propulsion configu- ration requires the use of a coupled aeropropulsive model across the flight envelope. In this work, we analyze the off-design performance of 18 optimized designs using an aeropropulsive model that is built with the OpenMDAO framework to couple 3-D RANS CFD simulations to 1-D thermodynamic cycle analyses. The designs are created via high-fidelity aeropropulsive design optimizations that span a range of fan pressure ratio and thrust values at the cruise conditions for the STARC-ABL concept, which was chosen as the aerodynamic design point for the propulsor. Performance analyses we present herein are then performed at a range of off-design flight conditions that span the flight envelope, including low-speed and low-altitude flight conditions. This study provides the first set of high-fidelity data for the STARC–ABL configuration at off-design conditions, and the results quantify the power savings through BLI compared to a traditional propulsion system across the entire mission profile.

optimization↗

Impacts of the Convective Transport Algorithm on Atmospheric Composition and Ozone-Climate Feedbacks in GEOS-CCM

Convective transport is one of the dominant factors in determining the composition of the troposphere. It is the main mechanism for lofting constituents from near-surface source regions to the middle and upper troposphere, where they can subsequently be advected over large distances. Gases reaching the upper troposphere can also be injected through the tropopause and play a subsequent role in the lower stratospheric ozone balance. Convection codes in climate models remain a great source of uncertainty for both the energy balance of the general circulation and the transport of constituents. This study uses the Goddard Earth Observing System Chemistry-Climate Model (GEOS CCM) to perform a controlled experiment that isolates the impact of convective transport of constituents from the direct changes on the atmospheric energy balance. Two multi-year simulations are conducted. In the first, the thermodynamic variable, moisture, and all trace gases are transported using the multi-plume Relaxed-Arakawa-Schubert (RAS) convective parameterization. In the second simulation, RAS impacts the thermodynamic energy and moisture in this standard manner, but all other constituents are transported differently. The accumulated convective mass fluxes (including entrainment and detrainment) computed at each time step of the GCM are used with a diffusive (bulk) algorithm for the vertical transport, which above all is less efficient at transporting constituents from the lower to the upper troposphere. Initial results show the expected differences in vertical structure of trace gases such as carbon monoxide, but also show differences in lower stratospheric ozone, in a region where it can potentially impact the climate state of the model. This work will investigate in more detail the impact of convective transport changes by comparing the two simulations over many years (1996-2010), focusing on comparisons with observed constituent distributions and similarities and differences of patterns of inter-annual variability caused by the convective transport algorithm. In particular, the impact on lower stratospheric composition will be isolated and the subsequent feedbacks of ozone on the climate forcing and tropopause structure will be assessed.

Pawson, S.↗

A Damage Model for the Simulation of Delamination in Advanced Composites under Variable-Mode Loading

A thermodynamically consistent damage model is proposed for the simulation of progressive delamination in composite materials under variable-mode ratio. The model is formulated in the context of Damage Mechanics. A novel constitutive equation is developed to model the initiation and propagation of delamination. A delamination initiation criterion is proposed to assure that the formulation can account for changes in the loading mode in a thermodynamically consistent way. The formulation accounts for crack closure effects to avoid interfacial penetration of two adjacent layers after complete decohesion. The model is implemented in a finite element formulation, and the numerical predictions are compared with experimental results obtained in both composite test specimens and structural components.

Turon, A.↗

Numerical Simulation of North Atlantic Sea Ice Variability, 1951 - 1980

A two-level dynamic-thermodynamic sea ice model is used to simulate the growth, drift and decay of sea ice in the Northern Hemisphere during a 30-year period, 1951 to 1980. The model is run with a daily timestep on a 222 km grid and is forced by interanually varying fields of geostrophic wind and temperature-derived thermodynamic fluxes. The objective is a quantitative description of large-scale sea ice variability in terms of the dynamic and thermodynamic processes responsible for the fluctuations, especially in the North Atlantic where sea ice represents a substantial input of fresh water. The fields of ice velocity and thickness contain strong seasonal as well as interannual variability. The mean drift pattern results in thicknesses of 4 to 5 m offshore of northern Canada and Greenland, while winter thicknesses of approximately 2 m are typical of Alaskan. Eurasian and East Greenland coastal waters. The 30-year mean fields are characterized by ecessive ice in the North Atlantic during winter and by a summer retreat that is more rapid than observed.

Walsh, J. E.↗

An Introduction to Thermodynamic Performance Analysis of Aircraft Gas Turbine Engine Cycles Using the Numerical Propulsion System Simulation Code

This document is intended as an introduction to the analysis of gas turbine engine cycles using the Numerical Propulsion System Simulation (NPSS) code. It is assumed that the analyst has a firm understanding of fluid flow, gas dynamics, thermodynamics, and turbomachinery theory. The purpose of this paper is to provide for the novice the information necessary to begin cycle analysis using NPSS. This paper and the annotated example serve as a starting point and by no means cover the entire range of information and experience necessary for engine performance simulation. NPSS syntax is presented but for a more detailed explanation of the code the user is referred to the NPSS User Guide and Reference document (ref. 1).

Jones, Scott M.↗

Aspects of Supercritical Turbulence: Direct Numerical Simulation of O2/H2 and C7H16/N2 Temporal Mixing Layers

Results from Direct Numerical Simulations of temporal, supercritical mixing layers for two species systems are analyzed to elucidate species-specific turbulence aspects. The two species systems, O2/H2 and C7HG16/N2, have different thermodynamic characteristics; thus, although the simulations are performed at similar reduced pressure (ratio of the pressure to the critical pressure), the former system is dose to mixture ideality and has a relatively high solubility with respect to the latter, which exhibits strong departures from mixture ideality Due to the specified, smaller initial density stratification, the C7H16/N2 layers display higher growth and increased global molecular mixing as well as larger turbulence levels. However, smaller density gradients at the transitional state for the O2/H2 system indicate that on a local basis, the layer exhibits an enhanced mixing, this being attributed to the increased solubility and to mixture ideality. These thermodynamic features are shown to affect the irreversible entropy production (i.e. the dissipation), which is larger for the O2/H2 layer and is primarily concentrated in high density-gradient magnitude regions that are distortions of the initial density stratification boundary. In contrast, the regions of largest dissipation in the C7H16/N2 layer are located in high density-gradient magnitude regions resulting from the mixing of the two fluids.

Bellan, J.↗

Introduction and Overview of the Object-Oriented Turbomachinery Analysis Code (OTAC)

This is an overview of the Object-Oriented Turbomachinery Analysis Code (OTAC). OTAC is an object-oriented code developed to perform 2-D meridional flowthrough design and analysis of turbomachines. OTAC uses an implicit formulation of the governing equations relevant to meanline and streamline simulations to solve for the thermodynamic and velocity state of the fluid at the exit of each blade row. OTAC is written using the interpreted language capability of the Numerical Propulsion System Simulation (NPSS) code; as such, OTAC requires a version of NPSS to run. A series of models and test cases have been developed to demonstrate the capability and features of OTAC.

Turbomachinery↗

Simulation of Delamination Under High Cycle Fatigue in Composite Materials Using Cohesive Models

A new thermodynamically consistent damage model is proposed for the simulation of high-cycle fatigue crack growth. The basis for the formulation is an interfacial degradation law that links Fracture Mechanics and Damage Mechanics to relate the evolution of the damage variable, d, with the crack growth rate da/dN. The damage state is a function of the loading conditions (R and (Delta)G) as well as the experimentally-determined crack growth rates for the material. The formulation ensures that the experimental results can be reproduced by the analysis without the need of additional adjustment parameters.

Camanho, Pedro P.↗

Kinematic and Thermodynamic Structures of Hurricane Bob (1991) Determined From A 1.3-KM Resolution Numerical Simulation

Hurricane Bob (1991) is simulated using the Penn State/NCAR mesoscale model MM5. The simulation is conducted for a 24-h period at 4-km resolution and for a 6-h period at 1.3-km resolution. The 4-km simulation is able to fairly realistically capture the intensity and structure of the storm. The 1.3-km simulation depicts very small-scale convective structures and produces a convective rain band outside of the eye wall that did not occur in the 4-km simulation. The 1.3-km results are used to characterize several kinematic structures in the storm, including low-level outflow in the eye wall, multiple outflows at upper levels, and the convective rain band structure outside of the eye wall. Thermodynamic characteristics of air parcels flowing into and rising within the eye wall will also be examined through trajectory calculations.

Braun, Scott A.↗