Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle-in-cell”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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 NTRS1992

Effects of magnetospheric electrons on polar plasma outflow - A semikinetic model

The effect or hot magnetospheric electrons on the polar-plasma outflow was investigated, using a semikinetic model developed by Wilson et al. (1990) and Ho et al. (1991) to simulate the effect. The model is based on a hybrid particle-in-cell approach, in which the H(+) and O(+) ions are treated as adiabatic parallel-drifting gyrocenters injected as the upgoing portions of drifting bi-Maxwellian distributions at 1.6 R(E), while the electrons are treated as a massless neutralizing fluid. The results show that, in order to simulate the polar outflow under the influence of hot magnetospheric electrons, it is necessary to consider the effect of the electron temperature gradient.

Ho, C. W.↗
NASA NTRS1992

Further study of the effect of the downstream plasma condition on accelerator grid erosion in an ion thruster

Further numerical results are presented of earlier particle-in-cell/Monte Carlo calculations of accelerator grid erosion in an ion thruster. A comparison between numerical and experimental results suggests that the accelerator grid impingement is primarily due to ions created far downstream from the accelerator grid. In particular, for the same experimental conditions as those of Monheiser and Wilbur at Colorado State University, it is found that a downstream plasma density of 2 x 10 exp 14/cu m is required to give the same ratio of accelerator grid impingement current to beam current (5 percent). For this condition, a potential hill is found in the downstream region of 2.5 V.

Peng, Xiaohang↗
NASA NTRS1992

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗
NASA NTRS1991

A new approach to plasmasphere refilling: Anomalous plasma effects

During the last 10 months of the grant, both laminar and anomalous plasma processes occurring during the refilling of the outer plasmasphere after magnetic storms are investigated. Theoretical investigations were based on two types of models: (1) two-stream hydrodynamic model in which plasma flows from the conjugate ionospheres are treated as separate fluids and the ion temperature anisotropies are treated self-consistently; and (2) large-scale particle-in-cell code.

Singh, N.↗
NASA NTRS1993

Density-conserving shape factors for particle simulations in cylindrical and spherical coordinates

It is established that a particle-to-grid weighting which is based on both conventional particle-in-cell (PIC) and cloud-in-cell (CIC) shape factors will yield nonuniform grid densities, even in the case of uniform particle distributions, in both cylindrical and spherical coordinates. In view of this, alternative density-conserving weighting schemes are discussed which encompass modified PIC and CIC weighting. These alternative shape factors are mathematically no more complicated than the conventional nondensity-conserving expressions.

Ruyten, Wilhelmus M.↗
NASA NTRS1993

Lower hybrid resonance acceleration of electrons and ions in solar flares and the associated microwave emission

The particle acceleration processes here studied are driven by the relaxation of unstable ion ring distributions; these produce strong wave activity at the lower hybrid resonance frequency which collapses, and forms energetic electron and ion tails. The results obtained are applied to the problem posed by the production of energetic particles by solar flares. The numerical simulation results thus obtained by a 2 1/2-dimensional particle-in-cell code show a simultaneous acceleration of electrons to 10-500 keV energies, and of ions to as much as the 1 MeV range; the energy of the latter is still insufficient to account for gamma-ray emission in the 4-6 MeV range, but furnish a seed population for further acceleration.

Mcclements, K. G.↗
NASA NTRS1992

Arcing in space structures in low Earth orbit

This report describes results of an experimental and theoretical program to investigate arcing of structures containing dielectric and conducting materials when they are biased negatively with respect to a plasma. An argon ion source generated Ar(+) ions of directed energy 20 to 40 eV and density approximately 10(exp 7) cm(exp -3) that impinged upon samples containing a dielectric material on top of a negatively biased Kovar plate. Arcing events were studied for bias voltages between -300 and -1000V with respect to the ion beam. The samples were Dow Corning 93-500 adhesive on Kovar, fused silica cover slips bonded on Kovar, and silicon solar cells mounted on Kovar. Measurements of discharge current, Kovar plate voltage, and radiation from the arc versus time were carried out. Microsecond duration exposure images and optical spectra in the 0.24 to 0.40 micron band were also acquired during arcing events. Arcing events were found to be associated with exposed adhesive and means were found to eliminate arcing altogether. The charging of a silica cover plate and the fields around the plate were calculated using a particle-in-cell code. Models were developed to explain the ignition of the arc and the physical processes occurring during the discharge.

Upschulte, B. L.↗
NASA NTRS1994

Extending HPF for advanced data parallel applications

The stated goal of High Performance Fortran (HPF) was to 'address the problems of writing data parallel programs where the distribution of data affects performance'. After examining the current version of the language we are led to the conclusion that HPF has not fully achieved this goal. While the basic distribution functions offered by the language - regular block, cyclic, and block cyclic distributions - can support regular numerical algorithms, advanced applications such as particle-in-cell codes or unstructured mesh solvers cannot be expressed adequately. We believe that this is a major weakness of HPF, significantly reducing its chances of becoming accepted in the numeric community. The paper discusses the data distribution and alignment issues in detail, points out some flaws in the basic language, and outlines possible future paths of development. Furthermore, we briefly deal with the issue of task parallelism and its integration with the data parallel paradigm of HPF.

Chapman, Barbara↗
NASA NTRS1994

Computer simulation of Alfven waves and double layers along auroral magnetic field lines

A plasma simulation has been developed to model interactions between inertial Alfven waves and double layers and to investigate their relative contributions at auroral particle acceleration. We use a novel one-dimensional particle-in-cell code, with periodic boundary conditions, to model the nonlinear excitation of current-driven weak double layers via free energy supplied by an inertial Alfven wave. Analysis of the simulation output shows that double layers are not the agent primarily responsible for electron acceleration. Rather, the inertial Alfven wave accelerates groups of electrons into a steepening beam as it encounters them. As the beam electrons reenter the main distribution, decelerated by anomalous resistive effects, they are replaced by electrons farther downstream. Hence, the particles do not free-stream over the length of the channel. Furthermore, this wave action persists even when the system is linearly stable to ion-acoustic modes, precluding the possibility that this behavior is brought about by the formation of ion-acoustic double layers.

Silberstein, M.↗
NASA NTRS1994

A splitting algorithm for Vlasov simulation with filamentation filtration

A Fourier-Fourier transformed version of the splitting algorithm for simulating solutions of the Vlasov-Poisson system of equations is introduced. It is shown that with the inclusion of filamentation filtration in this transformed algorithm it is both faster and more stable than the standard splitting algorithm. It is further shown that in a scalar computer environment this new algorithm is approximately equal in speed and far less noisy than its particle-in-cell counterpart. It is conjectured that in a multiprocessor environment the filtered splitting algorithm would be faster while producing more precise results.

Klimas, A. J.↗
NASA NTRS1994

Three-dimensional numerical simulation of current collection by a probe in a magnetized plasma

A three-dimensional numerical model for current collection in a magnetized plasma is reported. The model is based on an electrostatic particle-in-cell code. The model yields self-consistent sheath structure including distributions of plasma and the electric potential around the body and the collection of electrons. The analytical theory of current collection by a body in a magnetized plasma yields an upper bound for the collected current determined by the conservation of energy and canonical angular momentum. The theory shows that the collected charged particles come from a cylindrical volume aligned with the magnetic shadow of the body; the maximum radius r(sub o) of this volume is determined by the body size, body potential, and the ambient magnetic field. This theory does not deal with the sheath structure around the body. The condition for the actual current to approach the upper-bound value has been a matter of debate. Our simulations reveal when and why the collected current becomes equal to its upper-bound value. Sheath size in the radial direction perpendicular to the axial ambient magnetic field is determined by the current-limiting radius r(sub o). Our simulation yields time-average current in good agreement with its upper bound. This feature of the current collection is explained as follows: Once electrons enter the sheath, some of them are freely accelerated perpendicular to the magnetic field because they are demagnetized by the large gradients in the perpendicular electric fields. Simulations show a large perpendicular acceleration, producing perpendicular energy as large as that determined by the potential on the body, especially in the region where perpendicular electric fields (E perpendicular) are the strongest. An analysis shows that the demagnetization of electrons occurs above a threshold potential on the body. This threshold condition follows from the breakdown of the adiabaticity of the electron dynamics inside the sheath.

Singh, Nagendra↗
NASA NTRS1994

Relativistic particle acceleration in plerions

We discuss recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is either pure electrons and positrons or primarily electrons and positrons with an admixture of heavy ions. Particle-in-cell simulation techniques as well as analytic theory have been used to show that such shocks in pure pair plasmas are fully thermalized -- the downstream particle spectra are relativistic Maxwellians at the temperature expected from the jump conditions. On the other hand, shocks containing heavy ions which are a minority constituent by number but which carry most of the energy density in the upstream medium do put approximately 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) varies as E(exp -2), where N(E)dE is the number of particles with energy between E and E+dE. The mechanism of thermalization and particle acceleration is found to be synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front. The synchrotron maser modes radiated by the heavy ions are absorbed by the pairs at their (relativistic) cyclotron frequencies, allowing the maximum energy achievable by the pairs to be gamma(sub +/-)m(sub +/-)c squared = m(sub i)c squared gamma(sub 1)/Z(sub i), where gamma(sub 1) is the Lorentz factor of the upstream flow and Z(sub i) is the atomic number of the ions. The shock's spatial structure is shown to contain a series of 'overshoots' in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. This shock model is applied to an interpretation of the structure of the inner regions of the Crab Nebula, in particular to the 'wisps,' surface brightness enhancements near the pulsar. We argue that these surface brightness enhancements are the regions of magnetic overshoot, which appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure gamma(sub 1) and a number of other properties of the pulsar's wind. We also discuss applications of the shock theory to the termination shocks of the winds from rotation-powered pulsars embedded in compact binaries. We show that this model adequately accounts for (and indeed predicted) the recently discovered X-ray flux from PSR 1957+20, and we discuss several other applications to other examples of these systems.

Arons, Jonathan↗
NASA NTRS1994

Early time evolution of negative ion clouds and electron density depletions produced during electron attachment chemical release experiments

Two-dimensional electrostatic particle-in-cell simulations are used to study the early time evolution of electron depletions and negative ion clouds produced during electron attachment chemical releases in the ionosphere. The simulation model considers the evolution in the plane perpendicular to the magnetic field and a three-species plasma that contains electrons, positive ions, and also heavy negative ions that result as a by-product of the electron attachment reaction. The early time evolution (less than the negative ion cyclotron period) of the system shows that a negative charge surplus initially develops outside of the depletion boundary as the heavy negative ions move across the boundary. The electrons are initially restricted from moving into the depletion due to the magnetic field. An inhomogenous electric field develops across the boundary layer due to this charge separation. A highly sheared electron flow velocity develops in the depletion boundary due to E x B and Delta-N x B drifts that result from electron density gradients and this inhomogenous electric field. Structure eventually develops in the depletion boundary layer due to low-frequency electrostatic waves that have growth times shorter than the negative ion cyclotron period. It is proposed that these waves are most likely produced by the electron-ion hybrid instability that results from sufficiently large shears in the electron flow velocity.

Scales, W. A.↗
NASA NTRS1991

Fountain-Like Flow of Heavy Oxygen Ions From the Earth's Ionosphere in Response to Transverse Heating

Normally the gravitationally bound heavy O(+) ions in the Earth's ionosphere are in a diffusive equilibrium. However, when energized to superthermal energies of a few eV transverse to the geomagnetic field, the combined effects of the downward gravitational and the upward electric and mirror forces produce interesting flow patterns in the vertical direction like in a pulsating fountain. This flow pattern is studied by means of a particle-in-cell code.

Singh, Nagendra↗
NASA NTRS1994

Evolution of a localized Langmuir packet in the solar wind and on auroral field lines

Langmuir emissions in space are reported to be clumpy and intermittent. The high-frequency wave power appears concentrated in spatial packets, whether amidst the solar wind or on auroral field lines. Due to the plasma motion relative to the spacecraft, determining the source for the wave free energy in the three-dimensional electron distribution function has always been difficult, since the unstable features pass by the detector in presumably too short time to be measured. The range of unstable phase velocities and growth rates have generally been estimated rather than determined by unequivocal measurements. The analysis of wave-particle interactions in a space environment has taken recently a new turn with the development of wave correlators on board rockets and satellites. Such instruments seek to identify correlations between the phase of the wave-field and the fluxes of energetic particles. The data interpretation is complex, however, it must be backed by a detailed theoretical understanding of the wave-particle interaction, including the phase relation for inhomogeneous packets. To this end Langmuir packets interacting with fast electrons can be studied in the appropriate regime by means of particle-in-cell simulations, provided that one succeeds in reducing the level of the fluctuations, enhancing the signal-to-noise ratio, and incorporating the appropriate boundary conditions. The first results of such simulations are presented here as a test and expansion of previous analysis.

Roth, I.↗
NASA NTRS1989

Ion propulsion

An ion engine is a plasma thruster which produces thrust by extracting ions from the plasma and accelerating them to high velocity with an electrostatic field. The ions are then neutralized and leave the engine as high velocity neutral particles. The advantages of ion engines are high specific impulse and efficiency and their ability to operate over a wide range of input powers. In comparison with other electric thrusters, the ion engine has higher efficiency and specific impulse than thermal electric devices such as the arcjet, microwave, radio frequency, and laser heated thrusters and can operate at much lower current levels than the MPD thruster. However, the thrust level for an ion engine may be lower than a thermal electric thruster of the same operating power, consistent with its higher specific impulse, and therefore ion engines are best suited for missions which can tolerate longer duration propulsive phases. The critical issue for the ion engine is lifetime, since the prospective missions may require operation for several thousands of hours. The critical components of the ion engine, with respect to engine lifetime, are the screen and accelerating grid structures. Typically, these are large metal screens that must support a large voltage difference and maintain a small gap between them. Metallic whisker growth, distortion, and vibration can lead to arcing, and over a long period of time ion sputtering will erode the grid structures and change their geometry. In order to study the effects of long time operation of the grid structure, we are developing computer codes based on the particle-in-cell (PIC) technique and laser induced fluorescence (LIF) diagnostic techniques to study the physical processes which control the performance and lifetime of the grid structures.

Meserole, J. S.↗
NASA NTRS1993

(abstract) 3D Electromagnetic Plasma Particle Simulations

A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.

electromagnetic plasma particles simulations 3D al↗