Error analysis of a divergence-preserving mixed finite element scheme for the incompressible Hall--magnetohydrodynamic equations

技术摘要：保散度混合有限元格式对不可压缩霍尔磁流体动力学方程的误差分析

问题陈述
不可压缩霍尔磁流体动力学（Hall–MHD）系统描述了导电流体中流体运动与磁场之间的相互作用，其中包含了在离子惯性长度尺度下变得显著的霍尔效应。该效应在感应方程中引入了一个刚性且高度非线性的项（$\eta \text{curl}(\text{curl} B \times B)$），导致色散行为、陡峭梯度以及高频振荡（哨波）。这些特征带来了巨大的分析和计算挑战。具体而言，三维空间中强解的整体存在性仍是一个未解决的问题，而标准数值方法往往难以兼顾稳定性以及磁场严格保持无散（$\text{div } B = 0$）的要求。违反该无散条件会诱导出非物理的磁单极子并引发严重的数值不稳定性。此外，现有的数值方法（例如基于标量辅助变量（SAV）的方法或谱技术）往往耗散修正的能量泛函而非真实的物理能量，无法在离散层面精确强制执行无散条件，或缺乏严格的误差分析和几何灵活性。

方法论
为应对这些挑战，作者提出了一种针对 Voigt 正则化霍尔–MHD 系统的保结构、线性、全离散混合有限元方法（FEM）。Voigt 正则化将动量方程中的时间导数替换为空间滤波版本 $(I - \alpha_1 \Delta)\partial_t u$，并在广义欧姆定律中保留电子惯性项 $\alpha_2 \partial_t J$。这一改进在保持原系统基本物理结构和能量定律的同时，提高了适定性和数值稳定性。

数值格式由以下组件构建：

1. 混合公式：问题被重新表述为以电场 $E$ 和电流密度 $J$ 为主要变量，从而得到动量方程、感应方程、欧姆定律和安培定律方程的耦合系统。
2. 空间离散化：该格式采用有限元外微积分（FEEC）以确保无散条件被精确满足。
   • 速度 $u$ 和压力 $p$ 使用 MINI 有限元对（最低阶 inf-sup 稳定）进行近似。
   • 磁场 $B$ 和电流密度 $J$ 分别使用最低阶 Raviart–Thomas（$RT_0$）和 Nédélec（$X_0$）空间进行近似。
   • 这些空间构成离散 de Rham 复形，确保 $\text{div } B_h = 0$ 在离散层面被精确保持。
3. 时间离散化：采用反对称、线性隐式时间步进格式。对流项和洛伦兹项利用前一时间步的值进行半隐式处理，以保持线性性和无条件能量稳定性。
4. 误差分析：作者通过将误差分解为近似误差（利用斯托克斯投影和正交投影）和离散误差，建立了最优收敛速率。他们利用离散 Gaffney 不等式和能量估计来处理非线性耦合项。

主要贡献

• 保结构格式：本文提出了一种线性、全离散有限元格式，精确强制执行磁场的无散约束，避免了可能降低精度或稳定性的投影步骤或惩罚方法。
• 无条件能量稳定性：所提出的格式无条件满足离散能量耗散律，即使在存在刚性霍尔项的情况下，也镜像了连续物理能量定律。
• 严格的误差估计：作者为 Voigt 正则化问题提供了完整的误差分析，在适当的正则性假设下建立了 $O(h^s + \tau)$ 阶的最优收敛速率。
• 推广至未正则化系统：至关重要的是，该分析扩展到了未正则化的霍尔–MHD 系统（其中 $\alpha_1 = \alpha_2 = 0$）。作者证明，只要电阻率 $\sigma > 0$，即可证明未正则化系统的收敛性。在非电阻机制（$\sigma = 0$）下，Voigt 正则化被证明在误差估计中起着关键作用。
• 2.5D 和 3D 实现：该框架在 2.5D（平移不变性）和完全 3D 设置中均得到验证，证明了该方法适用于复杂几何结构和物理现象。

结果
2.5D 和 3D 中的数值模拟证实了理论发现：

• 收敛性：该格式实现了速度、磁场和电流密度的预测时空收敛速率。虽然在某些情况下速度和电流密度的收敛速率略有降低（归因于霍尔项导致的解正则性有限），但磁场保持了最优收敛性。
• 散度保持：最大散度误差（$\max | \text{div } B_h |$）在整个模拟过程中保持在求解器容差水平，证实了无散约束的精确保持。
• 能量稳定性：系统总能量随时间单调衰减，与理论能量耗散律一致。
• 物理保真度：对 Arnold–Beltrami–Childress（ABC）流、Orszag–Tang 涡旋以及 Harris 片中的霍尔驱动重联的模拟，成功捕捉了复杂的物理现象，包括磁重联、磁岛形成、通量绳以及伴随霍尔–MHD 动力学的特征四极电流模式。即使在相对粗糙的网格上，该格式也能准确解析这些特征。

意义
本文声称，这项工作为模拟霍尔–MHD 系统提供了一个数学上严谨且计算上稳健的框架。通过将 Voigt 正则化与保散度混合有限元公式相结合，作者克服了霍尔项的非线性与严格的无散要求这双重挑战。该方法提供了一种物理上自洽的近似，具有稳定性、收敛性，并能处理复杂几何结构，填补了以往方法缺乏严格误差分析或精确散度保持的空白。结果表明，该方法是在霍尔效应占主导地位的天体物理和实验室等离子体现象研究中的可行工具。