A two-dimensional PIC/MCC model is developed to simulate the nitrogen radio frequency hollow cathode discharge（rf-HCD）.It is found that both the sheath oscillation heating and the secondary electron heating together...A two-dimensional PIC/MCC model is developed to simulate the nitrogen radio frequency hollow cathode discharge（rf-HCD）.It is found that both the sheath oscillation heating and the secondary electron heating together play a role to maintain the rf-HCD under the simulated conditions.The mean energy of ions（N＋_2,N＋）in the negative glow region is greater than the thermal kinetic energy of the molecular gas（N2）,which is an important characteristic of rf-HCD.During the negative portion of the hollow electrode voltage cycle,electrons mainly follow pendulum movement and produce a large number of ionization collisions in the plasma region.During the positive voltage of the rf cycle,the axial electric field becomes stronger and its direction is pointing to the anode（substrate）,therefore the ions move toward the anode（substrate）via the axial electric field acceleration.Compared with dc-HCD,rf-HCD is more suitable for serving as a plasma jet nozzle at low pressure.展开更多
Numerical simulations by the code of Object-Oriented PIC （Particle-in-Cell） and the Monte Carlo Collision （MCC） method were carried out in order to obtain an insight into the characteristics of plasmas generated b...Numerical simulations by the code of Object-Oriented PIC （Particle-in-Cell） and the Monte Carlo Collision （MCC） method were carried out in order to obtain an insight into the characteristics of plasmas generated by glow discharges in low pressure helium in a four-anode DC glow discharge device. The results show that, the pressure, the external mirror magnetic field, and the virtual breadth of the annular electrode affect the radial distribution of the plasma density and temperature. The simulations are instructive for further experiments.展开更多
Radio frequency capacitively coupled plasma source(RF-CCP)with a hollow electrode can increase the electron density through the hollow cathode effect(HCE),which offers a method to modify the spatial profiles of the pl...Radio frequency capacitively coupled plasma source(RF-CCP)with a hollow electrode can increase the electron density through the hollow cathode effect(HCE),which offers a method to modify the spatial profiles of the plasma density.In this work,the variations of the HCE in one RF period are investigated by using a two-dimensional particle-in-cell/Monte-Carlo collision(PIC/MCC)model.The results show that the sheath electric field,the sheath potential drop,the sheath thickness,the radial plasma bulk width,the electron energy distribution function(EEDF),and the average electron energy in the cavity vary in one RF period.During the hollow electrode sheath's expansion phase,the secondary electron heating and sheath oscillation heating in the cavity are gradually enhanced,and the frequency of the electron pendular motion in the cavity gradually increases,hence the HCE is gradually enhanced.However,during the hollow electrode sheath's collapse phase,the secondary electron heating is gradually attenuated.In addition,when interacting with the gradually collapsed hollow electrode sheaths,high-energy plasma bulk electrons in the cavity will lose some energy.Furthermore,the frequency of the electron pendular motion in the cavity gradually decreases.Therefore,during the hollow electrode sheath's collapse phase,the HCE is gradually attenuated.展开更多
为了研究火花试验装置中电极在真空中放电的微观特性,本文建立了在真空环境下,以钨为阳极材料、镉为阴极材料的二维平行板放电模型.采用PIC/MCC(Particle-In-Cell/Monte Carlo Collision)方法对该模型进行仿真,研究了不同电子发射机制...为了研究火花试验装置中电极在真空中放电的微观特性,本文建立了在真空环境下,以钨为阳极材料、镉为阴极材料的二维平行板放电模型.采用PIC/MCC(Particle-In-Cell/Monte Carlo Collision)方法对该模型进行仿真,研究了不同电子发射机制下平行板电极放电的发展过程以及空间场强、阴极表面温度和场增强因子对空间电子变化的影响,得到在场致发射、热发射以及热-场致发射作用下放电过程中的电子浓度和阳极吸收电流的变化以及电子密度和电势的空间分布等.研究发现,场致发射是微间隙阴极电子发射的主导发射机制,当阴极表面温度在焦耳热的作用下达到镉金属的沸点1040K时将产生镉蒸汽,电流密度和电子浓度逐渐增大,此时热发射将开始作用于微间隙放电;当温度大于镉金属气化温度后,场强的影响将大于温度的影响;当场增强因子很小时,热发射几乎不起作用,随着场增强因子不断增大,热发射的作用逐渐增强,导致空间电子浓度明显增加,真空环境下微间隙放电是由热-场共同作用的.展开更多
电子反流失效模式是离子推力器关键失效模式之一,决定推力器工作寿命。为明确各参数对电子反流失效模式的影响程度,确定加速应力,为地面加速寿命实验验证方案和长寿命优化设计提供数据支持,采用Hybrid-PIC-MCC(Particle in Cell-Monte C...电子反流失效模式是离子推力器关键失效模式之一,决定推力器工作寿命。为明确各参数对电子反流失效模式的影响程度,确定加速应力,为地面加速寿命实验验证方案和长寿命优化设计提供数据支持,采用Hybrid-PIC-MCC(Particle in Cell-Monte Carlo Collision)方法,构建了三栅极系统数值仿真模型。采用模型研究了地面真空舱本底压力、屏栅电压、加速栅电压、屏栅与加速栅间距、屏栅上游等离子体密度和放电室工质利用率等参数的影响敏感度对比。研究结果显示,真空舱本底压力可以作为加速寿命试验的首选加速应力,在推力器结构和工作本征参数中工质利用率为最敏感应力,其次是屏栅电压、屏栅上游等离子体密度、加速栅电压、屏栅和加速栅间距。展开更多
In order to study the influence of three-grid assembly thermal deformation caused by heat accumulation on breakdown times and an ion extraction process,a hot gap test and a breakdown time test are carried out to obtai...In order to study the influence of three-grid assembly thermal deformation caused by heat accumulation on breakdown times and an ion extraction process,a hot gap test and a breakdown time test are carried out to obtain thermal deformation of the grids when the thruster is in 5 k W operation mode.Meanwhile,the fluid simulation method and particle-in-cell-Monte Carlo collision（PICMCC） method are adopted to simulate the ion extraction process according to the previous test results.The numerical calculation results are verified by the ion thruster performance test.The results show that after about 1.2 h operation,the hot gap between the screen grid and the accelerator grid reduce to 0.25–0.3 mm,while the hot gap between the accelerator grid and the decelerator grid increase from 1 mm to about 1.4 mm when the grids reach thermal equilibrium,and the hot gap is almost unchanged.In addition,the breakdown times experiment shows that 0.26 mm is the minimal safe hot gap for the grid assembly as the breakdown times improves significantly when the gap is smaller than this value.Fluid simulation results show that the plasma density of the screen grid is in the range 6？×10^（17）–6？×？10^（18） m^（13） and displays a parabolic characteristic,while the electron temperature gradually increases along the axial direction.The PIC-MCC results show that the current falling of an ion beam through a single aperture is significant.Meanwhile,the intercepted current of the accelerator grid and the decelerator grid both increase with the change in the hot gap.The ion beam current has optimal perveance status without thermal deformation,and the intercepted current of the accelerator grid and the decelerator grid are 3.65 m A and 6.26 m A,respectively.Furthermore,under the effect of thermal deformation,the ion beam current has over-perveance status,and the intercepted current of the accelerator grid and the decelerator grid are 10.46 m A and 18.24 m A,respectively.Performance test results indicate that the breakdown times increase obviously.The intercepted current of the accelerator grid and the decelerator grid increases to 13 m A and 16.5 m A,respectively,due to the change in the hot gap after 1.5 h operation.The numerical calculation results are well consistent with performance test results,and the error comes mainly from the test uncertainty of the hot gap.展开更多
The electron energy distribution function （EEDF）, predicted by the Boltzmann equation solver BOLSIG＋ based on the two-term approximation, is introduced into the fluid model for simulating the high-power microwave ...The electron energy distribution function （EEDF）, predicted by the Boltzmann equation solver BOLSIG＋ based on the two-term approximation, is introduced into the fluid model for simulating the high-power microwave （HPM） breakdown in argon, nitrogen, and air, and its validity is examined by comparing with the results of particle-in-cell Monte Carlo collision （PIC/MCC） simulations as well as the experimental data. Numerical results show that, the breakdown time of the fluid model with the Maxwellian EEDF matches that of the PIC/MCC simulations in nitrogen; however, in argon under high pressures, the results from the Maxwellian EEDF were poor. This is due to an overestimation of the energy tail of the Maxwellian EEDF in argon breakdown. The prediction of the fluid model with the BOLSIG＋ EEDF, however, agrees very well with the PIC/MCC prediction in nitrogen and argon over a wide range of pressures. The accuracy of the fluid model with the BOLSIG＋ EEDF is also verified by the experimental results of the air breakdown.展开更多
基金supported by Natural Science Foundation of Hebei Province,China(No.A2012205072)
文摘A two-dimensional PIC/MCC model is developed to simulate the nitrogen radio frequency hollow cathode discharge（rf-HCD）.It is found that both the sheath oscillation heating and the secondary electron heating together play a role to maintain the rf-HCD under the simulated conditions.The mean energy of ions（N＋_2,N＋）in the negative glow region is greater than the thermal kinetic energy of the molecular gas（N2）,which is an important characteristic of rf-HCD.During the negative portion of the hollow electrode voltage cycle,electrons mainly follow pendulum movement and produce a large number of ionization collisions in the plasma region.During the positive voltage of the rf cycle,the axial electric field becomes stronger and its direction is pointing to the anode（substrate）,therefore the ions move toward the anode（substrate）via the axial electric field acceleration.Compared with dc-HCD,rf-HCD is more suitable for serving as a plasma jet nozzle at low pressure.
文摘Numerical simulations by the code of Object-Oriented PIC （Particle-in-Cell） and the Monte Carlo Collision （MCC） method were carried out in order to obtain an insight into the characteristics of plasmas generated by glow discharges in low pressure helium in a four-anode DC glow discharge device. The results show that, the pressure, the external mirror magnetic field, and the virtual breadth of the annular electrode affect the radial distribution of the plasma density and temperature. The simulations are instructive for further experiments.
文摘Radio frequency capacitively coupled plasma source(RF-CCP)with a hollow electrode can increase the electron density through the hollow cathode effect(HCE),which offers a method to modify the spatial profiles of the plasma density.In this work,the variations of the HCE in one RF period are investigated by using a two-dimensional particle-in-cell/Monte-Carlo collision(PIC/MCC)model.The results show that the sheath electric field,the sheath potential drop,the sheath thickness,the radial plasma bulk width,the electron energy distribution function(EEDF),and the average electron energy in the cavity vary in one RF period.During the hollow electrode sheath's expansion phase,the secondary electron heating and sheath oscillation heating in the cavity are gradually enhanced,and the frequency of the electron pendular motion in the cavity gradually increases,hence the HCE is gradually enhanced.However,during the hollow electrode sheath's collapse phase,the secondary electron heating is gradually attenuated.In addition,when interacting with the gradually collapsed hollow electrode sheaths,high-energy plasma bulk electrons in the cavity will lose some energy.Furthermore,the frequency of the electron pendular motion in the cavity gradually decreases.Therefore,during the hollow electrode sheath's collapse phase,the HCE is gradually attenuated.
文摘为了研究火花试验装置中电极在真空中放电的微观特性,本文建立了在真空环境下,以钨为阳极材料、镉为阴极材料的二维平行板放电模型.采用PIC/MCC(Particle-In-Cell/Monte Carlo Collision)方法对该模型进行仿真,研究了不同电子发射机制下平行板电极放电的发展过程以及空间场强、阴极表面温度和场增强因子对空间电子变化的影响,得到在场致发射、热发射以及热-场致发射作用下放电过程中的电子浓度和阳极吸收电流的变化以及电子密度和电势的空间分布等.研究发现,场致发射是微间隙阴极电子发射的主导发射机制,当阴极表面温度在焦耳热的作用下达到镉金属的沸点1040K时将产生镉蒸汽,电流密度和电子浓度逐渐增大,此时热发射将开始作用于微间隙放电;当温度大于镉金属气化温度后,场强的影响将大于温度的影响;当场增强因子很小时,热发射几乎不起作用,随着场增强因子不断增大,热发射的作用逐渐增强,导致空间电子浓度明显增加,真空环境下微间隙放电是由热-场共同作用的.
文摘电子反流失效模式是离子推力器关键失效模式之一,决定推力器工作寿命。为明确各参数对电子反流失效模式的影响程度,确定加速应力,为地面加速寿命实验验证方案和长寿命优化设计提供数据支持,采用Hybrid-PIC-MCC(Particle in Cell-Monte Carlo Collision)方法,构建了三栅极系统数值仿真模型。采用模型研究了地面真空舱本底压力、屏栅电压、加速栅电压、屏栅与加速栅间距、屏栅上游等离子体密度和放电室工质利用率等参数的影响敏感度对比。研究结果显示,真空舱本底压力可以作为加速寿命试验的首选加速应力,在推力器结构和工作本征参数中工质利用率为最敏感应力,其次是屏栅电压、屏栅上游等离子体密度、加速栅电压、屏栅和加速栅间距。
基金supported by the National Key Laboratory Fund of Science and Technology on Vacuum Technology & Physics (Grant No.6142207030103)National Natural Science Foundation of China (Grant No.11702123)
文摘In order to study the influence of three-grid assembly thermal deformation caused by heat accumulation on breakdown times and an ion extraction process,a hot gap test and a breakdown time test are carried out to obtain thermal deformation of the grids when the thruster is in 5 k W operation mode.Meanwhile,the fluid simulation method and particle-in-cell-Monte Carlo collision（PICMCC） method are adopted to simulate the ion extraction process according to the previous test results.The numerical calculation results are verified by the ion thruster performance test.The results show that after about 1.2 h operation,the hot gap between the screen grid and the accelerator grid reduce to 0.25–0.3 mm,while the hot gap between the accelerator grid and the decelerator grid increase from 1 mm to about 1.4 mm when the grids reach thermal equilibrium,and the hot gap is almost unchanged.In addition,the breakdown times experiment shows that 0.26 mm is the minimal safe hot gap for the grid assembly as the breakdown times improves significantly when the gap is smaller than this value.Fluid simulation results show that the plasma density of the screen grid is in the range 6？×10^（17）–6？×？10^（18） m^（13） and displays a parabolic characteristic,while the electron temperature gradually increases along the axial direction.The PIC-MCC results show that the current falling of an ion beam through a single aperture is significant.Meanwhile,the intercepted current of the accelerator grid and the decelerator grid both increase with the change in the hot gap.The ion beam current has optimal perveance status without thermal deformation,and the intercepted current of the accelerator grid and the decelerator grid are 3.65 m A and 6.26 m A,respectively.Furthermore,under the effect of thermal deformation,the ion beam current has over-perveance status,and the intercepted current of the accelerator grid and the decelerator grid are 10.46 m A and 18.24 m A,respectively.Performance test results indicate that the breakdown times increase obviously.The intercepted current of the accelerator grid and the decelerator grid increases to 13 m A and 16.5 m A,respectively,due to the change in the hot gap after 1.5 h operation.The numerical calculation results are well consistent with performance test results,and the error comes mainly from the test uncertainty of the hot gap.
基金Project supported by the National Basic Research Program of China(Grant No.2013CB328904)the Fundamental Research Funds for the Central Universities,Chinathe Open Research Fund of Key Laboratory of Cognitive Radio and Information Processing of Ministry of Education of China
文摘The electron energy distribution function （EEDF）, predicted by the Boltzmann equation solver BOLSIG＋ based on the two-term approximation, is introduced into the fluid model for simulating the high-power microwave （HPM） breakdown in argon, nitrogen, and air, and its validity is examined by comparing with the results of particle-in-cell Monte Carlo collision （PIC/MCC） simulations as well as the experimental data. Numerical results show that, the breakdown time of the fluid model with the Maxwellian EEDF matches that of the PIC/MCC simulations in nitrogen; however, in argon under high pressures, the results from the Maxwellian EEDF were poor. This is due to an overestimation of the energy tail of the Maxwellian EEDF in argon breakdown. The prediction of the fluid model with the BOLSIG＋ EEDF, however, agrees very well with the PIC/MCC prediction in nitrogen and argon over a wide range of pressures. The accuracy of the fluid model with the BOLSIG＋ EEDF is also verified by the experimental results of the air breakdown.
