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 19 records

A photochemical equilibrium model for ionospheric conductivity

A photochemical equilibrium model of the high-latitude ionosphere has been developed. This model provides densities of the ionospheric constituents, N2(+), O2(+), O(+), and NO(+), from 85 km to approximately 220 km. These densities are then used to calculate Pedersen and Hall conductivities. A comparison of the model results with Arecibo and Chatanika radar observations was made, covering periods of solar minimum and solar maximum. The comparison showed the model to predict ionospheric densities to within 50 percent and conductivities to within 40 percent in the illuminated portion of the ionosphere. In regions of electron precipitation, the model showed good agreement with measurements. Results of this study indicate the following: (1) ionospheric conductivity increases by a factor of about 1.6 from solar minimum to solar maximum conditions, (2) the portion of the ionosphere above 170 km can contribute as much as 40 percent during daylight and 80 percent during nighttime to the total height-integrated Pedersen conductivity, and (3) the ratio of the height-integrated Hall to Pedersen conductivities is approximately 1.1-1.3 for sunlit conditions; this is appreciably lower than the value of 2 found in previous studies. These and other factors indicate that, under certain conditions, the height-integrated Pedersen conductivity may be as much as 2-3 times larger than previously reported.

Rasmussen, C. E.↗

A numerical code for a three-dimensional magnetospheric MHD equilibrium model

Two dimensional and three dimensional MHD equilibrium models were begun for Earth's magnetosphere. The original proposal was motivated by realizing that global, purely data based models of Earth's magnetosphere are inadequate for studying the underlying plasma physical principles according to which the magnetosphere evolves on the quasi-static convection time scale. Complex numerical grid generation schemes were established for a 3-D Poisson solver, and a robust Grad-Shafranov solver was coded for high beta MHD equilibria. Thus, the effects were calculated of both the magnetopause geometry and boundary conditions on the magnetotail current distribution.

Voigt, G.-H.↗

Particle orbits in two-dimensional equilibrium models for the magnetotail

Assuming that there exist an equilibrium state for the magnetotail, particle orbits are investigated in two-dimensional kinetic equilibrium models for the magnetotail. Particle orbits in the equilibrium field are compared with those calculated earlier with one-dimensional models, where the main component of the magnetic field (Bx) was approximated as either a hyperbolic tangent or a linear function of z with the normal field (Bz) assumed to be a constant. It was found that the particle orbits calculated with the two types of models are significantly different, mainly due to the neglect of the variation of Bx with x in the one-dimensional fields.

Karimabadi, H.↗

Radiative-convective equilibrium models of Jupiter and Saturn

Radiative-convective equilibrium models for Jupiter and Saturn have been produced in a study concentrating on the stratospheric energy balance and the possible role of aerosol heating. These models are compared directly with the thermal structure profiles obtained from Voyager radio occultation measurements. The method is based on a straightforward flux divergence formulation derived from earlier work. The balance between absorbed and emitted energies is computed iteratively at each level in the atmosphere, assuming local thermodynamic equilibrium and employing a standard treatment of opacities. Results for Jupiter indicate that a dust-free model furnishes a good mean thermal profile for the stratosphere when compared with the Voyager 1 radio occultation measurements. Observations of the equatorial region exhibit periodic vertical structure. The Saturn models are relatively simple and in good agreement with the Voyager 2 radio occultation temperature profiles at all levels. Aerosol heating played a minor role in Saturn's midlatitude stratospheric energy balance at the time of the Voyager 2 encounter.

Appleby, J. F.↗

He I lines in B stars - Comparison of non-local thermodynamic equilibrium models with observations

Profiles of He gamma-gamma 4026, 4387, 4471, 4713, 5876, and 6678 have been obtained in 17 stars of spectral type B0-B5. Parameters of the nonlocal thermodynamic equilibrium models appropriate to each star are determined from the Stromgren index and fits to H-alpha line profiles. These parameters yield generally good fits to the observed He I line profiles, with the best fits being found for the blue He I lines where departures from local thermodynamic equilibrium are relatively small. For the two red lines it is found that, in the early B stars and in stars with log g less than 3.5, both lines are systematically stronger than predicted by the nonlocal thermodynamic equilibrium models.

Heasley, J. N.↗

Radiative equilibrium model of Titan's atmosphere

The present global radiative equilibrium model for the Saturn satellite Titan is restricted to the two-stream approximation, is vertically homogeneous in its scattering properties, and is spectrally divided into one thermal and two solar channels. Between 13 and 33% of the total incident solar radiation is absorbed at the planetary surface, and the 30-60 ratio of violet to thermal IR absorption cross sections in the stratosphere leads to the large temperature inversion observed there. The spectrally integrated mass absorption coefficient at thermal wavelengths is approximately constant throughout the stratosphere, and approximately linear with pressure in the troposphere, implying the presence of a uniformly mixed aerosol in the stratosphere. There also appear to be two regions of enhanced opacity near 30 and 500 mbar.

Samuelson, R. E.↗

Ionization and thermal equilibrium models for O star winds based on time-independent radiation-driven wind theory

Ab initio ionization and thermal equilibrium models are calculated for the winds of O stars using the results of steady state radiation-driven wind theory to determine the input parameters. Self-consistent methods are used for the roles of H, He, and the most abundant heavy elements in both the statistical and the thermal equilibrium. The model grid was chosen to encompass all O spectral subtypes and the full range of luminosity classes. Results of earlier modeling of O star winds by Klein and Castor (1978) are reproduced and used to motivate improvements in the treatment of the hydrogen equilibrium. The wind temperature profile is revealed to be sensitive to gross changes in the heavy element abundances, but insensitive to other factors considered such as the mass-loss rate and velocity law. The reduced wind temperatures obtained in observing the luminosity dependence of the Si IV lambda 1397 wind absorption profile are shown to eliminate any prospect of explaining the observed O VI lambda 1036 line profiles in terms of time-independent radiation-driven wind theory.

Drew, J. E.↗

Radiative-convective equilibrium models of Uranus and Neptune

The present study of Uranus and Neptune radiative-convective equilibrium models gives emphasis to such aspects of the stratospheric energy balance as the influence of aerosol heating and convective penetration. The results obtained for Uranus imply that a 'continuum absorber' may be a significant factor in the stratosphere despite the great distance from the sun. The results obtained for Neptune show that such a continuum absorber could significantly contribute to the energy balance within a localized stratospheric region, although it probably cannot furnish sufficient power to account for the observed IR spectrum irrespective of its vertical distribution. Attention is accordingly given to the 'convective penetration' that could arise under such rapid vertical mixing that CH4's condensation cannot occur before the gas is carried above the condensation region.

Appleby, J. F.↗

Radiative-convective equilibrium models of Uranus and Neptune

The present study of Uranus and Neptune radiative-convective equilibrium models gives emphasis to such aspects of the stratospheric energy balance as the influence of aerosol heating and convective penetration. The results obtained for Uranus imply that a continuum absorber may be a significant factor in the stratosphere despite the great distance from the sun. The results obtained for Neptune show that such a continuum absorber could significantly contribute to the energy balance within a localized stratospheric region, although it probably cannot furnish sufficient power to account for the observed IR spectrum irrespective of its vertical distribution. Attention is accordingly given to the convective penetration that could arise under such rapid vertical mixing that CH4's condensation cannot occur before the gas is carried above the condensation region.

Appleby, J. F.↗

Infrared polar brightenings on Jupiter. V - A thermal equilibrium model for the north polar hot spot

Voyager IRIS instrument records of the IR hydrocarbon emissions from Jupiter's north polar region are presently studied to determine the spatial and other characteristics of the north polar hot spot. Attention is given to a thermal equilibrium model that exploits the asymmetry found in 7.8-micron emission of stratospheric methane with respect to system III longitude in order to estimate stratospheric zonal wind velocity. This model accurately predicts the observed asymmetry in acetylene's 13.6-micron emission; this requires, however, enhanced acetylene abundance in the hot spot, as well as ethane depletion. Energetic charged particles are suggested to be the most probable cause of these effects.

Halthore, Rangasayi↗

Partial equilibrium model for predicting concentration of CO in combustion

A simple kinetic model, based upon the concept of partial equilibrium, is developed for predicting carbon monoxide concentrations in such steady-flow hydrocarbon-air combustion systems as gas-turbine combustors. The only two kinetic constraints used in the model are on the rate of change of the total number of gaseous particles in the system, and on the rate of change of the CO concentration. The accuracy of the model is verified by comparison with experimental results obtained by burning kerosene with air in a 130,000 Btu/hr atmospheric-pressure steady-flow burner. Burnt-gas cooling rates of 1,000,000 K/sec, about the magnitude encountered in the dilution zone of gas-turbine combustors, were obtained with a compact water-cooled heat exchanger. Good agreement between measured CO concentrations and the values predicted by the partial-equilibrium model was obtained.

Morr, A. R.↗

An equilibrium model for the coupled ocean-atmosphere boundary layer in the tropics

An atmospheric convective boundary layer (CBL) model is coupled to an ocean mixed-layer (OML) model in order to study the equilibrium state of the coupled system in the tropics, particularly in the Pacific region. The equilibrium state of the coupled system is solved as a function of sea-surface temperature (SST) for a given surface wind and as a function of surface wind for a given SST. It is noted that in both cases, the depth of the CBL and OML increases and the upwelling below the OML decreases, corresponding to either increasing SST or increasing surface wind. The coupled ocean-atmosphere model is solved iteratively as a function of surface wind for a fixed upwelling and a fixed OML depth, and it is observed that SST falls with increasing wind in both cases. Realistic gradients of mixed-layer depth and upwelling are observed in experiments with surface wind and SST prescribed as a function of longitude.

Sui, C.-H.↗

Phase equilibrium modeling for high temperature metallization on GaAs solar cells

Recent trends in performance specifications and functional requirements have brought about the need for high temperature metallization technology to be developed for survivable DOD space systems and to enhance solar cell reliability. The temperature constitution phase diagrams of selected binary and ternary systems were reviewed to determine the temperature and type of phase transformation present in the alloy systems. Of paramount interest are the liquid-solid and solid-solid transformations. Data are being utilized to aid in the selection of electrical contact materials to gallium arsenide solar cells. Published data on the phase diagrams for binary systems is readily available. However, information for ternary systems is limited. A computer model is being developed which will enable the phase equilibrium predictions for ternary systems where experimental data is lacking.

Chung, M. A.↗

Modelling Equilibrium and Fractional Crystallization in the System MgO-FeO-CaO-Al2O3-SiO2

A mathematical modelling technique for use in petrogenesis calculations in the system MgO-FeO-CaO-Al2O3-SiO2 is reported. Semiempirical phase boundary and elemental distribution information was combined with mass balance to compute approximate equilibrium crystallization paths for arbitrary system compositions. The calculation is applicable to a range of system compositions and fractionation calculations are possible. The goal of the calculation is the computation of the composition and quantity of each phase present as a function of the degree of solidification. The degree of solidification is parameterized by the heat released by the solidifying phases. The mathematical requirement for the solution of this problem is: (1) An equation constraining the composition of the magma for each solid phase in equilibrium with the liquidus phase, and (2) an equation for each solid phase and each component giving the distribution of that element between that phase and the magma.

Herbert, F.↗

A variational approach to the question of temporal stability of equilibrium models of solar prominences. I - The formal theory

A Lagrangian approach to the equations describing small amplitude departures from equilibrium of solar prominences is presented. Seven quantities are derived which, by analogy to the concepts of energy, momentum, and angular momentum, are conserved under conditions corresponding to ignorable coordinates of classical mechanics. It is expected that these conservation laws will be useful as criteria of accuracy in obtaining eigenfrequencies for the perturbation equations when numerical techniques are used. The method is capable of handling complicated spatial dependences and has internal criteria for assessing the accuracy of any given numerical convergence scheme.

Lerche, I.↗

Aqueous Alteration on the Ryugu Parent Body Constrained by Chemical Equilibrium Models

The composition of solid samples returned from the asteroid 162173 Ryugu byHayabusa2 indicates aqueous alteration of rocks that formed a parent body of Ryugu [1, 2]. As on parent bodies of carbonaceous chondrites [3-5], Ryugu materials could have went through (1) accretion of anhydrous and reduced rocky grains together with water-rich ices, (2) melting of ices in the body’s interior through release of radiogenic heat [6], and (3)aqueous alteration of rocky and organic materials that produced secondary minerals, salt- and organic-bearing water solutions, gases, and altered organic matter. Although the bulk chemical composition and major mineralogy (phyllosilicates, magnetite, pyrrhotite, carbonates, phosphates, etc.) of Ryugu [1,2]are similar to those of CI carbonaceous chondrites [3-5, 7], the lack of sulfates and a deficiency of oxygen in the returned samples suggest a distinct alteration pathway. Chondritic materials [3-5] do not indicate complete chemical equilibration though low-temperature parent body processes. However, calculations of chemical equilibria in modeling rock-water-gas type systems[e.g., 8, 9], could assess deviations of chondritic materials from equilibrium conditions and constrain partitioning of chemical elements between solid, aqueous, and gas phase for different stages of alteration. Here we constrained speciation of solid, aqueous, and gaseous phases during aqueous alteration on the Ryugu parent body through chemical equilibrium calculations.

M Yu Zolotov↗

A two-phase restricted equilibrium model for combustion of metalized solid propellants

An Eulerian-Lagrangian two-phase approach was adopted to model the multi-phase reacting internal flow in a solid rocket with a metalized propellant. An Eulerian description was used to analyze the motion of the continuous phase which includes the gas as well as the small (micron-sized) particulates, while a Lagrangian description is used for the analysis of the discrete phase which consists of the larger particulates in the motor chamber. The particulates consist of Al and Al2O3 such that the particulate composition is 100 percent Al at injection from the propellant surface with Al2O3 fraction increasing due to combustion along the particle trajectory. An empirical model is used to compute the combustion rate for agglomerates while the continuous phase chemistry is treated using chemical equilibrium. The computer code was used to simulate the reacting flow in a solid rocket motor with an AP/HTPB/Al propellant. The computed results show the existence of an extended combustion zone in the chamber rather than a thin reaction region. The presence of the extended combustion zone results in the chamber flow field and chemical being far from isothermal (as would be predicted by a surface combustion assumption). The temperature in the chamber increases from about 2600 K at the propellant surface to about 3350 K in the core. Similarly the chemical composition and the density of the propellant gas also show spatially non-uniform distribution in the chamber. The analysis developed under the present effort provides a more sophisticated tool for solid rocket internal flow predictions than is presently available, and can be useful in studying apparent anomalies and improving the simple correlations currently in use. The code can be used in the analysis of combustion efficiency, thermal load in the internal insulation, plume radiation, etc.

Sabnis, J. S.↗

Errors in Sounding of the Atmosphere Using Broadband Emission Radiometry (SABER) Kinetic Temperature Caused by Non-Local Thermodynamic Equilibrium Model Parameters

The vast set of near global and continuous atmospheric measurements made by the SABER instrument since 2002, including daytime and nighttime kinetic temperature (T(sub k)) from 20 to 105 km, is available to the scientific community. The temperature is retrieved from SABER measurements of the atmospheric 15 micron CO2 limb emission. This emission separates from local thermodynamic equilibrium (LTE) conditions in the rarefied mesosphere and thermosphere, making it necessary to consider the CO2 vibrational state non-LTE populations in the retrieval algorithm above 70 km. Those populations depend on kinetic parameters describing the rate at which energy exchange between atmospheric molecules take place, but some of these collisional rates are not well known. We consider current uncertainties in the rates of quenching of CO2 (v2 ) by N2 , O2 and O, and the CO2 (v2 ) vibrational-vibrational exchange to estimate their impact on SABER T(sub k) for different atmospheric conditions. The T(sub k) is more sensitive to the uncertainty in the latter two and their effects depend on altitude. The T(sub k) combined systematic error due to non-LTE kinetic parameters does not exceed +/- 1.5 K below 95 km and +/- 4-5 K at 100 km for most latitudes and seasons (except for polar summer) if the Tk profile does not have pronounced vertical structure. The error is +/- 3 K at 80 km, +/- 6 K at 84 km and +/- 18 K at 100 km under the less favourable polar summer conditions. For strong temperature inversion layers, the errors reach +/- 3 K at 82 km and +/- 8 K at 90 km. This particularly affects tide amplitude estimates, with errors of up to +/- 3 K.

Garcia-Comas, Maya↗