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 451 records · Page 25

Calculation of Nonequilibrium Radiation in the Flight Regimes of Aeroassisted Orbital Transfer Vehicles

A computer code has been developed that calculates radiative properties of nonequilibrium air in the low-density regimes expected during the flight of aeroassisted, orbital transfer vehicles. From the given nonequilibrium thermodynamic state variables, the code calculates number densities of internal states and the accompanying emission and absorption characteristics. In addition, the code calculates the number density of the hypothetical gas in radiative equilibrium that produces the same radiation emission as the given nonequilibrium gas. Sample results are shown to demonstrate how the code predicts the non-linear variation of radiation with density at low densities.

Park, Chul↗

Cellular Spacing Selection During the Directional Solidification of Binary Alloys. A Numerical Approach

The evolution of cellular solid/liquid interfaces from an initially unstable planar front was studied by means of a two-dimensional computer simulation. The developed numerical model makes use of an interface tracking procedure and has the capability to describe the dynamics of the interface morphology based on local changes of the thermodynamic conditions. The fundamental physics of this formulation was validated against experimental microgravity results and the predictions of the analytical linear stability theory. The performed simulations revealed that in certain conditions, based on a competitive growth mechanism, an interface could become unstable to random perturbations of infinitesimal amplitude even at wavelengths smaller than the neutral wavelength, lambda(sub c), predicted by the linear stability theory. Furthermore, two main stages of spacing selection have been identified. In the first stage, at low perturbations amplitude, the selection mechanism is driven by the maximum growth rate of instabilities while in the second stage the selection is influenced by nonlinear phenomena caused by the interactions between the neighboring cells. Comparison of these predictions with other existing theories of pattern formation and experimental results will be discussed.

Catalina, Adrian V.↗

Using Satellite Data in Weather Forecasting: I

The GOES Product Generation System (GPGS) is a set of computer codes and scripts that enable the assimilation of real-time Geostationary Operational Environmental Satellite (GOES) data into regional-weather-forecasting mathematical models. The GPGS can be used to derive such geophysical parameters as land surface temperature, the amount of precipitable water, the degree of cloud cover, the surface albedo, and the amount of insolation from satellite measurements of radiant energy emitted by the Earth and its atmosphere. GPGS incorporates a priori information (initial guesses of thermodynamic parameters of the atmosphere) and radiometric measurements from the geostationary operational environmental satellites along with mathematical models of physical principles that govern the transfer of energy in the atmosphere. GPGS solves the radiative-transfer equation and provides the resulting data products in formats suitable for use by weather-forecasting computer programs. The data-assimilation capability afforded by GPGS offers the potential to improve local weather forecasts ranging from 3 hours to 2 days - especially with respect to temperature, humidity, cloud cover, and the probability of precipitation. The improvements afforded by GPGS could be of interest to news media, utility companies, and other organizations that utilize regional weather forecasts.

Jedlovec, Gary J.↗

A Model to Assess the Risk of Ice Accretion Due to Ice Crystal Ingestion in a Turbofan Engine and its Effects on Performance

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 were attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was one or more of the following anomalies: degraded engine performance, engine roll back, compressor surge and stall, and flameout of the combustor. The main focus of this research is the development of a computational tool that can estimate whether there is a risk of ice accretion by tracking key parameters through the compression system blade rows at all engine operating points within the flight trajectory. The tool has an engine system thermodynamic cycle code, coupled with 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 blade rows. Assumptions are made to predict the complex physics involved in engine icing. Specifically, the code does not directly estimate ice accretion and does not have models for particle breakup or erosion. Two key parameters have been suggested as conditions that must be met at the same location for ice accretion to occur: the local wet-bulb temperature to be near freezing or below and the local melt ratio must be above 10%. These parameters were deduced from analyzing laboratory icing test data and are the criteria used to predict the possibility of ice accretion within an engine including the specific blade row where it could occur. Once the possibility of accretion is determined from these parameters, the degree of blockage due to ice accretion on the local stator vane can be estimated from an empirical model of ice growth rate and time spent at that operating point in the flight trajectory. The computational tool can be used to assess specific turbine engines to their susceptibility to ice accretion in an ice crystal environment.

Jorgenson, Philip C. E.↗

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

Stirling Laboratory Research Engine: Preprototype configuration report

The concept of a simple Stirling research engine that could be used by industrial, university, and government laboratories was studied. The conceptual and final designs, hardware fabrication and the experimental validation of a preprototype stirling laboratory research engine (SLRE) were completed. Also completed was a task to identify the potential markets for research engines of this type. An analytical effort was conducted to provide a stirling cycle computer model. The versatile engine is a horizontally opposed, two piston, single acting stirling engine with a split crankshaft drive mechanism; special instrumentation is installed at all component interfaces. Results of a thermodynamic energy balance for the system are reported. Also included are the engine performance results obtained over a range of speeds, working pressures, phase angles and gas temperatures. The potential for a stirling research engine to support the laboratory requirements of educators and researchers was demonstrated.

Hoehn, F. W.↗

Characterization of real gas properties for space shuttle main engine fuel turbine and performance calculations

Real thermodynamic and transport properties of hydrogen, steam, the SSME mixture, and air are developed. The SSME mixture properties are needed for the analysis of the space shuttle main engine fuel turbine. The mixture conditions for the gases, except air, are presented graphically over a temperature range from 800 to 1200 K, and a pressure range from 1 to 500 atm. Air properties are given over a temperature range of 320 to 500 K, which are within the bounds of the thermodynamics programs used, in order to provide mixture data which is more easily checked (than H2/H2O). The real gas property variation of the SSME mixture is quantified. Polynomial expressions, needed for future computer analysis, for viscosity, Prandtl number, and thermal conductivity are given for the H2/H2O SSME fuel turbine mixture at a pressure of 305 atm over a range of temperatures from 950 to 1140 K. These conditions are representative of the SSME turbine operation. Performance calculations are presented for the space shuttle main engine (SSME) fuel turbine. The calculations use the air equivalent concept. Progress towards obtaining the capability to evaluate the performance of the SSME fuel turbine, with the H2/H2O mixture, is described.

Harloff, G. J.↗

A PC-based collapse factor analysis for high pressure cryogenic systems

A collapse factor simulation computer program, suitable for use on personal computers, for supercritical cryogen system pressurization is examined. The available literature is previewed, individual component models are developed, and the overall simulation program logic is examined. Extensive use is made of existing thermodynamic and transport property evaluation packages. Heat transfer is estimated using natural convection correlations, and transient wall heat conduction is considered. Examples of the simulation are demonstrated for two facilities. Partial validation is achieved by comparing the simulation output with data from an existing test facility. Preliminary results indicate that the collapse factors will seldom exceed 1.15 during supercritical operation of high-pressure cryogenic systems. Recommendations for improvements in the component models and the simulation are made.

Hodge, B. K.↗

The CO/SiO radiative instability in cool star atmospheres revisited

We revisit the formation of radiative instabilities in cool star atmospheres and compare our results with those given by Muchmore, Nuth, & Stencel. We have considered the combined influence of CO and SiO molecules and have computed models for a grid of effective temperatures and geometrical dilution factors for the stellar radiation. Our results are based on the analysis of the energy balance of gas elements with prescribed thermodynamic properties. Our results show that radiative instabilities are most likely primarily caused by CO, whereas SiO is expected to play only a minor role, except when the CO density is reduced compared to Local Thermodynamic Equilibrium (LTE) values or the CO band can be assumed to be optically thick. The onset of radiative instabilities is expected to be strongly modified when dynamic phenomena such as stochastic shocks are present. Our results provide strong evidence that dust formation can most likely occur via a radiative instability alone. Therefore, we present a revised version of the Muchmore et al. dust formation paradigm, which also considers hydrodynamic cooling. The new paradigm is particularly relevant in cases where dust is formed relatively close to the stellar photosphere.

Cuntz, M.↗

Flow Simulation in Secondary Flow Passages of a Rocket Engine Turbopump

This paper presents application of a Generalized Fluid System Simulation Program (GFSSP) to model fluid flow in a very complex secondary flow circuit of a rocket engine turbopump, SIMPLEX. SIMPLEX was a low cost turbopump designed, manufactured and tested to demonstrate the ability to reduce the overall cost and design cycle time of turbomachinery. GFSSP is a general purpose computer program for analyzing flowrate, pressure, temperature and concentration distribution in a complex flow network. The program employs a finite volume formulation of mass, momentum and energy conservation equations in conjunction with thermodynamic equation of state of a real fluid. The secondary flow passages modeled include flow between the impeller shroud and housing, impeller back face and housing, flow through labyrinth seal and bearing, return line flow to inducer and flow adjacent to front and back face of the turbine. Axial load on the bearing are calculated from the predicted pressure on all axial faces of the rotor. The predicted pressure and temperature distributions compared well with the test data.

Schallhorn, Paul↗

High Temperature Vaporization into Different Environments

Materials in extreme environments are often negatively affected by the process of vaporization. The maximum vaporization rate is determined by calculating a flux from the equilibrium vapor pressure(s) of each component. However, this maximum vaporization flux is decreased by the vaporization coefficient and the atmosphere above the material. This paper discusses each of these factors for common oxides, using SiO2 as an example. Vaporization into a vacuum is an ideal starting point, as the basic mechanisms for vaporization have been derived from these studies. Imposing an over-pressure of a static gas changes the vaporization rates from both a thermodynamic and kinetic aspect. Finally, a flowing gas is the most common situation encountered in many extreme environments. Laminar and turbulent flow effects are be treated with both analytical transport expressions and computational fluid dynamics (CFD). The example of a vaporizing SiO2 coupon in a laboratory furnace is examined with both the analytical expressions and with CFD.

Jacobson, Nathan↗

Modeling of A Large-Scale Liquid Oxygen Propellant Transfer Test

Current planned NASA missions to the Moon and Mars involve the transfer of cryogenic propellants in orbit. Correctly predicting the flow rate during the propellant transfer operation could help with planning these operations. In lieu of expensive tests conducted on-orbit, accurate predictive computational models of these propellant transfer processes can be used to predict the propellant transfer flow rate. To gain confidence in these computational models, appropriate anchoring and validation to experimental data can be performed. Ground tests conducted at SpaceX involving the transfer of liquid oxygen were carried out to investigate the effect of the thermodynamic conditions including inlet pressure and temperature and downstream pressure on the flow rate. The results from the test were compared against flow rate predictions from a commercially available two-phase flow software. On average the model predicted the flow rates within 10% of the experimentally observed flow rates for the full range of thermodynamic inlet and outlet conditions.

Daniel Hauser↗

Experimental Validation of Nozzle Flow Simulations for Rotating Detonation Rocket Engines

Rotating detonation engines (RDEs) promise increased thermodynamic performance that may significantly enhance the capabilities of current rocket platforms. Little work has been conducted thus far to characterize the effect of nozzle design on high chamber pressure RDEs. Previous computational work was completed to understand the nozzle performance of the Purdue methane/oxygen rocket RDE. A new experimental study on a similar kerosene/oxygen RDE was conducted to validate the results of the computational study. Both a nozzleless geometry and several aerospike designs were hot-fire tested. New pressure instrumentation on the different nozzle surfaces were included to better understand the flow physics unique to RDE chamber exit conditions. Both the previous computational study and the new experimental results confirmed that for a nozzleless geometry, the RDE cycle enhances suction on the base region, an important result to determining RDE engine performance separate from nozzle effects. While the aerospike experiments showed agreement with computational results, delay of flow separation due to the RDE cycle could not be confirmed. Two aerospike geometries with different nozzle pressure ratios were experimentally evaluated. This study showed both the strength and the necessity for 3D transient computations to better understand the RDE flow field with nozzle geometries.

Alexis J Harroun↗

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗

Simulating The Performance Of A Scramjet

SCRAM computer program determines, according to one-dimensional mathematical model, cycle performance for airframe-integrated subsonic or supersonic combustion Ramjets. Reliable, efficient, and speedy ramjet designer's software tool usable on all standard computers down to those compatible with IBM PC-AT's. Scramjet cycle analysis code simulates, from nose to tail, hydrogen-fueled, airframe-integrated scramjet in real gas flow with equilibrium thermodynamic properties. Enables rapid scramjet performance estimates. Written in FORTRAN 77.

Walton, James T.↗

Strongly coupled radiative transfer and Joule heating in the cathode of an arc heater

Radiation and Joule heating in the electrode region of an arc heater are discussed. Radiative transport equations for a true axisymmetric geometry are used. A subsonic code is developed to numerically solve the fluid equations with strongly coupled radiation and Joule heating representative of a high pressure and high current arc heater. Analytic expression for the divergence of radiative heat flux derived previously is used. Jacobians of the radiation term are derived. The Joule heating term is computed using a previously developed code. The equilibrium gas model consists of seven species. The fluxes are differenced using Van Leer flux splitting. Using this code, the effects of radiative cooling on the thermodynamic parameters of the arc core are discussed.

Durgapal, P.↗

Probabilistic Study of Fluid Structure Interaction

Probabilistic CFD design is needed because we are asked to do more with less. To cost effectively accomplish the design task, we need to formally quantify the effect of uncertainties (variables) in the design. Probabilistic design is one effective method to formally quantify the effect of uncertainties. Our objective is to establish a revolutionary new early design process, by developing non-deterministic physics-based probabilistic design tools, which will include all the life cycle processes. Breakthroughs will be sought in speed, accuracy, intelligence, and usability of the system. This paper is concerned with the usefulness of parametric optimization method coupled with a Navier-Stokes analysis code for the aero-thermodynamic design of turbomachinery combustor liner. The interconnection between the CFD code and NESSUS codes facilitated the coupling between the thermal profiles and structural design. We have developed new concepts for reducing the computational cost of unsteady, three-dimensional, compressible aerodynamic analyses for multistage turbomachinery flows. The flow was modeled by the three-dimensional Favre-Reynolds-averaged Navier-Stokes equations using the k-epsilon turbulence closure, which was integrated using an implicit third-order upwind solver. The methodology developed in this paper is expected to lead to the design optimization of turbomachinery blades.

Gorla, Rama S. R.↗

Multiphysics Thermal-Fluid Analysis of a Non-Nuclear Tester for Hot-Hydrogen Materials Development

The objective of this effort is to analyze the thermal field of a non-nuclear tester, as a first step towards developing efficient and accurate multiphysics, thermo-fluid computational methodology to predict environments for hypothetical solid-core, nuclear thermal engine thrust chamber design and analysis. The computational methodology is based on a multidimensional, finite-volume, turbulent, chemically reacting, radiating, unstructured-grid, and pressure-based formulation. The multiphysics invoked in this study include hydrogen dissociation kinetics and thermodynamics, turbulent flow, convective, radiative and conjugate heat transfers.

Wang, Ten-See↗