[论文解读] Error estimates for structure-preserving discretization of the incompressible MHD system

Yicong Ma, Jinchao Xu+1 more|arXiv (Cornell University)|Aug 10, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 25被引用 4

一句话总结

本文针对不可压缩磁流体动力学（MHD）系统，提出了一种保持结构的有限元离散化方法的最优误差估计，结合能量稳定性和物理守恒律。在适当的CFL条件下，该方法在空间和时间上均达到最优收敛阶，数值实验验证了收敛率与理论预测一致。

ABSTRACT

In this paper, we carry out the error analysis for the structure-preserving discretization of the incompressible MHD system. This system, as a coupled system of Navier-Stokes equations and Maxwell's equations, is nonlinear. We use its energy estimate and the underlying physical structure to facilitate the error analysis. Under certain CFL conditions, we prove the optimal order of convergence. To support the theoretical results, we also present numerical tests.

研究动机与目标

• 建立应用于不可压缩MHD系统的保持结构的有限元方法的严格误差估计。
• 分析基于H(curl)和H(div)协调元的混合有限元格式的收敛行为，该格式可保持关键物理定律（如磁场的高斯定律）。
• 在特定时间步长和网格尺寸约束下，证明非线性、时间依赖的MHD方程的最优收敛阶。
• 通过使用制造解和收敛率分析的数值实验验证理论结果。

提出的方法

• 该方法采用基于有限元外微分形式（FEEC）和离散外微分形式（DEC）的混合有限元格式，使用H(curl)-协调的Nédélec元模拟电场，使用H(div)-协调的Raviart-Thomas元模拟磁场。
• 该格式在离散层面精确保持磁场的散度为零，确保满足高斯定律。
• 采用半隐式时间格式，非线性对流项显式处理，线性项隐式处理，从而在CFL条件下保持稳定性。
• 通过基于MHD系统底层物理结构的能量估计进行误差分析，包括总能量守恒。
• 分析依赖于Sobolev范数中的先验估计以及有限元空间的逼近性质，通过速度和磁场梯度的有界性对非线性项进行仔细处理。
• 收敛框架源自演化鞍点问题，将现有理论扩展至具有结构保持特性的非线性、时间依赖系统。

实验结果

研究问题

• RQ1 在时间和空间离散化下，保持结构的有限元离散化方法对不可压缩MHD系统的最优收敛率是多少？
• RQ2 时间步长 $ k $ 和网格尺寸 $ h $ 如何影响速度、磁场、电场和压力变量的收敛阶？
• RQ3 能否利用能量稳定性和物理守恒律（如 $ \nabla \cdot \bm{B} = 0 $）推导出非线性、时间依赖MHD系统的严格误差估计？
• RQ4 所提出的方法是否在合理的CFL条件下实现最优收敛阶？与现有方法相比表现如何？
• RQ5 数值实验在多大程度上验证了误差分析所预测的理论收敛率？

主要发现

• 在条件 $ k \leq h^{1/2} $ 下，速度、磁场和电场在各自能量范数下的收敛阶达到最优，为 $ \mathcal{O}(k^2 + h^{2\hat{r}}) $。
• 当 $ k \leq h^{3/2} $ 时，压力的收敛阶达到最优 $ \mathcal{O}(k^2 + h^{2r}) $，表明压力精度对时间步长有更严格的限制。
• 电场误差满足 $ \| \bm{E} - \bm{E}_h \|^2 + k \| \nabla \times (\bm{E} - \bm{E}_h) \|^2 \leq C(k^2 + h^{2\hat{r}}) $，在相同CFL条件下确认了最优收敛性。
• 数值实验验证了预测的收敛率：对网格尺寸 $ h $ 和时间步长 $ k $，观测到的收敛阶均约为1，与理论分析一致。
• 该格式在离散层面精确保持磁场的散度为零，有助于提高解的稳定性和精度。
• 理论误差界依赖于精确解的正则性，包括 $ \| \bm{u} \|_{0,\infty} $、$ \| \nabla \bm{u} \|_{0,3} $、$ \| \bm{B} \|_{0,\infty} $ 和 $ \| \bm{j} \|_{0,\infty} $ 的有界性，以及区域几何形状。