“…For examples, Gunzburger et al
[12] showed the existence, uniqueness and finite element approximations for MHD equations in three‐dimension‐domain; Johanna
[22] investigated a finite difference scheme for incompressible magnetohydrodynamics equations in two‐dimensional case; Wang et al
[30] developed a new two‐level finite element algorithm based on the Newton iterative method for the two‐ and three‐dimensional cases; Gao and Qiu
[8] presented a semi‐implicit finite element method and established the optimal error estimates under low assumption on the exact solutions and domain geometries; Adler et al
[1] investigated multigrid preconditioners based on specialized relaxation schemes for MHD system; Zhang et al
[33, 34] developed fractional‐step algorithms to construct fully decoupled, linear and unconditionally energy stable time discretization scheme; Wang et al
[28, 29] proposed unconditionally energy stable schemes by using a modified Crank–Nicolson method and established rigorous error estimates for MHD system and Cahn–Hilliard‐MHD system. Other numerical schemes and analysis for the MHD system with constant density could be founded in
[15, 18] and references cited therein. However, in addition to the above mentioned difficulties, there are very few works to construct efficient numerical algorithm for MHD equations with variable density due to lack of the maximum principle for
in hyperbolic equation.…”