doi:10.1016/j.cpc.2004.06.045
Copyright © 2004 Elsevier B.V. All rights reserved.
Investigation of ion beam neutralization processes with 2D and 3D PIC simulations
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Adrian Wheelock
,
, a, David L. Cooke
, a and Nikolaos A. Gatsonis
, b
a Air Force Research Laboratory, Space Vehicles Directorate, Hanscom AFB, MA 01731, USA
b Mechanical Engineering Department, Worcester Polytechnic Institute, Worcester, MA 01609, USA
Available online 24 July 2004.
Abstract
While it is common knowledge that ion beams are easily neutralized for both current and charge density using a variety of means, the precise process of neutralization remains unknown. With the increasing importance of electric propulsion, and in particular micropropulsion systems, this question is of significant importance. Additionally, it has bearing on thruster design, space instrument calibration, electrodynamic tethers, and ionospheric research. A review of the present state of knowledge on this topic is presented as well as results from ion beam simulations using 2D and 3D Particle-in-Cell (PIC) codes. We investigate both the early “filling” problem of the beam starting to move away from the spacecraft and the steady state problem where the beam encounters a wall at an infinite distance from the spacecraft.
Author Keywords: Particle-in-Cell (PIC); Electric propulsion; Ion beam neutralization
52.65.Rr; 52.75.Di; 52.40.-w; 52.25.-b
Fig. 1. 2D simulation domain.
Fig. 2. Vx–x phase space plots for 100 eV “filling” beam. Left: Current neutralized, case 1.1. Right: Density neutralized, case 1.2. Top: Ion phase space. Bottom: Electron phase space. t=1.8e–6 s.
Fig. 3. Vx–x phase space plots for 1000 eV “filling” beam. Left: Current neutralized, case 1.3. Right: Density neutralized, case 1.4. Top: Ion phase space. Bottom: Electron phase space. t=1.2e–6 s.
Fig. 4. (A) 3D simulation domain; (B) 3D potential τ=2.2E–7 s.
Fig. 5. (A) 3D Vx–x phase space for ions. τ=2.2E–7 s; (B) 3D Vx–x phase space for electrons. τ=2.2E–7 s.
Table 1. 2D simulation parameters and current-density comparison results

Table 2. 2D velocity ratio simulation parameters and results

Table 3. Steady state simulation parameters and results

“Warm” ions indicate the ions were mobile and allowed to develop from a cold beam. “Frozen” ions indicate particles were place in the problem by the code in a perfectly cold beam.
Table 4. Current coupling simulation parameters and results

“Warm” ions indicate the ions were mobile and started from the developed case 3.2. “Cold” ions indicate the ions were immobile, but started from the developed case 3.2. “Frozen” ions indicates particles were placed in the problem by the code in a cold beam. The right electrode size is in multiples of the left electrode.
Corresponding author, Ph.D. Candidate WPI.