离子束从真空弧等离子体中引出的三维数值模拟.pdfVIP

离子束从真空弧等离子体中引出的三维数值模拟 蓝朝晖, 彭宇飞，杨 振，龙继东，董 攀，魏涛，李杰，王韬，王小胡，石金水 （中国工程物理研究院流体物理研究所， 四川 绵阳 621900 ） 摘 要：开发了无碰撞的三维粒子模拟程序用于模拟离子流从真空弧等离子体中引出，模拟的等离子体密度范围是 1016 m-3 ~1018 -3 m 。建立了含有全尺寸矩形栅网的等离子体扩张杯结构以及引出结构模型，对多种成分离子从扩张杯引出进行了数值模拟，程序 可以计算多种成分离子在 100kV 高压的引出全过程。数值模拟给出了不同等离子体密度情况下的空间电场分布，离子束发射面， 空间电荷分布，并且和无栅网情况下做了比较研究。数值模拟结果可用于引出结构以及栅网的优化设计。 关键词：真空弧等离子体；离子引出；粒子模拟；栅网 中图分类号：O463+.2 文献标志码：A 真空弧离子源是一种被深入研究的离子源，可以产生多种成分金属和非金属离子，广泛应用于离子注入，离 子束处理，重离子聚变以及离子加速器。这种类型离子源经过脉冲高压引出，可以产生强流脉冲离子束。基于强 流脉冲离子束（典型束流强度百 mA~A 量级，脉冲时间 μs 量级）的特种脉冲中子源（真空器件型）用途广泛， 例如可用于中子照相、爆炸物探测、聚变堆材料辐照损伤研究、半导体抗辐射加固和核裂变点火等工农业和国防 [1-2] [3] 研究领域 。国内外在器件整体及内部部件方面都开展了大量的研究工作，如美国的圣地亚实验室 、俄罗斯 [4] [5] [6] 的大电流所 、以色列 和中国工程物理研究院 等。 在特种脉冲中子源的研究中，离子光学系统的优化设计是一项重要内容。然而，由于探测手段、空间尺寸以 及高压环境的限制，实验研究不仅成本高、周期长，而且通常会面临很大的技术困难。数值模拟可以避免实验研 究的各种限制，极大地减小研究成本，提高研究效率。本文将采用自编的三维粒子模拟程序来研究离子的引出、 输运问题。首先建立三维的物理模型，然后计算不同密度等离子体情况下的离子引出，并对发射面的影响因素， 栅网的作用等进行分析和讨论 2 计算方法及物理模型 2.1 计算方法 粒子模拟方法是通过跟踪大量带电粒子在外加及自洽电磁场中的运动并统计平均而得到宏观特性及运动规 律的一种方法[7] 。离子引出计算可忽略瞬态磁场的影响，但是需要考虑空间电荷效应和相对论相应，因此使用静 电粒子模拟是合适的。整个控制方程如下： 2 d    m v q E q v B i i i dt （1） d 1 x =v E   2 dt 1 v c   以上方程按照顺序依次为泊松方程，带电粒子运动方程，带电粒子位移方程，电场计算方程以及相对论因子 计算公式。E 和 B 是电场和外磁场，电场根据电荷分布以及外加高压自洽求解，在程序中外磁场是恒定的，本文 设为 B=0 。 三维静电粒子模拟的计算步骤如下，首先给定一个初始密度和温度的均匀等离子体分布，接着求解泊松方程， 用超松弛迭代法计算整个计算空间网格节点的场分布。通过节点插值求得每个粒子所在