空洞共振器的电磁场仿真

完全相同。电子微波炉的磁控管是2.45GHz的空洞共振器。由金属壁完全围成的空间，在特定频率处形成驻波——这就是"共振"，是空洞共振器的基本原理。

应用非常广泛。通信卫星的带通滤波器、粒子加速器的加速空洞、雷达发振器、5G基站的滤波器——都是空洞共振器的应用。FEM仿真的主要目标是求得共振频率和Q值。Q值越高，共振越尖锐；滤波器需要高Q值，而宽带天线则需要低Q值。应用场景不同，设计目标也完全不同。

先从最简单的矩形空洞开始。对边长为 $a \times b \times d$ 的矩形空洞，将麦克斯韦方程在边界条件（壁面处切向电场为零）下求解，可以得到优美的分析解。

圆筒空洞（半径 $a$、高度 $d$）稍微复杂一点，涉及贝塞尔函数的零点。

不用背。查表或数值计算就行。实务中真正重要的是理解：不同的模式指数 $(m,n,p)$ 对应不同的共振频率和电磁场模式。低阶模式（TE₁₀₁、TM₀₁₀）的场分布简单易用，但频率升高后模式密集，这就是"模式拥挤"问题，是仿真的一大难点。

简单说就是"振动的持久力"。敲玻璃杯会"嗡嗡"长鸣，Q值高。敲纸箱只有"嘭"一声，Q值低。空洞共振器也一样，壁面损耗越小，Q值越高。

从表面阻抗 $R_s$ 和壁面切向磁场 $\mathbf{H}_t$ 来算。

实务中要区分"空载Q值"和"负载Q值"。不考虑外部耦合的是空载Q（$Q_0$），考虑耦合的是负载Q（$Q_L$）。耦合系数 $\beta$ 与之间有关系。

正是。反射为零的最优耦合状态。粒子加速器追求临界耦合，而滤波器设计中则故意用过耦合（$\beta > 1$）或欠耦合（$\beta < 1$）。典型的Q值参考：铜空洞用于通信滤波器时 $Q_0 \sim 5{,}000$～$20{,}000$，超导空洞（铌材）可达 $Q_0 \sim 10^{10}$，差异巨大。

沿空洞轴向（例如 $z$ 方向），电场无 $z$ 分量的是TE（横向电场）模，磁场无 $z$ 分量的是TM（横向磁场）模。矩形空洞的最低阶通常是TE₁₀₁模，也是微波炉的基本模式。

因为TE₀₁₁模的壁面电流不穿过接合部（角和边界），不受接触阻抗影响。而且这个模式有个独特性质：频率越高Q值越高（通常相反）。正因为这些特性，它被用在频率标准器中——包括世界时间标准的铯原子钟，内部也装有空洞共振器。

实际空洞有耦合孔、同轴探针、调谐螺钉、表面微小变形等。这些"理想形状的偏差"导致分析解失效。这时就用FEM——把连续体支配方程离散化成固有值问题求解。

$[S]$ 是curl-curl矩阵（刚度矩阵），$[T]$ 是质量矩阵。固有值 $k_0^2$ 给出共振频率 $f = k_0 c / (2\pi)$，固有矢量 $\{e\}$ 给出电磁场模式。和结构的固有振动分析框架完全一样。

好问题。麦克斯韦方程有 $\nabla \cdot \mathbf{D} = 0$ 的约束（空洞内无电荷）。用普通节点单元对电场做标量型插值时，这个散度零条件自动不满足，导致计算结果中混入"虚假的共振模"。

边界单元（Nedelec单元）将自由度分配在"边"上，保证电场的切向分量连续。法向分量允许在单元界面处不连续，这样离散层面上 $\nabla \cdot \mathbf{D} = 0$ 被准确取入，伪模就消失了。现在高频FEM分析实际上用边界单元是标准配置。

2阶边界单元拟合曲面形状性能更好。同样网格密度下，相比1阶单元精度显著提升。特别是圆筒空洞这样有曲面的形状，1阶单元的形状近似误差会控制精度。HFSS的自适应网格细化也是建立在2阶边界单元之上。

有3种主要方法。

设计初期用摄动法（快），详细设计用S参数法（含耦合结构影响）。超导空洞这样Q值极高的情况，复特征值法数值精度不足，必须用摄动法。

目的不同。

Q值10000的空洞，能量衰减到1/e需要 $Q/\pi \approx 3{,}000$ 个周期。10 GHz的空洞，这意味着0.3微秒的计算时间。时间步长通常几皮秒，要10万步以上。而FEM的特征值分析一次矩阵分解就搞定，效率差异巨大。

标准分析流程是这样的。

表皮深度远小于壁厚（铜壁表皮深度微米级、壁厚毫米级），电磁波无法透过壁，所以PEC近似精度极高。Q值计算时只需后处理时加入壁面表面阻抗 $R_s$ 即可。但对于超导薄膜空洞这样壁厚接近表皮深度的情况，需要直接建模壁厚。

空洞共振器网格最关键的是"每波长的单元数"。经验法则：2阶边界单元至少6个单元/波长；1阶单元需10～15个。2阶单元反而更经济

HFSS的自适应网格细化很强，共振频率变化达到收敛判据（默认0.02%）就停止。但初始网格太粗的话容易错误收敛。使用者有责任确保初始波长分辨率达到要求。

空洞分析主要用这4种。

不存在但能表示电磁场对称性。比如矩形空洞的TE₁₀₁模关于某平面电场对称，那个平面上可以加PMC条件，把模型缩成1/2。但要注意——在对称面加PEC和加PMC会得到不同模式。PEC得到反对称模，PMC得到对称模。两个都算起来能得全部模式。

矩形空洞有3个对称平面，能缩成1/8模型。网格单元数变成1/8，内存大约1/8，计算时间因求解器特性不同，大约缩到1/10～1/30。圆筒轴对称模（TM₀₁₀等 $m=0$）可降维成2D轴对称，计算量只有3D的百分之一以下。

但对称性模型无法看到与它对称性不符的模式。如果要全面了解所有模式，还是得用完整模型。在设计某个特定模式时，对称性利用很高效。

HFSS最大优势是Eigenmode Solver的完熟度。基于2阶边界单元FEM+自适应网格细化，把共振频率自动收敛到高精度。用法流程这样：

Eigenmode结果表格自动显示。内部用摄动法从壁面表面阻抗评价。把壁面材料从PEC改为有限导率（比如铜）时，$R_s$ 自动算出Q值。注意：PEC壁面显示 $Q = \infty$。要得到真实Q值，必须设定壁面材料。

CST也有Eigenmode Solver，基于JDM（Jacobi-Davidson法）的固有值求解，在FIT网格上运作。六面体网格自动生成速度快，简单形状的空洞能高效处理。

CST的杀手锏是"宽带S参数分析"。时间域求解器打入短脉冲，一次运算通过FFT得全带宽S参数。多模耦合空洞的宽带特性用CST更快，比HFSS的频率扫描高效得多。

RF Module的Eigenfrequency Study能做空洞分析。COMSOL特色是多物理耦合。粒子加速器的超导空洞，壁面损耗产热 → 温度升高 → 表面阻抗变化（BCS阻抗的温度关联） → Q值变化，这个全过程COMSOL内可一体求解。HFSS和CST需要配合外部热分析工具。

自适应网格细化能力弱。HFSS那样基于"固有值收敛"的自动细化COMSOL没有，网格收敛性得用户手动确认。大规模模型的固有值求解扩展性也不如HFSS。

可以。Elmer FEM是边界单元型固有值分析，用ARPACK做特征值求解，空洞模式能算。OpenEMS是FDTD型，从时间域振铃衰减或FFT得共振频率。学术研究级别足够，但比商用软件GUI弱，故障时需自力更生。粒子加速器社区有ACE3P（SLAC开发）和Superfish（LANL开发）等无费工具。

单空洞的共振频率和Q值求解，HFSS最高效可信。结合空洞设计有带通滤波特性的话，CST的宽带分析优势显著。需要热耦合用COMSOL。资源有限用Elmer+手工调。不同需求用不同工具，无一通用之选。

典型伪模问题。检查项：

把电场的散度 $\nabla \cdot \mathbf{E}$ 画图看。伪模里 $\nabla \cdot \mathbf{E} \neq 0$。另外，网格细化时，伪模频率漂变大；真实模有网格收敛性。

Q值偏差是最常遇到的实务问题。主要原因和对策：

$\Delta$ 是RMS（均方根）表面粗糙度，$\delta$ 是表皮深度。$\Delta / \delta \ll 1$ 时无视不计，但10 GHz铜壁 $\delta \approx 0.66$ μm，机械加工面 $\Delta \approx 0.5$～1.5 μm，差不多大小，影响明显。电解抛光把 $\Delta$ 降到0.05 μm，Q值就接近理论值。

高阶模密集帯厄介。対策：

固有値求解器不良大2。

1：Adaptive Mesh収束

2：固有値求解器本体失敗

鉄則3。深呼吸——焦設定問題複雑化。次最小再現案例作——空洞単純化、分析解一致確認（引算最効率）。最後一度1変再実行——複数同時変更何効分。

用形状函数 $N_i$ 来近似未知量:

用数学表达这样：

用数学表达这样：

连续体的支配方程离散化后，得这样的代数方程组：

其中 $[K]$ 是整体刚度矩阵（或等效的系统矩阵），$\{u\}$ 是未知节点变量矢量，$\{F\}$ 是外力矢量。

收敛速度：二次要素 $O(h^2)$ 的阶数改进误差（光滑解的情况）

空洞共振器分析的实务分析流程和注意要点来讲。

1. 前处理 (Pre-processing)

2. 求解 (Solving)

3. 后处理 (Post-processing)

Ansoft Corporation開発3D高周波電磁界。2008年AnsysAnsoft買収。

Computer Simulation Technology () 開発。2016年Dassault Systèmes買収SIMULIA統合。

、調子！実際手動一番的勉強。分聞。

空洞共振器分析対応主要商用CAE工具的機能比较、各製品的歴史的背景詳述。

Ansoft Corporation開発3D高周波電磁界。2008年AnsysAnsoft買収。

Computer Simulation Technology () 開発。2016年Dassault Systèmes買収SIMULIA統合。

1986年設立。MATLAB連携的FEMLAB開始、後COMSOL改名。多物理场強。

空洞共振器分析的工具選定以下考慮:

、調子！実際手動一番的勉強。分聞。

空洞共振器分析最新的研究動向先進的方法見。

症状：求解器指定反復回数内収束異常終了

考原因：

対策：