SEARCH · Engineering Papers

Results for “equilibrium modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Monte Carlo computer simulations of Venus equilibrium and global resurfacing models

Two models have been proposed for the resurfacing history of Venus: (1) equilibrium resurfacing and (2) global resurfacing. The equilibrium model consists of two cases: in case 1, areas less than or equal to 0.03 percent of the planet are spatially randomly resurfaced at intervals of less than or greater than 150,000 yr to produce the observed spatially random distribution of impact craters and average surface age of about 500 m.y.; and in case 2, areas greater than or equal to 10 percent of the planet are resurfaced at intervals of greater than or equal to 50 m.y. The global resurfacing model proposes that the entire planet was resurfaced about 500 m.y. ago, destroying the preexisting crater population and followed by significantly reduced volcanism and tectonism. The present crater population has accumulated since then with only 4 percent of the observed craters having been embayed by more recent lavas. To test the equilibrium resurfacing model we have run several Monte Carlo computer simulations for the two proposed cases. It is shown that the equilibrium resurfacing model is not a valid model for an explanation of the observed crater population characteristics or Venus' resurfacing history. The global resurfacing model is the most likely explanation for the characteristics of Venus' cratering record. The amount of resurfacing since that event, some 500 m.y. ago, can be estimated by a different type of Monte Carolo simulation. To date, our initial simulation has only considered the easiest case to implement. In this case, the volcanic events are randomly distributed across the entire planet and, therefore, contrary to observation, the flooded craters are also randomly distributed across the planet.

Dawson, D. D.↗

NASA NTRS1999

Model Atmospheres for Novae in Outburst: Summary of Research

This paper presents a final report and summary of research on Model Atmospheres for Novae in Outburst. Some of the topics include: 1) Detailed NLTE (non-local thermodynamic equilibrium) Model Atmospheres for Novae during Outburst: II. Modeling optical and ultraviolet observations of Nova LMC 1988 #1; 2) A Non-LTE Line-Blanketed Stellar Atmosphere Model of the Early B Giant epsilon CMa; 3) Spectroscopy of Low Metallicity Stellar atmospheres; 4) Infrared Colors at the Stellar/Substellar Boundary; 5) On the abundance of Lithium in T CrB; 6) Numerical Solution of the Expanding Stellar Atmosphere Problem; and 7) The NextGen Model Atmosphere grid for 3000 less than or equal to T (sub eff) less than or equal to 10000K.

Hauschildt, Peter H.↗

NASA NTRS2006

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters (or closwre models) for an initial thermai non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two CFD codes currently being used at Glenn Research Center (GRC) for Stirling engine modeling are Fluent and CFD-ACE. The porous-media models available in each of these codes are equilibrium models, which assmne that the solid matrix and the fluid are in thermal equilibrium at each spatial location within the porous medium. This is believed to be a poor assumption for the oscillating-flow environment within Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, we non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location end time during the cycle. A NASA regenerator research grant has been providing experimental and computational results to support definition of various empirical coefficients needed in defining a noa-equilibrium, macroscopic, porous-media model (i.e., to define "closure" relations). The grant effort is being led by Cleveland State University, with subcontractor assistance from the University of Minnesota, Gedeon Associates, and Sunpower, Inc. Friction-factor and heat-transfer correlations based on data taken with the NASAlSunpower oscillating-flow test rig also provide experimentally based correlations that are useful in defining parameters for the porous-media model; these correlations are documented in Gedeon Associates' Sage Stirling-Code Manuals. These sources of experimentally based information were used to define the following terms and parameters needed in the non-equilibrium porous-media model: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity (including themal dispersion and estimate of tortuosity effects}, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity (including the effect of tortuosity) was also estimated. Determination of the porous-media model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Convertor (TDC), which uses a random-fiber regenerator matrix. The non-equilibrium porous-media model presented is considered to be an initial, or "draft," model for possible incorporation in commercial CFD codes, with the expectation that the empirical parameters will likely need to be updated once resulting Stirling CFD model regenerator and engine results have been analyzed. The emphasis of the paper is on use of available data to define empirical parameters (and closure models) needed in a thermal non-equilibrium porous-media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates. However, it is anticipated that a thermal non-equilibrium model such as that presented here, when iacorporated in the CFD codes, will improve our ability to accurately model Stirling regenerators with CFD relative to current thermal-equilibrium porous-media models.

Tew, Roy↗

NASA NTRS1986

Diagnostic analysis of two-dimensional monthly average ozone balance with Chapman chemistry

Chapman chemistry has been used in a two-dimensional model to simulate ozone balance phenomenology. The similarity between regions of ozone production and loss calculated using Chapman chemistry and those computed using LIMS and SAMS data with a photochemical equilibrium model indicate that such simplified chemistry is useful in studying gross features in stratospheric ozone balance. Net ozone production or loss rates are brought about by departures from the photochemical equilibrium (PCE) condition. If transport drives ozone above its PCE condition, then photochemical loss dominates production. If transport drives ozone below its PCE condition, then photochemical production dominates loss. Gross features of ozone loss/production (L/P) inferred for the real atmosphere from data are also simulated using only eddy diffusion. This indicates that one must be careful in assigning a transport scheme for a two-dimensional model that mimics only behavior of the observed ozone L/P.

Stolarski, Richard S.↗

NASA NTRS2017

On the Development of a New Nonequilibrium Chemistry Model for Mars Entry

This paper represents a summary of results to date of an on-going effort at NASA Ames Research Center to develop a physics-based non-equilibrium model for hypersonic entry into the Martian atmosphere. Our approach is to first compute potential energy surfaces based on accurate solutions of the electronic Schroedinger equation and then use quasiclassical trajectory calculations to obtain reaction cross sections and rate coefficients based on these potentials. We have presented new rate coefficients for N2 dissociation and CO dissociation and exchange reactions. These results illustrate shortcomings with some of the rate coefficients in Parks original T-Tv model for Mars entries and with some of the 30-45 year old shock tube data. We observe that the shock tube experiments of CO + O dissociation did not adequately account for the exchange reaction that leads to formation of C + O2. This reaction is actually the primary channel for CO removal in the shock layer at temperatures below 10,000 K, because the reaction enthalpy for exchange is considerably lower than the comparable value for dissociation.

quasiclassical trajectory calculation↗

NASA NTRS1993

The embedded young stars in the Taurus-Auriga molecular cloud. II - Models for scattered light images

We describe NIR imaging observations of embedded young stars in the Taurus-Auriga molecular cloud. We find a large range in J-K and H-K colors for these class I sources. The bluest objects have colors similar to the reddest T Tauri stars in the cloud; redder objects lie slightly above the reddening line for standard ISM dust and have apparent K extinctions of up to 5 mag. Most of these sources also show extended NIR emission on scales of 10-20 arcsec which corresponds to linear sizes of 1500-3000 AU. The NIR colors and nebular morphologies for this sample and the magnitude of linear polarization in several sources suggest scattered light produces most of the NIR emission in these objects. We present modeling results that suggest mass infall rates that agree with predictions for cold clouds and are generally consistent with rates estimated from radiative equilibrium models. For reasonable dust grain parameters, the range of colors and extinctions require flattened density distributions with polar cavities evacuated by bipolar outflows. These results support the idea that infall and outflow occur simultaneously in deeply embedded bipolar outflow sources. The data also indicate fairly large centrifugal radii and large inclinations to the rotational axis for a typical source.

Kenyon, Scott J.↗

NASA NTRS1994

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗

NASA NTRS1990

Model of Jovian F region ionosphere (Jovian ionosphere model in offset dipole approximation)

The geomagnetic control of the Earth's atmosphere is well understood. In the F-region and the topside ionosphere, non-electrical forces transport plasma along the magnetic field lines only. In consequence, the worldwide distribution of ionization is strongly dependent on the dip angle. For example, the equatorial anomaly is roughly symmetrical about the dipole equator rather than the geographic. The same appears to be the case in the Jovian ionosphere (Mahajan, 1981). The influence of the magnetic field of Jupiter on its ionization pattern is one of several outstanding topics which need to be studied. Tan (1986) investigated the formation of the equatorial anomaly in the Jovian ionosphere under a centered dipole model. Tan (1988) further studied the effect of the tilt of the Jovian dipole. The results were in broad agreement with those of a diffusive equilibrium model (Tan and Wu, 1981). An off-centered dipole model is constructed and its effects on the ionization pattern are investigated.

Tan, A.↗

NASA NTRS2006

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters for an initial thermal non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two codes currently used at Glenn Research Center for Stirling modeling are Fluent and CFD-ACE. The codes porous-media models are equilibrium models, which assume solid matrix and fluid are in thermal equilibrium. This is believed to be a poor assumption for Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, use non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location and time during the cycle. Experimentally based information was used to define: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity was also estimated. Determination of model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Converter (TDC), which uses a random-fiber regenerator matrix. Emphasis is on use of available data to define empirical parameters needed in a thermal non-equilibrium porous media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates.

Tew, Roy C.↗

NASA NTRS1980

Linear analysis of poststall gyrations

The poststall gyrations of a high-performance aircraft are investigated by approximating the equations of motion with a linear mathematical model. Equilibrium poststall spin conditions are determined, which are not limited to small angles or spin rates and are used as reference conditions about which the linear model is developed. As a result, complete coupling between lateral and longitudinal motions is retained. Three aerodynamic models are evaluated using this linear analysis. These models include pure rotary derivatives obtained from VPI&SU stability tunnel, forced-oscillation data obtained from NASA tests, and a combination of the two which separates the unsteady and pure rotary derivatives. Results indicate that this last model gives better predictions of aircraft motion during poststall gyrations at high angles of attack and large sideslip angles. Comparisons are made with full-scale flight test results.

Hreha, M. A.↗

NASA NTRS1994

A model for studying the composition and chemical effects of stratospheric aerosols

We developed polynomial expressions for the temperature dependence of the mean binary and water activity coefficients for H2SO4 and HNO3 solutions. These activities were used in an equilibrium model to predict the composition of stratospheric aerosols under a wide range of environmental conditions. For typical concentrations of H2O, H2SO4, HNO3, HCl, HBr, HF, and HOCl in the lower stratosphere, the aerosol composition is estimated as a function of the local temperature and the ambient relative humidity. For temperatures below 200 K, our results indicate that (1) HNO3 contributes a significant mass fraction to stratospheric aerosols, and (2) HCl solubility is considerably affected by HNO3 dissolution into sulfate aerosols. We also show that, in volcanically disturbed periods, changes in stratospheric aerosol composition can significantly alter the microphysics that leads to the formation of polar stratospheric clouds. The effects caused by HNO3 dissolution on the physical and chemical properties of stratospheric aerosols are discussed.

Tabazadeh, Azadeh↗

NASA NTRS1993

A comparison of the predictive capabilities of several turbulence models using upwind and central-difference computer codes

Four turbulence models are described and evaluated for transonic flows using the upwind code CFL3D and the central-difference code TLNS3D. In particular, the effects of recent modifications to the half-equation model of Johnson-King are explored in detail, and different versions of the model are compared. This model can obtain good results for both two-dimensional (2D) and three-dimensional (3D) separated flows. The one-equation models of Baldwin-Barth and Spalart-Allmaras perform well for separated airfoil flows, but can predict the shock too far forward at the outboard stations of a separated wing. The equilibrium model of Baldwin-Lomax predicts the shock location too far aft for both 2D and 3D separated flows, as expected. In general, all models perform well for attached or mildly separated flows.

Rumsey, Christopher L.↗

NASA NTRS1993

A Comparison of the Predictive Capabilities of Several Turbulence Models using Upwind and Central-Difference Computer Codes

Rumsey, Christopher L.↗

NASA NTRS1981

Models of transition region and coronal plasma in solar loop structures

The theory of coronal loops was developed and the simple version of the theory and observations are compared. It is shown that the gross properties of the confined corona can be understood by considering simple hydrostatic equilibrium models. By observing the corona it is apparent that hydrostatic models are seriously inadequate. It is suggested that observations of intensity fluctuations, of persistent up and down flows, of relatively cool matter residing at coronal heights, of apparent spatial co-mingling of hot plasmas at quite different temperatures, need a more sophisticated modeling.

Raymond, J. C.↗

NASA NTRS1978

Stellar model chromospheres. VII - Capella /G5 III +/, Pollux /K0 III/, and Aldebaran /K5 III/

Data from high-resolution SEC vidicon spectroscopy with a ground-based telescope (for the Ca II K line) and from spectral scans made with the BUSS ultraviolet balloon spectrograph (for the Mg II h and k lines) are used to derive models of the chromospheres and upper photospheres of three G-K giants. The models are based on partial-redistribution analyses of the Ca II K line wings and cores and on the fluxes in the Mg II lines. The photospheres thus computed are hotter than predicted by radiative-equilibrium models. The minimum-to-effective temperature ratio is found to decrease with decreasing effective temperature, while the mass column density at the top of the chromosphere increases with decreasing stellar surface gravity. The computed pressure at the chromosphere top in the primary member of the Capella spectroscopic binary system is 70 times smaller than the transition-region pressure derived by Haisch and Linsky (1976), which suggests that additional terms must be included in the transition-region energy equations for giant stars. Estimates of the Ca II and hydrogen column densities are made for the circumstellar envelope of Aldebaran.

Kelch, W. L.↗

NASA NTRS1991

Production of organic compounds in plasmas: A comparison among electric sparks, laser-induced plasmas and UV light

In order to study the production of organic compounds in plasmas (and shocks), various mixtures of N2, CH4, and H2, modeling the atmosphere of Titan, were exposed to discrete sparks, laser-induced plasmas (LIP) and ultraviolet light. The yields of HCN and simple hydrocarbons were measured and compared to those calculated from a simple quenched thermodynamic equilibrium model. The agreement between experiment and theory was fair for HCN and C2H2. However, the yields of C2H6 and other hydrocarbons were much higher than those predicted by the model. Our experiments suggest that photolysis by ultraviolet light from the plasma is an important process in the synthesis. This was confirmed by the photolysis of gas samples exposed to the light, but not to the plasma or shock waves. The results of these experiments demonstrate that, in addition to the well-known efficient synthesis of organic compounds in plasmas, the yields of saturated species, e.g., ethane, may be higher than predicted by theory and that LIP provide a convenient and clean way of simulating planetary lightning and impact plasmas in the laboratory.

Scattergood, T. W.↗

NASA NTRS1981

On Hybrid and mixed finite element methods

Three versions of the assumed stress hybrid model in finite element methods and the corresponding variational principles for the formulation are presented. Examples of rank deficiency for stiffness matrices by the hybrid stress model are given and their corresponding kinematic deformation modes are identified. A discussion of the derivation of general semi-Loof elements for plates and shells by the hybrid stress method is given. It is shown that the equilibrium model by Fraeijs de Veubeke can be derived by the approach of the hybrid stress model as a special case of semi-Loof elements.

Pian, T. H. H.↗

NASA NTRS1997

On the Dramatic Spin-up/Spin-Down Torque Reversals in Accreting Pulsars

Dramatic torque reversals between spin-up and spin-down have been observed in half of the persistent X-ray pulsars monitored by the Burst and Transient Space Experiment (BATSE) all-sky monitor on the Compton Gamma Ray Observatory. Theoretical models developed to explain early pulsar timing data can explain spin-down torques via a disk-magnetosphere interaction if the star nearly corotates with the inner accretion disk. To produce the observed BATSE torque reversals, however, these equilibrium models require the disk to alternate between two mass accretion rates, with M+/- producing accretion torques of similar magnitude but always of opposite sign. Moreover, in at least one pulsar (GX 1+4) undergoing secular spin-down, the neutron star spins down faster during brief (approximately 20 day) hard X-ray flares-this is opposite the correlation expected from standard theory, assuming that BATSE pulsed flux increases with mass accretion rate. The 10 day to 10 yr intervals between torque reversals in these systems are much longer than any characteristic magnetic or viscous timescale near the inner disk boundary and are more suggestive of a global disk phenomenon. We discuss possible explanations of the observed torque behavior. Despite the preferred sense of rotation defined by the binary orbit, the BATSE observations are surprisingly consistent with an earlier suggestion for GX 1+4: the disks in these systems somehow alternate between episodes of prograde and retrograde rotation. We are unaware of any mechanism that could produce a stable retrograde disk in a binary undergoing Roche lobe overflow, but such flip-flop behavior does occur in numerical simulations of wind-fed systems. One possibility is that the disks in some of these binaries are fed by an X-ray-excited wind.

Nelson, Robert W.↗