基于流体模型和非平衡态电子能量分布函数的高功率微波气体击穿研究

用流体模型研究高功率微波气体击穿时, 电子能量分布函数常被假设为麦克斯韦分布形式, 此假设可能将给模拟结果带来较大的误差. 通过求解玻尔兹曼方程, 得到非平衡状态下的电子能量分布函数. 分别将上述两类分布函数引入到流体模型中, 对氩气击穿进行了数值模拟. 结果表明, 基于非平衡分布函数得到的击穿时间与粒子模拟结果符合得很好, 而当平均电子能量较低时, 麦克斯韦分布函数的高能尾部导致了较短的击穿时间. 最后, 采用非平衡分布函数计算了不同压强下的氩气击穿阈值, 发现其与实验结果基本符合. 关键词： 微波气体击穿 电子能量分布函数 流体模型 玻尔兹曼方程     

Short-pulse high-power microwave breakdown at high pressures

The fluid model is proposed to investigate the gas breakdown driven by a short-pulse(such as a Gaussian pulse) highpower microwave at high pressures.However,the fluid model requires specification of the electron energy distribution function(EEDF);the common assumption of a Maxwellian EEDF can result in the inaccurate breakdown prediction when the electrons are not in equilibrium.We confirm that the influence of the incident pulse shape on the EEDF is tiny at high pressures by using the particle-in-cell Monte Carlo collision(PIC-MCC) model.As a result,the EEDF for a rectangular microwave pulse directly derived from the Boltzmann equation solver Bolsig+ is introduced into the fluid model for predicting the breakdown threshold of the non-rectangular pulse over a wide range of pressures,and the obtained results are very well matched with those of the PIC-MCC simulations.The time evolution of a non-rectangular pulse breakdown in gas,obtained by the fluid model with the EEDF from Bolsig+,is presented and analyzed at different pressures.In addition,the effect of the incident pulse shape on the gas breakdown is discussed.     

沿面闪络流体模型电离参数粒子模拟确定方法

介绍了粒子模拟确定高功率微波介质沿面闪络击穿流体模型相关电离参数的方法.对粒子模拟方法 (包括带电粒子动力学方程、次级电子发射以及蒙特卡罗碰撞模型)和流体整体模型方法 (包括连续性方程和能量守恒方程)做了简介.基于自编的1D3V粒子模拟-蒙特卡罗碰撞程序给出了在高(低)气压、不同气体种类以及不同微波场强和微波频率下流体模型电离参数的粒子模拟结果,包括电离频率、击穿时间、平均电子能量、电子能量分布函数类型.研究结果表明:平均电子能量与电子能量分布函数类型关系不大;中低气压下,电子能量接近Maxwell分布,电子能量分布函数类型对电离参数几乎没有影响;中高气压下,电子能量分布函数类型对电离参数有重要影响,其依赖系数X趋于高阶形式.不同气体的电子能量分布函数类型不同,需要利用粒子模拟对电子能量分布函数类型进行标定.同时,电子能量分布函数依赖系数与微波场强和频率也有关系,其随微波场强增加而增大,随微波频率增加而减小.在给定考察范围(微波场强在7 MV/m以下,微波频率在40 GHz以内),中低气压下,平均电子能量随微波场强增加而迅速增大,电离频率随微波场强增加先增大后降低,平均电子能量随微波频率增加而降低,电离频率随微波频率增加先增加后降低;高气压下,平均电子能量随微波场强增加而缓慢增大,电离频率随微波场强增加而增大,微波频率对平均电子能量和电离频率影响不大.     

高功率微波单一气体及混合气体击穿特性

综合考虑高功率微波对电子的加速过程以及电子与气体分子的碰撞过程,建立了单一气体与混合气体击穿过程的蒙特卡罗仿真模型,编写了三维蒙特卡罗仿真程序(3D-MCC)。针对单一气体Ar和N2以及混合气体N2/O2展开研究,仿真了气体雪崩击穿电子云形成过程,对比分析了不同气体电子能量分布函数随压强的变化规律。发现了Ar击穿特性受电子能量分布函数影响较大,而N2击穿特性受电子能量分布函数影响较小。通过分析平均电子能量以及电子密度随时间的变化过程,得到了Ar和N2击穿时间,并通过与流体模型计算得到的击穿时间比对分析验证了3D-MCC模型的正确性。在真空腔体内开展了S波段高功率微波大气击穿实验,测量得到了场强为6.38 kV/cm时不同压强下的大气击穿时间。通过在辐射源与真空腔体之间增加聚焦透镜,大大减小了壁效应的影响,并且采用模型仿真得到的大气击穿时间与实验结果吻合较好。     

Effect of air breakdown on microwave pulse energy transmission

The energy transmission of the long microwave pulse for the frequency of 2.45 GHz and 5.8 GHz is studied by using the electron fluid model, where the rate coefficients are deduced from the Boltzmann equation solver named BOLSIG+. The breakdown thresholds for different air pressures and incident pulse parameters are predicted, which show good agreement with the experimental data. Below the breakdown threshold, the transmitted pulse energy is proportional to the square of the incident electric field amplitude. When the incident electric field amplitude higher than the breakdown threshold increases,the transmitted pulse energy decreases monotonously at a high air pressure, while at a low air pressure it first decreases and then increases. We also compare the pulse energy transmission for the frequency of 2.45 GHz with the case of 5.8 GHz.     

不同气压下介质表面高功率微波击穿的数值模拟

介绍了用于模拟介质表面高功率微波击穿的粒子模拟-蒙特卡罗碰撞方法,并采用该方法模拟研究了氩气环境不同气压下的介质表面高功率微波击穿过程,获得了该击穿过程中粒子数量和电子平均能量的时间变化图像,并得到了击穿延迟时间。数值模拟结果发现:在低气压下,次级电子倍增的作用比较明显,但电子数量在次级电子倍增饱和后的增速较低,击穿延迟时间较长;随着气压的升高,次级电子倍增的影响逐渐变小,气体电离逐渐占主导地位,击穿延迟时间逐渐变短;在高气压下,由于介质表面吸收沉积电子而呈负电性,次级电子倍增消失,击穿延迟时间由气体碰撞电离来决定。     

流注放电低温等离子体中电子输运系数计算的蒙特卡罗模型

流体或者粒子-流体混合数值仿真是研究流注放电基本物理机制的常用手段,而精确的电子输运系数是保证其仿真正确性的必要前提.鉴于现有电子输运系数求解工具存在一定缺陷,本文开发了采用蒙特卡罗方法求解低温等离子体中电子输运系数的仿真工具,测试表明其准确性和精确度均较高.研究了氮氧气体混合比及大气压下三体碰撞吸附对电子输运系数的影响.氮气中流注放电仿真表明,流体仿真中采用本模型改进后的电子输运系数可显著改善流注通道内部的等离子体参数分布.     

激光烧蚀金属靶时气体电离分析

我们使用Ｎｄ：ＹＡＧ激光器烧蚀金属Ａｌ靶获得等离子体,利用光谱时－空分辨技术,关于环境气体分解及其对等离子体辐射特征的影响进行了研究。使用的气体是Ａｒ,Ａｉｒ,Ｎ２。结果发现,在等离子体形成初期,等离子体时－空分辨谱中有大量的环境气体离子谱线。基于Ａｌ等离子体不同激光能量下的时间－空间分辨谱。对激光烧蚀金属靶时气体电离现象进行了分析,对气体分解造成的影响进行了简单的讨论。     

Substantiation of the two-temperature kinetic model by comparing calculations within the kinetic and fluid models of the positive column plasma of a dc oxygen discharge

Results from kinetic and fluid simulations of the positive column plasma of a dc oxygen discharge are compared using commercial CFDRC software (), which enables one to perform numerical simulations in an arbitrary 3D geometry with the use of both the fluid equations for all the components (fluid model) and the kinetic equation for the electron energy distribution function (kinetic model). It is shown that, for both the local and nonlocal regimes of the formation of the electron energy distribution function (EEDF), the non-Maxwellian EEDF can satisfactorily be approximated by two groups of electrons. This allows one to take into account kinetic effects within the conventional fluid model in the simplest way by using the proposed two-temperature approximation of the nonequilibrium and nonlocal EEDF (2T fluid model).     

�����Թ�ŵ����������̵�PIC/MCCģ��

采用粒子模拟与蒙特卡罗模型相结合(PIC/MCC)的方法,建立了氮分子气体辉光放电自洽的混合模型,其中带电粒子在电场中的运动及其产生的自洽场由PIC方法的静电模型描述,粒子间的碰撞过程由MC方法描述。模拟了放电过程及带电粒子(e、N2+、N+)在整个放电空间的行为及与之对应的自洽电场和电势的分布,通过计算带电粒子在负辉区的几率分布函数,讨论了带电粒子(e、N2+、N+)在负辉区的行为和分布特征。计算的分子离子密度与实验结果一致。     

Effect of molecular nitrogen on the electron mobility in a mixture of argon and optically excited sodium vapor

A parametric study of the electron energy distribution function (EEDF) and the electron mobility in the mixture Na + Ar + N₂ is carried out. An analysis is made of the conditions that obtain in a photoplasma when the detachment of the mean electron energy from the neutral gas temperature is due to superelastic collisions (collisions of the second kind) with excited sodium atoms. The case of low ionization of the medium at low vibrational temperatures of the ground state of the nitrogen molecules is considered. To find the EEDF a numerical solution of the Boltzmann transport equation is carried out. It is found that in the indicated mixture the presence of nitrogen leads to a depletion of the EEDF in the region of efficient vibrational excitation of the molecules and promotes the formation of inversion in the EEDF ∂f(ɛ)/∂ɛ>0 in the energy range corresponding to the Ramsauer minimum in the cross section of elastic collisions of electrons with the argon atoms. It is shown that the nonequilibrium character of the EEDF leads to a complicated dependence of the electron mobility on the partial ratios of the components of the mixture, the degree of ionization of the medium, and the population of the resonantly excited sodium atoms. Zh. Tekh. Fiz. 69, 14–19 (April 1999)     

Major elements analysis in bituminous coals under different ambient gases by laser-induced breakdown spectroscopy with PLS modeling

Three major elements, carbon, hydrogen, and nitrogen, in twenty-four bituminous coal samples, were measured by laser-induced breakdown spectroscopy. Argon and helium were applied as ambient gas to enhance the signals and eliminate the interference of nitrogen from surrounding air. The relative standard deviation of the related emission lines and the performance in the partial least squares (PLS) modeling were compared for different ambient environments. The results showed that argon not only improved the intensity, but also reduced signal fluctuation. The PLS model also had the optimal performance in multi-element analysis using argon as ambient gas. The root mean square error of prediction of carbon concentration decreased from 4.25% in air to 3.49% in argon, while the average relative error reduced from 4.96% to 2.98%. Hydrogen line demonstrated similar improvement. Yet, the nitrogen lines were too weak to be detected even in an argon environment which suggested the nitrogen signal measured in air come from the breakdown of nitrogen molecules in the atmosphere.     

Determination of the electron energy distribution in a plasma from the first and second derivatives of the probe current

An approach to the reconstruction of the electron energy distribution function (EEDF) as the sum of the second derivative of the electron current to the probe and some value proportional to its first derivative is proposed (hereafter we mean derivatives with respect to voltage). Solutions to model problems for typical electron distributions in a plasma show that this approach lowers the systematic error of EEDF reconstruction by several times in comparison with the conventional techniques using Langmuir and diffusion probes. The approach is applicable in a wide range of the ratio of the probe radius to the free path of electrons. It can be used for the determination of fast nonequilibrium electron distribution in neon discharge afterglow and also for the evaluation of the Maxwellian electron temperature in the low-energy range of nitrogen discharge afterglow.     

Enhancement of the electric field threshold for microwave breakdown by crossed magnetic field

The enhancement of breakdown threshold is of benefit to the high power microwave transmission. We propose a magnetic field in the transverse direction to the electric field to enhance the breakdown threshold. A theory of electric field threshold with crossed magnetic field for short pulse is developed, and verified by particle-in-cell/Monte Carlo collision (PIC/MCC) simulations. The result shows that the crossed magnetic field can enhance the breakdown threshold significantly.     

Particle in Cell/Monte Carlo Collision Method for Simulation of RF Glow Discharges: Effect of Super Particle Weighting

A parallel Particle in Cell/Monte Carlo Collision (PIC/MCC) numerical code for glow discharge plasma simulations is developed and verified. This method is based on simultaneous solution of the Lorentz equations of motion of super particles, coupled with the Poisson's equation for electric field. Collisions between the particles are modelled by the Monte Carlo method. Proper choice of particle weighting is critically important in order to perform adequate and efficient PIC simulations of plasma. Herein, effects of particle weighting on the simulations of capacitive radio‐frequency argon plasma discharges are studied in details. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

PIC/MCC Simulation of Electron and Ion Currents to Spherical Langmuir Probe

The Particle In Cell/Monte Carlo Collisions (PIC/MCC) simulation was used for the calculation of electron and ion currents to a spherical Langmuir (electrostatic) probe. This simulation took into account the collisions of collected charged particles with neutral gas particles around the probe and it can calculate the probe currents at higher neutral gas pressures. The improvements of usual simulation techniques enabled to speed up the simulation and to calculate the probe current even for neutral gas pressures above 1 kPa. The simulations were carried out for two cases: i) probe with radius of 0.5 mm in non‐thermal plasma with high electron temperature, ii) probe with radius of 10 µm in afterglow plasma with low electron temperature. The influence of probe radius on electron probe current was also studied. The simulations showed that thick sheath limit of OML theory provides incorrect values of probe current for probes with radii larger than 200 µm at plasma parameters considered even at very low neutral gas pressures. The probe characteristics were calculated for probe with 0.5 mm radius for pressures up to 500 Pa and for probe with 10 μm radius for pressures up to 3 kPa. The influence of collisions on electron and ion probe current was demonstrated and the procedure for determination of electron and ion densities from the probe measurement at higher pressures was developed. The results from PIC/MCC simulations were compared with results from continuum theory. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-Maxwellian electron energy distribution function in He,He/Ar,He/Xe/H<Subscript>2</Subscript> and He/Xe/D<Subscript>2</Subscript> low temperature afterglow plasma

Experimental studies of the electron energy distribution function “EEDF” under well defined conditions in flowing afterglow plasma, using a Langmuir probe are reported. The EEDF is measured in He₂ ⁺ and Ar⁺ dominated plasmas and in XeH⁺ and XeD⁺ dominated recombining plasmas. He is used as a buffer gas at medium pressures in all experiments (1600 Pa, 250 K). The deviation of the measured EEDF from Maxwellian distribution is shown to depend on plasma composition and on the processes governing the plasma decay. The influence of energetic electrons produced during the plasma decay on the body and tail of the EEDF is observed. The mechanism of energy balance in afterglow plasma is discussed.     

A Particle‐in‐Cell Simulation of the Breakdown Mechanism in Microdischarges with an Improved Secondary Emission Model

In this paper, the failure of the breakdown voltage from the Paschen's law at extremely small electrode separations is studied. The electrical breakdown in microgaps occurs at the voltages far below the Paschen curve minimum breakdown limit and the modified Paschen curve should be used. Offered explanation for the departure from the Paschen's law at small gap spacings is based on the increasing of the yield of the secondary electrons. The high electric fields existing in small gaps may enhance the secondary electron yield and this would lead to a lowering of the breakdown voltage and to the departure from the Paschen's law. Particlein‐cell/Monte‐Carlo (PIC/MCC) simulations with a new secondary emission model have been performed to estimate the importance of this mechanism in the discharge breakdown. Obtained simulation results suggest that deviations from the Paschen curve across the micron and submicorn gap spacing can be attributed to the ion‐enhanced field emissions. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)