In this paper, an Arrow-Hurwicz iterative finite element method based on charge-conservation for solving the stationary inductionless magnetohydrodynamic equations is proposed and analyzed. The current density and electric potential are discretized by \(\textbf{H}(\textrm{div},\varOmega )\times L^2(\varOmega )\) -conforming finite element pairs, respectively, which maintains charge-conservation. The hydrodynamic unknowns are discretized by stable velocity-pressure finite element pairs. Furthermore, we use the Arrow-Hurwicz iterative method to decouple the velocity and pressure, which leads to that there is no large-scale saddle point system that needs to be solved at each iteration step except for the determination of the initial functions. It has been proven that the proposed method converges geometrically with a contraction number independent of the mesh width. Numerical results are presented to verify the results of theoretical analysis and the effectiveness of the proposed method.
