1 基本概念
1.1 方法定义
时域有限差分法是一类通过时间步进方式求解波动问题的数值计算方法。其基本做法是将连续的时间与空间变量离散为网格点，再利用差分近似替代偏微分方程中的导数项，从而逐步计算场量在各个时刻的变化。
与直接求解稳态响应的方法不同，FDTD更强调“随时间演化”的过程。它能够在一次仿真中获得宽频信息，因此特别适合研究脉冲传播、瞬态响应以及复杂结构中的波动行为。
1.2 发展背景
时域有限差分法的形成与电磁数值计算的发展密切相关。早期研究者希望建立一种既能处理复杂几何，又能描述瞬态电磁现象的通用算法，FDTD由此逐渐成熟。
该方法之所以得到广泛应用，主要是因为它具有公式结构清晰、物理意义直观、易于程序实现等特点。随着计算机性能提升，FDTD从最初的基础电磁问题，扩展到光学、微波、声学及其他波动领域。
1.3 适用问题类型
FDTD适用于具有明确波动特征的偏微分方程问题，尤其适合分析瞬态过程、边界反射和散射行为。只要问题可被写成一组时间演化方程，通常都可以考虑采用这一方法。
1.3.1 波动方程的数值求解
对于一般波动方程，FDTD可通过对时间和空间的差分离散，追踪波包传播速度、波前形状以及反射与干涉现象。它常用于声波、弹性波和标量场传播等问题。
1.3.2 电磁场与麦克斯韦方程组
FDTD最经典的应用对象是麦克斯韦方程组。通过将电场与磁场交替更新，可直接模拟电磁波在自由空间、介质和复杂结构中的传播过程，并可用于分析反射、透射、谐振和辐射特性。
1.3.3 其他物理场中的扩展应用
除电磁学外，FDTD思想还可推广到热传导、流体波动、地震波传播以及某些量子与生物物理模型中。其核心仍是将连续场量离散化，并以显式迭代方式推进求解。
2 数学基础
2.1 偏微分方程离散化
FDTD的基础在于把偏微分方程转化为代数更新关系。离散化后，场量只在有限个网格点和离散时刻上被计算，从而将连续问题变为可迭代的数值问题。
2.1.1 时间离散
时间离散通常采用等间隔步长，将连续时间轴划分为一系列时刻。场量在相邻时间层之间通过更新公式递推，这种处理方式使算法具有明确的时间推进顺序。
2.1.2 空间离散
空间离散则是将计算区域划分为规则或非规则网格。每个网格点代表局部场值，空间导数通过相邻点之间的差分近似，使波动传播能够在网格上被模拟出来。
2.2 差分格式
差分格式决定了导数近似的具体方式，也影响算法的精度与稳定性。FDTD通常采用显式、中心型差分结构，以便实现简单且便于物理解释。
2.2.1 显式差分
显式差分的特点是下一时刻的场值可以直接由当前时刻和已知量计算得到，无需联立求解大规模方程组。这种形式计算效率较高，但稳定性通常需要满足特定条件。
2.2.2 中心差分
中心差分利用前后相邻网格点的信息近似导数，通常具有较好的对称性和较高精度。它在FDTD中被广泛使用，尤其适合构建电场与磁场的交错更新公式。
2.3 数值稳定性
稳定性是判断差分算法能否长期可靠运行的重要标准。若时间步长或网格设置不当，误差可能在迭代中放大，导致结果失真甚至发散。
2.3.1 CFL条件
CFL条件是FDTD中最常见的稳定性限制之一，它将时间步长与空间网格尺寸联系起来。一般而言，时间推进不能过快，否则波在一个时间步内传播的距离会超过网格允许范围，进而引发数值不稳定。
2.3.2 收敛性分析
收敛性分析用于判断网格细化后数值解是否逐渐逼近真实解。若算法满足一致性与稳定性，通常可获得较好的收敛表现；反之，则可能出现较大误差或不可靠结果。
3 经典FDTD算法
3.1 Yee网格
Yee网格是经典FDTD的核心结构，由电场与磁场在空间上交错存储而构成。该网格设计能够自然对应麦克斯韦方程组的旋度关系，并提高数值更新的稳定性。
3.1.1 电场与磁场交错布置
在Yee网格中，电场分量与磁场分量并不位于同一位置，而是错开半个网格单元。这种布置使得各分量之间的空间差分关系更加紧密，也便于实现时间上的半步交替更新。
3.1.2 网格单元结构
一个典型网格单元包含若干相互错位的场分量节点。该结构虽然看起来较复杂，但能够更准确地表达局部旋度和散度关系，是FDTD获得成功的重要原因之一。
3.2 更新方程
更新方程是FDTD的实际计算核心。它们把麦克斯韦方程离散后转化为一系列递推公式，使电磁场在每个时间步上交替演化。
3.2.1 电场更新
电场更新通常基于当前磁场的空间差分结果进行计算。通过这种方式，电场在下一时刻可由磁场旋度驱动得到，从而体现电磁耦合关系。
3.2.2 磁场更新
磁场更新则由当前电场的差分值决定。电场和磁场相互促进、交替推进，构成一个持续演化的数值循环。
3.3 时间推进策略
FDTD的时间推进具有典型的显式交替特征，电场与磁场分别在错开的时间层上更新。这一机制保证了算法简洁，同时维持了较好的物理对应关系。
3.3.1 半步交替更新
半步交替更新指电场与磁场在时间上相差半个步长。磁场先用电场更新，再由更新后的磁场推进电场，如此循环往复，实现整个场分布的逐步演化。
3.3.2 初始条件设定
仿真开始时必须给定初始场值，必要时还要设置源项、材料参数与边界条件。初始条件的合理性会影响早期瞬态响应，并决定数值计算是否符合物理情景。
4 边界与源处理
4.1 吸收边界条件
由于计算区域有限，若不加处理，波在边界处会发生非物理反射。吸收边界条件的目的就是尽量让波离开计算域时不产生明显回波。
4.1.1 一阶吸收边界
一阶吸收边界是较早且实现简单的一类方法，能够对近似正入射波产生一定吸收效果。它结构简洁，但对斜入射和复杂波形的抑制能力有限。
4.1.2 完美匹配层
完美匹配层是更常用也更有效的吸收技术。它在计算域外设置一层人工介质，使入射波逐渐衰减并尽量不被反射回内部区域，因此被广泛用于高精度仿真。
4.2 激励源建模
激励源用于在模型中引入电磁波或其他波动输入。不同源形式对应不同研究目标，源的设置直接影响频谱范围和场分布特征。
4.2.1 脉冲源
脉冲源通常具有短时、宽带的特点，适合一次性激发多个频率分量。它常用于获取系统的频率响应、共振特性以及瞬态传播过程。
4.2.2 连续波源
连续波源以单一频率或窄频带持续激励系统，便于分析稳态场分布、驻波结构及谐振增强效应。在天线和谐振腔研究中较为常见。
4.2.3 平面波入射
平面波入射常用于研究材料层、周期结构和散射体的电磁响应。通过设定合适的注入方式，可以模拟均匀波前照射目标结构的情形。
4.3 周期边界与对称边界
当结构本身具有周期性或对称性时，可利用特殊边界条件减少计算量，同时保持模型的物理代表性。
4.3.1 周期结构建模
周期边界适用于光子晶体、超材料和重复单元结构。只需计算一个基本单元及其边界关系，就能反映整个无限周期结构的行为。
4.3.2 镜像对称处理
镜像对称边界可用于具有对称几何的系统。借助电场或磁场的对称性，可以缩小仿真区域，并减少冗余计算。
5 材料建模
5.1 均匀介质
均匀介质是最基础的材料模型，参数在空间中保持不变。此时更新公式最为简单，常用于算法验证和基础物理分析。
5.2 色散材料
色散材料的电磁响应随频率变化，不能仅用常数参数描述。FDTD中通常借助附加微分方程或卷积形式来刻画其频散特性。
5.2.1 Drude模型
Drude模型常用于描述自由载流子主导的材料响应，例如某些金属在特定频段的行为。它能够反映导电电子对电磁场的延迟响应。
5.2.2 Lorentz模型
Lorentz模型适合描述带有共振特征的介质。通过引入谐振振子概念，可以模拟材料在特定频率附近出现的强吸收或色散现象。
5.3 各向异性材料
各向异性材料的物理参数随方向改变，场量之间的耦合关系也更复杂。处理这类材料时，更新方程通常需要使用张量形式的参数描述。
5.4 非线性材料
非线性材料中，响应不再与激励成简单线性比例，而可能依赖场强大小。此类材料可用于研究倍频、光学开关和强场效应等现象。
5.5 损耗介质
损耗介质会在波传播过程中吸收能量，使场强随传播距离逐渐衰减。其数值模型通常需要在更新方程中加入电导率或等效损耗项。
6 数值实现
6.1 网格剖分
网格剖分决定了仿真区域的空间离散方式，也影响计算精度与资源消耗。合理的剖分方案应兼顾几何表达能力与计算效率。
6.1.1 均匀网格
均匀网格在各方向采用固定步长，结构简单、实现方便。它适合规则几何或需要较强算法稳定性的场景。
6.1.2 非均匀网格
非均匀网格允许局部区域采用更细的分辨率，而在远离关注区的位置使用较粗网格。这种方法有助于降低总体计算量，但实现上更为复杂。
6.2 算法流程
FDTD的程序流程通常包括初始化、时间循环和结果保存等步骤。其整体结构较为规整，便于并行化与模块化实现。
6.2.1 初始化
初始化阶段需要设置网格、材料参数、边界条件、源项以及初始场分布。若参数配置准确，后续迭代会更加稳定可靠。
6.2.2 迭代更新
迭代更新是整个算法的主体。程序在每个时间步依次更新磁场、电场，再处理源、边界和材料响应，直到达到预设时间或收敛条件。
6.2.3 结果输出
结果输出包括场分布记录、监测点数据、频谱信息及派生物理量。良好的输出设计有助于后续分析反射、透射和谐振特性。
6.3 误差与精度控制
FDTD的精度受网格、时间步长、边界处理和材料模型等多种因素影响。为提高可信度，通常需要进行误差估计和参数收敛测试。
6.3.1 数值色散
数值色散是指离散算法中不同频率分量传播速度出现偏差的现象。它会导致波形展宽、相位失真，是FDTD中常见的误差来源之一。
6.3.2 网格分辨率要求
网格必须足够细致，才能准确刻画波长较短的场变化。一般而言，目标波长越短，所需网格越密，否则容易出现几何失真或计算误差增大。
6.3.3 时间步长选择
时间步长既要满足稳定性要求，也要保证足够的时间分辨率。过大的步长会破坏稳定性，过小则会显著增加计算耗时。
7 后处理与结果分析
7.1 场分布可视化
场分布可视化是理解仿真结果的重要手段。通过二维切片、三维云图或动态动画，可以直观看到波的传播路径、局域增强区域和反射位置。
7.2 频谱分析
频谱分析用于从时域数据中提取频率信息，判断系统对不同频段的响应情况。这对于识别谐振、带隙和滤波特性尤为重要。
7.2.1 傅里叶变换提取
傅里叶变换可将时域信号转换为频域谱线，从而获得各频率成分的幅度与相位信息。它是FDTD后处理中最常见的分析工具之一。
7.2.2 共振峰识别
共振峰识别主要用于寻找系统响应显著增强的频点。峰值位置可反映结构的固有频率，而峰宽则常与损耗和品质因子有关。
7.3 参数反演
参数反演是根据仿真或测量结果推断材料或结构参数的过程。FDTD可提供前向模型支持，使反演分析更具可操作性。
7.3.1 反射与透射计算
通过监测入射波、反射波和透射波的幅值变化，可以计算反射率与透射率。这类结果常用于评估薄膜、滤波器和界面结构的性能。
7.3.2 散射截面分析
散射截面描述目标对入射波的散射能力，是研究粒子、微结构和天线辐射特性的重要指标。FDTD能够较方便地从远场或近场数据中提取相关量。
8 典型应用
8.1 光学器件模拟
FDTD广泛用于模拟波导、谐振腔、透镜、滤光片和纳米光子结构。它能够较好地处理亚波长特征和复杂边界，因此在集成光学中十分常见。
8.2 天线与微波器件
在天线设计中，FDTD可用于分析辐射方向图、输入阻抗和带宽特性。对于滤波器、腔体和传输线等微波器件，它也能提供较全面的时域和频域信息。
8.3 超材料与光子晶体
超材料和光子晶体通常具有周期性或人工微结构特征，FDTD适合研究其带隙、负折射、局域共振及异常传播现象。由于可直接处理复杂单元结构，它在该领域应用非常普遍。
8.4 生物电磁问题
FDTD可用于分析电磁波在人体组织中的传播、吸收和局部能量沉积情况。相关研究常见于医疗成像、治疗评估和生物组织参数分析等方向。
8.5 传感与成像系统
在传感和成像中，FDTD可帮助设计高灵敏结构并评估信号响应。它能够模拟目标与波场的相互作用，为系统优化提供依据。
9 优势与局限
9.1 方法优势
FDTD之所以长期受到重视，源于其明确的物理含义和较强的通用性。对于许多瞬态波动问题，它都能提供直接、可解释的数值结果。
9.1.1 直观性
该方法的时间推进过程与真实波动演化相对应，便于理解场量如何随时间变化。对于初学者和工程应用人员而言，这种直观性尤其有价值。
9.1.2 宽频分析能力
只需一次宽带激励，就可获得多个频率点的信息，因此FDTD在谱分析方面很高效。相比逐频求解的方法，它更适合研究宽带响应。
9.1.3 复杂结构适应性
FDTD能够较容易地处理不规则边界、复杂材料分布和多层结构。只要网格足够细，几何和材料的建模灵活性通常较强。
9.2 主要局限
尽管应用广泛，FDTD仍存在一定缺点，尤其在高精度、大规模和长时程问题中更为明显。
9.2.1 计算资源消耗
当仿真区域较大或需要高分辨率时，存储和运算成本会迅速上升。三维问题尤其容易对内存和处理器性能提出较高要求。
9.2.2 网格伪影
由于离散化近似，曲面、斜边和细小结构可能被网格“台阶化”处理，从而引入伪影。若网格不够精细，这类误差会更明显。
9.2.3 长时间仿真成本
对于需要观察慢过程或高品质因子共振的情形，仿真时间往往较长。此时即便单步计算成本不高，总体耗时也可能相当可观。
10 相关方法与扩展
10.1 有限元法
有限元法通过将区域划分为单元，并在单元内构造近似函数来求解场问题。它对复杂几何的适应性较强，常用于与FDTD进行互补比较。
10.2 有限体积法
有限体积法强调守恒形式，常见于流体与传输问题。对于某些波动系统，它也可作为数值离散框架，但在电磁时域分析中不如FDTD常用。
10.3 频域方法
频域方法直接在频率空间求解稳态响应，适合单频或少数频点分析。与FDTD相比，它不强调时间演化，但在某些谐振和线性问题中效率较高。
10.4 高阶FDTD
高阶FDTD通过采用更高阶的差分近似来降低数值误差和色散效应。它在提升精度方面具有优势，但公式更复杂，对边界和稳定性处理也提出更高要求。
10.5 混合数值方法
混合数值方法将FDTD与其他求解技术结合，例如与有限元、解析展开或边界积分法耦合。此类方法常用于兼顾局部高精度与整体计算效率。