Comment on “Heat transfer in MHD viscoelastic boundary layer flow over a stretching sheet with thermal radiation and non-uniform heat source/sink”

•Solutions to energy equation in terms of Kummer’s function are not correct.•Homotopy analysis method is applied to viscoelastic boundary layer flow, heat transfer.•Graphical and tabular demonstration of convergence of the HAM solutions is provided.•Effects of parameters on skin friction coefficient and wall heat transfer are tabulated.•Homotopy perturbation method does not yield convergent solutions.

In this paper, we demonstrate that previously reported analytical solutions for the temperature field given in terms of Kummer’s function by Nandeppanavar et al. (2011) [1], are incorrect. We then provide valid solutions of the governing ordinary differential equations for the fluid flow and temperature field using the homotopy analysis method (HAM) for two general types of non-isothermal boundary conditions, namely, prescribed surface temperature and prescribed heat flux. Our analysis is supported by a graphical and tabular demonstration of convergence of the HAM solutions.