开槽波导3次谐波回旋行波放大管非线性理论与数值模拟
开槽波导3次谐波回旋行波放大管非线性理论与数值模拟
本文讨论了开槽圆柱波导的高频场分布，给出了注波互作用自洽非线性理论.在电子作大回旋运动与考虑速度零散的情况下，采用四阶龙格库塔法，对均匀截面开槽波导3次谐波回旋行波放大管注波互作用进行了数值计算，得出一些重要的互作用规律，为回旋行波放大管的进一步研究打下了基础.
关键词:回旋行波放大管；开槽波导；自洽非线性；高次谐波；速度零散
self-consistent nonlinear theory and simulation of a slotted third-harmonic gyro-twt amplifier
zhang hong-bin,li hong-fu,zhou xiao-lan,wang hua-jun,yu sheng,du pin-zhong
(inst.of high energy electronics,uest of china,chengdu 610054,china)
abstract:the distribution of rf field of the slotted cylindrincal wave guide is discussed and the self-consistent nonlinear theory of the beam-wave interaction is presented in this paper.the behavior of the slotted gyrotron travelling-wave amplifier (gyro-twt) with a uniform section is simulated by a runge-kutta algorithm code for a warm beam encircling around the axis of the wave guide.some important regulations are obtained.this work presents the bases to further studies of the gyro-twt.
key words:gyro-twt;slotted wave guide;self-consistent nonlinear;high harmonic wave;velocity spread
一、引言
回旋行波放大管属于毫米波放大器件，它以高功率、高效率、宽频带而著称，在雷达与通讯等领域有着极其重要的应用前景，自七十年代末以来，在理论和实验方面都取得了长足的进展［1～5］.
对于基次谐波回旋行波管，在毫米波波段需要很高的直流磁场，因而需要体积较大的超导系统或电磁铁系统来提供直流磁场.采用高次谐波互作用，便可大大降低管子对直流磁场的要求［2，3］，使采用永久磁铁成为可能，从而可大大减小管子的体积.由于开槽壁和光滑壁波导中高频场分布存在的差异，开槽波导更有利于注波互作用，对工作电压要求较低，工作效率比光滑壁波导要高，同时与光滑壁波导相比具有很好的模式竞争抑制能力［6］.本文以95ghz开槽3次谐波为例，对回旋行波放大管进行了数值模拟，得到了一些重要的互作用规律.
二、高频场模式和特性
图1所示为开槽波导结构以及电子注轨迹横截面图(虚圆表示电注横截面图).设n为开槽波导的槽数，θ0为间隙半张角，a、b分别为波导内外半径，r、φ、z为电子的柱坐标，v⊥为电子的横向速度，φ为动量空间角，即v⊥与x轴夹角.为了方便起见，将波导分为两个区域进行讨论，即：ⅰ区(0＜r＜a)和ⅱ区(a＜r＜b).由于在回旋行波管中电子注与波的有效互作用场为te波场，故仅需关心横电波高频场的分布情况［7～9］.这里只给出了高频电场分量的表达式，有关高频磁场分量的表达式可进一步能过电磁场分量关系求得.
图1 中空外开槽波导及电子注横截面示意图.虚圆为电子注横截面示意图
在ⅰ区(0＜r＜a)中
(1)
(2)
ez=0 (3)
在ⅱ区(a＜r＜b)中
ez=0 (4)
er=0 (5)
(6)
其中
(7)
b0=［-j′0(kcb)/y′0(kcb)］a0 (8)
在以上各式中，e0为高频场振幅，γ为角向谐波数，αγ为角向γ次谐波项的振幅系数，kc为截止波数，q为开槽序数(q=1,2,…,n)，m代表高频场的角向模式(m=0,1,2,…,n-1).aγ的值以及电路的色散关系可由电磁场在r=a处的边界条件确定.
(9)
色散关系为
(10)
式(9)表明，只有当空间谐波次数γ=m+ln时，非零空间谐波项才存在.角向模式决定相邻隙间高频场的相位差，对于每一具体模式，此相位差值为m2π/n.每一角向模式均由无数个角向谐波项组成，其谐波振幅系数由式(9)决定.在所有角向模式中有两个比较重要的模式，即π模式和2π模式，其角向谐波相对强弱分布情况见图2.由图2可知，2π模式的能量主要集中于零次谐波项中，而π模式的能量主要集中于±n/2次谐波项中.因此，π模式较2π模式更适合于高次回旋谐波互作用.如果电子注回旋谐波次数(用s表示)已经设定，那么槽数n的选择应保证最强非零次角向谐波项的次数г与回旋谐波次数s相等.如，对于π模式，槽数n应等于2s.
图2 角向谐波振幅对角向谐波数(γ)的相对分布示意图.(a)π模式(m=n/2,n=6,θ0=15°)，(b)2π模式(m=0,n=6,θ0=15°)
当角向模式m和槽深(即a/b的值)确定后，截止波数kc的值可由式(10)通过数值求解方法求得［6，8，9］.
三、自洽非线性理论
在热腔中，高频场沿轴向呈缓变分布状况，其对横坐标(r,φ)的分布函数与冷腔情况相同.下面给出ⅰ区中的热腔高频电场分量(te波)表达式.
(11)
(12)
ez=0 (13)
上述各式中，cmn为电场归一化系数，f(z)为一复函数，代表高频场沿z轴的缓变分布情况.cmn的值由下式求得
(14)
以下是自洽非线性注波互作用常微分方程组.
从洛伦兹公式
出发［8］，可推得电子在高频场(e,b)和直流磁场(b0)作用下的运动方程.每个电子有6个运动参量方程，这里仅给出了速度分量及动量空间角3个运动参量方程.
(15)
(16)
(17)
以上各式中，m0和γ分别为电子的静止质量和相对论因子，φ为动量空间角，u=γv,v为电子的速度，如图1所示.
从有源麦克斯韦方程出发，经过一系列复杂的推导并对电流进行离散化后得到非线性注波互作用场方程为
(18)
上式中，p为在一个高频场周期内所取的电子注批数，m为考虑电子注厚度因数而将电子注化分的圈数，n为每圈上所取的宏电子数，s为谐波次数.〈…〉表示对初始速度分布函数为g0(v⊥,vz)的速度空间进行平均.设电子注为单能电子注，速度零散主要来自于横纵向速度比值(v⊥/vz)的零散，这里按正态分布规律来处理速度零散，即初始速度分布函数为
式中k为归一化常数,△vz为平均纵向速度零散，δ为狄拉克函数．
边界条件
f(z)|z=0=f(0) (19)
(20)
式中f(0)为输入高频场电场幅值.
方程(15)～(18)为自洽非线性注波互作用方程组.将电子注离散为nt个宏电子，则一共有6nt+2个一阶非线性微分方程，结合边界条件(19)、(20)，利用四阶龙格库塔法对注波互作用进行数值计算，计算结果在下部分内容中给出并讨论.
四、结果与讨论
表1给出了互作用电路参数，各图表曲线相关参数见相应图表标注.图3给出了驱动功率为20w情况下，效率与电子速度比值α的关系.图中b0、bg分别为直流磁场和共振点磁场，ω为高频场频率，ωc为波导截止频率.由于在回旋行波管中波的能量取自于电子的横向能，又由于当α值增大，电子的横向能量以及回旋半径也随着增大，因此互作用效率也就随着α增大而增大.但当α增大到一定值后，注波互作用达到饱和，同时由于电子注回旋半径过大，电子在波导壁上产生截获，这样互作用效率又随α值增大而减小.
表1 数值模拟参数与结果
内半径 1.024mm
外半径 1.465mm
电路长度 87.9mm
注电压 60kv
注电流 6a
α 1.3
直流磁场 11.674kg
高频场模式 π
谐波次数 3
工作频率 95.08ghz
模拟结果
饱和效率 22.8%
饱和输出功率 82kw
饱和增益 36.15db
图3 效率与电子注速度比值α的关系(s=3,πmode,i=6a,v=60kv,ω/ωc=1.032，
b0/bg=0.99)
图4所示为饱和效率、饱和增益与b0/bg值之间的关系，虚线为饱和增益曲线.图中γz为纵向速度分量的相对论因子.图示表明，一方面，降低b0/bg值，有助于提高饱和互作用效率，但b0/bg值不能太低，否则失谐加重，注波互作用难以达到同步，饱和效率便会迅速降低；另一方面，增加b0/bg的值却有利于提高饱和增益.总的来说，磁场失谐率的选择应在效率和增益之间作优化折衷.
图4 饱和效率及增益与b0/bg值的关系(s=3,π mode,i=6a,v=60kv,ω/ωc=γz,
α=1.3)
图5所示电流分别为3a、6a和9a情况下(a)饱和效率、(b)饱和增益随频率变化的关系.可以看出饱和效率、饱和增益以及饱和带宽都随电流的增长而有所增加.在6a和图示情况下，饱和带宽为7%，电流为3a增大到9a时，饱和带宽从4.6%增大到8.3%.
图5 不同电流下，(a)饱和效率(b)饱和增益随频率变化的关系(s=3,π mode,v=60kv,α=1.3,b0/bg=0.99)
图6所示为几个不同磁场失谐率下饱和增益以及饱和效率随频率变化的关系.由图可见，磁场失谐率对饱和增益、饱和效率及饱和带宽都有较大影响，b0/bg值的提高有利于饱和增益及饱和带宽的提高，但饱和效率却有所降低.在图示条件下，当b0/bg值从0.983增大到0.998时，饱和带宽从4.8%增大到9.3%.
图6 不同磁场失谐率下，(a)饱和增益及(b)饱和效率随频率变化的关系(s=3,π mode,i=6a,v=60kv,α=1.3)
图7为在不同磁场失谐率下饱和效率随谐波次数的变化关系.由图表明，饱和效率随谐波次数的增大而降低，b0/bg值越低，谐波次数对饱和效率的影响越大.
图7 饱和效率随谐波次数的变化关系(π mode,i=6a,v=60kv,α=1.3,ω/ωc=γz,rl/a=0.7)
图8所示为不同谐波次数下饱和效率随频率的变化关系.图示表明谐波次数对饱和带宽有较大影响.在图示条件下，谐波次数从2增大到4时，饱和带宽从10.3%减小到5.7%.
图8 不同谐波次数下饱和效率随频率的变化关系(π mode,i=6a,v=60kv,α=1.3,b0/bg=0.99,ω/ωc=γz,rl/a=0.7)
五、结束语
本文在单模和未考虑空间电荷效应及波导壁损耗的情况下，对95ghz3次谐波回旋行波管注波互作用进行了数值计算，得出了一些重要的互作用规律.研究结果表明，开槽波导高次谐波回旋行波管能在较低的磁场下和较宽的频带范围内，获得较高的互作用效应.高次谐波互作用降低了管子对磁场的要求，使采用永久磁铁成为可能，但谐波次数的增高会削弱注波互作用效率和带宽.由于在π模式中，高频场能量主要集中于高次谐波项中，而在2π模式中，能量主要集中于零次谐波项中，因此，π模式较2π模式更有利于高次谐波放大.适当降低磁场失谐率b0/bg的值，有利于提高饱和互作用效率，但饱和增益及带宽却有所降低.适当提高横纵向速度比值(v⊥/vz)、电流值以及输入功率，有利于带宽、增益和互作用效应的提高.另外，速度零散对注波互作用亦有较大影响，随着速度零散的增加，注波互作用效率、增益都有所降低.
选择红外热像仪需要考虑那几点因素
英特尔发布全新视觉技术加速器解决方案 提升边缘设备深度学习和人工智能性能
当通信网络引入人工智能技术会怎么样
崇达技术宣布将收购普诺威35%的股权 将扩展进入集成电路领域培育新的业务增长点
中易云智慧能源管理平台可实现什么样的功能
开槽波导3次谐波回旋行波放大管非线性理论与数值模拟
宿舍安装人脸识别系统的优缺点是什么
如何对电动汽车动力电池散热方法在这！
什么是特锐祥贴片Y电容？一起来了解下！
直流变交流逆变器的工作原理及电路分享
比克如何应对补贴政策调整
程序员应该怎样保护自己的隐私
Shure宣布升级全线SE耳机
可调电源diy教程
C语言字符串转数字实现方法
电子芯闻早报：华为世界五百强排名上升99位
【回顾往年CES】LG展出Ultra HD 4K画质的OLED屏幕
运动蓝牙耳机选购指南：运动耳机什么款式比较好用
便携式氨气检测仪性能介绍
大部分VR开发者，在用哪个平台的硬件？
本文讨论了开槽圆柱波导的高频场分布，给出了注波互作用自洽非线性理论.在电子作大回旋运动与考虑速度零散的情况下，采用四阶龙格库塔法，对均匀截面开槽波导3次谐波回旋行波放大管注波互作用进行了数值计算，得出一些重要的互作用规律，为回旋行波放大管的进一步研究打下了基础.
关键词:回旋行波放大管；开槽波导；自洽非线性；高次谐波；速度零散
self-consistent nonlinear theory and simulation of a slotted third-harmonic gyro-twt amplifier
zhang hong-bin,li hong-fu,zhou xiao-lan,wang hua-jun,yu sheng,du pin-zhong
(inst.of high energy electronics,uest of china,chengdu 610054,china)
abstract:the distribution of rf field of the slotted cylindrincal wave guide is discussed and the self-consistent nonlinear theory of the beam-wave interaction is presented in this paper.the behavior of the slotted gyrotron travelling-wave amplifier (gyro-twt) with a uniform section is simulated by a runge-kutta algorithm code for a warm beam encircling around the axis of the wave guide.some important regulations are obtained.this work presents the bases to further studies of the gyro-twt.
key words:gyro-twt;slotted wave guide;self-consistent nonlinear;high harmonic wave;velocity spread
一、引言
回旋行波放大管属于毫米波放大器件，它以高功率、高效率、宽频带而著称，在雷达与通讯等领域有着极其重要的应用前景，自七十年代末以来，在理论和实验方面都取得了长足的进展［1～5］.
对于基次谐波回旋行波管，在毫米波波段需要很高的直流磁场，因而需要体积较大的超导系统或电磁铁系统来提供直流磁场.采用高次谐波互作用，便可大大降低管子对直流磁场的要求［2，3］，使采用永久磁铁成为可能，从而可大大减小管子的体积.由于开槽壁和光滑壁波导中高频场分布存在的差异，开槽波导更有利于注波互作用，对工作电压要求较低，工作效率比光滑壁波导要高，同时与光滑壁波导相比具有很好的模式竞争抑制能力［6］.本文以95ghz开槽3次谐波为例，对回旋行波放大管进行了数值模拟，得到了一些重要的互作用规律.
二、高频场模式和特性
图1所示为开槽波导结构以及电子注轨迹横截面图(虚圆表示电注横截面图).设n为开槽波导的槽数，θ0为间隙半张角，a、b分别为波导内外半径，r、φ、z为电子的柱坐标，v⊥为电子的横向速度，φ为动量空间角，即v⊥与x轴夹角.为了方便起见，将波导分为两个区域进行讨论，即：ⅰ区(0＜r＜a)和ⅱ区(a＜r＜b).由于在回旋行波管中电子注与波的有效互作用场为te波场，故仅需关心横电波高频场的分布情况［7～9］.这里只给出了高频电场分量的表达式，有关高频磁场分量的表达式可进一步能过电磁场分量关系求得.
图1 中空外开槽波导及电子注横截面示意图.虚圆为电子注横截面示意图
在ⅰ区(0＜r＜a)中
(1)
(2)
ez=0 (3)
在ⅱ区(a＜r＜b)中
ez=0 (4)
er=0 (5)
(6)
其中
(7)
b0=［-j′0(kcb)/y′0(kcb)］a0 (8)
在以上各式中，e0为高频场振幅，γ为角向谐波数，αγ为角向γ次谐波项的振幅系数，kc为截止波数，q为开槽序数(q=1,2,…,n)，m代表高频场的角向模式(m=0,1,2,…,n-1).aγ的值以及电路的色散关系可由电磁场在r=a处的边界条件确定.
(9)
色散关系为
(10)
式(9)表明，只有当空间谐波次数γ=m+ln时，非零空间谐波项才存在.角向模式决定相邻隙间高频场的相位差，对于每一具体模式，此相位差值为m2π/n.每一角向模式均由无数个角向谐波项组成，其谐波振幅系数由式(9)决定.在所有角向模式中有两个比较重要的模式，即π模式和2π模式，其角向谐波相对强弱分布情况见图2.由图2可知，2π模式的能量主要集中于零次谐波项中，而π模式的能量主要集中于±n/2次谐波项中.因此，π模式较2π模式更适合于高次回旋谐波互作用.如果电子注回旋谐波次数(用s表示)已经设定，那么槽数n的选择应保证最强非零次角向谐波项的次数г与回旋谐波次数s相等.如，对于π模式，槽数n应等于2s.
图2 角向谐波振幅对角向谐波数(γ)的相对分布示意图.(a)π模式(m=n/2,n=6,θ0=15°)，(b)2π模式(m=0,n=6,θ0=15°)
当角向模式m和槽深(即a/b的值)确定后，截止波数kc的值可由式(10)通过数值求解方法求得［6，8，9］.
三、自洽非线性理论
在热腔中，高频场沿轴向呈缓变分布状况，其对横坐标(r,φ)的分布函数与冷腔情况相同.下面给出ⅰ区中的热腔高频电场分量(te波)表达式.
(11)
(12)
ez=0 (13)
上述各式中，cmn为电场归一化系数，f(z)为一复函数，代表高频场沿z轴的缓变分布情况.cmn的值由下式求得
(14)
以下是自洽非线性注波互作用常微分方程组.
从洛伦兹公式
出发［8］，可推得电子在高频场(e,b)和直流磁场(b0)作用下的运动方程.每个电子有6个运动参量方程，这里仅给出了速度分量及动量空间角3个运动参量方程.
(15)
(16)
(17)
以上各式中，m0和γ分别为电子的静止质量和相对论因子，φ为动量空间角，u=γv,v为电子的速度，如图1所示.
从有源麦克斯韦方程出发，经过一系列复杂的推导并对电流进行离散化后得到非线性注波互作用场方程为
(18)
上式中，p为在一个高频场周期内所取的电子注批数，m为考虑电子注厚度因数而将电子注化分的圈数，n为每圈上所取的宏电子数，s为谐波次数.〈…〉表示对初始速度分布函数为g0(v⊥,vz)的速度空间进行平均.设电子注为单能电子注，速度零散主要来自于横纵向速度比值(v⊥/vz)的零散，这里按正态分布规律来处理速度零散，即初始速度分布函数为
式中k为归一化常数,△vz为平均纵向速度零散，δ为狄拉克函数．
边界条件
f(z)|z=0=f(0) (19)
(20)
式中f(0)为输入高频场电场幅值.
方程(15)～(18)为自洽非线性注波互作用方程组.将电子注离散为nt个宏电子，则一共有6nt+2个一阶非线性微分方程，结合边界条件(19)、(20)，利用四阶龙格库塔法对注波互作用进行数值计算，计算结果在下部分内容中给出并讨论.
四、结果与讨论
表1给出了互作用电路参数，各图表曲线相关参数见相应图表标注.图3给出了驱动功率为20w情况下，效率与电子速度比值α的关系.图中b0、bg分别为直流磁场和共振点磁场，ω为高频场频率，ωc为波导截止频率.由于在回旋行波管中波的能量取自于电子的横向能，又由于当α值增大，电子的横向能量以及回旋半径也随着增大，因此互作用效率也就随着α增大而增大.但当α增大到一定值后，注波互作用达到饱和，同时由于电子注回旋半径过大，电子在波导壁上产生截获，这样互作用效率又随α值增大而减小.
表1 数值模拟参数与结果
内半径 1.024mm
外半径 1.465mm
电路长度 87.9mm
注电压 60kv
注电流 6a
α 1.3
直流磁场 11.674kg
高频场模式 π
谐波次数 3
工作频率 95.08ghz
模拟结果
饱和效率 22.8%
饱和输出功率 82kw
饱和增益 36.15db
图3 效率与电子注速度比值α的关系(s=3,πmode,i=6a,v=60kv,ω/ωc=1.032，
b0/bg=0.99)
图4所示为饱和效率、饱和增益与b0/bg值之间的关系，虚线为饱和增益曲线.图中γz为纵向速度分量的相对论因子.图示表明，一方面，降低b0/bg值，有助于提高饱和互作用效率，但b0/bg值不能太低，否则失谐加重，注波互作用难以达到同步，饱和效率便会迅速降低；另一方面，增加b0/bg的值却有利于提高饱和增益.总的来说，磁场失谐率的选择应在效率和增益之间作优化折衷.
图4 饱和效率及增益与b0/bg值的关系(s=3,π mode,i=6a,v=60kv,ω/ωc=γz,
α=1.3)
图5所示电流分别为3a、6a和9a情况下(a)饱和效率、(b)饱和增益随频率变化的关系.可以看出饱和效率、饱和增益以及饱和带宽都随电流的增长而有所增加.在6a和图示情况下，饱和带宽为7%，电流为3a增大到9a时，饱和带宽从4.6%增大到8.3%.
图5 不同电流下，(a)饱和效率(b)饱和增益随频率变化的关系(s=3,π mode,v=60kv,α=1.3,b0/bg=0.99)
图6所示为几个不同磁场失谐率下饱和增益以及饱和效率随频率变化的关系.由图可见，磁场失谐率对饱和增益、饱和效率及饱和带宽都有较大影响，b0/bg值的提高有利于饱和增益及饱和带宽的提高，但饱和效率却有所降低.在图示条件下，当b0/bg值从0.983增大到0.998时，饱和带宽从4.8%增大到9.3%.
图6 不同磁场失谐率下，(a)饱和增益及(b)饱和效率随频率变化的关系(s=3,π mode,i=6a,v=60kv,α=1.3)
图7为在不同磁场失谐率下饱和效率随谐波次数的变化关系.由图表明，饱和效率随谐波次数的增大而降低，b0/bg值越低，谐波次数对饱和效率的影响越大.
图7 饱和效率随谐波次数的变化关系(π mode,i=6a,v=60kv,α=1.3,ω/ωc=γz,rl/a=0.7)
图8所示为不同谐波次数下饱和效率随频率的变化关系.图示表明谐波次数对饱和带宽有较大影响.在图示条件下，谐波次数从2增大到4时，饱和带宽从10.3%减小到5.7%.
图8 不同谐波次数下饱和效率随频率的变化关系(π mode,i=6a,v=60kv,α=1.3,b0/bg=0.99,ω/ωc=γz,rl/a=0.7)
五、结束语
本文在单模和未考虑空间电荷效应及波导壁损耗的情况下，对95ghz3次谐波回旋行波管注波互作用进行了数值计算，得出了一些重要的互作用规律.研究结果表明，开槽波导高次谐波回旋行波管能在较低的磁场下和较宽的频带范围内，获得较高的互作用效应.高次谐波互作用降低了管子对磁场的要求，使采用永久磁铁成为可能，但谐波次数的增高会削弱注波互作用效率和带宽.由于在π模式中，高频场能量主要集中于高次谐波项中，而在2π模式中，能量主要集中于零次谐波项中，因此，π模式较2π模式更有利于高次谐波放大.适当降低磁场失谐率b0/bg的值，有利于提高饱和互作用效率，但饱和增益及带宽却有所降低.适当提高横纵向速度比值(v⊥/vz)、电流值以及输入功率，有利于带宽、增益和互作用效应的提高.另外，速度零散对注波互作用亦有较大影响，随着速度零散的增加，注波互作用效率、增益都有所降低.
选择红外热像仪需要考虑那几点因素
英特尔发布全新视觉技术加速器解决方案 提升边缘设备深度学习和人工智能性能
当通信网络引入人工智能技术会怎么样
崇达技术宣布将收购普诺威35%的股权 将扩展进入集成电路领域培育新的业务增长点
中易云智慧能源管理平台可实现什么样的功能
开槽波导3次谐波回旋行波放大管非线性理论与数值模拟
宿舍安装人脸识别系统的优缺点是什么
如何对电动汽车动力电池散热方法在这！
什么是特锐祥贴片Y电容？一起来了解下！
直流变交流逆变器的工作原理及电路分享
比克如何应对补贴政策调整
程序员应该怎样保护自己的隐私
Shure宣布升级全线SE耳机
可调电源diy教程
C语言字符串转数字实现方法
电子芯闻早报：华为世界五百强排名上升99位
【回顾往年CES】LG展出Ultra HD 4K画质的OLED屏幕
运动蓝牙耳机选购指南：运动耳机什么款式比较好用
便携式氨气检测仪性能介绍
大部分VR开发者，在用哪个平台的硬件？