  A new MHD solver, based on the Nektar++ spectral/hp element framework, is
presented in this paper. The velocity and electric potential quasi-static MHD
model is used. The Hartmann flow in plane channel and its stability, the
Hartmann flow in rectangular duct, and the stability of Hunt's flow are
explored as examples. Exponential convergence is achieved and the resulting
numerical values were found to have an accuracy up to $10^{-12}$ for the state
flows compared to an exact solution, and $10^{-5}$ for the stability
eigenvalues compared to independent numerical results.
