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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibruim states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, Charles G.↗
NASA NTRS1992

Theoretical studies of the outer envelopes of young stellar objects

With the Monte Carlo code developed by Whitney and Hartmann, a series of models was computed of scattering in disks around young stellar objects. The code calculates scattering by dust, including polarization, in arbitrary geometries. By computing model images, it was found that disk, by themselves, around young stellar objects would be very difficult to detect with present day imaging techniques. In comparing these images to observations of young stellar objects which show diffuse structure, little resemblance was found. A flared disk system will only give high polarization when viewed edge-on, and the position angle is always oriented perpendicular to the disk plane. This suggests that an envelope, perhaps the remnant infalling envelope, must be present to scatter more stellar light than a disk can, and obscure the star at many inclinations. A grid was computed of models of scattering in a disk+envelope system. Evidence is presented that the wind of the pre-main sequence object FU Orionis arises from the surface of the luminous prostellar accretion disk. A disk wind model calculated assuming radiative equilibrium explains the differential behavior of the observed asymmetrical absorption line profiles. The model predicts that strong lines should be asymmetric and blueshifted, while weak lines should be symmetric and doubled peaked due to disk rotation, in agreement with observations.

Hartmann, Lee↗
NASA NTRS2018

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

This work presents a general framework for model reduction of non-equilibrium energy transfer and dissociation processes. The multi-group maximum-entropy method is coupled with the quasi-classical trajectory method to directly construct a reduced order model for chemical non-equilibrium. Kinetic data is calculated by applying detailed balance at a microscopic level, overcoming the limitation of calculating recombination kinetic data. This approach enables the construction of a reduced order model for kinetics which bypasses the need to compute state-to-state kinetic data. This physics based reduced order model ensures that an equilibrium distribution is reached given infinite time, while allowing for non-equilibrium distributions during the relaxation and dissociation processes. A proof-of-concept test case demonstrates the applicability of this model by comparison with state-to-state kinetic data for the N2(Chi (sup 1) Sigma ((sub g), (sup plus sign)) minus N ((exp 4) S (sub u)) system with excellent agreement.

Maximum Entropy↗
NASA NTRS1988

Phase equilibrium constraints on the howardite-eucrite-diogenite association

Model determinations of fractional crystallization and equilibrium partial melting in the 0-10-kbar range have been performed for a series of compositions relevant to diogenite and eucrite petrogenesis. Olivine is found to react with diogenite parent liquids along the plagioclase-absent olivine/low-Ca pyroxene liquidus boundary under conditions of both fractional and equilibrium crystallization up to about 2 kbar. Olivine also reacts with eucritic liquids saturated with plagioclase and low-Ca pyroxene to pressures in excess of 2 kbar. The ability of simple fractional crystallization at 2 kbar to account for mineralogical and chemical features of the diogenite-eucrite association is discussed.

Longhi, John↗
NASA NTRS1993

Combustion chamber analysis code

A three-dimensional, time dependent, Favre averaged, finite volume Navier-Stokes code has been developed to model compressible and incompressible flows (with and without chemical reactions) in liquid rocket engines. The code has a non-staggered formulation with generalized body-fitted-coordinates (BFC) capability. Higher order differencing methodologies such as MUSCL and Osher-Chakravarthy schemes are available. Turbulent flows can be modeled using any of the five turbulent models present in the code. A two-phase, two-liquid, Lagrangian spray model has been incorporated into the code. Chemical equilibrium and finite rate reaction models are available to model chemically reacting flows. The discrete ordinate method is used to model effects of thermal radiation. The code has been validated extensively against benchmark experimental data and has been applied to model flows in several propulsion system components of the SSME and the STME.

Przekwas, A. J.↗
NASA NTRS1983

Io's atmosphere - Pressure control by regolith cold trapping and surface venting

A new model for the basal pressure of Io's atmosphere is developed. This model takes into account the previously ignored fact that much of Io's surface has very high porosity, typically about 90 pct. Such porosity allows efficient subsurface cold trapping of atmospheric gases which tends to keep ambient surface pressures very low. SO2 is the only gas identified on Io, and the basal pressures for atmospheric models are usually pegged to local surface temperature via the SO2 vapor pressure equilibrium curve. Near Io's subsolar point the pressure in equilibrium with a surface SO2 frost deposit is about 1/10,000,000th bar. Porous surface models of the type developed invoke equilibrium with the colder, subsurface permafrost (at about 3-cm depth) and yield pressures of about 1/10 to the 12th bar. The subsurface cold trapping model explains many but not all observations relevant to Io's atmosphere. The new subsurface cold trapping model and the earlier surface frost equilibrium model, when taken together, provide lower and upper limits, respectively, on the basal SO2 pressure of Io's atmosphere.

Matson, D. L.↗
NASA NTRS1987

Multiple-Scale Turbulence Model

Turbulence in complex, recirculating flows not in spectral equilibrium. Report discusses mathematical fluid-flow model that includes multiple turbulence scales. Model developed for numerical prediction of confined, recirculating flows. Based on multiple-timescale concepts introduced in previous study, model takes into account turbulence in recirculating flow not in spectral equilibrium and different energy-transfer rates for eddies of different spatial scales treated separately.

Chen, C.↗
NASA NTRS1993

Mass loss from pre-main-sequence accretion disks. I - The accelerating wind of FU Orionis

We present evidence that the wind of the pre-main-sequence object FU Orionis arises from the surface of the luminous accretion disk. A disk wind model calculated assuming radiative equilibrium explains the differential behavior of the observed asymmetric absorption-line profiles. The model predicts that strong lines should be asymmetric and blueshifted, while weak lines should be symmetric and double-peaked due to disk rotation, in agreement with observations. We propose that many blueshifted 'shell' absorption features are not produced in a true shell of material, but rather form in a differentially expanding wind that is rapidly rotating. The inference of rapid rotation supports the proposal that pre-main-sequence disk winds are rotationally driven.

Calvet, Nuria↗
NASA NTRS2007

Quasi-Linear Parameter Varying Representation of General Aircraft Dynamics Over Non-Trim Region

For applying linear parameter varying (LPV) control synthesis and analysis to a nonlinear system, it is required that a nonlinear system be represented in the form of an LPV model. In this paper, a new representation method is developed to construct an LPV model from a nonlinear mathematical model without the restriction that an operating point must be in the neighborhood of equilibrium points. An LPV model constructed by the new method preserves local stabilities of the original nonlinear system at "frozen" scheduling parameters and also represents the original nonlinear dynamics of a system over a non-trim region. An LPV model of the motion of FASER (Free-flying Aircraft for Subscale Experimental Research) is constructed by the new method.

Shin, Jong-Yeob↗
NASA NTRS1989

A survey of upwind methods for flows with equilibrium and non-equilibrium chemistry and thermodynamics

Several versions of flux-vector split and flux-difference split algorithms were compared with regard to general applicability and complexity. Test computations were performed using curve-fit equilibrium air chemistry for an M = 5 high-temperature inviscid flow over a wedge, and an M = 24.5 inviscid flow over a blunt cylinder for test computations; for these cases, little difference in accuracy was found among the versions of the same flux-split algorithm. For flows with nonequilibrium chemistry, the effects of the thermodynamic model on the development of flux-vector split and flux-difference split algorithms were investigated using an equilibrium model, a general nonequilibrium model, and a simplified model based on vibrational relaxation. Several numerical examples are presented, including nonequilibrium air chemistry in a high-temperature shock tube and nonequilibrium hydrogen-air chemistry in a supersonic diffuser.

Grossman, B.↗
NASA NTRS1982

X-ray emission from young supernova remnants - Nonionization equilibrium abundances and emissivities

X-ray line emission from hot, low-density plasma in young supernova remnants is strongly enhanced by departures from ionization equilibrium. The X-ray emission from a Sedov blast wave has been calculated using a nonequilibrium evolutionary treatment of the ionization structure, and the resulting spectrum has been fitted to HEAO 2 SSS data for Tycho's remnant. These models yield dramatically different elemental abundances for heavy elements (Ne, Mg, Si, S, Ar, Ca, and Fe), compared with two-temperature component equilibrium models of Becker et al. (1979, 1980). Nonequilibrium broad-band X-ray emissivities result in lower mass determinations for the supernova ejecta. Areas of further improvement of remnant X-ray modeling are suggested.

Shull, J. M.↗
NASA NTRS1993

A wet-geology and cold-climate Mars model: Punctuation of a slow dynamic approach to equilibrium

It was suggested that Mars may have possessed a relatively warm humid climate and a vigorous hydrological cycle involving meteoric precipitation, oceans, and continental ice sheets. Baker hypothesized that these geologically active conditions may have been repeated several times; each of these dynamic epochs was followed by a collapse of the climate and hydrologic cycle of Mars into essentially current conditions, completing what is termed a 'Baker cycle'. The purpose is to present an endmember possibility that Martian glacial landscapes, including some that were previously considered to have formed under warm climatic conditions, might be explained by processes compatible with an extremely cold surface. Two aspects of hypothesized Martian glacial terrains were cited as favoring a warm climate during Baker cycles: (1) the formation of some landscapes, including possible eskers, tunnel channels, drumlins, and outwash plains, appears to have required liquid water, and (2) a liquid-surfaced ocean was probably necessary to feed the glaciers. The requirement for liquid water, if these features were correctly interpreted, is difficult to avoid; it is entirely possible that a comparatively warm climate was involved, but it is not clear that formation of landforms by wet-based glaciers actually requires a warm climate. Even less certain is the supposed requirement for liquid oceans. Formation of glaciers only requires a source of water or ice to supply an amount of precipitation that exceeds losses due to melting and sublimation. At Martian temperatures precipitation is very low, but so are melting and sublimation, so a large body of ice that is unstable with respect to sublimation may take the role of Earth's oceans in feeding the glaciers. Recent models suggest that even current Martian polar caps, long thought to be static bodies of ice and dust, might actually be slow-moving, cryogenic continental glaciers. Is it possible that subglacial processes beneath cryogenic (but wetbased) ice sheets formed the hypothesized Martian glacial landscapes?

Kargel, J. S.↗
NASA NTRS2019

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA's exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA's Entry Systems Modeling (ESM) Project, which serves as the Agency's only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specific models for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures.

Barnhardt, Michael D.↗
NASA NTRS1995

Non-LTE line-blanketed model atmospheres of hot stars. 1: Hybrid complete linearization/accelerated lambda iteration method

A new munerical method for computing non-Local Thermodynamic Equilibrium (non-LTE) model stellar atmospheres is presented. The method, called the hybird complete linearization/accelerated lambda iretation (CL/ALI) method, combines advantages of both its constituents. Its rate of convergence is virtually as high as for the standard CL method, while the computer time per iteration is almost as low as for the standard ALI method. The method is formulated as the standard complete lineariation, the only difference being that the radiation intensity at selected frequency points is not explicity linearized; instead, it is treated by means of the ALI approach. The scheme offers a wide spectrum of options, ranging from the full CL to the full ALI method. We deonstrate that the method works optimally if the majority of frequency points are treated in the ALI mode, while the radiation intensity at a few (typically two to 30) frequency points is explicity linearized. We show how this method can be applied to calculate metal line-blanketed non-LTE model atmospheres, by using the idea of 'superlevels' and 'superlines' introduced originally by Anderson (1989). We calculate several illustrative models taking into accont several tens of thosands of lines of Fe III to Fe IV and show that the hybrid CL/ALI method provides a robust method for calculating non-LTE line-blanketed model atmospheres for a wide range of stellar parameters. The results for individual stellar types will be presented in subsequent papers in this series.

Hubeny, I.↗
NASA NTRS1999

Tropical Convection and Climate Processes in a Cumulus Ensemble Model

Local convective-radiative equilibrium states of the tropical atmosphere are determined by the following external forcing: 1) Insolation, 2) Surface heat and moisture exchanges (primarily radiation and evaporation), 3) Heating and moistening induced by large-scale circulation. Understanding the equilibrium states of the tropical atmosphere in different external forcing conditions is of vital importance for studying cumulus parameterization, climate feedbacks, and climate changes. We extend our previous study using the Goddard Cumulus Ensemble (GCE) Model which resolves convective-radiative processes more explicitly than global climate models do. Several experiments are carried out under fixed insolation and sea surface temperature. The prescribed SST consists of a uniform warm pool (29C) surrounded by uniform cold SST (26C). The model produces "Walker"-type circulation with the ascending branch of the model atmosphere more humid than the descending part, but the vertically integrated temperature does not show a horizontal gradient. The results are compared with satellite measured moisture by SSM/I (Special Sensor Microwave/Imager) and temperature by MSU in the ascending and descending tropical atmosphere. The vertically integrated temperature and humidity in the two model regimes are comparable to the observed values in the tropics.

Sui, Chung-Hsiung↗
NASA NTRS1987

Statistical equilibrium in cometary C2. IV - A 10 level model including singlet-triplet transitions

Resonance fluorescence theory was used to calculate the population distribution in the energy states of the C2 molecule in comets. Ten electronic states, each with 14 vibrational states, were used in the calculations. These new calculations differ from earlier work in terms of additional electronic levels and the role of singlet-triplet transitions between the b and X levels. Since transition moments are not known, calculations are made of observable flux ratios for an array of possible values. Comparison with existing observations indicates that the a-X transition is very important, and there is marginal indication that the b-X transition is present. Swan band sequence flux ratios at large heliocentric distance are needed, as are accurate Mulliken/Swan and Phillips/Ballik-Ramsay (1963) observations.

Krishna Swamy, K. S.↗
NASA NTRS2022

LGM Paleoclimate Constraints Inform Cloud Parameterizations and Equilibrium Climate Sensitivity in CESM2

The Community Earth System Model version 2 (CESM2) simulates a high equilibrium climate sensitivity (ECS > 5°C) and a Last Glacial Maximum (LGM) that is substantially colder than proxy temperatures. In this study, we examine the role of cloud parameterizations in simulating the LGM cooling in CESM2. Through substituting different versions of cloud schemes in the atmosphere model, we attribute the excessive LGM cooling to the new CESM2 schemes of cloud microphysics and ice nucleation. Further exploration suggests that removing an inappropriate limiter on cloud ice number (NoNimax) and decreasing the time-step size (substepping) in cloud microphysics largely eliminate the excessive LGM cooling. NoNimax produces a more physically consistent treatment of mixed-phase clouds, which leads to an increase in cloud ice content and a weaker shortwave cloud feedback over mid-to-high latitudes and the Southern Hemisphere subtropics. Microphysical substepping further weakens the shortwave cloud feedback. Based on NoNimax and microphysical substepping, we have developed a paleoclimate-calibrated CESM2 (PaleoCalibr), which simulates well the observed twentieth century warming and spatial characteristics of key cloud and climate variables. PaleoCalibr has a lower ECS (∼4°C) and a 20% weaker aerosol-cloud interaction than CESM2. PaleoCalibr represents a physically more consistent treatment of cloud microphysics than CESM2 and is a valuable tool in climate change studies, especially when a large climate forcing is involved. Our study highlights the unique value of paleoclimate constraints in informing the cloud parameterizations and ultimately the future climate projection.

equilibrium climate sensitivity↗
NASA NTRS1991

Transient cryogenic liquid discharge in normal and micro-gravity

Transient discharge of liquid cryogens from partially filled tanks under microgravity is considered. The discharge is made to a low-pressure site through a straight tube. The ultimate objective of the study is to predict the transient mass flow rate. In addition, the pressure history in the tank is of interest, since the pressure falls below the triple point toward the end of the dump, causing solidification in the tank and/or the discharge line. As first steps, experiments have been performed in normal gravity and attempts are made to predict the transient rate and pressure history under this condition. Experiments are performed with liquid nitrogen dumped to sites with pressures below the triple point. Homogeneous equilibrium and nonequilibrium two-phase critical flow models are compared to the experimental results. Equilibrium and nonequilibrium tank pressure history models are also compared to experimental data.

Purwin, T. P.↗