The E×B drift instability is studied in Hall thruster using one-dimensional particle in cell(PIC)simulation method.By using the dispersion relation,it is found that unstable modes occur only in discrete bands in ...The E×B drift instability is studied in Hall thruster using one-dimensional particle in cell(PIC)simulation method.By using the dispersion relation,it is found that unstable modes occur only in discrete bands in k space at cyclotron harmonics.The results indicate that the number of unstable modes increases by increasing the external electric field and decreases by increasing the radial magnetic field.The ion mass does not affect the instability wavelength.Furthermore,the results confirm that there is an instability with short wavelength and high frequency.Finally,it is shown that the electron and ion distribution functions deviate from the initial state and eventually the instability is saturated by ion trapping in the azimuthal direction.Also for light mass ion,the frequency and phase velocity are very high that could lead to high electron mobility in the axial direction.展开更多
In order to understand the physical mechanism of multipactor discharge on dielectric window surface under high power microwave （HPM） excitation in vacuum, an electron movement simulation model based on the particle-...In order to understand the physical mechanism of multipactor discharge on dielectric window surface under high power microwave （HPM） excitation in vacuum, an electron movement simulation model based on the particle-in-cell （PIC） Monte Carlo （MC） is built in this paper. The influences of microwave electromagnetic field and electrostatic field from dielectric surface charging are simultaneously considered in this model. During the simulation, the emission velocity and angle distribution of secondary electrons from the dielectric surface are taken into account. The movement trajectories of electron clusters under complex field excitation are obtained. The influences of emergence angle and microwave electromagnetic parameters on the electron movement are analyzed. It is found that the emergence angle of electrons from the surface has significant effect on its movement, and both the impact energy and return time of electrons oscillate periodically with the phase of microwave field. The number of secondary electrons and induced electrostatic field from multipactoring are also investigated. The results reveal that both values oscillate periodically at twice the microwave frequency, which is due to the electron impact energy oscillating with microwave period. A schematic diagram is proposed to explain the periodical oscillation phenomena.展开更多
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.展开更多
In this study,we computationally examined the dynamics of dielectric barrier discharge in hydrogen sulfide.The simulations were performed with a 1 d3 v particle-in-cell/Monte Carlo collision model in which a parallel-...In this study,we computationally examined the dynamics of dielectric barrier discharge in hydrogen sulfide.The simulations were performed with a 1 d3 v particle-in-cell/Monte Carlo collision model in which a parallel-plate electrode geometry with dielectrics was used.Particle recombination process is represented in the model.The discharge mode was found to be initially Townsend discharge developing from the cathode to the anode,and at the peak of the current,a more stable glow discharge develops from the anode to the cathode.A higher applied voltage results in sufficient secondary electrons to trigger a second current peak,and then the current amplitude increases.As the frequency is increased,it leads to the advance of the phase and an increase in the amplitude of the current peak.A higher dielectric permittivity also makes the discharge occur earlier and more violently in the gap.展开更多
The main ion-atomic collision treatment methods based on Monte-Carlo simulation are considered and discussed. We have proposed an efficient scheme for simulation of time between collisions taking into account cross-se...The main ion-atomic collision treatment methods based on Monte-Carlo simulation are considered and discussed. We have proposed an efficient scheme for simulation of time between collisions taking into account cross-section dependence on ion velocity and random generation of ion velocities and scattering angles after collisions. The developed algorithm of simulation of interval between collisions takes into account the change of relative velocity of ion-atom pair as well as the change of cross-section of collision and atomic concentration. At the same time, unlike the widely used “null-collision” method, both the probability of collision and change of particles’ state which determines this probability are taken into consideration for each particle independently in time. The simulation results according to the techniques proposed are found to be close to the theoretical values of ion drift velocities. It is revealed that the “null-collision” method results in exceeding of drift velocity in strong and intermediate fields. At the same time the proposed method of accumulation of probability under the same conditions gives values close to theoretical ones. In weak fields calculated values of drift velocity in both methods exceed theoretical values to some small extent.展开更多
Plasma in the discharge channel of a pulsed plasma thruster(PPT) with flared electrodes is simulated by a self-developed two-dimensional code. The fully particle-in-cell method with Monte Carlo collision is employed t...Plasma in the discharge channel of a pulsed plasma thruster(PPT) with flared electrodes is simulated by a self-developed two-dimensional code. The fully particle-in-cell method with Monte Carlo collision is employed to model the particle movement and collisions and investigate the plasma properties and acceleration process. Temporal and spatial variations of the electron density distribution and the ion velocity between electrodes are calculated and analyzed in detail.The computational results of the electron number density, which is in the order of 1023 m-3,show good agreements with experimental results of a PPT named ADD SIMP-LEX. The ion velocity distributions along the center line of the channel lead to a comprehensive understanding of ions accelerated by electromagnetic field. The electron distributions of PPT with discharge voltages varying from 1300 to 2000 V are compared. The diffusion of electrons presents strong dependency on discharge voltage and implies higher degree of ionization for higher voltage.展开更多
Particle-in-cell (PIC) method has got much benefits from GPU-accelerated heterogeneous systems.However,the performance of PIC is constrained by the interpolation operations in the weighting process on GPU (graphic pro...Particle-in-cell (PIC) method has got much benefits from GPU-accelerated heterogeneous systems.However,the performance of PIC is constrained by the interpolation operations in the weighting process on GPU (graphic processing unit).Aiming at this problem,a fast weighting method for PIC simulation on GPU-accelerated systems was proposed to avoid the atomic memory operations during the weighting process.The method was implemented by taking advantage of GPU's thread synchronization mechanism and dividing the problem space properly.Moreover,software managed shared memory on the GPU was employed to buffer the intermediate data.The experimental results show that the method achieves speedups up to 3.5 times compared to previous works,and runs 20.08 times faster on one NVIDIA Tesla M2090 GPU compared to a single core of Intel Xeon X5670 CPU.展开更多
In this paper, we analytically explore the magnetic field and mass density evolutions obtained in particle-in-cell（PIC）and magnetohydrodynamics（MHD） simulations of a rarefied deuterium shell Z-pinch and compare th...In this paper, we analytically explore the magnetic field and mass density evolutions obtained in particle-in-cell（PIC）and magnetohydrodynamics（MHD） simulations of a rarefied deuterium shell Z-pinch and compare those results, and also we study the effects of artificially increased Spitzer resistivity on the magnetic field evolution and Z-pinch dynamic process in the MHD simulation. There are significant differences between the profiles of mass density in the PIC and MHD simulations before 45 ns of the Z-pinch in this study. However, after the shock formation in the PIC simulation,the mass density profile is similar to that in the MHD simulation in the case of using multiplier 2 to modify the Spitzer resistivity. Compared with the magnetic field profiles of the PIC simulation of the shell, the magnetic field diffusion has still not been sufficiently revealed in the MHD simulation even though their convergence ratios become the same by using larger multipliers in the resistivity. The MHD simulation results suggest that the magnetic field diffusion is greatly enhanced by increasing the Spitzer resistivity used, which, however, causes the implosion characteristic to change from shock compression to weak shock, even shockless evolution, and expedites the expansion of the shell. Too large a multiplier is not suggested to be used to modify the resistivity in some Z-pinch applications, such as the Z-pinch driven inertial confinement fusion（ICF） in a dynamic hohlraum. Two-fluid or Hall MHD model, even the PIC/fluid hybrid simulation would be considered as a suitable physical model when there exist the plasma regions with very low density in the simulated domain.展开更多
A rainstorm caused by mesoscale convective system (MCS) in Guizhou Province in June 25-26 in 2005 was simulated with the MM5 model. Based on the good simulated results of the MCS developing and the clouds physics proc...A rainstorm caused by mesoscale convective system (MCS) in Guizhou Province in June 25-26 in 2005 was simulated with the MM5 model. Based on the good simulated results of the MCS developing and the clouds physics process, and by means of reducing the height of Yunnan-Guizhou Plateau and cutting off the middle-east of the Yunnan-Guizhou Plateau on the simulated tests, the question as how the ladder terrain on the west of Yunnan-Guizhou Plateau impact on the rainstorm of Guizhou was studied. The analysis results showed that the second ladder terrain of Yunnan-Guizhou Plateau only affected the development of convective clouds on its backward position,and hardly affected the rain on its upward. The whole terrain of the Yunnan-Guizhou Plateau had a distinct impact not only on the windward slope rainfall of the west of the plateau, but also on the rainfall distribution, intensity and continuing time of the convective clouds on the middle-east of the plateau.展开更多
文摘The E×B drift instability is studied in Hall thruster using one-dimensional particle in cell(PIC)simulation method.By using the dispersion relation,it is found that unstable modes occur only in discrete bands in k space at cyclotron harmonics.The results indicate that the number of unstable modes increases by increasing the external electric field and decreases by increasing the radial magnetic field.The ion mass does not affect the instability wavelength.Furthermore,the results confirm that there is an instability with short wavelength and high frequency.Finally,it is shown that the electron and ion distribution functions deviate from the initial state and eventually the instability is saturated by ion trapping in the azimuthal direction.Also for light mass ion,the frequency and phase velocity are very high that could lead to high electron mobility in the axial direction.
基金supported in part by the National High Technology Research and Development Program of China
文摘In order to understand the physical mechanism of multipactor discharge on dielectric window surface under high power microwave （HPM） excitation in vacuum, an electron movement simulation model based on the particle-in-cell （PIC） Monte Carlo （MC） is built in this paper. The influences of microwave electromagnetic field and electrostatic field from dielectric surface charging are simultaneously considered in this model. During the simulation, the emission velocity and angle distribution of secondary electrons from the dielectric surface are taken into account. The movement trajectories of electron clusters under complex field excitation are obtained. The influences of emergence angle and microwave electromagnetic parameters on the electron movement are analyzed. It is found that the emergence angle of electrons from the surface has significant effect on its movement, and both the impact energy and return time of electrons oscillate periodically with the phase of microwave field. The number of secondary electrons and induced electrostatic field from multipactoring are also investigated. The results reveal that both values oscillate periodically at twice the microwave frequency, which is due to the electron impact energy oscillating with microwave period. A schematic diagram is proposed to explain the periodical oscillation phenomena.
基金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.
基金supported by National Natural Science Foundation of China(Nos.21503279,51777164,and 52077169)the Fundamental Research Funds for the Central Universities,China(Nos.xtr042019009 and PY3A083)HPC Platform,Xi’an Jiaotong University。
文摘In this study,we computationally examined the dynamics of dielectric barrier discharge in hydrogen sulfide.The simulations were performed with a 1 d3 v particle-in-cell/Monte Carlo collision model in which a parallel-plate electrode geometry with dielectrics was used.Particle recombination process is represented in the model.The discharge mode was found to be initially Townsend discharge developing from the cathode to the anode,and at the peak of the current,a more stable glow discharge develops from the anode to the cathode.A higher applied voltage results in sufficient secondary electrons to trigger a second current peak,and then the current amplitude increases.As the frequency is increased,it leads to the advance of the phase and an increase in the amplitude of the current peak.A higher dielectric permittivity also makes the discharge occur earlier and more violently in the gap.
文摘The main ion-atomic collision treatment methods based on Monte-Carlo simulation are considered and discussed. We have proposed an efficient scheme for simulation of time between collisions taking into account cross-section dependence on ion velocity and random generation of ion velocities and scattering angles after collisions. The developed algorithm of simulation of interval between collisions takes into account the change of relative velocity of ion-atom pair as well as the change of cross-section of collision and atomic concentration. At the same time, unlike the widely used “null-collision” method, both the probability of collision and change of particles’ state which determines this probability are taken into consideration for each particle independently in time. The simulation results according to the techniques proposed are found to be close to the theoretical values of ion drift velocities. It is revealed that the “null-collision” method results in exceeding of drift velocity in strong and intermediate fields. At the same time the proposed method of accumulation of probability under the same conditions gives values close to theoretical ones. In weak fields calculated values of drift velocity in both methods exceed theoretical values to some small extent.
基金supported by National Natural Science Foundation of China (Grant No. 11602016)
文摘Plasma in the discharge channel of a pulsed plasma thruster(PPT) with flared electrodes is simulated by a self-developed two-dimensional code. The fully particle-in-cell method with Monte Carlo collision is employed to model the particle movement and collisions and investigate the plasma properties and acceleration process. Temporal and spatial variations of the electron density distribution and the ion velocity between electrodes are calculated and analyzed in detail.The computational results of the electron number density, which is in the order of 1023 m-3,show good agreements with experimental results of a PPT named ADD SIMP-LEX. The ion velocity distributions along the center line of the channel lead to a comprehensive understanding of ions accelerated by electromagnetic field. The electron distributions of PPT with discharge voltages varying from 1300 to 2000 V are compared. The diffusion of electrons presents strong dependency on discharge voltage and implies higher degree of ionization for higher voltage.
基金Projects(61170049,60903044)supported by National Natural Science Foundation of ChinaProject(2012AA010903)supported by National High Technology Research and Development Program of China
文摘Particle-in-cell (PIC) method has got much benefits from GPU-accelerated heterogeneous systems.However,the performance of PIC is constrained by the interpolation operations in the weighting process on GPU (graphic processing unit).Aiming at this problem,a fast weighting method for PIC simulation on GPU-accelerated systems was proposed to avoid the atomic memory operations during the weighting process.The method was implemented by taking advantage of GPU's thread synchronization mechanism and dividing the problem space properly.Moreover,software managed shared memory on the GPU was employed to buffer the intermediate data.The experimental results show that the method achieves speedups up to 3.5 times compared to previous works,and runs 20.08 times faster on one NVIDIA Tesla M2090 GPU compared to a single core of Intel Xeon X5670 CPU.
基金Projected supported by the National Natural Science Foundation of China(Grant Nos.11675025,11135007,and 11405012)
文摘In this paper, we analytically explore the magnetic field and mass density evolutions obtained in particle-in-cell（PIC）and magnetohydrodynamics（MHD） simulations of a rarefied deuterium shell Z-pinch and compare those results, and also we study the effects of artificially increased Spitzer resistivity on the magnetic field evolution and Z-pinch dynamic process in the MHD simulation. There are significant differences between the profiles of mass density in the PIC and MHD simulations before 45 ns of the Z-pinch in this study. However, after the shock formation in the PIC simulation,the mass density profile is similar to that in the MHD simulation in the case of using multiplier 2 to modify the Spitzer resistivity. Compared with the magnetic field profiles of the PIC simulation of the shell, the magnetic field diffusion has still not been sufficiently revealed in the MHD simulation even though their convergence ratios become the same by using larger multipliers in the resistivity. The MHD simulation results suggest that the magnetic field diffusion is greatly enhanced by increasing the Spitzer resistivity used, which, however, causes the implosion characteristic to change from shock compression to weak shock, even shockless evolution, and expedites the expansion of the shell. Too large a multiplier is not suggested to be used to modify the resistivity in some Z-pinch applications, such as the Z-pinch driven inertial confinement fusion（ICF） in a dynamic hohlraum. Two-fluid or Hall MHD model, even the PIC/fluid hybrid simulation would be considered as a suitable physical model when there exist the plasma regions with very low density in the simulated domain.
文摘A rainstorm caused by mesoscale convective system (MCS) in Guizhou Province in June 25-26 in 2005 was simulated with the MM5 model. Based on the good simulated results of the MCS developing and the clouds physics process, and by means of reducing the height of Yunnan-Guizhou Plateau and cutting off the middle-east of the Yunnan-Guizhou Plateau on the simulated tests, the question as how the ladder terrain on the west of Yunnan-Guizhou Plateau impact on the rainstorm of Guizhou was studied. The analysis results showed that the second ladder terrain of Yunnan-Guizhou Plateau only affected the development of convective clouds on its backward position,and hardly affected the rain on its upward. The whole terrain of the Yunnan-Guizhou Plateau had a distinct impact not only on the windward slope rainfall of the west of the plateau, but also on the rainfall distribution, intensity and continuing time of the convective clouds on the middle-east of the plateau.
