[论文解读] Multisymplectic Geometry Method for Maxwell's Equations and Multisymplectic Scheme

Hongling Su, Mengzhao Qin|ArXiv.org|Feb 24, 2003

Numerical methods for differential equations参考文献 5被引用 3

一句话总结

本文通过将变系数介质中麦克斯韦方程组的拉格朗日量经勒让德变换导出多辛哈密顿形式，建立了一套用于非均匀、各向同性、无损耗介质中麦克斯韦方程组的多辛几何框架。提出了一种九点普雷萨曼格式，精确保持多辛守恒律，在电磁波传播的长时间数值模拟中展现出优异的稳定性和准确性。

ABSTRACT

In this paper we discussed the self-adjointness of the Maxwell's equations with variable coefficients $ε$ and $μ$. Three different Lagrangian are attained. By the Legendre transformation, a multisymplectic Bridge's (Hamilton) form is obtained. Based on the multisymplectic structure, the multisymplectic conservation law of the system is derived and a nine-point Preissman multisymplectic scheme which preserve the multisymplectic conservation law is given for the Maxwell's equations in an inhomogeneous, isotropic and lossless medium. At last a numerical example is illustrated.

研究动机与目标

• 建立具有可变 ε 和 μ 系数的麦克斯韦方程组的自伴变分形式。
• 通过拉格朗日密度的勒让德变换，推导出多辛哈密顿形式。
• 构建一种保持结构的数值格式，以守恒多辛守恒律。
• 通过具有精确解的波传播数值模拟验证该格式。

提出的方法

• 利用逆变分问题技术，为变系数介质中的麦克斯韦方程组推导出三种不同的拉格朗日密度。
• 应用勒让德变换，将拉格朗日系统转换为多辛哈密顿形式，揭示其内在几何结构。
• 从哈密顿系统的几何结构中推导出多辛守恒律。
• 通过直接对哈密顿形式进行离散化，构造一种保持多辛守恒律的九点普雷萨曼多辛格式。
• 在中间阶段引入辅助变量（U, V, P, Q），推导出最终的九点格式，随后以矩阵形式表达，便于数值实现。
• 在空间和时间方向均采用交错的跃迁型离散化，以保持对称性与守恒性质。

实验结果

研究问题

• RQ1 在何种条件下，变系数介质中的麦克斯韦方程组具有自伴性，从而可建立变分形式？
• RQ2 如何从麦克斯韦方程组的拉格朗日量中推导出多辛哈密顿结构？
• RQ3 在空间变化的 ε 和 μ 条件下，麦克斯韦方程组的多辛守恒律具有何种形式？
• RQ4 能否构造一种保持多辛结构的数值格式，同时准确模拟电磁波传播？
• RQ5 在长时间模拟中，九点多辛格式的性能与精确解相比如何？

主要发现

• 本文建立了麦克斯韦方程组在一阶形式下具有可变 ε 和 μ 时的自伴性必要且充分条件，从而支持变分形式的建立。
• 根据场表示的选择，推导出三种不同的拉格朗日量。
• 通过勒让德变换获得复数多辛哈密顿形式，揭示了系统内在的几何结构。
• 证明九点普雷萨曼格式能精确保持多辛守恒律，确保长时间模拟的稳定性。
• 当 ε = μ = 1 且 J = 0 时，数值结果与精确解高度一致，数值解在 1000 个时间步长内保持波状行为，误差增长极小。
• 在 t = 1000 时的绝对误差呈现出规则的振荡模式，与精确解的形状高度吻合，表明具有高精度与结构保真度。