波导模式分析
理论与物理
TE/TM模式分类
波导的模式有很多种吗？哪种会被使用呢？
矩形波导中TE₁₀是基本模式。只有在截止频率以上才能传播。在雷达和卫星通信中使用的Ku波段（12~18 GHz）标准波导WR-62中，TE₁₀的截止频率是9.49 GHz。
TE和TM是什么的缩写呢？
TE是Transverse Electric（横电波），指传播方向的电场分量 $E_z = 0$ 的模式。TM是Transverse Magnetic（横磁波），指 $H_z = 0$ 的模式。像同轴电缆那样存在TEM模式（$E_z = H_z = 0$）的结构也有，但单导体的波导中不存在TEM模式。这是因为波导内部的电场必须在横向形成“驻波”。
形成驻波，具体是指什么呢？
粗略地说，可以想象电磁波在波导壁之间反复反射，呈锯齿状前进。只有满足“恰好容纳整数个半波长”条件的模式才能存在。TEmn中的 m 和 n 分别表示宽度方向和高度方向容纳了多少个半波长。TE₁₀就是“宽度方向1个半波长，高度方向0个”。
模式的下标 $m, n$ 决定了横截面方向的驻波图案。矩形波导（宽度 $a$，高度 $b$，$a > b$）中的模式分类如下所示。
| 模式 | $E_z$ | $H_z$ | 最低阶 | 应用示例 |
|---|---|---|---|---|
| TEmn | 0 | $\neq 0$ | TE₁₀ | 标准传输（最常用） |
| TMmn | $\neq 0$ | 0 | TM₁₁ | 高次模式滤波器 |
| TEM | 0 | 0 | — | 仅同轴线（波导中不可） |
截止频率
截止频率是指，低于它的话电波就完全无法通过吗？
是的。低于截止频率 $f_c$ 时，电磁波会变成倏逝波，呈指数衰减。波导起到高通滤波器的作用。
矩形波导中TEmn模式和TMmn模式的截止频率由下式给出。
其中 $c$ 是介质中的光速，$a$ 是宽度（长边），$b$ 是高度（短边）。TE₁₀模式（$m=1, n=0$）时形式最简单。
那么，截止频率只由波导的宽度决定吗？例如实际数值会怎样呢？
例如在X波段（8~12 GHz）使用的WR-90波导，$a = 22.86$ mm，$b = 10.16$ mm。TE₁₀的截止频率是 $f_{c,10} = 3 \times 10^8 / (2 \times 0.02286) = 6.56$ GHz。下一个模式TE₂₀的截止频率是 $13.12$ GHz。所以在6.56~13.12 GHz的频带内，只有TE₁₀能传播，形成“单模带”。设计时通常使用截止频率的1.25倍到1.9倍范围。
以下是代表性波导（$a = 2b$ 的情况）最初几个模式按截止频率顺序的列表。
| 排名 | 模式 | $f_c / f_{c,10}$ | 备注 |
|---|---|---|---|
| 1 | TE₁₀ | 1.000 | 基本模式（设计频带） |
| 2 | TE₂₀ | 2.000 | 单模带的上限 |
| 3 | TE₀₁ | 2.000 | $a=2b$时与TE₂₀简并 |
| 4 | TM₁₁ / TE₁₁ | 2.236 | 高次模式的开始 |
传播常数与相速度
在截止频率以上传播时，波的速度和自由空间一样吗？
不一样。波导内的相速度比自由空间的光速快。但是，信息和能量传播的群速度比光速慢。这与相对论并不矛盾。
在截止频率以上（$f > f_c$）电磁波传播时，传播常数 $\beta$ 由下式给出。
由此可求得相速度 $v_p$ 和群速度 $v_g$。
两者之积始终满足 $v_p \cdot v_g = c^2$ 的关系。当 $f \to f_c$ 时，$v_g \to 0$（能量无法传播），$v_p \to \infty$。
在接近截止频率的频率下使用，群速度几乎为零，信号完全无法前进吗？
正是如此。所以在实际工作中，使用频率要高于截止频率的1.25倍。低于这个值则色散过大，脉冲信号会失真。在雷达脉冲压缩等应用中，如果不能准确把握群速度的频率依赖性，距离精度会急剧下降。
波阻抗
波导内的波阻抗因模式种类（TE/TM）而异。
其中 $\eta = \sqrt{\mu/\varepsilon} \approx 377 \, \Omega$（自由空间的固有阻抗）。TE模式在截止频率附近阻抗急剧增加，TM模式则趋近于零。
波阻抗随频率变化的话，波导和同轴电缆的连接处阻抗匹配似乎会很困难呢。
这正是微波电路设计的关键所在。在波导-同轴转换器（例如探针耦合器）中，通过调整探针的插入长度和位置，在宽频带内将VSWR（电压驻波比）抑制在1.2以下是实际工作的目标。在CAE中，通过S参数分析对反射特性S₁₁进行频率扫描，如果低于-20 dB则基本合格。
圆形波导的模式
圆形波导（半径 $a$）的模式用贝塞尔函数描述。TE模式的截止频率由 $J_n'(x_{nm}') = 0$ 的根决定，TM模式由 $J_n(x_{nm}) = 0$ 的根决定。
| 模式 | $x_{nm}'$ 或 $x_{nm}$ | 应用示例 |
|---|---|---|
| TE₁₁ | $x_{11}' = 1.841$ | 基本模式。最常用 |
| TM₀₁ | $x_{01} = 2.405$ | 轴对称。适合旋转对称结构 |
| TE₀₁ | $x_{01}' = 3.832$ | 低损耗模式。有利于长距离传输 |
圆形波导的TE₀₁模式低损耗是什么意思？
TE₀₁模式的壁面电流只沿圆周方向流动。频率越高，壁面电流密度越低，损耗越小。1970年代AT&T曾研究用于长距离通信的TE₀₁圆形波导，但后来被光纤取代。不过现在，它作为太赫兹频段的传输线再次受到关注。
波导的“截止”——与隧道中声音共鸣原理相同
波导的截止现象类似于隧道内的声音共鸣。如果隧道横截面宽度小于声波的半波长，就无法形成横向驻波，声音无法传播。波导也是如此，宽度 $a$ 小于半波长 $\lambda/2$ 的频率的电磁波无法传播。也就是说 $a > \lambda/2$ 即 $f > c/(2a)$ 是TE₁₀模式的传播条件。这种“高通滤波器”特性也有自然滤除设计外低频噪声的优点，有助于在雷达接收机前端排除不必要的低频干扰信号。
从麦克斯韦方程组推导波导模式
- 出发点：无源的麦克斯韦方程组可写为 $\nabla \times \mathbf{E} = -j\omega\mu\mathbf{H}$、$\nabla \times \mathbf{H} = j\omega\varepsilon\mathbf{E}$。假设时间因子为 $e^{j\omega t}$。
- $z$方向传播的假设：代入 $\mathbf{E}(x,y,z) = \mathbf{E}_t(x,y) e^{-j\beta z}$，则横向分量 $E_x, E_y, H_x, H_y$ 都可以用 $E_z$ 和 $H_z$ 表示。
- TE模式（$E_z = 0$）：$H_z$ 满足亥姆霍兹方程 $\nabla_t^2 H_z + k_c^2 H_z = 0$（其中 $k_c^2 = k^2 - \beta^2$），并在壁面施加诺伊曼条件 $\partial H_z / \partial n = 0$。
- TM模式（$H_z = 0$）：$E_z$ 满足亥姆霍兹方程 $\nabla_t^2 E_z + k_c^2 E_z = 0$，并在壁面施加狄利克雷条件 $E_z = 0$（理想导体）。
截止频率以下的倏逝波
- 当 $f < f_c$ 时，$\beta$ 为纯虚数 $\beta = -j\alpha$，电磁波以 $e^{-\alpha z}$ 指数衰减。这就是倏逝波。
- 衰减常数：$\alpha = k_c \sqrt{1 - (f/f_c)^2}$。在截止频率的一半（$f = f_c/2$）时，$\alpha \approx 0.866 k_c$，在一个波长的距离上衰减约50 dB。
- 倏逝波没有功率传播（坡印廷矢量的时间平均为零）。但由于近场耦合会有能量泄漏，因此波导法兰盘的间隙是导致性能劣化的原因。
波导参数一览
| 参数 | 符号 | SI单位 | 备注 |
|---|---|---|---|
| 截止频率 | $f_c$ | Hz | 由模式、截面尺寸决定 |
| 传播常数 | $\beta$ | rad/m | $f > f_c$ 时为实数 |
| 衰减常数 | $\alpha$ | Np/m | $f < f_c$ 时（倏逝波） |
| 相速度 | $v_p$ | m/s | 总是 $> c$ |
| 群速度 | $v_g$ | m/s | 总是 $< c$ |
| 波阻抗 | $Z_{TE}, Z_{TM}$ | $\Omega$ | 因模式类型而异 |
| 管壁损耗 | $\alpha_c$ | Np/m | 由壁面有限电导率引起 |
数值解法与实现
亥姆霍兹方程的推导
波导的模式分析，在计算机上是解什么方程呢？
从麦克斯韦方程组出发，归结为波导横截面上的二维亥姆霍兹方程。假设沿 $z$ 方向以 $e^{-j\beta z}$ 传播，就变成了仅包含横截面内分量的特征值问题。
在波导横截面上需要求解的亥姆霍兹方程形式如下。
其中 $\psi$ 在TE模式中是 $H_z$，在TM模式中是 $E_z$。$k_c$ 是截止波数，$k_c^2 = k^2 - \beta^2 = (\omega/c)^2 - \beta^2$。边界条件如下。
- TE模式：$\partial \psi / \partial n = 0$（壁面诺伊曼条件）
- TM模式：$\psi = 0$（壁面狄利克雷条件）
矩形波导存在解析解，但对于圆角或带脊等非标准截面，数值解法是必需的。
FEM公式化与边单元
听说波导的FEM分析中有“边单元”，和普通的FEM单元不同吗？
这是电磁场分析中最重要的点。如果用普通的基于节点的单元（标量单元）求解电场矢量，会产生大量物理上不存在的虚假模式——伪模式。这是因为无法在单元间保证麦克斯韦方程组中 $\nabla \cdot \mathbf{E} = 0$ 的条件。
那很麻烦啊。边单元就不会产生伪模式吗？
边单元（Nedelec单元）以“沿边的切向分量”作为自由度。这保证了电场切向分量的连续性，同时也能自然地表现法向分量的不连续性（在异种介质界面是物理上正确的行为）。结果，$\nabla \cdot \mathbf{E} = 0$ 被近似满足，伪模式被排除。现代的商用软件（HFSS、CST FEM、COMSOL RF Module）都基于边单元。
基于边单元的弱形式（伽辽金法）如下所示。
离散化后归结为广义特征值问题。
其中 $[S]$ 是旋度-旋度矩阵，$[T]$ 是质量矩阵，$\{e\}$ 是边上电场切向分量的未知向量。
作为特征值问题的公式化
特征值问题，意思是特征值对应截止频率吗？
正是如此。从特征值 $k_c^2$ 可以得到截止频率 $f_c = c k_c / (2\pi)$，特征向量则代表对应模式的电场分布。HFSS的Eigenmode Solver就是直接求解这类问题的。
实际应用中需要注意，特征值求解器的选择对计算效率有很大影响。
| 特征值求解器 | 特点 | 适用场景 |
|---|---|---|
| Lanczos法 | 最适合对称问题。高效计算少量特征值 | 波导模式分析的标准 |
| Arnoldi法 | 也可处理非对称问题。适用于有损介质的复特征值 | 有损波导的Q值计算 |
| FEAST法 | 并行计算指定频带内的全部特征值 | 密集模式的全面搜索 |
| 迭代子空间法 | 面向超大规模问题。内存效率高 | 超过100万自由度（DOF）的3D问题 |
FDTD方法
除了FEM，还有其他分析波导模式的方法吗？
FDTD（时域有限差分法）也可以使用。方法是在波导内激励脉冲，对时间响应进行FFT以确定谐振频率。CST的Transient Solver采用这种方法。与FEM相比，其优点是可以一次计算获得宽频带特性，但缺点是在高Q结构中需要很长的计算时间步长。
FDTD中，空间和时间在Yee网格上离散化，需要满足CFL稳定性条件 $c \Delta t \leq (\Delta x^{-2} + \Delta y^{-2} + \Delta z^{-2})^{-1/2}$。波导壁面作为理想电导体（PEC）边界 $\mathbf{n} \times \mathbf{E} = 0$ 实现。
频域 vs 时域的选用
频域分析（FEM特征值分析）类似于“用调谐器对准特定电台”——可以精确地逐个求解模式。时域分析（FDTD）类似于“同时录制所有频道”——可以一次获得宽频带响应。
なった
詳しく
報告