MFEM项目中求解耦合形式麦克斯韦方程的技术探讨

引言

在电磁场数值模拟领域，麦克斯韦方程组的求解是一个核心问题。本文基于MFEM项目中的一个典型问题，探讨了如何正确处理耦合形式的麦克斯韦方程组，特别是当电流密度与电场强度之间存在复杂关系时的数值求解策略。

问题描述

考虑以下耦合形式的麦克斯韦方程组：

1. 旋度方程： ∇×μ⁻¹∇×E = -∂J/∂t - ∂Js/∂t

2. 本构关系（欧姆定律）： J = σ̅E

其中Js表示外加激励源。这种耦合形式在考虑色散介质时尤为重要，因为直接代入欧姆定律会限制后续扩展的灵活性。

有限元离散化方案

在MFEM框架下，通常采用以下离散化方案：

• 电场E：使用H(curl)有限元空间（Nédélec元）
• 电流密度J：理论上适合使用H(div)空间

然而，实际计算表明，当电流密度J采用H(curl)空间离散时能得到正确结果，而采用H(div)空间则会出现问题。这一现象引发了深入的数值分析思考。

数值分析关键点

1. 质量矩阵等价性问题： 设M11(σ)表示H(curl)空间中与σ相关的质量矩阵，M12是H(curl)和H(div)空间的混合双线性形式，M22是H(div)空间的单位质量矩阵。一个重要的问题是：M11(σ)是否等于M12M22⁻¹M21(σ)？如果不等，则说明分离求解会引入投影/重插值误差。

2. 时间离散策略： 采用后向欧拉法进行时间离散时，分离求解的方式实际上将线性问题转化为"非线性"问题，这要求使用更小的时间步长，相当于显式方法的限制。

改进求解策略

针对这类耦合问题，推荐以下解决方案：

1. 整体求解方案： 构建块算子/矩阵系统，同时求解电场和电流密度，避免分离求解带来的误差积累。

2. 替代公式： 考虑使用磁扩散方程形式： ∇×σ⁻¹∇×B = μ∂B/∂t 这种形式避免了电导率的时间导数和质量项中可能的各向异性因子，计算电流密度时可通过B的旋度离散获得。

3. 非线性处理： 当本构关系确实呈现非线性时，可利用MFEM的NonlinearForm功能，结合非线性求解器进行处理。

结论

在MFEM框架下求解耦合形式的麦克斯韦方程组时，需要注意：

1. 谨慎选择各物理量的离散空间，H(div)空间对电流密度的离散可能引入不可忽视的误差
2. 优先考虑整体耦合求解而非分离迭代
3. 对于时变问题，采用隐式时间离散方法更有利于稳定性
4. 考虑问题物理特性，选择最合适的方程形式

这些经验对于处理更复杂的电磁问题，如色散介质中的波传播等，具有重要的参考价值。