PICPMCC 数值模拟四阳极直流辉光放电

利用OOPIC 代码并结合MCC 方法考察四阳极低气压直流辉光放电中粒子间的碰撞, 并对其进行了数值模拟。首先分析OOPIC 代码的计算流程, 给出模拟的边界条件和时间、空间步进; 然后对所得到的等离子体密度、温度的分布情况进行了分析。结果表明, 放电气压、电流以及环电极的有效宽度均影响等离子体密度和温度的径向分布。当气压为300Pa 时, 电子密度在径向较大范围内是均匀分布的, 而当气压升高到800Pa 时, 其径向不均匀增大, 且总的密度和电子的平均温度下降。放电电流的增大使得等离子体密度和温度均相应升高。但环电极有效宽度的减小使得等离子体密度的分布更不均匀, 并使总的密度和平均温度均下降。     

电子回旋共振放电电离特性的PIC/MCC模拟

A theoretical and computational model is presented to study the ionization of the argon electron cyclotron resonance(ECR)microwave discharge using a quasi-three-dimensional electromagnetic particle-in- cell plus Monte Carlo collision method.The interaction between the charged particles and microwave fields are described by the electromagnetic mode of particle-in-cell method.The collision processes are treated with Monte Carlo method.The simulation code is the original work.The results of the particle simulation for the ECR discharge of argon gas which include the microscopic features of charged particles and the electromagnetic characteristics of the ECR discharge plasma,and also the transient phenomena have been presented.     

Particle‐in‐cell/Monte Carlo simulation of electron and ion currents to cylindrical Langmuir probe

Electron and ion currents to a cylindrical Langmuir (electrostatic) probe were calculated using the particle‐in‐cell/Monte Carlo (PIC/MC) self‐consistent simulation for a neutral gas in the pressure range 2–3,000 Pa. The simulation enables us to calculate the probe currents even at high neutral gas pressures when the collisions of collected charged particles with neutral gas particles near the probe are important. The main aim of this paper is the calculation of probe currents at such high gas pressures and the comparison of the results with experimentally measured probe currents. Simulations were performed for two cases: (a) probes with varying radii in a non‐thermal plasma with high electron temperature at low neutral gas pressure of 2 Pa (in order to verify the correctness of our simulations), and (b) probe with the radius of 10 μm in the afterglow plasma with low electron temperature and a higher neutral gas pressure (up to 3,000 Pa). The electron probe currents obtained in case (a) show good agreement with those predicted by the orbital motion limited current (OMLC) theory for probes with radii up to 100 μm for the given plasma conditions. At larger probe radii and/or at higher probe voltages, the OMLC theory incorrectly predicts too high an electron probe current for the plasma parameters studied. Additionally, a formula describing the spatial dependence of the electron density in the presheath in the collisionless case is derived. The simulation at higher neutral gas pressures, i.e. case (b), shows a decrease of the electron probe current with increasing gas pressure and the creation of a large presheath around the probe. The simulated electron probe currents are compared with those of measurements by other authors, and the differences are discussed.     

Electron Beam Optics System PIC Simulation for MMW TWTs

The basic feature and design method of MMW TWTs electron optics are analyzed in this paper. With the increasing area convergence ratio of the electron gun in TWTs, the initial thermal velocity of electrons at cathode and the non-magnetic shielded effect can no longer be neglected. It brings new challenge for the design and simulation of MMW TWTs electron optics. A method of investigation the beam DC characteristics is described using particle-in-cell (PIC) simulation and incorporate electron beam optics model integrated electron gun with periodic permanent magnetic focusing system. It is valuable for achieving the goal of first-pass design success and the electron optics engineering design and optimization.     

Potential Formation in Front of an Electrode Close to the Plasma Potential Studied by PIC Simulation

Potential formation in front of an electrode that has the potential which is close to the plasma potential is studied by particle‐in‐cell (PIC) simulations. The code BIT1 [D. Tskhakaya, R. Schneider, J. Comp. Phys., 225 , 829 (2007)] is used for this purpose. It is shown that this code is very appropriate for our analysis because of its ability to create charged particles by uniform volume production in the entire system and to maintain in the same time the Maxwellian electron velocity distribution function prescribed in the input file. It turns out that some modifications of the code are necessary in order to achieve small sheath potential drops and to detect the cutoff in the electron velocity distribution function. The modifications of the code are described. The small sheath potential drop is achieved by introducing a finite reflectivity of the electrons at the left electrode, while the expected cutoffs of the electron velocity distribution function are found by tracking the reflected electrons as separate particle species. The simulation results are compared to the theoretical model of Jeli? [N. Jeli?, Phys. Plasmas, 18 , 113504 (2011)]. The matching is very good and this is a sign that the PIC simulations are very appropriate tool for the analysis of this type of problems. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

N2微空心阴极放电特性及其阴极溅射的PIC/MC模拟

采用两维PIC/MCC模型模拟了氮气微空心阴极放电以及轰击离子 (N₂⁺,N⁺) 的钛阴极溅射. 主要计算了氮气微空心阴极放电离子 (N₂⁺,N⁺) 及溅射原子Ti的行为分布, 并研究了溅射Ti 原子的热化过程. 结果表明: 在模拟条件下, 空心阴极效应是负辉区叠加的电子震荡; 在对应条件下, 微空心较传统空心放电两种离子 (N₂⁺,N⁺) 密度均大两个量级, 两种离子的平均能量的分布及大小几乎相同; 在放电空间N⁺的密度约为N₂⁺的1/6, 最大能量约大2倍; 在不同参数 (P, T, V)下, 轰击阴极内表面的氮离子(N₂⁺,N⁺)的密度近似均匀, 其平均能量几乎相等; 从阴极溅射出的Ti原子的初始平均能量约6.8 eV, 离开阴极约0.15 mm处几乎完全被热化. 模拟结果为N₂微空心阴极放电等离子体特性的认识提供了参考依据. 关键词： 微空心阴极放电 PIC/MC模拟 2等离子体')" href="#">N₂等离子体     

Three-dimensional PIC/MCC simulation of electron deposition in JAEA 10 A ion sources

A systematic research on the electron deposition process in the JAEA 10 A ion source is carried out by using a particle-in-cell/Monte Carlo collision simulation, which is based on a full three-dimensional self-developed code. Two parts are studied. One is the space and energy distribution of fast and slow electrons, the other is the vibration excitation collisions between electrons and hydrogen moleculars. The results show that the inhomogeneity of electrons comes from the Y direction drift of the fast electrons(Te 25 eV) due to the action of the magnetic fields. This drift also increases the number of vibration excitation collisions in the-Y direction, and results in the increase of Hain the-Y direction,eventually leading to the-Y drift of H-. It explains the spatial non-uniformity in the JAEA 10 A ion source.     

3 mm波段回旋速调管放大器设计与粒子模拟

设计了一支3 mm 波段基波回旋速调管,该回旋速调管工作在低损耗的TE01模式,包含四个谐振腔。首先使用线性理论确定工作参数的大致范围, 然后采用HFSS软件设计单个谐振腔,通过调整谐振腔尺寸和腔壁介质层参数使谐振腔的谐振频率和Q值符合设计要求, 最后使用粒子模拟程序优化设计了回旋速调管的互作用电路,研究了谐振腔参差调谐方案, Q值对回旋速调管性能的影响, 互作用电路的稳定性以及电子注参数变化对注-波互作用性能的影响。PIC粒子模拟结果表明,在电子注电压65 kV, 电流6 A, α(V⊥/V∥)1.5, 工作磁场3.6 T时,回旋速调管的3 dB带宽约为600 MHz,在93.7 GHz获得139 kW 的峰值输出功率,效率为35.6%,增益为28.4 dB。模拟中没有考虑电子注速度零散的影响。     

Particle-in-Cell/Monte Carlo Collision simulation of planar DC magnetron sputtering

In this paper a numerical simulation of a planar DC magnetron discharge is performed with the Particle-in Cell/Monte Carlo Collision （PIC/MCC） method. The magnetic field used in the simulation is calculated with finite element method according to experimental configuration. The simulation is carried out under the condition of gas pressure of 0.665 Pa and voltage magnitude of 400V. Typical results such as the potential distribution, charged particle densities, the discharge current density and ion flux onto the target are calculated. The erosion profile from the simulation is compared with the experimental data. The maximum erosion position corresponds to the place where the magnetic field lines are parallel to the target surface.     

1.2MeV脱硫脱硝电子加速器的设计与模拟

对工业用1.2MeV脱硫脱硝电子加速器进行了结构设计与性能分析，并利用PIC(Particle-In-Cell)模拟方法得到优化后的计算结果. 结果显示，加速器的电子光学性能有所提高，基本达到优化设计的目的.     

气体压强及表面等离激元影响表面波等离子体电离发展过程的粒子模拟

基于表面等离激元 (SPP) 的表面波等离子体 (SWP) 源, 具有高密度、低温度及高产率等优异性能, 其应用在电子器件微纳加工、材料改性等领域. 但由于SPP激励SWP放电的电离过程难于用理论分析和实验测量描述, 因而SWP源均匀稳定产生的电离发展过程一直未研究清晰. 本文以SWP放电的数值模拟为研究手段, 采用等离子体与电磁波相互作用的粒子模拟 (PIC) 方法, 结合蒙特卡罗碰撞 (MCC) 方法处理碰撞效应的优势, 研究气体压强影响电离过程的电磁能量耦合机理. 模拟结果表明SWP的高效产生是SPP的局域增强电场致使, 气体压强能够改变波模共振转换的出现时刻而影响了SWP的电离发展过程. 本文的研究成果展示了SPP维持SWP放电的电离过程, 可为下一代米级SWP源的参数优化提供设计建议. 关键词： 表面波等离子体 表面等离激元 粒子模拟 电离过程     

Optimization of Operating Parameters of Plasma Antenna

In this work, the operating parameters of the plasma antenna are optimized using a kinetic model based on Particle in Cell‐Monte Carlo Collisions (PIC‐MCC) method. This optimization study is performed via the investigation of variations in the operating parameters of the plasma antenna, i. e., its dimensions, background gas pressure, and the applied voltage frequency and their consequent effects on the plasma frequency, kinetic energy of electrons and plasma current density of plasma antenna. While the antenna performance is improved at higher tube lengths and applied frequencies, it is optimized at a particular tube radius. Moreover, higher background pressures have increasing effects on the plasma antenna operation. Based on this parametric study, the optimum operating parameters of the plasma antenna are established. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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

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

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.     

螺旋波等离子体推进器中双层形成的PIC模拟

建立了螺旋波等离子体推进器中双层等离子体形成过程的一维全粒子PIC模型,讨论了放电过程中电子加热模型和等离子体膨胀模型,计算了不同条件下螺旋波等离子体推进器中双层等离子体的分布。结果表明：采用氩气工质,在射频电流90 A/m2和工作压强0.05~0.15 Pa的条件下,可产生双层等离子体,其密度为1.31015~2.81015 m-3,双层电势降为10~30 V,可加速氩离子束流,喷射速度为3~12 km/s,双层加速效应对推进器的性能具有一定影响。     

Simulation of a Planar Emissive Probe in a Mid‐Sized Tokamak Plasma

Planar emissive probe is studied for the first time using a massively parallel particle‐in‐cell code BIT1 [22]. The probe is immersed in a plasma similar to edge plasmas of mid‐sized tokamaks. Dependence of the floating potential on electron emission from the probe is studied. With increasing emission the floating potential increases, but then saturates ～2T_e below the actual plasma potential (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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

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

W波段带状注扩展互作用振荡器的粒子模拟

设计了一种工作在W波段的振荡器。采用宽高比值为2的薄矩形电子注降低空间电荷效应,周期耦合腔慢波线作为高频结构以增加功率容量,阶梯渐变矩形波导作为输出结构,使输出结构与高频结构匹配连接减少反射。利用三维PIC粒子模拟程序进行数值模拟,分析慢波结构的周期数及周期长度对输出功率、频率及效率的影响。结果表明：当阴极电压为13.0~16.1 kV时,谐振腔能正常工作;谐振腔周期数为9,阴极电压为15.7 kV时输出功率和转换效率最大,输出峰值功率大于1.9 kW,平均功率为980 W,频率为91.6 GHz,转换效率达25%,输出频率随谐振腔周期长度的减小而增大。     

完全匹配层吸收边界在2D3V等离子体粒子模拟中的实现

基于Berenger完全匹配层的基本思想,导出了高斯单位制下2维完全匹配层吸收介质的控制方程及参数设置方法,并进行了数值验证和参数优化。依据所得结论编制了相应的吸收边界模块加入到2D3V等离子体粒子模拟程序PLASIM中,取得了良好的吸收效果。基于完全匹配层吸收边界的实现,在小窗口内对超强激光预脉冲与低密度等离子体的长时间相互作用过程进行了模拟研究,结果表明:激光强度和作用时间是影响等离子体密度改变的重要因素,超强激光其长时间的预脉冲可以通过有质动力的累积效应在低密度等离子体中引起较大的密度变化。     

基于Kr工质的离子推力器栅极PIC仿真

Kr是新型低成本离子推力器的优选工质之一。采用三维PIC粒子模拟方法对Kr等离子体在栅极系统中的输运过程进行仿真,并与Xe工质仿真结果进行对比,分析工质类型对束流分布、离子通过率、发散角损失和推力等参数特性和分布的影响。结果表明：相比于Xe工质,Kr工质条件下,屏栅上游鞘层区域更大,屏栅离子通过率略有上升;Kr+在加速栅孔轴向速度较大,径向位移小,撞击加速栅几率小,加速栅电流相对较小;Kr工质和Xe工质的发散角损失相当,但Kr工质工作电流区域更高;相比于Xe工质,相同功率条件下Kr工质的推力值约降低20%。