TY - JOUR
T1 - Kinetic, 3-D, PIC-DSMC Simulations of Ion Thruster Plumes and the Backflow Region
AU - Jambunathan, Revathi
AU - Levin, Deborah A.
N1 - Funding Information:
Manuscript received May 21, 2019; revised September 30, 2019 and December 29, 2019; accepted March 18, 2020. Date of publication April 29, 2020; date of current version June 10, 2020. This work was supported in part by Air Force Office of Scientific Research (AFOSR) under Grant AF FA9550-16-1-0193, in part by the National Science Foundation through the Blue Waters sustained-petascale computing project under Award OCI-0725070 and Award ACI-1238993, and in part by the state of Illinois. Blue Waters is a joint effort of the University of Illinois, Urbana-Champaign, and its National Center for Supercomputing Applications. The review of this article was arranged by Senior Editor F. Taccogna. (Corresponding author: Revathi Jambunathan.) The authors are with the Department of Aerospace Engineering, University of Illinois at Urbana–Champaign, Urbana, IL 61820 USA (e-mail: rjambunathan@lbl.gov).
Funding Information:
This work was supported in part by Air Force Office of Scientific Research (AFOSR) under Grant AF FA9550-16-1-0193, in part by the National Science Foundation through the Blue Waters sustained-petascale computing project under Award OCI-0725070 and Award ACI-1238993, and in part by the state of Illinois. Blue Waters is a joint effort of the University of Illinois, Urbana-Champaign, and its National Center for Supercomputing Applications. The review of this article was arranged by Senior Editor F. Taccogna. (Corresponding author: Revathi Jambunathan.)
Publisher Copyright:
© 1973-2012 IEEE.
PY - 2020/6
Y1 - 2020/6
N2 - Three-dimensional, fully kinetic, coupled particle-in-cell direct simulation Monte Carlo (PIC-DSMC) simulations are performed to accurately predict the self-consistent electric field and plume and backflow characteristics of an ion thruster, and are compared with results based on the usual Boltzmann approximation. The fully kinetic calculations are performed for the first time using the operational thruster exit number densities without geometric scaling and the actual electron-to-ion mass ratio. It was found that kinetic treatment of both electrons and ions significantly affects the self-consistent radial electric field as well as the backflow ion characteristics. For the xenon ion plume with a colocated electron source, increase in the thruster exit number density strengthens the radial electric field which, in turn, increases the energy of the backflow ions by an order of magnitude. In addition, the spanwise variation of the ion flux on the solar panel showed that charge-exchange ions mostly impinge at the tip of the solar panel, furthest from the thruster geometry, for the fully kinetic cases. In comparison, the Boltzmann relation predicts a lower radial electric field for the operational high number density case, resulting in near-normal ion incidence angles in the backflow region causing significant flux closer to the thruster geometry.
AB - Three-dimensional, fully kinetic, coupled particle-in-cell direct simulation Monte Carlo (PIC-DSMC) simulations are performed to accurately predict the self-consistent electric field and plume and backflow characteristics of an ion thruster, and are compared with results based on the usual Boltzmann approximation. The fully kinetic calculations are performed for the first time using the operational thruster exit number densities without geometric scaling and the actual electron-to-ion mass ratio. It was found that kinetic treatment of both electrons and ions significantly affects the self-consistent radial electric field as well as the backflow ion characteristics. For the xenon ion plume with a colocated electron source, increase in the thruster exit number density strengthens the radial electric field which, in turn, increases the energy of the backflow ions by an order of magnitude. In addition, the spanwise variation of the ion flux on the solar panel showed that charge-exchange ions mostly impinge at the tip of the solar panel, furthest from the thruster geometry, for the fully kinetic cases. In comparison, the Boltzmann relation predicts a lower radial electric field for the operational high number density case, resulting in near-normal ion incidence angles in the backflow region causing significant flux closer to the thruster geometry.
KW - Backflow plasma characteristics
KW - charge-exchange collisions
KW - fully-kinetic PIC-DSMC
KW - ion thruster plume
UR - http://www.scopus.com/inward/record.url?scp=85087061318&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85087061318&partnerID=8YFLogxK
U2 - 10.1109/TPS.2020.2988010
DO - 10.1109/TPS.2020.2988010
M3 - Article
AN - SCOPUS:85087061318
SN - 0093-3813
VL - 48
SP - 2017
EP - 2034
JO - IEEE Transactions on Plasma Science
JF - IEEE Transactions on Plasma Science
IS - 6
M1 - 9082014
ER -