“…In order to investigate the physical representation of the problem, the numerical values of velocity (f / (η)), temperature (θ(η)) and species concentration (φ (η)) with the boundary layer have been computed for different parameters as Williamson fluid parameter, We, Grashof number, Gr, modified Grashof number, Gc, Hartmann number, M, velocity ratio parameter, a/c, Prandtl number, Pr, thermal radiation, R, heat source, Q, Eckert number, Ec, Schmidt number, Sc and Soret number, So. In order to assess the accuracy of the numerical results the present results are compared with the previous studies (Fang et al, 2009;Nadeem et al, 2013;Ramzan et al, 2017;Salahuddin et al, 2016;Sher Akbar et al, 2015) presented in Tables 1-3 The velocity profile for different values of Williamson fluid parameter, We, are given in Fig. 3 which shows a decreasing trend in momentum boundary layer due to increase in We.…”