Error Analysis of a PFEM Based on the Euler Semi-Implicit Scheme for the Unsteady MHD Equations

In this article, we mainly consider a first order penalty finite element method (PFEM) for the 2D/3D unsteady incompressible magnetohydrodynamic (MHD) equations. The penalty method applies a penalty term to relax the constraint “ [Formula: see text] ”, which allows us to transform the saddle point problem into two smaller problems to solve. The Euler semi-implicit scheme is based on a first order backward difference formula for time discretization and semi-implicit treatments for nonlinear terms. It is worth mentioning that the error estimates of the fully discrete PFEM are rigorously derived, which depend on the penalty parameter [Formula: see text] , the time-step size [Formula: see text] , and the mesh size h. Finally, two numerical tests show that our scheme is effective.