2. 电磁仿真基础
电磁仿真的所有算法，最终都追溯到同一组物理方程——麦克斯韦方程组。理解这组方程如何被“离散化”，是理解 FDTD、FEM、MoM、FIT 各自能力边界的关键。本节从方程出发，梳理时域与频域两条技术路线，并给出主流算法的原理、数据结构与适用场景，为后续 HPC 并行策略的选择打下基础。
一、核心原理：从麦克斯韦方程组到数值离散
麦克斯韦方程组以微分形式描述电场与磁场的相互激发关系：
∇×E = -∂B/∂t （法拉第定律）
∇×H = J + ∂D/∂t （安培-麦克斯韦定律）
∇·D = ρ （高斯定律）
∇·B = 0 （磁通连续性）
其中 D = εE，B = μH，J = σE解析解只对极少数规则几何（如无限大平板、理想波导、球体散射）成立。工程问题中的天线、封装、腔体、超材料几乎都无解析解，必须数值求解。数值方法的核心动作只有三步：离散化（把连续空间/时间切成网格或面元）、构造代数方程（把微分/积分算子变成矩阵或迭代更新式）、求解并后处理（得到 S 参数、场分布、RCS 等）。不同算法的差异，本质上是这三步的取舍不同。
二、时域方法：FDTD 与 FIT
FDTD（时域有限差分）由 Yee 于 1966 年提出，核心是把电场和磁场在空间上交错半个网格、在时间上交错半个时间步，形成“蛙跳”更新。电场分量位于网格棱边中点，磁场分量位于面心，这种 Yee 网格天然满足旋度方程的几何关系。更新公式完全显式，无需解线性方程组：
H^{n+1/2} = H^{n-1/2} - (Δt/μ) · ∇×E^n
E^{n+1} = E^n + (Δt/ε) · (∇×H^{n+1/2} - σE^n)稳定性受 Courant 条件约束：Δt ≤ 1/(c·√(1/Δx²+1/Δy²+1/Δz²))。FDTD 一次时域推进即可通过傅里叶变换得到宽频带响应，非常适合超宽带天线、生物电磁剂量、光子晶体等场景。其内存与网格量成正比，对电大尺寸问题压力大，但天然适合 GPU 并行。
FIT（有限积分技术）与 FDTD 同属时域网格法，但离散对象是网格边上的电压和面上的磁通，而非场分量本身。它从积分形式的麦克斯韦方程出发，用拓扑矩阵（离散旋度、散度算子）构造更新式，因此能严格保持电荷守恒和磁通连续性。FIT 对曲线网格和复杂边界的适应性优于标准 FDTD，是 CST 等商业软件的核心算法之一。
三、频域方法：FEM 与 MoM
FEM（有限元法）基于变分原理或伽辽金弱形式，将求解域剖分为四面体/六面体单元，用基函数展开场量，最终形成稀疏线性方程组 Ax = b。FEM 对复杂几何、各向异性材料、多尺度结构的适应能力最强，广泛用于微波器件、封装互连、电机电磁场。代价是需要求解大规模稀疏矩阵，频域扫描时每个频点都要重新求解，计算量随频率点数线性增长。
MoM（矩量法）基于表面积分方程，只离散金属表面或介质分界面，用格林函数自动满足辐射边界条件，因此不需要吸收边界和空气盒。未知量少、精度高，特别适合开放域辐射和散射问题（天线、RCS）。但阻抗矩阵是稠密的，直接求解复杂度 O(N³)，即使采用 MLFMA 加速，内存和计算仍随 N 快速增长，通常限制在电小到电中尺寸。
四、主流算法适用场景对比
| 算法 | 离散对象 | 方程类型 | 矩阵特性 | 典型场景 | HPC 瓶颈 |
|---|---|---|---|---|---|
| FDTD | 体网格场分量 | 显式时间步进 | 无矩阵 | 超宽带天线、生物电磁、光子晶体 | 内存带宽、网格量 |
| FIT | 边/面通量 | 显式或隐式 | 稀疏 | 通用微波器件、EMC | 网格生成、时间步 |
| FEM | 四面体单元 | 频域弱形式 | 稀疏 | 复杂几何、多物理场、封装 | 稀疏直接/迭代求解 |
| MoM | 表面面元 | 频域积分方程 | 稠密 | 天线辐射、RCS、开放域 | 稠密矩阵存储与求解 |
五、案例与误区复盘
- 认为“网格越密越准”。FDTD 网格过密会因 Courant 条件使时间步极小，仿真时间爆炸；FEM 网格过密则矩阵规模失控。
- 用 MoM 求解封闭腔体内部问题。MoM 的格林函数假设开放边界，封闭结构应优先选 FEM 或 FIT。
- 忽略色散材料建模。FDTD 处理色散需引入 ADE 或 PLRC 等辅助微分方程，直接套用普通更新式会得到错误结果。
- 把“算法精度”与“并行效率”混为一谈。FDTD 并行扩展性最好，但精度受数值色散限制；FEM 精度高，但并行受线性求解器制约。
掌握这些算法的原理与边界，才能在面对具体电磁问题时做出合理选择，也才能进一步判断：哪些计算适合 GPU 加速，哪些受限于内存带宽，哪些必须依赖分布式稀疏求解器。这正是后续 HPC 资源优化的起点。