We have proposed a general numerical framework for plasma simulations on graphics processing unit clusters based on microscopic kinetic equations with full collision terms.Our numerical algorithm consistently deals wi...We have proposed a general numerical framework for plasma simulations on graphics processing unit clusters based on microscopic kinetic equations with full collision terms.Our numerical algorithm consistently deals with both long-range(classical forces in the Vlasov term)and short-range(quantum processes in the collision term)interactions.Providing the relevant particle masses,charges and types(classical,fermionic or bosonic),as well as the external forces and the matrix elements(in the collisional integral),the algorithm consistently solves the coupled multi-particle kinetic equations.Currently,the framework is being tested and applied in the field of relativistic heavy-ion collisions;extensions to other plasma systems are straightforward.Our framework is a potential and competitive numerical platform for consistent plasma simulations.展开更多
The paper builds the high-current plasma beams model under different dimensions (1D, 2D, and 3D) by continuum (magnetohydrodynamics MHD) and statistical (Monte Carlo MC) mechanics under conditions of low pressures (10...The paper builds the high-current plasma beams model under different dimensions (1D, 2D, and 3D) by continuum (magnetohydrodynamics MHD) and statistical (Monte Carlo MC) mechanics under conditions of low pressures (10<sup>-3</sup> Pa). After detailed presentation of the model, two methods firstly have been analyzed in terms of plasma beam properties. Then, we compare the simulation results of MHD numerical simulation with MC stochastic particles simulation. Finally, through further analysis, it is demonstrated that integrated hybrid MHD and MC method (IMHDMC) provides an innovative practical tool to capture essential properties of high-current plasma beams.展开更多
A high-energy electron beam generator is used to generate a plasma in atmosphere. Based on a Monte Carlo toolkit named GEANT4, a model including complete physics processes is established to simulate the passage of the...A high-energy electron beam generator is used to generate a plasma in atmosphere. Based on a Monte Carlo toolkit named GEANT4, a model including complete physics processes is established to simulate the passage of the electron beam in air. Based on the model, the characteristics of the electron beam air plasma are calculated. The energy distribution of beam electrons （BEs） indicates that high-energy electrons almost reside in the centre region of the beam, but low-energy electrons always live in the fringe area. The energy deposition is calculated in two cases, i.e., with and without secondary electrons （SEs）. Analysis indicates that the energy deposition of SEs accounts for a large part of the total energy deposition. The results of the energy spectrum show that the electrons in the inlet layer of the low-pressure chamber （LPC） are monoenergetic, but the energy spectrum of the electrons in the outlet layer is not pure. The SEs are largely generated at the outlet of the LPC. Moreover, both the energy distribution of BEs and the magnitude of the density of SEs are closely related to the pressure of LPC. Thus, a conclusion is drawn that a low magnitude of LPC pressure is helpful for reducing the energy loss in the LPC and also useful for greatly increasing the secondary electron density in dense air.展开更多
A volume-weighting cloud-in-cell (VW-CIC) model is developed to implement the particle-in-cell (PIC) simulation in axially symmetric systems. This model gives a first-order accuracy in the cylindrical system, and it i...A volume-weighting cloud-in-cell (VW-CIC) model is developed to implement the particle-in-cell (PIC) simulation in axially symmetric systems. This model gives a first-order accuracy in the cylindrical system, and it is incorporated into a PIC code. A planar diode with a finite-radius circular emitter is simulated with the code. The simulation results show that the VW-CIC model has a better accuracy and a lower noise than the conventional area-weighting cloud-in-cell (AW-CIC) model, especially on those points near the axis. The two-dimensional (2-D) space-charge-limited current density obtained from VW-CIC model is in better agreement with Lau’s analytical result. This model is more suitable for 2.5-D PIC simulation of axially symmetric plasmas.展开更多
For exploiting advantages of electron beam air plasma in some unusual applications, a Monte Carlo （MC） model coupled with heat transfer model is established to simulate the characteristics of electron beam air plasm...For exploiting advantages of electron beam air plasma in some unusual applications, a Monte Carlo （MC） model coupled with heat transfer model is established to simulate the characteristics of electron beam air plasma by considering the self-heating effect. Based on the model, the electron beam induced temperature field and the related plasma properties are investigated. The results indicate that a nonuniform temperature field is formed in the electron beam plasma region and the average temperature is of the order of 600 K. Moreover, much larger volume pear-shaped electron beam plasma is produced in hot state rather than in cold state. The beam ranges can, with beam energies of 75 keV and 80 keV, exceed 1.0 m and 1.2 m in air at pressure of 100 torr, respectively. Finally, a well verified formula is obtained for calculating the range of high energy electron beam in atmosphere.展开更多
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 order to study the extraction and acceleration mechanism of the dual-stage grid,a three-dimensional model based on the Particle-In-Cell/Monte Carlo Collision(PIC/MCC)method is performed.Dual-stage grid ion thruster...In order to study the extraction and acceleration mechanism of the dual-stage grid,a three-dimensional model based on the Particle-In-Cell/Monte Carlo Collision(PIC/MCC)method is performed.Dual-stage grid ion thruster is a new type of electrostatic ion thruster,which can break through the limitations of traditional gridded ion thrusters,and greatly improve the specific impulse.The high performance also makes the grid sensitive to operating parameters.In this paper,the influence of grid parameters on xenon ion thruster’s performance in a wide range is systematically simulated,and the optimal operating condition is given.Both the over-focusing of the plume,and the transparency of the screen grid are improved,and the grid corrosion is reduced through simulation optimization.The specific impulse under the given working conditions is 9877.24 s and the thrust is 7.28 mN.Based on the simulation optimization,the limitation of the dual-stage grid is discussed.The grid performs well under high voltage conditions(>3000 V)but not well under low voltage conditions(<2000 V).Finally,since argon is cheaper and more advantageous in future engineering applications,the plasma distribution and grid extraction ability under xenon and argon are analyzed and compared to study the flexibility of the dual-stage grid ion thruster.The simulation results show that a set of optimal parameters is only applicable to the corresponding propellant,which needs to be optimized for different propellant types.展开更多
Chirped wideband pump and seed pulses are usually considered for backward Raman amplification（BRA） in plasmas to achieve an extremely high-power laser pulse. However, current theoretical models only contain either a...Chirped wideband pump and seed pulses are usually considered for backward Raman amplification（BRA） in plasmas to achieve an extremely high-power laser pulse. However, current theoretical models only contain either a chirped pump or a chirped seed. In this paper, modified three-wave coupling equations are proposed for the BRA in the plasmas with both chirped wideband pump and seed. The simulation results can more precisely describe the experiments, such as the Princeton University experiment. The optimized chirp and bandwidth are determined based on the simulation to enhance the output intensity and efficiency.展开更多
In atmospheric radio-frequency （rf） discharges, the plasma parameters, such as electron density, sheath thickness and sheath voltage, are not easy to be probed experimentally, while the electrical characteristics, s...In atmospheric radio-frequency （rf） discharges, the plasma parameters, such as electron density, sheath thickness and sheath voltage, are not easy to be probed experimentally, while the electrical characteristics, such as impedance, resistance and reactance, are relatively convenient to be measured. In this paper we presented a simple theoretical model derived from the fluid description of generated plasmas without considering the circuit model, to investigate the relationship between the plasma impedance and plasma parameters. By introducing a relaxation frequency, the plasma impedance could be predicted by formulas presented in this study, and the mean electron density and sheath thickness can also be calculated from the measured or simulated impedance and reactance, respectively.展开更多
A numerical study of the effect of water content on OH production in a pulsed-dc atmospheric pressure helium-air plasma jet is presented. The generation and loss mechanisms of the OH radicals in a positive half-cycle ...A numerical study of the effect of water content on OH production in a pulsed-dc atmospheric pressure helium-air plasma jet is presented. The generation and loss mechanisms of the OH radicals in a positive half-cycle of the applied voltage are studied and discussed. It is found that the peak OH density increases with water content in air （varying from 0% to 1%） and reaches 6.3 x 10^18 m-3 when the water content is 1%. Besides, as the water content increases from 0.01% to 1%, the space-averaged reaction rate of three-body recombination increases dramatically and is comparable to those of main OH generation reactions.展开更多
Atomic and molecular processes relevant to the volumetric recombinationphenomena were investigated in a linear divertor plasma simulator MAP-Ⅱ. Volumetric recombinationis induced in He plasma by puffing of He or H_2....Atomic and molecular processes relevant to the volumetric recombinationphenomena were investigated in a linear divertor plasma simulator MAP-Ⅱ. Volumetric recombinationis induced in He plasma by puffing of He or H_2. In the He puffing case, the reduction of the ionflux is dominated by the electron-ion recombination. In the H_2 puffing case, however, it isdominated by the molecule-assisted recombination (MAR), which is characterized by the disappearanceof the Helium Rydberg spectra and by the existence of the hydrogen negative ions. Currentachievement and the future prospect are described.展开更多
Along with the introduction of the concept of dual-channel communication,we utilized the finite-difference time-domain（FDTD） method to simulate and measure the radiation pattern under certain plasma densities and pl...Along with the introduction of the concept of dual-channel communication,we utilized the finite-difference time-domain（FDTD） method to simulate and measure the radiation pattern under certain plasma densities and plasma collision frequencies.Results show that under certain settings,the radiation pattern of a plasma antenna resembles that of a metallic antenna.In contrast to a metallic antenna,a plasma antenna possesses other functionalities,such as dynamic reconfiguration and digital controllability.The data from simulation are similar to the measurement results,indicating that column plasma antenna can realize dual-channel communication.This work confirms the viability of realizing dual-channel communication by column plasma antenna,which adds a new but promising method for modern intelligent communication.展开更多
Plasma disruption is often an unavoidable aspect of tokamak operations. It may cause severe damage to in-vessel components such as the vacuum vessel conductors, the first wall and the divertor target plates. Two types...Plasma disruption is often an unavoidable aspect of tokamak operations. It may cause severe damage to in-vessel components such as the vacuum vessel conductors, the first wall and the divertor target plates. Two types of disruption, the hot-plasma vertical displacement event and the major disruption with a cold-plasma vertical displacement event, are simulated by the DINA code for HL-2M. The time evolutions of the plasma current, the halo current, the magnetic axis, the minor radius, the elongation as well as the electromagnetic force and eddy currents on the vacuum vessel during the thermal quench and the current quench are investigated. By comparing the electromagnetic forces before and after the disruption, we find that the disruption causes great damage to the vacuum vessel conductors. In addition, the hot-plasma vertical displacement event is more dangerous than the major disruption with the cold-plasma vertical displacement event.展开更多
The capacitively coupled radio frequency（CCRF）plasma has been widely used in various fields.In some cases,it requires us to estimate the range of key plasma parameters simpler and quicker in order to understand the ...The capacitively coupled radio frequency（CCRF）plasma has been widely used in various fields.In some cases,it requires us to estimate the range of key plasma parameters simpler and quicker in order to understand the behavior in plasma.In this paper,a glass vacuum chamber and a pair of plate electrodes were designed and fabricated,using 13.56 MHz radio frequency（RF）discharge technology to ionize the working gas of Ar.This discharge was mathematically described with equivalent circuit model.The discharge voltage and current of the plasma were measured atdifferent pressures and different powers.Based on the capacitively coupled homogeneous discharge model,the equivalent circuit and the analytical formula were established.The plasma density and temperature were calculated by using the equivalent impedance principle and energy balance equation.The experimental results show that when RF discharge power is 50–300 W and pressure is 25–250 Pa,the average electron temperature is about 1.7–2.1 e V and the average electron density is about 0.5？×10^17–3.6？×10^17m^-3.Agreement was found when the results were compared to those given by optical emission spectroscopy and COMSOL simulation.展开更多
Explicit structure-preserving geometric particle-in-cell(PIC)algorithm in curvilinear orthogonal coordinate systems is developed.The work reported represents a further development of the structure-preserving geometric...Explicit structure-preserving geometric particle-in-cell(PIC)algorithm in curvilinear orthogonal coordinate systems is developed.The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research.The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms,discrete exterior calculus,and explicit non-canonical symplectic integration.In addition to the truncated infinitely dimensional symplectic structure,the algorithm preserves exactly many important physical symmetries and conservation laws,such as local energy conservation,gauge symmetry and the corresponding local charge conservation.As a result,the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics.The algorithm has been implemented in the Sym PIC code,which is designed for highefficiency massively-parallel PIC simulations in modern clusters.The code has been applied to carry out whole-device 6 D kinetic simulation studies of tokamak physics.A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor.The state also admits a steady-state subsonic ion flow in the range of 10 km s-1,agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated.It is shown that high-n ballooning modes have larger growth rates than low-n global modes,and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2%level by the E×B flow generated by the instability.These results are consistent with early and recent electromagnetic gyrokinetic simulations.展开更多
The structural characteristics of zonal flows and their roles in the nonlinear interaction of multi-scale multi-mode turbulence are investigated numerically via a self-consistent Landau-fluid model.The multi-mode turb...The structural characteristics of zonal flows and their roles in the nonlinear interaction of multi-scale multi-mode turbulence are investigated numerically via a self-consistent Landau-fluid model.The multi-mode turbulence here is composed of a shorter wavelength electromagnetic(EM)ion temperature gradient(ITG)mode and a Kelvin-Helmholtz(KH)instability with long wavelengths excited by externally imposed small-scale shear flows.For strong shear flow,a prominent periodic intermittency of fluctuation intensity except for dominant ITG component is revealed in turbulence evolution,which onset time depends on the ion temperature gradient and the shear flow amplitudes corresponding to different KH instabilities.It is identified that the intermittency phenomenon results from the zonal flow dynamics,which is mainly generated by the KH mode and back-reacts on it.It is demonstrated that the odd symmetric components of zonal flow(same symmetry as the external flow)make the radial parity of the KH mode alteration through adjusting the drift velocities at two sides of the resonant surface so that the KH mode becomes bursty first.Afterwards,the ITG intermittency follows due to nonlinear mode coupling.Parametric dependences of the features of the intermittency are elaborated.Finally,associated turbulent heat transport is evaluated.展开更多
Turbulent transport resulting from drift waves,typically,the ion temperature gradient(ITG)mode and trapped electron mode(TEM),is of great significance in magnetic confinement fusion.It is also well known that turbulen...Turbulent transport resulting from drift waves,typically,the ion temperature gradient(ITG)mode and trapped electron mode(TEM),is of great significance in magnetic confinement fusion.It is also well known that turbulence simulation is a challenging issue in both the complex physical model and huge CPU cost as well as long computation time.In this work,a credible turbulence transport prediction model,extended fluid code(ExFC-NN),based on a neural network(NN)approach is established using simulation data by performing an ExFC,in which multi-scale multi-mode fluctuations,such as ITG and TEM turbulence are involved.Results show that the characteristics of turbulent transport can be successfully predicted including the type of dominant turbulence and the radial averaged fluxes under any set of local gradient parameters.Furthermore,a global NN model can well reproduce the radial profiles of turbulence perturbation intensities and fluxes much faster than existing codes.A large number of comparative predictions show that the newly constructed NN model can realize rapid experimental analysis and provide reference data for experimental parameter design in the future.展开更多
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.展开更多
This paper describe a numerical simulation method for the interaction between laser pulses and low density plasmas based on hydrodynamic approximation. We investigate Backward Raman Amplifying （BRA） experiments and ...This paper describe a numerical simulation method for the interaction between laser pulses and low density plasmas based on hydrodynamic approximation. We investigate Backward Raman Amplifying （BRA） experiments and their variants. The numerical results are in good agreement with experiments.展开更多
As an obstacle in high-performance discharge in future fusion devices,disruptions may do great damages to the reactors through causing strong electromagnetic forces,heat loads and so on.The drift waves in tokamak are ...As an obstacle in high-performance discharge in future fusion devices,disruptions may do great damages to the reactors through causing strong electromagnetic forces,heat loads and so on.The drift waves in tokamak are illustrated to play essential roles in the confinement performance as well.Depending on the plasma parameters and mode perpendicular wavelength,the mode phase velocity is either in the direction of electron diamagnetic velocity(namely,typical trapped electron mode)or in the direction of ion diamagnetic velocity(namely,the ubiquitous mode).Among them,the ubiquitous mode is directly investigated using gyro-fluid simulation associating with gyro-fluid equations for drift waves in tokamak plasmas.The ubiquitous mode is charactered by the short wavelength and propagates in ion diamagnetic direction.It is suggested that the density gradient is essential for the occurrence of the ubiquitous mode.However,the ubiquitous mode is also influenced by the temperature gradients and other plasma parameters including the magnetic shear and the fraction of trapped electrons.Furthermore,the ubiquitous mode may play essential roles in the turbulent transport.Meanwhile,the relevant parameters are scanned using a great number of electrostatic gyro-fluid simulations.The stability map is taken into consideration with the micro-instabilities contributing to the turbulent transport.The stability valley of the growth rates occurs with the assumption of the normalized temperature gradient equaling to the normalized density gradient.展开更多
基金supported by National Natural Science Foundation of China(No.12105227)。
文摘We have proposed a general numerical framework for plasma simulations on graphics processing unit clusters based on microscopic kinetic equations with full collision terms.Our numerical algorithm consistently deals with both long-range(classical forces in the Vlasov term)and short-range(quantum processes in the collision term)interactions.Providing the relevant particle masses,charges and types(classical,fermionic or bosonic),as well as the external forces and the matrix elements(in the collisional integral),the algorithm consistently solves the coupled multi-particle kinetic equations.Currently,the framework is being tested and applied in the field of relativistic heavy-ion collisions;extensions to other plasma systems are straightforward.Our framework is a potential and competitive numerical platform for consistent plasma simulations.
文摘The paper builds the high-current plasma beams model under different dimensions (1D, 2D, and 3D) by continuum (magnetohydrodynamics MHD) and statistical (Monte Carlo MC) mechanics under conditions of low pressures (10<sup>-3</sup> Pa). After detailed presentation of the model, two methods firstly have been analyzed in terms of plasma beam properties. Then, we compare the simulation results of MHD numerical simulation with MC stochastic particles simulation. Finally, through further analysis, it is demonstrated that integrated hybrid MHD and MC method (IMHDMC) provides an innovative practical tool to capture essential properties of high-current plasma beams.
文摘A high-energy electron beam generator is used to generate a plasma in atmosphere. Based on a Monte Carlo toolkit named GEANT4, a model including complete physics processes is established to simulate the passage of the electron beam in air. Based on the model, the characteristics of the electron beam air plasma are calculated. The energy distribution of beam electrons （BEs） indicates that high-energy electrons almost reside in the centre region of the beam, but low-energy electrons always live in the fringe area. The energy deposition is calculated in two cases, i.e., with and without secondary electrons （SEs）. Analysis indicates that the energy deposition of SEs accounts for a large part of the total energy deposition. The results of the energy spectrum show that the electrons in the inlet layer of the low-pressure chamber （LPC） are monoenergetic, but the energy spectrum of the electrons in the outlet layer is not pure. The SEs are largely generated at the outlet of the LPC. Moreover, both the energy distribution of BEs and the magnitude of the density of SEs are closely related to the pressure of LPC. Thus, a conclusion is drawn that a low magnitude of LPC pressure is helpful for reducing the energy loss in the LPC and also useful for greatly increasing the secondary electron density in dense air.
文摘A volume-weighting cloud-in-cell (VW-CIC) model is developed to implement the particle-in-cell (PIC) simulation in axially symmetric systems. This model gives a first-order accuracy in the cylindrical system, and it is incorporated into a PIC code. A planar diode with a finite-radius circular emitter is simulated with the code. The simulation results show that the VW-CIC model has a better accuracy and a lower noise than the conventional area-weighting cloud-in-cell (AW-CIC) model, especially on those points near the axis. The two-dimensional (2-D) space-charge-limited current density obtained from VW-CIC model is in better agreement with Lau’s analytical result. This model is more suitable for 2.5-D PIC simulation of axially symmetric plasmas.
基金supported by National Natural Science Foundation of China (No.10905044)
文摘For exploiting advantages of electron beam air plasma in some unusual applications, a Monte Carlo （MC） model coupled with heat transfer model is established to simulate the characteristics of electron beam air plasma by considering the self-heating effect. Based on the model, the electron beam induced temperature field and the related plasma properties are investigated. The results indicate that a nonuniform temperature field is formed in the electron beam plasma region and the average temperature is of the order of 600 K. Moreover, much larger volume pear-shaped electron beam plasma is produced in hot state rather than in cold state. The beam ranges can, with beam energies of 75 keV and 80 keV, exceed 1.0 m and 1.2 m in air at pressure of 100 torr, respectively. Finally, a well verified formula is obtained for calculating the range of high energy electron beam in atmosphere.
基金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.
基金co-supported by the National Key R & D Program for Intergovernmental International Scientific and Technological Innovation Cooperation,China(No. 2021YFE0116000)the National Natural Science Foundation of China (Nos. 12175032, 12102082 and 12275044),the National Natural Science Foundation of China+5 种基金the Belarusian Republican Foundation for Fundamental Research (No. 12211530449)the Fundamental Research Funds for the Central Universities of China (No. DUT22QN232)the S&T Program of Hebei, China (Nos. YCYZ202201 and 216Z1901G)the S&T Innovation Program of Hebei, China (Nos. SJMYF2022X18 and SJMYF2022X06)the Science and Technology Project of Hebei Education Department, China (No. ZC2023144)S&T Program of Langfang, China (No. 2022011039)
文摘In order to study the extraction and acceleration mechanism of the dual-stage grid,a three-dimensional model based on the Particle-In-Cell/Monte Carlo Collision(PIC/MCC)method is performed.Dual-stage grid ion thruster is a new type of electrostatic ion thruster,which can break through the limitations of traditional gridded ion thrusters,and greatly improve the specific impulse.The high performance also makes the grid sensitive to operating parameters.In this paper,the influence of grid parameters on xenon ion thruster’s performance in a wide range is systematically simulated,and the optimal operating condition is given.Both the over-focusing of the plume,and the transparency of the screen grid are improved,and the grid corrosion is reduced through simulation optimization.The specific impulse under the given working conditions is 9877.24 s and the thrust is 7.28 mN.Based on the simulation optimization,the limitation of the dual-stage grid is discussed.The grid performs well under high voltage conditions(>3000 V)but not well under low voltage conditions(<2000 V).Finally,since argon is cheaper and more advantageous in future engineering applications,the plasma distribution and grid extraction ability under xenon and argon are analyzed and compared to study the flexibility of the dual-stage grid ion thruster.The simulation results show that a set of optimal parameters is only applicable to the corresponding propellant,which needs to be optimized for different propellant types.
基金Project supported by the National Natural Science Foundation of China(Grant No.11305157)the Development Foundation of China Academy of Engineering Physics Laboratory(CAEPL)(Grant No.2013A0401019)
文摘Chirped wideband pump and seed pulses are usually considered for backward Raman amplification（BRA） in plasmas to achieve an extremely high-power laser pulse. However, current theoretical models only contain either a chirped pump or a chirped seed. In this paper, modified three-wave coupling equations are proposed for the BRA in the plasmas with both chirped wideband pump and seed. The simulation results can more precisely describe the experiments, such as the Princeton University experiment. The optimized chirp and bandwidth are determined based on the simulation to enhance the output intensity and efficiency.
基金supported by National Natural Science Foundation of China(No.11375107)Independent Innovation Foundation of Shandong University of China(No.2012TS067)
文摘In atmospheric radio-frequency （rf） discharges, the plasma parameters, such as electron density, sheath thickness and sheath voltage, are not easy to be probed experimentally, while the electrical characteristics, such as impedance, resistance and reactance, are relatively convenient to be measured. In this paper we presented a simple theoretical model derived from the fluid description of generated plasmas without considering the circuit model, to investigate the relationship between the plasma impedance and plasma parameters. By introducing a relaxation frequency, the plasma impedance could be predicted by formulas presented in this study, and the mean electron density and sheath thickness can also be calculated from the measured or simulated impedance and reactance, respectively.
基金supported by the National Natural Science Foundation of China(Grant No.11465013)the Natural Science Foundation of Jiangxi Province,China(Grant No.20151BAB212012)the International Science and Technology Cooperation Program of China(Grant No.2015DFA61800)
文摘A numerical study of the effect of water content on OH production in a pulsed-dc atmospheric pressure helium-air plasma jet is presented. The generation and loss mechanisms of the OH radicals in a positive half-cycle of the applied voltage are studied and discussed. It is found that the peak OH density increases with water content in air （varying from 0% to 1%） and reaches 6.3 x 10^18 m-3 when the water content is 1%. Besides, as the water content increases from 0.01% to 1%, the space-averaged reaction rate of three-body recombination increases dramatically and is comparable to those of main OH generation reactions.
基金The project partly supported by TEPCO Research Foundation
文摘Atomic and molecular processes relevant to the volumetric recombinationphenomena were investigated in a linear divertor plasma simulator MAP-Ⅱ. Volumetric recombinationis induced in He plasma by puffing of He or H_2. In the He puffing case, the reduction of the ionflux is dominated by the electron-ion recombination. In the H_2 puffing case, however, it isdominated by the molecule-assisted recombination (MAR), which is characterized by the disappearanceof the Helium Rydberg spectra and by the existence of the hydrogen negative ions. Currentachievement and the future prospect are described.
文摘Along with the introduction of the concept of dual-channel communication,we utilized the finite-difference time-domain（FDTD） method to simulate and measure the radiation pattern under certain plasma densities and plasma collision frequencies.Results show that under certain settings,the radiation pattern of a plasma antenna resembles that of a metallic antenna.In contrast to a metallic antenna,a plasma antenna possesses other functionalities,such as dynamic reconfiguration and digital controllability.The data from simulation are similar to the measurement results,indicating that column plasma antenna can realize dual-channel communication.This work confirms the viability of realizing dual-channel communication by column plasma antenna,which adds a new but promising method for modern intelligent communication.
基金Supported by the China ITER Plan Project Foundation under Grant Nos 2013GB113001 and 2014GB110004
文摘Plasma disruption is often an unavoidable aspect of tokamak operations. It may cause severe damage to in-vessel components such as the vacuum vessel conductors, the first wall and the divertor target plates. Two types of disruption, the hot-plasma vertical displacement event and the major disruption with a cold-plasma vertical displacement event, are simulated by the DINA code for HL-2M. The time evolutions of the plasma current, the halo current, the magnetic axis, the minor radius, the elongation as well as the electromagnetic force and eddy currents on the vacuum vessel during the thermal quench and the current quench are investigated. By comparing the electromagnetic forces before and after the disruption, we find that the disruption causes great damage to the vacuum vessel conductors. In addition, the hot-plasma vertical displacement event is more dangerous than the major disruption with the cold-plasma vertical displacement event.
基金supported by National Natural Science Foundation of China(Grant No.61378037)the Fundamental Research Funds for the Central Universities(Nos.2013B33614,2017B15214)+1 种基金the Research Funds of Innovation and Entrepreneurship Education Reform for Chinese Universities(No.16CCJG01Z004)the Changzhou Science and Technology Program(No.CJ20160027)
文摘The capacitively coupled radio frequency（CCRF）plasma has been widely used in various fields.In some cases,it requires us to estimate the range of key plasma parameters simpler and quicker in order to understand the behavior in plasma.In this paper,a glass vacuum chamber and a pair of plate electrodes were designed and fabricated,using 13.56 MHz radio frequency（RF）discharge technology to ionize the working gas of Ar.This discharge was mathematically described with equivalent circuit model.The discharge voltage and current of the plasma were measured atdifferent pressures and different powers.Based on the capacitively coupled homogeneous discharge model,the equivalent circuit and the analytical formula were established.The plasma density and temperature were calculated by using the equivalent impedance principle and energy balance equation.The experimental results show that when RF discharge power is 50–300 W and pressure is 25–250 Pa,the average electron temperature is about 1.7–2.1 e V and the average electron density is about 0.5？×10^17–3.6？×10^17m^-3.Agreement was found when the results were compared to those given by optical emission spectroscopy and COMSOL simulation.
基金supported by the the National MCF Energy R&D Program(No.2018YFE0304100)National Key Research and Development Program(Nos.2016YFA0400600,2016YFA0400601 and 2016YFA0400602)+1 种基金National Natural Science Foundation of China(Nos.11905220 and 11805273)supported by the U.S.Department of Energy(DE-AC02-09CH11466)。
文摘Explicit structure-preserving geometric particle-in-cell(PIC)algorithm in curvilinear orthogonal coordinate systems is developed.The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research.The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms,discrete exterior calculus,and explicit non-canonical symplectic integration.In addition to the truncated infinitely dimensional symplectic structure,the algorithm preserves exactly many important physical symmetries and conservation laws,such as local energy conservation,gauge symmetry and the corresponding local charge conservation.As a result,the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics.The algorithm has been implemented in the Sym PIC code,which is designed for highefficiency massively-parallel PIC simulations in modern clusters.The code has been applied to carry out whole-device 6 D kinetic simulation studies of tokamak physics.A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor.The state also admits a steady-state subsonic ion flow in the range of 10 km s-1,agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated.It is shown that high-n ballooning modes have larger growth rates than low-n global modes,and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2%level by the E×B flow generated by the instability.These results are consistent with early and recent electromagnetic gyrokinetic simulations.
基金supported by the National Key Research and Development Program of China(Grant Nos.2017YFE0301200 and 2017YFE0301201)partially by the National Natural Science Foundation of China(Grant Nos.11775069 and 11925501)the Liaoning Revitalization Talents Program(Grant No.XLYC1802009)。
文摘The structural characteristics of zonal flows and their roles in the nonlinear interaction of multi-scale multi-mode turbulence are investigated numerically via a self-consistent Landau-fluid model.The multi-mode turbulence here is composed of a shorter wavelength electromagnetic(EM)ion temperature gradient(ITG)mode and a Kelvin-Helmholtz(KH)instability with long wavelengths excited by externally imposed small-scale shear flows.For strong shear flow,a prominent periodic intermittency of fluctuation intensity except for dominant ITG component is revealed in turbulence evolution,which onset time depends on the ion temperature gradient and the shear flow amplitudes corresponding to different KH instabilities.It is identified that the intermittency phenomenon results from the zonal flow dynamics,which is mainly generated by the KH mode and back-reacts on it.It is demonstrated that the odd symmetric components of zonal flow(same symmetry as the external flow)make the radial parity of the KH mode alteration through adjusting the drift velocities at two sides of the resonant surface so that the KH mode becomes bursty first.Afterwards,the ITG intermittency follows due to nonlinear mode coupling.Parametric dependences of the features of the intermittency are elaborated.Finally,associated turbulent heat transport is evaluated.
基金supported by the National Key R&D Program of China (Nos. 2017YFE0301200 and 2017YFE0301201)partially by National Natural Science Foundation of China (Nos. 11775069 and 11925501)+1 种基金the Fundamental Research Funds for the Central Universities (No. DUT21GJ205)the Liao Ning Revitalization Talents Program (No. XLYC1802009)
文摘Turbulent transport resulting from drift waves,typically,the ion temperature gradient(ITG)mode and trapped electron mode(TEM),is of great significance in magnetic confinement fusion.It is also well known that turbulence simulation is a challenging issue in both the complex physical model and huge CPU cost as well as long computation time.In this work,a credible turbulence transport prediction model,extended fluid code(ExFC-NN),based on a neural network(NN)approach is established using simulation data by performing an ExFC,in which multi-scale multi-mode fluctuations,such as ITG and TEM turbulence are involved.Results show that the characteristics of turbulent transport can be successfully predicted including the type of dominant turbulence and the radial averaged fluxes under any set of local gradient parameters.Furthermore,a global NN model can well reproduce the radial profiles of turbulence perturbation intensities and fluxes much faster than existing codes.A large number of comparative predictions show that the newly constructed NN model can realize rapid experimental analysis and provide reference data for experimental parameter design in the future.
基金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.
文摘This paper describe a numerical simulation method for the interaction between laser pulses and low density plasmas based on hydrodynamic approximation. We investigate Backward Raman Amplifying （BRA） experiments and their variants. The numerical results are in good agreement with experiments.
基金Project partially supported by the National Natural Science Foundation of China(Grant Nos.12205035 and 11925501)also partially by the National Key Research and Development Program of China(Grant Nos.2017YFE0301200 and 2017YFE0301201).
文摘As an obstacle in high-performance discharge in future fusion devices,disruptions may do great damages to the reactors through causing strong electromagnetic forces,heat loads and so on.The drift waves in tokamak are illustrated to play essential roles in the confinement performance as well.Depending on the plasma parameters and mode perpendicular wavelength,the mode phase velocity is either in the direction of electron diamagnetic velocity(namely,typical trapped electron mode)or in the direction of ion diamagnetic velocity(namely,the ubiquitous mode).Among them,the ubiquitous mode is directly investigated using gyro-fluid simulation associating with gyro-fluid equations for drift waves in tokamak plasmas.The ubiquitous mode is charactered by the short wavelength and propagates in ion diamagnetic direction.It is suggested that the density gradient is essential for the occurrence of the ubiquitous mode.However,the ubiquitous mode is also influenced by the temperature gradients and other plasma parameters including the magnetic shear and the fraction of trapped electrons.Furthermore,the ubiquitous mode may play essential roles in the turbulent transport.Meanwhile,the relevant parameters are scanned using a great number of electrostatic gyro-fluid simulations.The stability map is taken into consideration with the micro-instabilities contributing to the turbulent transport.The stability valley of the growth rates occurs with the assumption of the normalized temperature gradient equaling to the normalized density gradient.
