## 1 Introduction  
The boundary layer theory has been a subject of interest for many researchers because of its significance in many industrial and engineering applications, i.e., airfoil design of of airplanes, the automobile industry and friction drag of a ship. The analytical solution on the boundary layer flow of a viscous fluid was introduced by [PERSON] [1] and [PERSON] and [PERSON] [2] for stretching and shrinking cases, respectively. A stretching sheet solution would generate a far field suction towards the sheet but opposite flow behaviors were observed for the shrinking case. The shrinking sheet seemed to distract the velocity away from the sheet. An adequate wall mass suction is needed to induce the flow due to the unconfined vorticity within the boundary layer. [PERSON] [3] analyzed the stagnation point flow over a shrinking sheet and discovered that the solutions were non-unique and did not exist for larger shrinking rates. [PERSON] et al. [4] analytically investigated viscous flow over shrinking rate with a slip condition and also concluded that the imposition of wall mass suction would induce two solutions. [PERSON] et al. [5] also concluded that the application of suction at certain values of the involving parameters contributed to the appearance of the second solution. Very recently, [PERSON] et al. [6] obtained dual solutions for mixed convective flow of a viscous fluid with copper nanoparticles over a permeable shrinking cylinder. The findings showed that the mixed effects of opposing flow and wall mass suction might generate two solutions. However, energy generation becomes a major issue in the industrial requirements. Attention was focused on enhancing the heat transfer rate of the systems, which include power stations, chemical plants, air conditioning and petrochemical industry. Nanofluid is created by dispersing the nanoparticles (size \(<\) 100 nm) into a base fluid such as water, ethylene or propylene glycol to inflate the thermal conductivity. Traditional Newtonian fluids such as water and air have been utilized as cooling fluids in many industrial and manufacturing processes but due to the final inspection of the product quality, it does not seem that the viscous fluid is a good choice for the cooling fluid. [PERSON] et al. [7] reported that the intensification of the working fluid's thermal conductivity can increase the thermal and electrical as well as the heat transfer efficiencies. The invention of nanofluids with great thermophysical properties can improve the performance for massive engineering operations such as in nuclear cooling systems, solar thermal energy, biomedical applications, lubrication, coolant in automobile radiators and others [8; 9; 10; 11; 12; 13; 14; 15]. Two remarkable nanofluid models were developed by [PERSON] [16] and [PERSON] and [PERSON] [17]. [PERSON] [16] highlighted seven slip mechanisms which can contribute relative velocity between the nanoparticles and the working fluids; however, only Brownian motion and thermophoresis are significant to model the nanofluids. [PERSON] and [PERSON] [17] proposed a model that investigated the effect of nanoparticles volume fraction on the enhancement of the thermal properties. [PERSON] and [PERSON] [18] initiated the investigation on the flow of a nanofluid towards a stretching sheet by using [PERSON]'s model. Latterly, [PERSON] et al. [19] explored the nanoparticles shape effect (brick, blade, cylindrical, platelet and spherical) on the boehmite alumina nanofluid flow characteristics. There are also a few boundary layer studies that modeled the non-Newtonian fluid as the operating fluid with the presence of nanoparticles as reported by [PERSON] et al. [20], [PERSON] and [PERSON] [21], [PERSON] and [PERSON] [22] and [PERSON] et al. [23]. Most of the boundary layer flow and heat transfer problems considered the classical Fourier model which represents the thermal conduction heat transfer.