[论文解读] A Mixed DG method and an HDG method for incompressible magnetohydrodynamics

Weifeng Qiu, Ke Shi|arXiv (Cornell University)|Feb 6, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 20被引用 3

一句话总结

本文提出了一种混合间断伽辽金（DG）方法和一种混合型间断伽辽金（HDG）方法，用于求解具有两种边界条件类型的稳态不可压缩磁流体动力学（MHD）问题。通过引入针对不连续多项式的新颖离散Sobolev嵌入估计，作者在最小正则性假设下，为所有变量——速度、磁场、压力和拉格朗日乘子——建立了最优先验误差估计，这是首次针对非线性MHD系统应用DG方法的此类分析。

ABSTRACT

In this paper we propose and analyze a mixed DG method and an HDG method for the stationary Magnetohydrodynamics (MHD) equations with two types of boundary (or constraint) conditions. The mixed DG method is based a recent work proposed by Houston et. al. for the linearized MHD. With two novel discrete Sobolev embedding type estimates for the discontinuous polynomials, we provide a priori error estimates for the method on the nonlinear MHD equations. In the smooth case, we have optimal convergence rate for the velocity, magnetic field and pressure in the energy norm, the Lagrange multiplier only has suboptimal convergence order. With the minimal regularity assumption on the exact solution, the approximation is optimal for all unknowns. To the best of our knowledge, this is the first a priori error estimates of DG methods for nonlinear MHD equations. In addition, we also propose and analyze the first divergence-free HDG method for the problem with several unique features comparing with the mixed DG method.

研究动机与目标

• 为具有两种不同边界条件或约束条件的稳态不可压缩MHD系统开发一种混合DG方法。
• 利用间断伽辽金方法，建立非线性MHD方程的先验误差估计，解决流体与磁场之间非线性耦合的挑战。
• 在精确解正则性要求最低的假设下，为所有未知量（速度、压力、磁场和拉格朗日乘子）实现最优收敛率。
• 提出并分析首个满足散度自由条件的HDG方法用于不可压缩MHD，确保局部守恒性并提升稳定性。
• 通过针对离散磁场$L^3$-范数的新离散Sobolev型估计，将现有线性化MHD DG分析推广至完整的非线性系统。

提出的方法

• 将先前针对线性化MHD工作的内部惩罚DG（IP-DG）格式，拓展至非线性稳态MHD系统。
• 提出两种针对不连续多项式逼近的新颖离散Sobolev嵌入型估计，以控制磁场的$L^3$-范数，这对处理非线性项至关重要。
• 采用混合有限元弱形式，引入混合变量：速度$\mathbf{u}$、压力$p$、磁场$\mathbf{b}$，以及用于散度自由约束的拉格朗日乘子$r$。
• 对于HDG方法，通过混合化消除单元内自由度，仅在网格骨架上保留全局耦合系统，从而提升计算效率。
• 对磁场逼近采用边、面和体积自由度，并重构一个全局类似连续的场$\tilde{\mathbf{b}}_h$，以促进误差分析。
• 应用形状正则性假设和柯西-施瓦茨不等式，界定离散磁场与其重构版本之间的差异，从而实现稳定性与收敛性分析。

实验结果

研究问题

• RQ1 在精确解正则性要求最低的假设下，能否在非线性不可压缩MHD系统的DG方法中为所有未知量实现最优收敛率？
• RQ2 如何为不连续多项式逼近构造离散Sobolev嵌入估计，以控制非线性MHD中磁场的$L^3$-范数？
• RQ3 对于具有两种不同磁场和拉格朗日乘子边界条件的非线性MHD混合DG方法，其收敛性特征如何？
• RQ4 能否为不可压缩MHD构造一种散度自由的HDG方法？其在稳定性与收敛性方面与混合DG方法相比如何？
• RQ5 DG与HDG格式中速度与磁场的局部守恒性对MHD模拟的精度与物理一致性有何影响？

主要发现

• 在精确解正则性要求最低的假设下，混合DG方法在能量范数下对速度、压力和磁场均实现了最优收敛率。
• 与磁场散度约束相关的拉格朗日乘子收敛阶次为次最优，但方法仍保持稳定且准确。
• 所提出的针对不连续多项式的离散Sobolev嵌入估计，能够有效控制磁场的$L^3$-范数，这对分析非线性耦合项至关重要。
• 首次提出适用于不可压缩MHD的散度自由HDG方法，保持了局部守恒性，并为混合DG方法提供了具有竞争力的替代方案。
• 分析结果确认，在相同正则性假设下，HDG方法对所有未知量均实现最优收敛，其收敛率与混合DG方法一致。
• 本工作首次为DG方法应用于非线性稳态不可压缩MHD方程提供了先验误差估计，填补了该领域的重要文献空白。