A convergent linearized lagrange finite element method for the magneto-hydrodynamic equations in two-dimensional nonsmooth and nonconvex domains

A new fully discrete linearized H1-conforming Lagrange finite element method is proposed for solving the two-dimensional magneto-hydrodynamics equations based on a magnetic potential formulation. The proposed method yields numerical solutions that converge in general domains that may be nonconvex, nonsmooth, and multiconnected. The convergence of subsequences of the numerical solutions is proved only based on the regularity of the initial conditions and source terms without extra assumptions on the regularity of the solution. Strong convergence in L2(0, T;L2(Ω)) was proved for the numerical solutions of both u and H without any mesh restriction.