摘要

This paper characterizes one penalty finite element method for the incompressible MHD equations. The method is an interesting combination of the classic iterative schemes (Stokes, Newton and Oseen iterations) with two different finite element pairs P(1)b-P-1-P(1)b and P-1-P-0-P-1. Moreover, the rigorous analysis of stability and error estimate for the proposed methods are given. Finally, the applicability and effectiveness of the presented schemes are illustrated in several numerical experiments.