Power transfer efficiency in an air-breathing radio frequency ion thruster

Due to a series of challenges such as low-orbit maintenance of satellites, the air-breathing electric propulsion has got widespread attention. Commonly, the radio frequency ion thruster is favored by low-orbit missions due to its high specific impulse and efficiency. In this paper, the power transfer efficiency of the radio frequency ion thruster with different gas compositions is studied experimentally, which is obtained by measuring the radio frequency power and current of the antenna coil with and without discharge operation. The results show that increasing the turns of antenna coils can effectively improve the radio frequency power transfer efficiency, which is due to the improvement of Q factor. In pure N_2 discharge,with the increase of radio frequency power, the radio frequency power transfer efficiency first rises rapidly and then exhibits a less steep increasing trend. The radio frequency power transfer efficiency increases with the increase of gas pressure at relatively high power, while declines rapidly at relatively low power. In N_(2)/O_(2) discharge, increasing the N_(2) content at high power can improve the radio frequency power transfer efficiency, but the opposite was observed at low power. In order to give a better understanding of these trends, an analytic solution in limit cases is utilized, and a Langmuir probe was employed to measure the electron density. It is found that the evolution of radio frequency power transfer efficiency can be well explained by the variation of plasma resistance, which is related to the electron density and the effective electron collision frequency.

A study of the influence of different grid structures on plasma characteristics in the discharge chamber of an ion thruster

The grid structure has significant effects on the discharge characteristics of an ion thruster.The discharge performances of a 30 cm diameter ion thruster with flat,convex and concave grids are studied.The analysis results show that the discharge chamber with a convex grid has a larger’magnetic-field free area’than the others,and the parallelism of the magnetic-field isopotential lines and anode is generally the same in the three models.Plasma densities of the three structures at the grid outlet are in the range of 3.1×1016-6.9×1017m-3.Along the thruster axis direction,the electron temperature in the chamber with the convex and concave grids is in the range of 3.3-3.5 eV,while that with a flat grid is lower,in the range of 3.1-3.5 eV.In addition,the convex and the concave grids have better uniform distribution of electron temperature.Moreover,the collision frequency ratios show that the axial degree of ionization of the three models is the highest,and the flat grid has the highest discharge efficiency,followed by the convex grid and the concave grid is the least efficient.The test and simulation results of the 30 cm diameter ion thruster with the convex grid show that the measurement and calculation results are 3.67 A and 3.44 A,respectively,and the error above mainly comes from the ignorance of the doubly charged ions and parameter settings in the model.The comparison error between the simulation and measurement of beam current density is mainly caused by the actual thermal deformation of the grids during the discharge process,which leads to the change in electric potential distribution and variation of the focusing characteristics of the grids.Upon consideration of discharge performance and the thermal grid gap variation,it can be concluded that the flat and concave grids are more suitable for small-diameter ion thrusters,while the convex grid is a more reasonable choice for the higher-power and larger-diameter thrusters.

Experimental study of a neutralizer-free gridded ion thruster using radio-frequency self-bias effect

An experimental study on the quasi-neutral beam extracted by a neutralizer-free gridded ion thruster prototype was presented.The prototype was designed using an inductively coupled plasma source terminated by a double-grid accelerator.The beam characteristics were compared when the accelerator was radio-frequency(RF)biased and direct-current(DC)biased.An RF power supply was applied to the screen grid via a blocking capacitor for the RF acceleration,and a DC power supply was directly connected to the screen grid for the DC acceleration.Argon was used as the propellant gas.Furthermore,the characteristics of the plasma beam,such as the floating potential,the spatial distribution of ion flux,and the ion energy distribution function(IEDF)were measured by a four-grid retarding field energy analyzer.The floating potential results showed that the beam space charge is compensated in the case of RF acceleration without a neutralizer,which is similar to the case of classical DC acceleration with a neutralizer.The ion flux of RF acceleration is 1.17 times higher than that of DC acceleration under the same DC component voltage between the double-grid.Moreover,there are significant differences in the beam IEDFs for RF and DC acceleration.The IEDF of RF acceleration has a widened and multipeaked profile,and the main peak moves toward the high-energy region with increasing the DC self-bias voltage.In addition,by comparing the IEDFs with RF acceleration frequencies of3.9 and 7.8 MHz,it is found that the IEDF has a more centered main peak and a narrower energy spread at a higher frequency.

Numerical study of the radio-frequency biased accelerating system in ion thrusters

A 2D-3V implicit immersed-finite-element particle-in-cell(IFE-PIC)model is introduced to investigate the radio-frequency(RF)self-bias accelerating system applied in the RF ion thruster.A set of holes in a two-grid system with slit apertures is simulated in Cartesian coordinates.The characteristics of the plasma plume,such as the ion density,the neutralization rate and the ion and electron current density were investigated for different RF voltage amplitudes(600-1200V)and frequencies(6-30 MHz).Furthermore,the performance of the thruster was also carefully studied.The simulation results show that a well-focused plasma beam can be formed when the voltage amplitude is larger than 900 V and the frequency exceeds the reciprocal of ion transit time(≥12 MHz)in our simulation cases.The performance of the system can be evidently improved by increasing the voltage amplitude and the frequency,and the losses of the particle and thrust are reduced correspondingly.The bulk region of the plasma beam downstream shows good quasi-neutrality,and the ions are dominant in the peripheral region when a well-focused state is achieved.The high ion density beamlet in the periphery of the ion beam is closer to the axis when the voltage amplitude is increasing,while it is expanded radially when increasing the frequency.Backstream electrons have been observed upstream,and this mainly occurs in the phase in which the electron cannot escape.

A study of grid failure mode drivers and methods for accelerated life testing of a 30cm diameter ion thruster

In view of the high cost caused by the 1:1 lifetime verification test of ion thrusters,the lifetime acceleration test should be considered.This work uses the PIC-MCC(Particle-in-Cell MonteCarlo Collision)method to analyze the five failure factors that lead to the failure of the accelerator grid of a 30 cm diameter ion thruster under the working mode of 5 k W.Meanwhile,the acceleration stress levels corresponding to different failure factors are obtained.The results show that background pressure has the highest stress level on the grid's erosion.The accelerator grid aperture's mass sputtering rate under the rated vacuum degree(1×10^(-4)Pa)of 5 k W work mode is 8.78 times that of the baseline vacuum degree(1×10^(-6)Pa),and the mass sputtering rate under worse vacuum degree(5×10^(-3)Pa)is 5.08 times that of 1×10^(-4)Pa.Under the influence of the other four failure factors,namely,the voltage of the accelerator grid,upstream plasma density,the screen grid voltage and mass utilization efficiency,the mass sputtering rates of the accelerator grid hole are 2.32,2.67,1.98 and 2.51 times those of the accelerator grid hole under baseline condition,respectively.The ion sputtering results of two 30 cm diameter ion thrusters(both installed with new grids assembly)after working for 1000 h show that the mass sputtering rate of the accelerator grid hole under vacuum conditions of 5×10^(-3)Pa is 4.54 times that under the condition of 1×10^(-4)Pa,and the comparison error between simulation results and test results of acceleration stress is about 10%.In the subsequent ion thruster lifetime verification,the working vacuum degree can be adjusted according to the acceleration stress level of background pressure,so as to shorten the test time and reduce the test cost.

Particle Simulation of Three-Grid ECR Ion Thruster Optics and Erosion Prediction

Ion optics is a critical component of ion thrusters. A two-dimensional axisymmetric model is developed to study the characteristics of three-grid electron cyclotron resonance ion thruster optics. The code is based on a particle-in-cell combined with the Monte Carlo collision method to simulate ion dynamics and charge-exchange processes in the grid region. The simulation results show that the mode can give a reasonable estimate of the physics characteristics of the ion optics. The design of the ion optics satisfies the requirement of preventing electron backstreaming. Charge-exchange ions can cause damage to the grids, especially to the accelerator grid. ＇Barrel＇ erosion can increase the accelerator grid aperture radius at a rate of 1.91~ 10-11 m/s, while the decelerator grid plays an important role in reducing ＇pits-and-grooves＇ erosion.

Study of the key factors affecting the triple grid lifetime of the LIPS-300 ion thruster

In order to ascertain the key factors affecting the lifetime of the triple grids in the LIPS-300 ion thruster,the thermal deformation,upstream ion density and component lifetime of the grids are simulated with finite element analysis,fluid simulation and charged-particle tracing simulation methods on the basis of a 1500 h short lifetime test.The key factor affecting the lifetime of the triple grids in the LIPS-300 ion thruster is obtained and analyzed through the test results.The results show that ion sputtering erosion of the grids in 5 kW operation mode is greater than in the case of 3 kW.In 5 kW mode,the decelerator grid shows the most serious corrosion,the accelerator grid shows moderate corrosion,and the screen grid shows the least amount of corrosion.With the serious corrosion of the grids in 5 kW operation mode,the intercept current of the acceleration and deceleration grids increases substantially.Meanwhile,the cold gap between the accelerator grid and the screen grid decreases from 1 mm to 0.7 mm,while the cold gap between the accelerator grid and the decelerator grid increases from 1 mm to 1.25 mm after 1500 h of thruster operation.At equilibrium temperature with 5 k W power,the finite element method（FEM）simulation results show that the hot gap between the screen grid and the accelerator grid reduces to 0.2 mm.Accordingly,the hot gap between the accelerator grid and the decelerator grid increases to 1.5 mm.According to the fluid method,the plasma density simulated in most regions of the discharge chamber is 1?×?10;-8?×?10;m;.The upstream plasma density of the screen grid is in the range 6?×?10;-6?×?10;m;and displays a parabolic characteristic.The charged particle tracing simulation method results show that the ion beam current without the thermal deformation of triple grids has optimal perveance status.The ion sputtering rates of the accelerator grid hole and the decelerator hole are 5.5?×?10;kg s;and 4.28?×?10;kg s;,respectively,while after the thermal deformation of the triple grids,the ion beam current has over-perveance status.The ion sputtering rates of the accelerator grid hole and the decelerator hole are 1.41?×?10;kg s;and 4.1?×?10;kg s;,respectively.The anode current is a key factor for the triple grid lifetime in situations where the structural strength of the grids does not change with temperature variation.The average sputtering rates of the accelerator grid and the decelerator grid,which were measured during the 1500 h lifetime test in5 k W operating conditions,are 2.2?×?10;kg s;and 7.3?×?10;kg s;,respectively.These results are in accordance with the simulation,and the error comes mainly from the calculation distribution of the upstream plasma density of the grids.

Numerical research of a 2D axial symmetry hybrid model for the radio-frequency ion thruster

Since the high efficiency discharge is critical to the radio-frequency ion thruster(RIT), a 2D axial symmetry hybrid model has been developed to study the plasma evolution of RIT. The fluid method and the drift energy correction of the electron energy distribution function(EEDF) are applied to the analysis of the RIT discharge. In the meantime, the PIC-MCC method is used to investigate the ion beam current extraction character for the plasma plume region. The beam current simulation results, with the hybrid model, agree well with the experimental results, and the error is lower than 11%, which shows the validity of the model. The further study shows there is an optimal ratio for the radio-frequency(RF) power and the beam current extraction power under the fixed RIT configuration. And the beam extraction efficiency will decrease when the discharge efficiency beyond a certain threshold(about 87 W). As the input parameters of the hybrid model are all the design values, it can be directly used to the optimum design for other kinds of RITs and radio-frequency ion sources.

2D hybrid-PIC simulation of the two and three-grid system of ion thruster

A 2D hybrid-PIC simulation model is proposed to investigate the beam extraction phenomena of the ion thruster. In which the electrons of the plasma sheath upstream the accelerator grid are assumed as particles while the downstream are fluid for improving the calculation efficiency. The ion transparency, plasma sheath formation, ion beam extraction characteristic of a two- and three-grid system have been compared in detail in this paper. From the comparison of the appearing time of the under-perveance phenomena in the two- and three-grid system, it illustrated that the two grid system has the wider operation range of the plasma densities than the three-grid one.

Influence of anode temperature on ignition performance of the IRIT4-2D iodine-fueled radio frequency ion thruster

This paper reports the ignition performance of the iodine-fueled radio frequency(RF)ion thruster(IRIT)at different anode temperaturesT_(a0).The experimental results show that the anode temperature plays important role on the ignition process of the IRIT.There were two characteristic temperatures related to the anode:the minimum ignition temperatureT_(i0)and the stable ignition temperatureT_(is),which were much lower than the pipeline temperature and the storage tank temperature.AtT_(a0)<T_(i0),it failed to discharge.WhenT_(i0)≤T_(a0)<T_(is),it was ignited with dramatical oscillations.AtT_(a0)≥T_(is),the discharge was stable in a large anode temperature range.The analysis showed that the different discharge phenomena at different anode temperatures were attributed to the change of iodine flow rate during the process of the iodine deposition-clogging and sublimation-clearing inside the thruster.The research helps improve the preheating design of the iodine-fueled electric thruster.

Three-dimensional particle simulation of ion thruster plume flows with EX-PWS

Ion thruster plumes from a multi-thruster array of different working configurations are simulated by a hybrid fluid-particle software. The particle in cell method is employed to model the transports of ions. The direct simulation Monte Carlo method is used to model momentum and charge exchange （CEX） collisions. The software is based on unstructured grids which make it easy to handle with complex geometry. The results of chamber simulation are compared with experimental data in ion current density and number density, which show good agreements. The maximum difference of current density along the thruster centerline is less than 9.30%. The interaction effects of plumes when multiple thrusters are operating in vacuum are predicted. Distributions of single charged xenon ions are significantly different in the near-field plume flow, however, merge into one in the far downstream region. Moreover, the interaction effect on the spatial distribution of CEX xenon ions is displayed as well.

Investigation of variable aperture on the performance and lifetime of ion thruster

Beam flatness is an important parameter that determines the performance and the lifetime of a gridded ion thruster.To improve the beam flatness of the 30 cm(LIPS-300)ion thruster,variable aperture ion optics that adapts to the decreasing ion density as the radius increases is proposed.It is the ion optics that the screen grid surface is divided into several zones,where the aperture diameter in each zone is determined by the ion density and the electron temperature upstream of the screen grid.The beam current density in the central area is artificially reduced.A particle in cell-Monte Carlo collision model is applied in this work to investigating the effect of variable aperture on the perveance and the maximum beam current per aperture by simulating the extraction,focusing and acceleration processes of ions.Taking into account the engineering implementability,the screen grid surface is divided into four zones.The hole diameter in each zone is decreased from 1.95 mm to 1.8 mm,1.9 mm,1.8 mm and 1.7 mm,respectively.The simulation results show that the maximum ion density in the center area of grid is decreased by10.6%and 6.99%,while it is increased by 6.49%and 22.3%in the edge region,respectively.The beam flatness of the variable aperture ion optics is improved from 0.69 to 0.88.The erosion rate is decreased by 31.9%,but the total beam current is also decreased by 7.15%.The simulation results can provide a valuable reference of the development of the ion thruster.

Numerical Simulation of Ion Extraction Through Ion Thruster Optics

Based on the particle-in-cell （PIC） method, a two-dimensional numerical scheme was developed to investigate the ion beam extraction phenomena through the ion thruster optics. According to the calculated results, the plasma sheath upstream of the screen grid, the electric field in the calculation domain, and the ion and electron spatial distributions are obtained for different accelerator grid voltages. The results indicate that the accelerator grid voltage affects the plasma sheath upstream of the screen grid significantly. It is found that a moderate accelerator grid voltage results in an ion optical performance better than either a higher or lower voltage, from a point of ion extraction from the discharge chamber and erosion mitigation of the accelerator grid due to the direct ion impingement.

Three Dimensional Simulation of Ion Thruster Plume-Spacecraft Interaction Based on a Graphic Processor Unit

Based on the three-dimensional particle-in-cell （PIC） method and Compute Unified Device Architecture （CUDA）, a parallel particle simulation code combined with a graphic processor unit （GPU） has been developed for the simulation of charge-exchange （CEX） xenon ions in the plume of an ion thruster. Using the proposed technique, the potential and CEX plasma distribution are calculated for the ion thruster plume surrounding the DS1 spacecraft at different thrust levels. The simulation results are in good agreement with measured CEX ion parameters reported in literature, and the CPU＇s results are equal to a CPU＇s. Compared with a single CPU Intel Core 2 E6300, 16-processor GPU NVIDIA GeForce 9400 GT indicates a speedup factor of 3.6 when the total macro particle number is 1.1 × 10^6. The simulation results also reveal how the back flow CEX plasma affects the spacecraft floating potential, which indicates that the plume of the ion thruster is indeed able to alleviate the extreme negative floating potentials of spacecraft in geosynchronous orbit.

An equivalent model of discharge instability in the discharge chamber of Kaufman ion thruster

The industrial application of the Kaufman ion thruster in its arc stage is limited owing to the instability of the discharge pulse.Presently,a complete prediction model that can predict the discharge pulse in the high-current stage does not exist.In this study,a complete prediction model for the pulse in the ion thruster is established using the zero-dimensional plasma discharge model and equivalent circuit model.The zero-dimensional plasma discharge model is used to obtain the corresponding plasma parameters by calculating the beam current,discharge current,voltage,and gas flow under actual working conditions.The input parameters of the equivalent circuit model are calculated using empirical formulae to acquire the estimated discharge waveforms.The pulse waveforms obtained using the model are found to be consistent with the experimental results.The model is used to evaluate the process of rapid changes in plasma density.Additionally,this model is employed to predict changes in the pulse waveforms when the volume of the discharge chamber and grid plate transmittance are changed.

Investigation into the thermal effect of the LIPS-200 ion thruster plume

The distribution of the thermal effects of the ion thruster plume are essential for estimating the influence of the thruster plume, improving the layout of the spacecraft, and for the thermal shielding of critical sensitive components. In order to obtain the heat flow distribution in the plume of the LIPS-200 xenon ion thruster, an experimental study of the thermal effects of the plume has been conducted in this work,with a total heat flow sensor and a radiant heat flow sensor over an axial distance of 0.5–0.9 m and a thruster angle of 0°–60°. Combined with a Faraday probe and a retarding potential analyzer, the thermal accommodation coefficient of the sensor surface in the plume is available. The results of the experiment show that the xenon ion thruster plume heat flow is mainly concentrated within a range of15°. The total and radial heat flow of the plume downstream of the thruster gradually decreases along the axial and radial directions, with the corresponding values of 11.78 k W m^(-2) and 0.3 k W m^(-2) for the axial 0.5 m position, respectively. At the same position, the radiation heat flow accounts for a very small part of the total heat flow, approximately 3%–5%. The thermal accommodation factor is0.72–0.99 over the measured region. Furthermore, the PIC and DSMC methods based on the Maxwell thermal accommodation coefficient model(EX-PWS) show a maximum error of 28.6% between simulation and experiment for LIPS-200 ion thruster plume heat flow, which, on the one hand, provides an experimental basis for studying the interaction between the ion thruster and the spacecraft, and on the other hand provides optimization of the ion thruster plume simulation model.

The plume diagnostics of 30 cm ion thruster with an advanced plasma diagnostics system

In order to achieve a better understanding of plume characteristics of LIPS-300 ion thruster,the beam current density,ion energy and electron number density of LIPS-300 ion thruster plume are studied with an Advanced Plasma Diagnostics System(APDS)which allows for simultaneous in situ measurements of various properties characterizing ion thruster,such as plasma density,plasma potential,plasma temperature and ion beam current densities,ion energy distribution and so on.The results show that the beam current density distribution has a double‘wing'shape.The high energy ions were found in small scan angle,while low energy ions were found in greater scan angle.Electron number density has a similar shape with the beam current density distribution.

Estimate of Lifetime of Ion Thruster Optics Based on Particle Simulation

A three-dimensional particle simulation of ion thruster optics with charge-exchange collision was developed in this study. The simulation code was based on tracking ions using the particle-in-cell method, and the Monte Carlo technique was used to model the charge-exchange collision. Simulations were performed for a 20 cm ion thruster optics. The results were compared with the corresponding experimental data from a test of the ion thruster optics for a duration of 800 hours. The Depth-From-Focus （DFF） method was used to measure the erosion depth of the downstream surface of the accelerator grid. The predicted erosion depth of the accelerator grid was consistent reasonably with the corresponding experimental data. The simulation results showed that the accelerator grid would be burned through after 1333 hours.

Wide-range characteristics of beam perveance and saddle point potential of LIPS-200 ion thruster

Both the long-life and multi-mode versions of LIPS-200 ion thruster are under investigation in LIP(Lanzhou Institute of Physics).To confirm the feasible ranges of the beam current and accel(abbreviation for accelaration)grid potential to apply to the thruster,the wide-range beam perveance(the state of beam focus)and saddle point potential(the lowest potential along beamlet centerline)characteristics of LIPS-200 are studied with a test-verified PIC-MCC(Particle in Cell-Monte Carlo Collisions)model.These characteristics are investigated with both the initial and the eroded states of the accel grid aperture diameter.The results show that the feasible ranges of these parameters with respect to perveance/crossover(overfocused)limit extend as the operating time accumulates,while the feasible range of accel grid potential narrows due to a reduced EBSF(electron backstreaming failure)margin.The feasible ranges determined by the initial condition are:(i)the beam current up to 0.981 A,and(ii)the accel grid potential up to−85 V.A 23%enlargement of the aperture diameter would bring up to 48 V of EBSF margin loss.

Numerical simulation of the start-up process of a miniature ion thruster

A particle-in-cell Monte Carlo collision model of a discharge chamber is established to investigate the start-up process of a miniature ion thruster.We present the discharge characteristics at different stages(the initial stage,development stage,and stable stage)according to the trend of the discharge current with time.The discharge current is the sum of the sidewall current and the backplate current.During the start-up process,the sidewall current lags behind the backplate current.The variation and distribution characteristics of the discharge current over time are determined by the electron density distribution and electric potential distribution.