Homotopy perturbation method for MHD non-Newtonian Williamson fluid over exponentially stretching sheet with viscous dissipation and convective boundary condition

This research is aimed at presenting the two-dimensional steady fluid flow, represented by Williamson constitutive model past a nonlinear exponential stretching sheet theoretically. The system of ODEs describing the physical problem is successfully solved numerically with the help of the homotopy perturbation method (HPM). Special attention is given to study the convergence analysis of the proposed method. The influences of the physical governing parameters acting on the fluid velocity and the fluid temperature are explained with the help of the figures and tables. Further, the presented numerical method is employed to calculate both the rate of heat transfer and the drag force for the Williamson fluid flow. In particular, it is observed that both the Eckert number and the dimensionless convective parameter have the effect of enhancing the temperature of the stretching surface, while the inverse was noted for the dimensionless mixed convection parameter. Finally, the comparison with previous numerical investigations of other authors at some special cases which is reported here proves that the results obtained via homotopy perturbation method are accurate and the numerical method is reliable.