A Hybrid DEC-SIE Framework for Potential-Based Electromagnetic Analysis of Heterogeneous Media

技术摘要：一种用于非均匀介质基于势能电磁分析的混合 DEC-SIE 框架

问题陈述
在复杂、多材料及无界环境中进行电磁场分析面临着显著的计算挑战。传统的微分方程（DE）方法（如有限元法 FEM）擅长模拟非均匀材料，但需要人工领域截断（例如通过完美匹配层 PML 或吸收边界条件 ABC）来处理开放区域，这增加了开放区域问题的复杂性。相反，积分方程（IE）方法能够自然地满足无界区域的辐射条件，但在处理不均匀介质时表现不佳。此外，标准的 DE 求解器通常会遭受“低频崩溃”问题，即随着频率趋于零，由于电能与磁能尺度之间的失衡，系统变得病态。虽然基于势能的公式（使用磁矢量势 $\mathbf{A}$ 和电标量势 $\Phi$）具有天生的免疫低频崩溃的特性，但将这些公式与用于开放区域的表面积分方程（SIE）相结合在历史上一直存在困难。针对可渗透介质的标准基于势能的 SIE 通常涉及复杂的算子结构，包含多达 14 个不同的积分算子，从而造成了沉重的数值和实现负担。

方法论
作者提出了一种混合数值框架，将用于内部非均匀区域的离散外微分（DEC）与用于外部均匀区域的简化 SIE 相耦合。

1. 基于势能的公式化： 该方法利用洛伦兹规范下的麦克斯韦方程组，以磁矢量势 $\mathbf{A}$ 和电标量势 $\Phi$ 进行表述。这种方法本质上避免了低频崩溃，并能自然地与多物理场耦合对齐。
2. 领域分解：
   • 内部（DEC）： 使用 DEC 对非均匀有界区域进行离散化。与传统的 FEM 不同，DEC 是一种坐标无关的框架，它保留了拓扑不变性，并允许在非结构化网格上进行保结构离散。
   • 外部（SIE）： 使用 SIE 处理均匀无界外部区域，这能自然满足辐射条件并消除对人工截断的需求。
3. 核心创新：标量分量重构：
   • SIE 简化： 作者通过将问题表示为矢量势的笛卡尔分量（$A_x, A_y, A_z$）及其法向导数，对 SIE 进行了重构。这使得所需的表面积分算子数量从传统矢量势公式中的十四个减少到了仅有的两个：经典的单层（$S$）和双层（$D$）势。
   • DEC 适配： 为了确保与简化的 SIE 兼容，对 DEC 框架进行了适配。不同于将矢量势 $\mathbf{A}$ 视为 1-形式（与网格边相关），本文将每个笛卡尔分量（$A_x, A_y, A_z$）视为独立的离散 0-形式（与网格节点相关）。这模仿了 DEC 中已建立的标量场处理方式，并与界面处的标量分量 SIE 表示法自然对齐。
4. 耦合机制： 混合系统通过强制执行势能及其法向导数的连续性，将内部 DEC 方程与边界 SIE 耦合。生成的线性系统具有稀疏且良置的矩阵。

主要贡献

• 算子缩减： 主要贡献是将表面积分算子集从 14 个减少到 2 个。这是通过限制分析对象为非磁性材料（具有空间均匀的磁导率）并利用矢量势的标量分量表示法来实现的。
• 保结构混合化： 本研究通过将 DEC 的 $\mathbf{A}$ 表示法适配为 0-形式，成功统一了 DEC 和 SIE。这避免了将基于边的 1-形式量与表面积分算子进行匹配的复杂性，在简化界面耦合的同时，保留了 DEC 框架底层的几何与拓扑结构。
• 低频稳定性： 通过采用基于势能的公式，该方法在广泛的频率范围内（包括传统方法失效的低频区域）保持了稳定性和准确性。

数值结果
所提方法通过三个数值实验进行了验证：

1. 介质球体： 该方法展示了随网格细化而收敛的特性，在细网格下实现了低于 1% 的相对 $L_2$ 误差。它表现出对高介电常数对比度（高达 $\epsilon_r = 45$）的鲁棒性，并保持了良好的数值条件（条件数 $\sim 10^3$）。
2. 多层球体： 对于具有多个介电层的球体，观察到了类似的收敛性和准确性。
3. 复杂 3D 散射体： 分析了一个由四个椭圆介质圆柱体组成介质板的结构。通过消去定理（Extinction Theorem）验证了准确性，结果证实了从边界数据重建的场能够正确抵消散射体内部的入射场，相对消去误差在 $10^{-4}$ 到 $10^{-5}$ 数量级。
4. 宽带稳定性： 研究表明，当频率趋于零（$k_0 a \to 0$）时，该方法仍能保持稳定，条件数保持有界，证实了其对低频崩溃的免疫力。

意义与主张
本文声称提供了一个物理一致、保结构且计算高效的框架，用于解决复杂几何形状和非均匀材料中的电磁散射与辐射问题。通过规避传统矢量势 SIE 的复杂性以及对人工领域截断的需求，该混合 DEC-SIE 方法为解决大规模问题提供了一种统一的方法。作者强调，该方法具有扩展到多物理场应用以及在超表面和开放共振结构中直接应用的潜力。这项工作强调，标量分量适配并非背离 DEC 原则，而是一种战略性的对齐，旨在增强兼容性和数值条件，同时不牺牲底层的几何结构。