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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2-D, two-temperature, single fluid magnetohydrodynamic code which incorporates classical plasma transport coefficients and Hall effects was developed to predict steady-state, self-field magnetoplasmadynamic (MPD) thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects, the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries, assuming a fully ionized argon propellant. A maximum predicted specific impulse of 1680 s is obtained for a thruster with an anode radius of 2.5 cm, a cathode radius of 0.5 cm, and equal electrode lengths of 2.5 cm. A scaling relation is developed to predict, within limits, the onset of cylindrical, self-field thruster instability as a function of geometry and operating condition.

Lapointe, M.↗
NASA NTRS1992

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2D, two-temperature, single fluid MHD code which incorporates classical plasma transport coefficients and Hall effects has been developed to predict steady-state, self-field MPD thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries assuming a fully ionized argon propellant.

Lapointe, M.↗
NASA NTRS1994

A two-dimensional particle-in-cell model of a dusty plasma

Dusty plasmas are present in comets, in the ring systems of the outer planets, and in the interstellar medium. A two-dimensional particle-in-cell (PIC) model of a dusty plasma is presented in this paper. The PIC code is best suited for modeling the plasma-dust interaction for large grains, with diameters of the order of a centimeter. We have modeled the charging process for an individual dust grain and the associated potential pattern in the surrounding plasma. We have also considered the case of a large number of grains in a plasma, with intergrain separations of the order of the Debye length, and have shown that the plasma becomes depleted and the charge on a dust grain is reduced, as other workers in this field have predicted (cf. C. K. Goertz, 1989). We examine the electron and ion distribution functions in the vicinity of a charged grain and demonstrate that the ions near a grain have clearly been accelerated by the electrostatic potential.

Young, B.↗
NASA NTRS1994

Landau damping of magnetospherically reflected whistlers

Unducted VLF signals produced by lightning activity can form a population of magnetospherically reflected (MR) whistlers in the inner magnetosphere. It has been suggested recently that in the absence of significant attenuation such waves could merge into a broadband continuum with sufficient intensity to account for plasmaspheric hiss. To test this conjecture we have evaluated the path-integrated attenuation of MR whistlers along representative ray paths using the HOTRAY code. Using a realistic plasma distribution modeled on in-situ data, we find that the majority of MR waves experience significant damping after a few transits across the equator. This is primarily due to Landau resonance with suprathermal (0.1-1 keV) electrons. The attenuation is most pronounced for waves that propagate through the outer plasmasphere; this can readily account for the infrequent occurrence of multiple-hop MR waves for L greater than or equal to 3.5. Selected waves that originate at intermediate latitudes (15 deg is less than or equal to lambda is less than or equal to 35 deg) and whose ray paths are confined to the inner plasma- sphere may experience up to 10 magnetospheric reflections before substantial attentuation occurs. These waves should form the population of observed MR waves. Wave attenuation becomes more pronounced at higher frequencies; this can account for the absence of multiple-hop waves above 5 kHz. Weakly attenuated MR waves tend to migrate outward to the L shell, where their frequency is comparable to the equatorial lower hybrid frequency. The enhanced concentration of waves due to a merging of ray paths would produce a spectral feature that rises in frequency at lower L. This is quite distinct from the reported properties of plasmaspheric hiss, which maintains a constant frequency band throughout the entire plasmasphere. Furthermore, in the absence of mode conversion, waves below 500 Hz, which often form an important if not dominant part of the spectral properties of hiss, are unable to escape from the topside ionosphere in the whistler mode. Consequently, we conclude that unducted lightning signals cannot account for the origin of plasmaspheric hiss.

Thorne, Richard M.↗
NASA NTRS2020

Atomic Data for Plasma Spectroscopy: The CHIANTI Database, Improvements and Challenges

CHIANTI is an atomic database and software package for modeling emission lines and continua from hot astrophysical plasmas. It is freely available to all researchers and has been widely used in the Heliophysics and Astrophysics communities for almost 25 years. In this review, we summarize the properties of the current version of the database and give an overview of the relevant atomic processes. We also discuss progress towards a complete implementation of collisional-radiative modeling, simultaneously solving for atomic level and ion populations for individual elements.

A&M databases↗
NASA NTRS1992

Ion-cyclotron waves at Jupiter - Possibility of detection by Ulysses

Recent thermal plasma data and a computer code by Horne (1989) are employed to evaluate the linear-path-integrated gain of waves propagating through the Io to predict the Jovian plasma-wave environment. Estimates of the nonlinear saturation amplitudes are utilized with the thermal plasma data from two frequency bands to study the convective growth of the ion-cyclotron (IC) waves. Strong cyclotron resonant damping is theorized to prevent wave propagation to the lower latitudes, and the thermal plasma and cyclotron resonant energetic ions are expected to further confine the IC waves. L-mode waves below the O(+) gyrofrequency in the equatorial region of the torus are shown to inhabit an unstable region. The IC waves probably achieve nonlinear amplitudes regardless of plasma properties due to the rapid amplification in this region. It is suggested that the Ulysses data complicate the identification of the waves because the magnetometer is not adequately sensitive and because of the low frequency of the plasma-wave detector.

Mei, YI↗
NASA NTRS1989

PIC code modeling of spacecraft charging potential during electron beam injection into a background of neutral gas and plasma, part 1

Injections of nonrelativistic electron beams from an isolated equipotential conductor into a uniform background of plasma and neutral gas were simulated using a 2-D electrostatic particle code. The ionization effects on 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 spacecraft produce 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 spacecraft charging potential measured during the SEPAC experiments from Spacelab 1.

Koga, J. K.↗
NASA NTRS2005

Dusty Plasmas in Planetary Magnetospheres Award

This is my final report for the grant Dusty Plasmas in Planetary Magnetospheres. The funding from this grant supported our research on dusty plasmas to study: a) dust plasma interactions in general plasma environments, and b) dusty plasma processes in planetary magnetospheres (Earth, Jupiter and Saturn). We have developed a general purpose transport code in order to follow the spatial and temporal evolution of dust density distributions in magnetized plasma environments. The code allows the central body to be represented by a multipole expansion of its gravitational and magnetic fields. The density and the temperature of the possibly many-component plasma environment can be pre-defined as a function of coordinates and, if necessary, the time as well. The code simultaneously integrates the equations of motion with the equations describing the charging processes. The charging currents are dependent not only on the instantaneous plasma parameters but on the velocity, as well as on the previous charging history of the dust grains.

Horanyi, Mihaly↗
NASA NTRS1987

Simulations of auroral plasma processes - Electric fields, waves and particles

Plasma processes driven by current sheets of finite thicknesses in an ambient magnetized plasma are studied using a 2 1/2 dimensional particle-in-cell code, and similarities are found between simulated plasma processes and those observed in the auroral plasma. Current sheets are shown to be bounded by large perpendicular electric fields occurring near their edges above the conducting boundary. Shaped potential structures form when the current sheets are narrow, and when the current sheets are wide, potential structures develop a significant parallel potential drop such that the electrons are accelerated upwards. Downward parallel electric fields of variable strength are noted in the downward current region, and double layer formation is seen in both narrow and wide current sheets. High frequency oscillations near the electron plasma frequency and its harmonic are seen, and low frequency waves are observed.

Singh, Nagendra↗
NASA NTRS2020

A Double Disturbed Lunar Plasma Wake

Under nominal solar wind conditions, a tenuous wake forms downstream of the lunar nightside. However, the lunar plasma environment undergoes a transformation as the Moon passes through the Earth's magnetotail, with hot subsonic plasma causing the wake structure to disappear. We investigate the lunar wake response during a passing coronal mass ejection (CME) on March 8, 2012 while crossing the Earth's magnetotail using both a magnetohydrodynamic (MHD) model of the terrestrial magnetosphere and a three-dimensional hybrid plasma model of the lunar wake. The CME arrives at 1 AU around 10:30 UT and its impact is first detected inside the geomagnetic tail after 11:10 UT by the Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun (THEMIS-ARTEMIS) satellites in lunar orbit. A global magnetospheric MHD simulation using Wind data for upstream conditions with the OpenGGCM model reveals the magnetosheath compression to the lunar position from 11:20–12:00 UT, accompanied by multiple flux rope or plasmoid-like features developing and propagating tailward. MHD results support plasma changes observed by the THEMIS-ARTEMIS satellites. Lunar-scale simulations using the Amitis hybrid code show a short and misaligned plasma wake during the Moon's brief entry into the magnetosheath at 11:20 UT, with plasma expansion into the void being aided by the higher plasma temperatures. Sharply accelerated flow speed and a compressed magnetic field lead to an enhanced electric field in the lunar wake capable of generating sudden changes to the nightside near-surface electric potential.

A. P. Rasca↗
NASA NTRS1993

Cosmological gamma-ray burst model

A cosmological gamma-ray burst model that reproduces the observed gamma-ray spectra was developed. This model, which is an outgrowth of work on synchrotron emission from cosmological sources, creates the observed spectra from a power law spectrum through Compton attenuation in dense molecular clouds. It restricts the burst source to the centers of galaxies, and it is easily tested through comparisons with time dependent burst spectra. Dr. Brainerd continued development of a Monte Carlo code that tracks the random walk of a gamma-ray through a plasma in a strong magnetic field and used this code to show that a two temperature plasma in a strong magnetic field can suppress soft x-ray emission.

Brainerd, J. J.↗
NASA NTRS2007

Development of Numerical Tools for the Investigation of Plasma Detachment from Magnetic Nozzles

A multidimensional numerical simulation framework aimed at investigating the process of plasma detachment from a magnetic nozzle is introduced. An existing numerical code based on a magnetohydrodynamic formulation of the plasma flow equations that accounts for various dispersive and dissipative processes in plasmas was significantly enhanced to allow for the modeling of axisymmetric domains containing three.dimensiunai momentum and magnetic flux vectors. A separate magnetostatic solver was used to simulate the applied magnetic field topologies found in various nozzle experiments. Numerical results from a magnetic diffusion test problem in which all three components of the magnetic field were present exhibit excellent quantitative agreement with the analytical solution, and the lack of numerical instabilities due to fluctuations in the value of del(raised dot)B indicate that the conservative MHD framework with dissipative effects is well-suited for multi-dimensional analysis of magnetic nozzles. Further studies will focus on modeling literature experiments both for the purpose of code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗
NASA NTRS1989

Development and simulation study of a new inverse-pinch high Coulomb transfer switch

The inverse-pinch plasma switch was studied using a computer simulation code. The code was based on a 2-D, 2-temperature magnetohydrodynamic (MHD) model. The application of this code was limited to the disk-type inverse-pinch plasma switch. The results of the computer analysis appear to be in agreement with the experimental results when the same parameters are used. An inverse-pinch plasma switch for closing has been designed and tested for high-power switching requirements. An azimuthally uniform initiation of breakdown is a key factor in achieving an inverse-pinch current path in the switch. Thus, various types of triggers, such as trigger pins, wire-brush, ring trigger, and hypocycloidal-pinch (HCP) devices have been tested for uniform breakdown. Recently, triggering was achieved by injection of a plasma-ring (plasma puff) that is produced separately with hypocycloidal-pinch electrodes placed under the cathode of the main gap. The current paths at switch closing, initiated by the injection of a plasma-ring from the HCP trigger are azimuthally uniform, and the local current density is significantly reduced, so that damage to the electrodes and the insulator surfaces is minimized. The test results indicate that electron bombardment on the electrodes and the insulator surfaces is minimized. The test results indicate that electron bombardment on the electrodes is four orders of magnitude less than that of a spark-gap switch for the same switching power. Indeed, a few thousand shots with peak current exceeding a mega-ampere and with hold-off voltage up to 20 kV have been conducted without showing measurable damage to the electrodes and insulators.

Choi, Sang H.↗
NASA NTRS1996

NASA's Technical Handbook for Avoiding On-Orbit ESD Anomalies Due to Internal Charging Effects

This paper describes NASA-HDBK-4002, "Avoiding Problems Caused by Spacecraft On-Orbit Internal Charging Effects". The handbook includes a description of internal charging and why it is of concern to spacecraft designers. It also suggests how to determine when a project needs to consider internal spacecraft charging, it contains an electron penetration depth chart, rationale for a critical electron flux criterion, a worst-case geosynchronous electron plasma spectrum, general design guidelines, quantitative design guidelines, and a typical materials characteristics list. Appendices include a listing of some environment codes, electron transport codes, a discussion of geostationary electron plasma environments, a brief description of electron beam and other materials tests, and transient susceptibility tests. The handbook will be in the web page, hftp://standards.nasa.gov. A prior document, NASA TP2361 "Design Guidelines for Assessing and controlling Spacecraft Charging Effects", 1984, is in use to describe mitigation techniques for the effects of surface charging of satellites in space plasma environments. HDBK-4002 is meant to complement 2361 and together, the pair of documents describe both cause and mitigation designs for problems caused by energetic space plasmas.

Whittlesey, Albert↗
NASA NTRS1998

A Computational Chemistry Database for Semiconductor Processing

The concept of 'virtual reactor' or 'virtual prototyping' has received much attention recently in the semiconductor industry. Commercial codes to simulate thermal CVD and plasma processes have become available to aid in equipment and process design efforts, The virtual prototyping effort would go nowhere if codes do not come with a reliable database of chemical and physical properties of gases involved in semiconductor processing. Commercial code vendors have no capabilities to generate such a database, rather leave the task to the user of finding whatever is needed. While individual investigations of interesting chemical systems continue at Universities, there has not been any large scale effort to create a database. In this presentation, we outline our efforts in this area. Our effort focuses on the following five areas: 1. Thermal CVD reaction mechanism and rate constants. 2. Thermochemical properties. 3. Transport properties.4. Electron-molecule collision cross sections. and 5. Gas-surface interactions.

Jaffe, R.↗
NASA NTRS1981

Boundary conditions which lead to excitation of instabilities in plasma simulations

Two examples of two-dimensional electrostatic particle-code simulations are shown in which one exhibits characteristics of a stable plasma while the other exhibits unstable, long wavelength plasma oscillations. The only difference between the two simulations is a change in the boundary condition on the electrostatic potential. An energy theorem is derived which shows that the rate of change of field and particle energy within a closed volume is related to a surface integral involving the electrostatic potential and the normal component of the electric current. An analytic theory is developed for a one-dimensional plasma to show how boundary effects can excite spurious plasma instabilities. The theory is tested with a series of one-dimensional plasma simulations. Finally, practical considerations on means of avoiding the non-physical instabilities in simulation plasmas are given.

Swift, D. W.↗
NASA NTRS1988

Simulations of plasma structure evolution in the high-latitude ionosphere

The Naval Research Laboratory has recently developed a two-dimensional inertial, electrostatic code which has been successfully applied to the development and evolution of ionospheric structure driven by plasma instabilities. This code models the ionosphere and magnetosphere as a set of horizontal two-dimensional layers which are coupled by the vertical magnetic field lines at high latitudes. It is shown that the development of instability-generated structure can be strongly dependent on this coupling. For example, the influence of magnetospheric coupling on the E x B gradient drift instability is to retard the instability's growth and to isotropize density irregularities. The influence of ionospheric coupling on the Kelvin-Helmholtz instability is to retard its growth and to suppress vortex formation. This paper presents the results of numerical simulations of these instabilities and discusses their application to high-latitude ionospheric structure.

Huba, J. D.↗
NASA NTRS2002

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗