Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma simulation”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
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 NTRS1983

Recent advances in numerical simulation of space-plasma-physics problems

Computer simulations have become an increasingly popular, important and insightful tool for studying space plasmas. This review describes MHD and particle simulations, both of which treat the plasma and the electromagnetic field in which it moves in a self consistent fashion but on drastically different spatial and temporal scales. The complementary roles of simulation, observations and theory are stressed. Several examples of simulations being carried out in the area of magnetospheric plasma physics are described to illustrate the power, potential and limitations of the approach.

Birmingham, T. J.↗
NASA NTRS1993

Magnetosheath-ionspheric plasma interactions in the cusp/cleft. 2: Mesoscale particle simulations

Ionospheric plasma flowing out from the cusp can be an important source of plasma to the magnetosphere. One source of free energy that can drive this outflow is the injection of magnetosheath plasma into the cusp. Two-dimensional (three velocity) mesoscale particle simulations are used to investigate the particle dynamics in the cusp during southward interplanetary magnetic field. This mesoscale model self-consistently incorporates (1) global influences such as the convection of plasma across the cusp, the action of the mirror force, and the injection of the magnetosheath plasma, and (2) wave-particle interactions which produce the actual coupling between the magnetosheath and ionospheric plasmas. It is shown that, because the thermal speed of the electrons is higher than the bulk motion of the magnetosheath plasma, an upward current is formed on the equatorward edge of the injection region with return currents on either side. However, the poleward return currents are the stronger due to the convection and mirroring of many of the magnetosheath electrons. The electron distribution in this latter region evolves from upward directed streams to single-sided loss cones or possibly electron conics. The ion distribution also shows a variety of distinct features that are produced by spatial and/or temporal effects associated with varying convection patterns and wave-particle interactions. On the equatorward edge the distribution has a downflowing magnetosheath component and an upflowing cold ionospheric component due to continuous convection of ionospheric plasma into the region. In the center of the magnetosheath region, heating from the development of an ion-ion streaming instability causes the suppression of the cold ionospheric component and the formation of downward ionospheric streams. Further poleward there is velocity filtering of ions with low pitch angles, so that the magnetosheath ions develop a ring-beam distribution and the ensuing wave instabilities generate downward ionospheric conics. These downward ionospheric components are eventually turned by the mirror force, leading to the production of upward conics at elevated energies throughout the region.

Winglee, R. M.↗
NASA NTRS1992

Relativistic, perpendicular shocks in electron-positron plasmas

One-dimensional particle-in-cell plasma simulations are used to examine the mechanical structure and thermalization properties of collisionless relativistic shock waves in electron-positron plasmas. Shocks propagating perpendicularly to the magnetic field direction are considered. It is shown that these shock waves exist, and that they are completely parameterized by the ratio of the upstream Poynting flux to the upstream kinetic energy flux. The way in which the Rankine-Hugoniot shock jump conditions are modified by the presence of wave fluctuations is shown, and they are used to provide a macroscopic description of these collisionless shock flows. The results of a 2D simulation that demonstrates the generality of these results beyond the assumption of the 1D case are discussed. It is suggested that the thermalization mechanism is the formation of a synchrotron maser by the coherently reflected particles in the shock front. Because the downstream medium is thermalized, it is argued that perpendicular shocks in pure electron-positron plasmas are not candidates as nonthermal particle accelerators.

Gallant, Yves A.↗
NASA NTRS1989

Monte Carlo simulation of reentry plasmas

Attention is given to the treatment of ionization and plasma effects in the direct simulation Monte Carlo method. The requirements for accurate modeling of reentry plasmas are discussed along with the difficulties these requirements present. The current method for modeling such plasmas is reviewed and an alternative method is presented. Both methods are applied to the flow of a 10 km/s shock wave in air at 0.1 torr; a flowfield directly relevant to the projected aeroassisted orbital transfer vehicle. The results are compared and the differences between the methods are discussed.

Carlson, Ann B.↗
NASA NTRS1990

Perpendicular heating of electrons by upper hybrid waves generated by a ring distribution

Satellite observations of electron conical distributions with enhanced fluxes just outside the loss cone suggest that telectrons have been heated perpendicularly to the magnetic field in the mid-altitude polar magnetosphere. To understand electron conical distributions, plasma simulations are conducted to examine an upper hybrid wave instability of a ring electron distribution perpendicular to the magnetic field in a cold electron background. The simulations indicate that both the cold and ring distributions are heated perpendicularly during the saturation stage. From the plasma data, a ring distribution can be identified as a trapped distribution function with an enhancement near 90-deg pitch angle in the phase space density plot. It is suggested that the ring distribution might provide an additional free energy source for generating upper hybrid waves associated with electron conical events.

Lin, C. S.↗
NASA NTRS1992

The dynamics of low-beta plasma clouds as simulated by a three-dimensional, electromagnetic particle code

A 3D electromagnetic and relativistic particle simulation code is used to investigate the dynamics of low-beta plasma clouds moving perpendicular to an ambient magnetic field in vacuum and in a background plasma. The simulations show the formation of the space charge sheaths at the sides of the cloud with the associated polarization electric field which facilitate the cross-field propagation, as well as the sheaths at the front and rear end of the cloud caused by the larger ion Larmor radius, which allows ions to move ahead and lag behind the electrons as they gyrate. It is found that, in a background plasma, electron and ion sheaths expand along the magnetic field at the same rate, whereas in vacuum the electron sheath expands much faster than the ion sheath. Sheath clouds are accelerated up to relativistic energies, which indicates that artificial plasma clouds released in the ionosphere or magnetosphere may generate optical emissions as energetic sheath electrons scatter in the upper atmosphere.

Neubert, T.↗
NASA NTRS1990

Injection of an overdense electron beam in space

A three-dimensional particle simulation model is used to study the injection and propagation of an overdense electron beam in the vicinity of a conductor in space. Beam electrons with a density of more than 100 times the ambient electron density are modeled using large-scale plasma simulations; in these simulations the surface area of the conductor is several thousands times that of the beam cross section at the injection point. The parameters of the simulations are chosen to allow the realistic simulation of active space experiments such as the Space Shuttle's Spacelab 2 electron beam mission. These simulations confirm space observations that an overdense electron beam can at least partially escape the near vicinity of the spacecraft, even in a fully ionized plasma. Once they have escaped from the vicinity of the spacecraft, these beam electrons should be able to propagate away freely until their energy is dissipated. It is suggested that such large-scale simulations could be useful in interpreting data from space experiments.

Okuda, Hideo↗
NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗
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 NTRS1987

Plasma physics abstracts, 1 January - 31 December, 1986

Topics addressed include: ion-cyclotron waves; plasma waves; solar wind lithium releases; bow shock; Pi2 wave bursts; auroral kilometric radiation; ion energization; magnetic field corrections; electric fields; magnetospheric processes; electron acceleration; inner heliosphere; nightside auroral zone; computerized simulation; plasma wave turbulence; and magnetohydrodynamic waves in plasma sheets.

Gurnett, D. A.↗
NASA NTRS1971

Experimental investigations of a uranium plasma pertinent to a self-sustaining plasma source

The research is pertinent to the realization of a self-sustained fissioning plasma for applications such as nuclear propulsion, closed cycle MHD power generation using a plasma core reactor, and heat engines such as the nuclear piston engine, as well as the direct conversion of fission energy into optical radiation (nuclear pumped lasers). Diagnostic measurement methods and experimental devices simulating plasma core reactor conditions are discussed. Studies on the following topics are considered: (1) ballistic piston compressor (U-235); (2) high pressure uranium plasma (natural uranium); (3) sliding spark discharge (natural uranium); (4) fission fragment interaction (He-3 and U-235); and (5) nuclear pumped lasers (He-3 and U-235).

Schneider, R. T.↗
NASA NTRS2003

Gyrophase-bunched Electrons: Cluster Observations and Simulations

The interaction of the solar wind with the Earth's bow shock is responsible for a number of phenomena both upstream and downstream. One of the least understood and studied of these is the generation of gyrophase-bunched electrons. We will describe initial work searching for and analyzing such events in (Plasma Electron And Current Experiment (PEACE) data from CLUSTER. Plasma simulations will complement the data analysis. Recent analyses suggest that gyrophase-bunched electrons are present rather frequently in the upstream region, due either to energization in the reflection off the shock front or, more likely, to phase trapping in locally produced whistler wave fields.

Gurgiolo, C.↗
NASA NTRS2004

Modeling of Sheath Ion-Molecule Reactions in Plasma Enhanced Chemical Vapor Deposition of Carbon Nanotubes

In many plasma simulations, ion-molecule reactions are modeled using ion energy independent reaction rate coefficients that are taken from low temperature selected-ion flow tube experiments. Only exothermic or nearly thermoneutral reactions are considered. This is appropriate for plasma applications such as high-density plasma sources in which sheaths are collisionless and ion temperatures 111 the bulk p!asma do not deviate significantly from the gas temperature. However, for applications at high pressure and large sheath voltages, this assumption does not hold as the sheaths are collisional and ions gain significant energy in the sheaths from Joule heating. Ion temperatures and thus reaction rates vary significantly across the discharge, and endothermic reactions become important in the sheaths. One such application is plasma enhanced chemical vapor deposition of carbon nanotubes in which dc discharges are struck at pressures between 1-20 Torr with applied voltages in the range of 500-700 V. The present work investigates The importance of the inclusion of ion energy dependent ion-molecule reaction rates and the role of collision induced dissociation in generating radicals from the feedstock used in carbon nanotube growth.

Hash, David B.↗
NASA NTRS1979

Numerical simulation of turbulent flows with a three-dimensional vortex-in-cell method

A three-dimensional vortex-in-cell method has been developed for the evaluation of local flow fields due to a family of vortex filaments which employs the principles and architecture of a code developed for magnetic field evaluation in plasma simulations. The computational effort in the new method, as compared to 'vortex pushing' by direct Biot-Savarat interaction, increases directly, rather than quadratically, with the number of vortex elements. The method is well suited for studying large number of vortex filaments or rings and can be used for simulating continuous vorticity.

Buneman, O.↗
NASA NTRS1990

Simulation studies of plasma waves in the electron foreshock - The transition from reactive to kinetic instability

Particle simulation experiments were used to analyze the electron beam-plasma instability. It is shown that there is a transition from the reactive state of the electron beam-plasma instability to the kinetic instability of Langmuir waves. Quantitative tests, which include an evaluation of the dispersion relation for the evolving non-Maxwellian beam distribution, show that a quasi-linear theory describes the onset of this transition and applies again fully to the kinetic stage. This stage is practically identical to the late stage seen in simulations of plasma waves in the electron foreshock described by Dum (1990).

Dum, C. T.↗
NASA NTRS1982

Threshold criterion for a space simulation beam-plasma discharge

An experimental and theoretical study of the threshold characteristics of a space simulation beam-plasma discharge with emphasis on density profiles and a density-dependent ignition criterion. The study included various beam-plasma conditions covering beam currents from 8 to 85 mA, beam energies from 0.8 to 2.0 keV, and magnetic fields at 0.9 and 1.5 G. The study included experimental determinations of radial profiles of electron density for each of the selected conditions extending from a low-density, pre-beam-plasma discharge state to a strong beam-plasma discharge condition. The experimental results are shown to agree with detailed model calculations, which consider the beam-plasma discharge to be produced by large-amplitude electron plasma waves resulting from the beam-plasma interaction.

Szuszczewicz, E. P.↗