A stabilized element-free Galerkin (EFG) method is proposed in this paper for numerical analysis of the generalized steady MHD duct flow problems at arbitrary and high Hartmann numbers up to \(10^{16}\) . Computational formulas of the EFG method for MHD duct flows are derived by using Nitsche’s technique to facilitate the implementation of Dirichlet boundary conditions. The reproducing kernel gradient smoothing integration technique is incorporated into the EFG method to accelerate the solution procedure impaired by Gauss quadrature rules. A stabilized Nitsche-type EFG weak formulation of MHD duct flows is devised to enhance the performance damaged by high Hartmann numbers. Several benchmark MHD duct flow problems are solved to testify the stability and the accuracy of the present EFG method. Numerical results show that the range of the Hartmann number Ha in the present EFG method is \(1\le Ha\le 10^{16}\) , which is much larger than that in existing numerical methods.
