The present paper studies the flow and heat transfer of the hybrid nanofluids flows induced by a permeable power-law stretching/shrinking surface modulated orthogonal surface shear. The governing partial differential equations were converted into non-linear ordinary differential equations by using proper similarity transformations. These equations were then solved applying a numerical technique, namely bvp4c solver in MATLAB. Results of the flow field, temperature distribution, reduced skin friction coefficient and reduced Nusselt number were deduced. It was found that increasing mass flux parameter slows down the velocity and, hence, decreases the temperature. Furthermore, on enlarging the stretching parameter, the velocity and temperature increases and decreases, respectively. In addition, that the radiation parameter can effectively control the thermal boundary layer. Finally, the temperature decreases when the values of the temperature parameter increases. We apply similarity transformation in order to transform the governing model into a system of ODEs (ordinary differential equations). Numerical solutions for particular values of involved parameters are in very good agreement with previous calculations. The most important and interesting result of this paper is that for both the cases of shrinking and stretching sheet flows exhibit dual solutions in some intervals of the shrinking and stretching parameter. In spite of numerous published papers on the flow and heat transfer over a permeable stretching/shrinking surface in nanofluids and hybrid nanofluids, none of the researchers studied the present problem. Therefore, we believe that the results of the present paper are new, and have many industrial applications.
Keywords: critical points; heat transfer; hybrid nanofluid; numerical methods; power law; stretching/shrinking sheet.