Performance investigation of a low-power Hall thruster fed on iodine propellant

The common propellants used for electric thrusters, such as xenon and krypton, are rare, expensive,and difficult to acquire. Solid iodine attracts much attention with the advantages of low cost,extensive availability, low vapor pressure, and ionization potential. The performance of a lowpower iodine-fed Hall thruster matched with a xenon-fed cathode is investigated across a broad range of operation conditions. Regulation of the iodine vapor's mass flow rates is stably achieved by using a temperature control method of the iodine reservoir. The thrust measurements are finished utilizing a thrust target during the tests. Results show that thrust and anode-specific impulse increase approximately linearly with the increasing iodine mass flow rate.At the nominal power of 200 W class, iodine mass flow rates are 0.62 and 0.93 mg/s, thrusts are7.19 and 7.58 m N, anode specific impulses are 1184 and 826 s, anode efficiencies are 20.8%and 14.5%, and thrust to power ratios are 35.9 and 37.9 m N/k W under the conditions of 250 V,0.8 A and 200 V, 1.0 A, respectively. The operating characteristics of iodine-fed Hall thruster are analyzed in different states. Further work on the measurements of plasma characteristics and experimental optimization will be carried out.

Growth mechanism and characteristics of electron drift instability in Hall thruster with different propellant types

The existence of a significant electron drift instability(EDI) in the Hall thruster is considered as one of the possible causes of the abnormal increase in axial electron mobility near the outlet of the channel. In recent years, extensive simulation research on the characteristics of EDI has been conducted, but the excitation mechanism and growth mechanism of EDI in linear stage and nonlinear stage remain unclear. In this work, a one-dimensional PIC model in the azimuthal direction of the thruster near-exit region is established to gain further insights into the mechanism of the EDI in detail, and the effects of different types of propellants on EDI characteristics are discussed. The changes in axial electron transport caused by EDI under different types of propellants and electromagnetic field strengths are also examined. The results indicate that EDI undergoes a short linear growth phase before transitioning to the nonlinear phase and finally reaching saturation through the ion Landau damping. The EDI drives a significant ion heating in the azimuthal direction through electron–ion friction before entering the quasi-steady state, which increases the axial mobility of the electrons. Using lighter atomic weight propellant can effectively suppress the oscillation amplitude of EDI, but it will increase the linear growth rate, frequency, and phase velocity of EDI. Compared with the classical mobility, the axial electron mobility under the EDI increases by three orders of magnitude, which is consistent with experimental phenomena. The change of propellant type is insufficient to significantly change the axial electron mobility. It is also found that the collisions between electrons and neutral gasescan significantly affect the axial electron mobility under the influence of EDI, and lead the strength of the electric field to increase and the strength of the magnetic field to decrease, thereby both effectively suppressing the axial transport of electrons.

Influence of acceleration stage electrode voltage on the performance of double-stage Hall effect thruster with adjustable zero magnetic point

The configuration of electrode voltage and zero magnetic point position has a significant impact on the performance of the double-stage Hall effect thruster. A 2D-3V model is established based on the two-magnetic peak type double-stage Hall thruster configuration, and a particle-in-cell simulation is carried out to investigate the influences of both acceleration electrode voltage value and zero magnetic point position on the thruster discharge characteristics and performances.The results indicate that increasing the acceleration voltage leads to a larger potential drop in the acceleration stage, allowing ions to gain higher energy, while electrons are easily absorbed by the intermediate electrode, resulting in a decrease in the anode current and ionization rate. When the acceleration voltage reaches 500 V, the thrust and efficiency are maximized, resulting in a 15%increase in efficiency. After the acceleration voltage exceeds 500 V, a potential barrier forms within the channel, leading to a decrease in thruster efficiency. Further study shows that as the second zero magnetic point moves towards the outlet of the channel, more electrons easily traverse the zero magnetic field region, participating in the ionization. The increase in the ionization rate leads to a gradual enhancement in both thrust and efficiency.

Thrust estimate method of an on-orbit Hall thruster using Hall drift current

In order to realize the thrust estimation of the Hall thruster during its flight mission,this study establishes an estimation method based on measurement of the Hall drift current.In this method,the Hall drift current is calculated from an inverse magnetostatic problem,which is formulated according to its induced magnetic flux density detected by sensors,and then the thrust is estimated by multiplying the Hall drift current with the characteristic magnetic flux density of the thruster itself.In addition,a three-wire torsion pendulum micro-thrust measurement system is utilized to verify the estimate values obtained from the proposed method.The errors were found to be less than 8%when the discharge voltage ranged from 250 V to 350 V and the anode flow rate ranged from 30 sccm to 50 sccm,indicating the possibility that the proposed thrust estimate method could be practically applied.Moreover,the measurement accuracy of the magnetic flux density is suggested to be lower than 0.015 mT and improvement on the inverse problem solution is required in the future.

Influence of the position relationship between the cathode and magnetic separatrix on the discharge process of a Hall thruster

Previous studies have shown that there is an obvious coupling relationship between the installation location of the external cathode and the magnetic separatrix in the plume region of a Hall thruster.In this paper,the particle-in-cell simulation method is used to compare the thruster discharge process under the conditions of different position relationships between the cathode and the magnetic separatrix.By comparing the distribution of electron conduction,potential,plasma density and other microscopic parameters,we try to explain the formation mechanism of the discharge difference.The simulation results show that the cathode inside and outside the magnetic separatrix has a significant effect on the distribution of potential and plasma density.When the cathode is located on the outer side of the magnetic separatrix,the potential above the plume region is relatively low,and there is a strong potential gradient above the plume region.This potential gradient is more conducive to the radial diffusion of ions above the plume,which is the main reason for the strong divergence of the plume.The distribution of ion density is also consistent with the distribution of potential.When the cathode is located on the outer side of the magnetic separatrix,the radial diffusion of ions in the plume region is enhanced.Meanwhile,by comparing the results of electron conduction,it is found that the traiectories of electrons emitted from the cathode are significantly different between the inner and outer sides of the magnetic separatrix.This is mainly because the electrons are affected by the magnetic mirror effect of the magnetic tip,which makes it difficult for the electrons to move across the magnetic separatrix.This is the main reason for the difference in potential distribution.In this paper,the simulation results of macroscopic parameters under several conditions are also compared,and they are consistent with the experimental results.The cathode is located on the inner side of the magnetic separatrix,which can effectively reduce the plume divergence angle and improve the thrust.In this paper,the cathode moves from R=50 mm to R=35 mm along the radial direction,the thrust increases by 3.6 mN and the plume divergence angle decreases by 23.77%.Combined with the comparison of the ionization region and the peak ion density,it is found that the main reason for the change in thrust is the change in the radial diffusion of ions in the plume region.

Experimental investigation of the polarityswitching process with different bipolar ionic liquid thruster operating frequencies

The bipolar ionic liquid thruster employs ionic liquid as a propellant to discharge positively and negatively charged high-energy particles under an alternating current(AC)power source,effectively suppressing electrochemical reaction and ensuring charge neutrality.Determining an optimal AC supply power source frequency is critical for sustained stable thruster operation.This study focuses on the emission characteristics of the ionic liquid thruster under varied AC conditions.The AC power supply was set within the frequency range of 0.5-64 Hz,with eight specific frequency conditions selected for experimentation.The experimental results indicate that the thruster operates steadily within a voltage range of±1470 to±1920 V,with corresponding positive polarity current ranging from 0.41 to 4.91μA and negative polarity current ranging from−0.49 to−4.10μA.During voltage polarity switching,an emission delay occurs,manifested as a prominent peak signal caused by circuit capacitance characteristics and a minor peak signal resulting from liquid droplets.Extended emission test was conducted at 16 Hz,demonstrating approximately 1 h and 50 min of consistent emission before intermittent discharge.These findings underscore the favorable impact of AC conditions within the 8-16 Hz range on the self-neutralization capability of the ionic liquid thruster.

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.

Performance of pulsed plasma thruster at low discharge energy

As the size of satellites scales down, low-power and compact propulsion systems such as the pulsed plasma thruster(PPT) are needed for stabilizing these miniature satellites in orbit. Most PPT systems are operated at 2 J or more of discharge energy. In this work, the performance of a PPT with a side-fed, tongue-flared electrode configuration operated within a lower discharge energy range of 0.5-2.5 J has been investigated. Ablation and charring of the polytetrafluoroethylene propellant surface were analyzed through field-effect scanning electron microscopy imaging and energy-dispersive X-ray spectroscopy. When the discharge energy fell below 2 J, inconsistencies occurred in the specific impulse and the thrust efficiency due to the measurement of the low mass bit. At energy ≥2 J, the performance parameters are compared with other PPT systems of similar configuration and discussed in depth.

Study on Characteristics of a High-Precision Cold Gas Micro Thruster

In order to improve the reliability of the spacecraft micro cold gas propulsion system and realize the precise control of the spacecraft attitude and orbit, a micro-thrust, high-precision cold gas thruster is carried out, at the same time due to the design requirements of the spacecraft, this micro-thrust should be continuous working more than 60 minutes, the traditional solenoid valve used for the thrusts can’t complete the mission, so a long-life micro latching valve is developed as the control valve for this micro thruster, because the micro latching valve can keep its position when it cuts off the outage. Firstly, the authors introduced the design scheme and idea of the thruster. Secondly, the performance of the latching valve and the flow characteristics of the nozzle were simulated. Finally, from the experimental results and compared with the numerical study, it shows that the long-life micro cold gas thruster developed in this paper meets the mission requirements.

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.

Experimental and numerical simulation study of the effect for the anode positions on the discharge characteristics of 300W class low power Hall thrusters

Low-power Hall thruster(LHT) generally has poor discharge efficiency characteristics due to the large surface-to-volume ratio.Aiming to further refine and improve the performance of 300 W class LHT in terms of thrust and efficiency,and to obtain the most optimal operating point,the experimental study of the discharge characteristics for three different anode positions was conducted under the operation of various discharge voltages(100-400 V) and anode mass flow rates(0.65 mg·s-1and 0.95 mg·s-1).The experimental results indicated that the thruster has the most excellent performance in terms of thrust and efficiency etc at a channel length of 27 mm for identical operating conditions.In addition,particle in cell simulations,employed to reveal the underlying physical mechanisms,show that the ionization and acceleration zone is pushed downwards towards the channel exit as the anode moves towards the exit.At the identical operating point,when the channel length is reduced from 32 to 27 mm,the ionization and acceleration zone moves towards the exit,and the parameters such as thrust and efficiency increase due to the high ionization rate,ion number density,and axial electric field.When the channel length is further moved to 24 mm,the parameters in terms of thrust(F) and efficiency(ηa) are reduced as a result of the decreasing ionization efficiency(ηm) and the larger plume divergence angle(α).In this paper,the results indicated that an optimum anode position(ΔL=27 mm) exists for the optimum performance.

Cathode position effects on microwave discharge cusped field thruster

The microwave discharge cusped field thruster is a novel concept of electric micropropulsion device,which operatesμN level thrust in low mass flow rate conditions,making use of a coaxial transmission line resonator.With its advantages of low thrust noise and high thrust resolution over a wide range of thrust,the thruster has emerged as a candidate thruster for the space-borne gravitational wave detection mission.The cathode effects commonly exist in many kinds of electric propulsion,and they are typically significant in micropropulsions.In order to find out the cathode position effects on a microwave discharge cusped field thruster,a thermionic cathode is mounted on a cross-slider for coupling.Under different cathode positions,the plume is analyzed by a Faraday probe and a retarding potential analyzer to analyze the performance and discharge characteristics.The results show that the magnetic mirror effect leads to significant degradation of anode current and an increase in low-energy ion ratio as the cathode moves away from the thruster exit.The electron conduction route also significantly impacts anode current efficiency,related to the cathode-exit distance and the thruster magnetic topology.

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.

Influence of neutral gas supply position on wall erosion of Hall thruster studied by particle-in-cell simulation

In this work,we have carried out a simulation study on the discharge process of Hall thrusters under the conditions of different neutral gas radial supply positions based on the particle-in-cell(PIC)and Monte Carlo collision(MCC)methods.This paper compares the two-dimensional(2D)distributions of neutral gas,plasma and wall erosion-related parameters under different neutral gas supply positions.The comparison results show that the change of the neutral gas supply position affects the radial distribution uniformity of the neutral gas and plasma in the channel.From the comparison of the density peaks,it can be found that the neutral gas density and the plasma density peak under the upper gas supply condition are relatively low,and the plasma density peak is 22.49%lower than the density peak under the middle gas supply condition.Meanwhile,as the radial position of the gas supply moves from the lower gas supply to the upper gas supply,the position of the ionization zone also gradually moves toward the anode.The results of erosion-related parameter distribution comparison show that the change of gas supply location has an obvious influence on erosion rate and erosion range.In terms of erosion rate,the wall erosion rate is relatively low under the upper gas supply condition,and the peak erosion rates of the inner and outer walls are 33.3%and 29.9%lower than those under the other two conditions.In terms of erosion range,as the gas supply position moves from the lower gas supply position to the upper gas supply position,the erosion range gradually increases from5 to 7.5 mm.

Numerical study of viscosity and heat flux role in heavy species dynamics in Hall thruster discharge

A two-and three-dimensional velocity space axisymmetric hybrid-PIC model of Hall thruster discharge called Hybrid2D has been developed.The particle-in-cell(PIC) method was used for neutrals and ions(heavy species),and fluid dynamics on a magnetic field-aligned(MFA) mesh was used for electrons.A time-saving method for heavy species moment interpolation on a MFA mesh was developed.The method comprises using regular rectangle and irregular triangle meshes,connected to each other on a pre-processing stage.The electron fluid model takes into account neither inertia terms nor viscous terms and includes an electron temperature equation with a heat flux term.The developed model was used to calculate all heavy species moments up to the third one in a stationary case.The analysis of the viscosity and the heat flux impact on the force and energy balance has shown that for the calculated geometry of the Hall thruster,the viscosity and the heat flux terms have the same magnitude as the other terms and could not be omitted.Also,it was shown that the heat flux is not proportional to the temperature gradient and,consequently,the highest moments should be calculated to close the neutral fluid equation system.At the same time,ions can only be modeled as a cold non-viscous fluid when the sole aim of modeling is the calculation of the operating parameters or distribution of the local parameters along the centerline of the discharge channel.This is because the magnitude of the viscosity and the temperature gradient terms are negligible at the centerline.However,when a simulation’s focus is either on the radial divergence of the plume or on magnetic pole erosion,three components of the ion temperature should be taken into consideration.The non-diagonal terms of ion pressure tensor have a lower impact than the diagonal terms.According to the study,a zero heat flux condition could be used to close the ion equation system in calculated geometry.

Study of beam divergence and thrust vector eccentricity characteristics of the Hall thruster based on dual Faraday probe array planes and its applications

The accurate knowledge of the thrust vector eccentricity and beam divergence characteristics of Hall thrusters are of significant engineering value for the beneficial integration and successful application of Hall thrusters on spacecraft.For the characteristics of the plume bipolar diffusion due to the annular discharge channel of the Hall thruster,a Gaussian-fitted method for thrust vector deviation angle and beam divergence of Hall thrusters based on dual Faraday probe array planes was proposed in respect of the Hall thruster beam characteristics.The results show that the ratios of the deviation between the maximum and minimum values of the beam divergence angle and the thrust vector eccentricity angle using a Gaussian fit to the optimized Faraday probe dual plane to the mean value are 1.4%and 11.5%,respectively.The optimized thrust vector eccentricity angle obtained has been substantially improved,by approximately 20%.The beam divergence angle calculated using a Gaussian fitting to the optimized Faraday probe dual plane is approximately identical to the non-optimized one.The beam divergence and thrust vector eccentricity angles for different anode mass flow rates were obtained by averaging the beam divergence and thrust vector eccentricity angles calculated by the dual-plane,Gaussian-fitted ion current density method for different cross-sections.The study not only allows for an immediate and effective tool for determining the design of thrust vector adjustment mechanisms of spacecraft with different power Hall thrusters but also for characterizing the 3D spatial distribution of the Hall thruster plume.

Endurance-test and theoretical prediction of a rare earth nanocathode for the applied field magnetoplasmadynamic thruster

The erosion loss of cathode is essential for the lifetime of magnetoplasmadynamic thruster(MPDT).In this work,an endurance test system for MPDT cathodes was designed and developed,and the erosion characteristics,erosion rate and erosion mechanism of the cathode were studied using the system under vacuum condition.The WCe20 hollow cathode was selected to carry out the long-term erosion of 540 h with the argon propellant supply flow rate of40 ml min^(-1),the input current of 25 A,and the central magnetic field intensity of 96 Gs.In order to predict the theoretical service life of cathode,a steady state erosion numerical model was established.The calculation results show that the total erosion rate of sputtering and evaporation is 11.58 mg h^(-1),which is slightly smaller than the test data of the average cathode corrosion rate of 12.70 mg h^(-1) in the experiment,because the experimental value includes start-up erosion rate.

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.

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.