SEARCH · Search NASA
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.
Solution of Boltzmann and rate equations for electron distribution function and state populations in nonequilibrium MHD plasmas
Steady state composition of low density nonequilibrium hydrogen plasma
A code for the simulation of nonequilibrium plasmas is being developed, with the capability to couple the plasma fluid-dynamics for a single fluid with a collisional-radiative model, where electronic states are treated as separate species. The model allows for non-Boltzmann distribution of the electronic states. Deviations from the Boltzmann distributions are expected to occur in the rapidly ionizing regime behind a strong shock or in the recombining regime during a fast expansion. This additional step in modeling complexity is expected to yield more accurate predictions of the nonequilibrium state and the radiation spectrum and intensity. An attempt at extending the code to molecular plasma flows is presented. The numerical techniques used, the thermochemical model, and the results of some numerical tests are described.
Solution of Boltzmann, and rate equations for bound electronic states in nonequilibrium magnetohydrodynamic plasma
A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.
A computer program called NEQRAP is described that calculates the radiative properties of nonequilibrium ionized hydrogen. From the given electron temperature, electron density, and atom density values (which do not necessarily satisfy the equilibrium relationship) and intensities of incident radiation, the non-Boltzmann populations of electronic states are computed by solving the equation of quasi-steady-state population distribution. Emission and absorption coefficients are determined as functions of wavelength by invoking the principle of detailed balance between the upper and lower states of each radiative transition. Radiative transport through the medium is computed assuming a one-dimensional uniform slab. The rate of ionic reaction is also computed. When used on a sample case, the program shows that there is a large difference between the calculated intensities of radiation emitted by a bulk of equilibrium and nonequilibrium hydrogen. The accuracy of the program is estimated to be better than 10%.
Initial results from the first high-resolution study of the X-ray spectrum of a supernova remnant are reported. The spectrum of Puppis A between 500 and 1100 eV has been surveyed with the Focal Plane Crystal Spectrometer on the Einstein Observatory and the flux in eight lines and three line blends from various transitions of highly ionized nitrogen, oxygen, neon, and iron has been measured. The spectrum resembles that of active regions in the solar corona, but the neon lines seem enhanced relative to the oxygen lines, and both are enhanced relative to iron. The observed line strengths serve as diagnostics of the physical conditions in the emitting material. Without a priori assumptions, the ionization temperatures of O (2.2 million K) and Ne (4 million K) and the dominant ionization stage of Fe (Fe XVII) are deduced. Isothermal equilibrium plasmas or homogeneous nonequilibrium plasmas as descriptions of the source are ruled out. It is concluded that Puppis A must contain plasma over a temperature range 2-million to 5-million K with an intervening absorption column of 4 x 10 to the 21st/sq cm.
Continuum intensity for electron-ion recombination in nonequilibrium plasma
Refractive index computation for free electrons and neutral atoms in cesium plasma at helium neon laser wavelengths
Godunov-type numerical solution of the 20 moment plasma transport equations. One of the centerpieces of our proposal was the development of a higher order Godunov-type numerical scheme to solve the gyration dominated 20 moment transport equations. In the first step we explored some fundamental analytic properties of the 20 moment transport equations for a low b plasma, including the eigenvectors and eigenvalues of propagating disturbances. The eigenvalues correspond to wave speeds, while the eigenvectors characterize the transported physical quantities. In this paper we also explored the physically meaningful parameter range of the normalized heat flow components. In the second step a new Godunov scheme type numerical method was developed to solve the coupled set of 20 moment transport equations for a quasineutral single-ion plasma. The numerical method and the first results were presented at several national and international meetings and a paper describing the method has been published in the Journal of Computational Physics. To our knowledge this is the first numerical method which is capable of producing stable time-dependent solutions to the full 20 (or 16) moment set of transport equations, including the full heat flow equation. Previous attempts resulted in unstable (oscillating) solutions of the heat flow equations. Our group invested over two man-years into the development and implementation of the new method. The present model solves the 20 moment transport equations for an ion species and thermal electrons in 8 domain extending from a collision dominated to a collisionless region (200 km to 12,000 km). This model has been applied to study O+ acceleration due to Joule heating in the lower ionosphere.
Nonequilibrium ionization in asymmetric plasma sheath determined for Apollo spacecraft at superorbital velocity reentry
The construction of second-order upwind schemes for nonequilibrium plasmas, for both one- and two-fluid formulations is demonstrated. Coupled relaxation processes, including ionization kinetics and radiative processes and their algorithms for nonequilibrium, multiple temperature conditions are described as well. The paper applies the numerical techniques on some simple test cases, points out critical problems and their solutions, and makes qualitative comparisons with known results, whenever possible.
Drift velocities and thermal flux vectors in five component nonequilibrium plasma in presence of electric and magnetic fields, using Grad method
Solution to modified Chapman-Enskog kinetic equation for calculating electrical conductivity of nonequilibrium plasma
The utilization of a cold or nonequilibrium plasma in chlorination processing is discussed. Titanium dioxide (TiO2) was successfully chlorinated at temperatures between 700 and 900 C without the aid of carbon. In addition to these initial experiments, a technique was developed for determining the temperature of a specimen in a plasma. Development of that technique has required evaluating the emissivity of TiO2, ZrO2, and FeOTiO2 and analyzing the specimen temperature in a plasma as a function of both power absorbed by the plasma and the pressure of the plasma. The mass spectrometer was also calibrated with TiCl4 and CCl4 vapor.
The purpose of this investigation is to examine the O VI emission from the X-ray bright ring of the supernova remnant 1 E0102.2-729 in the small Magellanic cloud. Three pointings were positioned tangent to the ring, north (N), northeast (NE), and southeast (SE), to examine a range of X-ray emitting regions overlapping a range of optical [O III] nebulosity and to examine the velocity structure. One background pointing was also obtained, but it was contaminated by a star. The background levels in the pointings on the remnant were low enough that the background pointing was not required for the remaining analysis. The SE pointing was reobserved in August, 2004, in order to bring the total exposure up to the originally requested 15 ks. The archive notified us of the data's availability in mid September. Significant broad O VI 1032 and O VI 1038 emission was found, brightest in the NE and SE pointings. In the NE and SE pointings, the FWHM of the broad O VI component is approx.800-1000 km/s, while in the N pointing, the line is approx.1500 km/s wide. The O VI is redshifted in the N (approx.380 km/s) and NE (approx.60 km/s) but is blueshifted in the SE (approx. -160 km/s). These FUSE O VI velocity dispersions can be compared to the X-ray gas velocities inferred from Doppler distortions in the Chandra X-ray data as reported by Flanagan et al. 2004 (ApJ 605, 230). The bulk velocities in the X-ray bright ring of order +/- 1000 km/s, comparable to the velocity dispersion seen in the FUSE data. However, the X-ray data indicates a redshift of approx.1000 km/s in the SE, while the FUSE data show a blueshift of approx.160 km/s, underscoring the complex velocity structure in this remnant. The O VI fluxes estimated from the fits to the FUSE data were combined with X-ray (XMM- Newton) O VI1 and O VIII fluxes and compared with predictions from a plasma nonequilibrium ionization model in a "line-based" analysis. We found that the plasma departs significantly from collisional ionization equilibrium, particularly in the SE, and that the plasma excitation conditions vary among the pointings: the O plasma sees different conditions than the Ne and Mg plasmas.
An investigation of the dynamics of shock waves in weakly ionized argon plasmas has been performed using a pressure ruptured shock tube. The velocity of the shock is observed to increase when the shock traverses the plasma. The observed increases cannot be accounted for by thermal effects alone. Possible mechanisms that could explain the anomalous behavior include a vibrational/translational relaxation in the nonequilibrium plasma, electron diffusion across the shock front resulting from high electron mobility, and the propagation of ion-acoustic waves generated at the shock front. Using a turbulence model based on reduced kinetic theory, analysis of the observed results suggest a role for turbulence in anomalous shock dynamics in weakly ionized media and plasma-induced hypersonic drag reduction.
Nonequilibrium excitation in recombining N plasma nozzle flows