Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibrium programming”

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

A General Method for Automatic Computation of Equilibrium Compositions and Theoretical Rocket Performance of Propellants

A general computer program for chemical equilibrium and rocket performance calculations was written for the IBM 650 computer with 2000 words of drum storage, 60 words of high-speed core storage, indexing registers, and floating point attachments. The program is capable of carrying out combustion and isentropic expansion calculations on a chemical system that may include as many as 10 different chemical elements, 30 reaction products, and 25 pressure ratios. In addition to the equilibrium composition, temperature, and pressure, the program calculates specific impulse, specific impulse in vacuum, characteristic velocity, thrust coefficient, area ratio, molecular weight, Mach number, specific heat, isentropic exponent, enthalpy, entropy, and several thermodynamic first derivatives.

Gordon, Sanford↗

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗

Computer program to solve two-dimensional shock-wave interference problems with an equilibrium chemically reacting air model

The computer program EASI, an acronym for Equilibrium Air Shock Interference, was developed to calculate the inviscid flowfield, the maximum surface pressure, and the maximum heat flux produced by six shock wave interference patterns on a 2-D, cylindrical configuration. Thermodynamic properties of the inviscid flowfield are determined using either an 11-specie, 7-reaction equilibrium chemically reacting air model or a calorically perfect air model. The inviscid flowfield is solved using the integral form of the conservation equations. Surface heating calculations at the impingement point for the equilibrium chemically reacting air model use variable transport properties and specific heat. However, for the calorically perfect air model, heating rate calculations use a constant Prandtl number. Sample calculations of the six shock wave interference patterns, a listing of the computer program, and flowcharts of the programming logic are included.

Glass, Christopher E.↗

NNEPEQ: Chemical equilibrium version of the Navy/NASA Engine Program

The Navy NASA Engine Program, NNEP, currently is in use at a large number of government agencies, commercial companies and universities. This computer code has bee used extensively to calculate the design and off-design (matched) performance of a broad range of turbine engines, ranging from subsonic turboprops to variable cycle engines for supersonic transports. Recently, there has been increased interest in applications for which NNEP was not capable of simulating, namely, high Mach applications, alternate fuels including cryogenics, and cycles such as the gas generator air-turbo-rocker (ATR). In addition, there is interest in cycles employing ejectors such as for military fighters. New engine component models had to be created for incorporation into NNEP, and it was found necessary to include chemical dissociation effects of high temperature gases. The incorporation of these extended capabilities into NNEP is discussed and some of the effects of these changes are illustrated.

Fishbach, Laurence H.↗

NNEPEQ - Chemical equilibrium version of the Navy/NASA Engine Program

The Navy NASA Engine Program, NNEP, currently is in use at a large number of government agencies, commercial companies and universities. This computer code has been used extensively to calculate the design and off-design (matched) performance of a broad range of turbine engines, ranging from subsonic turboprops to variable cycle engines for supersonic transports. Recently, there has been increased interest in applications for which NNEP was not capable of simulating, namely, high Mach applications, alternate fuels including cryogenics, and cycles such as the gas generator air-turbo-rocker (ATR). In addition, there is interest in cycles employing ejectors such as for military fighters. New engine component models had to be created for incorporation into NNEP, and it was found necessary to include chemical dissociation effects of high temperature gases. The incorporation of these extended capabilities into NNEP is discussed and some of the effects of these changes are illustrated.

Fishbach, L. H.↗

Transpiration and film cooling boundary layer computer program. Volume 1: Numerical solutions of the turbulent boundary layer equations with equilibrium chemistry

A finite difference turbulent boundary layer computer program has been developed. The program is primarily oriented towards the calculation of boundary layer performance losses in rocket engines; however, the solution is general, and has much broader applicability. The effects of transpiration and film cooling as well as the effect of equilibrium chemical reactions (currently restricted to the H2-O2 system) can be calculated. The turbulent transport terms are evaluated using the phenomenological mixing length - eddy viscosity concept. The equations of motion are solved using the Crank-Nicolson implicit finite difference technique. The analysis and computer program have been checked out by solving a series of both laminar and turbulent test cases and comparing the results to data or other solutions. These comparisons have shown that the program is capable of producing very satisfactory results for a wide range of flows. Further refinements to the analysis and program, especially as applied to film cooling solutions, would be aided by the acquisition of a firm data base.

Levine, J. N.↗

Combustion of hydrogen-air jets in local chemical equilibrium: A guide to the CHARNAL computer program

A guide to a computer program, written in FORTRAN 4, for predicting the flow properties of turbulent mixing with combustion of a circular jet of hydrogen into a co-flowing stream of air is presented. The program, which is based upon the Imperial College group's PASSA series, solves differential equations for diffusion and dissipation of turbulent kinetic energy and also of the R.M.S. fluctuation of hydrogen concentration. The effective turbulent viscosity for use in the shear stress equation is computed. Chemical equilibrium is assumed throughout the flow.

Spalding, D. B.↗

Program Computes Performances Of Scramjets

SRGULL, a scramjet-cycle-analysis computer code, engineer's software tool capable of nose-to-tail simulation of hydrogen-fueled, scramjet engine integrated with air-frame. Simulated flow that of real gas with equilibrium thermodynamic properties. Program facilities initial estimates of performance of scramjet cycle by linking code for two-dimensional forebody, inlet, and nozzle with code for one-dimensional combustor. Five computer codes (SCRAM, SEAGUL, INLET, Program HUD, and GASH) integrated in program to provide capability for analysis of changing flow conditions. Written in FORTRAN 77.

Walton, J. T.↗

Computer program for calculating and fitting thermodynamic functions

A computer program is described which (1) calculates thermodynamic functions (heat capacity, enthalpy, entropy, and free energy) for several optional forms of the partition function, (2) fits these functions to empirical equations by means of a least-squares fit, and (3) calculates, as a function of temperture, heats of formation and equilibrium constants. The program provides several methods for calculating ideal gas properties. For monatomic gases, three methods are given which differ in the technique used for truncating the partition function. For diatomic and polyatomic molecules, five methods are given which differ in the corrections to the rigid-rotator harmonic-oscillator approximation. A method for estimating thermodynamic functions for some species is also given.

Mcbride, Bonnie J.↗

A computer program for two-dimensional and axisymmetric nonreacting perfect gas and equilibrium chemically reacting laminar, transitional and-or turbulent boundary layer flows

A computer program is described in detail for laminar, transitional, and/or turbulent boundary-layer flows of non-reacting (perfect gas) and reacting gas mixtures in chemical equilibrium. An implicit finite difference scheme was developed for both two dimensional and axisymmetric flows over bodies, and in rocket nozzles and hypervelocity wind tunnel nozzles. The program, program subroutines, variables, and input and output data are described. Also included is the output from a sample calculation of fully developed turbulent, perfect gas flow over a flat plate. Input data coding forms and a FORTRAN source listing of the program are included. A method is discussed for obtaining thermodynamic and transport property data which are required to perform boundary-layer calculations for reacting gases in chemical equilibrium.

Miner, E. W.↗

NASA Glenn Coefficients for Calculating Thermodynamic Properties of Individual Species

This report documents the library of thermodynamic data used with the NASA Glenn computer program CEA (Chemical Equilibrium with Applications). This library, containing data for over 2000 solid, liquid, and gaseous chemical species for temperatures ranging from 200 to 20,000 K, is available for use with other computer codes as well. The data are expressed as least-squares coefficients to a seven-term functional form for C((sup o)(sub p)) (T) / R with integration constants for H (sup o) (T) / RT and S(sup o) (T) / R. The NASA Glenn computer program PAC (Properties and Coefficients) was used to calculate thermodynamic functions and to generate the least-squares coefficients. PAC input was taken from a variety of sources. A complete listing of the database is given along with a summary of thermodynamic properties at 0 and 298.15 K.

Bonnie J McBride↗

Interaction of pollution abatement with world dynamics

The world dynamics model of Jay W. Forrester was modified to account for pollution abatement. In the modified model, it is assumed that as pollution increases, efforts are made to control pollution. There is a competition between food supply, material standard of living, and pollution abatement for capital, and time is required for diversion of capital toward pollution abatement. Inclusion of pollution abatement in the model drastically alters the response of the world system for the case in which depletion of natural resources is not considered. Instead of undergoing a pollution catastrophe, all system levels move more or less smoothly toward an equilibrium. A FORTRAN program listing of the modified world dynamics model is included.

Smith, G. L.↗

Sodium sulfate: Vaporization thermodynamics and role in corrosive flames

Gaseous species over liquid Na2SO4 were identified by the technique of molecular beam mass spectrometry. The heat and entropy of vaporization of the Na2SO4 molecule were measured directly. Comparisons of the experimental entropy with values calculated using various molecular parameters were used to estimate the molecular structure and vibrational frequencies. The thermodynamic properties of gaseous and condensed phase Na2SO4, along with additional pertinent species, were used in a computer program to calculate equilibrium flame compositions and temperatures for representative turbine engine and burner rig flames. Compositions were calculated at various fuel-to-oxidant ratios with additions of sulfur to the fuel and the components of sea salt to the intake air. Temperatures for condensation of Na2SO4 were obtained as a function of sulfur and sea salt concentrations.

Kohl, F. J.↗

Sodium sulfate - Vaporization thermodynamics and role in corrosive flames

Mass spectrometer experiments were conducted to determine the thermodynamic properties of gaseous Na2SO4, and these data were used in a computer program to calculate equilibrium flame compositions and temperatures for representative turbine engine and burner rig flames. The work is important in that sodium sulfate is the major phase recovered from turbine surfaces after instances of corrosion, due to the presence of sulfur in fuels and sodium chloride in intake air.

Kohl, F. J.↗