-
回旋振荡管产生高功率毫米波、亚毫米波,是受控核聚变等离子体加热和电流驱动,是等离子体诊断、太空垃圾处理、高密度数据通信、核磁共振波谱、高分辨率雷达和拒止武器等应用领域不可替代的高效电磁辐射源[1-4]。目前,100~200 GHz频段的兆瓦级回旋振荡管[1, 5]凭借在电子回旋谐振加热和电流驱动等方面的出色表现成为磁约束受控核聚变加热系统的首选波源,受到越来越多的专家学者的关注。随着工作频率和功率的进一步提高,回旋振荡管需要高次模工作,这将使工作模式与竞争模式间分割度减小,带来强烈的模式竞争,降低回旋振荡管的工作效率和稳定性;高电压、大电流工作,受空间电荷效应和枪区自激振荡等因素的影响,电子光学系统的主要指标参数零散进一步加剧,影响整管的工作性能。因此,有效抑制模式竞争和设计高质量电子注是兆瓦级回旋振荡管面临的突出问题。探索高性能、有效抑制非工作模式的谐振腔,设计高质量、大功率电子枪是开展核聚变用兆瓦级回旋振荡管研究的重要内容[6-7]。
本文以基于耦合波理论的回旋管自洽非线性理论为基础,编制回旋振荡管冷、热腔模拟代码,以此为工具,通过考察TE模式谱、注-波耦合系数和工作模式及其附近模式竞争模式的起振电流,设计出了有效抑制非工作模式、高效互作用的170 GHz兆瓦级回旋振荡管。此外,建立电子注速度零散模型,研究了速度零散对该只回旋振荡管工作性能的影响,为设计满足互作用需要的高质量电子光学系统提供理论参考。
-
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]:
$$ \left\{ \begin{array}{l} {{\mathit{\boldsymbol{E}}}_t}({\mathit{\boldsymbol{r}}}, \phi, z) = \sum\limits_{mn} {V_{mn}^{(1)}(z){\mathit{\boldsymbol{e}}}_{mn}^{(1)}({\mathit{\boldsymbol{r}}}, \phi )} + \sum\limits_{mn} {V_{mn}^{(2)}(z){\mathit{\boldsymbol{e}}}_{mn}^{(2)}} ({\mathit{\boldsymbol{r}}}, \phi ) \hfill \\ {{\mathit{\boldsymbol{H}}}_t}({\mathit{\boldsymbol{r}}}, \phi, z) = \sum\limits_{mn} {I_{mn}^{(1)}(z){\mathit{\boldsymbol{h}}}_{mn}^{(1)}({\mathit{\boldsymbol{r}}}, \phi )} + \sum\limits_{mn} {I_{mn}^{(2)}(z){\mathit{\boldsymbol{h}}}_{mn}^{(2)}({\mathit{\boldsymbol{r}}}, \phi )} \hfill \\ \end{array} \right. $$ (1) 式中,$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)的正交归一化矢量波函数,定义为:
$$\left\{ \begin{array}{l} {\mathit{\boldsymbol{e}}}_{_{mn}}^{(1)}(r, \phi ) = - {\nabla _t}\prod\nolimits_{mn}^{(1)} {} \hfill \\ {\mathit{\boldsymbol{h}}}_{_{mn}}^{(1)}(r, \phi ) = {\nabla _t}\Pi _{_{mn}}^{(1)} \times {{\bf{i}}_z} = {{\bf{i}}_z} \times {\mathit{\boldsymbol{e}}}_{_{mn}}^{(1)}(r, \phi ) \hfill \\ {\mathit{\boldsymbol{e}}}_{_{mn}}^{(2)}(r, \phi ) = {{\bf{i}}_z} \times {\nabla _t}\prod\nolimits_{mn}^{(2)} {} \hfill \\ {\mathit{\boldsymbol{h}}}_{_{mn}}^{(2)}(r, \phi ) = - {\nabla _t}\prod\nolimits_{mn}^{(2)} = {{\bf{i}}_z} \times {\mathit{\boldsymbol{e}}}_{_{mn}}^{(2)}(r, \phi ) \hfill \\ \end{array} \right. $$ (2) 式中,$\prod\nolimits_{mn}^{(i)} {} $是赫兹位函数,且满足:
$$ \left\{ \begin{array}{l} \nabla _{^ \bot }^2\prod\nolimits_{mn}^{(i)} + k_c^2\prod\nolimits_{mn}^{(i)} { = 0} \hfill \\ {\left. {\prod\nolimits_{mn}^{(1)} {} } \right|_s} = 0 \hfill \\ {\mathit{\boldsymbol{n}}} \cdot \nabla {\left. {\prod\nolimits_{mn}^{(2)} {} } \right|_s} = 0 \hfill \\ \end{array} \right. $$ (3) 把式(1)代入电、磁场满足的旋度方程:
$$ \left\{ \begin{array}{l} \nabla \times {\mathit{\boldsymbol{E}}} = - \frac{{\partial {\mathit{\boldsymbol{B}}}}}{{\partial t}} \hfill \\ \nabla \times {\mathit{\boldsymbol{H}}} = {\mathit{\boldsymbol{J}}} + \frac{{\partial {\mathit{\boldsymbol{D}}}}}{{\partial t}} \hfill \\ \end{array} \right. $$ (4) 整理后得到可以描述回旋振荡管中注-波互作用的有源广义传输线方程组:
$$ \left\{ \begin{array}{l} \frac{{\text{d}}}{{{\text{d}}z}}V_{mn}^{(i)} = - {\text{j}}\omega \mu I_{mn}^{(i)} + \frac{{\text{j}}}{{\omega \varepsilon }}\iint\limits_s {{\nabla _t}{J_z} \cdot {\mathit{\boldsymbol{e}}}_{mn}^{(i)}}{\text{d}}s + \sum\limits_{mp} {V_{mp}^{(i)}C_{pn}^{i, i'}\frac{1}{a}\frac{{{\text{d}}a}}{{{\text{d}}z}}} \hfill \\ \frac{{\text{d}}}{{{\text{d}}z}}I_{mn}^{(i)} = \frac{{\text{j}}}{{\omega \mu }}\gamma _{mn}^{(i)2}V_{mn}^{(i)} - \iint\limits_s {{{\mathit{\boldsymbol{J}}}_t}} \cdot {\mathit{\boldsymbol{e}}}_{mn}^{(i)}{\text{d}}s - \sum\limits_{mp} {I_{mp}^{(i)}C_{np}^{i, i'}\frac{1}{a}\frac{{{\text{d}}a}}{{{\text{d}}z}}} \hfill \\ \end{array} \right. $$ (5) 式中,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}$置零,上式即演变为描述回旋振荡管谐振腔特征参数的无源广义传输线方程组。数值求解无源广义传输线方程组,可得到某一频点不同模式间的相对幅值,为设计出有效抑制非工作模式的谐振腔提供理论依据。
$$ \iint\limits_s {{\mathit{\boldsymbol{e}}}_{mp}^{(j)}} \cdot \frac{{\partial {\mathit{\boldsymbol{e}}}_{mn}^{(i) * }}}{{\partial z}}{\text{d}}s = C_{pn}^{i, j}\frac{1}{a}\frac{{{\text{d}}a}}{{{\text{d}}z}} \;\;\;\; i, j=1, 2 $$ (6) 则,对于模式TEmn、TEmp、TMmn、TMmp有:
$$ C_{mn, mp}^{1, 1} = \left\{ \begin{array}{l} - \frac{{{\varepsilon _m}}}{2}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; (p = n) \hfill \\ \frac{{\gamma _{mn}^2{\varepsilon _m}}}{{\gamma _{mp}^2 - \gamma _{mn}^2}}\;{( - 1)^{p + n}}\;\;{\text{ }}(p \ne n) \hfill \\ \end{array} \right. $$ (7) $$ C_{mn, mp}^{1, 2} = 0 $$ (8) $$ C_{mn, mp}^{2, 1} = \frac{{m{\varepsilon _m}}}{{\sqrt {\mu _{mp}^2 - {m^2}} }}{( - 1)^{p + n}} $$ (9) $$ C_{mn, mp}^{2, 2} = \left\{ \begin{array}{l} - \frac{{{\varepsilon _m}}}{2}\frac{{{m^2}}}{{\mu _{mn}^2 - {m^2}}}\;\; (p = n) \hfill \\ - \frac{{{\varepsilon _m}\mu _{mn}^2}}{{\mu _{mn}^2 - \mu _{mp}^2}} \hfill \\ \times \frac{{{{( - 1)}^{p + n}}({m^2} - \mu _{mn}^2)}}{{\sqrt {\mu _{mn}^2 - {m^2}} \sqrt {\mu _{mp}^2 - {m^2}} }}{\text{ }}(p \ne n)\; \hfill \\ \end{array} \right. $$ (10) 式(5)中各模式的激励源来自有相对论效应的回旋电子注,电子在电、磁场中的运动方程为:
$$ \frac{{{\text{d}}{\mathit{\boldsymbol{v}}}}}{{{\text{d}}t}} = - \frac{e}{{{m_0}\gamma }}({\mathit{\boldsymbol{E}}} + {\mathit{\boldsymbol{v}}} \times {\mathit{\boldsymbol{B}}} - \frac{{\mathit{\boldsymbol{v}}}}{{{c^2}}}{\mathit{\boldsymbol{v}}} \cdot {\mathit{\boldsymbol{E}}}) $$ (11) 式中,${\mathit{\boldsymbol{v}}}$是电子的运动速度;${m_0}$为电子的静止质量;$\gamma $为相对论因子。
式(5)中各模式在谐振腔的左、右端口分别满足消失波和行波条件:
$$ \sum\limits_{i = 1}^2 {\sum\limits_{m, n} {{{\left| {\frac{{{\text{d}}V_{mn}^{(i)}}}{{{\text{d}}z}} \pm \gamma _{mn}^{(i)}V_{mn}^{(i)}} \right|}_{Z = 0, L}}} = 0} $$ (12) 式(5)、式(11)和式(12)形成以传输线形式描述的回旋振荡管自洽非线性理论。
由起振电流估算公式[11]:
$$ \begin{array}{l} \frac{{ - 1}}{{{I_{st}}}} = {\left( {\frac{\pi }{\lambda }{{\int_0^L {\left| {\overline f (z)} \right|} }^2}dz} \right)^{ - 1}}\left( {s + \frac{1}{2}\frac{{\omega {\beta ^2}_{ \bot 0}}}{{{v_{z0}}}}\frac{\partial }{{\partial \Delta s}}} \right)\frac{{{L^2}\pi }}{4}{{\text{e}}^{\left( { - \frac{{{{({\Delta _s}L)}^2}}}{8}} \right)}} \cdot \\ \left( {\frac{{Q{Z_0}e}}{{8{\gamma _0}{m_e}{c^2}}}} \right){\left( {\frac{{{k_{mp}}{C_{mp}}{G_{mp}}}}{{{\beta _{z0}}(s - 1)!}}} \right)^2}{\left( {\frac{{c{k_{mp}}{\gamma _0}{\beta _{ \bot 0}}}}{{2{\Omega _0}}}} \right)^{2(s - 1)}} \\ \end{array} $$ (13) 式中,${R_b}$表示电子注引导中心半径;$\overline f (z)$是谐振腔中TE模式的纵向场幅值分布函数;L是谐振腔参与互作用的有效长度;s是谐波次数;${\gamma _0}, \lambda, {\Omega _0}, Q, {Z_0}$分别代表电子的相对论因子初值,模式波长,非相对论回旋频率,品质因素和真空中的波阻抗。进一步考察工作模式附近的竞争模式分布,为注-波互作用提供合适的引导磁场参数,验证谐振腔参数(品质因素)和电子注引导中心半径设置是否合适。其中,电子注引导中心半径[4]通过考察注波耦合系数来分析确定。
$$ {C_{BF}} = \frac{{\mu _{mn}^2J_{m \pm s}^2({\mu _{mn}}{R_b}/a)}}{{\pi a_{}^2(\mu _{mn}^2 - {m^2})J_m^2({\mu _{mn}})}} $$ (14) -
采用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的TE32, 12模式),证实了电子回旋脉塞理论。
-
电子注参数的一致性对回旋振荡管的工作性能会产生一定的影响。本文基于统计学规律,将单一电子注能量时的速度分布函数(正态分布)引入回旋振荡管注-波互作用模型,模拟研究了不同电子注速度零散下的注-波互作用效率。该模型中电子的速度函数分布为:
$$ {f_e}({v_ \bot }) = \frac{1}{{\sigma \sqrt {2\pi } }}\exp \left[{-\frac{{{{({v_ \bot }-{v_{ \bot 0}})}^2}}}{{2{\sigma ^2}}}} \right] $$ (15) 式中,${v_ \bot }$是电子的真实速度;${v_{ \bot {\text{0}}}}$是电子的平均速度;$\sigma $是方差。并定义速度零散:
$$ \delta {v_ \bot } = \frac{1}{2}\frac{{{v_{ \bot, \max }} - {v_{ \bot, \min }}}}{{{v_{ \bot, {\text{center}}}}}} $$ (16) 则互作用效率随速度零散的变化情况如图 8所示。可见,速度零散在10%范围内,回旋管的性能几乎没有变化,但速度零散大于20%时,互作用效率下降约5%。此结果和参考文献[13]利用MAGY软件考察速度零散时的得出的结论有较好的一致性。这将为下一步电子光学系统设计提供理论依据。
Design and Study of 170 GHz Gyrotron with MW Power
-
摘要: 该文基于广义传输线理论和电子回旋脉塞理论,考察模式谱、深入研究注-波耦合系数和不同模式的起振电流,选择TE32,12模式为170 GHz兆瓦级回旋振荡管的工作模式。通过合理调整电子注工作参数有效抑制高次模谐振腔中的暂态模式,优化谐振腔的光滑缓变段来抑制寄生振荡模式。当电子注的横、纵速度比为1.3,引导中心半径为9.4 mm,电压为80 kV,电流为45 A,互作用磁场约为6.715 T时,该只优化设计的回旋振荡管在170.18 GHz频点输出功率超过1.8 MW,效率约为50%。此外,建立切合实际的电子注速度零散模型,讨论了速度零散对该只管子的性能影响。Abstract: Based on the theories of generalized transmission line and electron cyclotron resonance maser, the TE32, 12 mode is selected as the operating mode of a 170 GHz megawatt gyrotron by investigating the mode spectrum and studying the beam-wave coupling coefficient and the start-current of the modes carefully. Through appropriately arranging the parameters of the beam, the transient modes in the high order mode resonator are suppressed effectively. And through designing the gradual sections of a kind of taped cavity, the parasitic modes are suppressed furthest. As a result, an output power of 1.8 MW, corresponding to 50% efficiency, and an oscillation frequency of 170.18 GHz have been achieved with a pitch factor 1.3, radius of guiding center 9.4mm, 80 kV, 45A helical electron beam at a guiding magnetic field of 6.715T. At the same time, we have also set up a velocity-spread model suitable to actual circumstances and discussed the performance impact to the gyrotron.
-
Key words:
- efficiency /
- gyrotron /
- mode competition /
- starting current
-
[1] THUMM M, ARNOLD A, BORIE E, et al. Frequency step-tunable (114-170 GHz) megawatt gyrotrons for physics application[J]. Fusion Engineering and Design, 2001, 53(1-4):407-421. doi: 10.1016/S0920-3796(00)00519-6 [2] DAMMERTZ G, ALBERTI S, ARNOLD A, et al. Development of multimegawatt gyrotrons for fusion plasma heating and current drive[J]. IEEE Trans on Electron Devices, 2005, 52(2):808-817. http://ieeexplore.ieee.org/document/1424363/ [3] PRATER R. Heating and current driving by electron cyclotron waves[J]. Physics of Plasmas, 2004, 11(5):2349- 2376. doi: 10.1063/1.1690762 [4] NOTAKE T, SAITO T, TATEMATSU Y, et al. Development of a novel high power sub-THz second harmonic gyrotron[J]. Physical Review Letters, 2009, 103(22):225002. doi: 10.1103/PhysRevLett.103.225002 [5] HILLESHEIM J C, DEBOO J C, PEEBLES W A, et al. Observation of a critical gradient threshold for electron temperature fluctuations in the DⅢ-D Tokamak[J]. Physical Review Letters, 2013, 110(4):045003. doi: 10.1103/PhysRevLett.110.045003 [6] THUMM M K. Recent Advances in the worldwide fusion gyrotron development[J]. IEEE Trans on Plasma Scinece, 2014, 42(3):590-599. doi: 10.1109/TPS.2013.2284026 [7] ESPOSITO B, GRANUCCI G, SMEULDERS P, et al. Disruption avoidance in the Frascati Tokamak upgrade by means of magneto hydrodynamic mode stabilization using electron-cyclotron-resonannce heating[J]. Physical Review Letters, 2008, 100(4):045006(4). http://www.ncbi.nlm.nih.gov/pubmed/18352292 [8] LI Hong-fu, THUMM M. Mode coupling in corrugated waveguides with varying wall impedance and diameter change[J]. International Journal of Electronics, 1991, 71(5):827-844. doi: 10.1080/00207219108925527 [9] 雷朝军, 喻胜, 李宏福, 等.缓变回旋管谐振腔研究[J].物理学报, 2012, 61(18):180202. doi: 10.7498/aps.61.180202 LEI Chao-jun, YU Sheng, LI Hong-fu, et al. Study on gradually varying cavity for a gyrotron[J]. Acta Phys Sin, 2012, 61(18):180202. doi: 10.7498/aps.61.180202 [10] HUANG Y, LI H F, YANG S W. Study of a 35 GHz third-harmonic low-voltage complex cavity gyrotron[J]. IEEE Trans Plasma Science, 1999, 27(2):368-373. doi: 10.1109/27.772263 [11] KARTIKEYAN M V, BORE E, THUMM M K A. Gyrotrons high power microwave and millimeter wave technology[M]. New York, USA:Springer Press, 2004. [12] 雷朝军, 刘迎辉, 喻胜, 等.高次模回旋管中的非工作模式抑制研究[J].电子科技大学学报, 2015, 38(1):79- 82. http://manu50.magtech.com.cn/dzkjdx/CN/abstract/abstract161.shtml LEI Chao-jun, LIU Ying-hui, YU Sheng, et al. Study on the mode suppression of a high-order mode Gyrotron[J]. Journal of University of Electronic Science and Technology of China, 2015, 38(1):79-82. http://manu50.magtech.com.cn/dzkjdx/CN/abstract/abstract161.shtml [13] PU R F, NUSINOVICH G S, SINITSYN O V, et al. Numercial study of efficiency for a 670 GHz gyrotron[J]. Physics of Plasmas, 2011, 18:023107. doi: 10.1063/1.3554656