AbstractThis study analyzed the MHD boundary layer flow of Eyring-Powell nanofluid past stretching cylinder with Cattaneo-Christov heat flux model. The governing non-linear partial differential equations corresponding to the momentum, energy and concentration have been formulated and transformed into a set of non-linear ordinary differential equations by using similarity transformations. Then the resulting non-linear high order ordinary differential equations of momentum, energy and concentration, subjected to boundary conditions were solved numerically by utilizing the second-order accurate implicit finite difference method known as Keller-Box which is programmed in the MATLABR2017b software. The results indicated that the velocity profile increases as the Eyring-Powell fluid parameterMand the curvature parameterγincrease but it decreases as the magnetic parameterHaincreases. Both the temperature and the concentration profiles have revealed an increment pattern for large values of the magnetic parameterHaand the thermophoresis parameterNtbut a decrement manner with increasing values of the Eyring-Powell fluid parameterM. The Brownian motion parameterNbhas shown an opposite influence on the temperature and the concentration profiles. The results also depicted that the local skin friction coefficient increases with increasing in Eyring-Powell fluid parameterM, magnetic parameterHa. Besides, it is found that both the local Nusselt numberNuxand the local Sherwood numberShxare higher for large vales of Eyring-Powell fluid parameterMand curvature parameterγ. Furthermore, the present results for the local skin friction coefficient, the local Nusselt number and the local Sherwood number are validated with the data of previously published literature for various limiting conditions where a very sound agreement has been attained.