In this paper, Optimal Homotopy Asymptotic Method (OHAM) is combined with Galerkin
method to investigate the MHD Jeffery-Hamel flows in non-parallel walls for strongly nonlinear ordinary
differential equations. Results for velocity profiles in divergent and convergent channels for various values
of Hartmann and Reynolds numbers are compared with numerical method (the Runge-Kutta method of
order 4). The results confirm the notion that the OHAM method is a powerful and efficient technique for
finding exact solutions for differential equations which have great significance in different fields of science
and engineering. Also, it can be found that increasing Hartmann number will lead to backflow reduction.
