Swirling flow of fluid containing (SiO 2 ) and (MoS 2 ) nanoparticles analyze via Cattaneo-Christov theory

Abstract

Investigation of heat transport mechanism in swirling flow of viscous fluid containing silicon dioxide $\left(Si{O}_2\right)$ and molybdenum disulfide $\left(Mo{S}_2\right)$ nanoparticles is performed. The flow is engendered due to stretchable rotating cylinder which immersed in infinite fluid. The boundary layer assumption is applied to simplify the governing equations of the problem. The theory of Cattaneo-Christov for thermal energy transportation is employed in the present phenomenon under the heat and mass constraints. The flow is also influenced by Lorentz force. The results for flow field, temperature, and concentration field are produced by employing the bvp4c numerical technique to the similar differential equations. According to the observations, it is noted that in the presence of Lorentz force the reduction in velocity field of the nanofluid occurs. The thermal and solutal relaxation phenomena also declines the energy transport in nanofluid flow. The outcomes are validated through the comparison with previous published studies.

Introduction

Researchers have been interested in the fluid flow characteristics and heat transfer aspects of a stretching/shrinking surfaces for a few years. There is no doubt that the theory of fluid flow and heat transfer at a stagnation point may be applied to a wide range of industrial processes, such as plastic extrusion, wire coating, and so on. Sakiadis ¹ proposed the concept of a boundary layer formed by a moving flow on a rigid sheet with constant velocity from a slit through a fluid at rest. Based on the outcomes, with increasing distance from the slit, the thickness of the boundary layer increases. The concept of boundary on a moving surface was investigated by Crane. ² Following Crane, ² Carragher and Crane ³ examined the thermal approach to this problem. McLeod and Rajagopal ⁴ address the existence and uniqueness of these stretched flow solutions. Wang ⁵ has provided many references on stretching surfaces in a viscous fluid. Miklavčič and Wang ⁶ examined the shrinking case of a viscous fluid initially, followed by Fang ⁷ and Tie-Gang et al. ⁸ As Goldstein ⁹ discusses, a shrinking sheet flow is fundamentally a backward boundary-layer flow that exhibits physical features distinct from stretching sheet flows. Very recently, the numerical simulation for the nanofluid flow due to a stretching/shrinking sheet was discussed by Ferdows et al. ¹⁰ and Mousavi et al. ¹¹

There are several applications of swirling flow due to rotating surfaces in industry such as paper drying, plastic manufacturing and wire coating. Many scientists put their effort to investigate the mechanism of swirling flow of viscous fluid with heat transport in the recent decade. Swirl motion of Newtonian fluid was investigated by Fang and Yao. ¹² In this article authors established the suitable flow similarities for the conversion of governing three-dimensional flow equations. Naumov et al. ¹³ investigated the thin layer in interface of two immiscible fluids flow of swirl motion with the help of numerical technique. Hayat et al. ¹⁴ analyzed the 3D rotating flow of nanofluid in porous medium by using Darcy-Forchheimer model. Flow between two rotating disks of electro-viscous fluid in the presence of porous medium with Hall effect was studied by Krishna et al. ¹⁵ In this study the variation in pressure gradient was taken as time dependent. MHD flow of non-Newtonian fluid with help of Darcy-Forchheimer model was scrutinized by Ali et al. ¹⁶ by using well-known finite element method. The studies in this area can be found here (see Refs^17–22).

Production of best quality products in polymers industry is mainly depends upon the flow as well as heat transport rate. As the conventional liquid have low thermal conduction so many engineering procedures were adopted to overcomes this problem. Such as fins (extended surfaces), surface roughness, use of interfacial forces etc. Scientists discover nanofluids (a new category of fluid) which comprises the suspended in nano-sized particles $(1-100nm)$ in low conductivity fluid. The concept of nanofluids initially was given by Choi ²³ . A lot of focus in recent years put on the thermal characteristics of nanofluids and their flow properties. Mahanthesh et al. ²⁴ examined numerical solutions for three dimensional electrically conducting fluid flow across a non-linear stretching surface. The flow caused by a rotating disk with varying thickness of the the water-based nanofluid with silver nanoparticles was explored by Hayat et al. ²⁵ Zheng et al. ²⁶ presented the heat transfer and fluid flow properties of different nanofluids using a curved plate heat exchanger in solar energy systems. Maleki et al. ²⁷ studied the thermal conductivity of nanofluids containing $ZnO$ particles by three strategies. The similar studies can be found in Refs.^28–33

The above comprehensive literature review prove that there is space for the investigation of nanofluid swirling flow. The focus of this article is to explore the thermal transportation characteristics of flow of viscous fluid containing $Si{O}_2$ / $Mo{S}_2$ particles. The flow produced via a stretched and rotating cylinder with constant surface temperature. The realistic approach of Catteneo-Christov theory is employed to predict the heat and mass transport mechanism in the nanofluid. The source/sink and chemical reaction parameters are also utilized as heat and mass constraint. The MATLAB built-in solver bvp4c is used to achieve numerical results of the governing problem. The results of the velocity, temperature, and concentration field are illustrated graphically for various values of significant parameters.

Mathematical formulation

Consider an elastic cylinder of radius $R_1$ which is immersed in an infinite viscous fluid which contains $Si{O}_2$ and $Mo{S}_2$ solid nanoparticles. The axisymmetric 3D flow produce due the constant rotation $E$ and stretching with rate $a$ of the cylinder. The flow field is assumed as $V=[u,v,w]$ in the $\left(z,\varphi ,r\right)$ directions. The magnetic field $B=\left(0,0,B_0\right)$ is applied normal to the cylinder axis (see Figure 1). Moreover, the heat source/sink and chemical reaction are also considered. Furthermore, it is supposed that the temperature at the cylinder surface $T_w$ and the ambient fluid temperature is $T$ _∞ are constant and $T_w>T_{\infty }$ . In view of conservation laws and Cattaneo-Christov theory for heat flux the governing equation for the present problem are ¹² :

$$\frac{\partial u}{\partial z}+\frac{w}{r}+\frac{\partial w}{\partial r}=0,$$

(1)

$$u\frac{\partial u}{\partial z}+w\frac{\partial u}{\partial r}=\frac{{\mu }_{nf}}{{\rho }_{nf}}\left[\frac{{\partial }^{2}u}{\partial r^2}+\frac{1}{r}\frac{\partial u}{\partial r}\right]-\frac{{\sigma }_{nf}}{{\rho }_{nf}}{B}_0^{2}u,$$

(2)

$$u\frac{\partial v}{\partial z}+w\frac{\partial v}{\partial r}+\frac{wv}{r}=\frac{{\mu }_{nf}}{{\rho }_{nf}}\left[\frac{{\partial }^{2}v}{\partial r^2}-\frac{v}{r^2}+\frac{1}{r}\frac{\partial v}{\partial r}\right]-\frac{{\sigma }_{nf}}{{\rho }_{nf}}{B}_0^{2}v,$$

(3)

$$\begin{array}{c}u\frac{\partial T}{\partial z}+w\frac{\partial T}{\partial r}+{\lambda }_t\left[\begin{array}{l}{u}^2\frac{{\partial }^{2}T}{\partial z^2}+2uw\frac{{\partial }^{2}T}{\partial z\partial r}\\ +w^2\frac{{\partial }^{2}T}{\partial r^2}+u\frac{\partial u}{\partial z}\frac{\partial T}{\partial z}\\ +w\frac{\partial u}{\partial r}\frac{\partial T}{\partial z}+u\frac{\partial w}{\partial z}\frac{\partial T}{\partial r}\\ +w\frac{\partial w}{\partial r}\frac{\partial T}{\partial r}\end{array}\right]\\ \begin{array}{l}=\frac{k_{nf}}{{\left(\rho C_p\right)}_{nf}}\left[\frac{{\partial }^{2}T}{\partial r^2}+\frac{1}{r}\frac{\partial T}{\partial r}\right]+\frac{Q_0}{{\left(\rho C_p\right)}_{nf}}\left(T-T_{\infty }\right)\\ +{\lambda }_t\frac{Q_0}{{\left(\rho C_p\right)}_{nf}}\left(u\frac{\partial T}{\partial z}+w\frac{\partial T}{\partial z}\right),\end{array}\end{array}\}$$

(4)

$$\begin{array}{c}u\frac{\partial C}{\partial z}+w\frac{\partial C}{\partial r}+{\lambda }_c\left[\begin{array}{l}{u}^2\frac{{\partial }^{2}C}{\partial z^2}+2uw\frac{{\partial }^{2}C}{\partial z\partial r}\\ +w^2\frac{{\partial }^{2}C}{\partial r^2}+u\frac{\partial u}{\partial z}\frac{\partial C}{\partial z}\\ +w\frac{\partial u}{\partial r}\frac{\partial C}{\partial z}+u\frac{\partial w}{\partial z}\frac{\partial C}{\partial r}\\ +w\frac{\partial w}{\partial r}\frac{\partial C}{\partial r}\end{array}\right]\\ \begin{array}{l}=D_B\left[\frac{{\partial }^{2}C}{\partial r^2}+\frac{1}{r}\frac{\partial C}{\partial r}\right]-k_r\left(C-C_{\infty }\right)\\ -{\lambda }_c{k}_r\left(u\frac{\partial C}{\partial z}+w\frac{\partial C}{\partial z}\right),\end{array}\end{array}\}$$

(5)

The following BC’s are utilized to analyze the above problem.

$$u=2az,v=E,w=0,T=T_{\infty },C=C_{\infty }atr=R_1,$$

(6)

$$u\to 0,v\to 0,T\to T_{\infty },C\to C_{\infty }asr\to \infty .$$

(7)

Figure 1.

Flow configuration.

Here, $(u,v,w)$ are the axial, azimuthal, and radial velocity components, respectively, ${\mu }_{nf}$ the dynamic viscosity of nanofluid, ${\rho }_{nf}$ density of nanofluid, ${\sigma }_{nf}$ electric conductivity of nanofluid, $k_{nf}$ thermal conductivity of nanofluid, ${\left(\rho c_p\right)}_{nf}$ specific heat of nanofluid, $B_0$ magnetic field strength, $\left({\lambda }_t,{\lambda }_c\right)$ the relaxation time for heat and mass transport, respectively, $Q_0$ heat generation source, $D_B$ the Brownian diffusion coefficient, and $(T,C),$ the liquid temperature and concentration, respectively.

Presenting the following flow similarities and with the help of Table 1

$$\begin{array}{c}u=2az{f}^{\prime }(\eta ),v=Eg\left(\eta \right),w=\frac{-a{R}_{1}f\left(\eta \right)}{{\eta }^{\frac{1}{2}}},\\ \theta (\eta )=\frac{T-T_{\infty }}{T_w-T_{\infty }},\varphi (\eta )=\frac{C-C_{\infty }}{C_w-C_{\infty }},\eta =\frac{r^2}{R_1^2}\end{array}\},$$

(8)

the equation (1) is satisfied automatically and equations (2−5) yield

$$\frac{A1}{A2}\left({{f^{\prime }}^{\prime }}^{\prime }+\frac{{f^{\prime }}^{\prime }}{\eta }\right)+\frac{Re}{\eta }f{f^{\prime }}^{\prime }\frac{Re}{\eta }{f}^{\prime }-\frac{A5MRe}{A2}\frac{f^{\prime }}{2\eta }=0,$$

(9)

$$\begin{array}{l}\frac{A1}{A2}\left(2{\eta }^2{g}^{″}+2\eta g^{\prime }-\frac{g}{2}\right)+2Re\eta f{g}^{\prime }\\ +Refg-\frac{A5}{A4}MRe\eta =0,\end{array}$$

(10)

$$\begin{array}{l}\left(\eta {\theta }^{″}+{\theta }^{\prime }\right)+\frac{RePr}{A4}\left(\frac{\delta }{2}\theta +A3f{\theta }^{\prime }\right)\\ -2{\beta }_t\frac{A3}{A4}RePr\left({\theta }^{″}{f}^2+f{f}^{\prime }{\theta }^{\prime }\right)=0,\end{array}$$

(11)

$$\begin{array}{l}\left(\eta {\varphi }^{″}+{\varphi }^{\prime }\right)+RePrLe\left(f{\varphi }^{\prime }-\frac{\gamma }{2}\varphi \right)\\ -2{\beta }_{c}RePrLe\left(f^2{\varphi }^{″}+f{f}^{\prime }{\varphi }^{\prime }\right)=0,\end{array}$$

(12)

with transformed boundary conditions

$$\begin{array}{c}f(1)=0,f^{\prime }(1)=1,g\left(1\right)=1,\theta \left(1\right)=\varphi (1)=1\\ f^{\prime }(\infty )=0,g\left(\infty \right)=0,\theta (\infty )=0,\varphi (\infty )=0\end{array}\},$$

(13)

Where $Re\left(=\frac{a{R}_1^2}{2v}\right)$ is the Reynolds number, $M\left(=\frac{{\delta }_f{B}_0^2}{{\rho }_{f}a}\right)$ the Magnetic field parameter, $Pr\left(=\frac{\nu }{{\alpha }_1}\right)$ the Prandtl number, $Le\left(=\frac{{\alpha }_1}{D_B}\right)$ the Lewis number, $\delta \left(=\frac{Q_0}{a{\left(\rho c_p\right)}_f}\right)$ the source/sink parameter, $\gamma \left(=\frac{k_r}{a}\right)$ the chemical reaction parameter, ${\beta }_t\left(={\lambda }_{t}a\right)$ and ${\beta }_c\left(={\lambda }_{c}a\right)$ the thermal and mass relaxation parameters. Moreover the $A1,A2,A3,A4,A5$ can be stated as follows

$$\begin{array}{l}A1=\left(\frac{{\mu }_{nf}}{{\mu }_f}\right),A2=\left(\frac{{\rho }_{nf}}{{\rho }_f}\right),A3=\left(\frac{{\left(\rho c_p\right)}_{nf}}{{\left(\rho c_p\right)}_f}\right),\\ A4=\left(\frac{k_{nf}}{k_f}\right),A5=\left(\frac{{\sigma }_{nf}}{{\sigma }_f}\right).\end{array}$$

(14)

A further transformation $\eta =e^x$ is employed to accelerate convergence of the solution. Thus equations (8) − (11) become

$$\begin{array}{l}\frac{A1}{A2}\left(f_{xxx}-2f_{xx}+f_x\right)-Re\left(f_x^2+f{f}_x-f{f}_{xx}\right)\\ -\frac{A5}{A4}\frac{MRe}{2\left(f_x{e}^x\right)}=0,\end{array}$$

(15)

$$2\frac{A1}{A2}{g}_{xx}-\frac{g}{2}+Re\left(2f{g}_x+fg\right)-\frac{A5}{A4}MRe\left(g{e}^x\right)=0,$$

(16)

$$\begin{array}{l}{\theta }_{xx}{e}^x+\frac{RePr}{A4}\left(e^{-2x}\frac{\delta }{2}\theta +A3f{\theta }_x{e}^x\right)\\ -2{\beta }_{t}RePr\frac{A3}{A4}\left(f^2{\theta }_{xx}-f^2{\theta }_x+2f{f}_x{\theta }_x\right)=0\end{array}$$

(17)

,

$$\begin{array}{l}{\varphi }_{xx}{e}^x+RePrLe\left(f{\varphi }_x{e}^x-\frac{\gamma \varphi }{2}\right)\\ -2{\beta }_{c}RePrLe\left(f^2{\varphi }_{xx}-f^2{\varphi }_x+f{f}_x{\varphi }_x\right)=0,\end{array}$$

(18)

with BCs

$$\begin{array}{c}f\left(0\right)=0,f_x\left(0\right)=1,g\left(0\right)=1,\theta \left(0\right)=1,\varphi \left(0\right)=1\\ li{m}_{x\to \infty }{e}^{-x}{f}_x=0,g\left(\infty \right)=0,\theta \left(\infty \right)=0,\varphi \left(\infty \right)=0\end{array}\},$$

(19)

Table 1.

Mathematical relations for thermophysical properties of nanofluids.

Properties	Nanofluid
Density	${\rho }_{nf}={\rho }_f\left(1-\varphi +\varphi \frac{{\rho }_f}{{\rho }_s}\right)$
Viscosity	${\mu }_{nf}=\frac{{\mu }_f}{{\left(1-\phi \right)}^{2.5}}$
Heat Capacity	${\left(\rho C_p\right)}_{nf}={\left(\rho C_p\right)}_f\left(1-\varphi +\varphi \frac{{\left(\rho C_p\right)}_s}{{\left(\rho C_p\right)}_f}\right)$
Thermal Conductivity	$\frac{k_{nf}}{k_f}=\frac{k_s+2k_f-2\varphi \left(k_f-k_s\right)}{k_s+2k_f+2\varphi \left(k_f-k_s\right)}$
Electric Coductivity	$\frac{{\sigma }_{nf}}{{\sigma }_f}=1+\frac{3\left(\sigma -1\right)\varphi }{\left(\sigma +2\right)-\left(\sigma -1\right)\varphi }$

Physical reasoning of the outcomes

In this portion of research, the graphical results which are computed numerically describe the impact of physical parameters such as Reynolds number Re, magnetic field parameter M, volume fraction coefficient $\varphi$ , Prandtl number $Pr$ , Lewis number $Le$ , sourcing parameter $\delta$ , chemical reaction parameter $\gamma$ , and thermal and mass relaxation parameters $\left({\beta }_t,{\beta }_c\right)$ for both nanoflu-ids $\left(Si{O}_2/\mathrm{water}\mathrm{and}Mo{S}_2/\mathrm{water}\right)$ with the provided thermophysical properties in Table 2 on flow field and temperature and concentration distributions. When we go through the numerical computation, we specify the values of relevant parameters as $Re=1.0,M=2.0,\varphi =0.025,\mathrm{Pr}=6.2,Le=1,\gamma =\delta =1.5,{\beta }_t=0.1,{\beta }_c=0.1$ . Figure 2(a) to (c) represents the results for velocity field, temperature, and solutal distributions for various values of Reynolds number $Re.$ It is noted that the axial velocity $f^{\prime }\left(\eta \right)$ , swirling velocity $g\left(\eta \right)$ , and radial velocity $\frac{f}{{\eta }^{\frac{1}{2}}}$ are significantly decrease when Re boosted. Physically, the opposing flow agent in the form of inertial force increase when the Reynolds number enhance with fixed viscosity of the fluid. Moreover, it is also observed that there is algebraic decay in flow field for higher Re and the disturbance in the fluid occurs very far away to the free stream from surface of cylinder when Re have small magnitude. As the transport mechanism of thermal and solutal energy in the present problem is exist due to the forced convection in the system and the fluid motion produced due cylinder rotation with stretching. The outcomes for flow field with higher Re prove that fluid motion limited to the surface of cylinder in this case. Thus, the forced convection also decreases for higher Re and diminishing temperature, concentration profiles are observed in Figure 3(a) and (b). The objective of Figure 4 is to demonstrate the effects of Prandtl number Pr on temperature. These results tell us that the thermal diffusivity of the fluid decreases as Pr values rise. Figure 5(a) to (c) demonstrate the effect of the magnetic field parameter M on the velocity profile. The velocity profile is shown to declines as the magnetic field parameter is raised. Because the transverse magnetic field in the flow system of electrically conducting fluid generates a drag-like Lorentz force. This force tends to oppose the motion of nanofluid. Figure 6(a) to (c) depicts the decreasing behavior of axial velocity, swirling velocity, and radial velocity as nanoparticles volume fractions is enhanced. Physically, the higher concentration of nanoparticles in the viscous fluid produce the resistance to the fluid motion as fluid particles require more kinetic energy to flow. Thus, in results flow field decreases. The thermal and solutal conduction enhancement of base fluid is booted due to higher concentration of solid particles in it. These results are shown in Figure 7(a) and (b). The volume fraction parameter is an important one that plays a role in improving the thermal characteristics of fluids. As a conclusion, it may deduce that in many industrial processes, temperature can be adjusted by changing the nanoparticles volume fraction. Figure 8(a) and (b) displays the impact of Lewis number $Le$ and chemical reaction parameter $\gamma$ on concentration field. It is noted the mass diffusivity of the fluid decrease with more values of $Le$ and a chemical reaction the fluid declines the mass transport. The temperature and concentration patterns are displayed versus the thermal relaxation time parameter ${\beta }_t$ and the mass relaxation time parameter ${\beta }_c$ , accordingly, in Figure 9(a) and (b). The temperature and concentration fields decline as the thermal relaxation time and thermal relaxation mass parameter values improve. In Cattaneo-Christov theory, higher values of the relaxation time parameter control the immediate propagation of heat waves in a given medium. As a consequence, the fluid with increasing relaxation parameter values required more time to carry the thermal and solutal energy in the base fluid. Thus, the temperature and concentration fields drop as a result. The results of the present study validated through the comparison Table 3. The values of thermal gradient at the cylinder surface are given in Table 4 as well.

Table 2.

Thermophysical characteristics of nanoparticles Silicon dioxide $\left(Si{O}_2\right)$ and Molybdenum disulfide $\left(Mo{S}_2\right)$ and base fluid water $\left(H_{2}O\right)$ .

	$\rho /kg{m}^{-3}$	$c_p/k{g}^{-1}{K}^{-1}$	$\sigma /S{m}^{-1}$	$k/Wm{K}^{-1}$
Water	997.1	4179	0.163	5.5x10−6
$\left(Mo{S}_2\right)$	5060	397.746	34.5	17.9x10−6
$\left(Si{O}_2\right)$	2650	730	1.5	10−12

Figure 2.

(a–c) Impact of Reynolds number Re on $f^{\prime }\left(\eta \right),g\left(\eta \right)$ and $f/{\eta }^{\frac{1}{2}}$

Figure 3.

(a and b) Impact of Reynolds number Re on $\theta \left(\eta \right)$ and $\varphi \left(\eta \right)$ .

Figure 4.

Impact of Prandtl number Pr on θ $\varphi \left(\eta \right)$ .

Figure 5.

(a–c) Impact of magnetic field parameter M on $f^{\prime }\left(\eta \right),g\left(\eta \right)$ , and $f/{\eta }^{\frac{1}{2}}$ .

Figure 6.

(a–c) Impact of volume fraction coefficient $\varphi$ on $f^{\prime }\left(\eta \right),g\left(\eta \right)$ , and $f/{\eta }^{\frac{1}{2}}$ .

Figure 7.

(a and b) Impact of volume fraction coefficient $\varphi$ on $\theta \left(\eta \right)$ and $\varphi \left(\eta \right)$ .

Figure 8.

(a and b) Impact of Lewis number $Le$ and chemical reaction parameter $\gamma$ on $\varphi \left(\eta \right)$ .

Figure 9.

(a and b) Impact of thermal ${\beta }_t$ and mass relaxation parameter ${\beta }_c$ on $\theta \left(\eta \right)$ and $\varphi \left(\eta \right)$ respectively.

Table 3.

The comparison table for numerical values of $f^{″}\left(1\right)$ and $g^{\prime }\left(1\right)$ .

Re		0.1	0.2	0.5	1.0	2.0	5.0	10.0
${f^{\prime }}^{\prime }\left(1\right)$	Fang and Yao 12	−0.48180	−0.61748	−0.8822	−1.17775	−1.59389	−2.41743	−3.34446
${f^{\prime }}^{\prime }\left(1\right)$	Present	−0.48252	−0.6161	−0.884	−1.17821	−1.59785	−2.41817	−3.34445
$g^{\prime }\left(1\right)$	Fang and Yao 12	−0.51019	−0.52605	−0.58488	−0.68772	−0.87263	−1.29788	−1.81006
$g^{\prime }\left(1\right)$	Present	−0.5102	−0.52595	−0.5859	−0.68796	−0.87266	−1.29759	−1.81019

Table 4.

Numerical values of $-{\theta }^{\prime }(1)$ for various $Re,Pr,\delta$ , and ${\beta }_t.$

$Re$	$Pr$	$\delta$	${\beta }_t$	$-{\theta }^{\prime }(1)$
0.5	6.2	2	0.1	2.242906
0.7				2.382541
0.9				2.552958
	0.5			1.751631
	0.7			1.747774
	0.9			1.744803
		0.5		2.842302
		0.7		2.613992
		0.9		2.242906
			0.5	2.341122
			0.7	1.740708
			0.9	0 368323

Concluding remarks

Thermal energy transport in the nanofluid flow caused by stretchable rotating cylinder investigated in this article. The water base fluid contains the two different nanoparticles ( $Si{O}_2/Mo{S}_2$ ). The mechanism of thermal and solutal energy transport predicted via well-known Cattaneo-Christov theory. Thermal conduction controlled with heat source and Lorentz force in the system. Governing similar equations solved numerically. Some major points of the entire analysis are as follows:

• The flow confined to the surface of cylinder when Reynolds number $Re$ was increased.

• Forced convection declined in case of higher inertial force.

• Flow of nanofluid feels resistance when Lorentz force and solid particles concentration were boosted.

• The temperature of the fluid rises as $\varphi$ increased.

• Thermal and solutal relaxation phenomenon limited the energy transportation in the nanofluid flow

Future scope

The present work can be extend for the rheology of non-Newtonian fluids with heat transport. Moreover, entropy generation analysis can also be performed in the present investigated phenomenon.