  We study the convergence analysis of a finite element method for the
approximation of solutions to a seven-fields formulation of a
magnetohydrodynamics model, which preserves the energy of the system, and the
magnetic and cross helicities on the discrete level.
