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 re...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.展开更多
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...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.展开更多
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 sys...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.展开更多
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 ana...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.展开更多
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 adap...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.展开更多
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...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.展开更多
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 sen...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.展开更多
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 ...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.展开更多
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 accelarati...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.展开更多
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...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.展开更多
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 ini...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.展开更多
In a 2.45 GHz electron cyclotron resonance(ECR)ion thruster powered with rod antenna under a cross magnetic field,abnormal behaviours such as sudden drop of ion beam current(Ib)and larger increasing-rate of Ibin the h...In a 2.45 GHz electron cyclotron resonance(ECR)ion thruster powered with rod antenna under a cross magnetic field,abnormal behaviours such as sudden drop of ion beam current(Ib)and larger increasing-rate of Ibin the high microwave power(Pw)discharges at high gas flow rates were observed.A differential method was proposed to reveal the changes in the radial profiles of gray values extracted from the end-view discharge images.The increasing-rate of Ibwith respect to Pwwas used to evaluate efficiencies of ion production and transport.Analyses indicate that discharges are dominantly sustained by ordinary wave via electron heating in the electron plasma resonance layer that can shift along the rod-antenna,and extraordinary wave can only ignite a discharge in the ECR layer in the low gas flow rate regime.In terms of the confinement region defined by the magnetic field lines intercepting with the screen grid,the confinement region of the optimized 2.45 GHz cross magnetic field takes the shape of hourglass,enabling the high increasing-rate of Ib with respect to Pw in high power discharges at high gas flow rates.Correlated with the accompanied bright boundary layer appearing in the differentiated image,the sudden drop of Ib in the low gas flow rate regime is attributed to the discharge ignited by the enhanced extraordinary wave in the ECR layer neighbouring the narrowest confinement region,where the produced ions can promptly enter the loss region.展开更多
In this project,we investigated how ion thrusters produce propulsion and how the design of ion thrusters affects the performance of the thruster.In the experiment,we build a high voltage power supply(0-50 kV)and foil ...In this project,we investigated how ion thrusters produce propulsion and how the design of ion thrusters affects the performance of the thruster.In the experiment,we build a high voltage power supply(0-50 kV)and foil rings to produce ion wind.When considering the design of the thruster,we focus on three variables:the volume of the space,where ions are produced and the electric field intensity.Thus,to investigate the first variable we made foil rings with different radius and change the distance between the ring and positive cathode.To determine the propulsion produced we use a speed sensor to determine the magnitude of the wind produced.展开更多
Ion optics is a critical component of ion thrusters.A two-dimensional axisymmetricmodel is developed to study the characteristics of three-grid electron cyclotron resonance ionthruster optics.The code is based on a pa...Ion optics is a critical component of ion thrusters.A two-dimensional axisymmetricmodel is developed to study the characteristics of three-grid electron cyclotron resonance ionthruster optics.The code is based on a particle-in-cell combined with the Monte Carlo collisionmethod to simulate ion dynamics and charge-exchange processes in the grid region.The simulationresults show that the mode can give a reasonable estimate of the physics characteristics of the ionoptics.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 thedecelerator grid plays an important role in reducing 'pits-and-grooves' erosion.展开更多
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 t...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.展开更多
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...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.展开更多
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 direc...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.展开更多
Based on the particle-in-cell (PIC) method, a two-dimensional numerical schemewas developed to investigate the ion beam extraction phenomena through the ion thruster optics.According to the calculated results, the pla...Based on the particle-in-cell (PIC) method, a two-dimensional numerical schemewas 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 electricfield in the calculation domain, and the ion and electron spatial distributions are obtained fordifferent accelerator grid voltages. The results indicate that the accelerator grid voltage affectsthe plasma sheath upstream of the screen grid significantly. It is found that a moderate acceleratorgrid voltage results in an ion optical performance better than either a higher or lower voltage, froma point of ion extraction from the discharge chamber and erosion mitigation of the accelerator griddue to the direct ion impingement.展开更多
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 ...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 GPU'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.展开更多
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 ...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^(18)-8?×?10^(18)m^(-3).The upstream plasma density of the screen grid is in the range 6?×?10^(17)-6?×?10^(18)m^(-3)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^(-14)kg s^(-1)and 4.28?×?10^(-14)kg s^(-1),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^(-13)kg s^(-1)and 4.1?×?10^(-13)kg s^(-1),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^(-13)kg s^(-1)and 7.3?×?10^(-13)kg s^(-1),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.展开更多
基金National Natural Science Foundation of China(No.61901202)Key Laboratory Funds for the Science and Technology on Vacuum Technology and Physics Laboratory,Lanzhou Institute of Physics(No.HTKJ2019KL510003)。
文摘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.
基金supported by Shenzhen Technology Projects(No.ZDSYS201707280904031)the China Postdoctoral Science Foundation(No.2022M710977)。
文摘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.
基金supported by the China Postdoctoral Science Foundation(No.2022M710977)National Natural Science Foundation of China(No.51907039)+1 种基金the Natural Science Foundation of Guangdong Province(Nos.2022A1515110215 and 2023A1515010137)Shenzhen Technology Projects(Nos.JCYJ20190806142603534 and ZDSYS201707280904031)。
文摘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.
基金supported by Key Laboratory Funds for the Science and Technology on Vacuum Technology and Physics Laboratory,Lanzhou Institute of Physics(Nos.HTKJ2022KL510003 and 6142207210303)Independent project of Hangzhou Institute for Advanced Study(No.2022ZZ01009)Science and Technology Project Affiliated to the Education Department of Chongqing Municipality(No.KJZD-K202101506)。
文摘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.
基金supported by Top Young Talents of China Aerospace ScienceTechnology Corporation,Innovation Group Project of Gansu ProvinceNational Natural Science Foundation of China(No.11802111)。
文摘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.
基金financial support from National Natural Science Foundation of China(No.11805265)Key Laboratory of Micro-Satellites,Chinese Academy of Sciences(No.KFKT201903)is gratefully acknowledged。
文摘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.
基金National Natural Science Foundation of China (No. 12005087)the Science and Technology Program of Gansu Province (Nos. 2006ZCTF0054, HTKJ2019KL510003, and 20JR10RA478)。
文摘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.
文摘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.
基金Group independent research and development projects(No.YF-ZZYF-2021-132).
文摘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.
基金the financial support from National Natural Science Foundation of China(Nos.11402025,11475019,and 11702123)the National Key Laboratory of Science and Technology on Vacuum Technology&Physics(No.ZWK1608)+1 种基金the Advanced Space Propulsion Laboratory of BICEBeijing Engineering Research Center of Efficient and Green Aerospace Propulsion Technology(No.Lab ASP-2018-03)。
文摘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.
文摘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.
基金supported by National Natural Science Foundation of China(No.11975070)Open Funds for Science and Technology on Vacuum Technology and Physics Laboratory,Lanzhou Institute of Physics(No.zwk1609)。
文摘In a 2.45 GHz electron cyclotron resonance(ECR)ion thruster powered with rod antenna under a cross magnetic field,abnormal behaviours such as sudden drop of ion beam current(Ib)and larger increasing-rate of Ibin the high microwave power(Pw)discharges at high gas flow rates were observed.A differential method was proposed to reveal the changes in the radial profiles of gray values extracted from the end-view discharge images.The increasing-rate of Ibwith respect to Pwwas used to evaluate efficiencies of ion production and transport.Analyses indicate that discharges are dominantly sustained by ordinary wave via electron heating in the electron plasma resonance layer that can shift along the rod-antenna,and extraordinary wave can only ignite a discharge in the ECR layer in the low gas flow rate regime.In terms of the confinement region defined by the magnetic field lines intercepting with the screen grid,the confinement region of the optimized 2.45 GHz cross magnetic field takes the shape of hourglass,enabling the high increasing-rate of Ib with respect to Pw in high power discharges at high gas flow rates.Correlated with the accompanied bright boundary layer appearing in the differentiated image,the sudden drop of Ib in the low gas flow rate regime is attributed to the discharge ignited by the enhanced extraordinary wave in the ECR layer neighbouring the narrowest confinement region,where the produced ions can promptly enter the loss region.
文摘In this project,we investigated how ion thrusters produce propulsion and how the design of ion thrusters affects the performance of the thruster.In the experiment,we build a high voltage power supply(0-50 kV)and foil rings to produce ion wind.When considering the design of the thruster,we focus on three variables:the volume of the space,where ions are produced and the electric field intensity.Thus,to investigate the first variable we made foil rings with different radius and change the distance between the ring and positive cathode.To determine the propulsion produced we use a speed sensor to determine the magnitude of the wind produced.
文摘Ion optics is a critical component of ion thrusters.A two-dimensional axisymmetricmodel is developed to study the characteristics of three-grid electron cyclotron resonance ionthruster optics.The code is based on a particle-in-cell combined with the Monte Carlo collisionmethod to simulate ion dynamics and charge-exchange processes in the grid region.The simulationresults show that the mode can give a reasonable estimate of the physics characteristics of the ionoptics.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 thedecelerator grid plays an important role in reducing 'pits-and-grooves' erosion.
基金supported by National Natural Science Foundation of China (No. 11702123)Technology Based Research Projects of SASTIND (No. JSZL2017203B008)Fund for Distinguished Young Talents of CAST (No. 9140C550206130C55003)
文摘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.
基金supported by National Natural Science Foundation of China under Grant No. 11702123
文摘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.
文摘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.
基金supported by the China Postdoctoral Science Foundation (CPSF) (No. 20090450275)National Natural Science Foundation of China (No. 10805004)
文摘Based on the particle-in-cell (PIC) method, a two-dimensional numerical schemewas 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 electricfield in the calculation domain, and the ion and electron spatial distributions are obtained fordifferent accelerator grid voltages. The results indicate that the accelerator grid voltage affectsthe plasma sheath upstream of the screen grid significantly. It is found that a moderate acceleratorgrid voltage results in an ion optical performance better than either a higher or lower voltage, froma point of ion extraction from the discharge chamber and erosion mitigation of the accelerator griddue to the direct ion impingement.
基金supported by National Natural Science Foundation of China (No. 10805004)Foundation of National Key Lab. of Science and Technology on Vacuum & Cryogenic of China (No. 9140C550404100C55)
文摘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 GPU'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.
基金supported by the fund of National Key Laboratory of Science and Technology on Vacuum Technology and Physics (Grant No. 6142207030103)
文摘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^(18)-8?×?10^(18)m^(-3).The upstream plasma density of the screen grid is in the range 6?×?10^(17)-6?×?10^(18)m^(-3)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^(-14)kg s^(-1)and 4.28?×?10^(-14)kg s^(-1),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^(-13)kg s^(-1)and 4.1?×?10^(-13)kg s^(-1),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^(-13)kg s^(-1)and 7.3?×?10^(-13)kg s^(-1),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.
