[论文解读] Curl-free positive definite form of time-harmonic Maxwells equations well-suitable for iterative numerical solving

V.Е. Moiseenko, O. Ågren|arXiv (Cornell University)|Mar 28, 2020

Electromagnetic Simulation and Numerical Methods参考文献 10被引用 2

一句话总结

本文提出了一种无旋度、正定的时间谐 Maxwell 方程形式，采用矢量势（A）和标量势（Φ）、磁场（H）以及电位移（D）。通过用拉普拉斯算子和散度梯度项替代旋度算子，并消除主导微分算子中的空间变化系数，该方法确保了离散系统的对角占优性和正定性，从而实现了共轭梯度法的高效收敛。该方法在三维问题中实现了最优迭代次数与 N^{1/3} 成比例的性能。

ABSTRACT

A new form of time-harmonic Maxwells equations is developed and proposed for numerical modeling. It is written for the magnetic field strength, electric displacement, vector potential and the scalar potential. There are several attractive features of this form. The first one is that the differential operator acting on these quantities is positive. The second is absence of curl operators among the leading order differential operators. The Laplacian stands for the leading order operator in the equations for the magnetic field strength, vector potential and the scalar potential, while the gradient of divergence stands for the electric displacement. The third feature is absence of space varied coefficients in the leading order differential operators that provides diagonal domination of the resulting matrix of the discretized equations. A simple example is given to demonstrate the applicability of this new form of time-harmonic Maxwells equations.

研究动机与目标

• 开发一种时间谐 Maxwell 方程的形式，以避免标准形式中旋度算子引入的虚假模态。
• 消除所得线性系统的符号不定性，从而可使用快速迭代求解器（如共轭梯度法）。
• 通过从主导微分算子中移除空间变化系数，确保离散系统的对角占优性。
• 通过为所有场量使用同一网格，简化有限差分离散化过程，并便于边界条件的实现。
• 通过三维数值基准测试，展示该方法在迭代求解器中更优的收敛行为。

提出的方法

• 采用归一化的矢量势 A = a / (ik₀) 和标量势 Φ = -iφ / (k₀)，以平衡场量大小，实现一致的离散化。
• 从标准时间谐 Maxwell 方程出发，利用势函数推导方程，并应用库仑规范条件 ∇·A = 0。
• 将系统表示为二次型 Q = ∫[|LA|² + |LH|² + |LD|² + |LΦ|²] dV，通过最小化该型推导弱形式。
• 采用标准三维有限差分法在均匀结构化网格上对所得弱形式进行离散化，使用体积加权的节点贡献。
• 通过在边界上设置 Φ = 0，并在边界和棱边上使用一致的五点模板，施加边界条件。
• 使用共轭梯度法求解线性系统 Mx = b，不使用预条件，通过归一化残差和误差监控收敛性。

实验结果

研究问题

• RQ1 能否基于矢量势和标量势构造出无旋度、正定的时间谐 Maxwell 方程形式？
• RQ2 主导项中若无旋度算子和空间变化系数，是否能确保离散系统具有对角占优性？
• RQ3 该形式是否能实现共轭梯度法的高效收敛，避免符号不定系统中常见的停滞现象？
• RQ4 在该形式下，三维问题中迭代次数如何随网格规模变化？
• RQ5 能否为所有场分量实现一致、稳定且精确的有限差分离散化，且仅使用单一网格？

主要发现

• 所提出的格式生成了正定的离散系统，使得共轭梯度法能够保证收敛。
• 主导微分算子为拉普拉斯算子（用于 H、A、Φ）和散度的梯度（用于 D），从而从主导项中消除了旋度算子。
• 主导算子中无空间变化系数，确保了所得线性系统矩阵的对角占优性。
• 数值结果表明，共轭梯度求解器收敛平稳且迅速，残差随迭代次数持续下降。
• 达到一半离散化误差时的最优迭代次数与总网格节点数 N 的关系近似为 N^{1/3}，表明具有高效的迭代性能。
• 三维电流源的解析解（含正弦分量）被数值解成功匹配，验证了方法的准确性和一致性。