## 4 Discussion  
The impacts of emerging parameters have been presented in this section. The emerging parameters are stretching parameter \((c)\), Williamson parameter \((We)\), Porosity parameter \((\gamma)\), thermal radiation parameter \((Rd)\), Prandtl number \((\Pr)\), thermophoresis parameter \((Nt)\), Brownian motion parameter \((Nb)\) and Biot numbers \((Bi_{1},Bi_{2})\) on velocities profiles \((f^{\prime}(\eta),g^{\prime}(\eta))\), temperature profile \(\theta(\eta)\) and concentration profile \(\phi(\eta)\).Figure 1 shoe the physical sketch of the flow. Figures 2-9 are plotted to observe the impact of these emerging parameters. Figure 2a is presented to see the impression of stretching parameter \(c\) on \(f^{\prime}(\eta),g^{\prime}(\eta)\), \(\theta(\eta)\) and \(\phi(\eta)\). Physically the stretching parameter is the ratio of the fluid velocity along the x-direction to the fluid velocity along the y-direction. So this parameter is directly related to the velocity in the y-direction and related indirectly to the velocity in the x-direction. Because of this reason, the stretching parameter reduces \(f^{\prime}(\eta)\) in the x-direction, as shown in Figure 2a, and the stretching parameter enhances \(g^{\prime}(\eta)\) in the y-direction, as shown in Figure 2b. The influences of stretching parameter \(c\) on \(\theta(\eta)\) and \(\phi(\eta)\) are shown in Figure 2c,d. The increasing values of stretching parameter \(c\) show a reduction in \(\theta(\eta)\) and \(\phi(\eta)\). Figure 3a,b shows the impact of \((We)\) on \(f^{\prime}(\eta),g^{\prime}(\eta)\) in the x- and y-directions respectively. This parameter is the ratio of the relaxation to retardation time. The relaxation time is enhanced by increasing the Williamson parameter. Because of this, the fluid viscosity increases and as a result the velocity of fluid particles falls. Figure 4a,b shows the impact of \((\gamma)\) on \(f^{\prime}(\eta),g^{\prime}(\eta)\). Physically, the porous medium influences the boundary layer, which leads to acceleration of the fluid. Due to this fact, the velocity profile reduces with the enhancement of the porosity parameter. Figure 5 shows the impact of \((Rd)\) on \(\theta(\eta)\). We see from the figure that the enhancement of \((Rd)\) upsurges \(\theta(\eta)\). \((Rd)\) plays a significant role in the heat transmission process. When we increase thermal radiation, the temperature of the boundary layer flow enhances. Figure 6 depicts the impact of \((\mathrm{Pr})\) on \(\theta(\eta)\). Physically, with the increasing in the Prandtl number, the thermal diffusion diminishes, and also a small amplitude of oscillation in temperature occurs. Therefore, \(\theta(\eta)\) reduces with increasing \((\mathrm{Pr})\). Figure 7a, b depicts the impact of \((Nt)\) on \(\theta(\eta)\) and \(\phi(\eta)\). Here we see that the enhancement of \((Nt)\) reduces \(\theta(\eta)\) and \(\phi(\eta)\). In addition, thermophoresis is high in Newtonian fluids as compared to non-Newtonian fluids. Figure 8a,b depicts the impact of \((Nb)\) on \(\theta(\eta)\) and \(\phi(\eta)\). In Figure 8a, we see that \(\theta(\eta)\) enhances with the enhancement of \((Nb)\), while the opposite behavior of \((Nb)\) on \(\phi(\eta)\) can be seen in Figure 8b. Figure 9a-c depicts the impact of \((Bi_{1},Bi_{2})\) on \(\theta(\eta)\) and \(\phi(\eta)\). From the figure, we can see that the increase of Biot numbers \((Bi_{1},Bi_{2})\) causes an increase in the temperature and concentration profiles. The numerical values of the surface drag force and impact of the different parameters; \(We\), \(\gamma\) and \(c\), are shown in Table 1. It is observed that increasing \(c\) and \(\gamma\) reduces the surface drag force \(f^{\prime\prime}(0)\) and \(g^{\prime\prime}(0)\), while increasing \(We\) augments the surface drag force \(f^{\prime\prime}(0)\) and \(g^{\prime\prime}(0)\). Table 2 displays the numerical values of the heat transfer rate for the different embedded parameters; \(c,\mathrm{Pr},Rd,Nb\) and \(Nt\).