A Framework Based on hp-Finite Elements for Solving the Time-Harmonic Maxwell Equations

技术摘要：一种基于 $hp$-有限元的时谐麦克斯韦方程组求解框架

问题陈述
精确计算求解时谐麦克斯韦方程组在计算电磁学领域面临着显著挑战。这些困难主要源于双旋算子（double-curl operator）的不定性、伪模（spurious modes）的存在，以及在处理高频振荡解时伴随的高昂计算成本。虽然标准有限元方法在处理不定波问题时经常难以克服“污染误差”（pollution error），但本研究旨在开发一个既具备数学严谨性又兼顾计算效率的稳健框架，特别是在处理无界区域和复杂场行为时。

方法论
作者开发了一种稳定的 $hp$-有限元方法（FEM），旨在求解时谐麦克斯韦方程组中的双旋算子。该方法围绕以下几个核心组件构建：

1. 数学公式化： 研究始于通过对瞬态方程进行拉普拉斯变换导出的时谐麦克斯韦系统。作者引入了一种稳定的弱形式，该形式通过弱方式引入散度自由条件，从而允许使用标准的节点单元。该公式利用 Engquist-Majda 吸波边界条件（一阶）作用于特定的边界子集，以处理无界区域，并在远程边界上应用诺依曼（Neumann）条件。
2. 稳定性与适定性： 研究通过包含加权 $L^2$、$H^1$ 和边界范数的三范数进行稳定性分析。作者通过应用 Lax–Milgram 定理证明了离散解的存在性与唯一性，确立了相关双线性型的强制性和连续性。
3. 误差分析： 本文建立了先验（a priori）和后验（a posteriori）误差估计。
   • 先验估计提供了理论上的离散误差界限，展示了依赖于网格尺寸 $h$ 和多项式阶数 $p$ 的最优衰减率。
   • 后验估计基于残差概念，定义了由连续形式与近似解之间的差异所衍生的误差指示器。这些估计将误差限制在局部残差和边界跳跃之内。
4. 自适应策略： 基于残差误差估计，开发了一套完整的自适应实现方案。该策略的核心是一个自动选择机制，用于在 $h$-细化（网格剖分）与 $p$-细化（多项式增强）之间进行切换。这种选择由局部光滑度指示器驱动，该指示器通过评估勒让德（Legendre）展开系数的衰减情况来进行判断。如果解在局部是光滑的，则应用 $p$-细化；否则，则采用 $h$-细化。

主要贡献

• 稳定化 $hp$-框架： 本文提出了一种用于求解时谐麦克斯韦方程组的稳定化有限元公式，该公式能有效处理问题的不定性，并以弱方式强制执行散度自由条件。
• 理论保证： 该工作提供了严谨的数学基础，包括关于稳定性、唯一性以及针对先验和后验误差估计的最优收敛率的证明。
• 实用的自适应实现： 作者实现了一种实用的自适应策略，详细阐述了动态平衡 $h$-细化与 $p$-细化的算法。该方法旨在将计算资源集中在最需要的地方——即低正则性区域（需要网格细化）或高正则性区域（需要多项式增强）。
• 数值验证： 该框架通过三个数值算例进行了验证，这些算例涉及各种源项、材料系数（包括非均匀和奇异情况）以及区域配置。

结果
在具有不同材料系数（$\epsilon$ 和 $\sigma$）及源项的计算区域内进行了数值实验。跨多个测试用例的结果表明：

• 单调收敛： 在 30 次自适应迭代中，后验误差估计器（$\eta$）和绝对 $L^2$ 误差（$e$）均呈现单调下降趋势。
• 高效细化： 自适应策略成功地将细化过程集中在具有大解梯度或奇异性的区域。单元数量（$N_e$）和自由度数量（$N_{dof}$）虽呈递增趋势，但增长具有选择性，避免了均匀细化。
• 指数收敛： 误差衰减表现出相对于自由度立方根（$N_{dof}^{1/3}$）的指数收敛特性。这通过半对数图上的近线性衰减得到了证实，确认了该方法能够实现 $hp$-FEM 在解析解下的最优收敛率。
• **动态 $hp$-选择：** 最大多项式阶数（$\max(p)$）的演化呈现出阶梯式上升，证实了算法能够正确识别用于$p$-增强的光滑区域和用于$h$-细化的不规则区域。

意义
本文声称，与传统的用于时谐电磁场的有限元方法相比，所提出的方法在精度和效率方面都有显著提升。通过将稳定化公式与自适应 $hp$-策略相结合，该框架有效地缓解了污染误差，并能在不产生过度计算成本的情况下处理复杂的场行为。这项工作弥合了理论分析（稳定性与误差估计）与实际应用（自适应算法）之间的鸿沟，为求解涉及无界区域和非均匀材料的麦克斯韦方程组提供了一个稳健的工具。结果证实，$h$-细化与$p$-细化之间的平衡实现了相对于自由度的最优收敛，验证了该公式的数学严谨性。