170 GHz兆瓦级回旋振荡管设计研究

雷朝军，刘迎辉，张双狮，蒙 林，李宏福

(1. 中国人民警察大学 河北 廊坊 065000；2. 电子科技大学电子科学与工程学院太赫兹科学技术研究中心 成都 610054)

回旋振荡管产生高功率毫米波、亚毫米波，是受控核聚变等离子体加热和电流驱动，是等离子体诊断、太空垃圾处理、高密度数据通信、核磁共振波谱、高分辨率雷达和拒止武器等应用领域不可替代的高效电磁辐射源[1-4]。目前，100～200 GHz频段的兆瓦级回旋振荡管[1,5]凭借在电子回旋谐振加热和电流驱动等方面的出色表现成为磁约束受控核聚变加热系统的首选波源，受到越来越多的专家学者的关注。随着工作频率和功率的进一步提高，回旋振荡管需要高次模工作，这将使工作模式与竞争模式间分割度减小，带来强烈的模式竞争，降低回旋振荡管的工作效率和稳定性；高电压、大电流工作，受空间电荷效应和枪区自激振荡等因素的影响，电子光学系统的主要指标参数零散进一步加剧，影响整管的工作性能。因此，有效抑制模式竞争和设计高质量电子注是兆瓦级回旋振荡管面临的突出问题。探索高性能、有效抑制非工作模式的谐振腔，设计高质量、大功率电子枪是开展核聚变用兆瓦级回旋振荡管研究的重要内容[6-7]。

本文以基于耦合波理论的回旋管自洽非线性理论为基础，编制回旋振荡管冷、热腔模拟代码，以此为工具，通过考察TE模式谱、注-波耦合系数和工作模式及其附近模式竞争模式的起振电流，设计出了有效抑制非工作模式、高效互作用的170 GHz兆瓦级回旋振荡管。此外，建立电子注速度零散模型，研究了速度零散对该只回旋振荡管工作性能的影响，为设计满足互作用需要的高质量电子光学系统提供理论参考。

1 理论与方法

HFSS、CST等专业电磁仿真软件可以研究谐振腔的特征参数，然而这些专业仿真软件模拟得到的是某一频率时谐振腔内的总场分布，无法准确获得不同寄生模式间的相对幅值；很难优化设计出能有效抑制非工作模式的回旋振荡管谐振腔。为了解决这一问题，可以采用模式耦合理论，将谐振腔内的总场用正交模式展开，以各模式幅值满足的广义传输线方程组(即广义传输线理论)来研究谐振腔，准确得到各寄生模式的与工作模式间的相对幅值大小。针对不同的寄生模式特点，优化设计回旋振荡管谐振腔，则可以得到工作模式纯度较高的回旋管高频结构。

回旋振荡管采用柱形开放式谐振腔，柱坐标系下，将谐振腔中的电、磁场分解为横向电场Et、横向磁场Ht和纵向场电场Ezez、纵向磁场Hzez；横向电场 Et(z )、磁场 Ht( z)是纵向位置z的函数，并用规则圆波导中的矢量波函数展开[8-10]：

把式(1)代入电、磁场满足的旋度方程：

整理后得到可以描述回旋振荡管中注-波互作用的有源广义传输线方程组：

式中，i(i′)可取1和2，取1时式(5)是电波方程组，取2时式(5)是磁波方程组；a代表谐振腔内壁半径；、分别代表电、磁波的传播常数；J、J分tz别代表回旋电子注的横、纵向电流密度；将式(5)中Jt和Jz置零，上式即演变为描述回旋振荡管谐振腔特征参数的无源广义传输线方程组。数值求解无源广义传输线方程组，可得到某一频点不同模式间的相对幅值，为设计出有效抑制非工作模式的谐振腔提供理论依据。

根据模式间的耦合系数的定义[5-7]：

则，对于模式TEmn、TEmp、TMmn、TMmp有：

式(5)中各模式的激励源来自有相对论效应的回旋电子注，电子在电、磁场中的运动方程为：

式中，v是电子的运动速度；m0为电子的静止质量；γ为相对论因子。

式(5)中各模式在谐振腔的左、右端口分别满足消失波和行波条件：

式(5)、式(11)和式(12)形成以传输线形式描述的回旋振荡管自洽非线性理论。

由起振电流估算公式[11]：

式中，Rb表示电子注引导中心半径；(z)是谐振腔中TE模式的纵向场幅值分布函数；L是谐振腔参与互作用的有效长度；s是谐波次数；γ0, λ,Ω0, Q ,Z0分别代表电子的相对论因子初值，模式波长，非相对论回旋频率，品质因素和真空中的波阻抗。进一步考察工作模式附近的竞争模式分布，为注-波互作用提供合适的引导磁场参数，验证谐振腔参数(品质因素)和电子注引导中心半径设置是否合适。其中，电子注引导中心半径[4]通过考察注波耦合系数来分析确定。

2 模式选择和冷腔设计

兆瓦级回旋振荡管一般采用水冷式冷却系统，在目前的技术水平下，腔壁损耗一般不超过3～4 kW/cm2，这就需要170 GHz回旋振荡管尽可能采用高阶模式工作。高次模回旋振荡管中的模式分割度减小，模式竞争严重影响着其工作性能和稳定性，需要从工作模式选择、起振电流、热腔工作参数、不同时段的模式起振过程等方面综合考虑，深入研究各模式间的耦合系数和注-波耦合系数，有效抑制暂态模式和寄生模式，保证回旋振荡管高效率稳定工作[9,12]。

2.1 工作模式选择

回旋振荡管一般选用磁波作为工作模式，其高能电子注由磁控注入式电子枪提供。根据绝热压缩理论和动量守恒定律，在磁压缩比和工作电流一定时，电子注在互作用区的引导中心半径越大，阴极半径就会越大，同电流发射密度的情况下，发射带则越窄，越容易设计出满足要求的高质量回旋电子注；谐振腔半径越大，模式在腔壁的单位面积上欧姆损耗越小，越容易实现的兆瓦级功率输出。可见，选择高阶模作为工作模式能更利于大功率回旋振荡管设计。然而，高次模工作，参与注-波互作用的模式数量增多，工作模式与非工作模式分隔度变小，加大了模式竞争风险。精心选择工作模式也是大功率回旋振荡管设计需要考虑的主要内容之一。

综合以上因素，考察图1中TE模式谱，并考察图2中TE32,12模式及其附近的竞争模式的横截面电场分布。显然，TE47,7模式可提供较大的电子注引导中心半径，但场能在腔壁附近较集中，壁损耗过大；TE17,18、TE6,23模式会引起电子注引导半径过小，增加电子注设计难度。TE32,12模式可兼得兆瓦级回旋振荡管的要求。本文选择最有代表性的TE32,12模式作为170 GHz兆瓦级回旋振荡管的工作模式。

图1 与工作模式相关的部分TE模式谱

图2 TE32,12模式及其附近模式的横向场分布

2.2 注-波耦合系数

根据式(14)，可以确定TE32,12模式及其附近各模式的注-波耦合系数分布。为了使图像更清晰，只取最可能与TE32,12模式有强模式竞争的4个模式考察其注-波耦合系数分布，如图3所示。可见，注半径约为9.4 mm时，工作模式有较强的注-波耦合，并能有效抑制其他非工作模式。

图3 相关模式的注-波耦合系数随半径的变化曲线

2.3 起振电流

图4 相关模式的起振电流与磁场的关系曲线

考察模式的起振电流分布，可以大致确定回旋振荡管的参数范围，并能了解竞争模式的参数分布。由式(13)得到在注电压为80 kV，注半径为9.4 mm，速度比为1.3，品质因素为1 100时，170 GHz频率附近相关模式的引导磁场与电流分布如图4所示。磁场在6.77～6.93 T区间时如图4a所示，TE32,12模式可优先起振(Istart＜40 A)。潜在的竞争模式在此磁场区间的起振电流均大于55 A，如图4b所示。前部分模式谱分析中最可能存在模式竞争的TE47,7(Istart＞108A)和TE17,18(Istart＞ 280 A)，由于注-波耦合系数在引导中心半径处较小，很难起振。本文设计的高功率回旋振荡管工作电流为45 A，可保证在起振阶段有效抑制暂态模式竞争。

2.4 高频结构

改变腔体缓变段结构[9]，可以有效抑制寄生模式。在大量模拟计算的基础上，优化设计出了腔内工作模式纯度达67 dB的光滑缓变结构谐振腔，其沿轴向缓变的半径轮廓线如图5所示；腔内工作模式TE32,12和最大寄生模式TE32,11幅值的轴向分布曲线也在图5中示出。

图5 缓变腔中相关模式幅值沿轴向分布

3 热腔分析

3.1 自洽非线性理论计算

图6 互作用效率沿轴向分布

采用4阶龙格库塔法求解第一部分中提到的由传输线理论描述的回旋管自洽非线性理论，得到第二部分优化设计的缓变腔回旋振荡管工作参数。在注电压为80 kV，注电流为45 A，横纵速度比为1.3，磁场为6.715 T时，互作用效率约为50%，如图6所示。同时可以看到，起振电流中的磁场范围和热腔分析有所不同，约有0.06 T的偏差，这主要是由于起振电流是用高斯场近似真实模式幅值分布求解的。显然，从图5冷腔场分布可以看出，这种近似是比较粗略的。

3.2 电子横向群聚图

回旋电子注在注-波互作用过程中的相空间演变过程如图7所示。图7a是某圈宏电子的初始相空间分布(Z=0，角向均匀分布)。图7b是该圈宏电子在经历一段时间注-波互作用后的相空间分布(Z=900Δh，出现一个明显的群聚中心，其中Δh为数值求解步长)。图7c是该圈宏电子在注-波互作用末期的相空间分布(Z=2 530Δh，绝大部分宏电子回旋半径明显减小)。可知，在注-波互作用末期，大量回旋电子将能量交给了预先设计的高频场(170 GHz的TE32,12模式)，证实了电子回旋脉塞理论。

图7 某圈宏电子在不同位置的相空间分布

3.3 速度零散

电子注参数的一致性对回旋振荡管的工作性能会产生一定的影响。本文基于统计学规律，将单一电子注能量时的速度分布函数(正态分布)引入回旋振荡管注-波互作用模型，模拟研究了不同电子注速度零散下的注-波互作用效率。该模型中电子的速度函数分布为：

式中，v⊥是电子的真实速度；v⊥0是电子的平均速度；σ是方差。并定义速度零散：

则互作用效率随速度零散的变化情况如图8所示。可见，速度零散在10%范围内，回旋管的性能几乎没有变化，但速度零散大于20%时，互作用效率下降约5%。此结果和参考文献[13]利用MAGY软件考察速度零散时的得出的结论有较好的一致性。这将为下一步电子光学系统设计提供理论依据。

图8 效率随速度零散的变化

4 结 束 语

本文基于广义传输线理论，在编写“冷”“热”腔程序的基础上，通过合理选择工作模式，考察注-波耦合系数和起振电流，确定电子注大致参数范围；通过优化缓变谐振腔结构，提高工作模式纯度，使寄生模式得到有效抑制。通过合理设置电子注引导中心半径和工作磁场，最大可能地避免了高次模回旋振荡管复杂的模式竞争。经过优化设计，得到了互作用效率约50%，输出功率大于1.8 MW的170 GHz兆瓦级回旋振荡管。此外，结合工程实际，建立电子注速度零散模型，考察了这只回旋振荡管在考虑电子注参数零散时的工作性能。本文将为研制170 GHz兆瓦级回旋振荡管提供设计指导和理论依据。