FDTD与DGTD算法对比分析

1 引言

在现代电磁仿真领域，时域有限差分法（FDTD）与时域间断伽辽金法（DGTD）是两种主流的时域数值算法。FDTD凭借其简单高效的特点，长期占据工业仿真主流地位；而DGTD作为高阶方法，近年来在复杂几何建模和高精度计算中展现出独特优势。随着5G通信、自动驾驶、国防隐身技术等领域的快速发展，对电磁仿真的精度和适应性提出了更高要求。本文旨在从数学原理、计算特性、市场应用等维度对比分析两者的优劣势，展望其未来发展趋势。

2 算法原理

2.1 FDTD算法

FDTD基于Yee网格对麦克斯韦方程组进行时空离散，通过交替更新电场和磁场分量模拟电磁波传播。其核心特点包括：

• 结构化网格：规则六面体网格简化了计算流程，但难以精确拟合复杂曲面（易产生阶梯误差）。
• 显式时间推进：计算效率高，支持GPU加速，适合大规模问题。
• 局限性：低阶精度导致数值色散误差较大，这种误差会在宽带信号或长时仿真中积累，可能导致谐振频率偏移或远场辐射图失真。

更多关于FDTD算法的介绍，请见：什么是FDTD算法？

2.2 DGTD算法

DGTD融合有限元法（FEM）和间断伽辽金（DG）方法，采用非结构化网格与高阶基函数离散麦克斯韦方程，核心特点包括：

• 非结构化网格：可精确拟合任意几何形状（如曲面、纳米结构）。
• 高阶精度：通过多项式基函数（p≥3）显著降低数值色散，相位精度提升30~50%。
• 并行计算：单元级独立运算支持高效MPI/GPU加速，显式时间推进规避全局矩阵求解瓶颈。
  

3 算法对比分析

3.1 共同点

FDTD与DGTD都是时域电磁仿真方法的核心算法，均通过离散化麦克斯韦方程组直接模拟电磁场的动态传播过程。两者均具备宽频带分析能力，适用于瞬态电磁现象的仿真，例如脉冲信号传播、散射场计算及多物理场耦合问题。此外，二者均支持并行计算架构（如MPI和GPU加速），能够有效应对大规模计算任务，在5G通信、雷达系统等高频应用中展现技术通用性。然而，其底层数学框架和实现逻辑的差异，导致两者在精度、效率和应用场景上存在显著差别。

3.2 不同点

3.2.1 网格适应性

FDTD的核心限制源于其依赖的结构化Yee网格。这种六面体网格在规则几何（如矩形波导、平面天线）中计算效率极高，但其刚性结构在面对曲面、复杂材料界面或微纳尺度器件时，会因“阶梯近似”引入显著误差。例如，在模拟球形或曲边结构时（如图1：双直波导环形谐振器的网格构建），阶梯状网格无法精确描述几何形状，导致场分布畸变。

DGTD则通过非结构化网格彻底突破了几何限制。其网格单元可为四面体、六面体或混合形态，能够精确贴合任意复杂表面（如图3：飞机表面网格）。这种灵活性在纳米光子学中尤为重要——例如，在模拟亚波长光栅或等离子体共振结构时，DGTD的非结构化网格可避免FDTD的几何失真，确保近场增强效应的准确捕捉。此外，DGTD支持h-p自适应网格策略，即在强梯度场区域（如金属边缘或介质突变界面）动态加密网格（h-refinement）或提升基函数阶数（p-refinement），从而在保证精度的同时优化计算资源分配。

3.2.2 计算精度与数值误差

FDTD作为一种低阶方法，其精度高度依赖网格细化程度。在仿真高频电磁波时，即便网格尺寸满足Nyquist采样定理，仍会因中心差分近似引入数值色散误差——即不同频率分量传播速度的偏差。这种误差在宽带信号或长时仿真中累积，可能导致谐振频率偏移或远场方向图畸变。例如，在5G毫米波天线设计中，FDTD需通过超细网格（网格尺寸$\le \lambda /20$）抑制色散，显著增加计算成本。

DGTD则通过高阶多项式基函数（通常p≥3）实现“超收敛”特性。其场量在单元内以高阶多项式展开，能够以更粗的网格达到与FDTD细网格相当的精度。例如，Lumerical的案例表明，在米氏散射仿真中，DGTD的非结构化网格配合p=3基函数，可将相位误差降低30%-50%，且收敛速度更快（图3）。此外，DGTD的数值通量机制通过修正单元边界场跳跃，有效抑制了材料界面处的反射和耗散误差，使其在多介质、突变场问题中表现卓越。

下图展示了《Simulation of optical resonators using DGTD and FDTD》中DGTD与FDTD对于环形谐振器仿真中结果对比。DGTD仿真结果相对于参考结果吻合的较好，而FDTD仿真对于网格要求较高，运行成本高。

3.2.3 计算效率

FDTD的显式时间推进机制使其单步计算量极低，且易于实现GPU加速。例如，在均匀介质区域的仿真中，FDTD的电场-磁场交替更新仅涉及邻近网格点的简单代数运算，计算速度可达每秒数百万时间步。这种特性使其适合快速参数扫描和优化设计，如天线阻抗匹配或滤波器频响分析。

DGTD的单步计算复杂度显著高于FDTD。其每个单元需计算高阶基函数的内积与通量积分，且时间推进依赖多级Runge-Kutta方法，导致单步耗时比FDTD增加了2-5倍。但DGTD具有天然并行性——每个单元的场更新仅依赖相邻单元的通量数据，无需全局矩阵求解，因此在大规模分布式计算中（如千核级MPI集群），DGTD的并行效率可达90%以上，弥补了这一劣势。

3.2.4 技术生态与工业落地挑战

FDTD的成熟工业生态是其最大优势。 FDTD求解器已进入主流电磁仿真软件二十余年，应用生态较为成熟。此外，FDTD的算法逻辑简单，易于与机器学习结合——例如，利用神经网络预测最优网格参数，进一步压缩设计周期。

DGTD的工业应用仍面临生态短板。尽管Lumerical和COMSOL已推出DGTD模块，但其功能尚未完善，例如Lumerical DGTD缺乏GPU加速，光源类型单一；COMSOL中没有专门的DGTD算法接口，相关模块分布在声学模块、压电传感模块、以及波动光学模块。COMSOL官网于DGTD相关的案例基本都是声学以及压电模块的多物理场仿真，缺乏电磁学仿真案例。

3.3 总结

通过上述多维度对比可见，FDTD与DGTD并非简单的“替代”关系，而是基于应用场景的“互补”选择。在可预见的未来，两者将共同构成时域电磁仿真的技术双翼，推动从毫米波通信到量子器件的全频谱创新。

4 市场前景与未来趋势

4.1 市场概况

根据专业的市场研究机构（如MarketsandMarkets、Grand View Research、Allied Market Research等）发布的报告，2022年全球电磁仿真市场规模约50亿美元，市场将以年均7-10%的速度增长，到2030年市场规模能够达到100亿美元。

FDTD算法作为最早发展起来的时域电磁仿真方法，在光子学、微波工程、天线设计等领域占据主导地位，2022年全球电磁仿真市场中FDTD相关应用约占20~30%份额，其成熟的计算框架和高效GPU加速能力，使其在工业界保持稳定增长。

DGTD算法起步较晚，但凭借其高阶精度和非结构化网格优势，在航空航天、国防、纳米光子学等高端领域逐渐崭露头角。目前，DGTD市场份额较低（约1%），但根据市场研究机构预测，到2030年，DGTD市场份额有望提升至5~7%。

4.2 未来发展趋势

随着5G/6G通信、毫米波雷达、先进光子器件以及复杂电磁环境仿真需求的增长，对高精度、高效率电磁仿真算法的需求也在不断提高。未来，FDTD和DGTD可能会在不同应用场景下各展所长。

• FDTD：研究方向包括高阶FDTD方法、GPU优化、人工智能结合优化计算等，仍将在光子学、天线设计等领域保持重要地位。
• DGTD：受计算硬件进步驱动，在航空航天、国防、纳米光子学等高端应用领域将获得更大市场。

5 结论

综合来看，FDTD作为成熟的电磁仿真方法，在工业界仍具有较强的市场优势，预计在未来仍将保持稳定增长。而DGTD凭借其高精度和对复杂几何的适应性，虽然当前市场份额较低，但在高端应用领域具有较大的发展潜力。随着计算资源的提升和算法优化的推进，DGTD有望在未来占据更重要的位置，并与FDTD形成互补关系，共同推动电磁仿真技术的发展。

账号简介

我们（SimWorks）是一家国产自主仿真软件企业，致力于开发基于FDTD等算法的光电仿真软件“SimWorks有限差分”。
软件设计对标行业标杆，功能完备、上手容易，在仿真精度与速度上与国际领先水平相当。核心算法及代码均由核心团队自主开发，全面支持GPU加速与云计算。
软件目前公开提供免费云端试用，欢迎大家到官网查看以下说明进行试用。
试用说明：云计算试用说明

参考文献

[1] J N ,W P ,K B .Simulation of optical resonators using DGTD and FDTD[J].Journal of Optics A: Pure and Applied Optics,2009,11(11):114015-114015.
[2] Angulo L D, Alvarez J, Pantoja M F, et al. Discontinuous Galerkin time domain methods in computational electrodynamics: State of the art[C]//Forum Electromagn. Res. Methods Appl. Technol. 2015, 10(4): 1-24.
[3] MarketsandMarkets：Electromagnetic Simulation Software Market - Global Forecast to 2027
[4] Grand View Research：Computational Electromagnetics Market Size, Share & Trends Analysis Report
[5] Allied Market Research：Electromagnetic Simulation Software Market by Deployment Mode and Application