3. 电磁模型基础

在射频电磁模型与电路模型一致性验证中，全波数值方法是获取“电磁真值”的主要手段。它直接求解麦克斯韦方程组，不依赖等效电路假设，因此常被用作电路模型、行为模型或测量结果的参照基准。但全波方法并非单一算法，FEM、FDTD、MoM 在求解域、离散方式、未知量规模和适用结构上差异显著。选错方法或设置不当的网格与边界条件，会让“电磁真值”本身产生系统性偏差，后续一致性验证就失去意义。

三种主流方法的求解机制与适用边界

FEM（有限元法）基于微分形式的波动方程，将求解域剖分为四面体或六面体单元，在每个单元内用基函数展开场量，最终形成稀疏线性方程组。它的优势在于能自然处理复杂几何、非均匀介质和各向异性材料，适合波导不连续、腔体、滤波器、封装互连等封闭或半封闭结构。FEM 通常配合自适应网格细化，可在场变化剧烈处自动加密。其代价是开放区域需要人工截断，且频率扫描时若采用直接求解，内存和计算量随频率点增长明显。

FDTD（时域有限差分法）在 Yee 网格上交替采样电场和磁场，按时间步进显式更新。它天然适合宽带响应、瞬态问题和开放辐射结构，一次时域仿真经傅里叶变换即可得到宽频带 S 参数。FDTD 的网格必须足够细以解析最小几何特征和最短波长，时间步长受 Courant 条件限制，因此细长结构或高 Q 腔体需要很长的稳定时间。色散材料、薄层和曲面边界需要特殊处理，否则会引入阶梯近似误差。

MoM（矩量法）基于积分方程，只对金属表面或介质分界面离散，未知量集中在表面电流或等效面电流上。它不需要对自由空间剖分，因此对天线、散射体、线缆和开放金属结构非常高效，且自动满足辐射边界条件。但 MoM 产生的是稠密矩阵，直接求解复杂度约为 O(N³)，迭代求解也受条件数影响。对于复杂非均匀介质体，MoM 需要体等效原理或混合方法，否则效率会迅速下降。

网格质量如何决定结果可信度

网格是数值方法的“采样率”。在 FEM 中，单元尺寸通常取介质中波长的 1/10～1/6，但真正决定精度的是场梯度大的区域：金属边缘、介质交界面、馈电点、缝隙和薄层。若这些位置网格过粗，S 参数在谐振点附近会出现频率偏移和幅度压缩。实际操作中可先以 λ/8 全局网格试算，再在回波损耗变化剧烈处做局部加密，观察 S11 曲线是否随网格细化收敛。若两次加密后谐振频率变化小于 0.5%，可认为网格基本收敛。

FDTD 的网格约束更严格。空间步长 Δx 需同时满足最小几何特征和最短波长要求，通常取 λ_min/20 以上；时间步长必须满足 Courant 条件 Δt ≤ Δx/(c√3)（三维）。若结构中有 0.1 mm 的细缝而工作频率对应波长约 30 mm，全局用 0.1 mm 网格会导致网格量爆炸。此时可用亚网格或共形网格，但需注意稳定性。判断 FDTD 结果是否可信，可观察时域波形是否在仿真窗口内衰减到噪声以下，以及宽频 S 参数是否满足因果性和无源性。

MoM 的网格针对表面电流，通常要求每个波长至少 10 个基函数，曲面结构还需满足几何逼近精度。对于线天线，段长取 λ/10 左右；对于贴片天线，表面三角形边长取 λ/10～λ/15。MoM 的收敛性可通过增加基函数阶数或细化网格观察输入阻抗是否稳定。若输入阻抗随网格加密变化超过 2%，说明当前离散不足。

边界条件：把开放空间“关”进有限计算域

边界条件直接决定波是否被正确吸收或反射。FEM 和 FDTD 处理开放区域时，常用 PML（完美匹配层）或 ABC（吸收边界条件）。PML 通过坐标拉伸或各向异性材料吸收 outgoing 波，理论上反射为零，但实际性能依赖层数、电导率分布和入射角。对于宽角扫描或低频截止情况，PML 可能失效，此时可改用高阶 ABC 或边界元混合。

FDTD 中 PML 通常设置 8～16 层，电导率按多项式或几何级数渐变，理论反射系数可低于 -60 dB。若仿真中观察到端口反射异常或 S11 在低频段抬高，应检查 PML 是否离结构太近，一般要求 PML 内边界距辐射体至少 λ/4。FEM 的 PML 可做成球面或长方体，同样需要与辐射结构保持足够距离。

对于封闭金属腔体，边界条件为理想电导体（PEC），此时无需 PML，但需注意腔体 Q 值高时频率分辨率要足够细。对于周期结构如频率选择表面，可用周期边界条件（PBC）或 Floquet 端口，此时单元网格必须与周期方向对齐，否则相位关系会出错。

方法选择与验证流程

面对一个具体射频结构，可按以下逻辑选择方法：若结构以金属表面和开放辐射为主，优先 MoM；若包含复杂非均匀介质体且尺寸不超过几个波长，优先 FEM；若需要宽频带瞬态响应或结构电大且以规则网格为主，优先 FDTD。对于混合结构，可采用 FEM-MoM 或 FDTD-MoM 混合方法，但需注意交界面上场量连续性和网格匹配。

无论选用哪种方法，建议建立“网格—边界—收敛”三级检查：第一级，用粗网格和默认边界快速试算，确认端口设置、材料参数和频率范围无误；第二级，逐步加密网格并加厚 PML，观察关键指标（S11 谐振频率、S21 幅度、输入阻抗）是否收敛；第三级，与解析解、测量值或另一种数值方法交叉比对。只有通过收敛性检查的电磁结果，才能作为电路模型一致性验证的可靠基准。

网格与边界条件不是仿真前的“一次性设置”，而是需要与结果一起迭代评估的变量。在射频电磁模型与电路模型一致性验证中，只有先确认全波结果本身已收敛，后续的电路模型偏差分析才有意义。