The boundary elements method for magneto-hydrodynamic (MHD) channel flows at high Hartmann numbers