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 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.展开更多
A PIC （particle-in-cell）-MC （Monte Carlo） code to model electron beam transport into dense matter is developed. The background target is treated as a cold, stationary fluid and the fast electrons as particles with...A PIC （particle-in-cell）-MC （Monte Carlo） code to model electron beam transport into dense matter is developed. The background target is treated as a cold, stationary fluid and the fast electrons as particles with the relativistic motions. The process is described by a particle-in-cell method with consideration of the influence of both the self-generated electric and magnetic fields as well as collisions between the fast electrons and the target. The collisional part of the code is solved by the Monte Carlo-type method. Furthermore by assuming that the background current balances with the fast electron current, the electric field is given by the Ohm＇s law and the magnetic field is calculated from the Faraday＇s law. Both are solved in a two-dimensional cylindrical geometry. The algorithms implemented in the code are demonstrated and the numerical experiments are performed for monoenergy homogeneous fast electron beam transport in an aluminum target when the fields, collision and angular scattering are switched on and off independently.展开更多
We have developed a three dimensional （3D） PIC （particle-in-cell）-MC （Monte Carlo） code in order to simulate an electron beam transported into the dense matter based on our previous two dimensional code. The rel...We have developed a three dimensional （3D） PIC （particle-in-cell）-MC （Monte Carlo） code in order to simulate an electron beam transported into the dense matter based on our previous two dimensional code. The relativistic motion of fast electrons is treated by the particle-in-cell method under the influence of both a self-generated transverse magnetic field and an axial electric field, as well as collisions. The electric field generated by return current is expressed by Ohm＇s law and the magnetic field is calculated from Faraday＇s law. The slowing down of monoenergy electrons in DT plasma is calculated and discussed.展开更多
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.展开更多
By performing one-dimensional particle-in-cell simulations, the nonlinear effects of electronacoustic(EA) waves are investigated in a multispecies plasma, whose constituents are hot electrons, cold electrons, and beam...By performing one-dimensional particle-in-cell simulations, the nonlinear effects of electronacoustic(EA) waves are investigated in a multispecies plasma, whose constituents are hot electrons, cold electrons, and beam electrons with immobile neutralized positive ions. Numerical analyses have identified that EA waves with a sufficiently large amplitude tend to trap cold electrons. Because EA waves are dispersive, where the wave modes with different wavenumbers have different phase velocities, the trapping may lead to the mixing of cold electrons. The cold electrons finally get thermalized or heated. The investigation also shows that the excited EA waves give rise to a broad range of wave frequencies, which may be helpful for understanding the broadband-electrostatic-noise spectrum in the Earth’s auroral region.展开更多
Stimulated Raman particle-in-cell （PIC） simulations scattering （SRS） in a low-density The backward stimulated Raman plasma slab is investigated by scattering （B-SRS） dominates initially and erodes the head of th...Stimulated Raman particle-in-cell （PIC） simulations scattering （SRS） in a low-density The backward stimulated Raman plasma slab is investigated by scattering （B-SRS） dominates initially and erodes the head of the pump wave, while the forward stimulated Raman scattering （F-SRS） subsequently develops and is located at the rear part of the slab. Two-stage electron acceleration may be more efficient due to the coexistence of these two instabilities. The B-SRS plasma wave with low phase velocities can accelerate the background electrons which may be further boosted to higher energies by the F-SRS plasma wave with high phase velocities. The simulations show that the peaks of the main components in both the frequency and wave number spectra occur at the positions estimated from the phase-matching conditions.展开更多
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.展开更多
In this paper,a scheme of commonly-resonated extended interaction circuit system based on high order TMn,mode is proposed to lock the phases of two extended interaction oscillators(EIOs)for generating high power at G-...In this paper,a scheme of commonly-resonated extended interaction circuit system based on high order TMn,mode is proposed to lock the phases of two extended interaction oscillators(EIOs)for generating high power at G-band.Two separate EIOs are coupled through a specific single-gap coupling field supported by a designed gap waveguide with length Lg,which form the phase-locked EIOs based on the commonly-resonated system.As a whole system,the system has been focused on with mode analysis based on different single-gap coupling fields,mode hopping,which present the variation of phase difference between the two-beam-wave interactions when changing Lg.To demonstrate the effectiveness of the proposed circuit system in producing the phase locking,we conducted particle-in-cell(PIC)simulations to show that the interesting mode hopping occurs with the phase difference of O and r between the output signals from two output ports,corresponding to the excitation of the TMn mode with different n.Simulation results show that 1)the oscillator can deliver two times of the output power obtained from one single oscillator at 220 GHz,2)the two EIOs can still deliver output signals with phase difference of O and when the currents of the two beams are different or the fabrication errors of the two EIO cavities are taken into account.The proposed scheme is promising in extending to phase locking between multiple EIOs,and generating higher power at millimeter-wave and higher frequencies.展开更多
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.展开更多
An analytical model for describing the charged-particle transport in a wall-confined laser-induced decaying plasma is established under an external electrostatic field,focusing on the effects of the initial plasma-ele...An analytical model for describing the charged-particle transport in a wall-confined laser-induced decaying plasma is established under an external electrostatic field,focusing on the effects of the initial plasma-electrode gap(IPEG)that exists in applications such as laser isotope separation.This newly developed analytical model is validated by particle-in-cell simulations and the experimental scaling relation,and can also be reduced to its previously published counterpart that did not consider IPEGs.Based on this analytical model,the influences of different IPEG spacings on the characteristics of the whole ion extraction process are studied.The results show that the ion extraction ratios at the endpoints of the first and second stages both decrease with increasing IPEG spacing,while the corresponding time durations for the first two stages show a non-monotonous variation trend.The specific ion extraction time,defined as the ion extraction time per unit mass to comprehensively characterize the ion extraction efficiency,increases generally with the increase of IPEG spacing.This study not only provides further insight into the fundamental physical processes in a wall-bounded decaying plasma under an externally applied electrostatic field,but also offers useful theoretical guidance for optimal designs of geometrical and operating parameters in laser isotope separation processes.展开更多
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.展开更多
The irradiation of few-nm-thick targets by a finite-contrast high-intensity short-pulse laser results in a strong pre-expansion of these targets at the arrival time of the main pulse.The targets decompress to near and...The irradiation of few-nm-thick targets by a finite-contrast high-intensity short-pulse laser results in a strong pre-expansion of these targets at the arrival time of the main pulse.The targets decompress to near and lower than critical densities with plasmas extending over few micrometers,i.e.multiple wavelengths.The interaction of the main pulse with such a highly localized but inhomogeneous target leads to the generation of a short channel and further self-focusing of the laser beam.Experiments at the Glass Hybrid OPCPA Scaled Test-bed(GHOST)laser system at University of Texas,Austin using such targets measured non-Maxwellian,peaked electron distribution with large bunch charge and high electron density in the laser propagation direction.These results are reproduced in 2D PIC simulations using the EPOCH code,identifying direct laser acceleration(DLA)[1]as the responsible mechanism.This is the first time that DLA has been observed to produce peaked spectra as opposed to broad,Maxwellian spectra observed in earlier experiments[2].This high-density electrons have potential applications as injector beams for a further wakefield acceleration stage as well as for pump-probe applications.展开更多
The radiation reaction effects on electron dynamics in counter-propagating circularly polarized laser beams are investigated through the linearization theorem and the results are in great agreement with numeric soluti...The radiation reaction effects on electron dynamics in counter-propagating circularly polarized laser beams are investigated through the linearization theorem and the results are in great agreement with numeric solutions.For the first time,the properties of fixed points in electron phase-space were analyzed with linear stability theory,showing that center nodes will become attractors if the classical radiation reaction is considered.Electron dynamics are significantly affected by the properties of the fixed points and the electron phase-space densities are found to be increasing exponentially near the attractors.The density growth rates are derived theoretically and further verified by particle-in-cell sim-ulations,which can be detected in experiments to explore the effects of radiation reaction qualitatively.The attractor can also facilitate realizing a series of nanometer-scaled flying electron slices via adjusting the colliding laser frequencies.展开更多
With two-dimensional quantum electrodynamics（QED） particle-in-cell simulations, a dense electron-positron（e^-e~＋） pair generation from laser-solid interactions is demonstrated. When the interaction of two linearl...With two-dimensional quantum electrodynamics（QED） particle-in-cell simulations, a dense electron-positron（e^-e~＋） pair generation from laser-solid interactions is demonstrated. When the interaction of two linearly polarized laser pulses with a thin target enters into the relativistic transparency regime, a stable standing wave（SW） field can be formed by the overlap of the two counter-propagating laser pulses directly. The present study aims to clarify the effects of the S W field on the dynamics of e^-e^＋ pair plasmas. Our results indicate that under the combined effect of the SW field and radiation reaction（RR） effect, the created e^-e^＋ pairs can be trapped into the electric field nodes when the field strength is strong.The trapping effect contributes to the generation of γAV≥400 and ultra-dense pair plasmas in the two-side irradiation scheme. Despite different laser intensities, these pair plasmas have a Maxwellian spectral distribution with a peak energy of 150 MeV. Besides, the periodical modulation of the average energy, spatial, phase-space, and angular patterns of the e^-e^＋ pair plasmas can be triggered. In the angular patterns, as long as the SW field exists, pair plasmas can be pinched along the laser polarization direction. These results may offer a better understanding of the laser-solid interactions in the experiments when 10-PW laser facilities come into operation in the future.展开更多
It is known that ion channel can effectively limit the radial expansion of an artificial electron beam during its longrange propagation in the space plasma environment.Most prior studies discussed the focusing charact...It is known that ion channel can effectively limit the radial expansion of an artificial electron beam during its longrange propagation in the space plasma environment.Most prior studies discussed the focusing characteristics of the beam in the ion channel,but the establishment process and transient properties of the ion channel itself,which also plays a crucial role during the propagation of the relativistic electron beam in the plasma environment,were commonly neglected.In this study,a series of two-dimensional(2D)particle-in-cell simulations is performed and an analytical model of ion channel oscillation is constructed according to the single-particle motion.The results showed that when the beam density is higher than the density of plasma environment,ion channel can be established and always continues to oscillate periodically over the entire propagation.Multiple factors,including the beam electron density,initial beam radius,and the plasma density can affect the oscillation properties of ion channel.Axial velocity of the beam oscillates synchronously with the ion channel and this phenomenon will finally develop into a two-stream instability which can seriously affect the effective transport for relativistic electron beam.Choosing appropriate beam parameters based on various plasma environments may contribute to the improvement of the stability of ion channel.Additionally,radial expansion of the beam can be limited by ion channel and a stable long-range propagation in terrestrial atmosphere may be achieved.展开更多
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.展开更多
Space satellite observations in an electron phase-space hole(electron hole) have shown that bipolar structures are discovered at the parallel cut of parallel electric field, while unipolar structures spring from the p...Space satellite observations in an electron phase-space hole(electron hole) have shown that bipolar structures are discovered at the parallel cut of parallel electric field, while unipolar structures spring from the parallel cut of perpendicular electric field. Particle-in-cell(PIC) simulations have demonstrated that the electron bi-stream instability induces several electron holes during its nonlinear evolution. However, how the unipolar structure of the parallel cut of the perpendicular electric field formed in these electron holes is still an unsolved problem,especially in a strongly magnetized plasma(Ω_e >ω_(pe), where Ω_e is defined as electron gyrofrequency and ω_(pe) is defined as plasma frequency, respectively). In this paper, with two-dimensional(2D) electrostatic PIC simulations, the evolution of the electron two-stream instability with a finite width in strongly magnetized plasma is investigated. Initially, those conditions lead to monochromatic electrostatic waves, and these waves coalesce with each other during their nonlinear evolution. At last, a solitary electrostatic structure is formed. In such an electron hole, a bipolar structure is formed in the parallel cut. of parallel electric field, while a unipolar structure presents in the parallel cut of perpendicular electric field.展开更多
It is known that gas flow rate is a key factor in controlling industrial plasma processing. In this paper, a 2D PIC/MCC model is developed for an rf hollow cathode discharge with an axial nitrogen gas flow. The effect...It is known that gas flow rate is a key factor in controlling industrial plasma processing. In this paper, a 2D PIC/MCC model is developed for an rf hollow cathode discharge with an axial nitrogen gas flow. The effects of the gas flow rate on the plasma parameters are calculated and the results show that： with an increasing flow rate, the total ion（N＋2, N＋） density decreases, the mean sheath thickness becomes wider, the radial electric field in the sheath and the axial electric field show an increase, and the energies of both kinds of nitrogen ions increase;and, as the axial ion current density that is moving toward the ground electrode increases, the ion current density near the ground electrode increases. The simulation results will provide a useful reference for plasma jet technology involving rf hollow cathode discharges in N2.展开更多
Hollow cathodes serve as electron sources in Hall thrusters,ion thrusters and other electric propulsion systems.One of the vital problems in their application is the cathode erosion.However,the basic erosion mechanism...Hollow cathodes serve as electron sources in Hall thrusters,ion thrusters and other electric propulsion systems.One of the vital problems in their application is the cathode erosion.However,the basic erosion mechanism and the source of high-energy ions cause of erosion are not fully understood.In this paper,both potential measurements and simulation analyses were performed to explain the formation of high-energy ions.A high-speed camera,a single Langmuir probe and a floating emissive probe were used to determine the steady and oscillatory plasma properties in the near-field plume of a hollow cathode.The temporal structure,electron temperature,electron density,and both static and oscillation of plasma potentials of the plume have been obtained by the diagnostics mentioned above.The experimental results show that there exists a potential hill（about 30 V） and also severe potential oscillations in the near-plume region.Moreover,a simple 2 D particle-in-cell model was used to analyze the energy transition between the potential hill and/or its oscillations and the ions.The simulation results show that the energy of ions gained from the static potential background is about 20 e V,but it could reach to 60 e V when the plasma oscillates.展开更多
基金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 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 by the National High Technology ICF Committee of Chinathe National Natural Science Fund of China(Nos. 10335020, 10375011, 10576007)the Laboratory of Computational Physics (No. 51479050205ZW0905)
文摘A PIC （particle-in-cell）-MC （Monte Carlo） code to model electron beam transport into dense matter is developed. The background target is treated as a cold, stationary fluid and the fast electrons as particles with the relativistic motions. The process is described by a particle-in-cell method with consideration of the influence of both the self-generated electric and magnetic fields as well as collisions between the fast electrons and the target. The collisional part of the code is solved by the Monte Carlo-type method. Furthermore by assuming that the background current balances with the fast electron current, the electric field is given by the Ohm＇s law and the magnetic field is calculated from the Faraday＇s law. Both are solved in a two-dimensional cylindrical geometry. The algorithms implemented in the code are demonstrated and the numerical experiments are performed for monoenergy homogeneous fast electron beam transport in an aluminum target when the fields, collision and angular scattering are switched on and off independently.
基金National High Technology ICF Committee in ChinaNational Natural Science Fund of China(Nos.10675024,10335020,10375011,and 10576007)+1 种基金National Basic Research Program of China(973 Program)(No.2007CB815101)the Laboratory of Computational Physics(No.51479050205ZW0905)
文摘We have developed a three dimensional （3D） PIC （particle-in-cell）-MC （Monte Carlo） code in order to simulate an electron beam transported into the dense matter based on our previous two dimensional code. The relativistic motion of fast electrons is treated by the particle-in-cell method under the influence of both a self-generated transverse magnetic field and an axial electric field, as well as collisions. The electric field generated by return current is expressed by Ohm＇s law and the magnetic field is calculated from Faraday＇s law. The slowing down of monoenergy electrons in DT plasma is calculated and discussed.
基金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 support from Chinese Academy of Science(CAS)TWAS for his Ph.D studies at the University of Science and Technology of China in the category of a 2016 CAS-TWAS President’s Fellowship Awardee(Series No.2016-172)+1 种基金partially supported by National Natural Science Foundation of China(Nos.41331067,41774169,and 41527804)the Key Research Program of Frontier Sciences,CAS(QYZDJ-SSW-DQC010)
文摘By performing one-dimensional particle-in-cell simulations, the nonlinear effects of electronacoustic(EA) waves are investigated in a multispecies plasma, whose constituents are hot electrons, cold electrons, and beam electrons with immobile neutralized positive ions. Numerical analyses have identified that EA waves with a sufficiently large amplitude tend to trap cold electrons. Because EA waves are dispersive, where the wave modes with different wavenumbers have different phase velocities, the trapping may lead to the mixing of cold electrons. The cold electrons finally get thermalized or heated. The investigation also shows that the excited EA waves give rise to a broad range of wave frequencies, which may be helpful for understanding the broadband-electrostatic-noise spectrum in the Earth’s auroral region.
基金supported by National High Technology ICF Committee in Chinathe National Natural Science Fund of China(Nos.10675024,10335020,10375011 and 10576007)the Laboratory of Computational Physics(No.51479050205ZW0905)
文摘Stimulated Raman particle-in-cell （PIC） simulations scattering （SRS） in a low-density The backward stimulated Raman plasma slab is investigated by scattering （B-SRS） dominates initially and erodes the head of the pump wave, while the forward stimulated Raman scattering （F-SRS） subsequently develops and is located at the rear part of the slab. Two-stage electron acceleration may be more efficient due to the coexistence of these two instabilities. The B-SRS plasma wave with low phase velocities can accelerate the background electrons which may be further boosted to higher energies by the F-SRS plasma wave with high phase velocities. The simulations show that the peaks of the main components in both the frequency and wave number spectra occur at the positions estimated from the phase-matching conditions.
基金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.
基金Supported in part by the National Natural Science Foundation of China(62401125)the Natural Science Foundation of Sichuan Province(2023NSFSC1376)the Fundamental Research Funds for the Central Universities(ZYGX2024J008)。
文摘In this paper,a scheme of commonly-resonated extended interaction circuit system based on high order TMn,mode is proposed to lock the phases of two extended interaction oscillators(EIOs)for generating high power at G-band.Two separate EIOs are coupled through a specific single-gap coupling field supported by a designed gap waveguide with length Lg,which form the phase-locked EIOs based on the commonly-resonated system.As a whole system,the system has been focused on with mode analysis based on different single-gap coupling fields,mode hopping,which present the variation of phase difference between the two-beam-wave interactions when changing Lg.To demonstrate the effectiveness of the proposed circuit system in producing the phase locking,we conducted particle-in-cell(PIC)simulations to show that the interesting mode hopping occurs with the phase difference of O and r between the output signals from two output ports,corresponding to the excitation of the TMn mode with different n.Simulation results show that 1)the oscillator can deliver two times of the output power obtained from one single oscillator at 220 GHz,2)the two EIOs can still deliver output signals with phase difference of O and when the currents of the two beams are different or the fabrication errors of the two EIO cavities are taken into account.The proposed scheme is promising in extending to phase locking between multiple EIOs,and generating higher power at millimeter-wave and higher frequencies.
基金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.
基金supported by the National Key Laboratory of Particle Transport and Separation Technology(Grant No.WZKF-2024-2).
文摘An analytical model for describing the charged-particle transport in a wall-confined laser-induced decaying plasma is established under an external electrostatic field,focusing on the effects of the initial plasma-electrode gap(IPEG)that exists in applications such as laser isotope separation.This newly developed analytical model is validated by particle-in-cell simulations and the experimental scaling relation,and can also be reduced to its previously published counterpart that did not consider IPEGs.Based on this analytical model,the influences of different IPEG spacings on the characteristics of the whole ion extraction process are studied.The results show that the ion extraction ratios at the endpoints of the first and second stages both decrease with increasing IPEG spacing,while the corresponding time durations for the first two stages show a non-monotonous variation trend.The specific ion extraction time,defined as the ion extraction time per unit mass to comprehensively characterize the ion extraction efficiency,increases generally with the increase of IPEG spacing.This study not only provides further insight into the fundamental physical processes in a wall-bounded decaying plasma under an externally applied electrostatic field,but also offers useful theoretical guidance for optimal designs of geometrical and operating parameters in laser isotope separation processes.
基金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 NNSA cooperative agreement DE-NA0002008the Defense Advanced Research Projects Agency's PULSE program(12-63-PULSE-FP014)the Air Force Office of Scientific Research(FA9550-14-1-0045).
文摘The irradiation of few-nm-thick targets by a finite-contrast high-intensity short-pulse laser results in a strong pre-expansion of these targets at the arrival time of the main pulse.The targets decompress to near and lower than critical densities with plasmas extending over few micrometers,i.e.multiple wavelengths.The interaction of the main pulse with such a highly localized but inhomogeneous target leads to the generation of a short channel and further self-focusing of the laser beam.Experiments at the Glass Hybrid OPCPA Scaled Test-bed(GHOST)laser system at University of Texas,Austin using such targets measured non-Maxwellian,peaked electron distribution with large bunch charge and high electron density in the laser propagation direction.These results are reproduced in 2D PIC simulations using the EPOCH code,identifying direct laser acceleration(DLA)[1]as the responsible mechanism.This is the first time that DLA has been observed to produce peaked spectra as opposed to broad,Maxwellian spectra observed in earlier experiments[2].This high-density electrons have potential applications as injector beams for a further wakefield acceleration stage as well as for pump-probe applications.
基金The work has been supported by the National Basic Research Program of China(Grant No.2013CBA01502),NSFC(Grant No.11535001)National Grand Instrument Project(2012YQ030142).
文摘The radiation reaction effects on electron dynamics in counter-propagating circularly polarized laser beams are investigated through the linearization theorem and the results are in great agreement with numeric solutions.For the first time,the properties of fixed points in electron phase-space were analyzed with linear stability theory,showing that center nodes will become attractors if the classical radiation reaction is considered.Electron dynamics are significantly affected by the properties of the fixed points and the electron phase-space densities are found to be increasing exponentially near the attractors.The density growth rates are derived theoretically and further verified by particle-in-cell sim-ulations,which can be detected in experiments to explore the effects of radiation reaction qualitatively.The attractor can also facilitate realizing a series of nanometer-scaled flying electron slices via adjusting the colliding laser frequencies.
基金Project supported by the National Basic Research Program of China(Grant No.2013CBA01504)the National Natural Science Foundation of China(Grant Nos.11347028,11405083,and 11675075)+1 种基金the Natural Science Foundation of Hunan Province,China(Grant No.2018JJ2315)the Youth Talent Project of Hunan Province,China(Grant No.2018RS3096)
文摘With two-dimensional quantum electrodynamics（QED） particle-in-cell simulations, a dense electron-positron（e^-e~＋） pair generation from laser-solid interactions is demonstrated. When the interaction of two linearly polarized laser pulses with a thin target enters into the relativistic transparency regime, a stable standing wave（SW） field can be formed by the overlap of the two counter-propagating laser pulses directly. The present study aims to clarify the effects of the S W field on the dynamics of e^-e^＋ pair plasmas. Our results indicate that under the combined effect of the SW field and radiation reaction（RR） effect, the created e^-e^＋ pairs can be trapped into the electric field nodes when the field strength is strong.The trapping effect contributes to the generation of γAV≥400 and ultra-dense pair plasmas in the two-side irradiation scheme. Despite different laser intensities, these pair plasmas have a Maxwellian spectral distribution with a peak energy of 150 MeV. Besides, the periodical modulation of the average energy, spatial, phase-space, and angular patterns of the e^-e^＋ pair plasmas can be triggered. In the angular patterns, as long as the SW field exists, pair plasmas can be pinched along the laser polarization direction. These results may offer a better understanding of the laser-solid interactions in the experiments when 10-PW laser facilities come into operation in the future.
基金supported by the Joint Funds of the National Natural Science Foundation of China(Grant Nos.61372050 and U1730247).
文摘It is known that ion channel can effectively limit the radial expansion of an artificial electron beam during its longrange propagation in the space plasma environment.Most prior studies discussed the focusing characteristics of the beam in the ion channel,but the establishment process and transient properties of the ion channel itself,which also plays a crucial role during the propagation of the relativistic electron beam in the plasma environment,were commonly neglected.In this study,a series of two-dimensional(2D)particle-in-cell simulations is performed and an analytical model of ion channel oscillation is constructed according to the single-particle motion.The results showed that when the beam density is higher than the density of plasma environment,ion channel can be established and always continues to oscillate periodically over the entire propagation.Multiple factors,including the beam electron density,initial beam radius,and the plasma density can affect the oscillation properties of ion channel.Axial velocity of the beam oscillates synchronously with the ion channel and this phenomenon will finally develop into a two-stream instability which can seriously affect the effective transport for relativistic electron beam.Choosing appropriate beam parameters based on various plasma environments may contribute to the improvement of the stability of ion channel.Additionally,radial expansion of the beam can be limited by ion channel and a stable long-range propagation in terrestrial atmosphere may be achieved.
基金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.
基金Supported by the National Science Foundation of China(41474125,41331067,41421063)973 Program(2013CBA01503)Key Research Program of Frontier Sciences,CAS(QYZDJ-SSW-DQC010)
文摘Space satellite observations in an electron phase-space hole(electron hole) have shown that bipolar structures are discovered at the parallel cut of parallel electric field, while unipolar structures spring from the parallel cut of perpendicular electric field. Particle-in-cell(PIC) simulations have demonstrated that the electron bi-stream instability induces several electron holes during its nonlinear evolution. However, how the unipolar structure of the parallel cut of the perpendicular electric field formed in these electron holes is still an unsolved problem,especially in a strongly magnetized plasma(Ω_e >ω_(pe), where Ω_e is defined as electron gyrofrequency and ω_(pe) is defined as plasma frequency, respectively). In this paper, with two-dimensional(2D) electrostatic PIC simulations, the evolution of the electron two-stream instability with a finite width in strongly magnetized plasma is investigated. Initially, those conditions lead to monochromatic electrostatic waves, and these waves coalesce with each other during their nonlinear evolution. At last, a solitary electrostatic structure is formed. In such an electron hole, a bipolar structure is formed in the parallel cut. of parallel electric field, while a unipolar structure presents in the parallel cut of perpendicular electric field.
基金supported by the Natural Science Foundation of Hebei Province,China(No.A2012205072)
文摘It is known that gas flow rate is a key factor in controlling industrial plasma processing. In this paper, a 2D PIC/MCC model is developed for an rf hollow cathode discharge with an axial nitrogen gas flow. The effects of the gas flow rate on the plasma parameters are calculated and the results show that： with an increasing flow rate, the total ion（N＋2, N＋） density decreases, the mean sheath thickness becomes wider, the radial electric field in the sheath and the axial electric field show an increase, and the energies of both kinds of nitrogen ions increase;and, as the axial ion current density that is moving toward the ground electrode increases, the ion current density near the ground electrode increases. The simulation results will provide a useful reference for plasma jet technology involving rf hollow cathode discharges in N2.
基金financial support from National Natural Science Foundation of China under Grant Nos.11402025 and 11475019China Academy of Space Technology under Grant Nos.YJJ0701 and ZWK1608
文摘Hollow cathodes serve as electron sources in Hall thrusters,ion thrusters and other electric propulsion systems.One of the vital problems in their application is the cathode erosion.However,the basic erosion mechanism and the source of high-energy ions cause of erosion are not fully understood.In this paper,both potential measurements and simulation analyses were performed to explain the formation of high-energy ions.A high-speed camera,a single Langmuir probe and a floating emissive probe were used to determine the steady and oscillatory plasma properties in the near-field plume of a hollow cathode.The temporal structure,electron temperature,electron density,and both static and oscillation of plasma potentials of the plume have been obtained by the diagnostics mentioned above.The experimental results show that there exists a potential hill（about 30 V） and also severe potential oscillations in the near-plume region.Moreover,a simple 2 D particle-in-cell model was used to analyze the energy transition between the potential hill and/or its oscillations and the ions.The simulation results show that the energy of ions gained from the static potential background is about 20 e V,but it could reach to 60 e V when the plasma oscillates.
