[论文解读] Note on transformation to general curvilinear coordinates for Maxwell's curl equations (Is the magnetic field vector axial?)

Dzmitry M. Shyroki|ArXiv.org|Jul 3, 2003

Electromagnetic Simulation and Numerical Methods参考文献 7被引用 4

一句话总结

本文主张，在广义曲线坐标系中，当磁感应强度 **H** 和磁通密度 **B** 被视为极矢量（真实矢量）而非轴矢量（赝矢量）时，麦克斯韦方程组的表述更加一致。其核心贡献在于：计算电磁学中广泛使用的完美匹配层（PML）方法仅在这一‘纯极矢量’表述下才能正确工作，因为标准的轴矢量约定会导致PML中复坐标缩放不兼容，从而破坏波导和谐振腔模拟中的数值稳定性和精度。

ABSTRACT

Arguments for H and B to be considered `axial', or pseudo vectors, are revisited. As a point against, we examine the complex-coordinate method for numerical grid truncation and mode loss analysis proved very successful in computational electrodynamics. This method is not compatible with convention that H and B are axial.

研究动机与目标

• 解决在非正交坐标变换下 **H** 和 **B** 变换行为长期存在的模糊性问题。
• 研究 **H** 和 **B** 的标准轴矢量约定与计算电磁学中使用的复坐标方法之间的兼容性。
• 证明麦克斯韦方程组的‘纯极矢量’表述对于一致实现完美匹配层（PML）是必要的。
• 表明PML方法——模拟开放边界电磁问题的关键工具——在传统的轴矢量处理 **H** 和 **B** 的方式下会失效。
• 为麦克斯韦方程组在曲线坐标系中的变换提供一个一致的框架，以在数值模拟中保持物理精度。

提出的方法

• 使用张量形式的核指标记法推导广义曲线坐标系中的麦克斯韦方程组，通过列维-奇维塔符号区分真实矢量与赝矢量。
• 提出一种‘纯极矢量’表述，其中 **H** 和 **B** 在坐标反演下按普通矢量变换，与标准的轴矢量约定形成对比。
• 对PML区域应用复坐标缩放，推导PML区域所需的介电常数和磁导率张量。
• 比较‘纯极矢量’与‘极-轴’约定下的PML公式，表明仅前者能导出标准PML张量（18s）。
• 通过在笛卡尔和极坐标网格上使用有限差分频域方法模拟阶跃折射率光纤的基模，对方法进行数值验证。
• 利用源自纯极矢量表述的边界条件，正确分类对称波导中的偶模与奇模，结果符合物理预期。

实验结果

研究问题

• RQ1 在广义曲线坐标变换下，**H** 和 **B** 的传统轴矢量处理是否仍保持一致性？
• RQ2 若将 **H** 和 **B** 视为轴矢量，能否从麦克斯韦方程组中一致推导出广泛使用的复坐标PML方法？
• RQ3 **H**/**B** 矢量类型约定对波导与谐振腔中数值模拟精度有何影响？
• RQ4 是否存在一种曲线坐标系中麦克斯韦方程组的表述形式，既能保持与PML方法的兼容性，又能避免赝矢量带来的复杂性？
• RQ5 在对称系统中，极-轴与纯极矢量表述在模场分类与场行为上存在哪些差异？

主要发现

• 标准PML公式（18s）仅与麦克斯韦方程组的‘纯极矢量’表述兼容，而与传统的‘极-轴’表述不兼容。
• 在轴矢量约定下，PML需要一个与标准PML根本不同的修正张量（18v），导致数值不稳定和模场错误。
• 数值模拟表明，基于纯极矢量表述的PML在非笛卡尔网格上产生的模有效折射率误差显著低于基于极-轴约定的PML。
• 纯极矢量表述能正确再现非正交坐标系中理想电导体的物理边界条件，如在 x¹ = 0 处 E₁ = 0 且 H₂ = H₃ = 0。
• 在对称波导与谐振腔中，纯极矢量表述与传统方法在模分类上结果一致，证明尽管数学处理方式不同，其物理预测具有等价性。
• 当 **H** 和 **B** 被视为轴矢量时，PML中使用的复坐标缩放方法无法恢复正确的张量形式，从而在该约定下从根本上破坏了该方法的理论基础。