Based on the analysis of the physical mechanism of the Stationary Plasma Thruster （SPT）, an integral equation describing the ion density of the steady SPT and the ion velocity distribution function at an arbitrary a...Based on the analysis of the physical mechanism of the Stationary Plasma Thruster （SPT）, an integral equation describing the ion density of the steady SPT and the ion velocity distribution function at an arbitrary axial position of the steady SPT channel are derived. The integral equation is equivalent to the Vlasov equation, but the former is simpler than the latter. A one dimensional steady quasineutral hybrid model is established. In this model, ions are described by the above integral equation, and neutrals and electrons are described by hydrodynamic equations. The transferred equivalency to the differential equation and the integral equation, together with other equations, are solved by an ordinary differential equation （ODE） solver in the Matlab. The numerical simulation results show that under various circumstances, the ion average velocity would be different and needs to be deduced separately.展开更多
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.展开更多
Ablation excited by current pulses is a very critical physical process in pulse plasma thrusters(PPT).Its effects on wall-plasma interaction directly determine the PPT performances.In order to reveal the process of th...Ablation excited by current pulses is a very critical physical process in pulse plasma thrusters(PPT).Its effects on wall-plasma interaction directly determine the PPT performances.In order to reveal the process of the ablated wall interaction with the discharge plasma in PPT,ablation models formulated by three different boundary conditions at the wall-plasma interface are studied.These are the two widely used high-speed evaporation models(Model-L and Model-M),and the recently developed Keida-Zaghloul model(Model-K)of the Knudsen layer that takes into account the internal degrees of freedom on the energy flux conservation.First,fundamental mechanisms of the three ablation models are clarified by comparative analysis in order to gain a comprehensive understanding of the wall-plasma interaction.Then,the applicability of different ablation models with the numerical solutions of LES-6 PPT is investigated in detail using magnetohydrodynamic(MHD)modeling.Results show that Model-L and Model-M are actually special cases of Model-K when a simplified jump conditions limited by high velocity at the vapor/plasma interface is used;A ratio of ablation rate in Model-L to that in Model-M is about 0.8at the same wall surface temperature,while it rises to 1 at different surface temperature determined by Model-L and Model-M in PPT.Even though Model-K solution requires significant computational time,it shows more accurate ablation feature for the wall-plasma interaction and possesses better computing precision of impulse bit during post-pulse which is useful for future studies of the late time ablation.展开更多
基金The project supported by National Fundamental Science Research Fundation of China (No. K1403060719)
文摘Based on the analysis of the physical mechanism of the Stationary Plasma Thruster （SPT）, an integral equation describing the ion density of the steady SPT and the ion velocity distribution function at an arbitrary axial position of the steady SPT channel are derived. The integral equation is equivalent to the Vlasov equation, but the former is simpler than the latter. A one dimensional steady quasineutral hybrid model is established. In this model, ions are described by the above integral equation, and neutrals and electrons are described by hydrodynamic equations. The transferred equivalency to the differential equation and the integral equation, together with other equations, are solved by an ordinary differential equation （ODE） solver in the Matlab. The numerical simulation results show that under various circumstances, the ion average velocity would be different and needs to be deduced separately.
基金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.
基金Project supported by Ph.D.Programs Foundation of Ministry of Education of China(20121101120004)Basic Research Foundation of Beijing Institute of Technology(20120142015)
文摘Ablation excited by current pulses is a very critical physical process in pulse plasma thrusters(PPT).Its effects on wall-plasma interaction directly determine the PPT performances.In order to reveal the process of the ablated wall interaction with the discharge plasma in PPT,ablation models formulated by three different boundary conditions at the wall-plasma interface are studied.These are the two widely used high-speed evaporation models(Model-L and Model-M),and the recently developed Keida-Zaghloul model(Model-K)of the Knudsen layer that takes into account the internal degrees of freedom on the energy flux conservation.First,fundamental mechanisms of the three ablation models are clarified by comparative analysis in order to gain a comprehensive understanding of the wall-plasma interaction.Then,the applicability of different ablation models with the numerical solutions of LES-6 PPT is investigated in detail using magnetohydrodynamic(MHD)modeling.Results show that Model-L and Model-M are actually special cases of Model-K when a simplified jump conditions limited by high velocity at the vapor/plasma interface is used;A ratio of ablation rate in Model-L to that in Model-M is about 0.8at the same wall surface temperature,while it rises to 1 at different surface temperature determined by Model-L and Model-M in PPT.Even though Model-K solution requires significant computational time,it shows more accurate ablation feature for the wall-plasma interaction and possesses better computing precision of impulse bit during post-pulse which is useful for future studies of the late time ablation.
