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 793 records · Page 44

A combined optical/X-ray study of the Galaxy cluster Abell 2256

The dynamics of Abell 2256 is investigated by combining X-ray observations of the intracluster gas with optical observations of the galaxy distribution and kinematics. Magnitudes and positions are presented for 172 galaxies and new redshifts for 75. Abell 2256 is similar to the Coma Cluster in its X-ray luminosity, mass, and galaxy density. Both the X-ray surface brightness and the galaxy surface density distributions exhibit an elliptical morphology. The radial galaxy distribution is steeper than the density profile of the X-ray-emitting gas, yet the galaxy velocity dispersion is higher than the equivalent value for the gas. Under the simplest assumptions that the galaxy velocity distribution is isotropic and the gas is isothermal, the galaxies and gas cannot be in hydrostatic equilibrium in a common gravitational potential. Models consistent with available data have mass-to-light ratios which increase with radius and galaxy orbits that are anisotropic with a radial bias.

Fabricant, Daniel G.↗

Computational equations for radiating and ablating shock layers

The computational equations governing the shock-layer flowfield over a blunt body dominated by radiation and ablation are derived. The flowfield considered is that expected around the vehicles entering the atmosphere of the planet Mars or the earth on return from Mars on a sprint-type manned mission. To determine the appropriate method of approach, the radiative transfer calculations are made first using the existing computer code RASLE (Nicolet et al., 1978), which uses the assumption of equilibrium and a coarse spectral band model for radiative transfer, and the code NONEQ which accounts for thermochemical nonequilibrium and carries out line-by-line calculation of radiation. By comparing the results, it is concluded that a thermochemical nonequilibrium model and a line-by-line spectral model must be used.

Park, Chul↗

Condensation in hypersonic nitrogen wind tunnels

Experimental observations and a theoretical model for the onset and disappearance of condensation are given for hypersonic flows of pure nitrogen at M = 10, 14 and 18. Measurements include Pitot pressures, static pressures and laser light scattering experiments. These measurements coupled with a theoretical model indicate a substantial non-equilibrium supercooling of the vapor phase beyond the saturation line. Typical results are presented with implications for the design of hypersonic wind tunnel nozzles.

Lederer, Melissa A.↗

Pressure-coupled vaporization and combustion responses of liquid-fuel droplets in high-pressure environments

The dynamic responses of liquid-fuel droplet vaporization and combustion to ambient pressure oscillations are examined. The analysis is based on the complete sets of conservation equations for both gas and liquid phases, and accommodates detailed treatments of finite-rate chemical kinetics and variable properties. With a full account of thermodynamic phase equilibrium at the droplet surface, the model enables a systematic examination of the effects of ambient flow conditions on the droplet behavior. The responses of hydrocarbon fuel droplets in both sub- and super-critical environments are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗

Monte Carlo analysis of dissociation and recombination behind strong shock waves in nitrogen

Computations are presented for the relaxation zone behind strong, 1D shock waves in nitrogen. The analysis is performed with the direct simulation Monte Carlo method (DSMC). The DSMC code is vectorized for efficient use on a supercomputer. The code simulates translational, rotational and vibrational energy exchange and dissociative and recombinative chemical reactions. A model is proposed for the treatment of three body-recombination collisions in the DSMC technique which usually simulates binary collision events. The model improves previous models because it can be employed with a large range of chemical-rate data, does not introduce into the flow field troublesome pairs of atoms which may recombine upon further collision (pseudoparticles) and is compatible with the vectorized code. The computational results are compared with existing experimental data. It is shown that the derivation of chemical-rate coefficients must account for the degree of vibrational nonequilibrium in the flow. A nonequilibrium-chemistry model is employed together with equilibrium-rate data to compute the flow in several different nitrogen shock waves.

Boyd, I. D.↗

Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus

Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.

Parmentier, E. M.↗

Analysis of the upper atmosphere CO2(nu-2) vibrational temperatures retrieved from ATMOS/Spacelab 3 observations

An analysis of the upper atmospheric (80-116 km) CO2(nu-2) vibrational temperatures retrieved from atmospheric trace molecule spectroscopy (ATMOS) experiment Spacelab 3 spectra by using a nonlocal thermodynamic equilibrium (non-LTE) radiative transfer model is presented. Thermal collisions with atmospheric atomic oxygen keep this vibrational state very close to LTE up to around 100 km. Above this height, the different deviations from LTE retrieved from ATMOS/Spacelab 3 spectra for the Northern and the Southern Hemispheres are explained in terms of this collisional process and in terms of the different kinetic temperature profiles measured at those locations. From these simultaneous observations of the kinetic and CO2(nu-2) vibrational temperatures, a deactivation rate of CO2(nu-2) by O(3P) has been derived which leads to a rate coefficient value between 3 and 6 x 10 exp -12 cu cm/s and favors an independent or negative temperature dependence rate constant for the atmospheric temperature range. Cooling rates induced by the CO2 15-micron fundamental band in the upper mesosphere and lower thermosphere were derived from the simultaneous kinetic temperature, CO2(nu-2) vibrational temperature, and CO2 concentration, as measured by ATMOS/Spacelab 3, and found to be a factor of between 5 and 10 times larger than those generally accepted until very recently.

Lopez-Puertas, M.↗

Emission lines from X-ray-heated accretion disks in low-mass X-ray binaries

We investigate the structure of accretion disks illuminated by X-rays from a central compact object in a binary system. X-rays can photoionize the upper atmosphere of the disk and form an accretion disk corona (ADC) where emission lines can form. We construct a model to calculate the vertical structure and the emission spectrum of the ADC with parameters appropriate to low-mass X-ray binaries. These models are made by nonlocal thermodynamic equilibrium calculations of ion and level populations and include a large number of atomic processes for 10 cosmically abundant elements. Transfer of radiation is treated by using the escape probability formalism. The vertical temperature profile of the ADC consists of a Compton-heated region and a mid-T zone where the temperature is approximately 10(exp 6) K. A thermal instability occurs close to the disk photosphere and causes the temperature of the ADC to drop abruptly from 10(exp 6) K to several times 10(exp 4) K. The emission spectrum in the optical, ultraviolet, extreme ultraviolet, and X-ray range is discussed and compared with the observations.

Ko, Yuan-Kuen↗

A Hubble Space Telescope ultraviolet spectrum of SN 1993J

We obtained a Hubble Space Telescope (HST) ultraviolet spectrum of the Type II supernova SN 1993J in M81 on 1993 April 15. The approximately 1650-2900 A region is smoother than observed for SN 1987A and SN 1992A and lacks strong P Cygni absorptions caused by iron peak element lines. Synthetic spectra calculated using a parameterized Local Thermodynamic Equilibrium (LTE) procedure and a simple model atmosphere do not fit the UV region. Radio observations suggest that SN 1993J is embedded in a thick circumstellar envelope. The UV spectra of other supernovae that are believed to have thick circumstellar envelopes also have approximately 1650-2900 A regions lacking in strong P Cygni absorptions. Interaction of supernova ejecta and circumstellar matter may cause the smooth UV spectrum. If so, UV observations of supernovae will provide insight into the circumstellar environment of the supernova progenitors.

Jeffery, David J.↗

Observing stellar coronae with the Goddard High Resolution Spectrograph. 1: The dMe star AU microscopoii

We report on an observation of AU Mic taken with the Goddard High Resolution Spectrograph (GHRS) aboard the Hubble Space Telescope. The data consist of a rapid sequence of spectra covering the wavelength range 1345-1375 A with a spectral resolution of 10,000. The observations were originally intended to search for spectral variations during flares. No flares were detected during the 3.5 hr of monitoring. A method of reducing the noise while combining the individual spectra in the time series is described which resulted in the elimination of half of the noise while rejecting only a small fraction of the stellar signal. The resultant spectrum was of sufficient quality to allow the detection of emission lines with an integrated flux of 10(exp -15) ergs/sq cm(sec) or greater. Lines of C I, O I, O V, Cl I, and Fe XXI were detected. This is the first indisputable detection of the 1354 A Fe XXI line, formed at T approximately = 10(exp 7) K, on a star other than the Sun. The line was well resolved and displayed no significant bulk motions or profile asymmetry. From the upper limit on the observed line width, we derive an upper limit of 38 km/s for the turbulent velocity in the 10(exp 7) K plasma. An upper limit is derived for the flux of the 1349 A Fe XII line, formed at T approximately = 1.3 x 10(exp 6) K. These data are combined with contemporaneous GHRS and International Ultraviolet Explorer (IUE) data to derive the volume emission measure distribution of AU Mic over the temperature range 10(exp 4)-10(exp 7) K. Models of coronal loops in hydrostatic equilibrium are consistent with the observed volume emission measures of the coronal lines. The fraction of the stellar surface covered by the footprints of the loops depends upon the loop length and is less than 14% for lengths smaller than the stellar radius. From the upper limit to the estimated width of the Fe XXI line prfile we find that the we cannot rule out Alfven wave dissipation as a possible contributor to the required quiescent loop heating rate.

Maran, S. P.↗

EUVE spectroscopy of epsilon Canis Majoris (B2 II) from 70 to 730 A

We present spectra of the brightest stellar source of extreme ultraviolet (EUV) radiation longward of 400 A, the B2 II star, epsilon CMa. These data were taken with the three spectrometers aboard the NASA Extreme Ultraviolet Explorer satellite (EUVE) during the first cycle of pointed observations. We report on our initial studies of the continuum and line spectrum of the stellar photosphere in the 320 to 730 A region, and on the wind emission lines observed in the 170-375 A region. This is the first EUV spectrum of an early-type star, and thus makes epsilon CMa the most comprehensively observed B star from the X-ray to infrared regimes. The radiation in both the H Lyman continuum and He I continuum (shortward of 504 A) are found to be significantly greater than predicted by both Local Thermodynamic Equilibrium (LTE) and non-LTE model atmospheres. Since epsilon CMa also exhibits a mid-infrared excess, this points to the outer layers being warmer than the models indicate. The anomalously large Lyman continuum flux, combined with the very low column density measured in the direction toward this star implies that it is the dominant source of hydrogen ionization of the local interstellar medium in the immediate vicinity of the sun. All of the lines predicted to be strong from model atmospheres are present and several wind absorption features are also identified. We have detected emission lines from highly ionized iron that are consistent with the ROSAT Position Sensitive Proportional Counter (PSPC) observations if a multi-temperature emission model is used, and the assumption is made that there is significant absorption beyond that of the neutral phase of the ISM. The spectrum shows strong O III 374 A line emission produced by the Bowen flourescence mechanism, which has not previously been observed in the spectra of hot stars.

Cassinelli, J. P.↗

X-ray induced ionization in the winds of near-main-sequence O and B stars

We investigate the effects of X-rays on the ionization state of stellar winds for O and early-B stars along the main sequence. In our investigation, detailed statistical equilibrium, radiative transfer, and atomic physics models are used to compute ionization distributions for H, He, C, N, O, and Si. X-rays are modeled as a spatially distributed source within the wind, with parameters constrained by ROSAT and Einstein observations. Our results indicate that the ionization balance in the winds of stars with spectral type B0 V and later is significantly altered by the X-ray radiation field. Unlike the case of denser O star winds, where the X-rays tend to perturb the level of ionization, the ionization state of the bulk wind of early-B stars can be significantly increased by soft X-rays. We examine in detail the case of tau Sco (B0 V), which has been well-studied at UV and X-ray wavelengths. Comparisons are made between calculated ionization fractions and those deduced from UV observations. In addition, we address the sensitivity of our results to the X-ray source characteristics, the wind temperature, and the photospheric extreme ultraviolet (EUV) flux. Our results suggest the possibility that for early-B stars X-rays play a critical role in both influencing the radiation line driving force, as well as ionizing and heating the wind all the way down to the top of the photosphere.

Macfarlane, J. J.↗

AKN 564: An Unusual Component in the X-ray Spectra of NLSy1 Galaxies

We present an ASCA observation of the narrow-line Seyfert l (NLSyl) Arakelian 564. The X-ray light curve shows rapid variability, but no evidence for energy-dependence to these variations, within the 0.6 - 10 keV bandpass. A strong (EW approximately 70 eV) spectral feature is observed close to 1 keV. A similar feature has been observed in TON S180, another member of the NLSyl class of objects, but has not been observed in broad-line Seyfert galaxies. The feature energy suggests a large contribution from Fe L-shell lines but its intensity is difficult to explain in terms of emission and/or absorption from photoionized gas. The models which predict most emission at 1 keV are characterized by extreme values of column density, Fe abundance and ionization parameter. Models based on gas in thermal equilibrium with kT approximately l keV provide an alternative parameterization of the soft spectrum. The latter may be interpreted as the hot intercloud medium, undergoing rapid cooling and producing strong Fe L-shell recombination lines. In all cases the physical conditions are rather different from those observed in broad-line Seyferts. The hard X-ray spectrum shows a broad and asymmetric Fe K(alpha) line of large equivalent width (approximately 550 eV) suggestive of significant emission from the inner accretion disk. The profile can be explained by a neutral disk viewed at approximately 60 deg to the line-of-sight, contrary to the hypothesis that NLSyls are viewed pole-on. The large EW of this line, the strong 1 keV emission and the strong optical Fe emission lines all suggest an extreme Fe abundance in this and perhaps other NLSyls.

Turner, T. J.↗

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗

The Isothermal Dendritic Growth Experiment (IDGE)

Dendritic solidification is one of the simplest examples of pattern formation where a structureless melt evolves into a ramified crystalline microstructure; it is a common mode of solidification in many materials, but especially so in metals and alloys. There is considerable engineering interest in dendrites because of the role dendrites play in the determination of microstructure, and thereby in influencing the physical properties of cast metals and alloys. Dendritic solidification provides important examples of non-equilibrium physics, pattern formation dynamics, and models for computational condensed matter and material physics. Current theories of dendritic growth generally couple diffusion effects in the melt with the physics introduced by the interface. Unfortunately, in terrestrial based experiments, convective effects in the melt alter the growth process in such a manner as to prevent definitive analysis of convective, diffusive or interfacial effects. Thus, the effective elimination of convection in the melt by operating experiments on orbit were required to produce high-fidelity data needed for achieving further progress. This simple fact comprised the scientific justification for the IDGE.

Glicksman, M. E.↗

The Diffuse EUV and X-Ray Background as a Probe of the Interstellar Medium

We have used the Deep Survey telescope of the Extreme Ultraviolet Explorer to investigate shadows in the diffuse EUV/soft X-ray background that are cast by clouds in the interstellar medium. We confirm the existence of a shadow previously reported and provide evidence for two new shadows. We used IRAS data to identify the clouds producing these shadows and to determine their optical depth to EUV radiation. The EUV-absorbing clouds are optically thick in the EUV, and all EUV emission detected in the direction of these shadows must be produced from material in front of the clouds. We obtained new optical data to determine the distance to these clouds. We use a new differential cloud technique to obtain the pressure of the interstellar medium. These results do not depend on any zero-level calibration of the data. Our results provide evidence that the pressure of the hot interstellar gas is the same in three different directions in the local interstellar medium and is at least 8 times higher than that derived for the Local Cloud surrounding our Sun. This provides new evidence for large thermal pressure imbalances in the local ISM and directly contradicts the basic assumption of thermal pressure equilibrium used in almost all present models of the interstellar medium.

Lieu, Richard↗

Mass Anomalies on Ganymede

Radio Doppler data from two Ganymede encounters (G1 and G2) on the first two orbits in the Galileo mission have been analyzed previously for gravity information . For a satellite in hydrostatic equilibrium, its gravitational field can be modeled adequately by a truncated spherical harmonic series of degree two. However, a fourth degree field is required in order to fit the second Galileo flyby (G2). This need for a higher degree field strongly suggests that Ganymede s gravitational field is perturbed by a gravity anomaly near the G2 closest approach point (79.29 latitude, 123.68 west longitude). In fact, a plot of the Doppler residuals , after removal of the best-fit model for the zero degree term (GM) and the second degree moments (J2 and C22), suggests that if an anomaly exists, it is located downtrack of the closest approach point, closer to the equator.

Schubert, G.↗

Boundary Layer Protuberance Simulations in Channel Nozzle Arc Jet

Various protuberance heights and shapes were modeled in the channel nozzle of the NASA Johnson Space Center Atmospheric Reentry Materials and Strictures Facility with the Data- Parallel Line Relaxation computational fluid dynamics code. The heating on the protuberance was compared to baseline (no protuberance) heating at a single fixed arc jet condition in order to obtain heating augmentation factors that will be used for flight traceability in the Boundary Layer Transition Flight Experiment on Space Shuttle Orbiter flights STS-119 (completed) and STS-128 (future flight). The arc jet simulations were performed in conjunction with the actual ground tests performed on the flight version (selected height and shape) of the protuberance. Thearc jet simulations for the final (flight version) protuberance included non-uniform inflow conditions beginning at the channel nozzle throat. The 2D inflow condition was modeled based on the current best practices methodology and used variable enthalpy and mass flow rate across the throat. Channel walls were modeled as fully catalytic isothermal surfaces, while the test section (consisting of Reaction Cured Glass tiles) was modeled as a partially catalytic radiative equilibrium wall. The results of the protuberance and baseline simulations were compared to the applicable ground test results. In addition, the obtained heating augmentation factors were compared to the factors derived from the STS-119 flight data. The effects of the protuberance shock on the opposite channel wall were also investigated.

Larin, M. E.↗