基于混合二维有限元/Leontovich方法的非理想导体平面波导结广义散射矩阵精确建模

《IEEE Transactions on Microwave Theory and Techniques》：Complete Modal Analysis of Planar Waveguide Junctions With a 2-D Hybrid FEM/Leontovich Framework Including Conductor Losses

本文针对毫米波和太赫兹频段平面波导结中导体损耗严重影响器件性能的问题，提出了一种混合二维有限元/Leontovich边界条件方法，用于精确计算包含侧壁和盖板损耗的广义散射矩阵。该方法通过变分公式和模态展开，高效地集成了导体损耗效应，为高性能无源器件的设计提供了关键工具，相关成果发表于《IEEE Transactions on Microwave Theory and Techniques》。

　　

在毫米波和太赫兹技术飞速发展的今天，高性能平面电路与波导系统成为推动通信、传感和成像等领域进步的核心。然而，随着频率升高至毫米波甚至太赫兹波段，导体的趋肤效应加剧，导致波导壁的欧姆损耗显著增加，这严重制约了滤波器、耦合器、功分器等关键无源器件的性能，如插入损耗恶化、品质因数下降。传统上，在分析此类平面波导结（如微带线、鳍线、基片集成波导中的不连续性结构）时，通常假设导体为理想电导体（PEC），这在高频情况下会引入不可忽视的误差，无法准确预测电路的实际性能。因此，发展一种能够精确计入导体损耗的电磁建模方法，对于设计高性能毫米波/太赫兹器件至关重要。

为了应对这一挑战，研究人员在《IEEE Transactions on Microwave Theory and Techniques》上发表论文，致力于解决非理想导体壁条件下平面波导结的电磁建模问题。本研究的核心目标是建立一种高效的数值方法，用于精确计算包含损耗的波导结的广义散射矩阵（GSM），从而为电路设计提供可靠依据。

本研究采用的关键技术方法主要包括：1）混合二维有限元法（2-D FEM）与模态分析：利用二维FEM处理波导结横向截面的场分布，结合解析的模态展开描述波导端口的激励与响应，将三维问题降维为二维问题，显著提升计算效率。2）Leontovich（或阻抗）边界条件（IBC）的应用：通过引入表面阻抗Z_s来模拟非理想导体（如金属）的损耗特性，将其作为边界条件融入变分公式中，分别处理波导侧壁和上下盖板的损耗效应。3）变分公式的构建与求解：针对横磁波（TM^z）和横电波（TE^z）两种基本激励模式，分别推导了包含损耗项的弱形式波动方程，并通过Galerkin方法进行离散化，最终形成可用于求解GSM的线性系统。该方法避免了复杂的三维体积网格剖分，适用于分析任意形状的平面波导结。

TM^z模激励下的二维FEM/Leontovich分析

研究人员首先分析了由TM^z模（磁场纵向分量为零）激励的波导结。通过将电场表示为矢量位A的旋度，并利用波导内场的纵向变化关系，将问题转化为求解二维标量亥姆霍兹方程。在推导弱形式时，将Leontovich边界条件（切向电场等于表面阻抗乘以切向磁场的叉乘）应用于波导的侧壁和盖板。对于侧壁损耗，通过边界条件关联了横向电场与磁场的切向分量，并最终体现为对横向波数k_tp的修正，引入了与表面阻抗相关的复数值修正项（k?_{p = kp - (1-j)δk²/(h kp)，其中δ为趋肤深度）。对于盖板（lid-wall）损耗，采用了类似的推导，但针对的是z方向的边界，也得到了对传播常数的修正形式。最终，通过离散化得到了一个线性系统，其解即为考虑损耗后的广义阻抗矩阵（GIM），进而通过公式S = (Z - I)(Z + I)^-1转换为广义散射矩阵（GSM）。}

TE^z模激励下的二维FEM/Leontovich分析

对于TE^z模（电场纵向分量为零）激励的情况，研究过程与TM^z模类似，但基于矢量位F进行展开。场表达式同样分离了横向和纵向依赖性。然而，在TE^z模情况下，将Leontovich边界条件应用于侧壁时，发现边界条件要求切向电场同时存在横向（t?）和纵向（?）分量，而基于理想导体假设的腔体场表达式无法自然满足这一条件。为了解决这个问题，研究人员从更一般的矢量三维亥姆霍兹方程出发，应用加权余量法（Galerkin方法）并利用格林恒等式，将面积分转化为体积分和边界积分。在边界积分中，分别处理端口边界和损耗性侧壁边界。通过精心选择矢量测试函数，并利用三角函数的正交性进行z方向的解析积分，最终成功地将三维矢量公式简化为二维标量公式，同时融入了侧壁损耗的影响。此公式包含了额外的矩阵项（如矩阵J，其元素涉及测试函数和基函数切向导数的乘积的线积分），这些项直接来源于对损耗性侧壁上边界条件的处理。

TM^z与TE^z激励的叠加

为了完整表征波导结在任意激励下的特性，需要将TM^z和TE^z两种极化模式的解结合起来。这意味着需要分别求解对应于两种极化的线性系统，然后将结果组合，以构成完整的广义散射矩阵。该矩阵按不同的纵向模态指数p进行分块对角化，反映了场沿波导高度方向（z方向）的不同变化模式。

研究结论表明，本文所提出的混合二维FEM/Leontovich方法能够有效地分析具有非理想导体壁的平面波导结。该方法通过引入表面阻抗边界条件，在二维有限元框架内精确地模拟了导体损耗，避免了繁琐的三维全波仿真，在计算效率和精度之间取得了良好平衡。讨论部分强调了该方法的普适性，它适用于任意形状的波导结，并且可以处理复数介电常数和磁导率的填充介质，从而扩展了其应用范围。这项研究的重要意义在于为毫米波和太赫兹频段高性能平面波导器件的精确设计和优化提供了强有力的理论工具和数值实现方案，对推进高频集成电路技术的发展具有重要价值。