Stable finite element methods preserving ∇ · B= 0 exactly for MHD models

Kaibo Hu, Yicong Ma, Jinchao Xu

Research output: Contribution to journalArticlepeer-review

49 Scopus citations


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 ∇ · B= 0. In contrast to most existing approaches that eliminate the electrical field variable E and give a direct discretization of the magnetic field, our new approach discretizes the electric field E by Nédélec type edge elements for H(curl) , while the magnetic field B by Raviart–Thomas type face elements for H(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.

Original languageEnglish (US)
Pages (from-to)371-396
Number of pages26
JournalNumerische Mathematik
Issue number2
StatePublished - Feb 1 2017
Externally publishedYes

Bibliographical note

Funding Information:
This material is based upon work supported in part by the US Department of Energy Office of Science, Office of Advanced Scientific Computing Research, Applied Mathematics program under Award Number DE-SC-0014400 and by Beijing International Center for Mathematical Research of Peking University, China.

Publisher Copyright:
© 2016, Springer-Verlag Berlin Heidelberg.


  • 65M60
  • 65N30
  • 76W05


Dive into the research topics of 'Stable finite element methods preserving ∇ · B= 0 exactly for MHD models'. Together they form a unique fingerprint.

Cite this