Theoretical study of non-Newtonian micropolar nanofluid flow over an exponentially stretching surface with free stream velocity

Abstract

The computational analysis of the second-order micropolar stagnation point flow of nanofluid over an exponentially permeable stretching sheet is considered. The freestream velocity with the thermal slip effects is taken into account in this analysis. This model is developed on the basis of flow assumptions and reduced into partial differential equations before applying the boundary layer approximations. The governing equations as a mathematical model are simplified with the help of suitable transformations. The differential system is further solved by using the bvp4c. Both graphs and tables are used to report observations. The skin friction and Nusselt number are reported for both weak and strong concentrations. The magnitude of skin friction is noticed greater for strong concentration in comparison with weak concentration. Subject to couple stress, the values of skin friction are relatively high for the case of weak concentration in comparison with strong concentration. Micropolar profile admit the direct relation toward micropolar parameter and micro-gyration parameter. Both Sherwood number and Nusselt number admits higher values for strong concentration as compared to weak concentration.

Keywords

Non-Newtonian fluid micropolar fluid exponential stretching numerical technique slip effects

Introduction

The study of non-Newtonian fluids acknowledged notable attention by the researchers because of a wide range of applications. Owning to such importance as yet numerous non-Newtonian fluid models are introduced to encounter the flow field properties of fluids like the flow of an incompressible second-order fluid past a stretched sheet was considered by Rajagopal et al.¹

The issue involved continuous extrusion of a polymer sheet from a die, which is a polymer processing application. Bujurke et al.² investigated second order fluid flow along with heat transfer aspects by use of Noll and Coleman constitutive equations . This flow’s boundary layer properties were determined. Both thermal and velocity boundary layer thickness were closely examined. Pontrelli³ examined the steady fluid flow of a second-grade homogeneous incompressible fluid across a stretching sheet. The flow partial differential equations are solved by using the collocation method. The solution dependency on the non-Newtonian parameter was explored, and the shear stress was provided along with several velocity distribution. The two dimensional flow of a viscoelastic fluid past a moving sheet was studied by Sadeghy and Sharifi.⁴ They utilized the boundary layer approach for their investigation. HAM solution for second grade fluid over a porous plate was solved by Hayat and Khan.⁵ Fetecau and Fetecau⁶ discussed the solutions of second grade fluid for unsteady and unidirectional flows. Khan et al.⁷ investigated the physical aspects of third grade nanoparticles with optimization of entropy. Flow of a viscoelastic fluid toward a fixed plate with non-axisymmetric Homann stagnation point was studied by Mahapatra and Sidui.⁸ Shamshuddin et al.⁹ examined the impacts of power law fluid model under the viscous dissipation with Hall current at exponential stretching sheet. Shamshuddin et al.¹⁰ considered the impacts of Joule heating and viscous dissipation on magneto power law fluid model over exponential stretching sheet. Ferdows et al.¹¹ investigated the effects of radiative micropolar fluid flow with non Darcy porous medium having homogeneous–heterogeneous reactions. Dawar et al.¹² analyzed the effects of MHD micropolar nanofluid chemically over stretching surface. They also highlighted the effects of variable heat sink/source and velocity slips effects for micropolar nanofluid. The recent developments in this regard can be assess in References.^13–15

The fluid with nanosized particles is termed as nanofluid. Owning to importance of nanofluids various authors reported their finds like Rahman et al.¹⁶ discussed the suspended nanoparticles in the fluid flow past a porous exponentially stretching/shrinking by using the Buongiorno’s model. Thermal radiation, magnetic field effects with slip flow over a stretching/shrinking sheet with nanofluid was concluded by Abdul Hakeem et al.¹⁷ Viscous dissipation and second order slip of melting heat transfer on magneto-hydrodynamics convective flow of a nanofluid past a stretching sheet was studied by Mabood and Das.¹⁸ Kadir et al.¹⁹ investigated the effects of von Karman swirling bioconvection over a rotating disk with the nanofluid flow. Beg et al.²⁰ offered the experimental results about the lubricity and rheology of drilling liquids enriched having nanoparticles. Shamshuddin et al.²¹ examined the time dependent micropolar fluid flow over rotating desk under the multi-physico-chemical magnetic field. Numerical outcomes of nanofluid effects discussed by Kumar and Sokhal.²² Khan et al.²³ studied the features of activation energy and Wu’s slip of bioconvection with magnetized couple stress nanofluid. Nadeem et al.²⁴ examined the impacts of CNTs with base fluid at rotating disk. Abbas et al.²⁵ reported the impact of induced MHD micropolar fluid at curved surface. The recent studies on the MHD nanofluid flow under the various flow assumptions can be reached out in References.^26–29

We studied the non-Newtonian micropolar nanofluid flow with viscous dissipation over an exponentially stretching sheet. The thermal slip effects and suction/injection are also considered. The mathematical model is developed using the boundary layer approximation under the flow assumptions. The developed differential model is solved through the numerical technique namely bvp4c. Before our study, no one highlighted the effects on the boundary layer flow discussed. The ultimate results are shared by using tables and graphs.

Mathematical formulation

Non-Newtonian stagnation point micropolar fluid flow over an exponential stretching surface is considered. Both suction and injection are considered and specifically, $V_{w} < 0$ signifies the injection flow and $V_{w} > 0$ symbolizes the suction flow. The nanoparticles are considered suspended in the flow regime. Buongiorno fluid model is used in this analysis. The geometry of problem is given in Figure 1. The flow field is mathematically modeled and the ultimate flow differential system for micropolar fluid (Rehman et al.,³⁰ Dero et al.³¹) is concluded as follows:

\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0,

(1)

\begin{array}{l} u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} = U_{e} \frac{\partial U_{e}}{\partial x} + (\frac{μ + k}{ρ}) \frac{\partial^{2} u}{\partial y^{2}} \\ + \frac{α_{1}}{ρ} (u \frac{\partial^{3} u}{\partial x \partial y^{2}} + \frac{\partial u}{\partial x} \frac{\partial^{2} u}{\partial y^{2}} + v \frac{\partial^{3} u}{\partial y^{3}} + \frac{\partial u}{\partial y} \frac{\partial^{2} u}{\partial x \partial y}) + \frac{k}{ρ} \frac{\partial N}{\partial y} \end{array}

(2)

\begin{matrix} u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y} = α \frac{\partial^{2} T}{\partial y^{2}} + \frac{ρ C_{p}}{ρ C_{f}} [D_{B} \frac{\partial C}{\partial y} \frac{\partial T}{\partial y} + \frac{D_{T}}{T_{\infty}} {(\frac{\partial T}{\partial y})}^{2}] \\ + (\frac{μ + k}{ρ C_{f}}) {(\frac{\partial T}{\partial y})}^{2}, \end{matrix}

(3)

u \frac{\partial C}{\partial x} + v \frac{\partial C}{\partial y} = D_{B} (\frac{\partial^{2} C}{\partial y^{2}}) + \frac{D_{T}}{T_{\infty}} \frac{\partial^{2} T}{\partial y^{2}}

(4)

u \frac{\partial N}{\partial x} + v \frac{\partial N}{\partial y} = \frac{γ^{*}}{j ρ} \frac{\partial^{2} N}{\partial y^{2}} - \frac{k}{j ρ} \frac{\partial u}{\partial y} - \frac{2 k}{j ρ} N,

(5)

Figure 1.

Flow of micropolar nanofluid pattern.

here, equation (1) represents the continuity equation, equation (2) represents the momentum equation, equation (3) represents the energy equation, equation (4) represents the concentration equation, and equation (5) represents the micropolar equation respectively. The endpoint conditions of the fluid flow are defined as

\begin{matrix} v = V_{w}, u = U_{w}, T = T_{w} + L \frac{\partial T}{\partial y}, D_{B} \frac{\partial C}{\partial y} \\ + \frac{D_{T}}{T_{\infty}} \frac{\partial T}{\partial y} = 0, N = - n \frac{\partial u}{\partial y}, at y \to 0, \\ u = U_{\infty}, T = T_{\infty}, C = C_{\infty}, N = 0 as y \to \infty . \end{matrix}

(6)

The suitable transformations are as follows:

\begin{array}{l} u = u_{0} e^{x / L} F^{'} (η), v = - \sqrt{\frac{u_{0} ν}{2 L}} e^{x / 2 L} (F (η) + η F^{'} (η)), \\ T = T_{\infty} + (T_{w} - T_{\infty}) θ (η), C = C_{\infty} + (C_{w} - C_{\infty}) ϕ (η), \\ N = - \sqrt{\frac{{u_{0}}^{3}}{2 ν L}} e^{3 x / 2 L} H (η), η = \sqrt{\frac{u_{0}}{2 ν L}} e^{x / 2 L} y . \end{array}

(7)

The suitable transformations are applied on equations (1)–(5) and we have,

\begin{matrix} (1 + K) F^{″'} + β_{1} (2 η F^{″} F^{″'} + 5 F^{'} F^{″'} + 3 F ″ F^{″} - F F^{iv}) \\ + K H^{'} - 2 F^{'} F^{'} + F F^{″} = 0, \end{matrix}

(8)

\frac{1}{P_{r}} θ^{″} + N_{b} {θ^{'}}^{2} + N_{t} θ^{'} ϕ^{'} - 2 θ F^{'} + θ^{'} F + Ec (1 + K) F ″ F ″ = 0,

(9)

ϕ^{″} + \frac{N_{t}}{N_{b}} θ^{″} - Sc (2 F' ϕ - ϕ^{'} F) = 0,

(10)

(1 + \frac{K}{2}) H^{″} + F H^{'} - 3 F^{'} H - K (2 H + F ″) = 0 .

(11)

The reduced conditions are

\begin{matrix} F (0) = S, θ (0) = 1 + λ θ^{'} (0), \\ F^{'} (0) = γ, N_{b} ϕ^{'} (0) + N_{t} θ^{'} (0) = 0, \\ H (0) = - n F^{″} (0), \\ F^{'} (\infty) = 1, θ (\infty) = 0, ϕ (\infty) = 0, H (\infty) = 0, \end{matrix}

(12)

where, $F (η)$ and $θ (η)$ are the velocity and temperature of the fluid. $ϕ (η)$ and $H (η)$ are micropolar and concentration profiles. $S > 0$ presented the suction and $S < 0$ represents the injection. The shear stress and wall heat stress on the permeable sheet are defined as follows:

\begin{matrix} C_{fx} = \frac{τ_{wx}}{ρ_{f} U_{w}^{2}}, C_{m} = \frac{τ_{m}}{ρ_{f} u_{w}^{2}}, N u_{x} = \frac{x q_{w}}{k_{f} (T_{w} - T_{\infty})}, \\ S h_{x} = \frac{q_{m}}{D_{B} (C_{w} - C_{\infty})}, \\ τ_{m} = - ρ_{f} (1 + \frac{K}{2}) {(\frac{\partial N}{\partial y})}_{y = 0}, \\ q_{w} = - k_{f} {(\frac{\partial T}{\partial y})}_{y = 0}, q_{m} = - D_{B} {(\frac{\partial C}{\partial y})}_{y = 0} \\ τ_{w} = {([(μ + k) \frac{\partial u}{\partial y} + \propto_{1} (2 \frac{\partial u}{\partial x} \frac{\partial u}{\partial y} + u \frac{\partial^{2} u}{\partial x \partial y} + v \frac{\partial^{2} u}{\partial y^{2}})] + KN)}_{y \to 0}, \end{matrix}

(13)

After dimensionalized the above equations, we have

\begin{matrix} {Re}_{x}^{1 / 2} C_{f} = (1 + (1 - n) K) F ″ (0) + β_{1} (7 F ″ (0) F ″ (0) - F (0) F ″' (0)) \\ {Re}_{x}^{\frac{1}{2}} C_{m} = (1 + K / 2) H' (0) \\ {Re}_{x}^{- \frac{1}{2}} N u_{x} = - θ^{'} (0), \\ {Re}_{x}^{- \frac{1}{2}} S h_{x} = - ϕ^{'} (0) . \end{matrix}

(14)

Numerical procedure

Equations (8)–(12) are nonlinear and hence cannot be solved exactly. Therefore we seek a numerical solution by using bvp4c through Matlab. The initial value problem subject to the above dimensionless system is as follows:

\begin{matrix} y (1) = F (η); y (2) = F^{'} (η); y (3) = F^{″} (η); \\ y (4) = F^{″'} (η); yy 1 = F^{″ ″} (η); \\ θ (η) = y (5); θ' (η) = y (6); θ ″ (η) = yy 2; \\ ϕ (η) = y (7); ϕ' (η) = y (8); ϕ ″ (η) = yy 3; \\ H (η) = y (9); H' (η) = y (10); H ″ (η) = yy 4; \end{matrix}

\begin{matrix} yy 1 = (\frac{1}{β_{1} y (1)}) ((1 + K) y (4) + β_{1} (2 xy (3) y (4) \\ + 5 y (2) y (4) + 3 y (3) y (3)) \\ + Ky (10) - 2 y (2) y (2) + y (1) y (3)), \end{matrix}

(15)

\begin{matrix} yy 2 = - \Pr (N_{b} y (6) y (6) + N_{t} y (6) y (8) - 2 y (5) y (2) \\ + y (1) y (6) + Ec (1 + K) y (3) y (3)), \end{matrix}

(16)

yy 3 = - (\frac{N_{t}}{N_{b}} yy 2 - Sc (2 y (2) y (7) - y (1) y (8))),

(17)

\begin{matrix} yy 4 = {(1 + \frac{K}{2})}^{- 1} (y (1) y (10) - 3 y (2) y (9) \\ - K (2 y (9) + y (3))) . \end{matrix}

(18)

Dimensionless boundary conditions takes the form

\begin{matrix} y 0 (1) - S; y 0 (5) - 1 - λ y 0 (6); \\ y 0 (2) - γ; N_{b} y 0 (8) + N_{t} y 0 (6); \\ y 0 (9) + ny 0 (3); \\ yinf (2) - 1; yinf (5); yinf (7); yinf (9) . \end{matrix}

(19)

It is important to note that the following values are used of default parameters

\begin{matrix} β_{1} = 0.1, K = 0.1, \Pr = 1.5, N_{b} = 0.1, N_{t} = 0.1, S = 0.1, \\ Ec = 0.1, Sc = 0.2, n = 0.5, λ = 0.1, and γ = 0.1 . \end{matrix}

Results and discussion

The dimensionless system of nonlinear differential equations are solved through numerical technique. The physical effects of involving parameters under the flow assumptions are discussed in the form of tables and figures. The impact of flow variables on $F^{″} (0)$ , $H' (0)$ , $ϕ' (0)$ , and $θ' (0)$ are highlighted in Tables 1 and 2 for both cases of strong and weak concentration. The variations of the couple stress and skin friction for different values of $S$ for both cases of weak and strong concentrations are presented in Table 1. It is seen that higher values of $S$ reduced the values of skin friction and the couple stress in both cases of strong and weak concentration. The value of skin friction is significant for the case of strong concentration as compared to weak concentration. The impact of stretching parameter on couple stress and skin friction is offered in Table 1 for both cases of strong and weak concentration. The increasing stretching parameter decline the skin friction and enhances the couple stress. The impacts of the material parameter on couple stress and skin friction is also reveals in Table 1. The couple stress and skin friction enhances due to increase in material parameter. The impacts of the micropolar parameter are reveals in Table 1. Both skin friction and couple stress admits significant variations toward micropolar parameter. The impact of variations of the suction parameter on Sherwood and Nusselt numbers are reported in Table 2. The increasing suction parameter improved the Nusselt number and Sherwood number in both cases for strong and weak concentration. The impact of increasing values of Brownian motion parameter is also offered in Table 2. The increasing values of Brownian motion parameter cause decline in the Nusselt number and Sherwood number in both cases for weak and strong concentrations. The influence of thermophoresis parameter on the Sherwood and Nusselt numbers is presented in Table 2. The increasing values of thermophoresis parameter improved the Sherwood and Nusselt numbers in both cases for strong and weak concentration. The influence of stretching parameter on both numbers is also offered in Table 2. Significant impact of the stretching parameter on the Nusselt number and Sherwood number is noticed for both cases namely strong and weak concentration. The impact of variation of material parameter on both numbers is given in Table 2. In absolute sense, the material parameter cause decline in Nusselt number and Sherwood number in case of weak concentration but Nusselt number and Sherwood number shows opposite trends for the case of strong concentration. The impact of thermal slip parameter on Nusselt and Sherwood number is also examined.

Table 1.

Numerical results of skin friction coefficient and couple stress for different values of involving physical parameters in both cases strong (n = 0.0) and weak (n = 0.5) concentrations.

				$n = 0.0$		$n = 0.5$
$S$	$γ$	$β_{1}$	$K$	$F ″ (0)$	$H' (0)$	$F ″ (0)$	$H' (0)$
0.1	0.2	0.1	0.2	1.8241	−0.0571	1.8194	1.3828
0.2	−	−	−	1.8204	−0.0590	1.8157	1.4230
0.3	−	−	−	1.8166	−0.0608	1.8119	1.4639
0.2	0.1	−	−	2.0192	−0.0474	2.0140	1.5175
−	0.2	−	−	1.8241	−0.0571	1.8194	1.3828
−	0.3	−	−	1.6287	−0.0668	1.8194	1.3828
−	0.2	0.1	−	1.8241	−0.0571	1.8194	1.3828
−	−	0.2	−	1.8297	−0.0581	1.8211	1.3841
−	−	0.3	−	1.8345	−0.0592	1.8225	1.3852
−	−	0.1	0.1	1.8429	−0.0875	1.8403	1.3006
−	−	−	0.2	1.8241	−0.0571	1.8194	1.3828
−	−	−	0.3	1.8079	−0.0305	1.8015	1.4550

Table 2.

Numerical results of Sherwood Number and Nusselt Number for different values of involving physical parameters in both cases of strong (n = 0.0) and weak (n = 0.5) concentrations.

										$n = 0.0$		$n = 0.5$
$S$	$N_{b}$	$N_{t}$	$γ$	$β_{1}$	$Sc$	$λ$	$Ec$	$K$	$\Pr$	$θ' (0)$	$φ' (0)$	$θ' (0)$	$φ' (0)$
0.1	0.3	0.2	0.3	0.1	0.2	0.1	0.1	0.2	6.0	−1.3491	−0.4497	−1.2924	−0.4308
0.2	–	–	–	–	–	–	–	–	–	−1.4591	−0.4864	−1.4052	−0.4684
0.3	–	–	–	–	–	–	–	–	–	−1.5651	−0.5217	−1.5144	−0.5048
–	0.1	–	–	–	–	–	–	–	–	−1.8359	−0.6120	−1.7615	−0.5872
–	0.2	–	–	–	–	–	–	–	–	−1.5567	−0.5189	−1.4922	−0.4974
–	0.3	–	–	–	–	–	–	–	–	−1.3491	−0.4497	−1.2924	−0.4308
–	–	0.1	–	–	–	–	–	–	–	−1.3055	−0.4352	−1.2468	−0.4156
–	–	0.2	–	–	–	–	–	–	–	−1.3491	−0.4497	−1.2924	−0.4308
–	–	0.3	–	–	–	–	–	–	–	−1.3912	−0.4637	−1.3383	−0.4461
–	–	–	0.1	–	–	–	–	–	–	−1.4025	−1.4025	−1.2937	−1.2937
–	–	–	0.2	–	–	–	–	–	–	−1.4013	−0.7007	−1.2930	−0.6465
–	–	–	0.3	–	–	–	–	–	–	−1.4002	−0.4667	−1.2930	−0.6465
–	–	–	–	0.1	–	–	–	–	–	−1.3491	−0.4497	−1.2924	−0.4308
–	–	–	–	0.2	–	–	–	–	–	−1.3492	−0.4498	−1.2923	−0.4307
–	–	–	–	0.3	–	–	–	–	–	−1.3493	−0.4499	−1.2922	−0.4306
–	–	–	–	–	–	–	–	–	–	−1.3509	−0.4503	−1.2924	−0.4308
–	–	–	–	–	–	–	–	–	–	−1.3491	−0.4497	−1.2903	−0.4301
–	–	–	–	–	–	–	–	–	–	−1.3472	−0.4491	−1.2881	−0.4294
–	–	–	–	–	–	0.1	–	–	–	−1.3509	−0.4503	−1.3489	−0.4496
–	–	–	–	–	–	0.2	–	–	–	−1.3216	−0.4405	−1.3196	−0.4399
–	–	–	–	–	–	0.3	–	–	–	−1.2923	−0.4308	−1.2903	−0.4301
–	–	–	–	–	–	–	0.1	–	–	−1.3509	−0.4503	−1.3489	−0.4496
–	–	–	–	–	–	–	0.2	–	–	−1.3216	−0.4405	−1.3196	−0.4399
–	–	–	–	–	–	–	0.3	–	–	−1.2923	−0.4308	−1.2903	−0.4301
–	–	–	–	–	–	–	–	0.1	–	−1.2924	−0.4308	−1.3497	−0.4499
–	–	–	–	–	–	–	–	0.2	–	−1.2923	−0.4308	−1.3489	−0.4496
–	–	–	–	–	–	–	–	0.3	–	−1.2922	−0.4308	−1.3482	−0.4494
–	–	–	–	–	–	–	–	–	6.0	−1.2923	−0.4308	−1.3489	−0.4496
–	–	–	–	–	–	–	–	–	6.5	−1.3158	−0.4386	−1.3753	−0.4584
–	–	–	–	–	–	–	–	–	7.0	−1.3377	−0.4459	−1.4000	−0.4667

The influence of Schmidt number on the both numbers is offered in Table 2. The increasing values of Schmidt number cause decline in Nusselt number and Sherwood number for both cases for strong and weak concentration. The influence of Eckert number on both numbers is offered in Table 2. Increasing values of Eckert number cause decline in the Nusselt number and Sherwood number for both cases namely strong and weak concentration. The influence of micropolar parameter on the Nusselt number and Sherwood number is presented in Table 2. Trifling variations are noticed in this case. The influence of Prandtl number on the both numbers is presented in Table 2. The increasing values of Prandtl number enhance the Sherwood and Nusselt numbers for both cases namely strong and weak concentration.

The impacts of the involving parameter on velocity, temperature, concentration, and micropolar profile are highlighted in the Figures 2 to 12. Figure 2 shows the impacts of K on the velocity. It is noted that velocity decreases with increase in micropolar parameter. Physically, vertex viscosity enhances which declines the velocity profile. The influence of material parameter on the velocity profile is given in Figure 3. Figure 3 indicates that velocity of the fluid declines by increasing the values of the material parameter $β_{1}$ of the fluid. The second grade fluid parameter enhances which enhance the viscosity of fluid and as a result the fluid velocity decline. Figure 4 is presented the effects of stretching parameter on the velocity function. The velocity function declines due to enhancing the values of stretching parameter. The variation of the suction parameter with velocity function which reveals in Figure 5. The velocity function found to be decline behavior when the values of suction parameter enhances. Figure 6 indicates the relation between the temperature profile and Prandtl number. The temperature function reduced when values of Prandtl number increase. Physically the thermal diffusivity decreases as the Prandtl number increases. Figure 7 presented the impact of variation of Eckert number on temperature function. The increases in Eckert number cause decrease in the temperature of the fluid. Eckert number indicates the kinetic energy of the flow relative to the boundary layer enthalpy differences. Eckert number shows an important role in high speed flows. Figure 8 shows influence of thermophoresis on temperature function. The temperature function curves increasing due to higher values of thermophoresis parameter. Figure 9 illustrates the effects of thermal slip parameter on the temperature function. The increasing thermal slip parameter cause decline in the temperature function. Figure 10 shows that influence of Brownian motion parameter on temperature function. The temperature function enhances due to higher values of Brownian motion parameter because the fluid temperature enhances due to the increase in the kinetic energy. Figure 11 reveals micropolar profile enhancement due to higher values of micro-gyration. Figure 12 shows the variation of the micropolar profile toward micropolar parameter. We found that micropolar profile is increasing function of micropolar parameter.

Figure 2.

Effect of K on the $F' (η)$ .

Figure 3.

$Effect of β_{1}$ on the $F' (η)$ .

Figure 4.

Effect of $γ$ on the $F' (η)$ .

Figure 5.

Effect of $S$ on the $F' (η)$ .

Figure 6.

Effect of $\Pr$ on the $θ (η)$ .

Figure 7.

Effect of $Ec$ on the $θ (η)$ .

Figure 8.

Effect of $N_{t}$ on the $θ (η)$ .

Figure 9.

Effect of $λ$ on the $θ (η)$ .

Figure 10.

Effect of $N_{b}$ on the $θ (η)$ .

Figure 11.

Effect of $n$ on the $θ (η)$ .

Figure 12.

Effect of $K$ on the $H (η)$ .

Conclusions

Non-Newtonian micropolar fluid flow along with suspended nanoparticles subject to an exponentially stretching surface under the stagnation point flow is considered in the current analysis. The key outcomes are as follows:

The values of skin friction are higher in magnitude for the case of weak concentration (n = 0) in comparison with strong concentration (n = 0.5).

The magnitude of couple stress is observed higher for the case of strong concentration (n = 0.5) as compared to weak concentration (n = 0).

Significant variations are noticed for Nusselt and Sherwood numbers toward involved flow parameters subject to strong and weak cases.

Fluid velocity is found as increasing function of $γ$ while opposite is the case for $S, K and β_{1} .$

Fluid temperature admits direct relation toward Ec, $N_{t}, and Nb$ but inverse is the case for $\Pr and λ .$

Micropolar profile shows direct relation toward $n and K .$

Footnotes

Appendix

Acknowledgements

Authors would like to thank Prince Sultan University for their support through the TAS research lab.

Handling Editor: Chenhui Liang

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iD

Nadeem Abbas

References

Rajagopal

Gupta

AS.

Flow of a viscoelastic fluid over a stretching sheet. Rheologica Acta 1984; 23: 213–215.

Bujurke

Biradar

Hiremath

PS.

Second-order fluid flow past a stretching sheet with heat transfer. Zeitschrift für angewandte Mathematik und Physik 1987; 38: 653–657.

Pontrelli

Flow of a fluid of second grade over a stretching sheet. Int J Non Linear Mech 1995; 30: 287–293.

Sadeghy

Sharifi

Local similarity solution for the flow of a “second-grade” viscoelastic fluid above a moving plate. Int J Non Linear Mech 2004; 39: 1265–1273.

Hayat

Khan

Homotopy solutions for a generalized second-grade fluid past a porous plate. Nonlinear Dyn 2005; 42: 395–405.

Fetecau

Starting solutions for some unsteady unidirectional flows of a second grade fluid. Int J Eng Sci 2005; 43: 781–789.

Khan

Culham

Makinde

OD.

Combined heat and mass transfer of third-grade nanofluids over a convectively-heated stretching permeable surface. Can J Chem Eng 2015; 93: 1880–1888.

Mahapatra

Sidui

Non-axisymmetric Homann stagnation-point flow of a viscoelastic fluid towards a fixed plate. Eur J Mech B/Fluids 2020; 79: 38–43.

Shamshuddin

Salawu

Ogunseye

, et al. Dissipative power-law fluid flow using spectral quasi linearization method over an exponentially stretchable surface with hall current and power-law slip velocity. Int Commun Heat Mass Transf 2020; 119: 104933.

10.

Shamshuddin

Khan

Anwar Bég

, et al. Hall current, viscous and Joule heating effects on steady radiative 2-D magneto-power-law polymer dynamics from an exponentially stretching sheet with power-law slip velocity: a numerical study. Therm Sci Eng Prog 2020; 20: 100732.

11.

Ferdows

Shamshuddin

Zaimi

Dissipative-radiative micropolar fluid transport in a non-Darcy porous medium with cross-diffusion effects. CFD Lett 2020; 12: 70–89.

12.

Dawar

Shah

Kumam

, et al. Chemically reactive MHD micropolar nanofluid flow with velocity slips and variable heat source/sink. Sci Rep 2020; 10: 20926.

13.

Mishra

Shamshuddin

Beg

, et al. Adomain computation of radiative-convective bi-directional stretching flow of a magnetic non-Newtonian fluid in porous media with homogeneous–heterogeneous reactions. Int J Mod Phys B 2020; 34: 2050165.

14.

Shah

Kumam

Selim

, et al. Impact of nanoparticles shape and radiation on the behavior of nanofluid under the Lorentz forces. Case Stud Therm Eng 2021; 26: 101161.

15.

Ferdows

Shamshuddin

Salawu

, et al. Numerical simulation for the steady nanofluid boundary layer flow over a moving plate with suction and heat generation. SN Appl Sci 2021; 3: 264.

16.

Rahman

Roşca

Pop

Boundary layer flow of a nanofluid past a permeable exponentially shrinking/stretching surface with second order slip using buongiorno’s model. Int J Heat Mass Transf 2014; 77: 1133–1143.

17.

Abdul Hakeem

Vishnu Ganesh

Ganga

Magnetic field effect on second order slip flow of nanofluid over a stretching/shrinking sheet with thermal radiation effect. J Magn Magn Mater 2015; 381: 243–257.

18.

Mabood

Das

Melting heat transfer on hydromagnetic flow of a nanofluid over a stretching sheet with radiation and second-order slip. Eur Phys J Plus 2016; 131: 1–12.

19.

Kadir

Mishra

Shamshuddin

, et al. Numerical simulation of von Karman swirling bioconvection nanofluid flow from a deformable rotating disk. In: ICHTFM 2018: 20th international conference on heat transfer and fluid mechanics, Istanbul, Turkey, 16–17 August 2018. http://usir.salford.ac.uk/id/eprint/47891

20.

Bég

Espinoza

DES

Kadir

, et al. Experimental study of improved rheology and lubricity of drilling fluids enhanced with nano-particles. Appl Nanosci 2018; 8: 1069–1090.

21.

Shamshuddin

Beg

Reddy

, et al. Rotating unsteady multi-physico-chemical magneto-micropolar transport in porous media: Galerkin finite element study. Comput Therm Sci 2018; 10: 167–197.

22.

Kumar

Sokhal

GS.

Numerical analysis of performance of water-based nanofluid flowing through tube bent at 90°. Heat Transf Asian Res 2020; 49: 18–32.

23.

Khan

Shehzad

Rauf

, et al. Thermally developed unsteady viscoelastic micropolar nanofluid with modified heat/mass fluxes: a generalized model. Physica A 2020; 550: 123986.

24.

Nadeem

Abbas

Elmasry

, et al. Numerical analysis of water based CNTs flow of micropolar fluid through rotating frame. Comput Methods Programs Biomed 2020; 186: 105194.

25.

Abbas

Malik

Nadeem

Transportation of magnetized micropolar hybrid nanomaterial fluid flow over a Riga curface surface. Comput Methods Programs Biomed 2020; 185: 105136.

26.

Islam

Khan

Kumam

, et al. Radiative mixed convection flow of maxwell nanofluid over a stretching cylinder with joule heating and heat source/sink effects. Sci Rep 2020; 10: 1–18.

27.

Shutaywi

Shah

Mathematical Modeling and numerical simulation for nanofluid flow with entropy optimization. Case Stud Therm Eng 2021; 26: 101198.

28.

Shamshuddin

Mabood

A numerical model for analysis of binary chemical reaction and activation energy of thermo solutal micropolar nanofluid flow through permeable stretching sheet: nanoparticle study. Phys Scr 2021; 96: 075206.

29.

Mabood

Shamshuddin

Mishra

SR.

Characteristics of thermophoresis and Brownian motion on radiative reactive micropolar fluid flow towards continuously moving flat plate: HAM solution. Math Comput Simul 2022; 191: 187–202.

30.

Rehman

Farooq

Ahmed

, et al. Boundary layer stagnation-point flow of second grade fluid over an exponentially stretching sheet. Am J Appl Math Stat 2015; 3: 211–219.

31.

Dero

Rohni

Saaban

MHD micropolar nanofluid flow over an exponentially stretching/shrinking surface: triple solutions. J Adv Res Fluid Mech Therm Sci 2019; 56: 165–174.