[论文解读] Computational Electromagnetism with Variational Integrators and Discrete Differential Forms

Ari Stern, Yiying Tong+2 more|arXiv (Cornell University)|Jul 30, 2007

Numerical methods for differential equations被引用 12

一句话总结

本文提出了一类基于时空离散微分形式的变分多辛数值方法，用于求解麦克斯韦方程组。证明了杨氏FDTD格式及其相关方法具有多辛性质，并可源自离散拉格朗日变分原理，从而在几何框架下统一了这些方法，实现了新型保结构数值方法的构建。

ABSTRACT

In this paper, we introduce a general family of variational, multisymplectic numerical methods for solving Maxwell’s equations, using discrete differential forms in spacetime. In doing so, we demonstrate several new results, which apply both to some well-established numerical methods and to new methods introduced here. First, we show that Yee’s finite-difference time-domain (FDTD) scheme, along with a number of related methods, are multisymplectic and derive from a discrete Lagrangian variational principle. Second, we generalize the

研究动机与目标

• 开发计算电磁学中尊重麦克斯韦方程组几何与辛结构的保结构数值方法。
• 利用离散微分形式，将杨氏FDTD等现有方法统一于同一变分与多辛框架之下。
• 推广这些方法，以构建求解时空麦克斯韦方程的新几何精确格式。

提出的方法

• 该方法采用离散微分形式对时空进行离散化，确保麦克斯韦方程组的几何结构在离散层面得以保持。
• 在时空网格上构建离散拉格朗日变分原理，确保所得格式具有多辛性质。
• 该方法在结构化网格上使用类似有限差分的模板，但其构造基于变分原理，而非启发式离散化。
• 离散形式定义于具有对偶网格的时空网格上，从而实现微分形式及其外导数的一致离散化。
• 该方法确保在空间与时间上均保持辛结构，从而实现长期稳定与高精度。
• 通过证明杨氏格式可自然地从离散变分原理导出，实现了对杨氏格式的推广，从而解释了其稳定性和结构保持性。

实验结果

研究问题

• RQ1 杨氏FDTD方法能否从离散变分原理导出，并证明其具有多辛性质？
• RQ2 经典电磁求解器（如FDTD）背后的几何结构是什么？
• RQ3 如何利用离散微分形式构造麦克斯韦方程的新多辛格式？
• RQ4 变分结构对电磁学中长期数值稳定性有何影响？
• RQ5 如何在保持辛与几何性质的前提下，对麦克斯韦方程的时空形式进行离散化？

主要发现

• 证明了杨氏FDTD格式具有多辛性质，并可源自离散拉格朗日变分原理，从而解释了其长期稳定性。
• 离散拉格朗日形式化提供了一个统一框架，用于理解与推广现有的电磁求解器。
• 使用离散微分形式可确保微分形式及其外导数的离散化具有一致性与结构保持性。
• 该方法自然导出新的多辛格式，其继承了连续方程的几何特性。
• 该方法表明，在时空离散化中实现几何结构保持是可行的，从而在长时间模拟中提升了数值精度。