Advances in fully-kinetic PIC simulations of near-vacuum Hall thruster and other plasma systems

Abstract

In recent years, many groups have numerically modeled the near-anode region of a Hall thruster in attempts to better understand the associated physics of thruster operation. Originally, simulations assumed a continuum approximation for electrons and used magnetohydrodynamic fluid equations to model the significant processes. While these codes were computationally efficient, their applicability to non-equilibrated regions of the thruster, such as wall sheaths, was limited, and their accuracy was predicated upon the notion that the energy distributions of the various species remained Maxwellian at all times. The next generation of simulations used the fully-kinetic particle-in-cell (PIC) model. Although much more computationally expensive than the fluid codes, the full-PIC codes allowed for non-equilibrated thruster regions and did not rely on Maxwellian distributions. However, these simulations suffered for two main reasons. First, due to the high computational cost, fine meshing near boundaries which would have been required to pro]perly resolve wall sheaths was often not attempted. Second, PIC is inherently a statistically noisy method and often the extreme tails of energy distributions would not be adequately sampled due to high energy particle dissipation. The current work initiates a third generation of Hall thruster simulation. A PIC-Vlasov hybrid model was implemented utilizing adaptive meshing techniques to enable automatically scalable resolution of fine structures during the simulation. The code retained the accuracy and versatility of a PIC simulation while intermittently recalculating and smoothing particle distribution functions within individual cells to ensure full velocity space coverage.

(cont.) In addition, this simulation extended the state of the art in Hall thruster anomalous diffusion modeling by adopting a "quench rule" which is able to predict the spatial and temporal structure of the cross-field transport without the aid of prior empirical data. Truly predictive computations are thus enabled. After being thoroughly tested and benchmarked, the simulation was then applied to the near vacuum Hall thruster recently constructed at MIT. Recommendations to improve that thruster's performance were made based on the simulation's results, and those optimizations are being experimentally implemented by other researchers. This work was conducted with the aid of Delta Search Labs' supercomputing facility and technical expertise. The simulation was fully-parallelized using MPI and tested on a 128 processor SGI Origin machine. We gratefully acknowledge that funding for portions of this work has been provided by the United States Air Force and the National Science Foundation.