In this paper, we develop and analyze a finite element projection method for magnetohydrodynamics equations in Lipschitz domain. A fully discrete scheme based on Euler semi-implicit method is proposed, in which continuous elements are used to approximate the Navier–Stokes equations and H(curl) conforming Nédélec edge elements are used to approximate the magnetic equation. One key point of the projection method is to be compatible with two different spaces for calculating velocity, which leads one to obtain the pressure by solving a Poisson equation. The results show that the proposed projection scheme meets a discrete energy stability. In addition, with the help of a proper regularity hypothesis for the exact solution, this paper provides a rigorous optimal error analysis of velocity, pressure and magnetic induction. Finally, several numerical examples are performed to demonstrate both accuracy and efficiency of our proposed scheme.


Mathematics and Statistics


National Natural Science Foundation of China, Grant 11871467

Keywords and Phrases

Error Analysis; Finite Element Method; Magnetohydrodynamics; Nédélec Edge Element; Projection Methods

International Standard Serial Number (ISSN)

1098-2426; 0749-159X

Document Type

Article - Journal

Document Version

Final Version

File Type





© 2023 Wiley, All rights reserved.

Publication Date

01 Mar 2023