Search NASASearch

Engineering topics

Koga, J.

Publications and source records attributed to Koga, J..

A simulation study of radial expansion of an electron beam injected into an ionospheric plasma

Injections of nonrelativistic electron beams from a finite equipotential conductor into an ionospheric plasma have been simulated using a two-dimensional electrostatic particle code. The purpose of the study is to survey the simulation parameters for understanding the dependence of beam radius on physical variables. The conductor is charged to a high potential when the background plasma density is less than the beam density. Beam electrons attracted by the charged conductor are decelerated to zero velocity near the stagnation point, which is at a few Debye lengths from the conductor. The simulations suggest that the beam electrons at the stagnation point receive a large transverse kick and the beam expands radially thereafter. The buildup of beam electrons at the stagnation point produces a large electrostatic force responsible for the transverse kick. However, for the weak charging cases where the background plasma density is larger than the beam density, the radial expansion mechanism is different; the beam plasma instability is found to be responsible for the radial expansion. The simulations show that the electron beam radius for high spacecraft charging cases is of the order of the beam gyroradius, defined as the beam velocity divided by the gyrofrequency. In the weak charging cases, the beam radius is only a fraction of the beam gyroradius. The parameter survey indicates that the beam radius increases with beam density and decreases with magnetic field and beam velocity. The beam radius normalized by the beam gyroradius is found to scale according to the ratio of the beam electron Debye length to the ambient electron Debye length. The parameter dependence deduced would be useful for interpreting the beam radius and beam density of electron beam injection experiments conducted from rockets and the space shuttle.

Koga, J.

A parallel particle-in-cell model for the massively parallel processor

The availability of the nearest-neighbor communication-incorporating Massively Parallel Processor has prompted the development of a two-dimensional, particle-in-cell algorithm which loads particles in a cell randomly onto a row of processors, filling only half of them with particles. Due to the simplification of communications among processors achieved in a row by the vacant processors and the random-particle sequence, the algorithm efficiently sorts particles and performs gather/scatter procedures for collecting charge density according to their cells. The algorithm calculates electric fields at the cells by FFT.

Lin, C. S.

Numerical investigation of a plasma beam entering transverse magnetic fields

Two-dimensional electrostatic and electromagnetic particle-in-cell codes are used to investigate plasma beam propagation across transverse magnetic fields. Both large- and small-ion-gyroradius beams are studied in the case of small beam momentum or energy (low drift kinetic beta). Large-ion-gyroradius beams with a large dielectric constant epsilon much greater than (M/m) to the 1/2 (where M/m is the ion-to-electron mass ratio) are found to propagate across the magnetic field via E x B drifts at nearly the initial injection velocity. Beam degradation and undulations are observed, in agreement with previous experimental and analytical results. When epsilon is of order (M/m) to the 1/2 the plasma beam propagates across field lines at only half its initial velocity and loses its coherent structure. When epsilon is much less than (M/m) to the 1/2 the beam particles decouple at the magnetic field boundary, scattering the electrons and deflecting the ions slightly.

Koga, J.

Excitation of low-frequency waves by auroral electron beams

The electron distribution functions measured by the Dynamics Explorer 1 satellite during an auroral pass in 1981 are used in a linear instability analysis of low-frequency electromagnetic and electrostatic waves near and below the hydrogen gyrofrequency. It is suggested that the low-frequency electric and magnetic noise in the auroral zone might be explained by O and H electromagnetic ion cyclotron waves excited by energetic electron beams. An instability analysis suggests that upward and downward streaming electrons throughout the central plasma sheet region provide the free energy for heating oxygen ion through oxygen electrostatic ion cyclotron waves.

Lin, C. S.

Simulation of radial expansion of an electron beam injected into a background plasma

A 2-D electrostatic particle code was used to study the beam radial expansion of a nonrelativistic electron beam injected from an isolated equipotential conductor into a background plasma. The simulations indicate that the beam radius is generally proportional to the beam electron gyroradius when the conductor is charged to a large potential. The simulations also suggest that the charge buildup at the beam stagnation point causes the beam radial expansion. From a survey of the simulation results, it is found that the ratio of the beam radius to the beam electron gyroradius increases with the square root of beam density and decreases inversely with beam injection velocity. This dependence is explained in terms of the ratio of the beam electron Debye length to the ambient electron Debye length. These results are most applicable to the SEPAC electron beam injection experiments from Spacelab 1, where high charging potential was observed.

Koga, J.

Plasma simulation using the massively parallel processor

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

Lin, C. S.