Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plasma codes”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

A DSMC Study of Low Pressure Argon Discharge

Work toward a self-consistent plasma simulation using the DSMC (Direct Simulation Monte Carlo) method for examination of the flowfields of low-pressure high density plasma reactors is presented. Presently, DSMC simulations for these applications involve either treating the electrons as a fluid or imposing experimentally determined values for the electron number density profile. In either approach, the electrons themselves are not physically simulated. Self-consistent plasma DSMC simulations have been conducted for aerospace applications but at a severe computational cost due in part to the scalar architectures on which the codes were employed. The present work attempts to conduct such simulations at a more reasonable cost using a plasma version of the object-oriented parallel Cornell DSMC code, MONACO, on an IBM SP-2. Due to availability of experimental data, the GEC reference cell is chosen to conduct preliminary investigations. An argon discharge is chosen to conduct preliminary investigations. An argon discharge is examined thus affording a simple chemistry set with eight gas-phase reactions and five species: Ar, Ar(+), Ar(*), Ar(sub 2), and e where Ar(*) is a metastable.

Hash, David B.↗
NASA NTRS1981

Simulation of charge exchange plasma propagation near an ion thruster propelled spacecraft

A model describing the charge exchange plasma and its propagation is discussed, along with a computer code based on the model. The geometry of an idealized spacecraft having an ion thruster is outlined, with attention given to the assumptions used in modeling the ion beam. Also presented is the distribution function describing charge exchange production. The barometric equation is used in relating the variation in plasma potential to the variation in plasma density. The numerical methods and approximations employed in the calculations are discussed, and comparisons are made between the computer simulation and experimental data. An analytical solution of a simple configuration is also used in verifying the model.

Robinson, R. S.↗
NASA NTRS2002

2-D Magnetohydrodynamic Modeling of A Pulsed Plasma Thruster

Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) MK-1 pulsed plasma thruster. Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.

Thio, Y. C. Francis↗
NASA NTRS2002

An Experimental Study of a Pulsed Electromagnetic Plasma Accelerator

Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) pulsed electromagnetic plasma accelerator (PEPA-0). Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.

Thio, Y. C. Francis↗
NASA NTRS1984

Plasma and magnetospheric research

Research developments in the following areas are discussed: (1) an ion trajectory computer code which plots the paths of ions ejected from the polar cusp ionosphere; (2) the response of plasmaspheric ion temperatures to geomagnetic activity; (3) spacecraft sheath effects; (4) plasma flow; (5) neutral gas temperatures; and (6) instrument fabrication, modification, and maintenance.

Comfort, R. H.↗
NASA NTRS1976

Energy distributions and radiation transport in uranium plasmas

An approximate analytic model, based on continuous electron slowing, has been used for survey calculations. Where more accuracy is required, a Monte Carlo technique is used which combines an analytic representation of Coulombic collisions with a random walk treatment of inelastic collisions. The calculated electron distributions have been incorporated into another code that evaluates both the excited atomic state densities within the plasma and the radiative flux emitted from the plasma.

Miley, G. H.↗
NASA NTRS1990

Current collection in a flowing magnetoplasma

Effects of plasma drift on the current collection by a long conducting cylinder in a magnetized plasma is studied by means of a 2 1/2 dimensional PIC code. It is found that for the drift velocity (V sub 0) perpendicular to the magnetic field (B sub 0), the electron current collected by a positive cylinder is considerably enhanced depending on the drift velocity. The distributions of plasma and the potential structure around the cylinder for several relative orientations between V sub 0 and B are presented along with the comparisons of current with and without the magnetic field. Simulations with the magnetic field in the simulation plane show that the potential structures around the cylinder are two-dimensional double layers with dimension (L) perpendicular to B much smaller than the dimension (L) parallel to B. In fact, L perpendicular is found to be approximately determined by the current limiting radius given by the Parker-Murphy model. However, it is found that the collected currents in the simulations are generally higher than those given by this model.

Singh, Nagendra↗
NASA NTRS1990

SEPAC data analysis in support of the environmental interaction program

Injections of nonrelativistic electron beams from an isolated equipotential conductor into a uniform background of plasma and neutral gas were simulated using a two dimensional electrostatic particle code. The ionization effects of spacecraft charging are examined by including interactions of electrons with neutral gas. The simulations show that the conductor charging potential decreases with increasing neutral background density due to the production of secondary electrons near the conductor surface. In the spacecraft wake, the background electrons accelerated towards the charged space craft produced an enhancement of secondary electrons and ions. Simulations run for longer times indicate that the spacecraft potential is further reduced and short wavelength beam-plasma oscillations appear. The results are applied to explain the space craft charging potential measured during the SEPAC experiments from Spacelab 1. A second paper is presented in which a two dimensional electrostatic particle code was used to study the beam radial expansion of a nonrelativistic electron beam injected from an isolated equipotential conductor into a background plasma. The simulations indicate that the beam radius is generally proportional to the beam electron gyroradius when the conductor is charged to a large potential. The simulations also suggest that the charge buildup at the beam stagnation point causes the beam radial expansion. From a survey of the simulation results, it is found that the ratio of the beam radius to the beam electron gyroradius increases with the square root of beam density and decreases inversely with beam injection velocity. This dependence is explained in terms of the ratio of the beam electron Debye length to the ambient electron Debye length. These results are most applicable to the SEPAC electron beam injection experiments from Spacelab 1, where high charging potential was observed.

Lin, Chin S.↗
NASA NTRS1989

Simulation of electrostatic ion instabilities in the presence of parallel currents and transverse electric fields

A spatially two-dimensional electrostatic PIC simulation code was used to study the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a field-aligned drift for L is much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field aligned drift and the E x B drift are individually subcritical. The simulation results show that the growing ion waves are associated with small vortices in the linear stage, which evolve to the nonlinear stage dominated by larger vortices with lower frequencies.

Nishikawa, K.-I.↗
NASA NTRS1990

Simulation of electrostatic turbulence due to sheared flows parallel and transverse to the magnetic field

A spatially two-dimensional electrostatic particle simulation code is used to examine the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a magnetic-field-aligned electron drift for L much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field-aligned drift and the E x B drift are individually subcritical. The simulation results indicate that a broadband turbulence is associated with such an equilibrium.

Nishikawa, K.-I.↗
NASA NTRS1991

LLNL contributions to MPD thrusters for SEI

Some of the topics covered with respect to the Lawrence Livermore National Laboratory's (LLNL's) contributions to Magnetoplasmadynamic (MPD) Thrusters for the Space Exploration Initiative (SEI) include: an IR camera, plasma-induced erosion/redeposition, the Mirror Fusion Test Facility-B (MFTF-B), the Thruster Lifetime Test Facility, the RACE Compact Torus Accelerator Facility, and a RACE program summary. Some of the other topics addressed include: flux contours for HAM simulation, comparison of RACE data of plasma ring formation with the HAM 2-D magnetohydrodynamic code, and the 2-D Ring Acceleration Code (TRAC).

Hooper, Edwin Bickford↗
NASA NTRS1994

Nonlinear evolution of a large-amplitude circularly polarized Alfven wave: High beta

The nonlinear dynamics following saturation of the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation in the case of high plasma beta and no wave dispersion. The magnetohydrodynamic (MHD) code permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Turbulent cascades develop after saturation of two coupled oblique three-wave parametric instabilities; one of which is an oblique filamentationlike instability reported earlier. Remnants of the parametric processes, as well as of the original Alfven pump wave, persist during late nonlinear times. Nearly incompressible MHD features such as spectral anisotropies appear as well.

Ghosh, S.↗
NASA NTRS1995

Measurement of photoemission and secondary emission from laboratory dust grains

The overall goal of this project is experimentally determine the emission properties of dust grains in order to provide theorists and modelers with an accurate data base to use in codes that predict the charging of grains in various plasma environments encountered in the magnetospheres of the planets. In general these modelers use values which have been measured on planar, bulk samples of the materials in question. The large enhancements expected due to the small size of grains can have a dramatic impact upon the predictions and the ultimate utility of these predictions. The first experimental measurement of energy resolved profiles of the secondary electron emission coefficient, 6, of sub-micron diameter particles has been accomplished. Bismuth particles in the size range of .022 to .165 micrometers were generated in a moderate pressure vacuum oven (average size is a function of oven temperature and pressure) and introduced into a high vacuum chamber where they interacted with a high energy electron beam (0.4 to 20 keV). Large enhancements in emission were observed with a peak value, delta(sub max) = 4. 5 measured for the ensemble of particles with a mean size of .022 micrometers. This is in contrast to the published value, delta(sub max) = 1.2, for bulk bismuth. The observed profiles are in general agreement with recent theoretical predictions made by Chow et al. at UCSD.

Hazelton, Robert C.↗
NASA NTRS2000

Charged Particle Environment Definition for NGST: L2 Plasma Environment Statistics

The plasma environment encountered by the Next Generation Space Telescope satellite in a halo orbit about L2 can include the Earth's magnetotail and magnetosheath in addition to the solar wind depending on the orbital radius chosen for the mission. Analysis of plasma environment impacts on the satellite requires knowledge of the average and extreme plasma characteristics to assess the magnitude of spacecraft charging and materials degradation expected for the mission lifetime. This report describes the analysis of plasma data from instruments onboard the IMP 8 and Geotail spacecraft used to produce the plasma database for the LRAD engineering-level phenomenology code developed to provide the NGST L2 environment definition.

Minow, Joseph I.↗
NASA NTRS2003

RF Manipulation and Detection of Protons in the High Performance Antiproton Trap (HiPAT) Experiment

HiPAT is a Penning-Malmberg trap designed to contain 10(exp 12) particles diagnostics of the plasma. Destructive particle detection is performed by extracting trapped ions from HiPAT and colliding them with a microchannel plate detector (providing number and energy information). The improved RF system has been used to detect various plasma modes for both electron and ion plasmas in the two traps at MSFC, including axd, cyclotron, and diocotron modes. New diagnostics are also being added to HiPAT to measure the axial density distribution of the trapped cloud to match measured RF plasma modes to plasma conditions. In addition, plasma modeling efforts have been started using the XOOPIC code.

Martin, James↗
NASA NTRS2012

Simulation Tool for Dielectric Barrier Discharge Plasma Actuators at Atmospheric and Sub-Atmospheric Pressures: SBIR Phase I Final Report

This report is the final report of a SBIR Phase I project. It is identical to the final report submitted, after some proprietary information of administrative nature has been removed. The development of a numerical simulation tool for dielectric barrier discharge (DBD) plasma actuator is reported. The objectives of the project were to analyze and predict DBD operation at wide range of ambient gas pressures. It overcomes the limitations of traditional DBD codes which are limited to low-speed applications and have weak prediction capabilities. The software tool allows DBD actuator analysis and prediction for subsonic to hypersonic flow regime. The simulation tool is based on the VORPAL code developed by Tech-X Corporation. VORPAL's capability of modeling DBD plasma actuator at low pressures (0.1 to 10 torr) using kinetic plasma modeling approach, and at moderate to atmospheric pressures (1 to 10 atm) using hydrodynamic plasma modeling approach, were demonstrated. In addition, results of experiments with pulsed+bias DBD configuration that were performed for validation purposes are reported.

Likhanskii, Alexandre↗
NASA NTRS1994

Target compressions by working fluids driven with solid liner implosions

Compression by a spherical solid liner of a gold target surrounded by a hydrogen plasma is simulated. Two-dimensional simulations that treat only a subset of the physics included in the one-dimensional code were performed in an attempt to assess multidimensional effects. A one-dimensional numerical code has been developed to study the effects of thermal radiation and conduction. Results of pressure, density, and energy deposited for different initial plasma conditions are presented and discussed. Results from both one- and two-dimensional codes show that the average target density at peak compression is 39-43 g/cu cm, using the SHIVA Star facility at 90 kV discharge.

Chiang, P.-R.↗