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.展开更多
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.展开更多
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.展开更多
This paper presents a numerical simulation of DC glow discharge at 2 Torr(1 Torr=1.33322×10~2 Pa)with a microhollow anode,using a two-dimensional(2D)PIC/MCC code to explore the impact of the hollow anode structur...This paper presents a numerical simulation of DC glow discharge at 2 Torr(1 Torr=1.33322×10~2 Pa)with a microhollow anode,using a two-dimensional(2D)PIC/MCC code to explore the impact of the hollow anode structure on discharge characteristics.Simulation results show that the plasma density in the anode channel decreases exponentially along the x-direction towards the outlet while the electric potential decreases linearly.The electron temperature,derived from the relationship between density and electric potential,shows a good agreement with the calculated temperature.The potential peak at the inlet,along with the grounded anode,forms a transverse potential well that draws electrons toward the center of the channel.The x-direction electric field,generated by the potential difference between the inlet and outlet of the anode,directs electrons toward the inlet of the anode.Low-energy electrons are confined within the potential well,while very few high-energy electrons outlet the channel.The hollow anode structure serves as a collimator for electrons,leading to an increment of moderate-velocity electron flux in the anode channel.Shortening the anode leads to a slower rate of density decrease,resulting in a higher density at the outlet.展开更多
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.展开更多
The deposition layer on the discharge channel wall of a Hall thruster after long-term operation occasionally detaches from the wall and interferes with the plasma inside the channel,resulting in current pulse and thre...The deposition layer on the discharge channel wall of a Hall thruster after long-term operation occasionally detaches from the wall and interferes with the plasma inside the channel,resulting in current pulse and threatening the power supply.To understand the generation mechanism of the current pulse,the interaction between the peeling material and the channel plasma was studied by the particle-in-cell and Monte Carlo collision(PIC/MCC)method.The plasma evolution and distribution in the Hall thruster channel when peeling material existed were simulated.Simulation results show that the peeling material changes the electron transportation and ionization,and the effects are related to the position of the peeling material.When the peeling material invades the zone where the ionization rate is originally the highest without peeling material,the ionization process is almost terminated.The ionization suppression by the peeling material will cause the propellant gas to refill the channel exit.As the peeling material moves away,gas discharge is induced near the channel exit,resulting in anode current pulse.展开更多
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.展开更多
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.展开更多
基金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.
文摘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.
文摘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.
基金Project partially supported by the National Natural Science Foundation of China(Grant Nos.12275060 and12075223)the Fund from Frontier Science Center for Interaction between Space Environment and Matter(Grant No.5740401024)。
文摘This paper presents a numerical simulation of DC glow discharge at 2 Torr(1 Torr=1.33322×10~2 Pa)with a microhollow anode,using a two-dimensional(2D)PIC/MCC code to explore the impact of the hollow anode structure on discharge characteristics.Simulation results show that the plasma density in the anode channel decreases exponentially along the x-direction towards the outlet while the electric potential decreases linearly.The electron temperature,derived from the relationship between density and electric potential,shows a good agreement with the calculated temperature.The potential peak at the inlet,along with the grounded anode,forms a transverse potential well that draws electrons toward the center of the channel.The x-direction electric field,generated by the potential difference between the inlet and outlet of the anode,directs electrons toward the inlet of the anode.Low-energy electrons are confined within the potential well,while very few high-energy electrons outlet the channel.The hollow anode structure serves as a collimator for electrons,leading to an increment of moderate-velocity electron flux in the anode channel.Shortening the anode leads to a slower rate of density decrease,resulting in a higher density at the outlet.
基金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.
基金supported by National Natural Science Foundation of China(No.U23B20152).
文摘The deposition layer on the discharge channel wall of a Hall thruster after long-term operation occasionally detaches from the wall and interferes with the plasma inside the channel,resulting in current pulse and threatening the power supply.To understand the generation mechanism of the current pulse,the interaction between the peeling material and the channel plasma was studied by the particle-in-cell and Monte Carlo collision(PIC/MCC)method.The plasma evolution and distribution in the Hall thruster channel when peeling material existed were simulated.Simulation results show that the peeling material changes the electron transportation and ionization,and the effects are related to the position of the peeling material.When the peeling material invades the zone where the ionization rate is originally the highest without peeling material,the ionization process is almost terminated.The ionization suppression by the peeling material will cause the propellant gas to refill the channel exit.As the peeling material moves away,gas discharge is induced near the channel exit,resulting in anode current pulse.
基金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.
文摘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.
