搜索

x

留言板

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

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

磁绝缘线振荡器中模式竞争的物理分析和数值模拟

杨温渊 董烨 孙会芳 董志伟

磁绝缘线振荡器中模式竞争的物理分析和数值模拟

杨温渊, 董烨, 孙会芳, 董志伟
PDF
HTML
导出引用
  • 作为一种紧凑型高功率微波器件, 磁绝缘线振荡器在起振过程中容易出现模式竞争现象, 如果不能对其进行有效抑制, 可能导致器件的最终输出性能下降. 由于磁绝缘线振荡器中波束互作用区通常采用同轴盘荷波导作为其慢波结构, 因此本文从同轴盘荷波导中几个可能被相对论电子束激发的低阶本征模与电子束之间的色散关系入手, 分析了三种类型的模式竞争的特点、产生的可能原因以及削弱方法. 基于以上分析, 给出了一种高功率紧凑型L波段磁绝缘线振荡器的物理模型, 并利用全电磁三维粒子程序对其进行了冷腔和热腔的数值模拟. 结果表明, 由于结构不完全对称和电子发射可能存在一定的非均匀性, 器件运行初期互作用区有竞争模式HEM11模出现, 与理论分析一致;起振一小段时间后(10 ns左右), 互作用区基模增长加快, 高阶模被抑制. 进一步优化后器件在基模获得了高效率和高功率微波输出, 饱和时输出功率约为8.1 GW, 输出效率达到了18%, 模式纯度约为97%. 本文研究结果可为磁绝缘线振荡器运行过程中出现的竞争模式的识别和输出性能优化提供理论参考和依据.
      通信作者: 杨温渊, yang_wenyuan@iapcm.ac.cn
    • 基金项目: 国家自然科学基金(批准号: 11875094)资助的课题
    [1]

    Barker R J, Schamiloglu E 2001 High-Power Microwaves Sources and Technologies (New York: Institute of Electrical and Electronic Engineer, Inc.) pp43–53

    [2]

    Eastwood J W, Hawkins K C, Hook M P 1998 IEEE Trans. Plasma Sci. 26 698

    [3]

    Lemke R W, Calico S E, Clark M C 1997 IEEE Trans. Plasma Sci. 25 364

    [4]

    Yang W Y 2008 IEEE Trans. Plasma Sci. 36 2801

    [5]

    Fan Y W, Li S R, Wang X Y, Li A K, Yu Y Q, Liu Z Y 2019 Rev. Sci. Instrum. 90 044704

    [6]

    董志伟, 孙会芳, 杨郁林, 杨温渊, 周前红, 张芳, 董烨 2016 强激光与粒子束 28 033023

    Dong Z W, Sun H F, Yang Y L, Yang W Y, Zhou Q H, Zhang F, Dong Y 2016 High Pow. Las. Part. Beam. 28 033023

    [7]

    Wang X Y, Fan Y W, Shu T, Li A K, Yu Y Q, Liu Z Y 2019 IEEE Trans. Plasma Sci. 47 3974

    [8]

    Nallasamy V, Datta S K, Reddy S U, Jain P K 2017 J. Electromagn. Waves Appl. 31 1864

    [9]

    Kumar A, Dwivedi S, Jain P K 2019 IEEE Trans. Plasma Sci. 47 4642

    [10]

    Kim D H, Jung H C, Min S H, Shin S H, Rhee M J, Park G S, Kim C H, Yim D W 2006 7th IEEE International Vacuum Electronics Conference (IVEC)/6th IEEE International Vacuum Electron Sources Conference (IVESC) Monterey, CA, April 25–27, 2006 p352

    [11]

    Dixit G, Kumar A, Jain P K 2017 Phys. Plasma. 24 013113

    [12]

    Nallasamy V, Narasimhamurthy C, Geetha B, Gupta S K, Datta S K, Reddy S U, Jain P K 2017 J. Electromagn. Waves Appl. 31 375

    [13]

    Qin F, Wang D, Xu S, Zhang Y, Fan Z K 2016 Rev. Sci. Instrum. 87 044703

    [14]

    孙会芳, 董志伟, 杨郁林 2010 强激光与粒子束 22 303

    Sun H F, Dong Z W, Yang Y L 2010 High Pow. Las. Part. Beam. 22 303

    [15]

    Kim D H, Jung H C, Min S H, Shin S H, Park G S 2007 Appl. Phys. Lett. 90 124103

    [16]

    Lemke R W 1989 J. Appl. Phys. 66 1089

    [17]

    Cousin R, Larour J, Gardelle J, Cassany B, Modin P, Gouard P, Raymond P 2007 IEEE Trans. Plasma Sci. 35 1467

    [18]

    Cousin R, Larour J, Gouard P, Raymond P 2006 J. Appl. Phys. 100 084512

    [19]

    王冬, 陈代兵, 范植开, 邓景康 2008 物理学报 57 4875

    Wang D, Chen D B, Fan Z K, Deng J K 2008 Acta Phys. Sin. 57 4875

    [20]

    Jiang T, Zhang J D, He J T, Li Z Q, Ling J P 2016 IEEE Trans. Plasma Sci. 44 755

    [21]

    秦奋, 王冬, 陈代兵, 文杰 2012 物理学报 61 094101

    Qin F, Wang D, Chen D B, Wen J 2012 Acta Phys. Sin. 61 094101

    [22]

    董烨, 董志伟, 杨温渊, 周海京 2009 强激光与粒子束 21 1199

    Dong Y, Dong Z W, Yang W Y, Zhou H J 2009 High Pow. Las. Part. Beam. 21 1199

    [23]

    姜利辉, 李浩, 吴泽威 2014 强激光与粒子束 26 063009

    Jiang L H, Li H, Wu Z W 2014 High Pow. Las. Part. Beam. 26 063009

  • 图 1  同轴盘荷波导的结构示意图 (a) 纵向中心截面; (b) 横截面

    Fig. 1.  Schematic drawings of the coaxial disk-loaded waveguide in (a) Axial and (b) radial cross section

    图 2  同轴盘荷波导中TM00模、TM01模、HEM11模以及电子束和光束的色散关系图

    Fig. 2.  The uncoupled dispersion curves of the coaxial disk-loaded waveguide (TM00, TM01 and HEM11), and the dispersion lines of the electron and light beam.

    图 3  同轴盘荷波导中TM00(π模)的轴向电场在(a) 纵向中心截面和(b) 横截面的等高图

    Fig. 3.  The contour plots of Ez of the TM00 (π mode) in the coaxial disk-loaded waveguide in the (a) Axial and (b) radial cross section.

    图 4  同轴盘荷波导中TM01的轴向电场在横截面的等高图

    Fig. 4.  The contour plots of Ez of the TM01 mode in the coaxial disk-loaded waveguide in the radial cross section.

    图 5  同轴盘荷波导中HEM11(π模)的轴向电场在 (a) 纵向中心截面和(b)横截面的等高图

    Fig. 5.  The contour plots of Ez of the HEM11(π mode)in the coaxial disk-loaded waveguide in the (a) Axial and (b) radial cross section.

    图 6  紧凑型MILO纵向截面示意图

    Fig. 6.  Schematic drawing of the compact MILO in the axial cross section.

    图 7  短脉冲宽频信号激发后MILO互作用区中点D处电场的频谱图

    Fig. 7.  The oscillation frequency of the Ez stimulated in cold cavity at point D.

    图 8  频率为1.61 GHz的纵向电场在 (a) 纵向和(b) 横向截面等高图

    Fig. 8.  The contour plots of Ez with the frequency of 1.61 GHz in (a) The axial and (b) the radial cross section.

    图 9  频率为1.77 GHz的纵向电场在 (a) 纵向和(b) 横向截面的等高图

    Fig. 9.  The contour plots of Ez with the frequency of 1.77 GHz in (a) The axial and (b) the radial cross section.

    图 10  (a) 阳极腔内观察点D的电场Ez随时间的变化曲线及不同时间窗口(b) 0−50 ns, (c) 8−14 ns和(d) 14−50 ns的FFT变换图

    Fig. 10.  (a) Variations of Ez with time at the observation point D in the anode cavity; the corresponding Fourier transform with different time intervals: (b) 0−50 ns; (c) 8−14 ns; (d) 14−50 ns.

    图 11  (a) 阳极腔内观察点D磁场By随时间的变化曲线及不同时间窗口(b) 0−50 ns, (c) 8−14 ns和(d) 14−50 ns的FFT变换图

    Fig. 11.  (a) Variations of By with time at the observation point D in the anode cavity, and the corresponding Fourier transform with different time intervals: (b) 0−50 ns, (c) 8−14 ns (d) 14−50 ns.

    图 12  初始阶段不同时刻互作用区横截面轴向电场的等高图 (a) t = 12.171 ns; (b) t = 13.094 ns; (c) t = 14.944 ns; (d) t = 16.793 ns

    Fig. 12.  Contour plots of Ez in the interaction region at different times at the initial stage.

    图 13  饱和时互作用区轴向电场在(a) 纵向截面和(b) 横向截面的等高图, t = 45.727 ns

    Fig. 13.  Contour plots of Ez in the (a) Axial and (b) radial cross section in the interaction region at saturation.

    图 14  (a) MILO输出周期平均功率随时间的变化曲线; (b)稳定后输出口电场的FFT变换图

    Fig. 14.  (a) Time plots of the periodic-average output power of the MILO; (b) the Fourier transform of Ex at the output port.

    图 15  电磁场在输出口横向截面分布的箭矢图 (a) 电场; (b)磁场

    Fig. 15.  The vector plots of the (a) Electric field and (b) magnetic field in the radial cross section at the output port.

  • [1]

    Barker R J, Schamiloglu E 2001 High-Power Microwaves Sources and Technologies (New York: Institute of Electrical and Electronic Engineer, Inc.) pp43–53

    [2]

    Eastwood J W, Hawkins K C, Hook M P 1998 IEEE Trans. Plasma Sci. 26 698

    [3]

    Lemke R W, Calico S E, Clark M C 1997 IEEE Trans. Plasma Sci. 25 364

    [4]

    Yang W Y 2008 IEEE Trans. Plasma Sci. 36 2801

    [5]

    Fan Y W, Li S R, Wang X Y, Li A K, Yu Y Q, Liu Z Y 2019 Rev. Sci. Instrum. 90 044704

    [6]

    董志伟, 孙会芳, 杨郁林, 杨温渊, 周前红, 张芳, 董烨 2016 强激光与粒子束 28 033023

    Dong Z W, Sun H F, Yang Y L, Yang W Y, Zhou Q H, Zhang F, Dong Y 2016 High Pow. Las. Part. Beam. 28 033023

    [7]

    Wang X Y, Fan Y W, Shu T, Li A K, Yu Y Q, Liu Z Y 2019 IEEE Trans. Plasma Sci. 47 3974

    [8]

    Nallasamy V, Datta S K, Reddy S U, Jain P K 2017 J. Electromagn. Waves Appl. 31 1864

    [9]

    Kumar A, Dwivedi S, Jain P K 2019 IEEE Trans. Plasma Sci. 47 4642

    [10]

    Kim D H, Jung H C, Min S H, Shin S H, Rhee M J, Park G S, Kim C H, Yim D W 2006 7th IEEE International Vacuum Electronics Conference (IVEC)/6th IEEE International Vacuum Electron Sources Conference (IVESC) Monterey, CA, April 25–27, 2006 p352

    [11]

    Dixit G, Kumar A, Jain P K 2017 Phys. Plasma. 24 013113

    [12]

    Nallasamy V, Narasimhamurthy C, Geetha B, Gupta S K, Datta S K, Reddy S U, Jain P K 2017 J. Electromagn. Waves Appl. 31 375

    [13]

    Qin F, Wang D, Xu S, Zhang Y, Fan Z K 2016 Rev. Sci. Instrum. 87 044703

    [14]

    孙会芳, 董志伟, 杨郁林 2010 强激光与粒子束 22 303

    Sun H F, Dong Z W, Yang Y L 2010 High Pow. Las. Part. Beam. 22 303

    [15]

    Kim D H, Jung H C, Min S H, Shin S H, Park G S 2007 Appl. Phys. Lett. 90 124103

    [16]

    Lemke R W 1989 J. Appl. Phys. 66 1089

    [17]

    Cousin R, Larour J, Gardelle J, Cassany B, Modin P, Gouard P, Raymond P 2007 IEEE Trans. Plasma Sci. 35 1467

    [18]

    Cousin R, Larour J, Gouard P, Raymond P 2006 J. Appl. Phys. 100 084512

    [19]

    王冬, 陈代兵, 范植开, 邓景康 2008 物理学报 57 4875

    Wang D, Chen D B, Fan Z K, Deng J K 2008 Acta Phys. Sin. 57 4875

    [20]

    Jiang T, Zhang J D, He J T, Li Z Q, Ling J P 2016 IEEE Trans. Plasma Sci. 44 755

    [21]

    秦奋, 王冬, 陈代兵, 文杰 2012 物理学报 61 094101

    Qin F, Wang D, Chen D B, Wen J 2012 Acta Phys. Sin. 61 094101

    [22]

    董烨, 董志伟, 杨温渊, 周海京 2009 强激光与粒子束 21 1199

    Dong Y, Dong Z W, Yang W Y, Zhou H J 2009 High Pow. Las. Part. Beam. 21 1199

    [23]

    姜利辉, 李浩, 吴泽威 2014 强激光与粒子束 26 063009

    Jiang L H, Li H, Wu Z W 2014 High Pow. Las. Part. Beam. 26 063009

  • [1] 杜朝海, 李铮迪, 薛志浩, 刘濮鲲, 薛谦忠, 张世昌, 徐寿喜, 耿志辉, 顾伟, 粟亦农, 刘高峰. W波段损耗介质加载回旋返波振荡器中模式竞争的研究. 物理学报, 2012, 61(7): 070703. doi: 10.7498/aps.61.070703
    [2] 王 冬, 陈代兵, 范植开, 邓景康. HEM11模磁绝缘线振荡器的高频分析. 物理学报, 2008, 57(8): 4875-4882. doi: 10.7498/aps.57.4875
    [3] 王冬, 陈代兵, 秦奋, 范植开. 双频磁绝缘线振荡器二维周期结构研究. 物理学报, 2009, 58(10): 6962-6972. doi: 10.7498/aps.58.6962
    [4] 秦奋, 王冬, 陈代兵, 文杰. L波段高阶模抑制型磁绝缘线振荡器研究. 物理学报, 2012, 61(9): 094101. doi: 10.7498/aps.61.094101
    [5] 陈代兵, 王冬, 秦奋, 文杰, 金晓, 安海狮, 张新凯. 磁绝缘线振荡器的起振电压与注入电压关系的分析. 物理学报, 2012, 61(1): 012901. doi: 10.7498/aps.61.012901
    [6] 陈代兵, 张运俭, 张北镇, 王冬, 秦奋, 文杰, 金晓, 吴勇, 于爱民. 磁绝缘线振荡器阴极烧蚀与电压波形的关系研究. 物理学报, 2013, 62(1): 012901. doi: 10.7498/aps.62.012901
    [7] 黄丽萍, 洪斌斌, 刘畅, 唐昌建. 220GHz三次谐波光子带隙谐振腔回旋管振荡器的研究. 物理学报, 2014, 63(11): 118401. doi: 10.7498/aps.63.118401
    [8] 杨 浩, 郭 霞, 关宝璐, 王同喜, 沈光地. 注入电流对垂直腔面发射激光器横模特性的影响. 物理学报, 2008, 57(5): 2959-2965. doi: 10.7498/aps.57.2959
    [9] 连天虹, 王石语, 寇科, 刘芸. 离轴抽运厄米-高斯模固体激光器. 物理学报, 2020, 69(11): 114202. doi: 10.7498/aps.69.20200086
    [10] 赵啦啦, 刘初升, 闫俊霞, 蒋小伟, 朱艳. 不同振动模式下颗粒分离行为的数值模拟. 物理学报, 2010, 59(4): 2582-2588. doi: 10.7498/aps.59.2582
    [11] 陈石, 王辉, 沈胜强, 梁刚涛. 液滴振荡模型及与数值模拟的对比. 物理学报, 2013, 62(20): 204702. doi: 10.7498/aps.62.204702
    [12] 刘漾, 巩华荣, 魏彦玉, 宫玉彬, 王文祥, 廖复疆. 有效抑制光子晶体加载矩形谐振腔中模式竞争的方法. 物理学报, 2009, 58(11): 7845-7851. doi: 10.7498/aps.58.7845
    [13] 刘腊群, 刘大刚, 王学琼, 杨超, 夏蒙重, 彭凯. 磁绝缘传输线中心汇流区电子能量沉积及温度变化的数值模拟研究. 物理学报, 2012, 61(16): 162902. doi: 10.7498/aps.61.162902
    [14] 马理强, 苏铁熊, 刘汉涛, 孟青. 微液滴振荡过程的光滑粒子动力学方法数值模拟. 物理学报, 2015, 64(13): 134702. doi: 10.7498/aps.64.134702
    [15] 靳冬欢, 刘文广, 陈星, 陆启生, 赵伊君. 三股互击式喷注器及燃烧室流场的数值模拟. 物理学报, 2012, 61(6): 064206. doi: 10.7498/aps.61.064206
    [16] 梁慧敏, 杜惊雷, 王宏波, 王治华, 罗时荣, 杨经国, 郑万国, 魏晓峰, 朱启华, 黄晓军, 王晓东, 郭 仪. 不同波长激光激发下C6H12受激拉曼散射模式竞争. 物理学报, 2007, 56(12): 6994-6998. doi: 10.7498/aps.56.6994
    [17] 李锟影, 李璞, 郭晓敏, 郭龑强, 张建国, 刘义铭, 徐兵杰, 王云才. 利用光反馈多模激光器结合滤波器产生平坦混沌. 物理学报, 2019, 68(11): 110501. doi: 10.7498/aps.68.20190171
    [18] 高启, 张传飞, 周林, 李正宏, 吴泽清, 雷雨, 章春来, 祖小涛. Z箍缩Al等离子体X特征辐射谱线数值模拟及考虑叠加效应后的修正. 物理学报, 2014, 63(12): 125202. doi: 10.7498/aps.63.125202
    [19] 成玉国, 夏广庆. 感应式脉冲推力器中等离子体加速数值研究. 物理学报, 2017, 66(7): 075204. doi: 10.7498/aps.66.075204
    [20] 耿少飞, 唐德礼, 赵杰, 邱孝明. 圆柱形阳极层霍尔等离子体加速器的质点网格方法模拟. 物理学报, 2009, 58(8): 5520-5525. doi: 10.7498/aps.58.5520
  • 引用本文:
    Citation:
计量
  • 文章访问数:  274
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-14
  • 修回日期:  2020-05-24
  • 上网日期:  2020-06-12
  • 刊出日期:  2020-10-05

磁绝缘线振荡器中模式竞争的物理分析和数值模拟

    基金项目: 国家自然科学基金(批准号: 11875094)资助的课题

摘要: 作为一种紧凑型高功率微波器件, 磁绝缘线振荡器在起振过程中容易出现模式竞争现象, 如果不能对其进行有效抑制, 可能导致器件的最终输出性能下降. 由于磁绝缘线振荡器中波束互作用区通常采用同轴盘荷波导作为其慢波结构, 因此本文从同轴盘荷波导中几个可能被相对论电子束激发的低阶本征模与电子束之间的色散关系入手, 分析了三种类型的模式竞争的特点、产生的可能原因以及削弱方法. 基于以上分析, 给出了一种高功率紧凑型L波段磁绝缘线振荡器的物理模型, 并利用全电磁三维粒子程序对其进行了冷腔和热腔的数值模拟. 结果表明, 由于结构不完全对称和电子发射可能存在一定的非均匀性, 器件运行初期互作用区有竞争模式HEM11模出现, 与理论分析一致;起振一小段时间后(10 ns左右), 互作用区基模增长加快, 高阶模被抑制. 进一步优化后器件在基模获得了高效率和高功率微波输出, 饱和时输出功率约为8.1 GW, 输出效率达到了18%, 模式纯度约为97%. 本文研究结果可为磁绝缘线振荡器运行过程中出现的竞争模式的识别和输出性能优化提供理论参考和依据.

English Abstract

    • 磁绝缘线振荡器(magnetically insulated transmission line oscillator, MILO)是一种紧凑型高功率微波器件, 它不需要外加引导磁场聚焦电子束, 输出功率可达吉瓦量级, 在许多领域有着良好的应用前景[1]. 为进一步提高MILO输出性能, 如提高其输出功率和效率, 增加模式稳定性等, 研究人员一直在进行相关研究和探索[2-21].

      为了获得高效稳定的高功率微波输出, MILO工作点一般设计在基模(TM00模)的π模, 但器件运行初期仍然容易有不同程度的模式竞争现象出现[11-21]. 由于MILO中波束互作用区通常采用同轴盘荷波导作为其慢波结构, 因此, 从同轴盘荷波导中几个可能被相对论电子束激发的低阶模式(TM00, TM01和HEM11模)与电子束之间的色散关系入手, 简要分析了三种类型的模式竞争的特点、产生的可能原因以及削弱方法, 具体包括基模TM00模相邻纵向模式之间的竞争、高阶对称模TM01与基模之间的竞争以及非对称模式HEM11与基模之间的竞争, 为磁绝缘线振荡器运行过程中出现的竞争模式的识别和输出性能的优化提供参考. 在此基础上, 给出了一种稳定紧凑型L波段MILO的物理模型并利用全电磁三维粒子程序进行了数值模拟, 重点对起振阶段的模式竞争问题进行了分析.

    • 由于MILO中波束互作用区通常采用同轴盘荷波导作为其慢波结构, 下面从同轴盘荷波导中几个可能被电子束激发的低阶模式与电子束之间的色散关系入手, 对MILO中模式竞争问题进行物理分析.

      同轴盘荷波导的纵向中心截面(即过轴线的纵向截面, 以下简称纵向截面)和横截面如图1所示, 包括阴极、阳极叶片和阳极外壳, 它们均为金属导体, 其中p为慢波结构的周期长度, 整个结构为轴对称的圆柱结构.

      图  1  同轴盘荷波导的结构示意图 (a) 纵向中心截面; (b) 横截面

      Figure 1.  Schematic drawings of the coaxial disk-loaded waveguide in (a) Axial and (b) radial cross section

      利用电磁场模拟程序, 计算了上述同轴盘荷波导中TM00模、TM01模和HEM11模的色散曲线, 具体计算方法可参考文献[4,20,22], 其中阴极半径为4.0 cm, 阳极叶片和阳极外壳的内半径分别为6.0 cm和10.0 cm, 周期长度p为3.0 cm, 电子束电压为500 kV, 具体结果如图2所示. 图2中还给出了电子束和真空中光束的色散曲线, ${v_{\rm{e}}}$${v_{{\rm{drift}}}}$分别为电子可能达到的最大速度和达到磁绝缘时的漂移速度. 图2A, BC所对应的点分别为电子束与TM00模, HEM11模和TM01模色散曲线的交点. MILO工作点一般设计在基模(TM00模)的π模, 即点A处, 对应的纵向电场在纵向和横向截面分布的等高图如图3(a)图3(b)所示.

      图  2  同轴盘荷波导中TM00模、TM01模、HEM11模以及电子束和光束的色散关系图

      Figure 2.  The uncoupled dispersion curves of the coaxial disk-loaded waveguide (TM00, TM01 and HEM11), and the dispersion lines of the electron and light beam.

      图  3  同轴盘荷波导中TM00(π模)的轴向电场在(a) 纵向中心截面和(b) 横截面的等高图

      Figure 3.  The contour plots of Ez of the TM00 (π mode) in the coaxial disk-loaded waveguide in the (a) Axial and (b) radial cross section.

    • 对于来自基模(TM00)非π模的竞争, 由图2可以看到, 可能参与竞争的纵向模式一般是与π模(点A)较接近的模式, 如(N–1)π/N, N为互作用区阳极腔的个数, 即慢波结构的周期数. 由于是基模之间的竞争, 因此竞争模式的电磁场分布仍然具有轴对称性(图3(b)), 共振频率略低于工作点π模对应的频率. 该竞争模式初期可能先于π模被激发, 一段时间后, 可能被π模抑制.

      已经知道, 在MILO中, 从阴极发射的电子束的能量分配到了径向和轴向两个方向. MILO设计时, 要求电子的平均漂移速度约等于基模TM00模的相速度. 但由于电子的动能与径向位置相关, 因此, 磁绝缘尚未完全形成前, 电子速度沿径向连续变化, 电子束色散线可扫过如图2的阴影部分[15,16]. 电子束色散线与TM00模的色散曲线将产生一系列交点, 从而可能激发起TM00模的含π模在内的其他纵向模式. TM00模的纵向模式之间的竞争已经在实验上观察到且被证实[17,18].

      通过合理选取阳极腔的个数来控制模式隔离度, 或者选取合适的阳极叶片半径, 尤其是提取腔叶片的半径来调整束波互作用腔的品质因数, 可以削弱或者抑制TM00模的纵向模式之间的竞争.

    • 对于基模与高阶模TM01之间的竞争, 此时竞争模式TM01的电磁场分布也具有轴对称性, 只是轴向电场沿径向正负交替一次, 如图4所示. 由色散关系图2可以看到, 点C为电子束和TM01模色散曲线的交点, 其对应的共振频率远高于点A对应的频率, 且交点位置随电子束纵向速度(电压)变化相对较为敏感, 因此激发难度相对较大. 但由于此时束波共振点对应的本征模的群速度小于零, 如果器件品质因数较高或者互作用区较长, 器件运行初期有可能被激发.

      图  4  同轴盘荷波导中TM01的轴向电场在横截面的等高图

      Figure 4.  The contour plots of Ez of the TM01 mode in the coaxial disk-loaded waveguide in the radial cross section.

      通过调节阳极叶片的半径进而降低束波互作用腔的品质因数, 或者通过增加负载电流强度缩短MILO达到磁绝缘状态的时间, 都可能有效抑制TM01模式的起振和放大.

    • 对于非对称模式HEM11模与基模的竞争, 由图2可以看到, HEM11各纵向模式的频率仅略高于与之对应的TM00模的频率, π模附近的频率相差最小(图中点A和点B). 其中, HEM11模在点B处对应的纵向电场在纵向和横向截面分布的等高图如图5(a)图5(b)所示. 可以看到, HEM11模的轴向电场Ez不再具有轴对称性, 而是沿角向正负交替一次.

      图  5  同轴盘荷波导中HEM11(π模)的轴向电场在 (a) 纵向中心截面和(b)横截面的等高图

      Figure 5.  The contour plots of Ez of the HEM11(π mode)in the coaxial disk-loaded waveguide in the (a) Axial and (b) radial cross section.

      已有研究表明: 任何小的结构非对称或者电子发射的非均匀, 都可能在器件运行初期激发HEM11模, 且起振较快, 在器件运行初期很难完全抑制[19-21]. 由于MILO慢波结构包括反射区、互作用区和提取输出区, 在这些区域, 阳极叶片半径各有差异, 加上输出区还有支撑杆, 导致与TM00模的频率接近的HEM11模容易被激发

      尽量保证器件结构的对称性和阴极电子发射的均匀性, 尽管不能够完全消除这种竞争模式, 但可以最大程度地缩短模式竞争的时间和减小竞争模式的强度.

    • 为了获得稳定的高功率微波输出, 基于以上分析和认识, 结合已有工作基础[3-6,13,14], 设计了一种紧凑型L波段MILO, 其结构示意图如图6所示, 除支撑杆沿角向以90°夹角均匀排列外, 整个结构为轴对称结构, 其中波束互作用区采用的慢波结构为前面部分讨论过的同轴盘荷波导, 点D为后面的模拟计算中在互作用区中设置的观察点. 下面我们对上述MILO中的模式竞争问题和其输出特性进行模拟和分析.

      图  6  紧凑型MILO纵向截面示意图

      Figure 6.  Schematic drawing of the compact MILO in the axial cross section.

    • 为了计算波束互作用区中的共振频率, 将上述MILO的互作用区两端封闭, 同时沿z方向在第四和第五叶片之间加载一个短脉冲宽频信号来激励电磁场振荡. 图7给出了上述激励信号在MILO互作用区中点D的位置激发出的电场频谱图. 由图7可以看到: 在容易被激发的频率中, 较低的两个频率分别为1.61和1.77 GHz.

      图  7  短脉冲宽频信号激发后MILO互作用区中点D处电场的频谱图

      Figure 7.  The oscillation frequency of the Ez stimulated in cold cavity at point D.

      图8(a)图8(b)图9(a)图9(b)分别给出了频率分别为1.61和1.77 GHz时对应的纵向电场在XZXY两个截面分布的等高图. 结合前面的分析, 可知这两个频率对应的电磁场模式分别是TM00模和HEM11模的π模. 由此可见, 与基模TM00模相同, 对于HEM11模, 与电子束容易产生相互作用的模式也是π模.

      图  8  频率为1.61 GHz的纵向电场在 (a) 纵向和(b) 横向截面等高图

      Figure 8.  The contour plots of Ez with the frequency of 1.61 GHz in (a) The axial and (b) the radial cross section.

      图  9  频率为1.77 GHz的纵向电场在 (a) 纵向和(b) 横向截面的等高图

      Figure 9.  The contour plots of Ez with the frequency of 1.77 GHz in (a) The axial and (b) the radial cross section.

    • 下面, 我们利用全电磁三维粒子模拟程序对给出的如图6所示的L波段MILO的输出性能进行数值模拟, 同时考察初始阶段的模式竞争情况, 其中注入的电子束的平均电压约为510 kV, 对应的电功率约为32.8 GW.

      图10(a)(d)图11(a)(d)给出了阳极腔内D点电场和磁场随时间的变化曲线及针对不同时段的FFT变换图. 由图10(a)和11(a)可以看到, 辐射场8.0 ns左右开始起振, 大约22.0 ns达到饱和, 其中在8.0至14.0 ns之间, 存在一个小的峰值. 对电磁场按时间分段做FFT变换, 可以得到, 全时段(0—50 ns)的电磁场的中心频率约为1.6 GHz, 但是可以看到有略高于该频率的小峰值, 提示有一定强度的竞争模式; 初期(8—14 ns)的电磁场的频率约为1.7 GHz, 稳定后(> 14 ns), 电磁场的频率约为1.6 GHz, 频谱比较纯, 略高于该频率的小峰值没有出现. 可见, 在起振后的初始阶段, 存在明显的模式竞争; 大约过10 ns左右, 竞争模式被抑制.

      图  10  (a) 阳极腔内观察点D的电场Ez随时间的变化曲线及不同时间窗口(b) 0−50 ns, (c) 8−14 ns和(d) 14−50 ns的FFT变换图

      Figure 10.  (a) Variations of Ez with time at the observation point D in the anode cavity; the corresponding Fourier transform with different time intervals: (b) 0−50 ns; (c) 8−14 ns; (d) 14−50 ns.

      图  11  (a) 阳极腔内观察点D磁场By随时间的变化曲线及不同时间窗口(b) 0−50 ns, (c) 8−14 ns和(d) 14−50 ns的FFT变换图

      Figure 11.  (a) Variations of By with time at the observation point D in the anode cavity, and the corresponding Fourier transform with different time intervals: (b) 0−50 ns, (c) 8−14 ns (d) 14−50 ns.

      图12(a)(d)给出了初始阶段不同时刻互作用区D点所在横截面轴向电场的等高图. 可以看到, 轴向电场沿角向的分布由初期的非轴对称分布, 逐渐转变为轴对称分布. 由轴向电场的分布特征和共振频率, 结合冷腔计算结果, 可以确定, 在初始阶段, 非对称模式HEM11模与TM00模产生竞争; 由于加载了电子束, 因此与冷腔分析时的结果进行比较, TM00模和HEM11模的频率由冷腔时的1.61和1.77 GHz分别降为1.6和1.7 GHz.

      图  12  初始阶段不同时刻互作用区横截面轴向电场的等高图 (a) t = 12.171 ns; (b) t = 13.094 ns; (c) t = 14.944 ns; (d) t = 16.793 ns

      Figure 12.  Contour plots of Ez in the interaction region at different times at the initial stage.

      图13(a)图13(b)给出了饱和时互作用区中轴向电场在纵向和横向截面的等高图, 可以看到, 饱和时轴向电场在互作用区的阳极腔中正负交替排列(π模场分布), 沿角向呈轴对称分布, 与TM00模的场分布完全一致, 因此, 可以确定器件稳定后互作用区的工作模式为TM00模的π模.

      图  13  饱和时互作用区轴向电场在(a) 纵向截面和(b) 横向截面的等高图, t = 45.727 ns

      Figure 13.  Contour plots of Ez in the (a) Axial and (b) radial cross section in the interaction region at saturation.

      由于是在直角坐标系中建立的三维计算模型, 加上输出区还有支撑杆, 因此互作用区很难实现完全轴对称. 同时, 由于阴极电子发射面大, 电子发射也存在一定的非均匀性, 这些因素使得器件运行初期容易激发HEM11模, 且很难完全抑制. 起振一小段时间后(10 ns左右), 互作用区基模增长加快, 高阶模被抑制.

      基于以上分析, 为了提高器件的输出功率和效率, 对器件的参数进行了进一步的优化, 主要调节了互作用区的周期长度和电子束电压, 其他参数保持不变. 模拟结果表明, 当周期长度p为3.1 cm, 注入的电子束的平均电压约为600 kV, 对应电功率约为45 GW时, 器件的输出效率达到了18%.

      图14(a)图14(b)给出了MILO输出周期平均功率随时间的变化曲线以及稳定后输出口电场的FFT变换图. 由图14可以看到, 饱和时输出功率高达8.1 GW, 输出频谱比较纯, 中心频率约为1.552 GHz. 由于MILO的输出区为同轴波导, 输出模式应为同轴TEM模. 图15(a)图15(b)给出了输出口横截面电磁场分布的箭矢图, 可以看到, 计算所得的电磁场分布确与理想同轴TEM模的分布一致. 利用文献[23]给出的方法, 计算得到了模式纯度约为97%.

      图  14  (a) MILO输出周期平均功率随时间的变化曲线; (b)稳定后输出口电场的FFT变换图

      Figure 14.  (a) Time plots of the periodic-average output power of the MILO; (b) the Fourier transform of Ex at the output port.

      图  15  电磁场在输出口横向截面分布的箭矢图 (a) 电场; (b)磁场

      Figure 15.  The vector plots of the (a) Electric field and (b) magnetic field in the radial cross section at the output port.

    • 从同轴盘荷波导中几种可能被电子束激发的低阶模式 (TM00, TM01和HEM11模)的色散曲线与电子束色散曲线关系入手, 简要分析了三种类型的模式竞争的特点和产生的可能原因以及削弱方法, 具体包括基模TM00模相邻纵向模式之间的竞争、高阶对称模TM01与基模之间的竞争以及非对称模HEM11与基模之间的竞争. 在此基础上, 给出了一种稳定紧凑型L波段磁绝缘线振荡器的物理模型, 并利用全电磁三维粒子程序对其中的模式竞争问题进行了数值模拟. 冷腔分析表明: 在容易被激发的本征模式中, 对应频率较低的模式分别是TM00和HEM11的π模, 对应的频率分别为1.61 GHz和1.77 GHz. 热腔的模拟结果表明, 该器件运行初期主要竞争模式是HEM11模, 起振约十个ns后, 竞争模式基本被抑制. 进一步优化后器件在基模获得了高效率和高功率微波输出, 饱和时输出功率约8.1 GW, 输出效率达到了18%, 模式纯度约为97%. 研究结果可为磁绝缘线振荡器运行过程中出现的竞争模式的识别和输出性能优化提供理论参考和依据. 下一步将在控制好模式竞争的基础上改进和优化器件结构, 以进一步提高器件的输出功率和效率.

参考文献 (23)

目录

    /

    返回文章
    返回