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.展开更多
This paper presents the concept of a Dielectric-lined Multiwave Cerenkov Generator producing high power millimeter waves, which has been investigated with a two and onehalf dimensional electromagnetic relativistic Par...This paper presents the concept of a Dielectric-lined Multiwave Cerenkov Generator producing high power millimeter waves, which has been investigated with a two and onehalf dimensional electromagnetic relativistic Particle-in-Cell (PIC) plasma simulation code. Themodified device can operate in a lower diode-voltage regime with much higher radiation efficiencyand slight downshift of operation frequency. There exist the optima for the permittivity of thedielectric liner and for the magnitude of the guiding magnetic field. The required intensity of theguiding field is reduced by the introduction of the liner. The enhanced propagation of the electronbeam is studied in the presence of the liner.展开更多
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.展开更多
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.展开更多
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.展开更多
Some nonlinear behavior in electron cyclotron resonance plasma was investigated using a two-dimension hybrid-mode with self-consistent microwave absorption. The saturation,oscillations of plasma parameters (plasma den...Some nonlinear behavior in electron cyclotron resonance plasma was investigated using a two-dimension hybrid-mode with self-consistent microwave absorption. The saturation,oscillations of plasma parameters (plasma density, potential, electron temperature) versus operating conditions (pressure, power) are discussed. Our simulation results are consistent qualitatively with many experimental measurements.展开更多
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.展开更多
Transport simulation of ECRH H-mode experiments on HL-2A tokamak is carried out using ONETWO code, the GLF23 and PEDESTAL models, along with TORAY code for ECRH. It is found that the initial electron and ion temperatu...Transport simulation of ECRH H-mode experiments on HL-2A tokamak is carried out using ONETWO code, the GLF23 and PEDESTAL models, along with TORAY code for ECRH. It is found that the initial electron and ion temperature profiles affect L-H transition significantly, and larger initial temperature gradient at the edge plasma benefits the transition. The simulation results show that it is possible to achieve ECRH H-mode with appropriate initial electron and ion temperature profiles under present discharge conditions on HL-2A tokamak. In addition, the pedestal density, electron temperature and pedestal width are predicted, and the evolutions of electron and ion temperature profile are calculated.展开更多
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.展开更多
The existence of two diffe1：ent discharge modes has been verified in an rf （radio-frequency） atmospheric pressure glow discharge （APGD） by Shi [J. Appl. Phys. 97, 023306 （2005）]. In the first mode, referred to ...The existence of two diffe1：ent discharge modes has been verified in an rf （radio-frequency） atmospheric pressure glow discharge （APGD） by Shi [J. Appl. Phys. 97, 023306 （2005）]. In the first mode, referred to as a mode, the discharge current density is relatively low and the bulk plasma electrons acquire the energy due to the sheath expansion. In the second mode, termed γ mode, the discharge current density is relatively high, the secondary electrons emitted by cathodc under ion bombardment in the cathode sheath region play an important role in sustaining the discharge. In this paper, a one-dimensional self-consistent fluid model for rf APGDs is used to simulate the discharge mechanisms in the mode in helium discharge between two parallel metallic planar electrodes. The results show that as the applied voltage increases, the discharge current becomes greater and the plasma density correspondingly increases, consequentially the discharge transits from the a mode into the γ mode. The high collisionality of the APGD plasma results in significant drop of discharge potential across the sheath region, and the electron Joule heating and the electron collisional energy loss reach their maxima in the region. The validity of the simulation is checked with the available experimental and numerical data.展开更多
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.展开更多
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.展开更多
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.展开更多
基金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.
文摘This paper presents the concept of a Dielectric-lined Multiwave Cerenkov Generator producing high power millimeter waves, which has been investigated with a two and onehalf dimensional electromagnetic relativistic Particle-in-Cell (PIC) plasma simulation code. Themodified device can operate in a lower diode-voltage regime with much higher radiation efficiencyand slight downshift of operation frequency. There exist the optima for the permittivity of thedielectric liner and for the magnitude of the guiding magnetic field. The required intensity of theguiding field is reduced by the introduction of the liner. The enhanced propagation of the electronbeam is studied in the presence of the liner.
文摘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.
文摘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.
文摘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.
文摘Some nonlinear behavior in electron cyclotron resonance plasma was investigated using a two-dimension hybrid-mode with self-consistent microwave absorption. The saturation,oscillations of plasma parameters (plasma density, potential, electron temperature) versus operating conditions (pressure, power) are discussed. Our simulation results are consistent qualitatively with many experimental measurements.
基金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 Magnetic Confinement Fusion Science Program of China(Nos.2013GB111005,2012GB105001,2011GB105004)National Natural Science Foundation of China(No.10875037)
文摘Transport simulation of ECRH H-mode experiments on HL-2A tokamak is carried out using ONETWO code, the GLF23 and PEDESTAL models, along with TORAY code for ECRH. It is found that the initial electron and ion temperature profiles affect L-H transition significantly, and larger initial temperature gradient at the edge plasma benefits the transition. The simulation results show that it is possible to achieve ECRH H-mode with appropriate initial electron and ion temperature profiles under present discharge conditions on HL-2A tokamak. In addition, the pedestal density, electron temperature and pedestal width are predicted, and the evolutions of electron and ion temperature profile are calculated.
基金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.
基金Project supported by the National Natural Science Foundation of China（Grant Nos 50528707 and 50537020）.
文摘The existence of two diffe1：ent discharge modes has been verified in an rf （radio-frequency） atmospheric pressure glow discharge （APGD） by Shi [J. Appl. Phys. 97, 023306 （2005）]. In the first mode, referred to as a mode, the discharge current density is relatively low and the bulk plasma electrons acquire the energy due to the sheath expansion. In the second mode, termed γ mode, the discharge current density is relatively high, the secondary electrons emitted by cathodc under ion bombardment in the cathode sheath region play an important role in sustaining the discharge. In this paper, a one-dimensional self-consistent fluid model for rf APGDs is used to simulate the discharge mechanisms in the mode in helium discharge between two parallel metallic planar electrodes. The results show that as the applied voltage increases, the discharge current becomes greater and the plasma density correspondingly increases, consequentially the discharge transits from the a mode into the γ mode. The high collisionality of the APGD plasma results in significant drop of discharge potential across the sheath region, and the electron Joule heating and the electron collisional energy loss reach their maxima in the region. The validity of the simulation is checked with the available experimental and numerical data.
基金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.
基金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.
基金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.
