Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 361 records · Page 20

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependant flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. This Technical Memorandum illustrates the application and verification of the code through 25 demonstrated example problems.

Majumdar, A. K.↗

Component-specific modeling

Accomplishments are described for a 3 year program to develop methodology for component-specific modeling of aircraft hot section components (turbine blades, turbine vanes, and burner liners). These accomplishments include: (1) engine thermodynamic and mission models, (2) geometry model generators, (3) remeshing, (4) specialty three-dimensional inelastic structural analysis, (5) computationally efficient solvers, (6) adaptive solution strategies, (7) engine performance parameters/component response variables decomposition and synthesis, (8) integrated software architecture and development, and (9) validation cases for software developed.

Mcknight, R. L.↗

Performance Modeling of Experimental Laser Lightcrafts

A computational plasma aerodynamics model is developed to study the performance of a laser propelled Lightcraft. The computational methodology is based on a time-accurate, three-dimensional, finite-difference, chemically reacting, unstructured grid, pressure-based formulation. The underlying physics are added and tested systematically using a building-block approach. The physics modeled include non-equilibrium thermodynamics, non-equilibrium air-plasma finite-rate kinetics, specular ray tracing, laser beam energy absorption and refraction by plasma, non-equilibrium plasma radiation, and plasma resonance. A series of transient computations are performed at several laser pulse energy levels and the simulated physics are discussed and compared with those of tests and literatures. The predicted coupling coefficients for the Lightcraft compared reasonably well with those of tests conducted on a pendulum apparatus.

Wang, Ten-See↗

Generalized Fluid System Simulation Program (GFSSP) - Version 6

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, flow control valves and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. Users can introduce new physics, non-linear and time-dependent boundary conditions through user-subroutine.

Majumdar, Alok↗

Performance Modeling of an Experimental Laser Propelled Lightcraft

A computational plasma aerodynamics model is developed to study the performance of an experimental laser propelled lightcraft. The computational methodology is based on a time-accurate, three-dimensional, finite-difference, chemically reacting, unstructured grid, pressure- based formulation. The underlying physics are added and tested systematically using a building-block approach. The physics modeled include non-equilibn'um thermodynamics, non-equilibrium air-plasma finite-rate kinetics, specular ray tracing, laser beam energy absorption and equi refraction by plasma, non-equilibrium plasma radiation, and plasma resonance. A series of transient computations are performed at several laser pulse energy levels and the simulated physics are discussed and compared with those of tests and literature. The predicted coupling coefficients for the lightcraft compared reasonably well with those of tests conducted on a pendulum apparatus.

Wang, Ten-See↗

A non-local thermodynamic equilibrium, line-blanketed synthetic spectrum of Iota Herculis - C, Al, and Si lines

A non-LTE line-blanketed model stellar atmosphere is used to compute a model of I Herculis (B3 IV) with a Teff of 17,500 K and a log g of 3.75, following the conclusions of Peters and Polidan (1985). Detailed profiles of a number of lines of C, Al, and Si in the 1200-2000-A region are computed, including the resonance lines of C II, Al II, and Al III. These profiles are compared to observations obtained from the coaddition of eight IUE SWP images, using a technique developed by Leckrone and Adelman (1989). Comparison of carbon lines with a model that is underabundant in carbon by a factor of 2 relative to the sun indicates that the C abundance of Iota Her is at most one-half solar. Non-LTE effects are examined by comparing an LTE model possessing identical atmospheric parameters with the non-LTE model. Substantial differences in the populations of the model atomic states are found, but differences in the temperature structure of the two models often mask the non-LTE effects in the synthetic spectra.

Grigsby, James A.↗

Thermodynamic and fluid mechanic analysis of rapid pressurization in a dead-end tube

Three models have been applied to very rapid compression of oxygen in a dead-ended tube. Pressures as high as 41 MPa (6000 psi) leading to peak temperatures of 1400 K are predicted. These temperatures are well in excess of the autoignition temperature (750 K) of teflon, a frequently used material for lining hoses employed in oxygen service. These findings are in accord with experiments that have resulted in ignition and combustion of the teflon, leading to the combustion of the stainless steel braiding and catastrophic failure. The system analyzed was representative of a capped off-high-pressure oxygen line, which could be part of a larger system. Pressurization of the larger system would lead to compression in the dead-end line, and possible ignition of the teflon liner. The model consists of a large plenum containing oxygen at the desired pressure (500 to 6000 psi). The plenum is connected via a fast acting valve to a stainless steel tube 2 cm inside diameter. Opening times are on the order of 15 ms. Downstream of the valve is an orifice sized to increase filling times to around 100 ms. The total length from the valve to the dead-end is 150 cm. The distance from the valve to the orifice is 95 cm. The models describe the fluid mechanics and thermodynamics of the flow, and do not include any combustion phenomena. A purely thermodynamic model assumes filling to be complete upstream of the orifice before any gas passes through the orifice. This simplification is reasonable based on experiment and computer modeling. Results show that peak temperatures as high as 4800 K can result from recompression of the gas after expanding through the orifice. An approximate transient model without an orifice was developed assuming an isentropic compression process. An analytical solution was obtained. Results indicated that fill times can be considerably shorter than valve opening times. The third model was a finite difference, 1-D transient compressible flow model. Results from the code show the recompression effect but predict much lower peak temperatures than the thermodynamic model.

Leslie, Ian H.↗

An investigation of the energy balance of solar active regions using the ACRIM irradiance data

The detection of a significant correlation between the solar irradiance, corrected for flux deficit due to sunspots, and both the 205 nm flux and a photometric facular index were examined. A detailed analysis supports facular emission as the more likely source of correlation with the corrected radiance, rather then the error in sunspot correction. A computer program which simulates two dimensional convection in a compressible, stratified medium was investigated. Subroutines to calculate ionization and other thermodynamic variables were also completed.

Petro, L. D.↗

Directional solidification with heat losses

The upward directional solidification of a dilute binary alloy in a tall, thin, two-dimensional slot is considered. Heat losses from the sidewalls produce curved isotherms, interface deformations, and melt convection. An asymptotic approach is used which relates the small heat loss to the thinness of the solidification cell, and solutions are obtained for the resulting melt convection, interface deflection, and the distribution of heat and solute. The model is compared to the numerical results of Chang and Brown (1983) which exclude the effects of solute buoyancy and thermodynamic equilibrium at the interface. The radial segregation at the interface (consistent with both interface shapes and temperature fields) is computed, and it is found that the radial segregation is independent of interfacial shapes; the radial segregation cannot be calculated from only a knowledge of the interface shape.

Brattkus, K.↗

Liquefaction of Cryogenic Fluids for Production and Storage of Commodities on Extra-Terrestrial Surfaces

The Cryogenic Fluid In-situ Liquefaction for Landers (CryoFILL) activity has been investigating concepts for the liquefaction of cryogenic fluids produced on the surface of the Moon and Mars. CryoFILL consists of four activities: liquid nitrogen liquefaction testing that was completed in 2019, liquid oxygen liquefaction testing planned to begin in early 2022, fiber-optic sensor temperature measurement system development, and lightweight vacuum jacketed systems designed for the atmosphere of Mars. Additional activities include the development and validation of numerical modeling systems at three stages: thermodynamic models of the liquefaction process for basic energy/mass estimates, nodal models to predict liquefaction tank system level performance, and computational fluid dynamics to assess fluid phenomena occurring within the tank, specifically condensation and stratification within the ullage. These developments form a combined liquefaction and storage solution to support Lunar and Martian exploration. Key testing results from liquid nitrogen testing as well as plans for liquid oxygen testing will be discussed including key sensitivities from analytical evaluation of completed test results. Progress on the development and improvements on fiber optic sensor testing and developments will also be provided. Finally, work on the progress of the lightweight vacuum jacketed systems being co-developed by NASA and industry will be summarized.

Liquefaction↗

Development of a Computational Chemical Vapor Deposition Model: Applications to Indium Nitride and Dicyanovinylaniline

A computational chemical vapor deposition (CVD) model is presented, that couples chemical reaction mechanisms with fluid dynamic simulations for vapor deposition experiments. The chemical properties of the systems under investigation are evaluated using quantum, molecular and statistical mechanics models. The fluid dynamic computations are performed using the CFD-ACE program, which can simulate multispecies transport, heat and mass transfer, gas phase chemistry, chemistry of adsorbed species, pulsed reactant flow and variable gravity conditions. Two experimental setups are being studied, in order to fabricate films of: (a) indium nitride (InN) from the gas or surface phase reaction of trimethylindium and ammonia; and (b) 4-(1,1)dicyanovinyl-dimethylaminoaniline (DCVA) by vapor deposition. Modeling of these setups requires knowledge of three groups of properties: thermodynamic properties (heat capacity), transport properties (diffusion, viscosity, and thermal conductivity), and kinetic properties (rate constants for all possible elementary chemical reactions). These properties are evaluated using computational methods whenever experimental data is not available for the species or for the elementary reactions. The chemical vapor deposition model is applied to InN and DCVA. Several possible InN mechanisms are proposed and analyzed. The CVD model simulations of InN show that the deposition rate of InN is more efficient when pulsing chemistry is used under conditions of high pressure and microgravity. An analysis of the chemical properties of DCVA show that DCVA dimers may form under certain conditions of physical vapor transport. CVD simulations of the DCVA system suggest that deposition of the DCVA dimer may play a small role in the film and crystal growth processes.

Cardelino, Carlos↗

Metal Hydroxide and Oxyhydroxide Gas Database

In a high temperature combustion situation such as a gas turbine or rocket engine, combustion products include water vapor. This water vapor reacts with oxides in the hot-gas-path to form volatile hydroxides and oxyhydroxides. The oxides may be structural oxide components or protective oxides on alloys or non-oxide ceramics. In either case, this enhanced volatilization leads to degradation of components. Reliable thermochemical data for many metal hydroxides and oxyhydroxides has not been available. We review a database put together at NASA over many years based on experimental and ab initio computations for thermochemical data of Al, Si, Ti, Cr, Mn, Y, Zr, La, Ta, Gd, and Yb hydroxides and oxyhydroxides. Experimental methods for measuring thermodynamic quantities of metal hydroxides and oxyhydroxides are limited since the required high temperatures and high pressures preclude many conventional methods. Nonetheless, transpiration has proven a valuable tool in our laboratories for the Si-O-H, Cr-O-H, and Ti-O-H systems. Data are analyzed with second and third law methods. These data can be leveraged with quantum chemistry to provide reliable free energy functions for the third law calculation and comparison to the derived enthalpies of formation. Both wave function and DFT methods can be utilized for deriving the thermodynamic functions for the various metal hydroxides and oxyhydroxides. We have found that, in general, reasonable molecular geometries and spectroscopic data can be derived with a relatively low level of theory, i.e. DFT with the B3LYP functional. Accurate calculations of the enthalpy of formation require higher levels of theory, such as the Coupled Cluster, Singles, Doubles, and perturbative Triples CCSD(T) approach with complete basis sets (CBS). Also, corrections for hindered rotations and anharmonicity are necessary. For some symmetries, the anharmonic correction has proven to be a challenge and we have found that this is only treated properly by certain codes. The database lists the fHo(298), So(298), and Cp in the form Cp = A + BT + CT2 + D/T + E/T2, where A, B, C, D, and E are the fitting constants. Thus, it can be easily used in any of the common computational thermochemistry codes. To date, the database contains 51 gaseous species of the formula MOx(OH)y.

Metal Hydroxides↗

Thermodynamic Data to 20,000 K For Monatomic Gases

This report contains standard-state thermodynamic functions for 50 gaseous atomic elements plus deuterium and electron gas, 51 singly ionized positive ions, and 36 singly ionized negative ions. The data were generated by the NASA Lewis computer program PAC97, a modified version of PAC91 reported in McBride and Gordon. This report is being published primarily to document part of the data currently being used in several NASA Lewis computer programs. The data are presented in tabular and graphical format and are also represented in the form of least-squares coefficients. The tables give the following data as functions of temperature : heat capacity, enthalpy, entropy Gibbs energy, enthalpy of formation, and equilibrium constant. A brief discussion and a comparison of calculated results are given for several models for calculating ideal thermodynamic data for monatomic gases.

Gordon, Sanford↗

A Mass Computation Model for Lightweight Brayton Cycle Regenerator Heat Exchangers

Based on a theoretical analysis of convective heat transfer across large internal surface areas, this paper discusses the design implications for generating lightweight gas-gas heat exchanger designs by packaging such areas into compact three-dimensional shapes. Allowances are made for hot and cold inlet and outlet headers for assembly of completed regenerator (or recuperator) heat exchanger units into closed cycle gas turbine flow ducting. Surface area and resulting volume and mass requirements are computed for a range of heat exchanger effectiveness values and internal heat transfer coefficients. Benefit cost curves show the effect of increasing heat exchanger effectiveness on Brayton cycle thermodynamic efficiency on the plus side, while also illustrating the cost in heat exchanger required surface area, volume, and mass requirements as effectiveness is increased. The equations derived for counterflow and crossflow configurations show that as effectiveness values approach unity, or 100 percent, the required surface area, and hence heat exchanger volume and mass tend toward infinity, since the implication is that heat is transferred at a zero temperature difference. To verify the dimensional accuracy of the regenerator mass computational procedure, calculation of a regenerator specific mass, that is, heat exchanger weight per unit working fluid mass flow, is performed in both English and SI units. Identical numerical values for the specific mass parameter, whether expressed in lb/(lb/sec) or kg/(kg/sec), show the dimensional consistency of overall results.

Juhasz, Albert J.↗

A Mass Computation Model for Lightweight Brayton Cycle Regenerator Heat Exchangers

Based on a theoretical analysis of convective heat transfer across large internal surface areas, this paper discusses the design implications for generating lightweight gas-gas heat exchanger designs by packaging such areas into compact three-dimensional shapes. Allowances are made for hot and cold inlet and outlet headers for assembly of completed regenerator (or recuperator) heat exchanger units into closed cycle gas turbine flow ducting. Surface area and resulting volume and mass requirements are computed for a range of heat exchanger effectiveness values and internal heat transfer coefficients. Benefit cost curves show the effect of increasing heat exchanger effectiveness on Brayton cycle thermodynamic efficiency on the plus side, while also illustrating the cost in heat exchanger required surface area, volume, and mass requirements as effectiveness is increased. The equations derived for counterflow and crossflow configurations show that as effectiveness values approach unity, or 100 percent, the required surface area, and hence heat exchanger volume and mass tend toward infinity, since the implication is that heat is transferred at a zero temperature difference. To verify the dimensional accuracy of the regenerator mass computational procedure, calculation of a regenerator specific mass, that is, heat exchanger weight per unit working fluid mass flow, is performed in both English and SI units. Identical numerical values for the specific mass parameter, whether expressed in lb/(lb/sec) or kg/ (kg/sec), show the dimensional consistency of overall results.

Juhasz, Albert J.↗

Modeling Commercial Turbofan Engine Icing Risk With Ice Crystal Ingestion

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which are ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in flight. The computational tool was utilized to help guide a portion of the PSL testing, and was used to predict ice accretion could also occur at significantly lower altitudes. The predictions were qualitatively verified by subsequent testing of the engine in the PSL. The PSL test has helped to calibrate the engine icing computational tool to assess the risk of ice accretion. The results from the computer simulation identified prevalent trends in wet bulb temperature, ice particle melt ratio, and engine inlet temperature as a function of altitude for predicting engine icing risk due to ice crystal ingestion.

Jorgenson, Philip C. E.↗

Modeling of Commercial Turbofan Engine With Ice Crystal Ingestion: Follow-On

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which is ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in flight. The computational tool was utilized to help guide a portion of the PSL testing, and was used to predict ice accretion could also occur at significantly lower altitudes. The predictions were qualitatively verified by subsequent testing of the engine in the PSL. In a previous study, analysis of select PSL test data points helped to calibrate the engine icing computational tool to assess the risk of ice accretion. This current study is a continuation of that data analysis effort. The study focused on tracking the variations in wet bulb temperature and ice particle melt ratio through the engine core flow path. The results from this study have identified trends, while also identifying gaps in understanding as to how the local wet bulb temperature and melt ratio affects the risk of ice accretion and subsequent engine behavior.

Jorgenson, Philip C. E.↗

Modeling of Commercial Turbofan Engine with Ice Crystal Ingestion; Follow-On

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which is ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in flight. The computational tool was utilized to help guide a portion of the PSL testing, and was used to predict ice accretion could also occur at significantly lower altitudes. The predictions were qualitatively verified by subsequent testing of the engine in the PSL. In a previous study, analysis of select PSL test data points helped to calibrate the engine icing computational tool to assess the risk of ice accretion. This current study is a continuation of that data analysis effort. The study focused on tracking the variations in wet bulb temperature and ice particle melt ratio through the engine core flow path. The results from this study have identified trends, while also identifying gaps in understanding as to how the local wet bulb temperature and melt ratio affects the risk of ice accretion and subsequent engine behavior.

Jorgenson, Philip C. E.↗