Unsteady squeezing flow of magnetohydrodynamic carbon nanotube nanofluid in rotating channels with entropy generation and viscous dissipation

Abstract

The most favorable gift of modern science is nanofluid. The nanofluid can able to move freely through micro channels with the spreading of nanoparticles. Due to improved convection between the base liquid surfaces and nanoparticles, the nano suspensions express high thermal conductivity. Also, the benefits of suspending nanoparticles in base fluids are increased heat capacity, surface area, effective thermal conductivity, collision, and interaction among particles. This research aim to study squeezing flow of carbon nanotubes based on water (H₂O) in rotating channels. Entropy generation is evaluated and for this purpose, second law of thermodynamics is employed. The influences of thermal radiation, viscous dissipation, and applied magnetic field on nanofluid are taken into account. The flow of the nanofluid is considered in unsteady three dimensions. The transformed ordinary differential equations (ODEs) are solved by homotopy analysis method with the help of similarity variables. Results obtained for single and multi-wall carbon nanotubes are compared. Plots have been presented in order to examine how the velocities, temperature, and entropy profiles become affected by numerous physical parameters. Generally, the velocity profiles escalate when the upper plate of the channel moves toward the lower stretching one and reduces when the upper plate is moving away from the lower one. The velocity profile in y-direction escalates with the escalation in nanoparticle volume fraction and suction parameter while the rotation parameter bids dual behavior with the escalating values. The velocity profile in x-direction bids the oscillatory behavior with the enhancement in nanoparticle volume fraction, rotation parameter, and magnetic parameter. The physical properties of carbon nanotubes, thermo physical properties of carbon nanotubes and nanofluid of some base fluids, and thermal conductivity of carbon nanotubes with different volume fractions are shown through tables.

Keywords

Entropy generation squeezing flow single-walled carbon nanotubes and multi-walled carbon nanotubes magnetohydrodynamic rotating channels thermal radiation homotopy analysis method

Introduction

Presently, the goal of many investigators is to examine the nanotechnology which was presented by Choi¹ in 1995. He brings into being higher thermal conductivity of nanofluids than base fluids such as ethylene glycol and water. Nanofluid is one of the key parts of the nanotechnology which is potential heat transfer fluid. Carbon nanotube (CNT) nanofluids have many application in structural materials (textiles, body arms, concrete, polyethylene, sports equipment, bridges, flywheels, fire protections), electromagnetic (buckypaper, light bulb filament, magnets, solar cells, electromagnetic antenna), electroacoustic (loudspeaker), chemicals (air pollution filter, water filter, chemical nanowires, sensors), mechanicals (oscillator, waterproof), optical, electrical circuits, interconnects, transistors, cancer treatment, cardiac autonomic regulation, in drug delivery, platelet activation, tissues regeneration, super capacitors, thin film electronics, actuators, and so on. Due to these applications the researchers are busy to work on CNTs. Kang et al.² revealed the precise results of Choi by testing experimentally. The preparation impacts of CNT nanofluid stability has been inspected by Ganji et al.³ Haq et al.⁴ reported that single-walled carbon nanotubes (SWCNTs) have higher Nusselt number and skin friction than multi-walled carbon nanotubes (MWCNTs) by considering water as a base fluid. Liu et al.⁵ deliberated by studying ethylene glycol and synthetic engine oil in the existence of MWCNTs that ethylene glycol with CNTs have higher thermal conductivities than ethylene glycol suspension without CNTs. The theory of nanofluid paste over an exponentially stretching sheet has been presented by Nadeem and Lee.⁶ With Reynolds’ model and Vogel’s model, Ellahi et al.⁷ analytically probed the series solution of non-Newtonian nanofluids. Nadeem et al.⁸ numerically studied the magnetohydrodynamic boundary layer flow of Maxwell fluid over a stretching sheet in the presence of nanoparticles. The other related studies to boundary layer flow can be seen in Nadeem and Haq.^9,10 Sheikholeslami¹¹ numerically examined the nanofluid free convection in a porous enclosure under the impact of electric field. Shah et al.¹² investigated the micropolar nanofluid flow of Casson fluid between two rotating parallel plates with hall effect. Nasir et al.¹³ investigated the three-dimensional rotating flow of magnetohydrodynamic SWCNTs over a stretching sheet under the influence thermal radiation. Dawar et al.¹⁴ examined the magnetohydrodynamic CNTs Casson nanofluid and radiative heat transfer in rotating channels. Sheikholeslami¹⁵ analyzed the nanofluid flow magnetic field in a porous media considering magnetic field and Brownian motion. The other related studies of Sheikholeslami can be seen in the literature.^16–26

Entropy is a disorder of a system and surrounding. It occurs when heat transmission occurs, because some additional movements happen when it moves, for example, molecular vibration and friction, displacements of molecules, kinetic energy, and spin movements, due to which loss of useful heat occurs and thus heat cannot be transmitted fully into work. Chaos in a system and surrounding created due to these additional movements. This microscopic chaos results in macroscopic level occurs because of some unnecessary irreversibilities. For example, electric resistance, friction, mixing of fluids, unstained expansion, chemical reaction, inelastic deformation of solids, and unnecessary heat transmission in finite temperature difference. The entropy generation was originally formulated by Bejan²⁷ over an unsteady stretching surface. Sarojamma et al.²⁸ has examined the entropy generation on a thin film flow. Along with an inclined permeable surface, Soomro et al.²⁹ recently examined numerically the entropy generation in magnetohydrodynamic (MHD) water-based CNTs. Mansour and Sahin³⁰ examined the entropy generation rate in a laminar viscous flow in a circular flow and have deliberated that the entropy generation rate is high near the wall than that of the center of the pipe. Rashidi at al.³¹ examined entropy generation over a permeable stretching sheet. The related study about entropy generation can be seen in the literature.^32–35 The flow of nanoparticles in a rotating system through entropy generation has been examined by Hayat et al.³⁶ The flow of nanofluids with spherical heat source/sink through entropy generation has been examined by Nouri et al.³⁷ Over a stretching surface, the flow of Jeffery nanofluids through entropy generation has been examined by Dalir et al.³⁸ The unsteady squeezing flow of viscous fluid through entropy generation has been examined by Ahmed et al.³⁹ The related studies about entropy generation can be seen in Rashidi et al.^40–42 Afridi et al.⁴³ investigated the second law analysis of boundary layer flow with energy dissipation and variable fluid properties over an exponentially stretching sheet through entropy generation. They found that entropy generation can be minimizing by escalating viscosity parameter and reducing the operating temperature. Some important study of nanofluid with respect to their modern application can be seen in the literature.^44–51

The focal intention of this inspection is to examine the three-dimensional squeezing flow of CNTs nanofluids considering MHD and thermal radiation impacts through entropy generation. The influences of emerging parameters are deliberated through graphs. The formulated problem has been solved by homotopy analysis method (HAM).^52–58 The impacts of emerging parameters on velocities profiles, temperature profile, entropy generation, Bejan number, skin fraction coefficient, and local Nusselt number have been shown through graphs and discussed in detail as well. The formulation of the problem is presented in section “Problem statement.” In section “Entropy analysis and Bejan number,” the formulation of entropy generation and Bejan number is presented. In section “Solution by HAM,” the results and discussion of the graphs is presented. The analytical solution of the problem is presented in section “Results and discussion.” The central points of this investigation are presented in section “Conclusion.”

Problem statement

Consider three-dimensional squeezing flow of CNTs nanofluids (SWCNT, MWCNT-nanoparticles) based on water $(H_{2} O)$ in a rotating channel with fixed and porous lower wall. The magnetic field $B_{0}$ is taken along y-direction and the nanofluid flow taken under the influence of thermal radiation. The lower wall of the channel is stretched along $y = 0$ with velocity $U_{w} = ax / (1 - α t)$ along x-axis and $v = - (V_{0} / (1 - α t))$ along y-axis. The channel height is varying with stretching out and shriveling from top to bottom of the walls. This height change follows from $v_{h} = dh / dt$ . For the stated problem, all the assumption and conditions are expressed in equation (9). The modeling of this problem is as follows⁵²

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

(1)

\begin{matrix} \frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + \frac{2 β}{1 - α t} w + \frac{1}{ρ_{nf}} \frac{\partial p^{*}}{\partial x} \\ - \frac{μ_{nf}}{ρ_{nf}} (\frac{\partial^{2} u}{\partial x^{2}} + \frac{\partial^{2} u}{\partial y^{2}}) + \frac{σ_{nf} B_{0}^{2}}{ρ_{nf}} \frac{u}{1 - α t} = 0 \end{matrix}

(2)

\frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y} + \frac{1}{ρ_{nf}} \frac{\partial p^{*}}{\partial y} - \frac{μ_{nf}}{ρ_{nf}} (\frac{\partial^{2} v}{\partial x^{2}} + \frac{\partial^{2} u}{\partial y^{2}}) = 0

(3)

\begin{matrix} \frac{\partial w}{\partial t} + u \frac{\partial w}{\partial x} + v \frac{\partial w}{\partial y} - \frac{2 β}{1 - α t} w - \frac{μ_{nf}}{ρ_{nf}} (\frac{\partial^{2} w}{\partial x^{2}} + \frac{\partial^{2} w}{\partial y^{2}}) \\ + \frac{σ_{nf} B_{0}^{2}}{ρ_{nf}} \frac{w}{1 - α t} = 0 \end{matrix}

(4)

\begin{matrix} \frac{\partial T}{\partial t} + u \frac{\partial T}{\partial x} + v \frac{\partial w}{\partial y} - \frac{μ_{nf}}{{(ρ C_{p})}_{nf}} (\frac{\partial^{2} T}{\partial x^{2}} + \frac{\partial^{2} T}{\partial y^{2}}) \\ - \frac{1}{{(ρ C_{p})}_{nf}} \frac{\partial Qr}{\partial y} - \frac{μ_{nf}}{{(ρ C_{p})}_{nf}} \\ {4 {(\frac{\partial u}{\partial x})}^{2} + (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}) + {(\frac{\partial w}{\partial x})}^{2} + {(\frac{\partial w}{\partial y})}^{2}} \\ + \frac{σ_{nf} B_{0}^{2}}{ρ_{nf}} \frac{(u^{2} + w^{2})}{1 - α t} = 0 \end{matrix}

(5)

Here $u, v, and w$ are the components of velocity. Also $p^{*}, σ_{nf}, μ_{nf}, ρ_{nf}, (ρ C_{p})_{nf}, T, Qr$ indicate the modified pressure, electrical conductivity, dynamic viscosity, density, heat capacitance, temperature, and radiative heat flux, respectively (Figure 1).

Figure 1.

Geometrical representation of the problem.

The thermal flux is defined as follows

Qr = - \frac{4 σ^{*} \partial T^{4}}{3 k^{*} \partial y}

(6)

where $σ^{*}$ is Stefan-Boltzman constant and $k^{*}$ is the absorption coefficient. By Taylor series expansion, $T^{4}$ in terms of $T_{\infty}$ is defined as follows

T^{4} \approx 4 {TT}_{\infty}^{3} - 3 T_{\infty}^{3}

(7)

Substituting equations (6) and (7) in equation (5), we obtain

\begin{matrix} \frac{\partial T}{\partial t} + u \frac{\partial T}{\partial x} + v \frac{\partial w}{\partial y} - \frac{μ_{nf}}{{(ρ C_{p})}_{nf}} (\frac{\partial^{2} T}{\partial x^{2}} + \frac{\partial^{2} T}{\partial y^{2}}) \\ - \frac{1}{{(ρ C_{p})}_{nf}} \frac{16 σ^{*}}{3 k^{*}} \frac{\partial}{\partial y} (T^{3} \frac{\partial T}{\partial y}) - \frac{μ_{nf}}{{(ρ C_{p})}_{nf}} \\ {4 {(\frac{\partial u}{\partial x})}^{2} + (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}) + {(\frac{\partial w}{\partial x})}^{2} + {(\frac{\partial w}{\partial y})}^{2}} \\ + \frac{σ_{nf} B_{0}^{2}}{ρ_{nf}} \frac{(u^{2} + w^{2})}{1 - α t} = 0 \end{matrix}

(8)

The boundary conditions for the above stated problem are as follows

\begin{matrix} u = U_{w} = \frac{ax}{1 - α t}, v = - \frac{V_{0}}{1 - α t}, \\ w = 0, T = T_{w} at y = 0 \\ u = 0, v = V_{h} = - \frac{α}{2} {(\frac{υ}{a (1 - α t)})}^{1 / 2}, \\ w = 0, T = T_{h} as y = h (t) \end{matrix}

(9)

The density of the nanofluid is as follows

ρ_{nf} = (1 - φ) ρ_{f} + φ ρ_{CNT}

(10)

The heat capacity of the nanofluid is as follows

{(ρ C_{p})}_{nf} = (1 - φ) {(ρ C_{p})}_{f} + φ {(ρ C_{p})}_{CNT}

(11)

The dynamic viscosity of the nanofluid is as follows

μ_{nf} = \frac{μ_{f}}{{(1 - φ)}^{2.5}}

(12)

The thermal conductivity of the nanofluid is as follows

k_{nf} = \frac{1 - φ + 2 φ \frac{k_{CNT}}{k_{CNT} - k_{f}} \ln \frac{k_{CNT} + k_{f}}{2 k_{f}}}{1 - φ + 2 φ \frac{k_{f}}{k_{CNT} - k_{f}} \ln \frac{k_{CNT} + k_{f}}{2 k_{f}}} k_{f}

(13)

The electrical conductivity of the nanofluid is as follows

σ_{nf} = 1 + \frac{3 (\frac{σ_{CNT}}{σ_{f}} - 1) φ}{(\frac{σ_{CNT}}{σ_{f}} + 2) - (\frac{σ_{CNT}}{σ_{f}} - 1) φ} σ_{f}

(14)

Here, $ρ_{CNT}, ρ_{f}, φ, k_{f}, and k_{CNT}$ are the density of CNTs, the fluid density, volume fraction of nanoparticles, the thermal conductivity of base fluid and CNTs, respectively.

In order to simplify equations (1)–(4), (8), and (9), the following definitions are given

\begin{matrix} u = U_{w} f' (λ), v = {(\frac{a ν}{a (1 - α t)})}^{1 / 2} f (λ), \\ w = U_{w} g (λ), θ (λ) = \frac{T - T_{w}}{T_{w} - T_{h}}, λ = \frac{y}{h (t)} \end{matrix}

(15)

Here $f' (λ), f (λ), g (λ), and θ (λ)$ is the velocity function along x-axis, y-axis, rotational velocity function, and dimensionless temperature function, respectively, and $λ$ is the local similarity variable.

By substituting equation (15) into equations (1)–(4), (8), and (9), the final obtained equations are as follows

\begin{matrix} f^{iv} - {(1 - φ)}^{2.5} Z_{1} \\ [(f' f ″ - ff ″') + 2 Ω g' + \frac{S}{2} (3 f ″ + λ f ″')] - Z_{5} Mf ″ = 0 \end{matrix}

(16)

\begin{matrix} g ″ + {(1 - φ)}^{2.5} Z_{1} \\ [(fg' - gf') + 2 Ω f' - S (g + \frac{λ}{2} g')] - Z_{5} Mg = 0 \end{matrix}

(17)

\begin{matrix} (1 + \frac{R}{Z_{4}}) θ ″ + \frac{Z_{2} \Pr}{Z_{4} {(1 - φ)}^{2.5}} (θ f' - \frac{S}{2} λ θ') \\ + \frac{Z_{2}}{Z_{4}} \Pr {\begin{matrix} Ec (4 f^{2} + g^{2}) \\ + E c_{m} ({f ″}^{2} + {g'}^{2}) \end{matrix}} - Z_{5} M ({f'}^{2} + g^{2}) = 0 \end{matrix}

(18)

With the following boundary conditions

\begin{matrix} f (0) = A, f' (0) = 1, g (0) = 0, θ (0) = 1 \\ f (1) = \frac{S}{2}, f' (1) = 0, g (1) = 0, θ (1) = 0 \end{matrix}

(19)

In equations (16)–(19), the followings parameters and numbers are $Ω = β / a$ rotation parameter, $S = α / a$ squeezing parameter, $R = (16 σ^{*} T^{3}) / (3 k^{*} k_{f})$ radiation parameter, $M = σ_{f} B_{0}^{2} h / ρ_{f} ν_{f}$ magnetic parameter, $A = V_{0} / ah$ suction parameter, $\Pr = ((μ_{f} (C_{p})_{f}) / k)$ Prandtl number, $Ec = v^{2} / (h^{2} (Cp)_{f} (T_{w} - T_{h}))$ Eckert number, and $E c_{m} = U_{w}^{2} / ((Cp)_{f} (T_{w} - T_{h}))$ modified Eckert number. Also, $Z_{1} = ρ_{nf} = (1 - φ) ρ_{f} + φ ρ_{CNT},$ $Z_{2} = (ρ C_{p})_{nf} = (1 - φ) (ρ C_{p})_{f} + φ {(ρ C_{p})}_{CNT},$ $Z_{3} = μ_{nf} = (μ_{f} / (1 - φ)^{2.5}),$ $Z_{4} = (k_{nf} / k_{f}) = (1 - φ + 2 φ (k_{CNT} / (k_{CNT} - k_{f})) \ln ((k_{CNT} + k_{f}) / 2 k_{f})) / (1 - φ + 2 φ (k_{f} / (k_{CNT} - k_{f})) \ln ((k_{CNT} + k_{f}) / 2 k_{f})),$ and $Z_{5} = σ_{nf} = 1 + (3 ((σ_{CNT} / σ_{f}) - 1) φ / ((σ_{CNT} / σ_{f}) + 2) - ((σ_{CNT} / σ_{f}) - 1) φ) σ_{f}$ .

The skin fraction coefficient $(C f_{x})$ and local Nusselt number $(N u_{x})$ are defined as follows

{(Cf)}_{lower} = {(\frac{μ_{nf} (\frac{\partial u}{\partial y})}{ρ_{nf} U_{0}^{2}})}_{y = 0}, {(Cf)}_{upper} = {(\frac{μ_{nf} (\frac{\partial u}{\partial y})}{ρ_{nf} U_{0}^{2}})}_{y = h (t)}

(20)

\begin{matrix} {(Nu)}_{lower} = {(\frac{X (\frac{\partial T}{\partial y})}{k_{nf} (T_{w} - T_{h})})}_{y = 0} + Qr, \\ {(Nu)}_{upper} = {(\frac{X (\frac{\partial T}{\partial y})}{k_{nf} (T_{w} - T_{h})})}_{y = h (t)} + Qr \end{matrix}

(21)

After simplifying equations (20) and (21), $(C f_{x})$ and $(N u_{x})$ are defined as follows

\begin{matrix} {(Cf {Re}_{x}^{1 / 2})}_{lower} = (\frac{1}{Z_{1} {(1 - φ)}^{5 / 2}}) f ″ (0), {(Cf {Re}_{x}^{1 / 2})}_{upper} \\ = (\frac{1}{Z_{1} {(1 - φ)}^{5 / 2}}) f ″ (1) \end{matrix}

(22)

\begin{matrix} {(Nu {Re}_{x}^{- 1 / 2})}_{lower} = - (Z_{4} + R) θ' (0), {(Nu {Re}_{x}^{- 1 / 2})}_{upper} \\ = - (Z_{4} + R) θ' (1) \end{matrix}

(23)

where ${Re}_{x} = (U_{w} X / ν_{f})$ is the local Reynolds number.

Entropy analysis and Bejan number

Entropy is a disorder of a system and surrounding. It occurs when heat transmission occurs, because some additional movements happen when it moves. The loss of useful heat occurs due to kinetic energy, spin movements, and so on and thus heat cannot be transmitted fully into work. Due to these additional movements, Chaos in a system and surrounding created. For the above stated problem, the local entropy generation rate is⁵³

\begin{matrix} E_{G} = \frac{k_{nf}}{T_{w}^{2}} (1 + \frac{16 σ^{*} T^{3}}{3 k^{*} k_{f}}) {{(\frac{\partial T}{\partial x})}^{2} + {(\frac{\partial T}{\partial y})}^{2}} \\ + \frac{μ_{nf}}{T_{w}} {\begin{matrix} 4 {(\frac{\partial u}{\partial x})}^{2} + (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}) \\ + {(\frac{\partial w}{\partial x})}^{2} + {(\frac{\partial w}{\partial y})}^{2} \end{matrix}} + \frac{σ_{nf}}{T_{w}} B_{0}^{2} (u^{2} + w^{2}) \end{matrix}

(24)

Executing equation (15), the above equation becomes

\begin{matrix} Ns = {E_{G}}_{0} [(Z_{4} + R) {θ'}^{2} + \frac{γ \Pr Z_{3}^{2}}{Z_{1} Z_{4}} {Ec (4 {f'}^{2} + g^{2}) + \frac{Z_{1} E c_{m}}{Z_{3}} (f' + {g'}^{2})} + \frac{Z_{5} ME c_{m}}{Z_{3}} ({f'}^{2} + g^{2})] \end{matrix}

(25)

where $E_{G 0}, γ$ is the characteristics entropy generation and dimensionless temperature, respectively, whose expression are given by

{E_{G}}_{0} = \frac{k_{f} {(T_{w} - T_{h})}^{2}}{T_{h}^{2} h^{2}}, g = \frac{T_{h}}{T_{w} - T_{h}}

(26)

The non-dimensional form of equation (25) is as follows

Ns = \frac{E_{G}}{{E_{G}}_{0}} = [\begin{matrix} (Z_{4} + R) {θ'}^{2} + \frac{γ \Pr Z_{3}^{2}}{Z_{1} Z_{4}} {Ec (4 {f'}^{2} + g^{2}) + \frac{Z_{1} E c_{m}}{Z_{3}} (f' + {g'}^{2})} \\ + \frac{Z_{5} ME c_{m}}{Z_{3}} ({f'}^{2} + g^{2}) \end{matrix}]

(27)

The Bejan number Be is defined as follows

Be = \frac{(Z_{4} + R) {θ'}^{2}}{[\begin{matrix} (Z_{4} + R) {θ'}^{2} + \frac{γ \Pr Z_{3}^{2}}{Z_{1} Z_{4}} {Ec (4 {f'}^{2} + g^{2}) + \frac{Z_{1} E c_{m}}{Z_{3}} (f' + {g'}^{2})} \\ + \frac{Z_{5} ME c_{m}}{Z_{3}} ({f'}^{2} + g^{2}) \end{matrix}]}

(28)

It is obvious from equation (28) that the Bejan number takes values from 0 to 1. $Be = 0$ leads to the case of irreversibility, due to magnetic field and fluid fraction dictate over the heat transfer irreversibility, while $Be = 1$ shows domination of irreversibility due to heat transfer on entropy generation. $Be = 1 / 2$ leads to the case where the irreversibility due to heat transfer is same as the sum of the irreversibility due to magnetic fields and fluid fraction.

Solution by HAM

To solve equations (16)–(18) with the boundary condition (19), we use HAM with the succeeding process.

The initial suppositions are chosen as follows

\begin{matrix} f_{0} (λ) = A + λ + (\frac{3}{2} S - 2 - 3 A) λ^{2} \\ + (1 + 2 A - S) λ^{3}, g_{0} (λ) = 0, θ_{0} (λ) = 1 - λ \end{matrix}

(29)

The $L_{f}, L_{g}, and L_{θ}$ are taken as follows

L_{f} (f) = f^{iv}, L_{g} (g) = g ″, L_{θ} (θ) = θ ″

(30)

which have the succeeding properties

\begin{matrix} L_{\bar{F}} (r_{1} + r_{2} λ + r_{3} λ^{2} + r_{4} λ^{3}) = 0, \\ L_{\bar{G}} (r_{5} + r_{6} λ) = 0, L_{\bar{ϕ}} (r_{7} + r_{8} λ) = 0 \end{matrix}

(31)

where the general solution of the problem have $r_{i} (i = 1 - 8)$ constants:

The resultant non-linear operators $N_{f}, N_{g}, and N_{θ}$ are indicated as follows

\begin{matrix} N_{f} [f (λ; τ), g (λ; τ)] = \frac{\partial^{4} f (λ; τ)}{\partial λ^{4}} - {(1 - φ)}^{2.5} Z_{1} \\ [\begin{matrix} {\frac{\partial f (λ; τ)}{\partial λ} \frac{\partial^{2} f (λ; τ)}{\partial λ^{2}} - f (λ; τ) \frac{\partial^{3} f (λ; τ)}{\partial λ^{3}}} \\ + 2 Ω \frac{\partial g (λ; τ)}{\partial λ} - \frac{S}{2} {3 \frac{\partial f^{2} (λ; τ)}{\partial λ^{2}} + λ \frac{\partial^{3} f (λ; τ)}{\partial λ^{3}}} \end{matrix}] - Z_{5} M \frac{\partial f^{2} (λ; τ)}{\partial λ^{2}} = 0 \end{matrix}

(32)

\begin{matrix} N_{g} [f (λ; τ), g (λ; τ)] = \frac{\partial^{2} g (λ; τ)}{\partial λ^{2}} + {(1 - φ)}^{2.5} Z_{1} \\ [\begin{matrix} {f (λ; τ) \frac{\partial g (λ; τ)}{\partial λ} - g (λ; τ) \frac{\partial f (λ; τ)}{\partial λ}} \\ + 2 Ω \frac{\partial f (λ; τ)}{\partial λ} - S {g (λ; τ) + \frac{λ}{2} \frac{\partial g (λ; τ)}{\partial λ}} \end{matrix}] - Z_{5} Mg (λ; τ) = 0 \end{matrix}

(33)

\begin{matrix} N_{θ} [f (λ; τ), g (λ; τ), θ (λ; τ)] \\ = (1 + \frac{R}{Z_{4}}) \frac{\partial^{2} θ (λ; τ)}{\partial λ^{2}} + \frac{Z_{2}}{Z_{4}} \Pr {\begin{matrix} θ (λ; τ) \frac{\partial f (λ; τ)}{\partial λ} \\ - \frac{S}{2} λ \frac{\partial θ (λ; τ)}{\partial λ} \end{matrix}} \\ + {Ec (4 {(f (λ; τ))}^{2} + {(g (λ; τ))}^{2}) + E c_{m} ({(\frac{\partial^{2} f (λ; τ)}{\partial λ^{2}})}^{2} + {(\frac{\partial g (λ; τ)}{\partial λ})}^{2})} \\ - Z_{5} M ({(\frac{\partial f (λ; τ)}{\partial λ})}^{2} + {(g (λ; τ))}^{2}) \end{matrix}

(34)

The 0th-order problems from equations (16)–(18) are as follows

(1 - τ) L_{f} [f (λ; τ) - f_{0} (λ)] = τ ℏ_{f} N_{f} [f (λ; τ), g (λ; τ)]

(35)

(1 - τ) L_{g} [g (λ; τ) - g_{0} (λ)] = τ ℏ_{g} N_{g} [f (λ; τ), g (λ; τ)]

(36)

(1 - τ) L_{θ} [θ (λ; τ) - θ_{0} (λ)] = τ ℏ_{θ} N_{θ} [f (λ; τ), g (λ; τ), θ (λ; τ)]

(37)

The equivalent boundary conditions are as follows

\begin{matrix} {f (λ; τ) |}_{λ = 0} = A, {\frac{\partial f (λ; τ)}{\partial λ} |}_{λ = 0} = 1, \\ {f (λ; τ) |}_{λ = 1} = \frac{S}{2}, {\frac{\partial f (λ; τ)}{\partial λ} |}_{λ = 1} = 0 \\ {g (λ; τ) |}_{λ = 0} = 0, {g (λ; τ) |}_{λ = 1} = 0 \\ {θ (λ; τ) |}_{λ = 0} = 1, {θ (λ; τ) |}_{λ = 1} = 0 \end{matrix}

(38)

When $τ = 0 and τ = 1$ , we have

f (λ; 1) = f (λ), g (λ; 1) = g (λ), and θ (λ; 1) = θ (λ)

(39)

Expanding $f (λ; τ), g (λ; τ), and θ (λ; τ)$ by Taylor’s series

\begin{matrix} f (λ; τ) = f_{0} (λ) + \sum_{q = 1}^{\infty} f_{q} (λ) τ^{q} \\ g (λ; τ) = g_{0} (λ) + \sum_{q = 1}^{\infty} g_{q} (λ) τ^{q} \\ θ (λ; τ) = θ_{0} (λ) + \sum_{q = 1}^{\infty} θ_{q} (λ) τ^{q} \end{matrix}

(40)

where

\begin{matrix} f_{q} (λ) = {\frac{1}{q!} \frac{\partial f (λ; τ)}{\partial λ} |}_{τ = 0}, g_{q} (λ) = {\frac{1}{q!} \frac{\partial g (λ; τ)}{\partial λ} |}_{τ = 0}, \\ and θ_{q} (λ) = {\frac{1}{q!} \frac{\partial θ (λ; τ)}{\partial λ} |}_{τ = 0} \end{matrix}

(41)

The secondary constraints $ℏ_{f}, ℏ_{g}, and ℏ_{θ}$ are selected such that the series (40) converges at $τ = 1$ , changing $τ = 1$ in equation (40), we get

\begin{matrix} f (λ) = f_{0} (λ) + \sum_{q = 1}^{\infty} f_{q} (λ) \\ g (λ) = g_{0} (λ) + \sum_{q = 1}^{\infty} g_{q} (λ) \\ θ (λ) = θ_{0} (λ) + \sum_{q = 1}^{\infty} θ_{q} (λ) \end{matrix}

(42)

The $q th$ -order problem satisfies the following

\begin{matrix} L_{f} [f_{q} (λ) - χ_{q} f_{q - 1} (λ)] = ℏ_{f} U_{q}^{f} (λ), \\ L_{g} [g_{q} (λ) - χ_{q} g_{q - 1} (λ)] = ℏ_{g} U_{q}^{g} (λ), \\ L_{θ} [θ_{q} (λ) - χ_{q} θ_{q - 1} (λ)] = ℏ_{θ} U_{q}^{θ} (λ) . \end{matrix}

(43)

The equivalent boundary conditions are as follows

\begin{matrix} f_{q} (0) = {f'}_{q} (0) = f_{q} (1) = {f'}_{q} (1) = 0 \\ g_{q} (0) = g_{q} (1) = 0 \\ θ_{q} (0) = θ_{q} (1) = 0 \end{matrix}

(44)

Here

\begin{matrix} U_{q}^{f} (λ) = f_{q - 1}^{iv} - {(1 - φ)}^{2.5} Z_{1} \\ [{\sum_{k = 0}^{q - 1} {f'}_{q - 1 - k} {f ″}_{k} - \sum_{k = 0}^{q - 1} f_{q - 1 - k} {f ″'}_{k}} \\ + 2 Ω {g'}_{q - 1} - \frac{S}{2} {3 {f ″}_{q - 1} + λ {f ″'}_{q - 1}}] - Z_{5} {Mf}_{q - 1}^{″} \end{matrix}

(45)

\begin{matrix} U_{q}^{g} (λ) = g_{q - 1}^{″} + {(1 - φ)}^{2.5} Z_{1} \\ [{\sum_{k = 0}^{q - 1} f_{q - 1 - k} {g'}_{k} - \sum_{k = 0}^{q - 1} g_{q - 1 - k} {f'}_{k}} \\ + 2 Ω {f ″}_{q - 1} - S {g_{q - 1} + \frac{λ}{2} {g'}_{q - 1}} \end{matrix}] - Z_{5} M g_{q - 1}

(46)

\begin{matrix} U_{q}^{θ} (λ) = (1 + \frac{R}{Z_{4}}) {θ ″}_{q - 1} + \frac{Z_{2}}{Z_{4}} \Pr {\sum_{k = 0}^{q - 1} θ_{q - 1 - k} {f'}_{k} - \frac{S}{2} λ {θ'}_{q - 1}} \\ + \frac{Z_{2}}{Z_{4}} \Pr {Ec (4 \sum_{k = 0}^{q - 1} f_{q - 1 - k} f_{k} + \sum_{k = 0}^{q - 1} g_{q - 1 - k} g_{k}) + E c_{m} (\sum_{k = 0}^{q - 1} {f ″}_{q - 1 - k} {f ″}_{k} + \sum_{k = 0}^{q - 1} {g'}_{q - 1 - k} {g'}_{k})} \\ - Z_{5} M (\sum_{k = 0}^{q - 1} {f'}_{q - 1 - k} {f'}_{k} + \sum_{k = 0}^{q - 1} g_{q - 1 - k} g_{k}) \end{matrix}

(47)

where

χ_{q} = {\begin{matrix} 0, if τ \leq 1 \\ 1, if τ > 1 \end{matrix}

Results and discussion

In this segment, we have briefly studied the influence of emerging parameters and non-dimensional numbers on linear and angular velocity functions $(f (λ), f' (λ), g (λ))$ , temperature function $(θ (λ))$ , entropy generation $(Ns)$ , Bejan number $(Be)$ , skin fraction coefficient $(Cf {Re}_{x}^{1 / 2})$ and local Nusselt number $(Nu {Re}_{x}^{1 / 2})$ for both SWCNTs and MWCNTs based on water $(H_{2} O)$ . These emerging parameters and non-dimensional numbers are nanoparticle volume fraction $(φ)$ , rotation parameter $(Ω)$ , squeeze parameter $(S)$ , suction parameter $(A)$ , radiation parameter $(R)$ , magnetic parameter $(M)$ , Eckert number $(Ec)$ , and modified Eckert number $(E c_{m})$ , respectively.

Velocity functions $(f (λ), f' (λ),$ and $g (λ))$

Figures 2 –16 are plotted to see the influence of emerging parameters and non-dimensional numbers on linear and angular velocity functions $(f (λ), f' (λ), and g (λ))$ for both SWCNTs and MWCNTs based on water $(H_{2} O)$ . The impacts of nanoparticle volume fraction with increasing amount on $(f (λ), f' (λ), and g (λ))$ are shown in Figures 2 –4. Figure 2 indicates the impact of $(φ)$ on $f (λ)$ . From here, we see that the increasing values of $(φ)$ shows enhancement in velocity function $f (λ)$ in y-direction for both SWCNTs and MWCNTs. This leads us to the claim that $f (λ)$ permanently has a superior value at $(S > 0)$ that of $(S < 0)$ . This is due to the fact that the density of SWCNTs is greater than those of MWCNTs. Figure 3 is plotted to see the impact of $(φ)$ on $f' (λ)$ . In the interval $0 \leq λ \leq 4$ , the increasing rate of nanoparticle in velocity function along x-direction is declining, while in the interval $4 < λ \leq 1$ the variation in nanoparticle volume fraction along x-direction shows escalation. It should be illustrious that when $S > 0$ , the variations in velocity function along x-direction are higher. When $S < 0$ , then this leads to the enhancement in rotational velocity function $g (λ)$ . This result is shown in Figure 4. It is also observed that $g (λ)$ has bigger value for the case of MWCNTs than those of SWCNTs. It should be noteworthy that $g (λ)$ has larger value at $S > 0$ than that of $S < 0$ . The impacts of rotation parameter $(Ω)$ on $(f (λ), f' (λ))$ $and g (λ)$ are shown in Figures 5 –7. It can be understood from Figure 5 that $(Ω)$ gives dual behavior on $f (λ)$ . It can also be seen that the escalation in $(Ω)$ shows reducing behavior in the interval $0 \leq λ \leq 0.6$ for both cases of SWCNTs and MWCNTs, while in the interval $0.6 < λ \leq 1$ shows escalating behavior with the escalating values of rotation parameter. It is significant to mention that values of $f (λ)$ for MWCNTs is greater than those of SWCNTs. Also the velocity at $S > 0$ is higher than that of $S < 0$ . Figure 6 displays the impact of rotation parameter on $f' (λ)$ . From here, we see that $f' (λ)$ expresses an oscillatory behavior comparative to the variation in rotation parameter. At closed distances to the upper and lower wall of the channel, the increase in the rotation parameter, in the intervals $0 \leq λ \leq 0.3$ and $0.8 \leq λ \leq 1.0$ has decreased $f' (λ)$ in both cases of SWCNTs and MWCNTs. But in mean position of the channel $0.3 < λ < 0.8$ , the escalating behavior is observed. Figure 7 shows the impact of rotation parameter on $g (λ)$ . From here, we see that the variation in $(Ω)$ has reverse relationship with $g (λ)$ . We see from here that for the case of MWCNTs, the velocity function $g (λ)$ is always higher than those of SWCNTs, and when $S > 0$ the velocity enhances more quickly compared to when $S < 0$ . The impacts of squeezing parameter $(S)$ on $(f (λ), f' (λ), and g (λ))$ are shown in Figures 8 –10. Figure 8 depicts the impact of S on $f (λ)$ . From here, we see that the increase in the squeeze parameter $(S)$ enhances the velocity function in y-direction. Also, the velocity increases when the channel is contracted $(S > 0)$ and the velocity reduces when the channel is stretched out $(S < 0)$ . Similar impacts can be seen in Figures 9 and 10. In the dominant area of the channel, there is a reverse flow that produced a noteworthy decline in the velocity, while at the upper and lower walls of the channel, this reverse slowly vanishes. Figures 11 –13 are schemed to see the impact of suction parameter $(A)$ on $(f (λ), f' (λ), and g (λ))$ . In both $(S > 0)$ and $(S < 0)$ cases, the increasing suction parameter enhances the velocity in y-direction, which is shown in Figure 11. Due to the fluid suction in the area near to the lower wall of the channel, a significant enhancement in velocity is observed. It should be noted that the velocity is enhanced by MWCNTs more than those of SWCNTs. Figure 12 indicates the relationship between the distribution rate of nanofluid velocity $f' (λ)$ and $(A)$ . From here, we see the reverse relationship between $(A)$ and $f' (λ)$ . The rise in $(A)$ results decline in $f' (λ)$ . It is interesting to mention that the reverse flow is greater when $(S < 0)$ as compared to $(S > 0)$ . In addition $f' (λ)$ in the case of MWCNTs is larger than those of SWCNTs. The same result is plotted for the rotational velocity function $g (λ)$ in Figure 13. With the difference that $g (λ)$ in the case of SWCNTs has been reduced more those of MWCNTs. Figures 14 –16 are plotted to see the impact of magnetic parameter $(M)$ on velocity functions $(f (λ), f' (λ), and g (λ))$ . According to Lorentz force theory, the opposing of the flow is due to the existence of $(M)$ in electrically conducting fluid. This opposing force inclines to reduce the velocity of the flow. Figures 14 and 16 show reducing behavior in velocity functions $(f (λ) and g (λ))$ , while in Figure 15, the velocity function $f' (λ)$ shows reduction behavior in the interval $0.4 \leq λ \leq 1.0$ . This behavior is due to the stretching of the lower plate. In Figure 14, it interesting to mention that the reverse flow is greater when $(S > 0)$ as compared to $(S < 0)$ . This behavior is due to the opposing force of $(M)$ . Also, we see from the figures that for the case of MWCNTs, the velocity functions $(f (λ), f' (λ), and g (λ))$ are always higher than those of SWCNTs.

Figure 2.

Influence of $φ$ on $f (λ)$ .

Figure 3.

Influence of $φ$ on $f' (λ)$ .

Figure 4.

Influence of $φ$ on $g (λ)$ .

Figure 5.

Outcome of $Ω$ on $f (λ)$ .

Figure 6.

Influence of $Ω$ on $f' (λ)$ .

Figure 7.

Influence of $Ω$ on $g (λ)$ .

Figure 8.

Influence of S on $f (λ)$ .

Figure 9.

Influence of S on $f' (λ)$ .

Figure 10.

Influence of S on $g (λ)$ .

Figure 11.

Influence of A on $f (λ)$ .

Figure 12.

Influence of A on $f' (λ)$ .

Figure 13.

Influence A on $g (λ)$ .

Figure 14.

Influence of M on $f (λ)$ .

Figure 15.

Influence of M on $f' (λ)$ .

Figure 16.

Influence of M on $g (λ)$ .

Temperature function $(θ (λ))$

Figures 17 –22 are plotted to see the influence of emerging parameters and non-dimensional numbers on temperature function $(θ (λ))$ for both SWCNTs and MWCNTs based on water $(H_{2} O)$ . Figure 17 shows the relationship between $(φ)$ and $θ (λ)$ . From here, we see that the variation in $θ (λ)$ is reducing with the escalating values of $(φ)$ . The case where the channel is contracted or the case where the channel is stretched out, this result is identical for both cases of SWCNTs and MWCNTs. Figure 18 shows the impact of $(S)$ on $θ (λ)$ . From here, we see that growing $(S)$ in both cases $(S > 0)$ and $(S < 0)$ increases the nanofluid temperature for both SWCNTs and MWCNTs. Figure 19 shows the impact of M on $θ (λ)$ . From here, we see the increasing behavior in $θ (λ)$ with the increasing values of magnetic parameter. Figure 20 indicates the impact of $(R)$ on $θ (λ)$ . The temperature function $θ (λ)$ reduces with escalation in $(R)$ for both cases $(S > 0)$ and $(S < 0)$ . This behavior physically understands as the enhancement in heat transfer rate and the decline in $θ (λ)$ caused by enhancement in $(R)$ . It should be noted that $θ (λ)$ is more declined for the case of SWCNTs than those of MWCNTs when $S > 0$ , but this result is opposite when $S < 0$ . There is a direct relationship between $θ (λ)$ and $(Ec)$ , and $(E c_{m})$ . The influence of $(Ec)$ and $(E c_{m})$ on $θ (λ)$ for both $(S > 0)$ and $(S < 0)$ cases are shown in Figures 21 and 22, respectively. This direct relationship physically represents the relationship between the velocity and $(Ec)$ , and hence the kinetic energy. Escalating Eckert number escalates the intermolecular collision and this escalates the temperature function. In addition, the temperature function for MWCNTs increases more than those of SWCNTs in the case of $(S > 0)$ , it is due to the fact that MWCNTs has less amount of thermal conductivity than the SWCNTs which results the reverse impact obtained in the case of $(S < 0)$ .

Figure 17.

Influence of $φ$ on $θ (λ)$ .

Figure 18.

Influence of S on $θ (λ)$ .

Figure 19.

Influence of M on $θ (λ)$ .

Figure 20.

Influence of R on $θ (λ)$ .

Figure 21.

Influence of Ec on $θ (λ)$ .

Figure 22.

Influence of $E c_{m}$ on $θ (λ)$ .

Entropy analysis (Ns) and Bejan number (Be)

Figures 23 –32 are plotted to see the impact of emerging parameters and non-dimensional numbers on entropy generation rate $(Ns)$ and Bejan number $(Be)$ for both SWCNTs and MWCNTs. Figure 23 is plotted to see the impact of $(R)$ on $(Ns)$ for both cases of SWCNTs and MWCNTs. From Figure 23, we observed that as the radiative parameter $(R)$ becomes higher entropy generation rate across the channel escalates. This is due to the enhancement in emission of $(R)$ . The impact of $(R)$ on $(Be)$ is depicted in Figure 24. The rise in $(R)$ results the enhancement in $(Be)$ for both SWCNTs and MWCNTs. Figure 25 illustrates the influence of $(Ec)$ on $(Ns)$ . The escalating values of $(Ec)$ shows dual behavior in $(Ns)$ . In the region $(0 \leq λ \leq 0.5)$ , it shows increasing behavior, while in the region $(0.5 < λ \leq 1.0)$ , it shows decreasing behavior. Figure 26 depicts the escalating behavior in $(Be)$ with the escalation in $(Ec)$ . Figure 27 depicts the impact of $(M)$ on $(Ns)$ . From here, we see that $(M)$ gives dual behavior with the escalation in $(Ns)$ . While from Figure 28 we have seen that the escalation in $(M)$ gives decreasing behavior in $(Be)$ . Figure 29 illustrates the impact of $(E c_{m})$ on $(Ns)$ . The escalating $(E c_{m})$ gives dual behavior in $(Ns)$ . In the region $(0 \leq λ \leq 0.3)$ , the modified Eckert number gives increasing behavior, while in the region $(0.3 < λ \leq 1.0)$ , the modified Eckert number gives reducing behavior. The Influence of $(E c_{m})$ on $(Be)$ is shown in Figure 30. We have seen the increasing behavior in $(Be)$ with the escalation in $(E c_{m})$ . Figure 31 illustrates the impact of $(\Pr)$ on $(Ns)$ . From here, we have seen decreasing behavior in the region $(0 \leq λ \leq 0.39)$ , while in the region $(0.39 < λ \leq 1.0)$ , we have seen the escalating behavior in $(Ns)$ with the escalation in $(\Pr)$ . The escalating $(\Pr)$ gives decreasing behavior in $(Be)$ , which is shown in Figure 32.

Figure 23.

Influence of R on Ns.

Figure 24.

Influence of R on Be.

Figure 25.

Influence of Ec on Ns.

Figure 26.

Influence of Ec on Be.

Figure 27.

Influence of M on Ns.

Figure 28.

Influence of M on Be.

Figure 29.

Influence of $E c_{m}$ on Ns.

Figure 30.

Influence of $E c_{m}$ on Be.

Figure 31.

Influence of $\Pr$ on Ns.

Figure 32.

Influence of $\Pr$ on Be.

Skin fraction (Cf Re_x^1/2) and local Nusselt number (Nu Re_x^1/2)

Figures 33 –42 are plotted to see the impact of emerging parameters and non-dimensional numbers on skin fraction coefficient $(Cf {Re}_{x}^{1 / 2})$ and local Nusselt number $(Nu {Re}_{x}^{- 1 / 2})$ for both SWCNTs and MWCNTs. Figures 33 and 34 illustrate the influence of $(Ω)$ on $(Cf {Re}_{x}^{1 / 2})$ at both lower and upper walls of the channel. The surface drag force reduces with the escalation in $(Ω)$ at the lower wall of the channel in both $(S > 0)$ and $(S < 0)$ cases, which is shown in Figure 33. In addition, when the upper and lower walls are separated further more from each other, $(Cf {Re}_{x}^{1 / 2})$ reduces abundance more. Also in the case of MWCNTs, the skin fraction coefficient has reduced more than those of SWCNTs. While reverse impact can be seen at the upper wall of the channel, which is shown in Figure 34. Figures 35 –38 depict the impact of squeezing parameter $(S)$ on $(Cf {Re}_{x}^{1 / 2})$ and $(Nu {Re}_{x}^{- 1 / 2})$ at both upper and lower walls of the channel. Generally, escalating squeeze parameter, the $(Cf {Re}_{x}^{1 / 2})$ increase at the lower wall of the channel and reduces at the upper wall of the channel. In addition, at the lower wall of the channel, $(Cf {Re}_{x}^{1 / 2})$ has reduced in the case of SWCNTs than in that of MWCNTs. Also, there is a direct relationship between $(S)$ and $(Nu {Re}_{x}^{- 1 / 2})$ at the lower wall of the channel, while at the upper wall of the channel there is a reverse relationship between $(S)$ and $(Nu {Re}_{x}^{- 1 / 2})$ . Figures 39 and 40 depict the impact of $(M)$ on $(Cf {Re}_{x}^{1 / 2})$ at both upper and lower walls of the channel. It can be seen from both figures that the increase in $(M)$ increases $(Cf {Re}_{x}^{1 / 2})$ and reduces $(Nu {Re}_{x}^{- 1 / 2})$ at both lower and upper walls of the channel in both $(S > 0)$ and $(S < 0)$ cases. In addition, $(Cf {Re}_{x}^{1 / 2})$ has increased in the case of MWCNTs than those of SWCNTs, while $(Nu {Re}_{x}^{- 1 / 2})$ reduced in the case of SWCNTs than in that of MWCNTs. Figures 41 and 42 depict to see the effect of $(R)$ on $(Nu {Re}_{x}^{- 1 / 2})$ at both lower and upper walls of the channel. When the upper wall moves to the lower wall of the channel, the escalating $(R)$ escalates $(Nu {Re}_{x}^{- 1 / 2})$ . It is obvious from the figure that escalating $(Nu {Re}_{x}^{- 1 / 2})$ in the case of SWCNTs is more than those of MWCNTs.

Figure 33.

Influence of $Ω$ on $(Cf {Re}_{x}^{1 / 2})_{lower}$ .

Figure 34.

Influence of $Ω$ on $(Cf {Re}_{x}^{1 / 2})_{upper}$ .

Figure 35.

Influence of S on $(Cf {Re}_{x}^{1 / 2})_{lower}$ .

Figure 36.

Influence of S on $(Cf {Re}_{x}^{1 / 2})_{upper}$ .

Figure 37.

Influence of S on $(Nu {Re}_{x}^{- 1 / 2})_{lower}$ .

Figure 38.

Influence of S on $(Nu {Re}_{x}^{- 1 / 2})_{upper}$ .

Figure 39.

Influence of M on $(Cf {Re}_{x}^{1 / 2})_{lower}$ .

Figure 40.

Influence of M on $(Cf {Re}_{x}^{1 / 2})_{upper}$ .

Figure 41.

Influence of R on $(Nu {Re}_{x}^{- 1 / 2})_{lower}$ .

Figure 42.

Influence of R on $(Nu {Re}_{x}^{- 1 / 2})_{upper}$ .

Tables 1 –3 are schemed to study the Physical properties of CNTs, thermo physical properties CNTs and nanofluids of some base fluids, and thermal conductivity ( $k_{nf}$ ) of CNTs with different volume fraction $(φ)$ , respectively.

Table 1.

Physical properties of CNTs.

Materials	SWCNT	MWCNT
Thermal conductivity, $k_{nf}$ (W/mK)	3000	3000
Electrical conductivity, $φ_{nf}$ (S/m)	$10^{6} - 10^{7}$	$10^{6} - 10^{7}$
Tensile strength (GPa)	150	150
Young’s modulus (GPa)	1054	1200

CNTs: carbon nanotubes; SWCNT: single-walled carbon nanotube; MWCNT: multi-walled carbon nanotube.

Table 2.

The thermo physical properties of CNTs and nanofluids of some base fluids.

Physical properties		Specific heat, $c_{p} (J / K)$	Density, $ρ (kg / m^{3})$	Thermal conductivity, $k (W / mK)$
Base fluid	Water	$4.197 \times 10^{3}$	$9.97 \times 10^{2}$	$6.13 \times 10^{- 1}$
	Kerosene (lamp) oil	$2.090 \times 10^{3}$	$7.83 \times 10^{2}$	$1.45 \times 10^{- 1}$
	Engine oil	$1.910 \times 10^{3}$	$8.84 \times 10^{2}$	$1.44 \times 10^{- 1}$
Nanofluids	SWCNT	$4.25 \times 10^{2}$	$2.6 \times 10^{3}$	$6.6 \times 10^{3}$
	MWCNT	$7.96 \times 10^{2}$	$1.6 \times 10^{3}$	$3 \times 10^{3}$

CNTs: carbon nanotubes; SWCNT: single-walled carbon nanotube; MWCNT: multi-walled carbon nanotube.

Table 3.

Thermal conductivity $(k_{nf})$ of CNTs with different volume fractions (ϕ).

$φ$	$k_{nf}$ for SWCNT	$k_{nf}$ for MWCNT
0	$1.45 \times 10^{- 1}$	$1.45 \times 10^{- 1}$
0.01	$1.74 \times 10^{- 1}$	$1.72 \times 10^{- 1}$
0.02	$2.04 \times 10^{- 1}$	$2 \times 10^{- 1}$
0.03	$2.35 \times 10^{- 1}$	$2.28 \times 10^{- 1}$
0.04	$2.66 \times 1 0^{- 1}$	$2.57 \times 1 0^{- 1}$

CNTs: carbon nanotubes; SWCNT: single-walled carbon nanotube; MWCNT: multi-walled carbon nanotube.

Conclusion

Entropy generation investigation in three-dimensional squeezing flow of CNTs based on water $(H_{2} O)$ has been done. The influences of thermal radiation, viscous dissipation, and applied magnetic field on nanofluid are taken into account. SWCNTs and MWCNTs are used in this model. The impacts of embedded of parameters are shown graphically. The concluding remarks of this investigation are listed below:

The velocity profile $f (λ)$ in y-direction escalates with the escalation in $(φ)$ and $(A)$ , and $(Ω)$ gives dual behavior with the escalating values, while with the escalation in $(M)$ , the velocity function gives decreasing behavior.

The velocity profile $f' (λ)$ in x-direction reduces with the escalation in $(A)$ and gives the oscillatory behavior with the enhancement in $(φ)$ , $(Ω)$ , and $(M)$ .

The velocity profile $g (λ)$ escalates with the escalation in $(φ)$ and reduces with the escalation in $(Ω)$ , $(A)$ , and $(M)$ .

Generally, the velocity profiles escalate with the decrease in $(S > 0)$ and reduce with escalation in $(S < 0)$ .

The temperature function $θ (λ)$ escalates with the escalation in $(S)$ , $(Ec)$ , and $(E c_{m})$ and reduces with the escalation in $(φ)$ and $(R)$ .

The entropy generation $(Ns)$ escalates with the escalation in $(R)$ , while $(Ec)$ , $(E c_{m})$ , $(M)$ , and $(\Pr)$ have dual behavior on $(Ns)$ .

The Bejan number $(Be)$ escalates with the escalation in $(R)$ , $(Ec)$ , and $(E c_{m})$ , while declines with $(M)$ and $(\Pr)$ .

The skin fraction coefficient $(Cf {Re}_{x}^{1 / 2})_{lower}$ escalates with the escalation in $(S)$ , $(M)$ , and $(R)$ , while reduces with the escalation in $(Ω)$ .

The skin fraction coefficient $(Cf {Re}_{x}^{1 / 2})_{upper}$ enhances with the enhancement in $(Ω)$ while reduces with $(S)$ and $(M)$ .

The local Nusselt number $(Nu {Re}_{x}^{1 / 2})_{lower}$ escalates with the enhancement in $(R)$ , while reduces with $(S)$ .

The local Nusselt number $(Nu {Re}_{x}^{1 / 2})_{upper}$ increases with the escalation in $(R)$ , while reduces with $(S)$ .

Footnotes

Handling Editor: Oronzio Manca

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

Abdullah Dawar

Zahir Shah

Waris Khan

References

Choi

SUS

. Enhancing thermal conductivity of fluids with nanoparticles. In: Proceedings of the ASME international mechanical engineering congress and exposition, San Francisco, CA, 12–17 November 1995, pp.99–105. New York: ASME.

Kang

Kim

JM.

Estimation of thermal conductivity of nanofluid using experimental effective particle volume. Exp. Heat Transfer 2006; 19: 181–191.

Ganji

Hatami

Squeezing Cu–water nanofluid flow analysis between parallel plates by DTM-Padé method. J Mol Liq 2014; 193: 37–44.

Haq

Nadeem

Khan

et al . Convective heat transfer in MHD slips flow over a stretching surface in the presence of carbon nanotubes. Phys B Condens Matter 2015; 457: 40–47.

Liu

Lin

MCC

et al . Enhancement of thermal conductivity with carbon nanotube for nanofluids. Int Comm Heat Mass Trans 2005; 32: 1202–1210.

Nadeem

Lee

Boundary layer flow of nanofluid over an exponentially stretching surface. Nanoscale Res Lett 2012; 94: 7: 94.

Ellahi

Raza

Vafai

Series solutions of non-Newtonian nanofluids with Reynolds’ model and Vogel’s model by means of the homotopy analysis method. Math Comput Model 2012; 55: 1876–1891.

Nadeem

Haq

Khan

ZH.

Numerical study of MHD boundary layer flow of a Maxwell fluid past a stretching sheet in the presence of nanoparticles. J Taiwan Inst Chem Eng 2014; 45: 121–126.

Nadeem

Haq

RU.

MHD boundary layer flow of a nanofluid past a porous shrinking sheet with thermal radiation. J Aerospace Eng 2015; 28: 02995.

10.

Nadeem

Haq

RU.

Effect of thermal radiation for megnetohydrodynamic boundary layer flow of a nanofluid past a stretching sheet with convective boundary conditions. J Comput Theor Nanosci 2013; 11: 32–40.

11.

Sheikholeslami

Numerical investigation of nanofluid free convection under the influence of electric field in a porous enclosure. J Mol Liq 2018; 249: 1212–1221.

12.

Shah

Islam

Ayaz

et al . Radiative heat and mass transfer analysis of micropolar nanofluid flow of Casson fluid between two rotating parallel plates with effects of hall current. J Heat Trans 2018; 141: 02240.

13.

Nasir

Islam

Gul

et al . Three-dimensional rotating flow of MHD single wall carbon nanotubes over a stretching sheet in presence of thermal radiation. Appl Nanosci 2018; 8: 1361–1378.

14.

Dawar

Shah

Islam

et al . Magnetohydrodynamic CNTs Casson nanofluid and radiative heat transfer in a rotating channels. J Phys Res Appl 2018; 1: 17–32.

15.

Sheikholeslami

CuO-water nanofluid flow due to magnetic field inside a porous media considering Brownian motion. J Mol Liq 2018; 249: 921–929.

16.

Sheikholeslami

Numerical investigation for CuO-H2O nanofluid flow in a porous channel with magnetic field using mesoscopic method. J Mol Liq 2018; 249: 739–746.

17.

Sheikholeslami

Shehzad

SA.

Numerical analysis of Fe₃O₄–H₂O nanofluid flow in permeable media under the effect of external magnetic source. Int J Heat Mass Tran 2018; 118: 182–192.

18.

Sheikholeslami

Influence of Lorentz forces on nanofluid flow in a porous cavity by means of Non- Darcy model. Eng Comput 2017; 34: 2651–2667.

19.

Sheikholeslami

Rokni

HB.

Numerical simulation for impact of Coulomb force on nanofluid heat transfer in a porous enclosure in presence of thermal radiation. Int J Heat Mass Tran 2018; 118: 823–831.

20.

Sheikholeslami

Shehzad

et al . Numerical modeling for Alumina nanofluid magnetohydrodynamic convective heat transfer in a permeable medium using Darcy law. Int J Heat Mass Tran 2018; 127: 614–622.

21.

Sheikholeslami

Shafee

Lorentz forces effect on NEPCM heat transfer during solidification in a porous energy storage system. Int J Heat Mass Tran 2018; 127: 665–674.

22.

Sheikholeslami

Lattice Boltzmann method simulation of MHD non-Darcy nanofluid free convection. Physics B 2017; 516: 55–71.

23.

Sheikholeslami

Influence of magnetic field on nanofluid free convection in an open porous cavity by means of Lattice Boltzmann method. J Mol Liq 2017; 234: 364–374.

24.

Sheikholeslami

Magnetohydrodynamic nanofluid forced convection in a porous lid driven cubic cavity using Lattice Boltzmann method. J Mol Liq 2017; 231: 555–565.

25.

Sheikholeslami

CuO-water nanofluid free convection in a porous cavity considering Darcy law. Eur Phys J Plus 2017; 132: 55.

26.

Sheikholeslami

Numerical investigation of MHD nanofluid free convective heat transfer in a porous tilted enclosure. Eng Comput 2017: 34: 1939–1955.

27.

Bejan

A study of entropy generation in fundamental convective heat transfer. J Heat Trans 1979; 101: 718–725.

28.

Sarojamma

Vajravelu

Sreelakshmi

A study on entropy generation on thin film flow over an unsteady stretching sheet under the influence of magnetic field, thermocapillarity, thermal radiation and internal heat generation/absorption. Commun Numer Anal 2017; 2: 141–156.

29.

Soomro

Haq

Khan

et al . Numerical study of entropy generation in MHD water-based carbon nanotubes along an inclined permeable surface. Eur Phys J Plus 2017; 132: 412.

30.

Mansour

Sahin

Entropy generation in developing laminar fluid flow through a circular pipe with variable properties. Heat Mass Trans 2005; 42: 1–11.

31.

Rashidi

Mohammadi

Abbasbandy

et al . Entropy generation analysis for stagnation point flow in a porous medium over a permeable stretching surface. J Appl Fluid Mech 2015; 8: 753–765.

32.

Bég

Rashidi

Kavyani

et al . Entropy generation in hydromagnetic convective on Karman swirling flow. Int J Appl Math Mech 2013; 9: 37–65.

33.

Srinivasacharya

Bindu

KH.

Entropy generation in a micropolar fluid flow through an inclined channel. Alexand Eng J 2016; 55: 973–982.

34.

Qing

Bhatti

Abbas

et al . Entropy generation on MHD Casson nanofluid flow over a porous stretching/shrinking surface. Entropy 2016; 18: 123.

35.

Ishaq

Gohar

Shah

et al . Entropy generation on nanofluid thin film flow of Eyring–powell fluid with thermal radiation and MHD effect on an unsteady porous stretching sheet. Entropy 2018; 20: 412.

36.

Hayat

Rafiq

Ahmad

et al . Entropy generation analysis for peristaltic flow of nanoparticles in a rotating frame. Int J Heat Mass Tran 2017; 108: 1775–1786.

37.

Nouri

Pasandideh-Fard

Oboodi

et al . Entropy generation analysis of nanofluid flow over a spherical heat source inside a channel with sudden expansion and contraction. Int J Heat Mass Tran 2018; 116: 1036–1043.

38.

Dalir

Dehsara

Nourazar

SS.

Entropy analysis for magnetohydrodynamic flow and heat transfer of a Jeffrey nanofluid over a stretching sheet. Energy 2015; 79: 351–362.

39.

Ahmad

Khan

Hayat

et al . Entropy generation optimization and unsteady squeezing flow of viscous fluid with five different shapes of nanoparticles. Colloid Surf A: Physico Eng 2018; 554: 197–210.

40.

Rashidi

Abelman

Mehr

NF.

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

41.

Rashidi

Kavyani

Abelman

Investigation of entropy generation in MHD and slip flow over a rotating porous disk with variable. Int J Heat Mass Tran 2014; 70: 892–917.

42.

Rashidi

Ali

Freidoonimehr

et al . Parametric analysis and optimization of entropy generation in unsteady MHD flow over a stretching rotating disk using artificial neural network and particle swarm optimization algorithm. Energy 2013; 55: 497–510.

43.

Afridi

Qasim

Makinde

OD.

Second law analysis of boundary layer flow with variable fluid properties. J Heat Trans 2017; 10: 104505.

44.

Das

Banu

Jana

et al . Entropy analysis on MHD pseudo-plastic nanofluid flow through a vertical porous channel with convective heating. Alexand Eng J 2015; 54: 325–337.

45.

Seth

Sarkar

Makinde

OD.

Combined free and forced convection couette-hartmann flow in a rotating channel with arbitrary conducting walls and hall effects. J Mech 2016; 5: 613–629.

46.

Mkwizu

Makinde

Nkansah-Gyekye

Numerical investigation into entropy generation in a transient generalized Couette flow of nanofluids with convective cooling. Sadhana Acad Proc Eng Sci 2015; 40: 2073–2093.

47.

Makinde

Khan

Aziz

On inherent irreversibility in Sakiadis flow of nanofluids. Int J Exergy 2013; 2: 159–174.

48.

Adesanya

Makinde

OD.

Entropy generation in couple stress fluid flow through porous channel with fluid slippage. Int J Exerg 2014; 15: 344–362.

49.

Sheremet

Oztop

Pop

et al . Analysis of entropy generation in natural convection of nanofluid inside a square cavity having hot solid block: Tiwari and Das’ model. Entropy 2016; 18: 9.

50.

Sheremet

Pop

Öztop

et al . Natural convection of nanofluid inside a wavy cavity with a non-uniform heating: entropy generation analysis. Int J Nume Meth Heat Fluid Flow 2017; 27: 958–980.

51.

Khan

Syed

Mohyud-Din

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.

52.

Shah

Dawar

Islam

et al . Darcy-Forchheimer flow of radiative carbon nanotubes with microstructure and inertial characteristics in the rotating frame. Case Stud Therm Eng 2018; 12: 823–832.

53.

Khan

Shah

Islam

et al . Darcy–Forchheimer flow of micropolar nanofluid between two plates in the rotating frame with non-uniform heat generation/absorption. Adv Mech Eng 2018; 10: 1–16.

54.

Shah

Bonyah

Islam

et al . Radiative MHD thin film flow of Williamson fluid over an unsteady permeable stretching. Heliyon 2018; 4: e00825.

55.

Jawad

Shah

Islam

et al . Darcy-Forchheimer flow of MHD nanofluid thin film flow with Joule dissipation and Navier’s partial slip. J Phys Commun 2018; 2: 11.

56.

Khan

Zuhra

Shah

et al . Slip flow of Eyring-Powell nanoliquid film containing graphene nanoparticles. AIP Adv 2018; 8: 115302.

57.

Khan

Nie

Shah

et al . Three-dimensional nanofluid flow with heat and mass transfer analysis over a linear stretching surface with convective boundary conditions. Appl Sci 2018; 8: 2244.

58.

Khan

Shah

Islam

et al . Darcy–Forchheimer flow of micropolar nanofluid between two plates in the rotating frame with non-uniform heat generation/absorption. Adv Mech Eng 2018; 10: 1–16.

Unsteady squeezing flow of magnetohydrodynamic carbon nanotube nanofluid in rotating channels with entropy generation and viscous dissipation

Abstract

Keywords

Introduction

Problem statement

Entropy analysis and Bejan number

Solution by HAM

Results and discussion

Velocity functions ( f ( λ ) , f ′ ( λ ) , and g ( λ ) )

Temperature function ( θ ( λ ) )

Entropy analysis (Ns) and Bejan number (Be)

Skin fraction (Cf Rex1/2) and local Nusselt number (Nu Rex1/2)

Conclusion

Footnotes

Declaration of conflicting interests

Funding

ORCID iDs

References

Velocity functions $(f (λ), f' (λ),$ and $g (λ))$

Temperature function $(θ (λ))$

Skin fraction (Cf Re_x^1/2) and local Nusselt number (Nu Re_x^1/2)