搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

收集极释气对相对论返波管影响的粒子模拟

王洪广 柳鹏飞 张建威 李永东 曹亦兵 孙钧

收集极释气对相对论返波管影响的粒子模拟

王洪广, 柳鹏飞, 张建威, 李永东, 曹亦兵, 孙钧
PDF
HTML
导出引用
  • 相对论返波管被认为是最具有应用潜力的高功率微波器件之一. 随着输出微波功率的进一步提高, 相对论返波管内部包括收集极处的击穿现象越来越严重, 最终导致脉冲缩短, 成为器件向高功率、高能量方向发展中的最大障碍, 也是目前制约其发展的重要问题之一. 本文基于自主研发的2.5维粒子模拟软件UNIPIC-2D, 采用动态释气模型研究了不同释气系数下相对论返波管收集极释气与电离过程及引导磁场的影响. 粒子模拟结果表明, 随着电子不断轰击收集极, 收集极表面气压升高, 并发生气体电离, 产生的等离子体沿引导磁场进入慢波结构区域, 影响束-波相互作用过程, 使得输出功率下降; 随着释气系数的增大, 脉冲缩短现象越来越明显; 在低引导磁场情况下, 击穿以及脉冲缩短现象得到一定的缓解.
      通信作者: 王洪广, wanghg@xjtu.edu.cn
    • 基金项目: 国家自然科学基金委员会与中国工程物理研究院联合基金(批准号: U1530133)和陕西省自然基金面上项目(批准号: 2018JM6062)资助的课题
    [1]

    Barker R J, Schamiloglu E 2001 High-Power Microwave Sources and Technologies (Piscataway, New Jersey: IEEE Press) pp310–380

    [2]

    宫玉彬, 张章, 魏彦玉, 孟凡宝, 范植开, 王文祥 2004 物理学报 53 3990

    Gong Y B, Zhang Z, Wei Y Y, Meng F B, Fan Z K, Wang W X 2004 Acta Phys. Sin. 53 3990

    [3]

    Li X Z, Wang J G, Song Z M, Chen C H, Sun J, Zhang X W, Zhang Y C 2012 Phys. Plasmas 19 83111

    [4]

    李小泽, 王建国, 童长江, 张海 2008 物理学报 57 4613

    Li X Z, Wang J G, Tong C J, Zhang H 2008 Acta Phys. Sin. 57 4613

    [5]

    Benford J, Benford G 1997 IEEE Trans. Plasma Sci. 25 311

    [6]

    Insepov Z, Norem J, Vetizer S, Mahalingam S 2011 AIP Conf. Proc. 1406 523

    [7]

    Cao Y B, Song Z M, Wu P, Fan Z Q, Zhang Y C, Teng Y, Sun J 2017 Phys. Plasmas 24 033109

    [8]

    Korovin S D, Mesyats G A, Pegel I V, Polevin S D, Tarakanov V P 2000 IEEE Trans. Plasma Sci. 28 485

    [9]

    Xiao R Z, Chen C H, Deng Y Q, Cao Y B, Sun J, Li J W 2016 Phys. Plasmas 23 063114

    [10]

    Zhang J, Jin Z X, Yang J H, Zhong H H, et al. 2011 IEEE Trans. Plasma Sci. 39 1438

    [11]

    Kovalev N F, Nechaev V E, Petelin M I 1998 IEEE Trans. Plasma Sci. 26 246

    [12]

    梁玉钦, 邵浩, 孙钧, 等 2014 强激光与粒子束 26 26063010

    Liang Y Q, Shao H, Sun J, et al. 2014 High Power Laser and Particle Beams 26 26063010

    [13]

    Miao T Z, Bai X C, Sun J, Zhang X W, Cao Y B, Wu P, Shi Y C, Shao H 2017 Phys. Plasmas 24 123106

    [14]

    傅竹风, 胡友秋 1995 空间等离子体数值模拟 (合肥: 安徽科学技术出版社) 第433−476页

    Fu Z F, Hu Y Q 1995 Numerical Simulation of Space Plasma (Hefei: Anhui Science and Technology Publishers) pp433−476 (in Chinese)

    [15]

    Wang J G, Zhang D H, Liu C L, Li Y D, Wang Y, Wang H G, Qiao H L, Li X Z 2009 Phys. Plasmas 16 033108

    [16]

    Bird R B, Lightfoot E N, Stewart W E 1961 AIChE J. 7 5J

    [17]

    Birdsall C K 1991 IEEE Trans. Plasma Sci. 19 65

    [18]

    Wang H G, Li Y D, Liu C L, Zhou Y, Liu M Q 2010 IEEE Trans. Plasma Sci. 38 2062

    [19]

    IAEA http://www-amdis.iaea.org/ALADDIN [2019-4-15]

    [20]

    董烨, 董志伟, 周前红, 杨温渊, 周海京 2014 物理学报 63 027901

    Dong Y, Dong Z W, Zhou Q H, Yang W Y, Zhou H J 2014 Acta Phys. Sin. 63 027901

    [21]

    蔡利兵, 王建国 2011 物理学报 60 025217

    Cai L B, Wang J G 2011 Acta Phys. Sin. 60 025217

    [22]

    Vaughan R M 1988 IEEE Trans. Electron Dev. 35 1172

    [23]

    杨文晋, 李永东, 刘纯亮 2013 物理学报 62 087901

    Yang W J, Li Y D, Liu C L 2013 Acta Phys. Sin. 62 087901

    [24]

    李姝敏, 李永东, 刘震 2017 强激光与粒子束 29 29063001

    Li S M, Li Y D, Liu Z 2017 High Power Laser and Particle Beams 29 29063001

    [25]

    邵剑波, 马乔生, 谢鸿全, 李正红 2015 微波学报 31 62

    Shao J B, Ma Q S, Xie H Q, Li X H 2015 J. Microw. 31 62

  • 图 1  RBWO收集极释气的PIC/MCC模拟流程图

    Fig. 1.  PIC/MCC simulation flow chart for RBWO collecting extremely outgassing.

    图 2  3.4 T引导磁场RBWO示意图

    Fig. 2.  3.4 T-guide magnetic field RBWO schematic.

    图 3  无释气情况下3.4 T RBWO输出功率

    Fig. 3.  3.4 T RBWO output power without outgassing.

    图 4  不同释气系数λ下输出功率对比

    Fig. 4.  Output power comparison under different outgassing coefficients λ.

    图 5  无释气情况下的模拟结果 (a)电子实空间分布(b)电子z-Vz相空间分布

    Fig. 5.  Simulation results without outgassing: (a) Electronic real-time spatial distribution; (b) electronic z-Vz phase spatial distribution

    图 6  释气系数λ = 2情况下的模拟结果 (a)电子实空间分布; (b)电子相空间分布; (c)−(f)分别为4, 16, 28, 40 ns时离子实空间分布; (g)收集极表面气压随时间的变化

    Fig. 6.  Simulation result with outgassing coefficient λ of 2: (a) Electronic real-time spatial distribution; (b) electronic phase spatial distribution; (c)−(f) polar space distribution at 4, 16, 28, 40 ns; (g) surface pressure curve over time.

    图 7  λ = 2时慢电子占空间电子总数的比率

    Fig. 7.  Ratio of slow electrons to the total number of space electrons at λ = 2.

    图 8  微波输出功率以及脉宽随释气系数λ的变化

    Fig. 8.  Relationship of microwave output power and pulse width to outgassing coefficient λ.

    图 9  不同释气系数下40 ns时刻RBWO收集极表面气压(一个网格内的平均气压)的模拟结果

    Fig. 9.  Simulation results of RBWO collector surface pressure (average pressure in a grid) at 40 ns with different outgassing coefficients.

    图 10  不同引导磁场下无释气情况与λ = 2时平均输出功率的对比图

    Fig. 10.  Comparison of the average output power with no outgassing and λ = 2 under different guiding magnetic fields.

  • [1]

    Barker R J, Schamiloglu E 2001 High-Power Microwave Sources and Technologies (Piscataway, New Jersey: IEEE Press) pp310–380

    [2]

    宫玉彬, 张章, 魏彦玉, 孟凡宝, 范植开, 王文祥 2004 物理学报 53 3990

    Gong Y B, Zhang Z, Wei Y Y, Meng F B, Fan Z K, Wang W X 2004 Acta Phys. Sin. 53 3990

    [3]

    Li X Z, Wang J G, Song Z M, Chen C H, Sun J, Zhang X W, Zhang Y C 2012 Phys. Plasmas 19 83111

    [4]

    李小泽, 王建国, 童长江, 张海 2008 物理学报 57 4613

    Li X Z, Wang J G, Tong C J, Zhang H 2008 Acta Phys. Sin. 57 4613

    [5]

    Benford J, Benford G 1997 IEEE Trans. Plasma Sci. 25 311

    [6]

    Insepov Z, Norem J, Vetizer S, Mahalingam S 2011 AIP Conf. Proc. 1406 523

    [7]

    Cao Y B, Song Z M, Wu P, Fan Z Q, Zhang Y C, Teng Y, Sun J 2017 Phys. Plasmas 24 033109

    [8]

    Korovin S D, Mesyats G A, Pegel I V, Polevin S D, Tarakanov V P 2000 IEEE Trans. Plasma Sci. 28 485

    [9]

    Xiao R Z, Chen C H, Deng Y Q, Cao Y B, Sun J, Li J W 2016 Phys. Plasmas 23 063114

    [10]

    Zhang J, Jin Z X, Yang J H, Zhong H H, et al. 2011 IEEE Trans. Plasma Sci. 39 1438

    [11]

    Kovalev N F, Nechaev V E, Petelin M I 1998 IEEE Trans. Plasma Sci. 26 246

    [12]

    梁玉钦, 邵浩, 孙钧, 等 2014 强激光与粒子束 26 26063010

    Liang Y Q, Shao H, Sun J, et al. 2014 High Power Laser and Particle Beams 26 26063010

    [13]

    Miao T Z, Bai X C, Sun J, Zhang X W, Cao Y B, Wu P, Shi Y C, Shao H 2017 Phys. Plasmas 24 123106

    [14]

    傅竹风, 胡友秋 1995 空间等离子体数值模拟 (合肥: 安徽科学技术出版社) 第433−476页

    Fu Z F, Hu Y Q 1995 Numerical Simulation of Space Plasma (Hefei: Anhui Science and Technology Publishers) pp433−476 (in Chinese)

    [15]

    Wang J G, Zhang D H, Liu C L, Li Y D, Wang Y, Wang H G, Qiao H L, Li X Z 2009 Phys. Plasmas 16 033108

    [16]

    Bird R B, Lightfoot E N, Stewart W E 1961 AIChE J. 7 5J

    [17]

    Birdsall C K 1991 IEEE Trans. Plasma Sci. 19 65

    [18]

    Wang H G, Li Y D, Liu C L, Zhou Y, Liu M Q 2010 IEEE Trans. Plasma Sci. 38 2062

    [19]

    IAEA http://www-amdis.iaea.org/ALADDIN [2019-4-15]

    [20]

    董烨, 董志伟, 周前红, 杨温渊, 周海京 2014 物理学报 63 027901

    Dong Y, Dong Z W, Zhou Q H, Yang W Y, Zhou H J 2014 Acta Phys. Sin. 63 027901

    [21]

    蔡利兵, 王建国 2011 物理学报 60 025217

    Cai L B, Wang J G 2011 Acta Phys. Sin. 60 025217

    [22]

    Vaughan R M 1988 IEEE Trans. Electron Dev. 35 1172

    [23]

    杨文晋, 李永东, 刘纯亮 2013 物理学报 62 087901

    Yang W J, Li Y D, Liu C L 2013 Acta Phys. Sin. 62 087901

    [24]

    李姝敏, 李永东, 刘震 2017 强激光与粒子束 29 29063001

    Li S M, Li Y D, Liu Z 2017 High Power Laser and Particle Beams 29 29063001

    [25]

    邵剑波, 马乔生, 谢鸿全, 李正红 2015 微波学报 31 62

    Shao J B, Ma Q S, Xie H Q, Li X H 2015 J. Microw. 31 62

  • [1] 宫玉彬, 张 章, 魏彦玉, 孟凡宝, 范植开, 王文祥. 高功率微波器件中脉冲缩短现象的粒子模拟. 物理学报, 2004, 53(11): 3990-3995. doi: 10.7498/aps.53.3990
    [2] 黄华, 范植开, 谭杰, 马乔生, 甘延青, 常安碧. 长脉冲相对论速调管中束流脉冲缩短的研究. 物理学报, 2004, 53(4): 1129-1135. doi: 10.7498/aps.53.1129
    [3] 陈昭福, 黄华, 常安碧, 许州, 何琥, 雷禄容, 胡进光, 袁欢, 刘振帮. S波段长脉冲相对论速调管放大器中杂模振荡导致的脉冲缩短抑制研究. 物理学报, 2014, 63(23): 238402. doi: 10.7498/aps.63.238402
    [4] 刘艺, 杨佳, 李兴, 谷伟, 高志鹏. 微秒脉冲电场下Pb0.99(Zr0.95Ti0.05)0.98Nb0.02O3陶瓷击穿过程电阻变化规律. 物理学报, 2017, 66(11): 117701. doi: 10.7498/aps.66.117701
    [5] 黄华, 罗雄, 雷禄容, 罗光耀, 张北镇, 金晓, 谭杰. 长脉冲相对论扩展互作用腔振荡器的初步研究. 物理学报, 2010, 59(3): 1907-1912. doi: 10.7498/aps.59.1907
    [6] 刘振帮, 金晓, 黄华, 陈怀璧. X波段长脉冲同轴多注相对论速调管放大器的分析与设计. 物理学报, 2012, 61(12): 128401. doi: 10.7498/aps.61.128401
    [7] 李小泽, 王建国, 童长江, 张 海. 充填不同气体相对论返波管特性的PIC-MCC模拟. 物理学报, 2008, 57(7): 4613-4622. doi: 10.7498/aps.57.4613
    [8] 王辉辉, 蒙林, 刘大刚, 刘腊群, 杨超. 基于相对论返波管的全三维PIC/PSO数值优化研究. 物理学报, 2013, 62(13): 138401. doi: 10.7498/aps.62.138401
    [9] 王宇, 陈再高, 雷奕安. 等离子体填充0.14 THz相对论返波管模拟. 物理学报, 2013, 62(12): 125204. doi: 10.7498/aps.62.125204
    [10] 陈再高, 王建国, 王玥, 朱湘琴, 张殿辉, 乔海亮. 相对论返波管超辐射产生与辐射的数值模拟研究. 物理学报, 2014, 63(3): 038402. doi: 10.7498/aps.63.038402
    [11] 杨德文, 陈昌华, 史彦超, 肖仁珍, 滕雁, 范志强, 刘文元, 宋志敏, 孙钧. X波段高效率速调型相对论返波管研究. 物理学报, 2020, 69(16): 164102. doi: 10.7498/aps.69.20200434
    [12] 李 潇, 张海英, 尹军舰, 刘 亮, 徐静波, 黎 明, 叶甜春, 龚 敏. 磷化铟复合沟道高电子迁移率晶体管击穿特性研究. 物理学报, 2007, 56(7): 4117-4121. doi: 10.7498/aps.56.4117
    [13] 施卫, 田立强, 王馨梅, 徐鸣, 马德明, 周良骥, 刘宏伟, 谢卫平. 高压超大电流光电导开关及其击穿特性研究. 物理学报, 2009, 58(2): 1219-1223. doi: 10.7498/aps.58.1219
    [14] 马乔生, 金晓, 绪明, 李正红, 吴洋. X波段10 GW高功率返波管振荡器设计. 物理学报, 2011, 60(10): 105201. doi: 10.7498/aps.60.105201
    [15] 李爽, 常超, 王建国, 刘彦升, 朱梦, 郭乐田, 谢佳玲. 横磁模下介质表面二次电子倍增的抑制. 物理学报, 2015, 64(13): 137701. doi: 10.7498/aps.64.137701
    [16] 张力, 林志宇, 罗俊, 王树龙, 张进成, 郝跃, 戴扬, 陈大正, 郭立新. 具有p-GaN岛状埋层耐压结构的横向AlGaN/GaN高电子迁移率晶体管. 物理学报, 2017, 66(24): 247302. doi: 10.7498/aps.66.247302
    [17] 王辉辉, 刘大刚, 蒙林, 刘腊群, 杨超, 彭凯, 夏蒙重. 气体电离的全三维电磁粒子模拟/蒙特卡罗数值研究. 物理学报, 2013, 62(1): 015207. doi: 10.7498/aps.62.015207
    [18] 陈再高, 王建国, 王玥, 乔海亮, 郭伟杰, 张殿辉. 基于粒子模拟和并行遗传算法的高功率微波源优化设计. 物理学报, 2013, 62(16): 168402. doi: 10.7498/aps.62.168402
    [19] 郝建红, 丁 武, 董志伟. 磁绝缘传输线振荡器中的次级电子倍增现象. 物理学报, 2006, 55(9): 4789-4794. doi: 10.7498/aps.55.4789
    [20] 吴洋, 周自刚. S波段低磁场高效率相对论返波管振荡器研究. 物理学报, 2019, 68(19): 194101. doi: 10.7498/aps.68.20182155
  • 引用本文:
    Citation:
计量
  • 文章访问数:  1317
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-16
  • 修回日期:  2019-06-23
  • 上网日期:  2019-09-01
  • 刊出日期:  2019-09-20

收集极释气对相对论返波管影响的粒子模拟

  • 1. 西安交通大学, 电子物理与器件教育部重点实验室, 西安 710049
  • 2. 西安交通大学电子科学与工程学院, 西安 710049
  • 3. 西北核技术研究所, 西安 710024
  • 通信作者: 王洪广, wanghg@xjtu.edu.cn
    基金项目: 国家自然科学基金委员会与中国工程物理研究院联合基金(批准号: U1530133)和陕西省自然基金面上项目(批准号: 2018JM6062)资助的课题

摘要: 相对论返波管被认为是最具有应用潜力的高功率微波器件之一. 随着输出微波功率的进一步提高, 相对论返波管内部包括收集极处的击穿现象越来越严重, 最终导致脉冲缩短, 成为器件向高功率、高能量方向发展中的最大障碍, 也是目前制约其发展的重要问题之一. 本文基于自主研发的2.5维粒子模拟软件UNIPIC-2D, 采用动态释气模型研究了不同释气系数下相对论返波管收集极释气与电离过程及引导磁场的影响. 粒子模拟结果表明, 随着电子不断轰击收集极, 收集极表面气压升高, 并发生气体电离, 产生的等离子体沿引导磁场进入慢波结构区域, 影响束-波相互作用过程, 使得输出功率下降; 随着释气系数的增大, 脉冲缩短现象越来越明显; 在低引导磁场情况下, 击穿以及脉冲缩短现象得到一定的缓解.

English Abstract

    • 高功率微波(high power microwave, HPM)技术[1]在国防、航天、能源和通信等领域具有广阔的应用前景, 随着研究的不断深入, 正在从技术研究阶段向应用阶段加速前进. HPM源作为HPM技术的核心, 一般工作在中低真空状态(约10–2 Pa水平), 输出功率在GW量级, 电磁场功率流密度通常高达数百MW/cm2. 随着微波输出功率的不断提高, 器件内部强电磁场下的真空击穿导致微波输出功率下降以及脉冲缩短现象, 严重限制了高功率、高能量的输出, 极大地影响了器件的工作性能. 相对论返波管(relativistic backward wave oscillator, RBWO)是一种典型的HPM源, 利用电子束与慢波结构–1次空间谐波的相互作用, 在慢波结构里激励起电磁场并辐射出去. RBWO具有高功率、高效率、可重复频率运行等特点[24], 被普遍认为是最具有研究价值和应用潜力的HPM器件之一.

      在HPM的产生以及传输过程中, RBWO的高频结构、电子束收集极和介质窗三个区域是强电磁场真空击穿的薄弱环节. 强电磁场下真空击穿导致大量等离子体产生, 严重影响HPM的产生与传输条件, 造成微波功率下降以及脉冲缩短, 甚至导致结构破坏, 大幅降低HPM源的可靠性和寿命, 成为RBWO向更高功率发展的主要限制因素, 从而影响其在各个领域的推广及应用[5,6]. 电子束收集极是RBWO中的重要部分, 其主要作用是在外加引导磁场的约束下, 收集经过束-波互作用后的残余电子并且将其能量有效释放. 对于GW级输出功率的RBWO而言, 收集极不可避免地会产生大量等离子体, 其中包括解吸附气体的电子碰撞电离. 较高密度的等离子体沿着外加引导磁场向上游扩散进入高频结构, 可能导致RBWO输出功率下降甚至引发微波脉冲缩短[6,7]. 因此, 本文研究收集极等离子体的扩散及其对器件的影响机制, 对提升RBWO系统的功率、重复频率、长时间稳定工作的能力具有重要的现实意义.

      高密度等离子体影响了微波与强相对论电子束(intense relativistic electron beams, IREBs)之间的功率转换效率, 导致微波功率下降甚至脉冲缩短[711]. 同时有文献表明引导磁场的大小对收集极上的能量沉积分布也有较大的影响[12]. 但是对于RBWO中的收集极等离子体研究较少, 近些年研究人员才开始关注收集极击穿所产生等离子体对RBWO的影响. 然而目前已有的文献均采用固定气压的静态背景气体模型, 与实际情况存在差异[13]. 为了更加接近实际, 本文基于自主研发的2.5维粒子模拟(particle-in-cell, PIC)[14]软件UNIPIC-2D[15], 采用动态释气模型研究不同引导磁场、不同释气系数下的收集极释气与电离过程, 通过粒子模拟获得RBWO的输出功率、电子实空间和相空间分布等结果.

    • 本文研究释气现象对RBWO收集极击穿的影响问题, 考虑到在HPM产生过程中, 与高频结构相互作用后的IREBs仍携带很高能量不断轰击到收集极上, 并且将残余电子剩余的能量以热能的形式沉积在电子束收集极的收集面内, 并轰击出二次电子, 还引起收集极表面吸附气体分子的解吸附; 当沉积在金属中的能量达到临界值时, 会引起金属材料烧蚀气化. 由于烧蚀过程过于复杂, 在本文模拟中尚未考虑, 将其近似等效为材料表面释放中性气体分子, 通过设置释气系数λ为单个电子轰击收集极表面后释放气体分子的个数来实现. 气体解吸附后将发生扩散、碰撞电离等过程, 形成的等离子体对RBWO工作产生影响, 本文通过在PIC方法中考虑气体扩散和气体碰撞过程, 对这种现象进行模拟研究, 流程图如图1所示.

      图  1  RBWO收集极释气的PIC/MCC模拟流程图

      Figure 1.  PIC/MCC simulation flow chart for RBWO collecting extremely outgassing.

      对于解吸附气体的扩散过程, 采用Malthus反应[16]扩散方程描述:

      $\begin{split}\frac{{\partial {n_{{\rm{neutral}}}}\left( { {{r}} ,t} \right)}}{{\partial t}} =\, & \nabla \cdot \left( {{D_{{\rm{neutral}}}}\left( { {{r}} ,t} \right)\nabla {n_{{\rm{neutral}}}}} \right) \\ &- L\left( {{{r}} ,t} \right) + G\left( {{{r}} ,t} \right),\end{split}$

      其中, ${n_{{\rm{neutral}}}}$为中性气体分子密度, ${D_{{\rm{neutral}}}}$为中性气体分子扩散系数, L为损耗项(气体电离), G增益项(电子与离子复合). 该方程可以描述低气压气体的扩散过程, 其中扩散系数与气压成反比:

      $\begin{split}{D_{{\rm{neutral}}}}\left[ {\frac{{{\rm{c}}{{\rm{m}}^2}}}{{\rm{s}}}} \right] =\, & {C_1}{\left( {\frac{{T\left[ K \right]}}{{\sqrt {{T_{{\rm{c,A}}}}{T_{{\rm{c,B}}}}} }}} \right)^{{C_2}}}{\left( {{T_{{\rm{c,A}}}}{T_{{\rm{c,B}}}}} \right)^{{5 / {12}}}}\\ & \times{\left( {\frac{1}{{{M_{\rm{A}}}}} + \frac{1}{{{M_{\rm{B}}}}}} \right)^{{1 / 2}}}\frac{{{{\left( {{p_{{\rm{c,A}}}}{p_{{\rm{c,B}}}}} \right)}^{{1 / 3}}}}}{{p\left[ {{\rm{atm}}} \right]}},\end{split}$

      其中, ${C_1},{C_2}$为气体常数, 对于非极性气体一般都是定值; ${T_{{\rm{c,S}}}}$为中性气体分子S的临界温度(K); ${M_{\rm{S}}}$为中性气体分子S的分子质量(u); ${p_{{\rm{c,S}}}}$为中性气体分子S的临界气压(atm, 1 atm = 1.01325 × 105 Pa). 对于H2, 模拟中${p_{{\rm{c,S}}}}$取12.8, ${T_{{\rm{c,S}}}}$取33.3, ${M_{\rm{S}}}$取2.016.

      UNIPIC软件中采用蒙特卡罗碰撞(Monte Carlo collision, MCC)[17,18]模型来处理气体碰撞过程. 一般的MCC是采用碰撞时间随机的方法, 而粒子模拟方法中, 粒子推进和场推进的时间步长是固定的. 因此在PIC方法中, 采用了在一个时间步长内随机决定粒子间是否发生碰撞的方法来实现MCC. 目前, UNIPIC软件中的MCC模块能模拟带电粒子与中性原子之间的弹性碰撞、激发碰撞和电离碰撞等复杂的物理过程[4]. 在收集极等离子体中, 氢离子是主要成分之一, 并且由于质量小的离子运动快, 对器件的影响更为显著, 因此本文研究中解吸附气体采用H2, 在模拟中主要考虑了电子与H2分子的弹性碰撞、激发和电离碰撞, 所使用的碰撞截面数据来自国际原子能组织的核数据库[19]. 除此之外, 在实际过程中, 还有包括电荷交换等其他的碰撞类型, 但是因为碰撞截面较小, 故忽略不计[20,21].

      二次电子发射[18,22,23]是由于“初始”电子的碰撞从而导致材料表面发射出电子的现象. 当电子运动到收集极表面时, 有可能被反射、也可能沉积在材料中, 还有可能产生二次电子发射, 电子在表面的行为对气体电离现象有重要影响. 本文采用高能二次电子发射模型[23], 二次电子发射系数为

      $ \begin{split}\delta = \, & \dfrac{{785\dfrac{{{b_1}}}{\varepsilon } \cdot \dfrac{{\rho Z}}{A} \cdot \dfrac{{\left( {{e^3} - 1} \right)}}{{\beta {e^3}}}}}{{{E_0}\cos \theta }} \cdot \ln \left( {\dfrac{{1.166{E_0}}}{J}} \right) \\ & + {b_2} \cdot {{\rm{e}}^{ - {{\left| {{E_0} \cdot \cos \theta - 10^4} \right|} / {{(E_{{{\rm{1}} / {\rm{e}}}}} - 10^4)}}}},\end{split} $

      $J = \left\{ {\begin{aligned} & {\left( {9.76Z + \frac{{58.5}}{{{Z^{0.19}}}}} \right)\;({\rm eV}),\;Z \geqslant 13,} \\ & {11.5Z\;({\rm eV}),\;~~~~~~~~~~~~~~~~~~~~~ Z \leqslant 12,} \end{aligned}} \right.$

      其中, E0为入射能量(eV); $\varepsilon $为产生一个二次电子所需要消耗的能量, 对于同一种材料可认为$\varepsilon $是定值; $\rho $是材料密度(g/cm3); Z是原子序数; A是材料的原子量; J是材料的平均电离能(eV), $\beta $是材料对二次电子的吸收系数, b1$x = 0$时的二次电子逸出概率, $J,\;{b_1}{\rm{,}}\;\beta ,\;\varepsilon $由材料本身决定; b2E0 = 104 eV时的背散射电子发射系数, ${E_{{{\rm{1}} / {\rm{e}}}}}$为背散射电子发射系数降低到b2/e时所对应的入射电子能量值; $\theta $为入射角度.

    • 采用的3.4 T引导磁场RBWO模型结构[24]图2所示, 图中I为爆炸发射阴极, 在脉冲功率源驱动下发射强流相对论电子束; II, III, IV为高频结构区域, 其中II为半径2.8 cm、宽度1 cm的谐振腔反射器; III为慢波结构, 采用非均匀梯形结构以保证较高的工作效率; IV为半径3.2 cm、宽度0.7 cm的提取腔. 在模拟过程中的工作电压为700 kV, 电流为10 kA, 电子束流内外半径分别为1.9 cm和2.0 cm; 引导磁场3.4 T; 在收集极的局部区域(25 cm < z < 30 cm, r = 2.2 cm)收集电子束并产生释气和二次电子; 模拟时间为40 ns.

      图  2  3.4 T引导磁场RBWO示意图

      Figure 2.  3.4 T-guide magnetic field RBWO schematic.

    • 在无释气情况下, RBWO输出微波功率如图3所示, 输出的平均功率为1.94 GW左右, 输出微波功率比较稳定.

      图  3  无释气情况下3.4 T RBWO输出功率

      Figure 3.  3.4 T RBWO output power without outgassing.

      RBWO收集极产生释气与无释气情况下的模拟结果进行对比, 如图4所示. 相比无释气情况, 微波输出功率、脉冲宽度随着释气系数增加呈现减小趋势, 在释气系数为2时脉冲缩短现象明显, 在30 ns时器件的微波输出截止.

      图  4  不同释气系数λ下输出功率对比

      Figure 4.  Output power comparison under different outgassing coefficients λ.

      图5为无释气情况下40 ns时的电子分布, 其中图5(a)为电子z-r实空间分布, 图5(b)为电子z-Vz相空间分布. 从图5可以看出, 在无释气情况下, 慢波结构区域中电子束与高频电场相互作用, 发生能量转换和速度调制; 非均匀慢波结构的前段部分对电子束进行简单的预调制, 慢波结构的后半部分与电子束相互作用, 大部分电子通过慢波结构末端实现了电子与微波场的能量交换, 从而保证了器件的高效率运行. 从相空间图上可以看到有少部分电子反向运动, 这主要是由空间电荷效应和群聚作用过程引起的.

      图  5  无释气情况下的模拟结果 (a)电子实空间分布(b)电子z-Vz相空间分布

      Figure 5.  Simulation results without outgassing: (a) Electronic real-time spatial distribution; (b) electronic z-Vz phase spatial distribution

      收集极发生释气时, 在电子束的轰击下, 表面气压升高, 形成等离子体, 离子质量较大, 运动的距离较小; 电子则沿着引导磁场进入慢波结构区域, 影响束-波相互作用, 等离子体不断增加最终导致输出功率下降. 图6为释气系数λ = 2情况下的模拟结果, 图6(a)为电子实空间分布, 图6(b)为电子相空间分布, 图6(c)(f)为离子实空间分布. 从电子相空间分布图可以看到, 相比无释气情况, 此时的低速电子(慢电子)增多, 这是因为在收集极释气时, 电离产生的电子沿着引导磁场进入慢波结构, 在高频结构区域受到驻波场的强调制[13], 扩散速度较快. 由于上游二极管电场的作用, 慢电子被限制在互作用区域, 形成慢电子累积, 在一定程度上屏蔽了电场, 影响了正常的IREBs调制, 并且从微波场中获得能量导致微波输出功率降低. 图7为释气系数λ = 2时慢电子数占空间电子总数的比率, 随着时间的推移不断增大, 在10 ns时慢电子占约40%, 接近主电子束的密度, 从图4可以看出此时微波输出功率相比无释气情况已降低了约一半; 随着慢电子比例的进一步增大, 最终导致微波截止.

      图  6  释气系数λ = 2情况下的模拟结果 (a)电子实空间分布; (b)电子相空间分布; (c)−(f)分别为4, 16, 28, 40 ns时离子实空间分布; (g)收集极表面气压随时间的变化

      Figure 6.  Simulation result with outgassing coefficient λ of 2: (a) Electronic real-time spatial distribution; (b) electronic phase spatial distribution; (c)−(f) polar space distribution at 4, 16, 28, 40 ns; (g) surface pressure curve over time.

      图  7  λ = 2时慢电子占空间电子总数的比率

      Figure 7.  Ratio of slow electrons to the total number of space electrons at λ = 2.

      为了分析收集极产生等离子体对RBWO输出微波性能的影响规律, 在其他条件均保持不变的前提下改变释气系数, 得到最大输出功率、脉宽与释气系数的关系, 如图8所示. 结果表明, 当收集极产生释气时, RBWO的输出微波功率开始出现下降, 随着释气系数的增大, 单位时间内产生的气体分子增加, 等离子体密度上升, 输出功率迅速减小, 微波脉冲宽度缩短. 在3.4 T引导磁场下, 当释气系数为5时输出微波几乎完全截止, 可以推测此时器件已经发生严重的击穿.

      图  8  微波输出功率以及脉宽随释气系数λ的变化

      Figure 8.  Relationship of microwave output power and pulse width to outgassing coefficient λ.

      图9为不同释气系数下40 ns时刻RBWO收集极表面气压的模拟结果, 显示出随着λ的增大表面气压也随之增大, 结合图8可知当收集极表面气压大于0.0105 Torr (1 Torr = 1.33322 × 102 Pa)时, 对器件影响较大; 收集极表面气压气压小于0.0105 Torr时, 在40 ns内器件输出功率的影响较小.

      图  9  不同释气系数下40 ns时刻RBWO收集极表面气压(一个网格内的平均气压)的模拟结果

      Figure 9.  Simulation results of RBWO collector surface pressure (average pressure in a grid) at 40 ns with different outgassing coefficients.

      为了进一步分析释气对RBWO收集极击穿的影响, 在该器件的适用引导磁场范围内, 计算了不同引导磁场下该器件收集极产生释气后击穿的模拟情况. 图10所示为该器件在不同引导磁场下无释气情况与释气系数λ = 2时平均输出功率的对比. 可以看出在无释气情况下, 随着引导磁场的减小, 器件输出功率递减; 而在释气系数λ = 2时, 随着引导磁场的减小, 器件微波输出功率基本呈递增状态, 这可能是因为随着引导磁场的降低, 对电子束的约束变小, 轰击收集极的面积增大, 从而释气气压降低; 同时引导磁场对等离子体的约束变小, 等离子体分布范围变大, 而密度降低, 对主电子束的影响减小, 所以输出微波功率逐渐增加. 随着HPM技术的不断发展, 采用较低引导磁场是实现HPM产生系统小型化和低功耗的重要手段, 基于低磁场的永磁包装器件是当前发展方向之一[25]. 由模拟结果可知, 减小引导磁场也能够降低收集极释气对输出功率的影响.

      图  10  不同引导磁场下无释气情况与λ = 2时平均输出功率的对比图

      Figure 10.  Comparison of the average output power with no outgassing and λ = 2 under different guiding magnetic fields.

    • 本文采用2.5维UNIPIC软件, 基于动态释气与气体扩散模型研究了RBWO在不同引导磁场、不同释气系数下的收集极击穿过程, 模拟得到了RBWO的输出功率、电子实空间分布、相空间分布、等离子体分布和气压变化等结果. 模拟结果表明, 随着电子不断轰击收集极, 收集极表面气压升高, 轰击产生的二次电子以及收集极表面气体解吸附, 产生等离子体扩散到互作用区域, 使得RBWO发生输出功率下降和脉冲缩短现象; 随着收集极释气系数的不断增大, RBWO的脉冲缩短现象越来越严重; 低引导磁场在一定程度上可以缓解脉冲缩短现象.

参考文献 (25)

目录

    /

    返回文章
    返回