Stable finite element methods preserving \(\nabla \cdot \varvec{B}=0\) exactly for MHD models

详细信息    查看全文

• 作者：Kaibo Hu ; Yicong Ma ; Jinchao Xu
• 关键词：Mathematics Subject Classification65N30 ; 65M60 ; 76W05
• 刊名：Numerische Mathematik
• 出版年：2017
• 出版时间：February 2017
• 年：2017
• 卷：135
• 期：2
• 页码：371-396
• 全文大小：
• 刊物类别：Mathematics and Statistics
• 刊物主题：Numerical Analysis; Mathematics, general; Theoretical, Mathematical and Computational Physics; Mathematical Methods in Physics; Numerical and Computational Physics, Simulation; Appl.Mathematics/Comput
• 出版者：Springer Berlin Heidelberg
• ISSN：0945-3245
• 卷排序：135

文摘

This paper is devoted to the design and analysis of some structure-preserving finite element schemes for the magnetohydrodynamics (MHD) system. The main feature of the method is that it naturally preserves the important Gauss’s law, namely \(\nabla \cdot \varvec{B}=0\). In contrast to most existing approaches that eliminate the electrical field variable \(\varvec{E}\) and give a direct discretization of the magnetic field, our new approach discretizes the electric field \(\varvec{E}\) by Nédélec type edge elements for \(H(\mathrm {curl})\), while the magnetic field \(\varvec{B}\) by Raviart–Thomas type face elements for \(H(\mathrm {div})\). As a result, the divergence-free condition on the magnetic field holds exactly on the discrete level. For this new finite element method, an energy stability estimate can be naturally established in an analogous way as in the continuous case. Furthermore, well-posedness is rigorously established in the paper for the Picard linearization of the fully nonlinear systems by using the Brezzi theory. This well-posedness naturally leads to robust (and optimal) preconditioners for the linearized systems.