### Review Article

**Magnetohydrodynamic CNTs Casson Nanofluid and Radiative heat transfer in a Rotating Channels
**

**
**Abdullah Dawar^{1}, Zahir Shah^{2*}, Saeed Islam^{2}, Muhammad Idress^{3} and Waris Khan^{3}

Abdullah Dawar^{1}, Zahir Shah^{2*}, Saeed Islam^{2}, Muhammad Idress^{3} and Waris Khan^{3}

^{1}Department of Mathematics, *Q*urtuba University of Science and Information Technology, Peshawar 25000, Pakistan

^{2}Department of Mathematics, Abdul Wali Khan University, Mardan 23200, KP, Pakistan

^{3}Department of Mathematics, Islamia College University, Mardan 25000, KP, Pakistan

***Address for Correspondence:** Zahir Shah, Department of Mathematics, Abdul Wali Khan University, Mardan 23200, KP, Pakistan, Email: zahir1987@yahoo.com

**Dates:** **Submitted:** 27 July 2018; **Approved:** 16 August 2018; **Published:** 17 August 2018

**How to cite this article:** Dawar A, Shah Z, Islam S, Idress M, Khan W. Magnetohydrodynamic CNTs Casson Nanofluid and Radiative heat transfer in a Rotating Channels. J Phys Res Appl. 2018; 1: 017-032. DOI: 10.29328/journal.ijpra.1001002

**Copyright:** © 2018 Dawar A, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

**Keywords:** SWCNTs; MWCNTs; Nanoparticles; Casson fluids; Thermal radiation; MHD; Parallel Plates; Rotating System; HAM

### Abstract

The main purpose of this investigation is to inspect the innovative conception of the magneto hydrodynamic (MHD) nanoparticles of single wall carbon nanotubes base on the fluids (water, engine oil, and ethylene, glycol and kerosene oil) between two rotating parallel plates. Carbon nanotubes (CNTs) parade sole assets due to their rare structure. Such structure has significant optical and electronics features, wonderful strength and elasticity, and high thermal and chemical permanence. The heat exchange phenomena is deliberated subject to thermal radiation. Kerosene oil is taken as based nano fluids because of its unique attention due to their advanced thermal conductivities, exclusive features, and applications. The fluid flow is presumed in steady state. With the help of suitable resemblance variables, the fundamental leading equations have been converted to a set of differential equations. To obtain the solution of the modeled problem, the homotopic approach has been used. The influence of imbedded physical variables upon the velocities and temperature profiles are defined and deliberated through graphs. Moreover, for the several values of relevant variables, the skin fraction coefficient and local Nusselt number are tabulated. Plots have been presented in order to examine how the velocities and temperature profile get affected by various flow parameters.

### Introduction

A nanometer-scale tube-like structure is called nanotube. Carbon nanotubes (CNTs) were the first to be discovered in 1952. CNTs are commonly divided into two types called Single-walled carbon nanotubes (SWCNTs) and multi-walled carbon nanotubes (MWCNTs). The potential of CNTs are catalysts, flat panel display, absorption, and shielding of electro-magnetic waves, nanoelectrodes, sensors, supercapicator and energy conversion. Haq et al. [1], Reported that SWCNTs has higher Nusselt number and skin friction than MWCNTs by considering water as a base fluid. Liu et al. [2], Deliberated that ethylene glycol with CNTs have higher thermal conductivities than ethylene glycol suspension without CNTs by studying synthetic engine oil and ethylene glycol in the existence of MWCNTs. The theory of nanofluid past an exponentially stretching sheet has been presented by Nadeem and Lee [3]. Over stretching/shrinking surfaces the mentioned authors [4-7] examined the impacts of nanoparticles for boundary layer flow. Choi [8] was the pioneer who presented the idea of nanofluid by dipping the nanometer-sized particle into the base fluid. For the study of nanofluid flow, the mathematical model was established by Boungiorno [9]. The recent experimental and theoretical study of Sheikholeslami on nanofluids and its applications with different behaviour, properties and effects with uses of different numerical and analytical techniques can be seen in [10^{-1}5]. Sheikholeslami and Rokni [16^{-1}9] have recently investigated Simulation of nanofluids and CuOeH2O nanofluid in a curved porous enclosure by using mesoscopic approach. Elias et al. [20], have examined the influence of shape on heat transmission. Liquid metals, plasmas, slat water, electrolytes are the examples of such magnetofluids. MHD was presented by Hannes Alfven [21]. Cooling of devices, MHD instruments, crystal growth, electromagnetic casting, and magnetic drug targeting are the real-world uses of MHD in the field of engineering and technologies. Recently Shah et al. [22-25], have studied MHD nanofluid with rotating system.

The determination of current research is to study the nanoparticles of SWCNTs and MWCNTs Casson Fluid based on the fluids (water, engine oil, ethylene glycol and kerosene oil) between two rotating parallel plate under the influence of magnetic parameter and thermal radiation. Kerosene oil is taken as base nanofluid. Casson [26] was the first one who proposed the Casson fluid model, this model characterizes a shear retreating fluid which is presumed to obligate infinite viscosity at zero rate of shear stress. Mehmood et al. [27,28], recently investigated Casson micro polar fluid over a stretching sheet with internal heat transmission by using numerical techniques. Singh Megahe et al. [29], inspected the liquid film flow of Casson fluid in the existence of varied heat flux using slip velocity. Abolbashari et al. [30], studied the Casson nanofluid with entropy generation. The other related study about Casson fluid can be seen in [31-36]. The effect of all embedding parameters has been studied graphically. The analytical result for velocities and temperature profiles are obtained using the HAM technique [37-42].

**Some recent application of CNTs**

CNTS used for

• Cancer treatment

• Cardiac autonomic regulation

• In drug delivery

• Platelet activation

• Tissue regeneration

• In solar cell

**Mathematical modeling**

Consider the flows of CNTs nanofluid between two parallel plates. The distance between the lower and upper plates is labeled with h. SWCNTs and MWCNTs are used as nano-scale materials where kerosene oil is a base liquid. Through thermal radiation, the heat transportation mechanism is examined. Around the y-axis the plates are rotated with a constant angular velocity *γ*. It should be noted that *γ* > 0 indicates that both plates rotate in the same direction, *γ* > 0 indicates that both plates rotate in the opposite directions, *γ* > 0 is for the static case. The rotation of the lower plate which is moving with velocity U w = cx (c>0) is quicker than the upper plate. A coordinate system (*x, y, z*) is chosen in such a way that the x-axis is parallel to the plates, the y-axis is vertical to the plates, and the z-axis is normal to the xy plane. The plates are positioned at *y = 0* and *y = h*. With the help of two forces with same magnitude but opposite direction, the lower plate is being kept stretchable so the position (0, 0, 0) cannot changes. The fluid flow and heat transfer is supposed in steady state which is incompressible, laminar and stable. Along y direction, the magnetic field B0 is substituted with which the fluid is rotating as shown in Figure 1. The rheological model that illustrates the Casson fluid is known as [34-36]:

${{\rm T}}_{ij}=2\{\begin{array}{l}\left({\mu}_{{\rm B}}+\frac{{{\rm P}}_{y}}{\sqrt{2\kappa}}\right){m}_{ij},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\kappa >{\kappa}_{c},\\ \left({\mu}_{{\rm B}}+\frac{{{\rm P}}_{y}}{\sqrt{2{\kappa}_{c}}}\right){m}_{ij},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\kappa >{\kappa}_{c}.\end{array}\text{(1)}$

T* _{ij}* Denotes Cauchy stress tensor,

*µB*denotes the dynamic viscosity of the Casson fluid,

*k = m*is the square of components of strain rate P

_{ij}.m_{ij}_{y}is the yield stress of the fluid and

*k*represent critical value of

_{c}*k*

The governing equations for the state of problem are [22-24,37]

$\frac{\partial u}{\partial x}+\frac{\partial u}{\partial y}=0,\text{(2)}$

${\rho}_{nf}\left(u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}+2\gamma w\right)=-\frac{\partial {P}^{\ast}}{\partial x}+\left(1+\frac{1}{\beta}\right){\mu}_{nf}\left(\frac{{\partial}^{2}u}{\partial {x}^{2}}+\frac{{\partial}^{2}u}{\partial {y}^{2}}\right)-{\sigma}_{nf}{B}_{0}^{2}u,\text{(3)}$

${\rho}_{nf}\left(u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}\right)=-\frac{\partial {P}^{\ast}}{\partial y}+{\mu}_{nf}\left(1+\frac{1}{\beta}\right)\left(\frac{{\partial}^{2}u}{\partial {x}^{2}}+\frac{{\partial}^{2}u}{\partial {y}^{2}}\right),\text{(4)}$

${\rho}_{nf}\left(u\frac{\partial w}{\partial x}+v\frac{\partial w}{\partial y}-2\gamma u\right)={\mu}_{nf}\left(1+\frac{1}{\beta}\right)\left(\frac{{\partial}^{2}w}{\partial {x}^{2}}+\frac{{\partial}^{2}w}{\partial {y}^{2}}\right)+{\sigma}_{nf}{B}_{0}^{2}w,\text{(5)}$

Where ${P}^{\ast}=P-\frac{{\gamma}^{2}{x}^{2}}{2},{\sigma}_{nf},{\mu}_{nf}$
represents the meshes cross flow beside z-axis and *β* is the parameter of Casson fluid. Mathematically, the phenomenon of heat transmission can be indicated as

$u\frac{\partial {\rm T}}{\partial x}+v\frac{\partial {\rm T}}{\partial y}=\frac{{k}_{nf}}{{(\rho c)}_{nf}}\left(\frac{{\partial}^{2}{\rm T}}{\partial {x}^{2}}+\frac{{\partial}^{2}{\rm T}}{\partial {y}^{2}}\right)-\frac{\partial Qr}{\partial y}.\text{(6)}$

${\rm T},{\alpha}_{nf}=\frac{{k}_{nf}}{{(\rho c)}_{nf}},Qr$ Signifies the fluid temperature, thermal diffusivity and the radiative heat flux respectively. The radiative heat flux is defined as $Qr=-\frac{4{\sigma}^{\ast}}{3k}\frac{\partial {{\rm T}}^{4}}{\partial y},\text{(7)}$

Where ${\sigma}^{\ast},k$ indicates the absorption coefficient and Stefan Boltzmann constant respectively.

As ${{\rm T}}^{4}=4{{\rm T}}_{h}^{3}{\rm T}-3{{\rm T}}_{h}^{4}$ the equation (6) becomes

$u\frac{\partial {\rm T}}{\partial x}+v\frac{\partial {\rm T}}{\partial y}=\frac{{k}_{nf}}{{(\rho c)}_{nf}}\left(1+\frac{16{\sigma}^{\ast}{{\rm T}}_{h}^{3}}{3k}\right)\frac{{\partial}^{2}{\rm T}}{\partial {y}^{2}}.\text{(8)}$

The density, heat capacity, and dynamic viscosity of the nanofluid through mathematical equations are [40].

$\begin{array}{l}{(}_{\rho}=\left(1-\phi \right){\rho}_{f}+\phi {\rho}_{CNT},{\rho}_{nf}=\left(1-\phi \right){\left(\rho c\right)}_{f}+\phi {\left(\rho c\right)}_{CNT},\\ {\mu}_{nf}=\frac{{\mu}_{f}}{{\left(1-\phi \right)}^{2.5}},{k}_{nf}={k}_{f}\left(\frac{1-\phi +2\phi \frac{{k}_{CNT}}{{k}_{CNT-{k}_{f}}}\mathrm{ln}\left(\frac{{k}_{CNT+{k}_{f}}}{2{k}_{f}}\right)}{1-\phi +2\phi \frac{{k}_{f}}{{k}_{CNT-{k}_{f}}}\mathrm{ln}\left(\frac{{k}_{CNT+{k}_{f}}}{2{k}_{f}}\right)}\right).\text{(9)}\end{array}$

Where, ${\mu}_{f},\phi ,{k}_{nf}$
is the dynamic viscosity of base fluid, nanoparticle volumetric friction, and thermal conductivity respectively. The subscripts *CNT, f, nf* represent carbon nanotubes, base fluid, and nano fluid respectively. T_{0} is the lower temperature at upper wall and T_{h} is the higher temperature at lower wall. T_{h} is retained higher than T_{0} i.e. T_{h} > T_{0} . The boundary conditions for the system can be defined as [37].

$\begin{array}{l}\overrightarrow{u}={U}_{w}=cx,\overrightarrow{v}=0,\overrightarrow{w}=0,{\rm T}={{\rm T}}_{h}\text{at}y=0\text{,(10)}\\ \overrightarrow{u}\text{=0,}\overrightarrow{v}{\text{=-Z}}_{0}\text{,}\overrightarrow{w}=0,{\rm T}={{\rm T}}_{0}\text{at}y=h.\text{}\end{array}$

In equation (9) Z_{0} is the uniform suction/injection velocity at the upper wall. If (Z_{0} > 0) than it is called uniform suction velocity and if (Z_{0} < 0) than it is called uniform injection velocity.

Excluding the pressure gradient and introducing the similarity variables, we have [22-24].

$\overrightarrow{u}=cx{f}^{\prime}\left(\Gamma \right),\overrightarrow{v}=-chf\left(\Gamma \right),\overrightarrow{w}=cxg\left(\Gamma \right),\theta \left(\Gamma \right)=\frac{{\rm T}-{{\rm T}}_{0}}{{{\rm T}}_{h}-{{\rm T}}_{0}},\Gamma =\frac{y}{h}.\text{(11)}$

$\begin{array}{l}\left(1+\frac{1}{\beta}\right){f}^{iv}\left(\Gamma \right)-{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{1}\left({f}^{\prime}\left(\Gamma \right){f}^{\u2033}\left(\Gamma \right)-f\left(\Gamma \right){f}^{\u2034}\left(\Gamma \right)\right)\text{(12)}\\ -{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{2}{g}^{\prime}\left(\Gamma \right)-\left(\text{\hspace{0.17em}}\frac{{\sigma}_{nf}}{{\sigma}_{f}}\right){\rm M}{f}^{\u2033}\left(\Gamma \right)=0,\end{array}$

$\begin{array}{l}\left(1+\frac{1}{\beta}\right){g}^{\u2033}\left(\Gamma \right)-{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{1}\left(g\left(\Gamma \right){f}^{\prime}\left(\Gamma \right)-f\left(\Gamma \right){g}^{\prime}\left(\Gamma \right)\right)\text{(13)}\\ +2{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{2}{f}^{\prime}\left(\Gamma \right)+\left(\text{\hspace{0.17em}}\frac{{\sigma}_{nf}}{{\sigma}_{f}}\right){\rm M}g\left(\Gamma \right)=0,\end{array}$

$\left(1+Rd\right){\theta}^{\u2033}\left(\Gamma \right)+{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]\frac{\mathrm{Pr}{A}_{1}}{\left(\frac{{k}_{nf}}{{k}_{f}}\right)}f\left(\Gamma \right){\theta}^{\prime}\left(\Gamma \right)=0,\text{(14)}$

The relevant boundary conditions are $\begin{array}{l}f\left(0\right)=0,{f}^{\prime}\left(0\right)=1,g\left(0\right)=0,\theta \left(0\right)=1,\text{(15)}\\ f\left(1\right)=Q,{f}^{\prime}\left(1\right)=0,g\left(1\right)=0,\theta \left(1\right)=0.\end{array}$

In equation (14), $Q=\raisebox{1ex}{${Z}_{0}$}\!\left/ \!\raisebox{-1ex}{$ch$}\right.$ is the suction and injection parameter. If $\left(Q>0\right)$ then it is called suction parameter and if $\left(Q<0\right)$ then it is called injection parameter. The above non-dimensional system of equations has the following parameters.

$\begin{array}{l}{A}_{1}=\frac{c{h}^{2}}{{\nu}_{f}}\left(Reynold\text{'}snumber\right),{A}_{2}=\frac{\gamma {h}^{2}}{{\nu}_{f}}\left(rotationparameter\right),\\ \mathrm{Pr}=\frac{{\mu}_{f}{c}_{p}}{{k}_{f}}\left(prandtlnumber\right),Rd=\frac{16{\sigma}^{\ast}{{\rm T}}_{h}^{3}}{3k{k}^{*}}\left(radiationparameter\right),\\ {\rm M}=\frac{{\sigma}_{f}{B}_{0}^{2}{h}^{2}}{{\rho}_{f}{\nu}_{f}}\left(magneticparameter\right),\frac{{k}_{nf}}{{k}_{f}}=\frac{{k}_{s}+2{k}_{f}-2\varphi \left({k}_{f}-{k}_{s}\right)}{{k}_{s}+2{k}_{f}+2\varphi \left({k}_{f}-{k}_{s}\right)},\text{(16)}\\ \text{\hspace{0.17em}}\frac{{\sigma}_{nf}}{{\sigma}_{f}}=1+\frac{3\left(\text{\hspace{0.17em}}\frac{{\sigma}_{s}}{{\sigma}_{f}}-1\right)\varphi}{\left(\text{\hspace{0.17em}}\frac{{\sigma}_{s}}{{\sigma}_{f}}+2\right)-\left(\text{\hspace{0.17em}}\frac{{\sigma}_{s}}{{\sigma}_{f}}-1\right)\varphi},\end{array}$

The skin friction coefficient and local Nusselt number are defined as [38]:

${\tilde{C}}_{f}=\frac{{\mu}_{nf}}{{\mu}_{f}}{f}^{\u2033}(0)\text{(17)}$

$N{u}_{x}=-\frac{{k}_{nf}}{{k}_{f}}\left(1+Rd\right){\theta}^{\prime}\left(0\right).\text{(18)}$

**Solution by HAM**

We use HAM [39-42] to solve Eqs. (12, 13, and 14) with boundary conditions (15) by the succeeding process.

The primary suppositions are chosen as follows:

The linear operators are chosen as ${\text{L}}_{\overline{f}}{\text{,L}}_{\overline{g}}{\text{andL}}_{\overline{\theta}}$

${L}_{\overline{f}}(\overline{f})=\overline{f}\text{'}\text{'}\text{'},{\text{L}}_{\overline{g}}\text{( g \xaf )= g \xaf ''},{\text{L}}_{\overline{\theta}}\text{(}\overline{\theta}\text{)=}\overline{{\theta}^{\u2033}},$

Which have the succeeding properties:

${L}_{\overline{f}}({b}_{1}+{b}_{2}\Gamma +{b}_{3}{\Gamma}^{2})=0,{\text{L}}_{\overline{g}}({b}_{4}+{b}_{5}\Gamma )=0,{\text{L}}_{\overline{\theta}}({b}_{6}+{b}_{7}\Gamma )=0,\text{(21)}$

Where ${b}_{i}(i=1-7)$

The consequence non-linear operators ${N}_{\overline{f}},{N}_{\overline{g}}$ and ${N}_{\overline{\theta}}$ are indicated as:

$\begin{array}{l}{N}_{\overline{f}}\left[\overline{f}(\Gamma ;\psi ),\overline{g}(\Gamma ;\psi )\right]=\left(1+\frac{1}{\beta}\right)\frac{{\partial}^{4}\overline{f}(\Gamma ;\psi )}{\partial {\Gamma}^{4}}\\ -{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{1}\left\{\begin{array}{l}\frac{\partial \overline{f}(\Gamma ;\psi )}{\partial \Gamma}\frac{{\partial}^{2}\overline{f}(\Gamma ;\psi )}{\partial {\Gamma}^{2}}-\\ \overline{f}(\Gamma ;\psi )\frac{{\partial}^{3}\overline{f}(\Gamma ;\psi )}{\partial {\Gamma}^{3}}\end{array}\right\}\text{(22)}\\ -{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{2}\frac{\partial \overline{g}(\Gamma ;\psi )}{\partial \Gamma}-\left(\frac{{\sigma}_{nf}}{{\sigma}_{f}}\right)M\frac{{\partial}^{2}\overline{f}(\Gamma ;\psi )}{\partial {\Gamma}^{2}},\end{array}$

$\begin{array}{l}{N}_{\overline{g}}\left[\overline{f}(\Gamma ;\psi ),\overline{g}(\Gamma ;\psi )\right]=\left(1+\frac{1}{\beta}\right)\frac{{\partial}^{2}\overline{g}(\Gamma ;\psi )}{\partial {\Gamma}^{2}}\\ -{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{1}\left(\begin{array}{l}\overline{g}(\Gamma ;\psi )\frac{\partial \overline{f}(\Gamma ;\psi )}{\partial \Gamma}\\ -\overline{f}(\Gamma ;\psi )\frac{\partial \overline{g}(\Gamma ;\psi )}{\partial \Gamma}\end{array}\right)\text{(23)}\\ +2{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{2}\frac{\partial \overline{f}(\Gamma ;\psi )}{\partial \Gamma}+\left(\frac{{\sigma}_{nf}}{{\sigma}_{f}}\right)M\overline{g}(\Gamma ;\psi ),\end{array}$

${N}_{\overline{\theta}}\left[\overline{f}(\Gamma ;\psi ),\overline{\theta}(\Gamma ;\psi )\right]=(1+Rd)\frac{{\partial}^{2}\overline{\theta}(\Gamma ;\psi )}{\partial {\Gamma}^{2}}+\left[\left(1-\phi \right)+\phi \frac{{\left(\rho {c}_{p}\right)}_{CNT}}{{\left(\rho {c}_{p}\right)}_{f}}\right]\frac{\mathrm{Pr}{A}_{1}}{\left(\frac{{k}_{nf}}{{k}_{f}}\right)}\overline{f}\left(\Gamma ,\psi \right)\frac{\partial \overline{\theta}(\Gamma ;\psi )}{\partial \Gamma}.\text{(24)}$

The zero^{th} order problems from In organize to solve Eqs. (12, 13, and 14) are:

$(1-\psi ){L}_{\overline{f}}\left[\overline{f}(\Gamma ;\psi )-{\overline{f}}_{0}(\Gamma )\right]=\psi {\hslash}_{\overline{f}}{N}_{\overline{f}}\left[\overline{f}(\Gamma ;\psi ),\overline{g}(\Gamma ;\psi )\right],\text{(25)}$

$(1-\psi ){L}_{\overline{g}}\left[\overline{g}(\Gamma ;\psi )-{\overline{g}}_{0}(\Gamma )\right]=\psi {\hslash}_{\overline{g}}{N}_{\overline{g}}\left[\overline{f}(\Gamma ;\psi ),\overline{g}(\Gamma ;\psi )\right],\text{(26)}$

$(1-\psi ){L}_{\overline{\theta}}\left[\overline{\theta}(\Gamma ;\psi )-{\overline{\theta}}_{0}(\Gamma )\right]=\psi {\hslash}_{\overline{\theta}}{N}_{\overline{\theta}}\left[\overline{f}(\Gamma ;\psi ),\overline{\theta}(\Gamma ;\psi )\right].\text{(27)}$

The equivalent boundary conditions are:

$\begin{array}{l}{\overline{f}(\Gamma ;\psi )|}_{\Gamma =0}=0,\text{}{\frac{\partial \overline{f}(\Gamma ;\psi )}{\partial \Gamma}|}_{\Gamma =0}=1,\text{}{\overline{f}(\Gamma ;\psi )|}_{\Gamma =1}=Q,\text{}{\frac{\partial \overline{f}(\Gamma ;\psi )}{\partial \Gamma}|}_{\Gamma =1}=0,\\ \text{}{\overline{g}(\Gamma ;\psi )|}_{\Gamma =0}=0,\text{}{\overline{g}(\Gamma ;\psi )|}_{\Gamma =1}=0,\text{(28)}\\ {\overline{\theta}(\Gamma ;\psi )|}_{\Gamma =0}=1,\text{}{\overline{\theta}(\Gamma ;\psi )|}_{\Gamma =1}=0,\end{array}$

Where $\psi \in [0,1]$ is the emerging parameter? In case of $\psi =0\text{and}\psi =1$ we have:

Expanding $\overline{f}(\Gamma ;\psi ),\overline{g}(\Gamma ;\psi )\text{and}\overline{\theta}(\Gamma ;\psi )$ by Taylor’s series

$\begin{array}{l}\overline{f}(\Gamma ;\psi )={\overline{f}}_{0}(\Gamma )+{\displaystyle \sum _{q=1}^{\infty}{\overline{f}}_{q}(\Gamma ){\psi}^{q}},\text{}\\ \overline{g}(\Gamma ;\psi )={\overline{g}}_{0}(\Gamma )+{\displaystyle \sum _{q=1}^{\infty}{\overline{g}}_{q}(\Gamma ){\psi}^{q}},\text{(30)}\\ \overline{\theta}(\Gamma ;\psi )=\text{}{\overline{\theta}}_{0}(\Gamma )+{\displaystyle \sum _{q=1}^{\infty}{\overline{\theta}}_{q}(\Gamma ){\psi}^{q}}.\end{array}$

Where

As the series (31) converges at, $\psi =1$ changing $\psi =1$ in (31), we get:

$\begin{array}{l}\overline{f}(\Gamma )={\overline{f}}_{0}(\Gamma )+{\displaystyle \sum _{q=1}^{\infty}{\overline{f}}_{q}(\Gamma )},\text{}\\ \overline{g}(\Gamma )={\overline{g}}_{0}(\Gamma )+{\displaystyle \sum _{q=1}^{\infty}{\overline{g}}_{q}(\Gamma )},\text{(32)}\\ \overline{\theta}(\Gamma )={\overline{\theta}}_{0}(\Gamma )+{\displaystyle \sum _{q=1}^{\infty}{\overline{\theta}}_{q}(\Gamma )}.\end{array}$

The ${q}^{th}$ -order problem gratifies the following:

$\begin{array}{l}{L}_{\overline{f}}\left[{\overline{f}}_{q}(\Gamma )-{\chi}_{q}{\overline{f}}_{q-1}(\Gamma )\right]={\hslash}_{\overline{f}}{U}_{q}^{\overline{f}}(\Gamma ),\\ {L}_{\overline{g}}\left[{\overline{g}}_{q}(\Gamma )-{\chi}_{q}{\overline{g}}_{q-1}(\Gamma )\right]={\hslash}_{\overline{g}}{U}_{q}^{\overline{g}}(\Gamma ),\text{(33)}\\ {L}_{\overline{\theta}}\left[{\overline{\theta}}_{q}(\Gamma )-{\chi}_{q}{\overline{\theta}}_{q-1}(\Gamma )\right]={\hslash}_{\overline{\theta}}{U}_{q}^{\overline{\theta}}(\Gamma ).\end{array}$

The equivalent boundary conditions are:

$\begin{array}{l}{\overline{f}}_{q}(0)={\overline{{f}^{\prime}}}_{q}(0)={\overline{f}}_{q}(1)={\overline{{f}^{\prime}}}_{q}(1)=0,\text{}\\ {\overline{g}}_{q}(0)={\overline{g}}_{q}(1)=0,\text{(34)}\\ {\overline{\theta}}_{q}(0)={\overline{\theta}}_{q}(1)=0.\end{array}$

Here

$\begin{array}{l}{U}_{q}^{\overline{f}}(\Gamma )=\left(1+\frac{1}{\beta}\right){\overline{f}}_{q-1}^{iv}-{A}_{1}{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]\left({\displaystyle \sum _{k=0}^{q-1}{\overline{{f}^{\prime}}}_{q-1-k}}{\overline{{f}^{\u2033}}}_{k}-{\displaystyle \sum _{k=0}^{q-1}{\overline{f}}_{q-1-k}}{\overline{{f}^{\u2034}}}_{k}\right)\\ -{A}_{2}{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{\overline{g}}_{q-1}-\left(\frac{{\sigma}_{nf}}{{\sigma}_{f}}\right)M{\overline{{f}^{\u2033}}}_{q-1},\text{(35)}\end{array}$

$\begin{array}{l}{U}_{q}^{\overline{g}}(\Gamma )=\left(1+\frac{1}{\beta}\right){\overline{{g}^{\u2033}}}_{q-1}-{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{1}\left({\displaystyle \sum _{k=0}^{q-1}{\overline{g}}_{q-1-k}}{\overline{{f}^{\prime}}}_{k}-{\displaystyle \sum _{k=0}^{q-1}{\overline{f}}_{q-1-k}}{\overline{{g}^{\prime}}}_{k}\right)\\ +2{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \frac{{\rho}_{CNT}}{{\rho}_{f}}\right]{A}_{2}{\overline{{f}^{\prime}}}_{q-1}+\left(\frac{{\sigma}_{nf}}{{\sigma}_{f}}\right)M{\overline{g}}_{q-1},\text{(36)}\end{array}$

${U}_{q}^{\overline{\theta}}(\Gamma )=(1+Rd){\overline{{\theta}^{\u2033}}}_{q-1}+\left[\left(1-\phi \right)+\phi \frac{{\left(\rho {c}_{p}\right)}_{CNT}}{{\left(\rho {c}_{p}\right)}_{f}}\right]\frac{\mathrm{Pr}{A}_{1}}{\left(\frac{{k}_{nf}}{{k}_{f}}\right)}{\displaystyle \sum _{k=0}^{q-1}{\overline{f}}_{q-1-k}}{\overline{{\theta}^{\u2033}}}_{k}.\text{(37)}$

Where

${\chi}_{q}=\{\begin{array}{l}0,\text{if}\psi \le \text{1}\\ 1,\text{if}\psi \text{1}\end{array}\text{(38)}$

### Results and Discussion

In order to investigate the flow and heat transfer performance for both SWCNTs and MWCNTs based on kerosene nanoliquids between two rotating parallel plate, (Figure 3,6) are plotted. Figure 1 displays different structure of CNTs and Figure 2 shows physical shape of the flow. Figures 3-7 are plotted to see the impact of various parameters on velocity and temperature profiles for both SWCNTs and MWCNTs-kerosene nanoliquids. These parameters are nanoparticles volume friction *φ*, Reynolds number (), Rotation parameter (*A _{2}*), Magnetic parameter

*(M)*, Suction parameter (

*Q*> 0) and Injection parameter (

*Q*< 0) respectively. Figure 3 is plotted to see the comparison between the SWCNTs and MWCNTs with the escalating values of nanoparticle volume fraction

*φ*. From Figure 3 we observed that there is petty disparity in ${f}^{\prime}(\Gamma )$ with the escalating values of

*φ*. We also observed that the MWCNTs has comparatively much greater ${f}^{\prime}(\Gamma )$ as compared to SWCNTs. Figure 4 is plotted to see the impact of (

*A*) on ${f}^{\prime}(\Gamma )$ with the different escalating values. We observed from the figure, that with the different escalating values of (

_{1}*A*), the velocity profile changes ${f}^{\prime}(\Gamma )$ its behavior from rising to reducing at the medium of the two plates $(i.e.\Gamma =0.5)$ $(i.e.\Gamma =0.5)$. This is because of stretching of the lower plates. Figure 5 is plotted to see the impact of (

_{1}*A*) on ${f}^{\prime}(\Gamma )$. From figure 5, we observed that the disparity in ${f}^{\prime}(\Gamma )$ gives dual behavior. Within the region $(0\le \Gamma \le 0.5),$ due to numerous escalating values of (

_{2}*A*) the velocity profile ${f}^{\prime}(\Gamma )$ shows declining behavior, however within the region $(0.5<\Gamma \le 1)$, ${f}^{\prime}(\Gamma )$ gives escalating behavior for various escalating values of (

_{2}*A*). In addition, at the upper half of the channel (

_{2}*A*) gives more dominant variation for ${f}^{\prime}(\Gamma )$. Figure 6 is designed to see the influence of

_{2}*(M)*on ${f}^{\prime}(\Gamma )$. According to Lorentz force theory the magnetic field parameter has reverse effect on ${f}^{\prime}(\Gamma )$, that is ${f}^{\prime}(\Gamma )$ reduces with the escalation in

*(M)*. Figure 7 and 8 are plotted to see the impact of suction (

*Q*> 0) and injection (

*Q*< 0) parameters at the upper plate for ${f}^{\prime}(\Gamma ).$ Form figure 7, we see that ${f}^{\prime}(\Gamma ).$ escalates with suction (

*Q*< 0), that is ${f}^{\prime}(\Gamma ).$ increases with positive the values of (

*Q*), while form Figure 8, it is observed that ${f}^{\prime}(\Gamma ).$ reduces with injection (

*Q*< 0), that is ${f}^{\prime}(\Gamma ).$ reduces with the negative (

*Q*). It is fairly obvious that the presence of CNTs nanoparticles has improved the velocity function ${f}^{\prime}(\Gamma )$, while the velocity function ${f}^{\prime}(\Gamma )$ increases more when (

*Q*> 0) is present. But the decrement in velocity function ${f}^{\prime}(\Gamma )$ is due to the fact that (

*Q*< 0) absorbs the internal heat energy from the surface. Figure 9 is plotted to see the impact of (

*β*) on ${f}^{\prime}(\Gamma )$. From here we observed that the increasing values of (

*β*) shows reduction in ${f}^{\prime}(\Gamma )$. Figure 10 is plotted to see the impact of

*φ*on $g(\Gamma )$.

From Figure 10 we observed that there is petty disparity in $g(\Gamma )$ with the escalating values of *φ*. We also observed that the MWCNTs has comparatively much greater $g(\Gamma )$ as compared to SWCNTs. Figure 11 is plotted to see the impact of on (*A _{1}*) $g(\Gamma )$.
Figure 11 shows reduction in $g(\Gamma )$ between the two rotating plates with the increasing values of (

*A*) and the position of the maximum amount of $g(\Gamma )$ approaches to the stretching sheet. Figure 12 is plotted to see the impact of (

_{1}*A*) on $g(\Gamma )$. From Figure 12, we observed that the increasing values of (

_{2}*A*) gives decreasing behavior to $g(\Gamma )$. It is also observed that the disturbance in $g(\Gamma )$ is higher at the medium of the channels as compared to upper and lower surface of the channels. Figure 13 is plotted to see the impact of

_{2}*(M)*on $g(\Gamma )$. According to Lorentz force theory the magnetic field parameter has reverse effect on $g(\Gamma )$, that is $g(\Gamma )$ reduces with the escalation in

*(M)*. Figure 14 is plotted to see the comparison between the SWCNTs and MWCNTs with the escalating values of nanoparticle volume fraction

*φ*. From figure 14 we observed that there is petty disparity in $\theta (\Gamma )$ with the escalating values of

*φ*.

We also observed that the MWCNTs has comparatively much greater $\theta (\Gamma )$
as compared to SWCNTs. Figure 15 is plotted to see the impact of (*A _{1}*) on $\theta (\Gamma )$.
From here we observed that the escalating value of (

*A*) shows escalation in $\theta (\Gamma )$. As the distance from the surface escalates, $\theta (\Gamma )$ decreases. Figure 16 is plotted to see the impact of (Pr) on $\theta (\Gamma )$. From here we see that the escalating values of (Pr) shows reduction in $\theta (\Gamma )$. Physically, the nanofluids have a large thermal diffusivity with small (Pr), but this effect is revers for higher (Pr), therefore the temperature of liquid shows decreasing behavior. Figure 17 is plotted to see the impact of (

_{1}*Rd*) on $\theta (\Gamma )$ Thermal radiation has leading rule in heat transmission when the coefficient of convection heat transmission is small. From here we see that the escalating values of (

*Rd*) shows acceleration in $\theta (\Gamma )$.

### Discussion

Tables 1,2 are schemed to see the influences of different embedding parameters on skin fraction coefficient $\left({\tilde{C}}_{f}\right)$
and local Nusselt number (*Nu _{x}*). From Table 1, we see that the effect of

*A*on skin fraction coefficient. It is clear from the table that at 0.1, the table shows increasing behavior to 0.3 but from 0.3 to 0.5, the table values show decreasing behavior. At the medium of the channel the skin friction coefficient changes its behavior from increasing to decreasing. This is because of stretching of the lower plates. It is also clear form the table that the escalating values of

_{1}*A*shows reducing behavior. The escalating values of

_{2}*(M)*shows decreasing behavior in skin fraction coefficient $\left({\tilde{C}}_{f}\right)$. It is due to Lorentz force theory. The increasing and decreasing values of suction/injection parameter (

*Q*) show two different behaviors. The escalating values of injection parameter (

*Q*< 0) show decreasing behavior and also the escalating values of suction parameter (

*Q*> 0) show increasing behavior in skin fraction $\left({\tilde{C}}_{f}\right)$. From Table 2, we see that the increasing values of (

*A*) and positive

_{1}*Q*(

*i.e.*> 0) show increasing behavior in local Nusselt number (

*Q**Nu*), while (

_{x}*A*),

_{2}*(M)*and negative

*Q*(

*i.e.*> 0) show decreasing behavior in local Nusselt number (

*Q**Nu*). Tables 3-5 are schemed to study the Physical properties of CNTs, thermo physical properties CNTs and nanofluids of some base fluids, Thermal conductivity

_{x}*(k*of CNTs with different volume fraction

_{nf})*φ*respectively.

### Conclusion

The investigation for the two-dimensional flow of kerosene oil based nanofluid over an inclined stretching sheet with suction/injection, MHD and radiative heat flux effects are examined. SWCNTs and MWCNTs are used in this model. With the help of similarity variables, the system of governing partial differential equations is changed into ordinary differential equations. The impacts of embedded of parameters are shown graphically. The impacts of skin friction coefficient and local Nusselt number are shown through Table 1 and 2. On the achieved study, the key remarks are listed below.

1. The velocity function ${f}^{\prime}(\Gamma )$
escalates with the escalation in *φ*, (*β*) and suction parameter (*i.e. Q > 0*), and reduces with the escalation in

*A*,

_{1}*(M)*and injection parameter (

*i.e.*).

*Q*> 02. The velocity function ${f}^{\prime}(\Gamma )$
shows reducing behavior with escalating values of (*A _{2}*) within the region $(0\le \Gamma \le 0.5)$
and shows escalation behavior with escalating values of (

*A*) within the region $(0.5<\Gamma \le 1).$

_{2}3. The velocity function $g(\Gamma )$
shows escalation behavior with escalation in *φ* and show reducing behavior with escalation in (*A _{1}*),

*A*and

_{2}*(M)*.

4. The temperature function $\theta (\Gamma )$
show escalating behavior with the escalation in *φ* and (*Rd*) shows reducing behavior with the escalation in (*A _{1}*) and (Pr).

### References

- Haq RU, Nadeem S, Khan ZH, Noor NFM. 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. Ref.: https://tinyurl.com/ybhuvpbe
- Liu MS, Lin MCC, Te HI, Wang CC. Enhancement of thermal conductivity with carbon nanotube for nanofluids. Int Comm in Heat and Mass Transfer. 2015; 32: 1202-1210; Ref.: https://tinyurl.com/ycqbu7c9
- Nadeem S, Lee C. Boundary layer flow of nanofluid over an exponentially stretching surface. Nanoscale Res Lett. 2012; 7: 94. Ref.: https://tinyurl.com/ybp2rygp
- Ellahi R, Raza M, Vafai K. 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. Ref.: https://tinyurl.com/y9tn4s3b
- Nadeem S, Haq RU, Khan ZH. Numerical study of MHD boundary layer flow of a Maxwell fluid past a stretching sheet in the presence of nanopar-ticles. J TaiwanInst Chem Eng. 2014; 45: 121-126. Ref.: https://tinyurl.com/y9hsb38v
- Nadeem S, Haq RU. MHD boundary layer flow of a nano fluid past a porous shrinking sheet with thermal radiation, J Aerospace Eng. 1943-5525. Ref.: https://tinyurl.com/y97mpwdv
- Nadeem S, Haq RU. Effect of thermal radiation for megnetohydrody- namic boundary layer flow of a nanofluid past a stretching sheet with convective boundary conditions. J Comput Theor Nanosci. 2013; 11: 32-40. Ref.: https://tinyurl.com/y96y4ddt
- Choi SUS. Enhancing thermal conductivity of fluids with nanoparticles. Int mech eng congress and exposition. 1995; 66: 99-105. Ref.: https://tinyurl.com/ycq6by6d
- Buongiorno J. Convective transport in nanofluids. J Heat Transf. 2006; 128: 240-250. Ref.: https://tinyurl.com/yamfme6z
- Sheikholeslami M. Influence of magnetic field on nanofluid free convection in an open porous cavity by means of Lattice Boltzmann method. J of Molecular Liquids. 2017; 234: 364. Ref.: https://tinyurl.com/y9yxzq6f
- Sheikholeslami M. Magnetic field influence on nanofluid thermal radiation in a cavity with tilted elliptic inner cylinder. J of Molecular Liquids. 2017; 229: 137-147. Ref.: https://tinyurl.com/y7jrhyy3
- Sheikholeslami M. Magnetohydrodynamic nanofluid forced convection in a porous lid driven cubic cavity using Lattice Boltzmann method. J of Molecular Liquids. 2017; 231: 555-565. Ref.: https://tinyurl.com/y86kyh8f
- Sheikholeslami M. Numerical simulation of magnetic nanofluid natural convection in porous media. Physics Letters A. 2017; 381: 494-503. Ref.: https://tinyurl.com/ybllfpc2
- Sheikholeslami M. Numerical study of heat transfer enhancement in a pipe filled with porous media by axisymmetric TLB model based on GPU. Eur Phys J Plus. 2014; 129: 248.
- Sheikholeslami M. CVFEM for magnetic nanofluid convective heat transfer in a porous curved enclosure. Eur Phys J Plus. 2016; 131: 413. Ref.: https://tinyurl.com/y99bt4sh
- Sheikholeslami M, Rokni HB. Simulation of nanofluid heat transfer in presence of magnetic field: A review. Int J of Heat and Mass Transfer. 2017; 115: 1203-1233. Ref.: https://tinyurl.com/ybahhr9g
- Sheikholeslami M, Rokni HB. Free convection of CuOeH2O nanofluid in a curved porous enclosure using mesoscopic approach. Int J of Heat and Mass Transfer. 2017; 42: 15393. Ref.: https://tinyurl.com/y8jetqnk
- Sheikholeslami M. Houman BR. Numerical simulation for impact of Coulomb force on nanofluid heat transfer in a porous enclosure in presence of thermal radiation. International Journal of Heat and Mass Transfer. 2018; 118: 823-831. Ref.: https://tinyurl.com/y99vplln
- Sheikholeslami M, Ganji DD, Javed MY, Ellahi R. Effect of thermal radiation on magnetohydrodynamics nanofluid flow and heat transfer by means of two phase model. Journal of Magnetism and Magnetic Materials. 2015; 374: 36-43. Ref.: https://tinyurl.com/y8kd9kv6
- Elias MM, Miqdad M, Mahbubul IM, Saidur R, Kamalisarvestani M, et al. Effect of nanoparticle shape on the heat transfer and thermodynamic performance of a shell and tube heat exchanger. Int Commun Heat Mass Transfer. 2013; 44: 93-99. Ref.: https://tinyurl.com/y9owqk5n
- Alfven H. Existence of electromagnetic-hydrodynamic waves. Nature. 1942; 150: 405-406. Ref.: https://tinyurl.com/y9x9jb2p
- Shah Z, Gul T, Khan AM, Ali I, Islam S. Effects of hall current on steady three dimensional non-newtonian nanofluid in a rotating frame with brownian motion and thermophoresis effects. J Eng Technol. 2017; 6: 280-296. Ref.: https://tinyurl.com/yb28z25t
- Shah Z, Islam S, Gul T, Bonyah E, Khan MA. The electrical MHD and hall current impact on micropolar nanofluid flow between rotating parallel plates. Results Phys. 2018: 9: 1201-1214. Ref.: https://tinyurl.com/ybelycjt
- Shah Z, Islam S, Gul T, Bonyah E, Khan MA. Three dimensional third grade nanofluid flow in a rotating system between parallel plates with Brownian motion and thermophoresis effects. Results Phys. 2018; 10: 36-45. Ref.: https://tinyurl.com/yc24emj2
- Casson NA. A flow equation for pigment oil suspension of printing ink type. In: Mill CC (ed) Rheology of dispersed system. Pergamon Press. Oxford. 1959.
- Mehmood Z, Mehmood R, Iqbal Z. Numerical Investigation of Micropolar Casson Fluid over a Stretching Sheet with Internal Heating. Commun Theor Phys. 2017; 67: 443-448. Ref.: https://tinyurl.com/y7xbr4rv
- Iqbal Z, Mehmood R, Ehtsham Azhar, Mehmood Z. Impact of inclined magnetic field on micropolar Casson fluid using Keller box algorithm. Eur Phys J Plus. 2017; 132: 175. Ref.: https://tinyurl.com/yarhf5cl
- Megahe AM. Effect of slip velocity on Casson thin film flow and heat transfer due to unsteady stretching sheet in presence of variable heat flux and viscous dissipation. Open Journal of Fluid Dynamic. 2016; 6: 303-320. Ref.: https://tinyurl.com/y85wxkph
- Abolbashari MH, Freidoonimehr N, Rashidi MM. Analytical modeling of entropy generation for Casson nano-fluid flow induced by a stretching surface. Advanced Powder Technology. 2016; 2: 542-552. Ref.: https://tinyurl.com/y7fxzbae
- Attia HA, Ahmed MES. Transient MHD couette flow of a Casson fluid between parallel plates with heat transfer. Ital J Pure Appl Math. 2010; 27:
- Mukhopadhyay S. Casson fluid flow and heat transfer over a nonlinearly stretching surface. Chin Phys B. 2013; 27: 074701-074705. Ref.: https://tinyurl.com/yc6purfo
- Qasim M, Noreen S. Heat transfer in the boundary layer flow of a Casson fluid over a permeable shrinking sheet with viscous dissipation. Eur Phys J Plus. 2014;
**7:**129-137. Ref.: https://tinyurl.com/y9gm47zv - JC. Maxwell. Electricity and magnetism. 3rd ed. Clarendon, Oxford. 1904.
- Casson NA. Flow Equation for Pigment Oil Suspensions of the Printing Ink Type, Rheology of Disperse Systems. Pergamon Press: New York. NY, USA. 1959; 84-104. Ref.: https://tinyurl.com/ycb3o3z4
- Hussanan A, Salleh MZ, Tahar RM, Khan I. Unsteady boundary layer flow and heat transfer of a Casson fluid past an oscillating vertical plate with Newtonian heating. PLoS ONE. 2014; 9: e108763. Ref.: https://tinyurl.com/yaep7gfk
- Walicka A, Falicki J. Reynolds number effects in the flow of an electro rheological fluid of Casson type between fixed surface of revolution. Appl Math Comput. 2015; 250: 639-649. Ref.: https://tinyurl.com/y8yqa48z
- Hussain ST, Rizwan-ul-Haq, Khan ZH, Nadeem S. Water driven flow of carbon nanotubes in a rotating channel. Journal of Molecular Liquids. 2016; 214: 136-144. Ref.: https://tinyurl.com/ydbw4xpm
- Sheikholeslami M, Hatami D, Ganji DD. Nanofluid flow and heat transfer in a rotating system in the presence of a magnetic field. Journal of Molecular Liquids. 2014; 190: 112-120. Ref.: https://tinyurl.com/y8r7d73a
- Ishaq M, Ali G, Shah, Z, Islam S, Muhammad S. 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. Ref.: https://tinyurl.com/yb39tfy7
- Nasir S, Islam S, Gul T, Shah Z, Khan MA, et al. Three‑dimensional rotating flow of MHD single wall carbon nanotubes over a stretching sheet in presence of thermal radiation. Applied Nano science. 2018; 8: 1361-1378. Ref.: https://tinyurl.com/y9bpe8v9
- Hammed H, Haneef M, Shah Z, Islam S, Khan W, et al. The Combined Magneto hydrodynamic and electric field effect on an unsteady Maxwell nanofluid Flow over a Stretching Surface under the Influence of Variable Heat and Thermal Radiation. Appl Sci. 2018; 8: 160. Ref.: https://tinyurl.com/ycfnfqpn
- Ali A, M Sulaiman, S Islam, Z Shah, E Bonyah. Three-dimensional magnetohydrodynamic (MHD) flow of Maxwell nanofluid containing gyrotactic micro-organisms with heat source/sink. AIP Advances. 2018; 8: 085303. Ref.: https://tinyurl.com/yatjo8op