HL-2Aװ�����������˵��о�

用密度调制的方法研究了等离子体中粒子输运问题。采用了注入脉冲式补充送气和超声分子束两种不同的密度调制方法。在HL-2A装置常规欧姆放电的情况下,运用有限差分法和Nagashima矩阵技术,求解了粒子平衡方程。计算出了粒子的输运系数(对流速度v和扩散系数D)。研究了粒子输运系数与等离子体线平均密度之间的关系。实验结果表明,在欧姆放电的情况下,等离子体芯部的粒子对流速度方向始终是向内的,并且密度低时,粒子输运系数(粒子扩散系数D和对流速度v)较大;密度高时,粒子输运系数较小。     

纳秒脉冲下变压器油两相流注放电仿真研究

为揭示液体电介质击穿过程中形成的气体放电通道对液体电介质放电过程的影响,以针—板电极间隙变压器油为研究对象,基于等离子体流体力学模型,引入了液体电介质放电过程中气相放电通道对电离机制及自由电荷迁移率的影响,建立了用于模拟脉冲电压下液体电介质放电过程的两相流体模型,仿真研究了纳秒脉冲下针板电极流注放电的起始与发展过程。仿真结果表明：采用Heaviside方程可以在模型的不同区域同时实现气相物理过程和液相物理过程的模拟与计算。气相物理过程的引入导致流注尾部电场显著降低,流注头部电场进一步增强,使流注通道的发展速度要高于传统液相模型,有助于加深对纳秒脉冲下液体电介质中预击穿流注的起始、发展过程的认识和理解。     

射频辉光放电CH4等离子体中电子的能量分布

通过求解Lorentz简化的玻尔兹曼方程,得到射频放电CH₄等离子体中电子的能量分布函数.求解过程中使用一个简化的射频电场模型代替泊松方程求解放电电场.共计包含6类环境气体及27种电子碰撞反应.通过EEDF对等离子体中的电子反应率系数、电子平均能量、电子的传输率系数等进行求解分析.结果表明,在等离子体鞘层区域电子能量具有Maxwell分布形式,在正柱区域具有Druyvesteyn分布形式.最高电子能量和最大反应率系数出现在鞘层区域.电子的迁移率系数和扩散率系数随射频周期的演化时空分布不均匀.     

正极性纳秒脉冲电压下变压器油中流注放电仿真研究

建立了二维轴对称流体模型, 仿真研究了正极性纳秒脉冲电压下变压器油中针-板电极流注放电的起始与发展过程, 得到了不同的外施电压幅值、脉冲上升沿时间与电极间隙距离下油中流注放电的形貌、 电场强度与空间电荷密度分布等. 仿真结果表明: 空间电荷加强了流注头部前方电场, 使流注通道更易于向前推进, 形成"电离波"; 随着外施电压幅值升高, 流注发展的平均速度显著变大; 较陡的脉冲上升沿形成的放电半径较大, 对应的最大电场强度值变小; 随着电极间隙距离的增加, 流注发展平均速度变快. 仿真显示纳秒脉冲下放电中油温无明显升高, 表明此类放电过程没有明显的油气化现象. 我们认为, 场致电离是油中带电粒子产生的主导机制; 空间电荷效应增强流注前方电场使得电离进一步发展, 最终导致击穿. 本研究有助于加深对变压器油中放电起始、发展直至击穿过程的认识以及对液体电介质中电离机制的理解. 关键词： 变压器油 流体模型 流注放电 空间电荷效应     

常压下气体放电等离子体振荡的实验与理论研究

在实验上研究了高压交流电弧发生器电极间隙的气体放电及等离子体振荡, 观察到了气体放电过程中的纳秒脉冲.以电子的流体运动方程和麦克斯韦方程为理论基础, 利用δ函数来描述交变外电场作用下电极处的电子堆积现象,建立了常压下气体放电时等离子体在外电场中振荡的理论模型,通过Laplace变换求解出电极间的放电电压.理论与实验结果基本符合, 从而可估算出实验中等离子体的电子数密度为1.3× 10¹²/m³.     

光电离速率影响大气压空气正流注分支的机理研究

大气压空气中的流注放电有广泛的理论和应用研究价值,包括雷电机理、输变电系统空气绝缘理论以及材料表面改性等.流注是一个快速发展的强电离区域,在传播过程中存在着一种重要的特点—分支现象.光电离为正流注发展提供必要的自由电子,且实验结果表明分支特征与流注头部的光电离速率密切相关.本文基于新的流注分支判据,采用了粒子网格单元与蒙特卡罗碰撞相结合(PIC-MCC)的三维放电模型(Pamdi3D)进行数值仿真验证.为了研究光电离速率对正流注分支的影响,仿真了毫米尺度间隙针-板电极正流注发展,系统研究了不同光电离参数的影响.当减小氮气-氧气比例、光子吸收截面或光电离效率系数后,流注均更早地出现分支现象.这些计算结果表明大气压空气中流注头部光电离速率的降低将导致其发生分支的概率更高.     

脉冲电磁驱动低温螺旋流注放电机理

特定脉冲放电参数条件下,原来沿直线传播的低温等离子体射流放电通道会在脉冲电磁驱动下发生放电模式转换,形成螺旋形态的流注通道.与传统的螺旋波放电不同,实验中并没有外加恒定磁场等破坏介质管极向对称性的因素,而且螺旋的手性方向会随放电参数变化.为深入了解等离子体射流中螺旋结构的电磁作用机理,探索其导致螺旋形态和决定手征性的极向电场来源及影响因素,本文通过建立完备自洽的等离子体理论模型,对螺旋流注的手性方向、螺距变化、分叉等复杂放电特性和电磁机理进行分析.研究发现击穿时极向电场波模的初始相位对螺旋流注的手性选择存在影响,电子密度对螺旋流注的螺距参数存在影响,放电的脉冲重复频率对螺旋流注分叉的分叉点位置存在影响.以上放电特性及其影响因素对探明电磁波与等离子体的相互作用机制具有重要科学意义,同时也为低温等离子体的手性应用提供实验和理论支撑.     

电子散射和能量分配方式对电子输运系数的影响

电子输运系数是确保低温等离子体建模准确性的关键因素,通过模拟电子的输运过程可对其数值求解.在模拟电子输运时,电子和中性粒子碰撞后的散射和能量分配方式有多种处理方法.为了研究不同处理方法对电子输运系数的影响,本文基于蒙特卡罗碰撞方法,建立了电子输运系数的计算模型,模拟约化电场10—1000 Td (1 Td=10~(–21) V·m~2)氢原子气中的电子输运过程.计算结果表明,各向同性假设对电子输运系数的影响随电场强度增加而增加,但即使对于较低的约化电场(10 Td),各向异性散射假设下电子的平均能量、通量迁移率和通量扩散系数也分别比各向同性假设下的值高38.34%, 17.38%和119.18%.不同的能量分配方式对中高电场强度下( 200 Td)的电子输运系数影响较为显著.在高电场时,均分法计算得出的电子平均能量、通量迁移率和通量扩散系数均小于零分法对应的值,汤森电离系数则相反. Opal法得出的电子输运系数介于均分法和零分法之间.此外,考虑各向异性散射时,不同能量分配方式对输运系数的影响高于各向同性.本研究表明,在计算电子输运系数时需要考虑各向异性的电子散射,高电场条件下尤其要注意能量分配方式的选择.     

超热电子在稠密等离子体中输运的混合粒子模拟方法

相对论激光等离子体相互作用以及所产生的带电粒子束在高密度等离子体中的输运行为非常重要.该物理问题的数值模拟研究仍面临技术挑战.本文介绍一种粒子/流体混合模拟方法.该方法中超热电子采用动力学方法描述,背景冷的稠密等离子体采用简化的流体方程描述,适合于超热电子密度远小于背景电子密度,超热电子能量远大于背景电子温度.我们的三维并行混合模拟程序HEETS的模拟结果表明：背景材料的电离和电阻率模型至关重要,将严重影响高能电子输运过程的模拟.     

容性耦合射频氩等离子体放电诊断研究及仿真模拟

针对中等气压、中等功率下射频容性耦合(CCRF)等离子体的放电特性,采用基于流体模型的COMSOL软件仿真,建立一维等离子体放电模型,以Ar为工作气体,研究同一气压时不同射频输入功率下等离子体电子温度和电子密度的分布规律。同时依据仿真模型设计制作相同尺寸的密闭玻璃腔体和平板电极,实验测量了不同射频输入功率时放电等离子体的有效电流电压及发射光谱,进而计算等离子体的电子温度及电子密度;利用玻耳兹曼双线测温法,得到光谱法下等离子体的电子温度及电子密度。结果表明：当气体压强为250 Pa、输入功率为100~450 W时,等离子体电压电流呈线性关系,电子密度随功率的增大而增大,而电子温度并未随功率的变化而有明显变化,其与功率无关。运用仿真模拟验证了实验的准确性,通过比较,三种方法所得的结果相近。通过结合等效回路法、光谱法和数值模拟仿真法初步诊断出中等气压下等离子体的放电参数,提出了结合三种方法作为实验研究的方法,使实验结果更具说服力,证明其方法的可靠性,也为进一步的等离子体特性研究提供依据。     

Streamer discharge simulation in flue gas

Streamer discharges are formed in a dielectric-barrier discharge used for nonthermal plasma generation. The results of simulation of streamer type discharge in a flue gas mixture is reported. A Monte Carlo simulation is done to obtain the transport and appropriate rate coefficients. The transport and rate coefficients calculated from the Monte Carlo is used to solve the conservation equations for electron, positive and negative ions, together with the Poisson's equation. The G-factor (radicals produced per 100 eV of electrical energy input to the discharge) obtained for Townsend-type discharge is higher as compared to a streamer-type discharge. Also experimental results of the SO₂ removal efficiency is compared to theoretical predictions     

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.     

Benchmark calculations for Monte Carlo simulations of electrontransport

Benchmark calculations have been performed for electron transport coefficients with an aim to produce a body of data required to verify the codes used in plasma modeling. The present code for the time resolved Monte Carlo simulation (MCS) was shown to represent properly DC transport coefficients in a purely electric field, in crossed electric and magnetic fields, and in the presence of nonconservative collisions, ionization, and attachment. In addition, we have suggested tests of the time dependent solutions. Relaxation of the initial transport coefficient may serve as an accurate test of the code as well as the input data for some fluid codes. In this paper, we show only one example, but several different sets of conditions and cross sections should be used as well. Finally, we propose application of the quasi-steady state results in RF fields. As an example we suggest calculation of the components of diffusion tensor showing anomalous longitudinal diffusion and calculations made with nonconservative collisions (ionization in this case). We also check the application of approximate formulas to determine drift velocity on the basis of total collision frequency and to determine a diffusion coefficient by using the Einstein relation. Other tests required to verify the transport data calculations are discussed as well     

Transport Coefficients for Inelastic Maxwell Mixtures

The Boltzmann equation for inelastic Maxwell models (IMM) is used to determine the Navier–Stokes transport coefficients of a granular binary mixture in d-dimensions. The Chapman–Enskog method is applied to solve the Boltzmann equation for states near the (local) homogeneous cooling state. The mass, heat, and momentum fluxes are obtained to first order in the spatial gradients of the hydrodynamic fields, and the corresponding transport coefficients are identified. There are seven relevant transport coefficients: the mutual diffusion, the pressure diffusion, the thermal diffusion, the shear viscosity, the Dufour coefficient, the pressure energy coefficient, and the thermal conductivity. All these coefficients are exactly obtained in terms of the coefficients of restitution and the ratios of mass, concentration, and particle sizes. The results are compared with known transport coefficients of inelastic hard spheres (IHS) obtained analytically in the leading Sonine approximation and by means of Monte Carlo simulations. The comparison shows a reasonably good agreement between both interaction models for not too strong dissipation, especially in the case of the transport coefficients associated with the mass flux     

The nonequilibrium behavior of electron swarms in nonuniform fieldsin SF6

The behavior of electron swarms in nonuniform electric fields in SF₆ has been analyzed by Monte Carlo simulation. Decreasing and increasing electric field configurations have been studied. The spatial variations of the electron mean energy, drift velocity, and ionization and attachment coefficients are calculated for different field slopes and compared with the values calculated assuming equilibrium field conditions. It is found that the nonequilibrium behavior depends on the ratio of the field slope to pressure, and not on the field slope alone. The mean energy and drift velocity are almost the same as the equilibrium values, while the ionization and attachment coefficients show significant nonequilibrium features. The electron and ion distributions along the uniformly decreasing and increasing fields have also been studied     

低能电子在固体表面背散射系数的直接Monte Carlo方法模拟

应用单次碰撞的直接Monte Carlo方法计算能量范围从100 eV~10 keV的电子在固体Al,Si,Au表面的背散射系数,其中低能电子在固体中的弹性散射和非弹性散射截面分别应用Mott散射截面和Born近似下的广义振子强度计算模型得到.通过与压缩历史Monte Carlo方法的模拟计算结果及实验值的比较,结果表明,对于100 eV~10 keV范围的低能区电子,采用直接方法计算得到的电子背散射系数与实验值符合较好,直接方法比压缩历史方法更适合于能量在10 keV以下的电子输运计算.     

One-Particle Distribution Functions and Thermodynamics of Crystals with Many-Body Forces. II. Quantum Corrections

Crystalline state is described by non-symmetrical distribution functions in a self-consistent field approximation. The first quantum corrections to the HELMHOLTZ free energy and thermodynamic properties of strongly anharmonic crystals with many-body interactions are calculated in the quasi-classical case. The properties of crystalline inert gases are computed using the pair potentials of LENNARD -JONES , of BARKER and POMPE , of BARKER and BOBETIC , and of BARKER , FISHER and WATTS , together with the AXILROD -TELLER -MUTO three-body forces. The results are compared with both the experimental data and Monte Carlo calculations.     

Vorstadium des Funkens,untersucht mit dem Bildverstärker

Streamer development (near cathode, midgap, and near anode streamers) was investigated by means of image intensification and image deflection as to conditions for streamer onset and streamer velocities during the various stages of the avalanche-to-plasma channel transition in several gases, e.g. N₂, H₂, A, CH₄ as well as air, and some more mixtures of gases, in the pressure range 100 to 500 Torr. Quantitative results were obtained, and the influence of the discharge parameters upon streamer development is discussed. — The enhanced and accelerated streamer stage (II) was now observed in either direction; the significantly fast and steep ionization waves (stage III) were found to start at the cathode as well as at the anode on arrival of the respective streamers. Velocities up to 10⁹ cm/sec were measured.     

Backscattering coefficients and mean penetration depths of 1–4 keV electron scattering in solids

In the present study, we have performed Monte Carlo simulation of 1–4 keV electrons impinging on semi-infinite Al and Au to determine the transport cross-section, the backscattering coefficient and the mean penetration depth using a new approximation of the differential elastic scattering cross section. The mean number of the wide angle collisions suffered by the electron before slowing down to rest, and the backscattering coefficient are analytically calculated using Vicanek and Urbassek theory. The analytical results are compared with the numerical ones obtained from Monte Carlo simulation. The present results are found to reveal good agreement with experimental results. PACS 68.37.-d; 68.37.Nq