摘要：
离子推力器推力密度分布对航天器轨道维持和修正具有重要影响.采用粒子模拟-蒙特卡罗碰撞方法模拟束流等离子体输运过程,分析束流多组分粒子喷出数量和速度等微观参数,并计算得到单孔束流推力,结合放电室出口等离子体密度分布,进一步对推力密度分布特性分析,最后开展实验验证.研究结果显示:束流中单价离子、双荷离子以及交换电荷离子的推力贡献比分别为84.63%, 15.35% 和1.82%,可见推力主要来源于束流中的单价离子和双荷离子,交换电荷离子对推力贡献很小;推力密度分布具有较好的中心轴对称性,从推力器中心沿着径向先快速下降后趋于缓慢;与实验结果对比,经验模型相对误差约为4.1%,数值模型相对误差约为2.8%,相比经验模型,数值模型具有更好的准确性.研究结果可为离子推力器推力密度分布均匀性等优化提供参考.%Thrust density distribution of ion thruster is an important factor that affects the orbit correction and station keeping of the spacecraft. Current empirical models mainly concern themselves with the overall thrust of the ion thruster, yet the thrust density distribution has not been fully understood. Hence it is necessary to investigate the thrust density characteristics of the ion thruster to devise the approach to optimizing the thruster performances. In this study, the thrust density characteristics of the ion thruster is analyzed and discussed by combining the empirical and theoretical methods. An ion thruster utilizes biased grids to extract ions from discharge chamber and accelerate them to high velocities, thereby forming a beam and generating thrust. In this paper, we analyze the working process of the ion thruster. The thrust expression as a function of beam micro-particle parameters is presented. Meanwhile the transport process of the plasma in the beam stream is simulated by the particle in cell-Monte Carlo (PIC-MCC) method for two-grid optics. The motion behavior of ions is modeled by the PIC method,while the collisions of particles are modeled by the MCC method. In the simulation, the particle trajectories are traced and the micro information about ejected charged ions is recorded with respect to singly charged ion, doubly charged ion and charge exchanged (CEX) ion. By analyzing the density and axial velocity of the charged particles in the beam stream,the thrust of the beam from a single grid hole can be calculated, based on which the thrust distribution of the thruster can be inferred by considering the distribution of plasma density at the exit of discharge chamber. Moreover, the above theoretical analysis of the thrust density is tested experimentally. The studies show that the thrust contribution percentages of the singly charged ion, doubly charged ion and CEX ion in the beam current are 84.63%, 15.35%, and 1.82%, respectively. Apparently, the main contributions to the thrust are made by the singly charged ions and doubly charged ions in the beam plasma,while the CEX ions have a trivial effect on the variation of the thrust. The distribution of the thrust density shows good symmetry along the central axis and it levels off after a fast decline in the radial direction. Comparisons of empirical and numerical results with the experimental results show that the empirical results have an error of about 4.1% and the numerical results have an error of about 2.8%. This indicates that the computational accuracy of our numerical model is better than that of the empirical model This work provides a reference for optimizing the thrust density uniformity of an ion thruster.
