基于负氢离子源的全三维PIC/MCC模拟算法研究

在分析负氢离子源中等离子体物理机理基础下, 研究并优化粒子模拟算法, 设计高效的粒子存储方法. 研究并运用粒子碰撞蒙特卡罗方法, 考虑等离子体势以及带电粒子间库仑碰撞, 研制了全三维粒子模拟/蒙特卡罗算法(PIC/MCC). 采用磁荷模型, 运用时域有限差分方法计算多峰磁场, 并结合国外负氢离子源JT-60U, 考虑负氢离子源中主要反应, 对全三维PIC/MCC模拟算法模拟验证.     

潘宁源放电的全三维电磁粒子模拟/蒙特卡罗碰撞数值算法研究

深入研究潘宁放电的物理机制, 研制了全三维高品质算法粒子模拟软件(PIC), 设计并添加了相应物理情景的蒙特卡罗碰撞模块(MCC), 并对电子、氢分子离子(H₂⁺)、氢正离子(H⁺)、氢三正离子(H₃⁺)同时进行了跟踪, 成功研制了全三维电磁PIC/MCC数值算法. 结合国内研究较热的潘宁放电模型, 对该算法进行模拟验证. 模拟结果显示: 采用有效的滤波算法能抑制电磁数值噪声, 电子能量呈麦克斯韦分布, 由于电子的径向漂移和加速导致离子源顶端H₂⁺产量较大. 关键词： 潘宁离子源 高品质算法 粒子模拟/蒙特卡罗     

气体电离的全三维电磁粒子模拟/蒙特卡罗数值研究

深入研究电子束对中性气体电离的物理机理,在粒子模拟(PIC)软件CHIPIC的基础上,设计了蒙特卡罗(MCC)电离碰撞模块,并对电子与离子同时进行了跟踪,成功研制了全三维电磁PIC/MCC代码.通过该软件对填充氦气相对论返波管的模拟,对全三维电磁PIC/MCC代码进行验证.模拟结果显示:填充气体可以中和空间电荷限制效应,有效提高电流大小;填充适量的气体可提高功率峰值,扩展脉冲宽度;过量的气体则会缩短脉冲宽度、降低峰值功率.     

蒙特卡罗碰撞算法在CHIPIC软件中的数值实现

分析离子源中电子所参与的碰撞,利用电子与其它粒子间的空碰撞模型,研究电子之间的库仑碰撞,采用不同的存储和调用机制研制了两种碰撞模型的蒙特卡罗碰撞处理模块,将数值模块添加到粒子模拟软件CHIPIC中,对离子源JAERI 10A进行模拟验证,通过对模拟结果的分析,表明所设计的三维MCC算法正确.     

电子回旋共振放电的电离特性PIC/MCC模拟(Ⅰ)——物理模型与理论方法

采用粒子模拟与蒙特卡罗相结合(PIC/MCC)的方法对电子回旋共振(ECR)放电中的电离过程进行了模拟，其中带电粒子与微波的相互作用由PIC方法的电磁模型描述，粒子间的碰撞过程由MCC方法描述.考虑的碰撞类型有电子与中性粒子的弹性、激发、电离碰撞，离子与中性粒子的弹性、电荷交换碰撞，碰撞截面均依赖于能量而变化.阐述了理论分析的过程，为数值模拟ECR放电奠定了基础. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗 电离     

电子回旋共振放电中电子与微波互作用特性的粒子模拟和蒙特卡罗碰撞模拟

对电子回旋共振（ECR）放电电离过程中的电子与微波互作用特性进行了理论分析与数值模拟.采用粒子模拟（PIC）方法描述带电粒子与微波的互作用,采用蒙特卡罗碰撞（MCC）方法描述粒子间碰撞过程及带电粒子与边界的相互作用.编写了准三维的PIC/MCC数值模拟程序,并对放电过程中电子能量与微波场随时间、空间的演化进行了数值诊断. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗方法 电离     

15～20 mA多峰负氢离子源全三维数值模拟

对自主开发的全三维粒子模拟/蒙特卡罗(PIC/MCC)算法进行了阐述,包括算法流程设计、碰撞处理机制等。采用该算法对中国原子能科学研究院研制的15~20mA负氢离子源进行数值模拟研究,结果显示:该多峰会切离子源能得到空间分布均匀的离子束,采用优化后的虚拟过滤场,不仅能有效过滤高能电子,还将导致离子源下游区域低能电子增多,有效提高负氢离子体积产生效率。     

电子回旋共振放电的电离特性PIC/MCC模拟(Ⅱ)——数值模拟与结果讨论

采用粒子模拟与蒙特卡罗相结合(PIC/MCC)的方法，应用电磁模型，编写了准三维的电子回旋共振(ECR)放电电离过程的模拟程序，得到了ECR放电过程中电子与离子的相空间分布、电磁场分布.通过对这些分布随时间演化的分析，得出ECR加热发生在ω≈ω_c0且垂直于轴向的区域；ECR区域，微波能量几乎全部耦合给电子，获得能量的电子通过与中性粒子的电离碰撞产生了大量的带电粒子；随着放电的进行，大量带电粒子通过频繁的碰撞，分布由各向异性逐渐趋于各向同性. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗 电离     

负氢离子源中电子能量沉积三维数值模拟研究

本文理论计算了多峰离子源永磁体, 采用二体碰撞模型(binary collision)处理电子之间的库仑碰撞, 采用空碰撞(null-collision)方法处理电子与氢元素相关粒子, 开发了全三维PIC-MCC模拟算法, 并利用该算法模拟了多峰离子源中电子沉积过程, 分析了多峰磁场对电子的空间和能量分布影响, 结果显示: 电子的非均匀分布起源于高能电子在引出区的B× ▽B漂移.     

库仑碰撞截面在等离子体粒子模拟中的应用

在等离子体粒子模拟中,TA模型和NanBu模型被广泛用于处理库仑碰撞,这两种模型要求每个时间步长内全部粒子参与计算。为了降低参与碰撞的粒子数,提高库仑碰撞的计算效率,提出了一种基于截面的库仑碰撞模拟方法,并给出了库仑碰撞概率的计算公式。采用该方法对不同温度不同密度电子气的弛豫过程进行模拟,分别对比了电子速度分布函数、电子温度以及电子x、y方向上的温度与电子温度之比的模拟值与理论值,验证了该方法的准确性。在相同的小时间步长上,该方法相比TA模型计算效率提升可达40%以上。对于较大的时间步长,该方法仍能得到与理论解近似的模拟结果,相比Nanbu模型,在相同的精度下可取更大的时间步长,计算效率也有所提升。研究表明,该方法同样适用于电子-离子碰撞。因此在提高库仑碰撞计算效率上,该方法具有碰撞粒子数少以及适用于大时间步长的优势。     

Particle simulation of Coulomb collisions: Comparing the methods of Takizuka & Abe and Nanbu

The interactions of charged particles in a plasma are governed by long-range Coulomb collision. We compare two widely used Monte Carlo models for Coulomb collisions. One was developed by Takizuka and Abe in 1977, the other was developed by Nanbu in 1997. We perform deterministic and statistical error analysis with respect to particle number and time step. The two models produce similar stochastic errors, but Nanbu’s model gives smaller time step errors. Error comparisons between these two methods are presented.     

Self-consistent modelling of X-ray preionized XeCl-laser discharges

In this theoretical work a 0-D model for a self-sustained X-ray preionized XeCl-laser discharge is presented. The model is self-consistent in the sense that it simultaneously solves, contrarily to the usual decoupling procedure, the Boltzmann equation for electrons, the kinetic equations for excited and ionic species, the equations for the electrical circuit and the laser photon density. It includes a rather complete kinetics of HCl(v) vibrational excitation, dissociation and dissociative attachment. The influence of electron collisions with excited species and of e-e Coulomb collisions on the plasma parameters and transport coefficients is discussed. Some evidence of the non-stationary equilibrium between the electron distribution and the reduced electric field E/N is given. Results of the model are compared with experimental ones corresponding to a XeCl-laser discharge driven by a L-C inversion circuit. The model predicts well the main trends for the variation of the laser energy in a large range of experimental conditions. The discrepancy between experiment and model for absolute values of the laser energy is discussed.     

Understanding the accuracy of Nanbu’s numerical Coulomb collision operator

We investigate the accuracy of and assumptions underlying the numerical binary Monte Carlo collision operator due to Nanbu [K. Nanbu, Phys. Rev. E 55 (1997) 4642]. The numerical experiments that resulted in the parameterization of the collision kernel used in Nanbu’s operator are argued to be an approximate realization of the Coulomb–Lorentz pitch-angle scattering process, for which an analytical solution for the collision kernel is available. It is demonstrated empirically that Nanbu’s collision operator quite accurately recovers the effects of Coulomb–Lorentz pitch-angle collisions, or processes that approximate these (such interspecies Coulomb collisions with very small mass ratio) even for very large values of the collisional time step. An investigation of the analytical solution shows that Nanbu’s parameterized kernel is highly accurate for small values of the normalized collision time step, but loses some of its accuracy for larger values of the time step. Careful numerical and analytical investigations are presented, which show that the time dependence of the relaxation of a temperature anisotropy by Coulomb–Lorentz collisions has a richer structure than previously thought, and is not accurately represented by an exponential decay with a single decay rate. Finally, a practical collision algorithm is proposed that for small-mass-ratio interspecies Coulomb collisions improves on the accuracy of Nanbu’s algorithm.     

Untersuchung der Prozesse in der aktiven Zone der Hohlkatodenbogenentladung I. Die Energieverteilungsfunktion der Elektronen,die Anregungs- und Ionisationsgeschwindigkeit

In this paper the processes of excitation and ionization inside the cathode of a hollow cathode are discharge will be studied. The electron energy distribution function is calculated from the kinetic equation. For this we take into account the following processes: elastic and inelastic collisions of electrons with atoms, Coulomb interaction between emitted fast electrons and electrons of the plasma.     

On the Kinetics of the Active Zone of the Stationary SF6 Beam Discharge Plasma

The paper deals with the impact of intensive electron attachment on the kinetics of the electrons in the active zone of the stationary band-like beam discharge plasma in SF₆ which is an alternative useful plasma medium for “dry etching”. The energy distribution of the electrons in this plasma was obtained by numerically solving the Boltzmann equation which includes apart from elastic collisions, different exciting collision processes, attachment in electron collisions, direct ionization, the ambipolar loss of electrons, Coulomb interaction between electrons and of electrons with ions and the power input to the electrons by the turbulent electric field. In particular, due to the needed fulfilment of the consistent electron particle balance, for an extended region of the turbulence energy density in this plasma a large impact on the electron kinetics of the intensive electron attachment, which is the prevailing electron loss process, was found enforcing independent of the turbulence energy density always a large power input to the electrons, smooth and only slowly decreasing energy distributions even in the energy region of direct ionization.     

Ion-induced auger electron emission of Mg,Al and Si as a function of ion energy

Experiments have been performed in which ion-induced Auger electron yields of Mg, Al and Si were measured as a function of ion energy. A computer simulation based on a binary collision cascade model yielded results that were in close agreement with the experimental findings. From this a model is proposed, where the creation of inner shell vacancies occurs during symmetrical collisions in the collision cascade generated by the ion. From the minimum energy required for Auger emission a distance of closest approach between the particles can be derived, using the Molière approximation of the screened Coulomb potential to describe the collisions. This distance can be correlated with the radii of the interacting orbitals derived from Hartree-Fock-Slater calculations.     

直流空心阴极放电中鞘层区电子的输运过程

采用蒙特－卡洛方法对空心阴极放电中鞘层区电子的输运过程进行了研究，电子在鞘层区被非均匀电场加速，两次碰撞之间的步长是由电子与中性粒子的碰撞频率确定。模型中三种碰撞截面积是由实验和理论数据拟合而来。研究了电子平均能量、电子密度和电离系数在径向的分布。电子能量的空间分布结果与实验很好地符合。 关键词：     

高振动激发LiCs与CO_2碰撞中的振动-转动能量转移

受激发射泵浦得到的LiCs高位振动态与CO_2碰撞,研究了LiCs不同能量对碰撞能量转移过程的影响.利用高分辨率瞬时激光感应荧光(LIF)测量得到CO_2(0000,J=2-74)原生转动态分布.转动布居数的半对数描绘给出了一个双指数结构,得到了一个T_a=660±73K[对LiCs(E=4300cm~(-1))]和550±61K[对LiCs(E=3700cm~(-1))]的低能分布,一个T_b=2380±261K[对LiCs(E=4300cm~(-1))]和1980±217K[对LiCs(E=3700cm~(-1))]的高能分布,低能分布属于弹性或弱非弹性碰撞,高能分布属于强非弹性碰撞.转动分布对LiCs能量是敏感的.但弹性与非弹性分支比基本是相同的.在一次碰撞的条件下,测量了各J态的出现和倒空速率系数,它们仅弱依赖于LiCs激发能.出现和倒空速率系数对不同LiCs能量基本一致,是振动-转动/平移弛豫轨道的重要证据.