磁绝缘线振荡器内气体碰撞电离数值模拟.pdfVIP

第24卷第3期 强 激 光 与 粒 子 束 V01．24。No．3 I。ASERANDPARTICLEBEAMS Mar．。2012 2012年3月 HlGHPOWER 文章编号： 1001—4322(2012)03—0732—07 磁绝缘线振荡器内气体碰撞电离数值模拟。 董 烨， 董志伟， 杨温渊， 陈 军， 周前红， 孙会芳， 周海京 (北京应用物理与计算数学研究所，北京100094) 摘要：针对气体碰撞电离过程，介绍了蒙特卡罗碰撞(MMC)的处理方法，利用MMC方法编写了气体 碰撞电离模块。将其移植到3维全电磁粒子模拟程序NEPTUNE之中，模拟了充有He气的磁绝缘线振荡器 (MILO)。模拟结果表明：当He气密度较低时，电离的正离子由于较重无法自由移动。形成了正离子通道，可 以有效中和电子束空间电荷场，有利于电子束传输和群聚，提高了柬波互作用效率，微波输出功率得到了明显 提高，起振时间也有所缩短；当进一步增加He气密度时，电离碰撞增强，电子和离子数目会雪崩式增长，电子 束由于碰撞增强而导致能散度增大，其负效应已经远大于中和空间电荷场的正效应。反而不利于电子柬的群聚 和共振，从而导致输出微波功率降低乃至截断，起振时间缩短是由于其在非雪崩阶段的正效应积累所致，但是 随着负效应的增强起振功率不能得以维持，二极管最终将闭合。另外，还模拟了MILO填充空气、水蒸气及二 氧化碳等多原子、多组分气体的碰撞电离物理过程。模拟结果显示，同压强情况下，填充空气、水蒸气及二氧化 碳的脉冲缩短现象要比填充He气等较低原子序数气体的情况严重得多。 关键词： 高功率微波；蒙特卡罗碰撞；粒子模拟程序；磁绝缘线振荡器；脉冲缩短 文献标志码： A 中图分类号：TN814，0411．3 射频结构电击穿[1]是高功率微波(HPM)产生和传输的根本限制因素，也是输出微波脉冲缩短的直接诱 因。初级或次级电子与中性气体分子发生碰撞电离形成等离子体，其密度到达一定程度便会击穿，破坏柬波互 作用机制，致使输出微波脉冲缩短。因此，研究气体碰撞电离对认识微波脉冲缩短形成机制具有重要意义。目 的MCC—PIC描述方法，有关结论尚待进一步考证。为满足我国高功率微波发展需求，我们研发了并行3维全 此模块对磁绝缘线振荡器(MIL0)内气体碰撞电离进行数值模拟研究。 1 MCC-PIC程序求解过程及气体碰撞类型简介 MCC—PIC程序的计算流程如图1所示【l喝]。 首先，将位于网格棱心和面心位置的电场与磁场 分量利用PIC插值方法权重到粒子位置计算出粒 子受力；而后根据粒子受力和初始位置，利用 Newton-Lorentz方程推进粒子，得到粒子速度和 新位置；接着判断粒子是否到达边界被吸收或是 否发射新粒子；若存在粒子碰撞，则进入MCC过 程，判断粒子是否产生某类碰撞，若产生则计算其 Fig．1 flowchartofMCC—PICcode 碰撞后的速度；而后，根据PIC方法插值权重各个 Computational 图l MCC-PIC程序计算流程图 粒子到网格结点和棱心位置，计算电荷密度和电 *收稿日期：2011-10—17；修订日期：2011-11-07 展计划项目 作者简介：董烨(198l一)。男。助理研究员，从事高功率微波技术方面的研究；dongye0682@sina．corn。 万方数据 第3期 董烨等：磁绝缘线振荡器内气体碰撞电离数值模拟 733 collisionbet、v嘲electronandatomcollisionbet、Ⅳ唧ionandatom 时间步的计算。如图2所示，MCC的过