Numerical approximation for the MHD equations with variable density based on the Gauge-Uzawa method

In this paper, we consider the numerical approximation of incompressible magnetohydrodynamic (MHD) system with variable density. Firstly, we provide first- and second-order time discretization schemes based on the convective form of the Gauge-Uzawa method. Secondly, we prove that the proposed schemes are unconditionally stable. We also provide error estimates through rigorous theoretical analysis. Then, we construct a fully-discrete first-order scheme with finite elements in space and provide its stability result. Finally, we present some numerical experiments to validate the effectiveness of the proposed schemes. Furthermore, we also present the conserved scheme and its numerical results.