we investigate the system of partial differential equations used in resistive magnetohydrodynamic modeling of fusion plasmas. this system couples the euler and maxwell equations for evolution of a charged fluid in an electromagnetic field, hence the magnetic field in the resulting pde system must evolve on a divergence-free constraint manifold. as traditional numerical solution approaches often violate these constraints, we investigate a reformulation of the resistive mhd system to allow for accurate evolution of the continuum-level equations, while simultaneously ensuring that the solution satisfies the solenoidal constraint.