Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equilibrium calculations”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Analytical study of striated nozzle flow with small radius of curvature ratio throats

An analytical method was developed which is capable of estimating the chamber and throat conditions in a nozzle with a low radius of curvature throat. The method was programmed using standard FORTRAN 4 language and includes chemical equilibrium calculation subprograms (modified NASA Lewis program CEC71) as an integral part. The method determines detailed and gross rocket characteristics in the presence of striated flows and gives detailed results for the motor chamber and throat plane with as many as 20 discrete zones. The method employs a simultaneous solution of the mass, momentum, and energy equations and allows propellant types, 0/F ratios, propellant distribution, nozzle geometry, and injection schemes to be varied so to predict spatial velocity, density, pressure, and other thermodynamic variable distributions in the chamber as well as the throat. Results for small radius of curvature have shown good comparison to experimental results. Both gaseous and liquid injection may be considered with frozen or equilibrium flow calculations.

Norton, D. J.↗
NASA NTRS1993

Shock-induced devolatization of calcium sulfate and implications for K-T extinctions

Calcium sulfate devolatization during the impact at Chicxulub, Mexico and dispersal in the stratosphere of the resultant sulfuric acid aerosol have been suggested as a possible mechanism for the Cretaceous-Tertiary extinctions. In this paper, we investigated two shock-induced devolatization reactions of calcium sulfate up to 42 GPa in the laboratory: CaSO4 + SiO2 yields CaSiO3 + SO3(degassed) and CaSO4 yields CaO + SO2(degassed) + 1/2 O2(degassed). We found both to proceed to a much less extent than calculated by equilibrium thermodynamic calculations. Reaction products are found to be 10(exp -2) times those calculated for equilibrium. Consequently our estimate of the amount of sulfur oxides degassed into the atmosphere from shock devolatization of CaS04 in the Chicxulub lithographic section (6x10(exp 15)-2x10(exp 16)g in sulfur mass) is lower by a factor of 70 to 400 than previous estimates; the related environmental stress arising from the resultant global cooling of approximately 4 K and fallout of acid rain does not appear to suffice to explain the widespread K-T extinctions.

Chen, Guangqing↗
NASA NTRS1994

Shock-induced devolatilization of calcium sulfate and implications for K-T extinctions

The devolatilization of calcium sulfate, which is present in the target rock of the Chicxulub, Mexico impact structure, and dispersal in the stratosphere of the resultant sulfuric acid aerosol have been suggested as a possible mechanism for the Cretaceous-Tertiary extinctions. We measured the amount of SO2 produced from two shock-induced devolatilization reactions of calcium sulfate up to 42 GPa in the laboratory. We found both to proceed to a much lower extent than calculated by equilibrium thermodynamic calculations. Reaction products are found to be approx. 10(exp -2) times those calculated for equilibrium. Upon modeling the quantity of sulfur oxides degassed into the atmosphere from shock devolatilization of CaSO4 in the Chicxulub lithographic section, the resulting 9 x 10(exp 16) to 6 x 10(exp 17) g (in sulfur mass) is lower by a factor of 10-100 than previous upper limit estimates, the related environmental stress arising from the resultant global cooling and fallout of acid rain is insufficient to explain the widespread K-T extinctions.

Chen, Guangqing↗
NASA NTRS1980

Measurements of quiescent prominence spectra - Ca/+/ double-lambda 8542, 8498 and He lambda 4471, Ti/+/ lambda 4468

High-dispersion measurements of the line profiles of the Ca(+) lines at 8542 and 8498 A and the He 4471-A, Ti(+) 4468-A line pair in the spectra of quiescent solar prominences are presented. The measured Ca(+) and He line profiles are found to be accurately described by the two-component model of Landman et al. (1977) and Landman and Mongillo (1979) in which the prominence along the line of sight consists of a low-temperature region surrounded by high-temperature material representing contributions from the prominence/corona interface; the Ti(+) lines are too weak to distinguish among different models. The Ca(+) data allow a constant infrared line source function for the low-temperature region of 50,000 ergs/sq cm per sec per ster per A to be calculated, which is in disagreement with the previous Ca(+) excitation equilibrium calculations but is consistent with generally accepted values of the electron density and temperature. The integrated intensity ratio of 0.048 between Ti(+) and He emission indicates a low-temperature region temperature of 6620 K and microturbulence of 5.1 km/sec and high-temperature microturbulence of 11 km/sec for Ca(+) and 20 km/sec for He. The results exclude macroturbulence as a possible explanation of the generally observed wing enhancements.

Landman, D. A.↗
NASA NTRS1989

Determining Equilibrium Position For Acoustical Levitation

Equilibrium position and orientation of acoustically-levitated weightless object determined by calibration technique on Earth. From calibration data, possible to calculate equilibrium position and orientation in presence of Earth gravitation. Sample not levitated acoustically during calibration. Technique relies on Boltzmann-Ehrenfest adiabatic-invariance principle. One converts resonant-frequency-shift data into data on normalized acoustical potential energy. Minimum of energy occurs at equilibrium point. From gradients of acoustical potential energy, one calculates acoustical restoring force or torque on objects as function of deviation from equilibrium position or orientation.

Barmatz, M. B.↗
NASA NTRS1970

Thermodynamics of the internal combustion engine

The ubiquitous automobile has played both the role of hero and villain during the course of its history. At the present time it is being cast as the villain with increasing frequency because it has been indicted as a major source of air pollution. It is just this characterization that is responsible for our renewed interest in the automobile. When we speak of the automobile as a source of air pollution, we really mean the heart of the automobile, the internal combustion engine. The prototype of the modern internal combustion engine was first successfully operated by Nicholas August Otto in 1876. In spite of its long history, the internal combustion engine is still poorly understood in terms of the details of the physical processes that take place during its operations. The internal combustion engine operates on what is known as the Otto cycle, and Otto cycle calculations have been with us for many years. The emphasis in earlier calculations was usually on the performance rather than on the pollution aspects. Performance characteristics can often be estimated reasonably accurately even with a considerable degree of simplification. We are interested in the pollution aspects of the internal combustion engine, and therefore the calculations must be performed with more care. Our extensive experience with equilibrium calculations places us in a position to carry out such careful calculations. We have just recently begun such an analysis, and its objective is the prediction of exhaust gas properties. I will describe our progress by first reviewing the Otto cycle, and then I will illustrate the effect of various assumptions on the calculated results.

Frank J. Zeleznik↗
NASA NTRS2011

The lagRST Model: A Turbulence Model for Non-Equilibrium Flows

This study presents a new class of turbulence model designed for wall bounded, high Reynolds number flows with separation. The model addresses deficiencies seen in the modeling of nonequilibrium turbulent flows. These flows generally have variable adverse pressure gradients which cause the turbulent quantities to react at a finite rate to changes in the mean flow quantities. This "lag" in the response of the turbulent quantities can t be modeled by most standard turbulence models, which are designed to model equilibrium turbulent boundary layers. The model presented uses a standard 2-equation model as the baseline for turbulent equilibrium calculations, but adds transport equations to account directly for non-equilibrium effects in the Reynolds Stress Tensor (RST) that are seen in large pressure gradients involving shock waves and separation. Comparisons are made to several standard turbulence modeling validation cases, including an incompressible boundary layer (both neutral and adverse pressure gradients), an incompressible mixing layer and a transonic bump flow. In addition, a hypersonic Shock Wave Turbulent Boundary Layer Interaction with separation is assessed along with a transonic capsule flow. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWTBLI flows assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗
NASA NTRS1994

The metallicity of M4: Accurate spectroscopic fundamental parameters for four giants

High-quality spectra, covering the wavelength range 5480 to 7080 A, have been obtained for four giant stars in the intermediate-metallicity CN-bimodal globular cluster M4 (NGC 6121). We have employed a model atmosphere analysis that is entirely independent from cluster parameters, such as distance, age, and reddening, in order to derive accurate values for the stellar parameters effective temperature, surface gravity, and microturbulence, and for the abundance of iron relative to the Sun, (Fe/H), and of calcium, Ca/H, for each of the four stars. Detailed radiative transfer and statistical equilibrium calculations carried out for iron and calcium suggest that departures from local thermodynamic equilibrium are not significant for the purposes of our analysis. The spectroscopically derived effective temperatures for our program stars are hotter by about 200 K than existing photometric calibrations suggest. We conclude that this is due partly to the uncertain reddening of M4 and to the existing photometric temperature calibration for red giants being too cool by about 100 K. Comparison of our spectroscopic and existing photometric temperatures supports the prognosis of a significant east-west gradient in the reddening across M4. Our derived iron abundances are slightly higher than previous high-resolution studies suggested; the differences are most probably due to the different temperature scale and choice of microturbulent velocities adopted by earlier workers. The resulting value for the metallicity of M4 is (Fe/H )(sub M4) = -1.05 + or - 0.15. Based on this result, we suggest that metallicities derived in previous high-dispersion globular cluster abundance analyses could be too low by 0.2 to 0.3 dex. Our calcium abundances suggest an enhancement of calcium, an alpha element, over iron, relative to the Sun, in M4 of (Ca/H) = 0.23.

Drake, J. J.↗
NASA NTRS1975

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.↗
NASA NTRS1986

Temperature fluctuations and infrared emission from dust particles in a hot gas

In this paper, the effect of temperature fluctuations in very small dust particles caused by electronic collisions with an ambient hot gas, is calculated. The dust-temperature distribution differs strongly from that derived on the basis of equilibrium heating models, peaking at lower temperatures in the stochastically heated case. The resulting infrared spectrum is broadened at short wavelengths, giving much more emission than expected from equilibrium calculations, and at long wavelengths exhibits a lower color temperature. Dust-temperature fluctuations are thus expected to play an important role in determining the infrared spectrum of dusty, X-ray-emitting plasmas. A power-law distribution of grain radii extended to very small grain sizes will have a distinct infrared spectrum characterized by an infrared excess at the Wien side of the spectrum. This excess emission represents a high-temperature component in the spectrum, resulting from the stochastic heating of very small dust particles. Its magnitude depends on the grain-size distribution and the temperature and density of the ambient hot gas. This excess emission may be observed by IRAS in supernova remnants and rich clusters of galaxies.

Dwek, E.↗
NASA NTRS1978

Grain formation in astronomical systems: A critical review of condensation processes

An analysis is presented of the assumption and the applicability of the three theoretical methods for calculating condensations in cosmic clouds where no pre-existing nuclei exist. The three procedures are: thermodynamic equilibrium calculations, nucleation theory, and a kinetic treatment which would take into account the characteristics of each individual collision. Thermodynamics provide detailed results on the composition temperature and composition of the condensate provided the system attains equilibrium. Because of the cosmic abundance mixture of elements, large supersaturations in some cases and low pressures, equilibrium is not expected in astronomical clouds. Nucleation theory, a combination of thermodynamics and kinetics, has the limitations of each scheme. Kinetics, not requiring equilibrium, avoids nearly all the thermodynamics difficulties but requires detailed knowledge of many reactions which thermodynamics avoids. It appears to be the only valid way to treat grain formation in space. A review of experimental studies is given.

Donn, B.↗
DOE OSTI2025

Reduced Model of Ionization Lag in Intense Laser-Produced Plasmas

A physics-based empirical formula is derived to predict the ionization lag in underdense plasmas generated by an intense laser. Time-dependent nonlocal thermodynamic equilibrium calculations demonstrate significantly delayed ionization, due to rapid changes in plasma conditions, which critically impacts plasma properties such as opacity, emissivity, and heat transport. The reduced model, based on these calculations, enables the estimation of ionization lag without requiring in-depth knowledge of nonlocal thermodynamic equilibrium modeling. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by photoionization—plays a crucial role in this ionization delay, with collisional excitation setting the timescale for ionization. Simulations across a range of elements, from beryllium to germanium, demonstrate that ionization lag is a widespread phenomenon, underscoring the importance of incorporating such processes into ionization modeling in radiation hydrodynamic simulations for various laser-plasma experiments.

direct drive↗
NASA NTRS1991

Computerized reduction of elementary reaction sets for combustion modeling

If the entire set of elementary reactions is to be solved in the modeling of chemistry in computational fluid dynamics, a set of stiff ordinary differential equations must be integrated. Some of the reactions take place at very high rates, requiring short time steps, while others take place more slowly and make little progress in the short time step integration. The goal is to develop a procedure to automatically obtain sets of finite rate equations, consistent with a partial equilibrium assumptions, from an elementary set appropriate to local conditions. The possibility of computerized reaction reduction was demonstrated. However, the ability to use the reduced reaction set depends on the ability of the CFD approach in incorporate partial equilibrium calculations into the computer code. Therefore, the results should be tested on a code with partial equilibrium capability.

Wikstrom, Carl V.↗
NASA NTRS1994

Computer program for calculation of complex chemical equilibrium compositions and applications. Part 1: Analysis

This report presents the latest in a number of versions of chemical equilibrium and applications programs developed at the NASA Lewis Research Center over more than 40 years. These programs have changed over the years to include additional features and improved calculation techniques and to take advantage of constantly improving computer capabilities. The minimization-of-free-energy approach to chemical equilibrium calculations has been used in all versions of the program since 1967. The two principal purposes of this report are presented in two parts. The first purpose, which is accomplished here in part 1, is to present in detail a number of topics of general interest in complex equilibrium calculations. These topics include mathematical analyses and techniques for obtaining chemical equilibrium; formulas for obtaining thermodynamic and transport mixture properties and thermodynamic derivatives; criteria for inclusion of condensed phases; calculations at a triple point; inclusion of ionized species; and various applications, such as constant-pressure or constant-volume combustion, rocket performance based on either a finite- or infinite-chamber-area model, shock wave calculations, and Chapman-Jouguet detonations. The second purpose of this report, to facilitate the use of the computer code, is accomplished in part 2, entitled 'Users Manual and Program Description'. Various aspects of the computer code are discussed, and a number of examples are given to illustrate its versatility.

Gordon, Sanford↗
DOE OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
NASA NTRS1986

Cosmogenic neutron-capture-produced nuclides in stony meteorites

The complete neutron-flux results and production rates for Cl-36, Ni-59, and Co-60 in stony meteorites of various radii and composition are presented. The relative neutron source strengths and neutron production-versus-depth profiles were determined by using calculated H-3 production rates. The absolute source strengths were normalized to that determined for the moon by Woolum et al. (1975). The energy spectrum of the source neutrons and the neutron transport calculations, which employed the ANISN computer code, were similar to those used for the moon by Lingenfelter et al. (1972). The production rates of the three radionuclides were determined as a function of depth in various spherical meteoroids from the calculated equilibrium neutron-flux distributions and from energy-dependent neutron-capture cross sections. Rates for producing these radionuclides by spallation reactions were also calculated.

Spergel, M. S.↗
NASA NTRS1982

Solar transition region response to variations in the heating rate

The response of a numerical model for the upper chromosphere, transition region, and corona to variations in the energy input has been examined. The numerical model solves the set of one-dimensional two-fluid hydrodynamic equations in a simple vertical magnetic flux tube. The atmosphere responds to both the increase and decrease in energy deposition by smoothly readjusting the temperature gradient and the amount of material in the region of peak radiating efficiency to radiate away energy being deposited. At no time during this readjustment is a departure from a thin laminar transition region structure seen. In addition, a time-dependent description of the nonequilibrium ionization of all of the ionization stages of oxygen has been included. This calculation is coupled with the self-consistent calculations of the dynamical variables. It is found that the nonequilibrium ionization balance calculations for both heating and cooling small loops in the quiet sun predict relative ionic abundances which differ substantially from those which would be predicted by an equilibrium calculation

Mariska, J. T.↗