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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

SEPAC data analysis in support of the environmental interaction program

The main objective is to conduct data analyses of SEPAC data and computer modeling to investigate spacecraft environment effects associated with injection of electron beam, plasma clouds, and neutral gas clouds from the Shuttle Orbiter. To understand the dependence of spacecraft charging potential on beam density and other plasma parameters, a two dimensional electrostatic particle code was used to simulate the injection of electron beams from an infinite conductor into a plasma. The ionization effects on spacecraft charging are examined by including interactions of electrons with neutral gases. A survey of the simulation results is presented and discussed.

Lin, C. S.↗
NASA NTRS1983

Potentials in a plasma over a biased pinhole

The NASCAP/LEO code is used to simulate measurements taken at Jet Propulsion Laboratory of potentials near a simulated pinhole. The insulator near the high-voltage pinhole obeys an electric field boundary condition resulting from secondary electron hopping conductivity. The code predictions are in good agreement with the measurements.

Mandell, M. J.↗
NASA NTRS1993

Strongly nonlinear evolution of low-frequency wave packets in a dispersive plasma

The evolution of strongly nonlinear, strongly modulated wave packets is investigated in a dispersive plasma using a hybrid numerical code. These wave packets have amplitudes exceeding the strength of the external magnetic field, along which they propagate. Alfven (left helicity) wave packets show strong steepening for p < 1, while fast (fight heIicity) wave packets hardly steepen for any beta. Substantial regions of opposite helicity form on the leading side of steepened Alfven wave packets. This behavior differs qualitatively from that exhibited by the solutions to the derivative nonlinear Schrodinger (DNLS) equation.

Vasquez, Bernard J.↗
NASA NTRS2006

Towards the Identification of the Keeper Erosion Cause(s): Numerical Simulations of the Plasma and Neutral Gas Using the Global Cathode Model OrCa2D-II

Numerical simulations with the time-dependent Orificed Cathode (OrCa2D-II) computer code show that classical enhancements of the plasma resistivity can not account for the elevated electron temperatures and steep plasma potential gradients measured in the plume of a 25-27.5 A discharge hollow cathode. The cathode, which employs a 0.11-in diameter orifice, was operated at 5.5 sccm without an applied magnetic field using two different anode geometries. It is found that anomalous resistivity based on electron-driven instabilities improves the comparison between theory and experiment. It is also estimated that other effects such as the Hall-effect from the self-induced magnetic field, not presently included in OrCa2D-II, may contribute to the constriction of the current density streamlines thus explaining the higher plasma densities observed along the centerline.

plasma modeling↗
NASA NTRS2013

MHD Simulations of the Plasma Flow in the Magnetic Nozzle

The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the magnetic field in both the near- and far-field regions. However, in that work the downstream domain was constrained to a channel of constant cross-sectional area. In the present work we seek to address this issue by modeling the downstream region with a domain that permits free expansion of the plasma, permitting a better evaluation of the downstream effects the applied field has on the plasma. The inlet boundary conditions and applied magnetic field values will also be varied to determine the effect the initial plasma energy content and applied magnetic field energy density have on the near- and far-field plasma properties on the MHD code. This will determine the effect of inlet boundary conditions on the results downstream and address issues related to the restrictive numerical domain previously used.

Smith, T. E. R.↗
NASA NTRS1992

Anomalous ion mixing within a Kelvin-Helmholtz vortex in a collisionless plasma

Anomalously fast ion mixing is observed in a hybrid code simulation (particle ions and a massless electron fluid) of the Kelvin-Helmholtz instability in a collisionless plasma. While the traditional view predicts that the ion mixing occurs on a time scale longer than the roll-up time T(v) of the Kelvin-Helmholtz vortex by a factor of (velocity shear scale length)/(ion Larmor radius), the simulation results show that the mixing occurs within a time scale comparable to T(v). It is shown that this fast mixing is due to the scattering of ions by unsteady electromagnetic fields.

Terasawa, T.↗
NASA NTRS1987

Plasma simulation using the massively parallel processor

Two dimensional electrostatic simulation codes using the particle-in-cell model are developed on the Massively Parallel Processor (MPP). The conventional plasma simulation procedure that computes electric fields at particle positions by means of a gridded system is found inefficient on the MPP. The MPP simulation code is thus based on the gridless system in which particles are assigned to processing elements and electric fields are computed directly via Discrete Fourier Transform. Currently, the gridless model on the MPP in two dimensions is about nine times slower that the gridded system on the CRAY X-MP without considering I/O time. However, the gridless system on the MPP can be improved by incorporating a faster I/O between the staging memory and Array Unit and a more efficient procedure for taking floating point sums over processing elements. The initial results suggest that the parallel processors have the potential for performing large scale plasma simulations.

Lin, C. S.↗
NASA NTRS1990

Concurrent Algorithm For Particle-In-Cell Simulations

Separate decompositions used for particle-motion and field calculations. General Concurrent Particle-in-Cell (GCPIC) algorithm used to implement motions of individual plasma particles (ions and electrons) under influence of particle-in-cell (PIC) computer codes on concurrent processors. Simulates motions of individual plasma particles under influence of electromagnetic fields generated by particles themselves. Performed to study variety of nonlinear problems in plasma physics, including magnetic and inertial fusion, plasmas in outer space, propagation of electron and ion beams, free-electron lasers, and particle accelerators.

Liewer, Paulett C.↗
NASA NTRS2022

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗
NASA NTRS1988

Evolution of photon and particle spectra in compact, luminous objects

The physics of high energy photons and particles (especially electrons and positrons) in the compact, high-energy-density of galactic nuclei and quasars was investigated. A numerical code was developed which follows the nonlinear spectral evolution of a pair/photon plasma, due to two-body scattering and interaction process, in an unmagnetized system. The code was applied both to static plasmas and to relativistic expanding winds.

Eilek, Jean A.↗
NASA NTRS1987

A magnetohydrodynamic model of the interaction of the solar wind with the Jovian magnetosphere and a magnetohydrodynamic simulation of the interaction of the solar wind with the out flowing plasma from a comet

A three-dimensional code for a rapidly rotating magnetosphere in which the MHD equations and the Maxwell equations were solved by using the two step Lax Endroff scheme, was developed. Preliminary results were presented at the Fall AGU meeting in San Francisco. The basic simulation model to study the solar wind interactions was adapted to other bodies in addition to Jupiter. Because of the recent comet flybys, a comet was chosen as the first model. The aim was to model the formation of the contact surface and the plasma tail. Later, work was begun on a three-dimensional model which would include the effects of mass loading. This model was designed to study the weak cometary bow shocks observed by the probes to comets Halley and Giacobini-Zinner. The model was successful in reproducing the position and shape of the bow shock which was determined by using observations from the Suisei spacecraft.

Walker, R. J.↗
NASA NTRS1989

Ion beam generation at the plasma sheet boundary layer by kinetic Alfven waves

A two-dimensional quasi-linear numerical code was developed for studying ion beam generation at the plasma sheet boundary layer by kinetic Alfven waves. The model assumes that the central plasma sheet is the particle source, and that the last magnetic field lines on which kinetic Alfven waves exist and diffusion occurs can be either open or closed. As the possible source for the excitement of the kinetic Alfven waves responsible for ion diffusion, the resonant mode conversion of the surface waves to kinetic Alfven waves is considered. It is shown that, depending on the topology of the magnetic field at the lobe side of the simulation system, i.e., on whether field lines are open or closed, the ion distribution function may or may not reach a steady state.

Moghaddam-Taaheri, E.↗
NASA NTRS1991

SEPAC data analysis in support of the environmental interaction program

Data analyses of the Space Experiments with Particle Accelerators (SEPAC) data and computer modeling were conducted to investigate spacecraft environmental effects associated with injection of electron beams, plasma clouds, and neutral gas clouds from the Shuttle orbiter. The data analysis indicates that Extremely Low Frequency oscillations from 150 to 200 Hz were seen in the Langmuir probe current when the beam was fired in a continuous mode. The strongest oscillations occurred when the ambient pressure was augmented by neutral gas releases from the SEPAC plasma accelerator magnetoplasma-dynamic (MPD) arcjet. To understand the dependence of spacecraft charging potential on beam density and other plasma parameters, a two-dimensional electrostatic particle code was used to simulate the injection of electron beams from an infinite conductor into a plasma. The simulations show that the conductor charging potential depends critically on the reflection coefficient of the conductor surface, which is defined as the percentage of incident particles reflected by the conductor. The ionization effects on spacecraft charging were 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. The simulations also indicate that the beam radius is generally proportional to the beam electron gyroradius when the conductor is charged to a large potential. It appears that the charge buildup at the beam stagnation point causes the beam radial expansion. A survey of the simulation results suggests 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. These results are useful for explaining the spacecraft charging phenomena observed during SEPAC experiments from Spacelab 1.

Lin, Chin S.↗
NASA NTRS1994

Hybrid simulations of whistler waves generation and current closure by a pulsed tether in the ionosphere

The dynamic response of a magnetized collisionless plasma to an externally driven, finite size, sudden switch-on current source across the magnetic field has been studied using a two dimensional hybrid code. It was found that the predominant plasma response was the excitation of whistler waves and the formation of current closure by induced currents in the plasma. The results show that the current closure path consists of: (a) two antiparallel field-aligned current channels at the end of the imposed current sheet; and (b) a cross-field current region connecting these channels. The formation of the current closure path occured in the whistler timescale much shorter than that of MHD and the closure region expanded continuously in time. The current closure process was accompanied by significant energy loss due to whistler radiation.

Chang, C. L.↗
NASA NTRS1989

AFE dynamic effects in inhomogeneous plasmas

The Microwave Reflectometer Ionization Sensor (MRIS) is an instrument on the Aeroassist Flight Experiment (AFE) satellite which will be deployed from the space shuttle. The flow characteristic around a hypersonic bluff reentry vehicle will be measured by the AFE. The general mission of the MRIS is to measure the electron density within the range from 10(exp12) to 10(exp 15) electrons per cu cm and determine the distance to the location of each measured from the surface of the aerobrake. These measurements will be compared with prior aerothermodynamic computer code predictions. Since a knowledge of plasma dynamic effects is important for MRIS design and post flight analysis, it is of interest to consider any possibility of plasma dynamic effects and especially in inhomogeneous plasmas. Of particular interest is the need to study plasma dynamic effects that may emerge from a flow field stationary state that was determined without regard to electric or magnetic fields.

Juang, Jeng-Nan↗
NASA NTRS1985

Space plasma simulations; Proceedings of the Second International School for Space Simulations, Kapaa, HI, February 4-15, 1985. Parts 1 & 2

Space plasma simulations, observations, and theories are discussed. Papers are presented on the capabilities of various types of simulation codes and simulation models. Consideration is given to plasma waves in the earth's magnetotail, outer planet magnetosphere, geospace, and the auroral and polar cap regions. Topics discussed include space plasma turbulent dissipation, the kinetics of plasma waves, wave-particle interactions, whistler mode propagation, global energy regulation, and auroral arc formation.

Ashour-Abdalla, M.↗
NASA NTRS1989

The POLAR code wake model - Comparison with in situ observations

Measurements made by the plasma diagnostics package during Spacelab 2 are discussed and compared with predictions made with the Air Force Geophysics Laboratory POLAR wake code. This code uses a complex geometric model of the orbiter and the self-similar solution of the expansion of a plasma into a vacuum as its model basis. Excellent qualitative and quantitative agreement is found at distances greater than about 30 m. At least to the first order, the model's approximations are justified.

Murphy, G.↗
NASA NTRS1992

Studies of dynamic processes related to active experiments in space plasmas

This is the final report for grant NAGw-2055, 'Studies of Dynamic Processes Related to Active Experiments in Space Plasmas', covering research performed at the University of Michigan. The grant was awarded to study: (1) theoretical and data analysis of data from the CHARGE-2 rocket experiment (1keV; 1-46 mA electron beam ejections) and the Spacelab-2 shuttle experiment (1keV; 100 mA); (2) studies of the interaction of an electron beam, emitted from an ionospheric platform, with the ambient neutral atmosphere and plasma by means of a newly developed computer simulation model, relating model predictions with CHARGE-2 observations of return currents observed during electron beam emissions; and (3) development of a self-consistent model for the charge distribution on a moving conducting tether in a magnetized plasma and for the potential structure in the plasma surrounding the tether. Our main results include: (1) the computer code developed for the interaction of electrons beams with the neutral atmosphere and plasma is able to model observed return fluxes to the CHARGE-2 sounding rocket payload; and (2) a 3-D electromagnetic and relativistic particle simulation code was developed.

Banks, Peter M.↗