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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

PLASIM: A computer code for simulating charge exchange plasma propagation

The propagation of the charge exchange plasma for an electrostatic ion thruster is crucial in determining the interaction of that plasma with the associated spacecraft. A model that describes this plasma and its propagation is described, together with a computer code based on this model. The structure and calling sequence of the code, named PLASIM, is described. An explanation of the program's input and output is included, together with samples of both. The code is written in ANSI Standard FORTRAN.

Robinson, R. S.↗
NASA NTRS2004

Conductance Effects on Inner Magnetospheric Plasma Morphology: Model Comparisons with IMAGE EUV, MENA, and HENA Data

Modeling results of the inner magnetosphere showing the influence of the ionospheric conductance on the inner magnetospheric electric fields during the April 17, 2002 magnetic storm are presented. Kinetic plasma transport code results are analyzed in combination with observations of the inner magnetospheric plasma populations, in particular those from the IMAGE satellite. Qualitative and quantitative comparisons are made with the observations from EW, MENA, and HENA, covering the entire energy range simulated by the model (0 to 300 keV). The electric field description, and in particular the ionospheric conductance, is the only variable between the simulations. Results from the data-model comparisons are discussed, detailing the strengths and weaknesses of each conductance choice for each energy channel.

Liemohn, M.↗
NASA NTRS1995

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗
NASA NTRS2005

Plasma Heating Simulation in the VASIMR System

The paper describes the recent development in the simulation of the ion-cyclotron acceleration of the plasma in the VASIMR experiment. The modeling is done using an improved EMIR code for RF field calculation together with particle trajectory code for plasma transport calculat ion. The simulation results correlate with experimental data on the p lasma loading and predict higher ICRH performance for a higher density plasma target. These simulations assist in optimizing the ICRF anten na so as to achieve higher VASIMR efficiency.

Ilin, Andrew V.↗
NASA NTRS2014

Grid-Free 2D Plasma Simulations of the Complex Interaction Between the Solar Wind and Small, Near-Earth Asteroids

We present results from a new grid-free 2D plasma simulation code applied to a small, unmagnetized body immersed in the streaming solar wind plasma. The body was purposely modeled as an irregular shape in order to examine photoemission and solar wind plasma flow in high detail on the dayside, night-side, terminator and surface-depressed 'pocket' regions. Our objective is to examine the overall morphology of the various plasma interaction regions that form around a small body like a small near-Earth asteroid (NEA). We find that the object obstructs the solar wind flow and creates a trailing wake region downstream, which involves the interplay between surface charging and ambipolar plasma expansion. Photoemission is modeled as a steady outflow of electrons from illuminated portions of the surface, and under direct illumination the surface forms a non-monotonic or ''double-sheath'' electric potential upstream of the body, which is important for understanding trajectories and equilibria of lofted dust grains in the presence of a complex asteroid geometry. The largest electric fields are found at the terminators, where ambipolar plasma expansion in the body-sized night-side wake merges seamlessly with the thin photoelectric sheath on the dayside. The pocket regions are found to be especially complex, with nearby sunlit regions of positive potential electrically connected to unlit negative potentials and forming adjacent natural electric dipoles. For objects near the surface, we find electrical dissipation times (through collection of local environmental solar wind currents) that vary over at least 5 orders of magnitude: from 39 Micro(s) inside the near-surface photoelectron cloud under direct sunlight to less than 1 s inside the particle-depleted night-side wake and shadowed pocket regions

Asteroids↗
NASA NTRS2003

Magnetized Plasma Expansion and its Interaction with a Plasma Stream

Expansion of magnetized plasma in the magnetic field of a solenoid is studied by means of simulations using a 3-dimensional hybrid code. The plasma expands against a high- density and slow plasma stream (PS). The expansion causes inflation of the magnetic field; near the solenoid the magnetic field variation with increasing distance (r) remains as B(alpha)r(sup -3), but at farther distances B(alpha)r(sup -p), where the exponent p is found in the range 0.5 approx. less than p approx. less than 1.2 forming a plateau in the magnetic field distribution B(r). At the start of injection of plasma from the ends of the solenoid, the PS interacts with the solenoid magnetic field and creates a bow shock at a distance where the Larmor radius (r(sub il)) of the PS ions in the solenoid magnetic field nearly equals the scale length (L) of B(r), that is, r(sub il) approx. L = (B(sup -1)dB L /dr)(sup -1). As the injected plasma accumulates in the solenoid magnetic field, it expands inflating the magnetic field and the bow shock moves outward. The speed of the expansion front and the shock progressively decreases and eventually a stand-off occurs when the PS dynamic pressure is eventually balanced by the magnetic and kinetic pressures of the expanding plasma and the inflating magnetic field. The inflating field shows wave-like behavior, with considerable structures in the field and current distributions.

Singh, Nagendra↗
NASA NTRS2003

Spectroscopic Determination of the Physical Conditions in Hot Optically Thin Sources

The Astrophysics Plasma Emission Code and Database (APEC/APED), developed in part under this grant, have been upgraded to Version 1.3.0 and are now beginning to find widespread applications to X-ray spectral data from Chandra and XMM-Newton. These models represent the best theoretical data currently available and are reasonably complete below about 30 A. We are making significant progress in improving the spectral models between 30 and 90 A, by the addition of HULLAC calculations and NIST wavelengths for L-shell ions of Ne, Mg, Si, S, Ar, and Ca. While the output models are now public, the code itself is not ready for public release. We have made significant progress toward this goal through programming improvements, including optimized error codes, a major (elusive) bug fix, and additions to our testing protocols. Secondarily, we are beginning to investigate limited uses of APEC for X-ray photoionized plasma, as for opacity modeling. Stellar coronae are being used to benchmark the atomic data in APED as part of the Emission Line Project. The models appear to be in good agreement with the observations for most of the strong lines; however, we have identified significant discrepancies in the 3s/3d line ratios not only for Fe XVII, but also for Fe XVIII and XIX. The careful analysis of these high quality X-ray spectra has also shown that there are serious problems with some of the analysis techniques in widespread use. There are several manifestations of the problem, but the main problem is that line-to-continuum ratios may be seriously mismeasured. We continue to collaborate widely with atomic physicists in order to find the best atomic data, and to solicit calculations and measurements that are needed.

Brickhouse, Nancy↗
NASA NTRS2013

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗
NASA NTRS2004

Spectroscopic Determination of the Physical Conditions in Hot Optically Thin Sources

Astrophysical Plasma Emission Database - APED consists of atomic data, primarily theoretical, needed to calculate X-ray through optical emission spectra of hot thermal plasmas. These data are supplemented by experimental values, such as laboratory plasma wave- lengths, and are validated by experiment only rarely. Thus the comparison of predicted spectra to astrophysical spectra serves as invaluable feedback on the quality of the models. Following up on unfavorable comparisons, we have communicated atomic data needs to atomic theorists and experimental plasma physicists so that improvements can be made to the database. The database is now reasonably complete from approx. 1 to 30 A, and from - 90 to 180 A, with significant gaps remaining in the range 30 to 90 A. Astrophysical Plasma Emission Code - APEC uses the APED data as inputs and computes level populations and line emissivities using a rate matrix solver. Tabulated models for a fine grid of temperatures and densities can be applied to observational data. The current public version of APEC output uses an input ionization balance model, assuming CIE; however, we have a working version to calculate the ionization state of the gas. We will begin testing non-equilibrium ionization (NEI) models over the next year - we are now finalizing the atomic database needed to perform the ionization state calculations. At that time, the code will be made public. We have developed optimized error codes, additional testing protocols and substantial documentation, in preparation for this public release. Secondarily, we have begun to investigate the use of APED and APEC for X-ray photoionized plasma, beginning with opacity modeling (e.g. warm absorbers). Application to bright AGN Chandra grating spectra suggests that the accurate wavelengths in APED are of great importance to fitting the data. We are exploring additional subroutines to APEC to extend its usefulness for X-ray photoionized plasma, e.g. to include collisional or photo-excitation of metastable levels for absorption.

Oliversen, Ronald J.↗
NASA NTRS1990

X-ray spectral signatures of photoionized plasmas

Plasma emission codes have become a standard tool for the analysis of spectroscopic data from cosmic X-ray sources. However, the assumption of collisional equilibrium, typically invoked in these codes, renders them inapplicable to many important astrophysical situations, particularly those involving X-ray photoionized nebulae. This point is illustrated by comparing model spectra which have been calculated under conditions appropriate to both coronal plasmas and X-ray photoionized plasmas. It is shown that the (3s-2p)/(3d-2p) line ratios in the Fe L-shell spectrum can be used to effectively discriminate between these two cases. This diagnostic will be especially useful for data analysis associated with AXAF and XMM, which will carry spectroscopic instrumentation with sufficient sensitivity and resolution to identify X-ray photoionized nebulae in a wide range of astrophysical environments.

Liedahl, Duane A.↗
NASA NTRS2003

Spectroscopic Determination of the Physical Conditions in Hot Optically Thin Sources

The Astrophysics Plasma Emission Code and Database (APEC/APED), developed in part under this grant, have been upgraded to ATOMDB Version 1.3.1: and are now beginning to find widespread applications t o X-ray spectral data from Chandra and XMM-Newton (37 citations in published work, according to the ADS, plus numerous other conference and prepublication papers).ATOMDB is now linked through the Plasma Gate website: http://plasma-gate.weizmann.ac.il/DBfAPP.html. The major difference from Version 1.3.0 is that the new models naw ex- tend t o 50 keV rather than stopping at 10 keV. This means that ATOMDB can be used with redshifted observations. There are minor differences in emissivities due t o radiative recombination and cascades. Stellar coronae are being used to benchmark the atomic data in APED as part of the Emission Line Project. The models appear to be in good agreement with the observations for most of the strong lines; however, we have identified significant discrepancies in the 3s/3d line ratios not only for Fe XVII, but also for Fe XVIII and XIX. The Fe XVII problem has been known from solar observations, and is currently being tested under EBIT laboratory conditions by two groups. The Fe XVIII problem is substantially worse, but perhaps will shed light on the relevant underlying theoretical issues. Ming-Feng Gu has recently published new calculations, which we are comparing with APEC and with the obser- vations. His calculations appear to improve the emissivities of lines affected by cascades, but other problems remain.

Brickhouse, Nancy S.↗
NASA NTRS1998

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗
NASA NTRS1996

Effects of plasmaspheric ion heating due to ionospheric and magnetospheric sources

In an initial study, the He(+) observations from the Retarding Ion Mass Spectrometer on Dynamics Explorer 1 (RIMS/DE 1) was examined for more than 120 transits of the plasmasphere in the fall of 1981. The He(+) to H(+) ratio was determined as it varied spatially over portions of the DE 1 orbit, and its variation with solar and magnetic activities and with local time, focusing specifically on the inner plasmasphere. These variations were compared along the L = 2 field line with calculations made by the Field Line Interhemispheric Plasma (FLIP) code. In a recently submitted paper, the He(+) to H(+) density ratio was examined for all the available data from 1981 to 1984 from the RIMS on DE 1. There are two basic characteristics of the ratio: one is that the ratio decreases with radial distance in the plasmasphere, and the other is the strong dependence of the density ratio on solar activity. In addition to the He(+)/H(+) ratio research, a phenomenon has been studied in the topside ionosphere which relates to the thermal coupling of the ionosphere to the plasmasphere. There is little or no correlation with magnetic and solar activity here. Another study has been directed toward the relation of plasma properties to the density gradients forming the plasmapause. The study has followed a two-pronged approach. First, the observations have been analyzed to determine what happens to the plasma properties across these boundary layers (density gradients). Second, comparisons were made with FLIP model calculations to determine how well the model is able to treat these conditions. Among the significant lessons learned in these studies are two that bear directly on the direction of future investigations in this area. First, composition cannot be viewed independently of thermal structure. Second, solar and magnetic activity effects are real; but the causal relationship between activity and effects is frequently quite complicated because several different processes appear to be operating in different ways and on different time scales. Under these circumstances, large correlation coefficients should not be expected and are not generally found.

Comfort, Richard H.↗
NASA NTRS1986

Predicting the magnetospheric plasma of weather

The prediction of the plasma environment in time, the plasma weather, is discussed. It is important to be able to predict when large magnetic storms will produce auroras, which will affect the space station operating in low orbit, and what precautions to take both for personnel and sensitive control (computer) equipment onboard. It is also important to start to establish a set of plasma weather records and a record of the ability to predict this weather. A successful forecasting system requires a set of satellite weather stations to provide data from which predictions can be made and a set of plasma weather codes capable of accurately forecasting the status of the Earth's magnetosphere. A numerical magnetohydrodynamic fluid model which is used to model the flow in the magnetosphere, the currents flowing into and out of the auroral regions, the magnetopause, the bow shock location and the magnetotail of the Earth is discussed.

Dawson, John M.↗
NASA NTRS1989

Development of a numerical solution to the time dependent kinetic equation

A numerical solution was developed for the time dependent Fokker-Planck equation for arbitrary distributions of electrons injected into a magnetized plasma. The code which includes energy loss and pitch angle scattering due to Coulomb collisions and changes in pitch angle due to inhomogeneous magnetic fields was calibrated and tested. The numerical method is versatile so that other scattering or radiation terms can be easily included. Using this code many processes associated with the impulsive phase of solar flares will be investigated.

Hamilton, Russell J.↗
NASA NTRS1991

Numerical studies of electron dynamics in oblique quasi-perpendicular collisionless shock waves

Linear and nonlinear electron damping of the whistler precursor wave train to low Mach number quasi-perpendicular oblique shocks is studied using a one-dimensional electromagnetic plasma simulation code with particle electrons and ions. In some parameter regimes, electrons are observed to trap along the magnetic field lines in the potential of the whistler precursor wave train. This trapping can lead to significant electron heating in front of the shock for low beta(e). Use of a 64-processor hypercube concurrent computer has enabled long runs using realistic mass ratios in the full particle in-cell code and thus simulate shock parameter regimes and phenomena not previously studied numerically.

Liewer, P. C.↗
NASA NTRS2001

Spectroscopic Determination of the Physical Conditions in Hot, Optically Thin Sources

This annual report covers the period from Oct 1, 2000 to Sep 30, 2001. The Astrophysics Plasma Emission Code and Database (APEC/APED), developed in part under this grant, have been upgraded to Version 1.1 and are now beginning to be used outside our research group in applications to X-ray spectral data from Chandra and XMM-Newton. These models represent the best theoretical data currently available and are reasonably complete below about 30 angstroms. Stellar coronae are being used to benchmark the atomic data in APED as part of the Emission Line Project. Initial results suggest that the models for most of the strongest lines are in good agreement with the observations, including the H-like and He-like emission from dominant elements and the Fe L shell emission, both near 1 keV and in the extreme ultraviolet near 100 eV. At this point in time, we define 'good agreement' at the level of accuracy expected from the atomic physics, approximately 20-30%. In order to benchmark the spectral models beyond the canonical theoretical accuracy, we are working closely with the Chandra gratings calibration group to ensure that we are using the optimal calibration, primarily effective areas and line response functions.

Brickhouse, N.↗