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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Low-energy X-ray line emission from IC 443

HEAO 1 observations of the spectrum of the supernova remnant IC 443 in the energy range of 0.4-3 keV reveal the presence of a complex structure suggestive of emission from Fe XVIII-XX, S XV-XVI, and Si XIII-XIV ions. The best electron temperature in the band for simple model fits is about (6-11) x 10 to the 6th K. Raymond and Smith collisional equilibrium emission models do not adequately fit the data. These results are discussed in terms of possible nonequilibrium effects in the remnant.

Charles, P. A.↗
NASA NTRS2013

Prominence Mass Supply and the Cavity

A prevalent but untested paradigm is often used to describe the prominence-cavity system; the cavity is under-dense because it it evacuated by supplying mass to the condensed prominence. The thermal non-equilibrium (TNE) model of prominence formation offers a theoretical framework to predict the thermodynamic evolutin of the prominence and the surrounding corona. We examine the evidence for a prominence-cavity connection by comparing the TNE model and diagnostics of dynamic extreme ultraviolet (EUV) emission surrounding the prominence, specifically prominence horns. Horns are correlated extensions of prminence plasma and coronal plasma which appear to connect the prominence and cavity. The TNE model predicts that large-scale brightenings will occur in the Solar Dynamics Observatory Atmospheric Imaging Assembly 171 A badpass near he prominence that are associated with the cooling phase of condensation formation. In our simulations, variations in the magnitude of footpoint heating lead to variations in the duration, spatial scale, and temporal offset between emission enhancements in the other EUV bandpasses. While these predictions match well a subset of the horn observations, the range of variations in the observed structures is not captured by the model. We discuss the implications of one-dimensional loop simulations for the three-dimensional time-averaged equilibrium in the prominence and the cavity. Evidence suggests that horns are likely caused by condensing prominence plasma, but the larger question of whether this process produces a density-depleted cavity requires a more tightly constrained model of heating and better knowledge of the associated magnetic structure.

PROMINENCES↗
NASA NTRS1995

A new criterion for bar-forming instability in rapidly rotating gaseous and stellar systems. 1: Axisymmetric form

We analyze previous results on the stability of uniformly and differentialy rotating, self-gravitating, gaseous and stellar, axisymmetric systems to derive a new stability criterion for the appearance of torodial, m = 2 intermediate or I-modes and bar modes. In the process, we demonstrate that the bar modes in stellar systems and the m = 2 I-modes in gaseous systems have many common physical characteristics and only one substantial difference: because of the anisotropy of the stress tensor, dynamical instability sets in at lower rotation in stellar systems. This difference is reflected also in the new stability criterion. The new stability parameter alpha equals (T(sub J))/(absolute value of W) is formulated first for uniformly rotating systems and is based on the angular momentum content rather than on the energy content of a system. (T(sub J) is defined as ((L)(Omega(sub J)))/2; L is the total angular momentum; Omega(sub J) is the Jeans frequency introduced by self-gravity; and W is the total gravitational potential energy.) For stability of stellar systems alpha less than or equal to 0.254-0.258 while alpha less than or equal to 0.341-0.354 for stability of gaseous systems. For uniform rotation, one can write alpha = ((ft)/2)(exp 1/2), where t is defined as T/(absolute value of W), T is the total kinetic energy due to rotation, and f is a function characteristic of the topology/connectedness and the geometric shape of a system. Equivalently, alpha equals t/(chi), where chi is defined as Omega/Omega(sub J) and Omega is the rotation frequency. Using these forms, alpha can be extended to and calculated for a variety of differentially rotating, gaseous and stellar, axisymmetric disk and spheroidal models whose equilibrium structures and stability characteristics are known. In this paper, we also estimate alpha for gaseous torodial models and for stellar disk systems embedded in an inert or responsive 'halo.' We find that the new stability criterion holds equally well for all these previously published axisymmetric models.

Christodoulou, Dimitris M.↗
NASA NTRS2019

Behavior of Triple Langmuir Probes in Non-Equilibrium Plasmas

A model of current collection in Langmuir probes is used to investigate the effect representative non- equilibrium plasmas under various conditions have on the electron temperature and number density that would be calculated through analysis of the probe collection characteristics. The model uses the distribution function to calculate the charged particle flux to a probe and then, for fixed applied voltages, the current values that satisfy continuity in the probes are determined. The triple probe is not scanned in voltage, so there is no practical way using experimental triple probe data to determine if the plasma is in equilibrium. As a consequence, a triple probe analysis typically relies on the assumption that the plasma is in equilibrium. Proceeding from this point, the numerically-generated non-equilibrium triple probe data are analyzed assuming that the plasma is in equilibrium, with the data compared to the initial distribution function inputs of plasma temperature and number density to determine the effect the non-equilibrium distribution has on plasma measurements. The temperature and number density are both significantly affected when a fraction of the particles in the distribution are shifted from the equilibrium configuration into the non-equilibrium part of the distribution function. For all instances studied, the computed electron temperature and number density are extremely sensitive to small deviations from equilibrium ( ≤ 5% of the plasma shifted into the non-equilibrium function). Shifting more of the plasma into the non-equilibrium distribution beyond this initial level does not produce a significant additional shift in the computed plasma properties.

Polzin, Kurt A.↗
NASA NTRS1974

Application of an invariant second-order closure model to compressible turbulent shear layers

A second-order closure model for two-dimensional, compressible shear flows is investigated using the invariant modeling technique developed by Donaldson. The invariant model parameters were originally selected by comparison of model predictions with critical experimental data on basic incompressible flows. Additional modeling for compressible flows has been introduced. A number of different shear flows such as the free shear layer, flat plate boundary layer and a simulation of the mixing region of a chemical laser are computed. The model is consistent with first-order closure turbulence models for equilibrium flows and is further capable of predicting nonequilibrium flows that cannot be correctly solved by eddy viscosity models. The results are in generally good agreement with experimental measurements but suggest the need for inclusion of a Mach number dependent model of the pressure diffusion terms in order to adequately represent high speed compressible flows.

Varma, A. K.↗
NASA NTRS1988

Computations of ideal and real gas high altitude plume flows

In the present work, complete flow fields around generic space vehicles in supersonic and hypersonic flight regimes are studied numerically. Numerical simulation is performed with a flux-split, time asymptotic viscous flow solver that incorporates a generalized equilibrium chemistry model. Solutions to generic problems at various altitude and flight conditions show the complexity of the flow, the equilibrium chemical dissociation and its effect on the overall flow field. Viscous ideal gas solutions are compared against equilibrium gas solutions to illustrate the effect of equilibrium chemistry. Improved solution accuracy is achieved through adaptive grid refinement.

Feiereisen, William J.↗
NASA NTRS1990

Laminar flamelet modeling of turbulent diffusion flames

In modeling turbulent combustion, decoupling the chemistry from the turbulence is of great practical significance. In cases in which the equilibrium chemistry model breaks down, laminar flamelet modeling (LFM) is a promising approach to decoupling. Here, the validity of this approach is investigated using direct numerical simulation of a simple chemical reaction in two-dimensional turbulence.

Mell, W. E.↗
NASA NTRSDate not supplied

WMLES for the Fifth High-Lift PredictionWorkshop Cases Using FUN3D

This paper presents solution assessments and grid convergence studies for the test cases outlined in the Fifth High-Lift Prediction Workshop (HLPW-5), focusing on the high-lift Common Research Models (CRM-HL). The study utilizes a wall-modeled large-eddy simulation (WMLES) methodology developed in the unstructured-grid, node-centered flow solver FUN3D. The second-order accurate simulations conducted in this study utilize a finite-volume spatial discretization and an implicit temporal scheme. Large-scale turbulent features are resolved away from the wall, with small-scale effects captured by the Vreman subgrid-scale model. An equilibrium wall function uses the first grid point off the wall serving as the critical interface between the wall model and the large-eddy simulation region, thus requiring careful placement in grid design. WMLES solutions are assessed for HLPW-5 cases, including a clean wing-body configuration and geometry-buildup configurations corresponding to the 5.1\% ONERA CRM-HL model. Grid-convergence studies are systematically conducted using uniformly refined grids. Moreover, simulation results and grid sensitivity are presented for the NASA 5.2\% CRM-HL configuration at both moderate and flight-scale Reynolds numbers. Overall, WMLES results are satisfactory and agree well with available experimental data, especially on sufficiently fine grids.

high-lift aerodynamics↗
NASA NTRS2013

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗
NASA NTRS1993

An experimental study of trace element partitioning between olivine, orthopyroxene and melt in chondrules - Equilibrium values and kinetic effects

Mineral/melt partition coefficients were measured using an ion microprobe for 32 elements in orthopyroxene and olivine in equilibrium and dynamic crystallization experiments on compositions corresponding to chondrules. The mineral/melt partition coefficients calculated from the measured concentrations for both olivine and orthopyroxene show very little change between equilibrium experiments and dynamic experiments with cooling rates of up to 100 C/h. The results provide a self-consistent set of partition coefficients that can be used in thermodynamic models of equilibrium and kinetic partitioning between olivine, orthopyroxene, and melt. These data can be used in models of partial melting and crystal fractionation in olivine- and orthopyroxene-rich systems, such as chondrules. The results may also be applicable to mantle peridotites, komatiitic and picritic lavas, and ultramafic intrusions.

Kennedy, A. K.↗
NASA NTRS1995

Modeling of turbulent chemical reaction

Viewgraphs are presented on modeling turbulent reacting flows, regimes of turbulent combustion, regimes of premixed and regimes of non-premixed turbulent combustion, chemical closure models, flamelet model, conditional moment closure (CMC), NO(x) emissions from turbulent H2 jet flames, probability density function (PDF), departures from chemical equilibrium, mixing models for PDF methods, comparison of predicted and measured H2O mass fractions in turbulent nonpremixed jet flames, experimental evidence of preferential diffusion in turbulent jet flames, and computation of turbulent reacting flows.

Chen, J.-Y.↗
NASA NTRS2008

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗
NASA NTRS1988

Two classes of models for temporal variability of the Io torus

Two classes of spatially uniform models for the Io torus are analyzed for the stability and temporal variability of the torus. It is shown that stable-equilibrium torus models exist in two classes of models, one with nonlinear radial transport loss of plasma (e.g., from centrifugally driven interchanges) and one with linear transport (e.g., from externally driven interchanges). The two classes of models make substantially different predictions regarding temporal variability of the Io torus. Observational tests are proposed that make it possible to discriminate between the two classes.

Cheng, Andrew F.↗
NASA NTRS1994

Analysis of helium-rich subdwarf O stars. 1: NLTE models, methods, and fits for 21 Palomar Green survey sdOs

Atmospheric parameters for 21 helium-rich hot subdwarf O stars from the Palomar Green survey are found from fits of non-local thermodynamic equilibrium (NLTE) models to optical spectra. About 250 new NLTE models in the parameter range T(sub eff) from 35,000 to 65,000 K, log (g) from 4.0 to 6.5, and epsilon(He) from 50% He to 99% He have been calculated. A fit for each object is presented. Estimated distances and luminosities are calculated, assuming a mass of 0.5 solar masses. Large distances above the Galactic plane are found, and this implies that the majority of the sdO stars belong to a different stellar population than the sdB, planetary nebulae and white dwarf stars. Possibilities may include the halo and a 'thick disk.' Kinematical data from the literature are also discussed in view of the found distances. Five of the stars have been observed, modeled, and analyzed by Dreizler et al., and significant differences exist between their results and ours for these stars. The reason for these differences is not yet known, but the problem does not alter the conclusions about population membership that we present.

Thejll, P.↗
NASA NTRS1975

Application of a general boundary layer analysis to turbulent boundary layers subjected to strong favorable pressure gradients

Theoretical predictions of turbulent boundary layer development under the influence of strong favorable pressure gradients made using a finite-difference calculation procedure are compared to experimental data. Comparisons are presented for low speed flows with and without wall heat transfer as well as for supersonic flows with adiabatic walls. The turbulence model used is governed by an integral form of the turbulence kinetic energy equation and the results are compared with predictions made using a conventional equilibrium turbulence model based upon Prandtl's mixing length, a Clauser-type eddy viscosity model used by Cebecci and Mosinskis, and a two-equation turbulence energy model of Launder and Jones.

Kreskovsky, J. P.↗
NASA NTRS1999

Carbon Monoxide Formation in SN 1987A

The formation of carbon monoxide in the supernova SN 1987A at 200 days is investigated. Rather than the one-zone that have previous been employed, this work utilizes a radially dependent radiative transport model to compute the temperature and ionization structure and the destruction of CO and CO+ by ultraviolet radiation. The CO profile is computed assuming chemical equilibrium. Two models for the density and abundance profiles are examined: model 10H (unmixed) and model 10HMM (mixed) of Pinto & Woosley. Sensitivity to adopted rates and temperature as it might be affected by CO cooling is examined. The models give a total mass of 10-4 M_sun at 200 days for the unmixed model, which is comparable with the estimated observed abundance by Syromilio et al. but nearly 2 orders of magnitude less than the estimated observed and that computed in thermal-chemical models by Liu & Dalgarno. The effect of different model assumptions and results concerning ionization structure and radiative transfer are discussed. We confirm that CO+ is not expected to be produced in significant amounts and that the amount of CO is sensitive to the degree of the mixing of the composition of various elements present in the ejecta.

Gearhart, Rob A.↗
NASA NTRS1992

A first-order global model of Late Cenozoic climatic change: Orbital forcing as a pacemaker of the ice ages

The development of a theory of the evolution of the climate of the earth over millions of years can be subdivided into three fundamental, nested, problems: (1) to establish by equilibrium climate models (e.g., general circulation models) the diagnostic relations, valid at any time, between the fast-response climate variables (i.e., the 'weather statistics') and both the prescribed external radiative forcing and the prescribed distribution of the slow response variables (e.g., the ice sheets and shelves, the deep ocean state, and the atmospheric CO2 concentration); (2) to construct, by an essentially inductive process, a model of the time-dependent evolution of the slow-response climatic variables over time scales longer than the damping times of these variables but shorter than the time scale of tectonic changes in the boundary conditions (e.g., altered geography and elevation of the continents, slow outgassing, and weathering) and ultra-slow astronomical changes such as in the solar radiative output; and (3) to determine the nature of these ultra-slow processes and their effects on the evolution of the equilibrium state of the climatic system about which the above time-dependent variations occur. All three problems are discussed in the context of the theory of the Quaternary climate, which will be incomplete unless it is embedded in a more general theory for the fuller Cenozoic that can accommodate the onset of the ice-age fluctuations. We construct a simple mathematical model for the Late Cenozoic climatic changes based on the hypothesis that forced and free variations of the concentration of atmospheric greenhouse gases (notably CO2), coupled with changes in the deep ocean state and ice mass, under the additional 'pacemaking' influence of earth-orbital forcing, are primary determinants of the climate state over this period. Our goal is to illustrate how a single model governing both very long term variations and higher frequency oscillatory variations in the Pleistocene can be formulated with relatively few adjustable parameters.

Saltzman, Barry↗