张治强,王克强,胡 标*,张庆元
(1. 西北核技术研究所高功率微波技术重点实验室 西安 710024;2. 电子科技大学电子科学与工程学院 成都 611731)
模式变换器是高功率微波链路系统中最为常见的导波器件,在军事、科研、民用等领域得到广泛使用[1-5]。自20 世纪70 年代以来,基于规则波导的变模技术获得了发展,国际上提出了波导内各模式传输特性的耦合波理论[6]。对于常规矩形波导、圆波导等规则截面波导而言,相应的耦合波方程也容易推导[7-8]。然而,针对诸如椭圆波导等异型截面波导,由于公式中的赫兹函数难以求出,给耦合波方程组的应用带来极大困难。为了解决此类问题,本文尝试运用矢量波形函数来推导非规则波导内的耦合波方程组。矢量波形函数可以直接根据时域有限元法进行求解,所以对于任何非规则截面波导都可以采用这种方法。本文初步设计给出的X 波段椭圆TM01-TE11波导模式变换器结果表明,变模效率高于95%的带宽为10%,最高转换效率高于99%,较经典的圆波导变模器件有更好的性能表现,充分说明了基于矢量波函数耦合波理论在非规则截面波导模式转换技术中应用的有效性。
为了分析弯曲波导中的场,首先引用无源麦克斯韦方程组的复量形式:
在弯曲波导中,电磁场可以分解为横向分量和纵向分量,即:
式中,t 代表横向截面分量; z代表纵向传播分量。作为场方程组(1)和纵向单位矢量iw的矢积,根据环形坐标系中的关系,可化简为:
为了验证推导的矢量波函数耦合波理论的正确性,本文设计一种高功率微波传输系统中重要的波导器件TM01-TE11模式变换器,其主要作用是将高功率微波源产生的不利于传输和辐射的TM01模式转换为适合传输与辐射的TE11模式[9-12]。
首先,利用时域有限元法对给定椭圆波导(中心频率为9 GHz,椭圆长半轴和短半轴长度分别为30 mm 和20 mm)内的本征模进行求解,得到各个模式的矢量波形函数;然后,利用编写好的计算程序和前面推导的公式对椭圆波导内耦合波方程组进行数值求解;再利用迭代算法对椭圆TM01-TE11模式变换器进行合成设计[13]。计算结果显示,当波导弧长为0.35 m 时,转换效率可以达到99%以上,达到迭代收敛要求。整个模式转换器整体为三弯曲对称互易结构,输入和输出端口平行且中心点在同一水平直线,便于系统前后级对接。为了清晰地展示波导内各个模式的耦合过程,TM01-TE11模式变换器的结构及各个模式沿纵向的能量变化情况如图1 所示,从图中可以看出TM01逐渐转变为TE11模式。另外,在波导前半段和后半段有少量能量分别转换为TM11和TE21模式,其主要原因是TM01-TM11和TE11-TE21模式对之间的耦合能力较强。
图1 TM01-TE11 模式变换器各个模式沿纵向的能量变化情况
将数值计算合成的椭圆TM01-TE11模式变换器在CST 内进行建模与仿真,电磁波在TM01-TE11模式变换器内部将沿椭圆波导传播方向的电场分布及输入/输出端口处的模式电场分布如图2 所示。显然,在输入端口注入的椭圆TM01模式经过变换器后在输出端口转换为椭圆TE11模。
图2 椭圆TM01-TE11 模式变换器结构及端口和传播方向处的电场分布图
通常,在高功率微波应用系统中,模式转换器在真空状态下(小于10-2Pa)的峰值击穿场强阈值通常取为1×108V/m。因此,根据式(15)可以计算得到变换器在各个频段的功率容量:
式中,Emax为真空条件下微波最大击穿场强阈值;E1W为当输入功率为1W 时的峰值场强。从图中还可看出,当注入功率为0.5 W 时,弯头内最大峰值电场强度为1 048 V/m,折合计算出该模式变换器的最大功率容量达到4.55 GW。
分别采用数值计算和CST 仿真的方法对椭圆变换器的带宽特性进行计算和比较,结果如图3所示。通过对比可以发现,该变换器在中心频点的转换效率达到99.5%,95%以上转换效率带宽超过0.9 GHz。另外,在靠近中心频点两侧,数值与仿真结果存在一定误差,这主要是由于偏离中心点后数值计算精度降低所致。但总体上数值与仿真结果具有良好的一致性,验证了上述对非规则截面波导内耦合波方程组推导的正确性。
图3 数值计算和CST 仿真的带宽特性比较
基于赫兹矢量函数的耦合方程及其表达形式是描述圆波导、方波导等规则结构的弯曲波导模式变换的经典理论。随着应用的不断扩展,此类变模技术已无法满足多频点、宽频带和紧凑性等技术要求。本文重新推导并给出了基于采用矢量波形函数的耦合波理论表达形式,结合数值求解的方法可以处理任意非规则截面轴向弯曲波导内的电磁特性,具有更强的求解能力和更广的应用范围。设计的X 波段椭圆TM01-TE11波导模式变换器证明了本文方案的有效性,为推动非规则截面变模技术的发展提出了新的努力方向。