Convergent finite element discretization of the multi-fluid nonstationary incompressible magnetohydrodynamics equations