Engineering Papers⌕ Search

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.

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

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

On Kinetics Modeling of Vibrational Energy Transfer

Two models of vibrational energy exchange are compared at equilibrium to the elementary vibrational exchange reaction for a binary mixture. The first model, non-linear in the species vibrational energies, was derived by Schwartz, Slawsky, and Herzfeld (SSH) by considering the detailed kinetics of vibrational energy levels. This model recovers the result demanded at equilibrium by the elementary reaction. The second model is more recent, and is gaining use in certain areas of computational fluid dynamics. This model, linear in the species vibrational energies, is shown not to recover the required equilibrium result. Further, this more recent model is inconsistent with its suggested rate constants in that those rate constants were inferred from measurements by using the SSH model to reduce the data. The non-linear versus linear nature of these two models can lead to significant differences in vibrational energy coupling. Use of the contemporary model may lead to significant misconceptions, especially when integrated in computer codes considering multiple energy coupling mechanisms.

Gilmore, John O.↗
NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗
NASA NTRS1987

Atmospheric reentry of the in-core thermionic SP-100 reactor system

Presumed end-of-life atmospheric reentry of the GA SP-100 system was studied to assess dispersal feasibility and associated hazards. Reentry was studied by sequential use of an orbital trajectory and a heat analysis computer program. Two heating models were used. The first model assumed a thermal equilibrium condition between the stagnation point aerodynamic heating and the radiative cooling of the skin material surface. The second model allowed for infinite conductivity of the skin material. Four reentering configurations were studied representing stages of increased SP-100 breakup: (1) radiator, shield and reactor, (2) shield and reactor, (3) reactor with control drums, and (4) reactor without control drums. Each reentering configuration was started from a circular orbit at 116 km having an inertial velocity near Mach 25. The assumed failing criterion was the attainment of melting temperature of a critical system component. The reentry analysis revealed breakup of the vessel in the neighborhood of 61 km altitude and scattering of the fuel elements. Subsequent breakup of the fuel elements was not predicted. Oxidation of the niobium skin material was calculated to cause an increase in surface temperature of less than ten percent. The concept of thermite analogs for enhancing reactor reentry dispersal was assessed and found to be feasible in principle. A conservative worst-case hazards analysis was performed for radioactive and nonradioactive toxic SP-100 materials assumed to be dispersed during end-of-life reentry. The hazards associated with this phase of the SP-100 mission were calculated to be insignificant.

Stamatelatos, M. G.↗
NASA NTRS1995

An extended structure-based model based on a stochastic eddy-axis evolution equation

We have proposed and implemented an extension of the structure-based model for weak deformations. It was shown that the extended model will correctly reduce to the form of standard k-e models for the case of equilibrium under weak mean strain. The realizability of the extended model is guaranteed by the method of its construction. The predictions of the proposed model were very good for rotating homogeneous shear flows and for irrotational axisymmetric contraction, but were seriously deficient in the case of plane strain and axisymmetric expansion. We have concluded that the problem behind these difficulties lies in the algebraic constitutive equation relating the Reynolds stresses to the structure parameters rather than in the slow model developed here. In its present form, this equation assumes that under irrotational strain the principal axes of the Reynolds stresses remain locked onto those of the eddy-axis tensor. This is correct in the RDT limit, but inappropriate under weaker mean strains, when the non-linear eddy-eddy interactions tend to misalign the two sets of principal axes and create some non-zero theta and gamma.

Kassinos, S. C.↗
NASA NTRS1989

The transient performance of a two-phase fluid reservoir

Thermal control of future large, high power spacecraft will require a two-phase fluid central bus. The two-phase fluid reservoir is a critical component in the two-phase fluid bus. It both controls the saturation temperature and provides a space for volumetric changes. A dynamic reservoir simulation model does not currently exist, but it is needed to expedite efforts and reduce risk. During 1989 an effort was made to develop a simulation model of the transient performance of a two-phase fluid reservoir. As a beginning, a preliminary model was developed. It is based upon component mathematical models in lumped parametric form and build upon five component mathematical models for calculating dynamic responses of two-phase fluid reservoirs, primary feedback elements, controller commands, heater actuators, and reservoir heaters. As much as possible, the model took advantage of the available SINDA'85/FLUINT thermal/fluid integrator. Additional calculation logic and computer subroutines were developed to complete implementation of the model. The model is capable of simulating dynamic response of an equilibrium two-phase fluid reservoir. Modification of the model to include the liquid/vapor nonequilibrium is required for applications of the model to simulate performance of reservoir in which the liquid and vapor phases of the reservoir fluid are not in equilibrium. In addition, the model in its present form, needs to be refined in several respects. More empirical data are needed to guide the model development. The model may then be used to conduct a full parametric study of two-phase fluid reservoirs. More complexities in two-phaes flow regions in laboratory and flight conditions may have to be considered eventually if empirical data cannot be simulated satisfactorily. System with other components arrangement also need to be simulated if optimization is ever to be attained. The present model does, however, preliminarily demonstrates that such analyses are quite possible and offers a far less expensive method to understand the transient of a two-phase fluid reservoir system than a totally headware approach.

Chi, Joseph↗
NASA NTRS2014

State-to-State Internal Energy Relaxation Following the Quantum-Kinetic Model in DSMC

A new model for chemical reactions, the Quantum-Kinetic (Q-K) model of Bird, has recently been introduced that does not depend on macroscopic rate equations or values of local flow field data. Subsequently, the Q-K model has been extended to include reactions involving charged species and electronic energy level transitions. Although this is a phenomenological model, it has been shown to accurately reproduce both equilibrium and non-equilibrium reaction rates. The usefulness of this model becomes clear as local flow conditions either exceed the conditions used to build previous models or when they depart from an equilibrium distribution. Presently, the applicability of the relaxation technique is investigated for the vibrational internal energy mode. The Forced Harmonic Oscillator (FHO) theory for vibrational energy level transitions is combined with the Q-K energy level transition model to accurately reproduce energy level transitions at a reduced computational cost compared to the older FHO models.

Liechty, Derek S.↗
NASA NTRS1987

Kinetic model for the vibrational energy exchange in flowing molecular gas mixtures

The present study is concerned with the development of a computational model for the description of the vibrational energy exchange in flowing gas mixtures, taking into account a given number of energy levels for each vibrational degree of freedom. It is possible to select an arbitrary number of energy levels. The presented model uses values in the range from 10 to approximately 40. The distribution of energy with respect to these levels can differ from the equilibrium distribution. The kinetic model developed can be employed for arbitrary gaseous mixtures with an arbitrary number of vibrational degrees of freedom for each type of gas. The application of the model to CO2-H2ON2-O2-He mixtures is discussed. The obtained relations can be utilized in a study of the suitability of radiation-related transitional processes, involving the CO2 molecule, for laser applications. It is found that the computational results provided by the model agree very well with experimental data obtained for a CO2 laser. Possibilities for the activation of a 16-micron and 14-micron laser are considered.

Offenhaeuser, F.↗
NASA NTRS1992

Ignition and combustion characteristics of metallized propellants

During this reporting period, theoretical work on the secondary atomization process was continued and the experimental apparatus was improved. A one-dimensional model of a rocket combustor, incorporating multiple droplet size classes, slurry combustion, secondary atomization, radiation heat transfer, and two-phase slip between slurry droplets and the gas flow was derived and a computer code was written to implement this model. The STANJAN chemical equilibrium solver was coupled with this code to yield gas temperature, density, and composition as functions of axial location. Preliminary results indicate that the model is performing correctly, given current model assumptions. Radiation heat transfer in the combustion chamber is treated as an optically-thick participating media problem requiring a solution of the radiative transfer equation. A cylindrical P sub 1 approximation was employed to yield an analytical expression for chamber-wall heat flux at each axial location. The code exercised to determine the effects of secondary atomization intensity, defined as the number of secondary drops produced per initial drop, on chamber burnout distance and final Al2O3 agglomerate diameter. These results indicate that only weak secondary atomization is required to significantly reduce these two parameters. Stronger atomization intensities were found to yield decreasing marginal benefits. The experimental apparatus was improved to reduce building vibration effects on the optical system alignment. This was accomplished by mounting the burner and the transmitting/receiving optics on a single frame supported by vibration-isolation legs. Calibration and shakedown tests indicate that vibration problems were eliminated and that the system is performing correctly.

Mueller, D. C.↗
NASA NTRS1980

Nonlinear transient analysis by energy minimization: A theoretical basis for the ACTION computer code

The formulation basis for establishing the static or dynamic equilibrium configurations of finite element models of structures which may behave in the nonlinear range are provided. With both geometric and time independent material nonlinearities included, the development is restricted to simple one and two dimensional finite elements which are regarded as being the basic elements for modeling full aircraft-like structures under crash conditions. Representations of a rigid link and an impenetrable contact plane are added to the deformation model so that any number of nodes of the finite element model may be connected by a rigid link or may contact the plane. Equilibrium configurations are derived as the stationary conditions of a potential function of the generalized nodal variables of the model. Minimization of the nonlinear potential function is achieved by using the best current variable metric update formula for use in unconstrained minimization. Powell's conjugate gradient algorithm, which offers very low storage requirements at some slight increase in the total number of calculations, is the other alternative algorithm to be used for extremely large scale problems.

Kamat, M. P.↗
NASA NTRS1987

Quasi-thermal models

Non-local thermodynamics equilibrium (LTE) effects in the photosphere; recent research on the chromosphere of the M and C stars; and elementary shock-waves and pulsation theories and their applications to Mira long-period variables are discussed.

Delareza, Ramiro↗
NASA NTRS1995

The Superorbital Expansion Tube concept, experiment and analysis

In response to the need for ground testing facilities for super orbital re-entry research, a small scale facility has been set up at the University of Queensland to demonstrate the superorbital expansion tube concept. This unique device is a free piston driven, triple diaphragm, impulse shock facility which uses the enthalpy multiplication mechanism of the unsteady expansion process and the addition of a secondary shock driver to further heat the driver gas. The pilot facility has been operated to produce quasi-steady test flows in air with shock velocities in excess of 13 km/s and with a usable test flow duration of the order of 15 micro sec. an experimental condition produced in the facility with total enthalpy of 108 MJ/kg and a total pressure of 335 MPa is reported. A simple analytical flow model which accounts for non-ideal rupture of the light tertiary diaphragm and the resulting entropy increase in the test gas is discussed. It is shown that equilibrium calculations more accurately model the unsteady expansion process than calculations assuming frozen chemistry. This is because the high enthalpy flows produced in the facility can only be achieved if the chemical energy stored in the test flow during shock heating of the test gas is partially returned to the flow during the process of unsteady expansion. Measurements of heat transfer rates to a flat plate demonstrate the usability of test flow for aerothermodynamic testing and comparison of these rates with empirical calculations confirms the usable accuracy of the flow model.

Neely, A. J.↗
NASA NTRS2013

The Metallicity Evolution of Low Mass Galaxies: New Contraints at Intermediate Redshift

We present abundance measurements from 26 emission-line-selected galaxies at z approx. 0.6-0.7. By reaching stellar masses as low as 10(exp 8) M stellar mass, these observations provide the first measurement of the intermediate-redshift mass-metallicity (MZ) relation below 10(exp 9)M stellar mass. For the portion of our sample above M is greater than 10(exp 9)M (8/26 galaxies), we find good agreement with previous measurements of the intermediate-redshift MZ relation. Compared to the local relation, we measure an evolution that corresponds to a 0.12 dex decrease in oxygen abundances at intermediate redshifts. This result confirms the trend that metallicity evolution becomes more significant toward lower stellar masses, in keeping with a downsizing scenario where low-mass galaxies evolve onto the local MZ relation at later cosmic times. We show that these galaxies follow the local fundamental metallicity relation, where objects with higher specific (mass-normalized) star formation rates (SFRs) have lower metallicities. Furthermore, we show that the galaxies in our sample lie on an extrapolation of the SFR-M∗ relation (the star-forming main sequence). Leveraging the MZ relation and star-forming main sequence (and combining our data with higher-mass measurements from the literature), we test models that assume an equilibrium between mass inflow, outflow, and star formation.We find that outflows are required to describe the data. By comparing different outflow prescriptions, we show that momentum, driven winds can describe the MZ relation; however, this model underpredicts the amount of star formation in low-mass galaxies. This disagreement may indicate that preventive feedback from gas heating has been overestimated, or it may signify a more fundamental deviation from the equilibrium assumption.

Low Mass↗
NASA NTRS1994

Implementation of a vibrationally linked chemical reaction model for DSMC

A new procedure closely linking dissociation and exchange reactions in air to the vibrational levels of the diatomic molecules has been implemented in both one- and two-dimensional versions of Direct Simulation Monte Carlo (DSMC) programs. The previous modeling of chemical reactions with DSMC was based on the continuum reaction rates for the various possible reactions. The new method is more closely related to the actual physics of dissociation and is more appropriate to the particle nature of DSMC. Two cases are presented: the relaxation to equilibrium of undissociated air initially at 10,000 K, and the axisymmetric calculation of shuttle forebody heating during reentry at 92.35 km and 7500 m/s. Although reaction rates are not used in determining the dissociations or exchange reactions, the new method produces rates which agree astonishingly well with the published rates derived from experiment. The results for gas properties and surface properties also agree well with the results produced by earlier DSMC models, equilibrium air calculations, and experiment.

Carlson, A. B.↗
NASA NTRSDate not supplied

Assessment of Mars 2020 Forebody Heating Predictions with Coupled Material Response

The objective of this work was to assess the impact of flow field - material response cou-pling on reproducing temperatures within the Mars 2020 forebody thermal protection system.Both equilibrium and nonequilibrium coupled material response approaches were considered.Turbulent transition estimates were included in the flow field modeling based on transitiontimes inferred from flight measurements. A coupled equilibrium ablation model was foundto overpredict the in-depth temperatures near the surface by over 100 K in the stagnationregion and by about 200 K near the leeside shoulder, which was subject to turbulent heatingaugmentation. For nonequilibrium ablation, two finite rate surface chemistry models wereused to model the chemical process on the vehicle surface. The first was a model consistingof N2catalysis and CO2recombination model, and the second was a model consisting of onlyN2and O2catalysis. Near the stagnation point, the N2and CO2model was found to be withinabout 10 K from the peak in-depth temperature near the surface; however, overpredictionof the in-depth temperatures was observed early in the trajectory, similar to the equilibriumablation model. The N2and O2catalysis model was found to significantly underpredict nearsurface in-depth temperatures.

Thomas K West↗
NASA NTRS1984

Origin of the Greenland Ice Sheet

A climate model resolving the seasonal cycle and the two horizontal dimensions was developed and applied to several problems of current interest. Models of this type are useful when for various reasons a general circulation model experiment is not warranted or not feasible. For example, in cases where the signal to natural variability is small it may be advantageous to first consider such a statistical dynamical model because extremely long runs may be necessary in the application. In this case the simpler statistical dynamical model serves as a pilot study device. The model developed is a thermodynamic model whose solution yields the equilibrium seasonal cycle for the surface temperature field over the globe. The model is essentially a statement of the conservation of heat energy for individual columns of the earth atmosphere system. Various terms such as the infrared radiation flux to space are parameterized with earth radiation budget data from satellites such as Nimbus 6. The primary agent modulating the seasonal cycle amplitude is the heat capacity per unit area which is a strong function of surface type -- ocean surface can store 60 times more heat per unit time than land. By adjusting its few empirical parameters the model can be brought into remarkable agreement with the observed seasonal cycle. The model is then very useful for looking at the dependence of the seasonal cycle of the temperature on such externally defined variables as the Earth's orbital elements (eccentricity, tilt, precession of equinoxes) or the configuration of land-sea geography which can be changed by continental drift.

North, G. R.↗
NASA NTRS1994

Magnetic mapping and Birkeland currents in the Toffoletto-Hill and Tsyganenko magnetosphere models

We investigate the mapping of magnetic flux tubes between the high-latitude ionosphere and the equatorial magnetosphere, using the Toffoletto and Hill (1989) theoretical model and the Tsyganenko (1987, 1989) empirical models. In all models, strong distortion of flux tubes occurs near the magnetopause and in the distant tail. These effects are illustrated for various interplanetary magnetic field directions in the Toffoletto and Hill (1989) model and for various Kp values in the Tsyganenko (1987, 1989) model. We calculate the Birkeland current that is actually present in each model (from Ampere's law), as well as the Birkeland current that would be implied if the model were in magnetostatic equilibrium with isotropic plasma pressure (from the Vasyliunas equation). Comparison of the actual and implied Birkeland currents indicates that no physically significant Birkeland current is included in any of the models tested. Each model contains a 'noise' component of Birkeland current that is, however, much smaller in magnitude than the 'real' Birkeland currents that remain to be included.

Ding, Cheng↗
NASA NTRS2018

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.

nonequilibrium modeling↗
NASA NTRS1986

A dynamic model of filament eruptions and two ribbon flares

Two basically different models for the filament equilibrium by Kippenhahn and Schluter (1957) and Kuperus and Raadu (1974) have appeared in the literature. A further analyses by van Tend and Kuperus (1978) added the force due to the horizontal component of the background field to the Kuperus and Raadu model. In order to obtain a better model which actually describes these phenomena, the evolution of the filament has to be considered in detail. A first attempt was recently presented by Kaastra. Kaastra did not formulate the precise energy balance equations for the problem, as is done in the present work. In the present model not only the force balance, but also the energy balance of the filament is taken into account. Thus a fully closed system of equations is obtained, that describes the evolution of the filament, first in force equilibrium during the current build-up phase, then in the non-equilibrium phase before the eruption, and the eruption itself. A neutral point appears above the photospheric surface in the non-equilibrium phase, but long before the eruption. It was found that although the filament itself may be in non-equilibrium, the evolution may still be slow up to the height where the eruption takes place. The eruption of the filament itself causes a large induced electric field at the neutral point which leads to the observed flare phenomena.

Kuin, N. Paul M.↗