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

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

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

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

Validity of the two-term Boltzmann approximation employed in the fluid model for high-power microwave breakdown in gas

The electron energy distribution function （EEDF）, predicted by the Boltzmann equation solver BOLSIG＋ based on the two-term approximation, is introduced into the fluid model for simulating the high-power microwave （HPM） breakdown in argon, nitrogen, and air, and its validity is examined by comparing with the results of particle-in-cell Monte Carlo collision （PIC/MCC） simulations as well as the experimental data. Numerical results show that, the breakdown time of the fluid model with the Maxwellian EEDF matches that of the PIC/MCC simulations in nitrogen; however, in argon under high pressures, the results from the Maxwellian EEDF were poor. This is due to an overestimation of the energy tail of the Maxwellian EEDF in argon breakdown. The prediction of the fluid model with the BOLSIG＋ EEDF, however, agrees very well with the PIC/MCC prediction in nitrogen and argon over a wide range of pressures. The accuracy of the fluid model with the BOLSIG＋ EEDF is also verified by the experimental results of the air breakdown.     

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.     

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

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

高功率微波与等离子体相互作用理论和数值研究

研究高功率微波与等离子体的相互作用,对于微波放电和电磁兼容研究均具有重要意义.基于波动方程、等离子体的流体力学方程以及波尔兹曼方程,建立高功率微波脉冲与等离子体相互作用的理论模型,并结合等离子体的特征参数,采用时域有限差分方法分析了等离子体电子密度和高功率微波传输特性的变化.结果表明,由于高功率微波的电子加热作用,等离子体中的非线性效应明显,发生击穿使得等离子体电子密度增大,从而导致微波的反射增强,透过率降低.所提出的模型和相关结果对于高功率微波和电磁脉冲防护具有指导意义.     

Numerical simulation of super-short pulsed discharge in Helium with particle-in-cell Monte--Carlo collisions technique

This paper reports that a simulation of glow discharge in pure helium gas at the pressure of 1.333×10³ Pa under a high-voltage nanosecond pulse is performed by using a one-dimensional particle-in-cell Monte Carlo collisions (PIC--MCC) model. Numerical modelling results show that the cathode sheath is much thicker than that of anode during the pulse discharge, and that there exists the phenomenon of field reversal at relative high pressures near the end of the pulse, which results from the cumulative positive charges due to their finite mobility during the cathode sheath expansion. Moreover, electron energy distribution function (EEDF) and ion energy distribution function (IEDF) have been also observed. In the early stage of the pulse, a large amount of electrons can be accelerated above the ionization threshold energy. However, in the second half of the pulse, as the field in bulk plasma decreases and thereafter the reverse field forms due to the excessive charges in cathode sheath, although the plasma density grows, the high energy part of EEDF decreases. It concludes that the large volume non-equilibrium plasmas can be obtained with high-voltage nanosecond pulse discharges.     

Microwave breakdown in slots

The properties of microwave-induced breakdown of air in narrow metallic slots are investigated, both theoretically and experimentally, with emphasis on factors important for protection against transmission of incident high-power microwave radiation. The key factors investigated are breakdown power threshold, breakdown time, peak-leakage power, and total transmitted energy, as functions of incident pulse shape and power density. The theoretical investigation includes estimates of the electric field intensification in narrow slots and basic breakdown plasma modeling. New results important for application to the high-power microwave field, such as the influence of pulse shape on breakdown time and peak-leakage power, are presented. The experimental investigation comprises a set of slot breakdown experiments at atmospheric pressure, which are analyzed to extract key parameters, such as transmission cross section, breakdown time, peak leakage power, and transmitted energy. The experimental data is compared and shown to be in good agreement with results obtained in the theoretical investigation.     

介质表面高功率微波击穿的数值模拟

研究了用于模拟高功率微波条件下介质表面击穿的静电PIC-MCC模型,并通过自行编写的数值模拟程序模拟了真空及不同气压条件下介质表面击穿过程中的次级电子倍增和气体电离等过程.模拟结果发现,在真空及低气压条件下,电子的主要来源是次级电子倍增,电子数量以两倍于入射场的频率振荡;在高气压情况下,电子的主要来源是气体电离. 关键词： 介质表面击穿 高功率微波 数值模拟 次级电子倍增     

高功率微波窗内外表面闪络击穿流体模拟研究

建立理论模型,将电磁场时域有限差分方法与等离子体流体模型结合,编制一维电磁场与等离子流体耦合程序,数值研究了3 GHz高功率微波窗内外表面闪络击穿的不同物理过程.研究结果表明:外表面闪络击穿中,输出微波脉宽缩短(未完全截止),窗体前均方根场强呈驻波分布,波节与波腹位置不变,窗体外表面形成有一层高密(约10~(21)·m~(-3)量级)极薄(约mm量级)等离子体(扩散缓慢),入射波可部分透过该薄层等离子体,脉宽缩短主要源于等离子体吸收效应;降低初始等离子体密度、厚度、入射波场强及缩短入射波脉宽等方式,可不同程度地改善输出脉宽缩短效应.内表面闪络击穿中,窗体前均方根场强亦出现驻波分布f但波节与波腹位置随时间变化),等离子体向波源方向运动;强释气下,输出脉宽缩短(未完全截止),形成多丝状高密(约10~(21)·m~(-3)量级)极薄(约mm量级)等离子体区域(扩散缓慢),间距1/4微波波长,脉宽缩短主要源于等离子体吸收效应;弱释气、低场强下,脉宽缩短有所改善(但最终截止),形成多带状致密(约10~(18)·m~(-3)量级)略厚(mm-cm量级)等离子体区域(扩散较快),间距1/4波长,脉宽缩短主要源于等离子体吸收效应;弱释气、高场强下,脉宽缩短严重(很快截止),形成块状高密(约10~(21)·m~(-3)量级)较厚(约cm量级)等离子体区域(扩散迅速),脉宽缩短主要源于等离子体反射效应.     

Similarity principle of microwave argon plasma at low pressure

In order to validate the similarity principle of microwave breakdown, a two-dimensional(2 D) fluid model of lowpressure microwave argon plasma is established and solved by the finite-element method. Proportional conditions are used in this model to build three different breakdown processes that meet the premise of a similarity principle, and these breakdown processes are called "similar cases" in this paper. Similar cases have proportionately sized breakdown regions,where the ratio of frequency of incident microwave f to gas pressure p( f/p), and the reduced field E/p in them are kept the same. All the important physical parameters such as electron density, electron temperature, and reduced electric field can be obtained from the simulation of this model. The results show that the parameters between similar cases are in constant ratio without changing with time, which means that the similarity principle is also valid in microwave breakdown.     

Theoretical Studies of the Breakdown Characteristics at Microwave Frequencies

This paper describes several theoretical approaches for studies of the breakdown characteristics for a few atomic and molecular gases at microwave frequencies in small gaps. Numerical studies have been carried out extensively over a wide range of pressures for a microwave frequency of 2.45 GHz using both the kinetic and fluid approaches and compared with the results obtained by an analytical-experimental approach. The obtained results illustrate the physical phenomena that occur during breakdown with emphasis on the determination of the breakdown voltage and its dependence on the gap size and gas pressure at microwave frequency. Good agreement between the results achieved by using different approaches suggests that they can be applied even to gaps of a few hundred micrometers and high pressures.     

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

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

高功率微波沿面闪络击穿机制粒子模拟

针对高功率微波介质沿面闪络击穿物理过程,首先建立了理论模型,包括:动力学方程、粒子模拟算法、二次电子发射, 以及电子与气体分子蒙特卡罗碰撞模型、电子碰撞介质表面退吸附气体分子机制;其次,基于理论模型,编制了1D3V PIC-MCC程序,分别针对真空二次电子倍增、高气压体电离击穿和低气压面电离击穿过程,运用该程序仔细研究了电子和离子随时间演化关系、电子运动轨迹、电子及离子密度分布、空间电荷场时空分布、电子平均能量、碰撞电子平均能量、碰撞电子数目随时间演化关系、电子能量分布函数、平均二次电子发射率以及能量转换关系。研究结果表明:真空二次电子倍增引发的介质表面沉积功率只能达到入射微波功率1%左右的水平,不足以击穿;气体碰撞电离主导的高气压体电离击穿,是由低能电子（eV量级）数目指数增长到一定程度导致的,形成位置远离介质表面,形成时间为s量级;低气压下的介质沿面闪络击穿,是在二次电子倍增和气体碰撞电离共同作用下,由于数目持续增长的高能电子（keV量级）碰撞介质沿面导致沉积功率激增而引发的,形成位置贴近介质沿面,形成时间在ns量级。     

微波输出窗内表面闪络击穿3维全电磁等离子体流体模拟北大核心CSCD

通过建立电磁场等离子体流体耦合物理模型,基于自主研发的3维全电磁粒子模拟大规模并行程序NEPTUNE3D,编制了3维电磁场与等离子流体耦合程序模块,对1.3GHz高功率微波窗内表面闪络击穿物理过程进行了数值模拟。研究结果表明:微波窗内侧表面形成的等离子体构型与初始种子电子分布形式密切相关。中心点源分布下,等离子体发展为"蘑菇"形状,输出微波脉冲缩短并不严重,等离子体吸收微波功率大于反射微波功率;面源分布下,等离子体发展为"帽子"形状,输出微波脉冲缩短严重,输出微波完全截断,开始阶段等离子体吸收微波功率占优,待等离子体密度增加到一定程度后,反射微波功率占优。通过降低窗体表面场强、表面释气率及初始种子电子密度等方法,可不同程度地延长输出微波脉冲宽度。窗体表面不同气体层厚度对闪络击穿下的输出微波脉冲宽度影响不大。     

Variations in the electron energy distribution function in discharge plasma of a unipolar gas breakdown

It is experimentally demonstrated that the electron energy distribution function (EEDF) and electron density in plasma of a little-studied pulse discharge of a unipolar gas breakdown (UGB) significantly vary in time. The revealed behavior of its plasma parameters is explained based on the proposed physical model of a UGB discharge.     

射频感应耦合Ar-N2等离子体物理特性的Langmuir探针测量及理论研究

分别通过Langmuir探针测量和动力学模型模拟方法研究了射频感应耦合Ar-N₂等离子体中电子能量分布、电子温度、电子密度等物理量随N₂含量的变化规律.实验研究结果表明：电子能量分布呈现出非Maxwell型分布,并由双温分布向三温分布过渡;电子温度在不同的气压下随N₂含量的增加呈现出不同的变化规律.在放电气压小于1.3 Pa时,电子温度随N₂含量的增加而下降;当气压大于1.3 Pa时,电子温度随N₂关键词： 感应耦合等离子体 2混合气体放电')" href="#">Ar-N₂混合气体放电 电子能量分布 Langmuir探针     

微波大气击穿阈值的理论研究

本文通过对使用有效场强(或均方根场强)得到的微波大气击穿阈值表达式进行讨论, 指出其推导中所做的假设及这些假设应用到微波大气击穿过程中存在的问题. 然后分别使用解析理论和数值模拟对微波大气击穿过程中的有效电子温度变化过程和击穿阈值进行研究, 并将其与直流电场进行比较. 分析发现在高气压下, 电子能量转移频率高, 有效电子温度随电场大幅振荡, 由于电离频率随有效电子温度的增长率大于电子能量损失随有效电子温度的增长率, 因此在高气压时, 微波大气击穿阈值低于使用有效场强的击穿阈值. 通过大量分析, 给出了理论推导和数值模拟得到的微波大气击穿阈值拟合表达式.     

开放空间重复频率微波脉冲击穿概率理论与实验研究

综合考虑有效初始电子产生理论、雪崩电子击穿理论等过程中的击穿延迟时间,探讨了开放空间微波脉冲的击穿延时概率分布,提出了重复频率微波脉冲击穿概率模型,定义了基于概率模型的微波脉冲击穿阈值。利用S波段微波准光学反射聚焦系统对一定气压大气击穿过程进行了模拟,监测击穿放电发光时刻作为击穿时间,分别在铯137放射源存在与否情况下开展了系列实验。研究结果表明,提高种子电子产生率相较于提高电离率是增大脉冲击穿概率更有效的方法;重复频率过程中,若存在累积效应,击穿延时概率分布曲线将左移并趋于稳定,击穿后的气体在短时间内容易再次击穿。     

高功率微波输出窗内侧击穿动力学的PIC/MCC模拟研究

高功率微波在受控热核聚变加热、微波高梯度加速器、高功率雷达、定向能武器、超级干扰机及冲击雷达等方面有着重要的应用.本文针对高功率微波输出窗内侧氩气放电击穿过程,建立了二次电子倍增和气体电离的一维空间分布、三维速度分布（1D3V）模型,并开发了相应的PIC/MC程序代码.研究了气压、微波频率、微波振幅对放电击穿的影响.结果表明：在真空情况下,介质窗放电击穿只存在二次电子倍增过程;在低气压和稍高气压时,二次电子倍增和气体电离共存;在极高气压时,气体电离占主导.给出了不同气压下电子、离子的密度和静电场的空间分布.此外还观察到,在500 mTorr时,随着微波振幅或微波频率的变化,气体电离出现的时刻和电离产生的等离子体峰值位置有较大差异,尤其是当微波频率（GHz）在数值上是微波振幅（MV/m）的2倍时,气体电离出现的较早.