Design and Study of 170 GHz Gyrotron with MW Power

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

回旋振荡管采用柱形开放式谐振腔，柱坐标系下，将谐振腔中的电、磁场分解为横向电场${{\mathit{\boldsymbol{E}}}_t}$、横向磁场${{\mathit{\boldsymbol{H}}}_t}$和纵向场电场${E_z}{{\mathit{\boldsymbol{e}}}_z}$、纵向磁场${H_z}{{\mathit{\boldsymbol{e}}}_z}$；横向电场${{\mathit{\boldsymbol{E}}}_t}(z)$、磁场${{\mathit{\boldsymbol{H}}}_t}(z)$是纵向位置$z$的函数，并用规则圆波导中的矢量波函数展开^[8-10]：

式中，$V_{mn}^{(1)}(z)$、$I_{mn}^{(1)}(z)$分别是电波的“电压”“电流”幅度函数；$V_{mn}^{(2)}(z)$、$I_{mn}^{(2)}(z)$分别为磁波的“电压”“电流”幅度函数，${\mathit{\boldsymbol{e}}}_{mn}^{(1)}({\mathit{\boldsymbol{r}}}, \phi )$、${\mathit{\boldsymbol{e}}}_{mn}^{(2)}({\mathit{\boldsymbol{r}}}, \phi )$分别代表电波(TM)、磁波(TE)的正交归一化矢量波函数，定义为：

式中，$\prod\nolimits_{mn}^{(i)} {} $是赫兹位函数，且满足：

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

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

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

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

则，对于模式TE_mn、TE_mp、TM_mn、TM_mp有：

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

式中，${\mathit{\boldsymbol{v}}}$是电子的运动速度；${m_0}$为电子的静止质量；$\gamma $为相对论因子。

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

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

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

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

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

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

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

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

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

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

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

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

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

式中，${v_ \bot }$是电子的真实速度；${v_{ \bot {\text{0}}}}$是电子的平均速度；$\sigma $是方差。并定义速度零散：

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

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