介质谐振器的电磁场仿真
理论与物理
什么是介质谐振器
介质谐振器就是一块陶瓷在谐振吗？和金属空腔有什么不同？
简单来说，它是一种电磁波被限制在高介电常数陶瓷（$\varepsilon_r = 20 \sim 90$）内部并发生谐振的器件。与用金属壁包围的空腔谐振器不同，它是利用介质与空气的介电常数差，通过“内部反射”来限制电磁波的。
诶，没有金属壁也能被限制住吗？
是的。这和光纤的全反射原理相同。$\varepsilon_r$ 越大，波长就越短——因为 $\lambda = \lambda_0 / \sqrt{\varepsilon_r}$，所以即使频率相同，谐振器也可以做得更小。例如，对于 $\varepsilon_r = 36$ 的陶瓷，尺寸可以做到金属空腔的 $1/6$。它被用于5G基站的滤波器和振荡器的频率稳定化。
我听说Q值也很高，大概是多少呢？
无载Q值（$Q_u$）极高，在1万到10万之间。金属空腔由于导体损耗，Q值通常停留在数千左右，但介质可以使用 $\tan\delta$ 在 $10^{-4}$ 以下的低损耗材料。所以基站滤波器才能干净利落地滤除带外信号。
| 项目 | 介质谐振器（DR） | 金属空腔谐振器 |
|---|---|---|
| 谐振原理 | 介电常数差引起的内部反射 | 金属壁的全反射 |
| 典型 Qu | 10,000 〜 100,000 | 2,000 〜 10,000 |
| 尺寸 | 小型（$\propto 1/\sqrt{\varepsilon_r}$） | 大型 |
| 频率稳定性 | 可实现 $\tau_f \approx 0$ ppm/K | 依赖于金属热膨胀 |
| 主要用途 | 基站滤波器、振荡器 | 粒子加速器、雷达 |
谐振模式（TE/TM/HE/EH）
DR的谐振模式，和空腔谐振器的TE/TM模式不同吗？
金属空腔中存在纯粹的TE模式和TM模式，但DR是开放结构，所以会变成轴向同时残留电场和磁场的“混合模式”。HE模式是以磁场成分为主的混合模式，EH模式是以电场成分为主的混合模式。
滤波器设计中最常用的是哪个模式？
应用最广泛的是 TE01δ 模式。这是轴向（z方向）没有电场分量 $E_z$ 的最低阶模式，电场沿圆周方向分布。下标 $\delta$ 表示轴向限制不完全——也就是说电磁场会略微渗透到介质外部。这个模式的Q值最高，与其他模式的频率分离度也好，是滤波器设计的标准选择。
在双模滤波器中，有时会利用 HE11δ 模式的双重简并（两个正交极化），用一个DR实现两个谐振。这样可以将元件数量减半，因此在卫星通信的输入输出滤波器中很受欢迎。
谐振频率的近似公式
谐振频率是怎么计算的？能解析求解吗？
对于圆柱形DR，虽然不存在精确解，但常用Courtney模型或Kajfez的近似公式。TE01δ 模式的一个代表性近似公式是：
$a$ 和 $h$ 分别代表什么？
$a$ 是圆柱的半径，$h$ 是高度。例如，对于 $\varepsilon_r = 36$、$a = 5\,\text{mm}$、$h = 4\,\text{mm}$ 的DR，$f_0 \approx 4.8\,\text{GHz}$。实际工作中，先用这个公式估算初始值，然后用FEM进行精确优化。
作为更精确的近似公式，还有Kishk等人提出的公式（在较宽的介电常数和形状比范围内，精度约为 $\pm 2\%$）：
Q值（品质因数）
我知道Q值高滤波器特性会变好，但Q值具体是衡量什么的呢？
Q值是“储存的能量”与“每周期损耗的能量”之比。用公式表示就是：
对于DR，损耗可分为三部分：介质损耗 $Q_d$、导体损耗 $Q_c$（如果有金属外壳）、辐射损耗 $Q_r$：
它们各自大概是多少值？
通常介质损耗占主导，可以用 $Q_d = 1/\tan\delta$ 近似。对于 $\tan\delta = 10^{-4}$ 的材料，$Q_d \approx 10{,}000$。如果放入金属屏蔽外壳，$Q_c$ 大约为数万，$Q_r$ 在外壳内实际上可以视为无穷大。所以通常 $Q_u \approx Q_d$。
| 损耗因素 | 公式 | 典型值 | 主导条件 |
|---|---|---|---|
| 介质损耗 $Q_d$ | $1 / \tan\delta$ | 5,000 〜 100,000 | 大多数情况下占主导 |
| 导体损耗 $Q_c$ | $\propto \sqrt{f} / R_s$ | 10,000 〜 50,000 | 外壳壁面较近时 |
| 辐射损耗 $Q_r$ | $\propto (\varepsilon_r)^{3/2}$ | $> 100{,}000$ | 仅开放结构 |
温度系数 $\tau_f$
基站放在室外，温度变化很大吧？谐振频率稳定吗？
用温度系数 $\tau_f$（单位: ppm/K）来评估温度变化引起的频率偏移：
这里 $\tau_\varepsilon$ 是介电常数的温度系数，$\alpha_L$ 是线膨胀系数。实用的DR材料通过调整成分，可以将 $\tau_f$ 控制在 $\tau_f \approx 0\,\text{ppm/K}$ 附近。例如Ba-Zn-Ta系（BZT）实现了 $\tau_f = 0 \pm 2\,\text{ppm/K}$，即使在 $-40°$C 到 $+60°$C 的环境变化下，频率偏移也能控制在 0.01% 以下。
精度真高！CAE也会一起做温度分析吗？
是的。使用COMSOL等多物理场工具进行热-电磁场耦合分析，可以评估由外壳内温度分布引起的局部 $\varepsilon_r$ 波动如何影响谐振频率。特别是对于大功率滤波器，自发热效应不可忽视，因此耦合分析变得必不可少。
陶瓷的烧制程度决定频率
介质谐振器的材料开发，正是与“陶瓷烧制程度”的斗争。通过以0.1mol%为单位调整BaTiO3-BaZrO3系的混合比，并以$\pm5°$C的精度控制烧成温度，可以改变 $\varepsilon_r$，而这直接关系到谐振频率。过去依赖资深技术人员的经验，但现在主流是结合CAE参数化研究和实验设计法（DoE）进行材料设计。这可以说是“匠人技艺”被数字孪生取代的一个好例子。
各项的物理意义
- 波长缩短效应 $\lambda = \lambda_0 / \sqrt{\varepsilon_r}$：在高介电常数介质中，光速降低为 $c/\sqrt{\varepsilon_r}$，因此波长变短。这是DR能够小型化的根本原理。$\varepsilon_r = 36$ 时，波长为自由空间的 $1/6$。
- 谐振条件：电磁波在DR内部形成驻波。对于圆柱形DR，电场分布在径向上用贝塞尔函数、轴向上用指数函数（渗透部分）描述。
- $\tan\delta$（损耗角正切）：介质损耗的度量。设 $\varepsilon_r = \varepsilon_r' - j\varepsilon_r''$，则 $\tan\delta = \varepsilon_r'' / \varepsilon_r'$。微波频段的DR材料要求 $10^{-4}$ 以下。
- $\tau_f$（温度系数）：温度每升高1K时谐振频率的变化（ppm）。由 $\tau_\varepsilon$（介电常数的温度变化）和 $\alpha_L$（热膨胀）的贡献决定。可以通过材料成分控制到 $\tau_f \approx 0$。
假设条件与适用范围
- Courtney近似假设为 孤立DR（周围无导体壁）。放入外壳后，与壁的相互作用会导致频率偏移
- $\varepsilon_r \gg 1$ 的假设：当 $\varepsilon_r < 10$ 时，电磁波渗透较大，近似精度变差
- $a/h$ 的范围：$0.5 < a/h < 2.5$ 左右为实用范围。极端扁平或细长形状需要其他近似
- 各向同性假设：部分DR材料在烧结时具有各向异性。单晶系需要定义各方向的 $\varepsilon_r$
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 谐振频率 $f_0$ | Hz | 实用上为GHz频段（1〜100 GHz）。$\omega_0 = 2\pi f_0$ |
| 相对介电常数 $\varepsilon_r$ | 无量纲 | DR材料：20〜90。具有温度、频率依赖性 |
| 损耗角正切 $\tan\delta$ | 无量纲 | 需要 $10^{-4}$ 以下。$Q_d = 1/\tan\delta$ |
| 温度系数 $\tau_f$ | ppm/K | 目标: $|\tau_f| < 5$ ppm/K |
| 半径 $a$, 高度 $h$ | m | 实用上为mm〜cm。微波频段为数mm |
数值解法与实现
特征值问题的公式化
既然有近似公式，为什么还需要FEM分析呢？
近似公式只适用于孤立的简单圆柱DR。实际产品中，DR可能放在金属外壳里，或者通过耦合孔与其他DR相连，或者带有支撑台和调谐螺丝。对于这种复杂形状，要准确计算谐振频率和Q值，就需要对麦克斯韦方程组进行全波求解。
怎么把麦克斯韦方程组变成特征值问题呢？
假设时谐场并分离出 $e^{j\omega t}$，就可以得到关于电场的亥姆霍兹型特征值方程：
这里 $k_0 = \omega/c$ 成为特征值。用边单元离散化后，就变成广义特征值问题：
$[S]$ 相当于旋度-旋度矩阵，$[T]$ 相当于质量矩阵。将 $k_0^2$ 作为特征值求解，就可以从每个特征值通过 $f_n = k_0 c / (2\pi)$ 得到谐振频率。Q值则是将特征值扩展到复数域，用 $Q = \text{Re}(k_0) / (2 \cdot \text{Im}(k_0))$ 来计算。
边单元（Nedelec单元）与节点单元
必须用边单元吗？用普通的节点单元不行吗？
这是电磁场FEM最重要的一点。用节点单元（标量单元）近似矢量场时，无法自动满足 $\nabla \cdot \mathbf{E} = 0$ 的散度条件。结果，大量物理上不存在的“伪模式”会混入特征值中，使得难以分辨真正的模式。
为什么边单元不会出现伪模式？
边单元（Nedelec单元）将自由度配置在边上，只保证电场切向分量的连续性。法向分量可以不连续——这正确地再现了在不同介电常数界面上电场法向分量不连续的物理现象。结果，散度条件自动得到满足，伪模式被排除。
| 单元类型 | 自由度配置 | 连续性保证 | 伪模式 | 对DR分析的适用性 |
|---|---|---|---|---|
| 节点单元（标量） | 节点 | 全分量连续 | 会产生 | 不适用 |
| 一阶边单元（CT/LN） | 边 | 仅切向分量 | 排除 | 基本OK |
| 二阶边单元（LT/QN） | 边+面 | 仅切向分量 | 排除 | 推荐（高精度） |
网格划分策略
DR的网格划分有什么特殊的技巧吗？
有三个要点。首先，介质与空气的界面要划分得足够细。因为这里的电场变化剧烈。标准是介质内的波长 $\lambda_d = \lambda_0 / \sqrt{\varepsilon_r}$ 的 $1/10$ 以下。其次，金属壁附近——需要表皮深度 $\delta_s$ 的 $1/3$ 以下的单元尺寸。最后，空气区域的外侧边界至少要距离介质 $\lambda_0 / 4$ 以上。
なった
詳しく
報告