FDTD法（有限差分时间域法） — CAE术语解说
FDTD法是什么
老师，在电磁波仿真中经常听到FDTD法，用一句话怎么解释呢？
简单说，就是将Maxwell方程在空间和时间上用有限差分直接离散化，在Yee单元上交替更新电场和磁场的显式解法。1966年由Kane Yee提出，正式名称是Finite-Difference Time-Domain法。因为不需要求解矩阵，只需按时间步长逐步迭代，所以算法极其简洁。
不用求矩阵，这和FEM的方法完全不同呢。具体使用Maxwell方程的哪种形式呢？
直接使用时间域的旋度方程。在没有源的情况下是这样的：
上面是法拉第定律，下面是安培定律（包含导电率 $\sigma$ 的损耗项）。这两个式子交替进行差分更新，这就是FDTD法的精髓。比如说隐形战机被雷达电磁波照射后的散射过程，可以像放电影一样逐帧追踪。
Yee单元的结构
听说"Yee单元"这种特殊的网格结构很关键，和普通网格有什么不同呢？
最核心的特点就是电场和磁场在空间上错开半个网格放置。考虑一个 $\Delta x \times \Delta y \times \Delta z$ 的长方体单元：
- 电场分量（$E_x, E_y, E_z$）放在单元边的中点
- 磁场分量（$H_x, H_y, H_z$）放在单元面的中心
这样的话，某个 $H$ 分量周围的4个 $E$ 分量恰好相隔半个网格，Maxwell方程中的旋度（curl）就能用中心差分自然、精确地表示出来。
不仅空间错开，时间上也错开吧？
没错！这叫蛙跳（Leapfrog）方式。时刻 $n$ 更新 $\mathbf{E}$，时刻 $n+\tfrac{1}{2}$ 更新 $\mathbf{H}$，再时刻 $n+1$ 更新 $\mathbf{E}$……循环往复。比如 $E_x$ 的一维更新式大概是这样的：
其中 $C_a, C_b$ 是由材料的誘電率 $\varepsilon$ 和导电率 $\sigma$ 决定的更新系数。每个单元只需要乘法和加法，所以和GPU并行计算的兼容性特别好。
单元是长方体的话，曲面就没法很好地表示了吧？
你观察得很敏锐。曲面会变成"阶梯近似"，这是FDTD法的一个弱点。比如圆形贴片天线的边缘用阶梯表示会导致共振频率偏移。解决办法有局部细化网格的次网格技术，或者Dey-Mittra法这样的曲面修正技巧。但如果整个结构都是曲面（比如球形雷达罩），FEM方法会更方便。
Courant稳定条件（CFL条件）
既然是显式法，时间步长肯定有限制吧？在结构分析中也有CFL条件…
完全一样的概念。FDTD法要求满足Courant稳定条件（CFL条件）。三维情况下是这样的：
这里 $c$ 是光速（真空中约 $3 \times 10^8$ m/s）。如果网格各向同性 $\Delta x = \Delta y = \Delta z = \Delta$，那么：
$$\Delta t \leq \frac{\Delta}{c\sqrt{3}}$$只要违反这个条件，解就会指数级发散，完全没法用。
实际工作中网格尺寸怎么定呢？相对于波长要多细？
经典的经验法则是最短波长的1/10以下。也就是最高频率 $f_{\max}$ 对应的网格尺寸：
比如28 GHz的毫米波5G天线，$\lambda_{\min} \approx 10.7$ mm，所以 $\Delta \leq 1.07$ mm左右。想要更高精度可以用λ/20。但问题是网格减半的话，三维情况下单元数变成8倍，加上Courant条件导致时间步数也翻倍，计算成本就会增加16倍。这个平衡的把握就是电磁解析工程师的真本领。
Courant数设到极限值也能稳定吗？
理论上是稳定的，但实务中一般设成0.9～0.95。这是为了防止浮点舍入误差导致偶尔超限。而且加上PML边界或分散性介质时，实际的稳定限界会略微降低，留点余量比较保险。
PML吸收边界条件
FDTD是有限计算域，域边界处的电磁波不会反射吧？
这正是个大问题。1994年Berenger发明的PML（Perfectly Matched Layer）就是革命性的解决方案。在计算域外周设置8～16个单元厚的"虚拟吸收层"，使进入其中的电磁波衰减并消失。
这和普通的有损介质（比如导体）放在那儿不一样吗？
差别很大！普通有损介质放在那儿会因为阻抗不匹配产生反射。而PML的妙处在于不管入射角、不管频率，理论反射系数都是零。物理上这种介质不存在，但在Maxwell方程的框架内完全匹配，所以可以被创造出来。
现在的主流是UPML（单向PML）和CPML（卷积PML）这种改进版。特别是CPML，对倏逝波的吸收性能更好，而且实现也更简单，大多数FDTD代码都采用它作为标准配置。
PML层数加得越多越好吗？
理论上是的，但实际上8～16层就能把PML表面反射控制在-60～-80 dB以下了。再多就是浪费内存和计算时间。更重要的是PML的导电率梯度分布的调校——一般用3～4阶多项式从PML内层向外逐步增加导电率，这是最佳实践。如果增长太急，离散化误差反而会增加反射。
宽带脉冲分析
听说FDTD法的强项是"一次计算就能得到宽带特性"，怎么做到的呢？
用高斯脉冲或调制高斯脉冲作为激励源。比如这样的波形：
$f_0$ 是中心频率，$\tau$ 是脉冲宽度。脉冲越短，频率带宽越宽。这个脉冲从输入端口打进去，计算时间演化，记下入射波、反射波和透射波的时间波形。计算完成后做FFT变换到频域，就一下子得到全频率范围的S参数（$S_{11}$, $S_{21}$等）了。
FEM的频域解法是逐频率求解吧？相比之下效率差太多了。
完全同意。比如要评估一个UWB天线的3～12 GHz频段，在100个频点上做评估。用频域FEM的话要求解100次线性方程组，但FDTD法一次计算就全搞定。5G毫米波或车载雷达（77 GHz频段）这种需要宽带评估的应用，FDTD法的计算效率优势压倒性。
反过来说，高Q值共振器或滤波器这种非常窄带的结构，时间域波形需要很久才能衰减，迭代步数会多到可怕，这种情况FEM会更高效。
天线的远场方向图也能用FDTD算吗？总不能把网格拉到很远吧…
好问题。用近场-远场变换（NTFF: Near-to-Far-Field Transformation）就解决了。在天线周围的虚拟闭合曲面上用FDTD计算近场，然后转换成等效电磁流，再解析地算出远场。所以不用把网格延伸到很远。这样就能高效地得到辐射方向图、指向性增益、偏振特性等等。
商用求解器对比：HFSS vs CST
实务中主要用CST Studio Suite和Ansys HFSS吧。FDTD法和它们的关系怎样？
大概是这样对应的：
| 项目 | CST Time Domain Solver | Ansys HFSS |
|---|---|---|
| 方法 | FIT（≒FDTD） | FEM（频域） |
| 网格 | 六面体（结构网格） | 四面体（非结构网格） |
| 擅长 | 宽带、过渡、大规模 | 高Q共振、曲面形状 |
| 带宽效率 | 一次全频域 | 逐频点求解 |
| GPU支持 | 非常强 | 有支持（有限） |
| 应用例 | 5G天线、EMC、雷达 | 导波管滤波器、连接器 |
CST的时间域求解器基于FIT（有限积分技术）体系，在六面体网格上数学上基本等价于FDTD法。所以CST能充分享受FDTD的宽带优势。
那HFSS是不是就没必要了？
绝对不是！导波管滤波器或空腔共振器这种Q值达到几千到几万的结构，FDTD法需要计算几十万个时间步才能让波形衰减，效率低得可怕。HFSS只需在共振频率周围求几个点。另外HFSS的自适应网格细分功能特别强，初学者甚至都不用太费劲调网格参数就能得到高精度结果，这在实务上很有价值。
最近HFSS也加入了Transient Solver，CST也有FEM频域求解器，两个工具在相互渗透。判断的关键是："频带宽度是否很大？""形状里有没有大量曲面？""要不要看过渡现象？"沿着这三个轴线来选择比较好。
实务技巧·常见陷阱
用FDTD法做仿真时有什么容易出错的地方吗？
常见的陷阱有这么几个：
- PML和结构间距太近 — 解析对象与PML边界之间至少要留λ/4以上的空间，否则倏逝波进入PML会增加反射。特别是贴片天线的接地面边缘要特别注意。
- 网格纵横比太大 — 如果 $\Delta x : \Delta z = 1 : 10$ 这样的扁平单元，Courant条件会让 $\Delta t$ 极小，导致计算时间爆增。一般纵横比控制在3以下比较好。
- 仿真停止时间过早 — 如果时间波形还没衰减到充分的水平就停止计算，FFT结果会有涟漪。要等到信号比峰值低30 dB以上才能停。
- 分散性介质的实现错误 — 频率相关的誘電率（Debye模型、Drude模型等）要是推导有问题，会破坏因果律，得到非物理的结果。用辅助微分方程（ADE）法来引入比较安全。
GPU计算的加速效果到底怎样？
FDTD法的单元更新完全独立，天生适合GPU的超并行架构。CST的GPU加速能达到CPU唯一计算的10～50倍。比如整个车体的EMC分析（10亿单元规模）在4块GPU上可以几小时内完成。开源的FDTD代码（gprMax、MEEP、openEMS等）也陆续支持CUDA/OpenCL，没有商业工具也能试用。
总结
把前面的内容归纳一下，FDTD法的关键点是什么？
浓缩成5点的话：
- Yee单元：电场和磁场空间错开半格子，用中心差分自然准确地离散化Maxwell的旋度
- 蛙跳时间迭代：显式法，无需求矩阵。最适合GPU并行
- Courant稳定条件：$\Delta t \leq \Delta/(c\sqrt{3})$。违反就立即发散。实务设成0.9～0.95倍
- PML吸收边界：开放空间用有限域模拟。CPML是现在的标配
- 宽带脉冲分析：高斯脉冲+FFT，一次计算得全频特性
FDTD法的另一大魅力是算法直观、易于调试，还能把时间演化做成动画看，物理现象"活生生"地呈现。电磁波从天线辐射出去、和散射体相互作用、产生衍射……这些过程可以像看电影一样观察。其他方法很难有这样的优势。
准确理解CAE术语是团队沟通的基石。 — NovaSolver项目致力于支持实务工作者的学习。
告诉我们您在FDTD实务中遇到的挑战
NovaSolver项目面对CAE工程师的日常难题——网格尺寸决策、PML参数调校、计算成本削减——致力于找到解决方案。您的实务经验将成为改进工具的原动力。
联系我们（筹备中）相关主题
参考
更详细
错误