High order H(div) discontinuous Galerkin methods for MHD equations