Impacts of various shaped Cu-nanomaterial on the heat transfer over a bilateral stretchable surface: Numerical investigation

Abstract

The heat transport in the nanofluids attained much interest of the researchers and engineers due to broad uses in medical sciences, paint industries, aerodynamics, wheel alignment and manufacturing of aircraft parts. Therefore, keeping in mind the paramount significance of the heat transfer, the study of Cu-nanomaterials based nanofluid is conducted. The governing nanofluid model transformed in dimensionless version via similarity transformations. For numerical simulation of the dynamics of Cu-H₂O, RK technique with shooting algorithm is employed and presented behavior of the fluid motion, temperature, wall shear stresses and local thermal performance rate via graphical aid. It is noted that the heat transfer augmented promptly by increasing $λ$ and volumetric fraction of Cu nanomaterial. Further, graphical and tabular comparison is also provided under certain assumptions which authenticate the study.

Keywords

Nanofluids Cu nanoparticle stretching sheet numerical scheme Hamilton and Crossers model

Introduction

The fundamental fluid mechanics problem in two dimensions due to a stretching slit, firstly explored by Crane.¹ Applications of such type of problems comprised in paper production, plastic films drawing, and polymer extrusion etc. Due to versatile uses, researcher’s community focused to analyze such type of problems under certain conditions. These are mixed convection, mass transfer, convective condition, variable slit temperature, suction and injection etc. The flow characteristics of Newtonian fluid in three dimensions past a bi-laterally stretchable sheet in three dimensions reported by Wang² in 1988. In 2014, Sheikholeslami and Ganji³ reported the variations in the nanofluid characteristics in rotating Cartesian coordinates system. For mathematical investigation of the model, they implemented numerical technique and examined the variations of flow quantities in the flow field. Implementation of two-phase model for nanofluid is also part of his work. In 2010, Rashidi and Keimanesh⁴ inspected magnetohydrodynamic flow of laminar liquid film which placed horizontally. The sheet is capable to stretching. They adopted analytical technique for mathematical analysis. The effects of Lorentz force on the nanofluid flow explored in Sheikholeslami et al.⁵ MHD flow of $A l_{2} O_{3}$ -H₂O examined in Sheikholeslami et al.⁶ They focused on the flow characteristics due to fluctuating various flow quantities.

A numerical investigation of thermal transport in the boundary layer flow of the nanofluid over a stretchable surface was conducted by Rana and Bhargava.⁷ They used the model comprising the thermophoretic and Brownian motion effects. The results for the flow regimes against the parameters were drawn by implementing variational finite element technique. In 2015, Das⁸ presented the analysis of nanofluid dynamics over a permeable surface. To find fascinating results of the study, they incorporated slip effects over the surface and tackled the resultant model numerically. The study was interesting in the context of the topic but, the more novel results could be achieved by ingrained the effects of thermal radiations, Lorentz forces and resistive heating.

The dynamics of non-Newtonian nanofluid under convective flow condition is a paramount research topic in the area of Fluid Dynamics. Keeping in view the significance of such flow, Sandeep et al.⁹ analyzed the flow characteristics over a stretching surface with suction/blowing properties of the surface. They observed that thermal performance in Oldroyd-B nanofluid is very effective compared to other fluid like Maxwell and Jeffery nanofluids. In order to make the results more effective, they plugged the magnetic field phenomena in the constitutive model. The heat transfer inspection in radiative and dissipative nanofluid over a stretchable geometry by taking the effects of Lorentz forces into account were found in Pal et al.¹⁰ The influences of suction and injection, Lorentz forces and Grashof number in the dynamics of fluid were core outcomes of the study. Another significant heat transfer investigation in the nanofluid over an inclined stretchable surface is conducted in Thumma et al.¹¹ The Keller’s box technique is implemented for mathematical analysis and displayed the results against the flow quantities. They validated the study by comparing the results with existing literature.

The effects of Lorentz forces on the dynamics of fluid are significant and broadly applicable in the accomplishment of various industrial products. Due to these applications, a heat transfer investigation in the nanofluid by taking the impacts of Lorentz forces into account over a stretching surface were examined in Ibrahim et al.¹² They found that the imposed magnetic field resists the fluid motion due to which maximum heat transfer at the surface is observed. Further, they validate the study by making a comparative analysis. The investigation of second law analysis in Jeffery nanofluid under the impacts of nonlinear thermal radiation and activation energy is presented in Hayat et al.¹³ They observed reverse variations in the temperature and concentration regimes against Deborah number. Further, the escalations in the entropy generation are inspected for stronger Lorentz forces and Bejan number. A 3D heat transport mechanism in Carreau fluid under thermal radiation and heat generation/absorption was determined in Khan et al.¹⁴ They imposed convective condition on the surface and explored the interesting results for the heat transport.

The nonlinear flow model in the presence of various shaped nanoparticles of copper between nonparallel walls was reported in Khan et al.¹⁵ They also discussed the model analytically and explored the velocity and thermal fields graphically. Recently, heat transfer in magneto-nanofluid with natural convection discussed in Sheikholeslami et al.¹⁶ They used Lattice Boltzmann technique for solution purpose. The effects of Ag-nanoparticles on non-Newtonian in the existence of Lorentz forces past a wavering surface studied in Mohd Zin et al.¹⁷ They considered the appropriate geometry of the vertically placed plate.

In 2015, Gul et al.¹⁸ examined the effects of various shaped nanoparticles in porous media. Ahmed et al.¹⁹ examined H₂O flow composed by SWCNT bounded between two Riga plates. They discussed the velocity and thermal profiles graphically. The flow of ferro-fluids over a plate positioned vertically described in Gul et al.²⁰ They also examined the stimulus of convective boundary condition. The characteristics of rotating nanofluid by encountering the influences of magnetic field and radiative heat flux presented in Sheikholeslami et al.²¹ Furthermore, they focused on numerical computation of the model. The effects of variable magnetic field and Cu-H₂O cavity examined in Sheikholeslami and Vajravelu²² and Sheikholeslami,²³ respectively. The study comprised in Mohyud-Din et al.²⁴ confined to H₂O flow composed by Cu nanoparticles in stretchable opening/narrowing channel. Further, they examined the impacts of Lorentz force on the flow properties. The flow of magneto-casson fluid past a stretchable sheet by considering of slip parameter discussed in Ullah et al.²⁵ For numerical and analytical investigation of Newtonian model between two oblique walls described in Adnan et al.^26,27 respectively. Rashidi et al.²⁸ discussed rotating nanofluid model between porous disks. Further, 3D rotating nanofluid squeezed flow discussed in Khan et al.²⁹ For the flow of nanofluids in various geometries under certain assumptions and boundary conditions, we can study^30–36 and references comprised therein.

From the inspection of literature, it is pointed that the heat transfer in the nanofluids under various shape effects over a bidirectional stretchable sheet is not discussed so far. Therefore, the analysis is conducted to fill this research opening. The nanofluid model is formulated via similarity transforms and then handled numerically. At the end, the variations in the velocity and temperature of the nanofluids against various plugged parameters are explored and discussed.

Formulation of the model

Considered the nanofluid flow over a plane stretchable sheet which is capable to stretching in opposite directions. The sheet is placed at $z = 0 .$ The flow over the sheet is due to stretching of the sheet. The velocity in the directions of x and y axes are described by U_w = a/(x + y)⁻ⁿ and V_w = b/(x + y)⁻ⁿ, respectively. Here, invariables are a, b, and n. The variable temperature at the sheet is described by $T_{w} (x, y) = T_{\infty} + A {(x + y)}^{n}$ , in which $A$ is constant. Figure 1 demonstrated the flow configuration:

Figure 1.

Flow configuration of the nanofluid.

In the light of above imposed restriction, the specific set of PDE’s representing the flow of nanofluid flow over a bi-directionally stretchable surface are³⁷:

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

(1)

\tilde{u} \frac{\partial \tilde{u}}{\partial x} + \tilde{v} \frac{\partial \tilde{u}}{\partial y} + \tilde{w} \frac{\partial \tilde{u}}{\partial z} = \frac{μ_{nf}^{*}}{ρ_{nf}^{*}} (\frac{\partial^{2} \tilde{u}}{\partial z^{2}}),

(2)

\tilde{u} \frac{\partial \tilde{v}}{\partial x} + \tilde{v} \frac{\partial \tilde{v}}{\partial y} + \tilde{w} \frac{\partial \tilde{v}}{\partial z} = \frac{μ_{nf}^{*}}{ρ_{nf}^{*}} (\frac{\partial^{2} \tilde{v}}{\partial z^{2}}),

(3)

\tilde{u} \frac{\partial \tilde{T}}{\partial x} + \tilde{v} \frac{\partial \tilde{T}}{\partial y} + \tilde{w} \frac{\partial \tilde{T}}{\partial z} = \frac{k_{nf}^{*}}{{(ρ C_{p})}^{*}_{nf}} (\frac{\partial^{2} \tilde{T}}{\partial z^{2}}),

(4)

Equation (1) presenting the continuity equation in three dimensions which gratify identically. The set of equations (2)–(4) described momentum and energy equations, respectively. Furthermore, thermal conductivity is $k_{nf}^{*}$ , dynamic viscosity and density of the nanofluid denoted by $μ_{nf}^{*}, and ρ_{nf}^{*}$ , respectively. The following theoretical models¹⁸ are used for thermo-physical characteristics of nanofluid:

\begin{matrix} μ_{nf}^{*} = μ_{f}^{*} (1 + a^{*} ϕ^{*} + b^{*} {ϕ^{*}}^{2}), \\ k_{nf}^{*} = k_{f}^{*} [\frac{k_{s}^{*} + (n - 1) k_{f}^{*} + (n - 1) (k_{s}^{*} - k_{f}^{*}) ϕ^{*}}{k_{s}^{*} + (n - 1) k_{f}^{*} - (k_{s}^{*} - k_{f}^{*}) ϕ^{*}}], \\ ρ_{nf}^{*} = \frac{ρ_{f}^{*}}{{(1 - ϕ^{*})}^{- 1}} + ϕ^{*} ρ_{s}^{*}, \\ {(ρ C_{p})}^{*}_{nf} = [(1 - ϕ^{*}) + \frac{ϕ^{*} {(ρ C_{p})}^{*}_{s}}{{(ρ C_{p})}^{*}_{f}}] {(ρ C_{p})}^{*}_{f}, \end{matrix}

In the above models, $a^{*} and b^{*}$ are constants. The empirical shape factor parameter is denoted by $n$ , $ϕ^{*}$ denotes the volumetric fraction of the nanofluid, densities of carrier fluid and nanofluid are $ρ_{f}^{*}$ and $ρ_{s}^{*}$ , respectively. The empirical shapes factor of the nanoparticles, Sphericity, thermal and physical properties of the nanofluid comprised in.¹⁸

The dimensional boundary conditions for the specific model are as under³⁷:

\begin{matrix} At z = 0 : \\ \begin{matrix} \tilde{u} = U_{w} = \frac{a}{{(x + y)}^{- n}} \\ \tilde{v} = V_{w} = \frac{b}{{(x + y)}^{- n}} \\ w = 0 \\ \tilde{T} = T_{w} = T_{\infty} + \frac{A}{{(x + y)}^{- n}} \end{matrix}}, \\ At z \to \infty : \\ \begin{matrix} \tilde{u} \to 0 \\ \tilde{v} \to 0 \\ \tilde{T} \to T_{\infty} \end{matrix}}, \end{matrix}

The similarity variables are³⁷:

\begin{matrix} \tilde{u} = \frac{F^{'} (η) a}{{(x + y)}^{- n}} \\ \tilde{v} = \frac{G^{' (η)} a}{{(x + y)}^{- 1}} \\ \tilde{w} = - \sqrt{a ν_{f}} [\frac{(\frac{n + 1}{2}) (F (η) + G (η))}{{(x + y)}^{\frac{1 - n}{2}}} + \frac{(\frac{n - 1}{2}) η (F^{'} (η) + G^{'} (η))}{{(x + y)}^{\frac{1 - n}{2}}}] \\ β (η) = \frac{T - T_{\infty}}{T_{w} - T_{\infty}} \\ η = \sqrt{\frac{a}{ν_{f}}} z {(x + y)}^{\frac{n - 1}{2}} \end{matrix}},

By entreating the suitable partial derivatives and above defined self-similar transformations in equations (1)–(4), the following is attained:

\begin{matrix} F ″' - \frac{(1 - ϕ^{*}) + ϕ^{*} [\frac{ρ_{s}^{*}}{ρ_{f}^{*}}]}{(1 + a^{*} ϕ^{*} + b^{*} {ϕ^{*}}^{2})} \\ (n F^{'} (F^{'} + G^{'}) + 0.5 (n + 1) F ″ (F + G)) 0, \end{matrix}

(5)

\begin{matrix} G ″' + \frac{(1 - ϕ^{*}) + ϕ^{*} [\frac{ρ_{s}^{*}}{ρ_{f}^{*}}]}{(1 + a^{*} ϕ^{*} + b^{*} {ϕ^{*}}^{2})} \\ (0.5 (n + 1) G ″ (F + G) - nG' (F^{'} + G')) = 0, \end{matrix}

(6)

\begin{matrix} β ″ + \Pr \frac{(1 - ϕ^{*}) + ϕ^{*} [\frac{{(ρ C_{p})}^{*}_{s}}{{(ρ C_{p})}^{*}_{f}}]}{[\frac{k_{s}^{*} + (n - 1) k_{f}^{*} + (n - 1) (k_{s}^{*} - k_{f}^{*}) ϕ^{*}}{k_{s}^{*} + (n - 1) k_{f}^{*} - (k_{s}^{*} - k_{f}^{*}) ϕ^{*}}]} \\ (0.5 (n + 1) β^{'} (F + G) - n β (F^{'} + G^{'})) = 0, \end{matrix}

(7)

The supporting self-similar boundary conditions are:

\begin{matrix} \begin{matrix} F (η) = G (η) = 0 \\ F^{'} (η) = 1 \\ G^{'} (η) = λ \\ β (η) = 1 \end{matrix}} at η = 0, \\ \begin{matrix} F^{'} (η) \to 0 \\ β (η) \to 0 \\ G' (η) \to 0 \end{matrix}} at η \to \infty, \end{matrix}

Here, $\Pr = \frac{μ_{f}^{*} {(c_{p})}^{*}_{f}}{k_{f}^{*}}$ .

Further, dimensional formulas for shear stresses and local Nusselt number are as follows:

C_{Fx} = \frac{τ_{zx}}{ρ_{nf} U_{w}^{2}}, C_{Fy} = \frac{τ_{zy}}{ρ_{nf} V_{w}^{2}} and N u_{x} = \frac{q_{w} (x + y)}{k (T_{w} - T_{\infty})},

Where, $τ_{zx}, τ_{zy}$ and radiative heat flux $q_{w}$ described as under:

\begin{matrix} τ_{zx} = μ_{nf} (\frac{\partial u}{\partial z} + \frac{\partial w}{\partial x}) ↓_{z = 0}, \\ τ_{zy} = μ_{nf} (\frac{\partial v}{\partial z} + \frac{\partial w}{\partial y}) ↓_{z = 0} and q_{w} = - k_{nf} (\frac{\partial T}{\partial z}) ↓_{z = 0}, \end{matrix}

After simplification, transformed version is given as:

\begin{matrix} C_{Fx} \sqrt{R e_{x}} = \frac{{\overset{⌣}{A}}_{2}}{{\overset{⌣}{A}}_{1}} F ″ (η) ↓_{η = 0}, \\ C_{Fy} \sqrt{R e_{y}} λ^{\frac{3}{2}} = \frac{{\overset{⌣}{A}}_{2}}{{\overset{⌣}{A}}_{1}} G ″ (η) ↓_{η = 0} and \\ N u_{x} {(R e_{x})}^{- \frac{1}{2}} = - \overset{β}{{\overset{⌣}{A}}_{4}}' (η) ↓_{η = 0} \\ {\overset{⌣}{A}}_{1} = \frac{1}{{(1 - ϕ^{*})}^{- 1}} + ϕ^{*} [\frac{ρ_{s}^{*}}{ρ_{f}^{*}}], {\overset{⌣}{A}}_{2} = 1 + a^{*} ϕ^{*} + b^{*} {ϕ^{*}}^{2}, \end{matrix}

$\overset{=}{{\overset{⌣}{A}}_{4}} [\frac{k_{s}^{*} + {gk}_{f}^{*} + g (k_{s}^{*} - k_{f}^{*}) ϕ^{*}}{k_{s}^{*} + {gk}_{f}^{*} - (k_{s}^{*} - k_{f}^{*}) ϕ^{*}}]$ , where, $g = n - 1$ .

Numerical analysis of the model

It is a paramount fact that the coupled nonlinear mathematical models fail to tackle in the form of closed solution. Therefore, numerical solutions are credible for such nature of the flow models. The formulated model is highly nonlinear and coupled in nature so, closed form solution does not hold. Thus, the model is then treated numerically with RK technique³⁷ with coupling of shooting technique. In this scheme, firstly the nonlinear mathematical model transformed into a system of first order initial value problem and then solved using the algorithm. Most importantly, $η_{\infty}$ is chosen in such a way that the velocity and thermal fields satisfy the asymptotic behavior far from the surface. In order to obtain the desired first order model, the following transformation are made for the velocity and thermal equations:

\begin{matrix} {\overset{⌣}{s}}_{1} = F, {\overset{⌣}{s}}_{2} = F^{'}, {\overset{⌣}{s}}_{3} = F ″, {\overset{⌣}{s}}_{4} = G, {\overset{⌣}{s}}_{5} = G^{'}, {\overset{⌣}{s}}_{6} = G ″, \\ {\overset{⌣}{s}}_{7} = β and {\overset{⌣}{s}}_{8} = β' . \end{matrix}

In the light of these substitutions, the nonlinear coupled model transformed into the following version:

[\begin{matrix} \overset{⌣}{s}'_{1} \\ \overset{⌣}{s}'_{2} \\ \overset{⌣}{s}'_{3} \\ \overset{⌣}{s}'_{4} \\ \overset{⌣}{s}'_{5} \\ \overset{⌣}{s}'_{6} \\ \overset{⌣}{s}'_{7} \\ \overset{⌣}{s}'_{8} \end{matrix}] = [\begin{matrix} {\overset{⌣}{s}}_{2} \\ {\overset{⌣}{s}}_{3} \\ \frac{(1 - ϕ^{*}) + ϕ^{*} [\frac{ρ_{s}^{*}}{ρ_{f}^{*}}]}{(1 + a^{*} ϕ^{*} + b^{*} {ϕ^{*}}^{2})} (\frac{n + 1}{2} ({\overset{⌣}{s}}_{1} + {\overset{⌣}{s}}_{4}) {\overset{⌣}{s}}_{3} + n ({\overset{⌣}{s}}_{2} + {\overset{⌣}{s}}_{5}) {\overset{⌣}{s}}_{2}) \\ s_{5} \\ s_{6} \\ - \frac{(1 - ϕ^{*}) + ϕ^{*} [\frac{ρ_{s}^{*}}{ρ_{f}^{*}}]}{(1 + a^{*} ϕ^{*} + b^{*} {ϕ^{*}}^{2})} (\frac{n + 1}{2} ({\overset{⌣}{s}}_{1} + {\overset{⌣}{s}}_{4}) {\overset{⌣}{s}}_{6} + n ({\overset{⌣}{s}}_{2} + {\overset{⌣}{s}}_{5}) {\overset{⌣}{s}}_{5}) \\ s_{8} \\ - \frac{(1 - ϕ^{*}) + ϕ^{*} [\frac{{(ρ C_{p})}^{*}_{s}}{{(ρ C_{p})}^{*}_{f}}]}{[\frac{k_{s}^{*} + (n - 1) k_{f}^{*} + (n - 1) (k_{s}^{*} - k_{f}^{*}) ϕ^{*}}{k_{s}^{*} + (n - 1) k_{f}^{*} - (k_{s}^{*} - k_{f}^{*}) ϕ^{*}}]} \Pr (\frac{n + 1}{2} ({\overset{⌣}{s}}_{1} + {\overset{⌣}{s}}_{4}) {\overset{⌣}{s}}_{8} - n s_{2} ({\overset{⌣}{s}}_{2} + {\overset{⌣}{s}}_{5}) {\overset{⌣}{s}}_{7} \end{matrix}],

The appropriate initial conditions are:

[\begin{matrix} {\overset{⌣}{s}}_{1} \\ {\overset{⌣}{s}}_{2} \\ {\overset{⌣}{s}}_{3} \\ {\overset{⌣}{s}}_{4} \\ {\overset{⌣}{s}}_{5} \\ {\overset{⌣}{s}}_{6} \\ {\overset{⌣}{s}}_{7} \\ {\overset{⌣}{s}}_{8} \end{matrix}] = [\begin{matrix} 0 \\ 1 \\ \overset{⌣}{g_{1}} \\ 0 \\ λ \\ \overset{⌣}{g_{2}} \\ 1 \\ \overset{⌣}{g_{3}} \end{matrix}]

Where, $\overset{⌣}{g_{1}} = F ″ (0), \overset{⌣}{g_{2}} = G ″ (0)$ and $\overset{⌣}{g_{3}} = β' (0)$ . After this, integration is carried out and Mathematica 10.0 is used for further computation of the results.

Results and discussion

This section deals with the graphical results and their comprehensive discussion. The results are presented against various flow parameters over the domain of interest. these results are displayed in Figures 2 to 10. The results are displayed for the velocities, temperature, shear stresses and local thermal performance of the nanofluids against the pertinent flow parameters.

Figure 2.

Effects of $λ$ on $F' (η)$ .

Figure 3.

Effects of $λ$ on $G' (η)$ .

Figure 4.

Effects of $λ$ on $β (η)$ .

Figure 5.

Effects of $λ$ on the shear stresses for (a) $n = 1$ and (b) $n = 2$ in $x$ -direction.

Figure 6.

Effects of $ϕ$ on the shear stresses for (a) $n = 1$ and (b) $n = 2$ in $x$ -direction.

Figure 7.

Effects of $λ$ on the shear stresses (a) $n = 1$ and (b) $n = 2$ in $y$ -direction.

Figure 8.

Effects of $ϕ$ on the shear stresses (a) $n = 1$ and (b) $n = 2$ in $y$ -direction.

Figure 9.

Effects of $λ$ on the heat transport rate for (a) $n = 1$ and (b) $n = 2$ .

Figure 10.

Effects of $ϕ$ on the heat transport rate for (a) $n = 1$ and (b) $n = 2$ .

Velocity and thermal distribution

The trends in horizontal velocity component $F' (η)$ against arising values of the parameter $λ$ are displayed in Figure 2. The velocity distribution is sketched for platelets, blades, cylinders, and bricks shaped based nanofluid. It is examined that the fluid velocity in the horizontal direction drops for higher values of $λ$ . The motion of the nanofluid comprising platelets nanomaterial is slow in comparison with rest of the nanofluids. The physical reason behind these trends is the thermophysical characteristics of the nanomaterial. The effective density of platelets nanomaterial is less compared to rest of the nanomaterials. Due to low dense material, the internal molecular forces become weaker which allows the fluid to move freely. On the other side, rest of the nanofluids become more denser due to high density of bricks, blades and cylinders nanomaterial. Ultimately internal fluid forces become dominant which opposes the fluid motion.

The vertical fluid motion against $λ$ for various nanofluids is examined in Figure 3. It is observed that the vertical velocity $G' (η)$ augmented for increasing values of $λ$ . The velocity rises slowly for the nanofluids composed by blades and bricks nanomaterial. The physical reason behind these trends is the effective density. Due to higher density of these materials, the resultant fluid becomes more denser which leads to slow increment in the vertical component of the velocity. Near the stretchable surface abrupt fluid motion is examined and its gradually slow down away from the it. The velocity of the nanofluids asymptotically vanishes at ambient level of the surface.

Thermal behavior of the nanofluids for the multiple ranges of the parameter $λ$ are portrayed in Figure 4. It is obvious that $λ$ opposes the fluid temperature under the particular scenario. The temperature promptly declines in the region $0.2 < η < 0.8$ . Beyond the region, fluid temperature asymptotically vanishes and obeys the ambient temperature condition. The thermal boundary layer region strated beyond $η > 1.2$ . A slow decrement in the temperature is observed for the nanofluid composed by blades nanomaterial in which thermophysical characteristics playing dominant role especially thermal conductivity.

Trends in shear stresses

The trends in stresses along horizontal and vertical directions due to varying flow parameters are displayed in Figures 5 to 8, respectively. The stresses along horizontal direction against $ϕ$ versus $λ$ are depicted in Figure 5. It is examined that the shear stresses drop inconsequently as the volumetric fraction tends to zero and abruptly drops for higher volumetric fraction. The physical reason for negligible decrement in the shear stresses for smaller volumetric fraction is that the fluid under goes to its conventional liquid state for zero $ϕ$ due to which the decrement is negligible. The more decreasing trends in the shear stresses are examined for stronger volumetric fraction of the nanomaterials. These trends are elucidated in Figure 5(a) and (b), respectively.

Figure 6 portraying the trends in shear stresses against $ϕ$ versus $λ$ . The volumetric fraction of the nanomaterials varies curve wise while $λ$ varies horizontally. The shear stresses promptly decline for higher $λ$ . These trends controlled by increasing the index parameter $n$ in the flow model. In both the cases, abrupt decrement is examined for the nanofluid composed by platelets shaped nanomaterial in which the effective density playing its significant role.

The behavior of vertical shear stresses by increasing the values of $λ$ and volumetirc fraction of the nanomaterial is presented in Figures 7 and 8, respectively. From Figure 7, it is observed that the vertical stresses drops very rapidly by increasing the parameter $λ$ and this decrement becomes almost similar as the volumetric fraction tends to zero. For the nanofluid containing bricks shaped nanomaterial, the vertical trends in the shear stresses becomes slow to rest of the nanofluids. Further, it is obvious that by increasing index parameter value, the vertical stresses declines abruptly. In Figure 8, the higher volumetric fraction causes the decrement in the shear stresses but it is slow compared to the trends displayed in Figure 7. For higher $λ$ and $ϕ$ , prompt decrement is examined for vertical shear stresses in Figure 8.

Trends in local thermal performance

Figures 9 and 10 elucidating the behavior of local thermal performance in the nanofluids under varying $λ$ and the volumetric fraction $ϕ$ , respectively. Form Figures 9 and 10, it is examined that the local thermal performance rate escalated for higher $λ$ and $ϕ$ . However, slow increment is observed as the volume fraction and $λ$ tends to zero. The maximum rate of thermal transport is investigated for the nanofluid suspended by blades shaped nanomaterial.

Authentication of the study

The under consideration nanofluid model is modified version of Khan et al.³⁷ In order to verify the presented results with existing science literature under the same boundary conditions and the same flow parameters, we need to impose some restriction ( $ϕ = 0$ ) on our model mentioned in Figure 11(a) and (b), respectively. From Figure 11(a) and (b), it is evident that the presented study has an excellent agreement with the previous literature Khan et al.³⁷ which validate the study.

Figure 11.

Comparison of present study with previous work Khan et al.³⁷ for (a) F′(η) and (b) β(η)

Moreover, a comparative analysis for $F ″ (0)$ and $G ″ (0)$ is provided for zero volumetric fraction of the nanomaterial. The results are computed and given in Table 1. From the table, the present results are in very good agreement with existing science literature by taking $η_{\infty} = 15$ . This proves the validation of the analysis.

Table 1.

Authentication of the study for $F ″ (0)$ and $G ″ (0)$ , $ϕ = 0$ .

$n$	$λ$	Present		Khan et al.³⁷
		$F ″ (0)$	$G ″ (0)$	$F ″ (0)$	$G ″ (0)$
$1.0$	$0.0$	$- 1.0$	$0$	$- 1.0$	$0$
	$0.5$	$- 1.22474$	$- 0.612372$	$- 1.224745$	$- 0.612372$
	$1.0$	$- 1.41421$	$- 1.41421$	$- 1.414214$	$- 1.414214$

Conclusions

The thermal transport mechanism is examined under various shaped nanomaterial based nanofluids. The colloidal model is treated numerically and displayed the results against the flow quantities over the feasible region. It is observed that the vertical motion of the nanofluids significantly escalated by increasing $λ$ . The temperature of the nanofluids containing bricks, blades, cylinders and platelets shaped nanomaterial drop for higher $λ$ . Further, the shear stresses along horizontal direction declines for higher volumetric fraction and it could be controlled by increasing the value of index parameter $n$ . The local thermal performance of the nanofluids improved against higher $ϕ$ and significant contribution of blades nanomaterial is examined regarding thermal rate analysis. Finally, the comparative investigation is provided which showing good agreement with existing literature.

Footnotes

Appendix

Handling Editor: Chenhui Liang

Author Note

Syed Tauseef Mohyud-Din is also affiliated to University of Multan (UoM), Multan, Pakistan.

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 iDs

Adnan

Umar Khan

Syed Tauseef Mohyud-Din

References

Crane

LJ.

Flow past a stretching plate. Z Angew Math Phys 1970; 21: 645–647.

Wang

CY.

Stretching a surface in a rotating fluid. Z Angew Math Phys 1988; 39: 177–185.

Sheikholeslami

Ganji

DD.

Numerical investigation for two phase modeling of nanofluid in a rotating system with permeable sheet. J Mol Liq 2014; 194: 13–19.

Rashidi

Keimanesh

Using differential transform method and padé approximant for solving MHD flow in a laminar liquid film from a horizontal stretching surface. Math Probl Eng 2010; 2010: 491319.

Sheikholeslami

Ganji

Gorji-Bandpy

, et al. Magnetic field effect on nanofluid flow and heat transfer using KKL model. J Taiwan Inst Chem Eng 2014; 45: 795–807.

Sheikholeslami

Gorji-Bandpy

Ganji

, et al. Magnetohydrodynamic free convection of Al₂O₃-water nanofluid considering thermophoresis and Brownian motion effects. Comput Fluids 2014; 94: 147–160.

Rana

Bhargava

Flow and heat transfer of a nanofluid over a nonlinearly stretching sheet: a numerical study. Commun Nonlinear Sci Numer Simul 2012; 17: 212–226.

Das

Nanofluid flow over a non-linear permeable stretching sheet with partial slip. J Egypt Math Soc 2015; 23: 451–456.

Sandeep

Rushi Kumar

Jagadeesh Kumar

MS.

A comparative study of convective heat and mass transfer in non-Newtonian nanofluid flow past a permeable stretching sheet. J Mol Liq 2015; 212: 585–591.

10.

Pal

Mandal

Vajravelu

MHD convection–dissipation heat transfer over a non-linear stretching and shrinking sheets in nanofluids with thermal radiation. Int J Heat Mass Transf 2013; 65: 481–490.

11.

Thumma

Anwar Bég

Kadir

Numerical study of heat source/sink effects on dissipative magnetic nanofluid flow from a non-linear inclined stretching/shrinking sheet. J Mol Liq 2017; 232: 159–173.

12.

Ibrahim

Shankar

Nandeppanavar

MM.

MHD stagnation point flow and heat transfer due to nanofluid towards a stretching sheet. Int J Heat Mass Transf 2013; 56: 1–9.

13.

Hayat

Kanwal

Qayyum

, et al. Entropy generation optimization of MHD Jeffrey nanofluid past a stretchable sheet with activation energy and non-linear thermal radiation. Phys A Stat Mech Appl 2020; 544: 123437.

14.

Khan

Sardar

Hashim . Heat generation/absorption and thermal radiation impacts on three-dimensional flow of Carreau fluid with convective heat transfer. J Mol Liq 2018; 272: 474–480.

15.

Khan

Ahmed

Mohy-ud-Din

ST.

Analysis of magnetohydrodynamic flow and heat transfer of Cu–water nanofluid between parallel plates for different shapes of nanoparticles. Neural Comput Appl 2018; 29: 695–703.

16.

Sheikholeslami

Gorji-Bandpy

Ganji

DD.

Lattice Boltzmann method for MHD natural convection heat transfer using nanofluid. Powder Technol 2014; 254: 82–93.

17.

Mohd Zin

Khan

Shafie

The impact silver nanoparticles on MHD free convection flow of Jeffrey fluid over an oscillating vertical plate embedded in a porous medium. J Mol Liq 2016; 222: 138–150.

18.

Gul Khan

Shafie

. Energy transfer in mixed convection MHD flow of nanofluid containing different shapes of nanoparticles in a channel filled with saturated porous medium. Nanoscale Res Lett 2015; 10: 490.

19.

Ahmed

Adnan Khan

, et al. Influence of thermal radiation and viscous dissipation on squeezed flow of water between Riga plates saturated with carbon nanotubes. Colloids Surf A Physicochem Eng Asp 2017; 522: 389–398.

20.

Gul

Khan

Shafie

, et al. Heat transfer in MHD mixed convection flow of a ferrofluid along a vertical channel. PLoS One 2015; 10: e0141213.

21.

Sheikholeslami

Hatami

Ganji

DD.

Nanofluid flow and heat transfer in a rotating system in the presence of a magnetic field. J Mol Liq 2014; 190: 112–120.

22.

Sheikholeslami

Vajravelu

Nanofluid flow and heat transfer in a cavity with variable magnetic field. Appl Math Comput 2017; 298: 272–282.

23.

Sheikholeslami

CuO-water nanofluid free convection in a porous cavity considering Darcy law. Eur Phys J Plus Epub ahead of Print 30 January 2017. DOI: 10.1140/epjp/i2017-11330-3.

24.

Mohyud-Din

Khan

Ahmed

, et al. Magnetohydrodynamic flow and heat transfer of nanofluids in stretchable convergent/divergent channels. Appl Sci 2015; 5: 1639–1664.

25.

Ullah

Shafie

Khan

Effects of slip condition and Newtonian heating on MHD flow of casson fluid over a nonlinearly stretching sheet saturated in a porous medium. J King Saud Univ Sci 2017; 29: 250–259.

26.

Adnan Khan

Ahmed

, et al. Thermo-diffusion and diffusion-thermo effects on flow of second grade fluid between two inclined plane walls. J Mol Liq 2016; 224: 1074–1082.

27.

Adnan Asadullah

Khan

, et al. Analytical and numerical investigation of thermal radiation effects on flow of viscous incompressible fluid with stretchable convergent/divergent channels. J Mol Liq 2016; 224: 768–775.

28.

Rashidi

Abelman

Freidooni Mehr

Entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid. Int J Heat Mass Transf 2013; 62: 515–525.

29.

Khan

Ahmed

Mohy-ud-Din

ST.

Numerical investigation for three dimensional squeezing flow of nanofluid in a rotating channel with lower stretching wall suspended by carbon nanotubes. Appl Therm Eng 2017; 113: 1107–1117.

30.

Ullah

Khan

Shafie

MHD natural convection flow of Casson nanofluid over nonlinearly stretching sheet through porous medium with chemical reaction and thermal radiation. Nanoscale Res Lett 2016; 11: 527.

31.

Ali

Gohar

Khan

MHD flow of water-based Brinkman type nanofluid over a vertical plate embedded in a porous medium with variable surface velocity, temperature and concentration. J Mol Liq 2016; 223: 412–419.

32.

Sheikholeslami

CVFEM for magnetic nanofluid convective heat transfer in a porous curved enclosure. Eur Phys J Plus 2016; 131: 413.

33.

Sheikholeslami

Hayat

Alsaedi

Numerical study for external magnetic source influence on water based nanofluid convective heat transfer. Int J Heat Mass Transf 2017; 106: 745–755.

34.

Aman

Khan

Ismail

, et al. Impacts of gold nanoparticles on MHD mixed convection Poiseuille flow of nanofluid passing through a porous medium in the presence of thermal radiation, thermal diffusion and chemical reaction. Neural Comput Appl 2018; 30: 789–797.

35.

Choi

Enhancing thermal conductivity of fluids with nanoparticles in developments and applications of non-Newtonians flows. ASME 1995; 66: 99–105.

36.

Sheikholeslami

Shehzad

SA.

Magnetohydrodynamic nanofluid convection in a porous enclosure considering heat flux boundary condition. Int J Heat Mass Transf 2017; 106: 1261–1269.

37.

Khan

Mustafa

Hayat

, et al. On three-dimensional flow and heat transfer over a non-linearly stretching sheet: analytical and numerical solutions. PLoS One 2014; 9: e107287.