LAUSR

Abstract The MHD equation has some applications in physics and engineering. The main aim of the current paper is to propose a new numerical algorithm for solving the MHD equation. At first, the temporal direction has been discretized by the Crank–Nicolson scheme. Also, the unconditional stability and convergence of the time-discrete scheme have been investigated by using the energy method. Then, an improvement of element free Galerkin (EFG) i.e. the interpolating element free Galerkin method has been employed to discrete the spatial direction. Furthermore, an error estimate is presented for the full discrete scheme based on the Crank–Nicolson scheme by using the energy method. We prove that convergence order of the numerical scheme based on the new numerical scheme is O ( τ 2 + δ m ) . In the considered method the appeared integrals are approximated using Gauss Legendre quadrature formula. Numerical examples confirm the efficiency and accuracy of the proposed scheme.