An analysis has been carried out to study a problem of the boundary-layer flow of Carreau nanofluid over a non-linearly stretched sheet with chemical reaction and the heat generation/absorption in a porous medium. A power-law model includes a two-phase model for Carreau nanofluid with a convective condition. The governing PDEs with the corresponding boundary conditions are modified to a system of non-linear ODEs with the appropriate boundary conditions by picking local similarity conversions and solved numerically by using Runge–Kutta–Fehlberg 4th–5th order numerical method (RKF45) on based shooting technique. This investigation discusses the effects of study parameters like the porosity parameter K 1 , the heat source λ > 0 or sink λ 0 parameter, the chemical reaction parameter γ 1 and the Biot number B i on flow velocity, temperature and nanofluid volume fraction in addition to the heat and mass transfer rates tabular and graphically. A comparative study is likewise revealed showing the comparison of current results with previously published data.