Shape Effect of Nanosize Particles on Magnetohydrodynamic Nanofluid Flow and Heat Transfer over a Stretching Sheet with Entropy Generation

Abstract

Magnetohydrodynamic nanofluid technologies are emerging in several areas including pharmacology, medicine and lubrication (smart tribology). The present study discusses the heat transfer and entropy generation of magnetohydrodynamic (MHD) Ag-water nanofluid flow over a stretching sheet with the effect of nanoparticles shape. Three different geometries of nanoparticles—sphere, blade and lamina—are considered. The problem is modeled in the form of momentum, energy and entropy equations. The homotopy analysis method (HAM) is used to find the analytical solution of momentum, energy and entropy equations. The variations of velocity profile, temperature profile, Nusselt number and entropy generation with the influences of physical parameters are discussed in graphical form. The results show that the performance of lamina-shaped nanoparticles is better in temperature distribution, heat transfer and enhancement of the entropy generation.

1. Introduction

The entropy generation is related to thermo-dynamic irreversibility; entropy generation is common in every category of heat transfer processes. The entropy generation technique is introduced by Bejan [1,2]. The quantity of irreversibility equal to that of entropy can be measured during a process. Entropy analysis is an influential means to find the efficiency of a system or process. Limited studies on the entropy generation of turbulent nanofluid flow have been analyzed in the microchannel, stretching sheet and channels; the entropy generation in the microchannel and channel is studied more closely although the studies described can apply to nanofluid, fluid turbulent, or laminar flow [3]. Some researchers have worked on entropy generation by considering different types of geometry. Ko and Cheng [4] numerically studied the entropy generation in a channel wavy wall. According to the results of Ko and Cheng, amplitude of the wavy wall is increasing function of the entropy generation. Mahian et al. [5] presented the entropy generation of water-TiO₂ and ethylene glycol-Al₂O₃ nanofluids between two circular cylinders. The obtained results displayed that the entropy generation has inverse relation with the volume fraction of nanofluid. Tshehla and Makinde [6] presented the entropy generation of steady flow between two concentric pipes. They demonstrated that entropy generation is directly relational to Brinkman number. Govindaraju et al. [7] analyzed the impact of slip magnetohydrodynamic nanofluid on entropy generation in a stretching sheet. The authors found that the entropy generation reduces in the flow system in the presence of slip parameter and volume fraction. Govindaraju et al. [1] presented the entropy generation study of an incompressible magnetohydrodynamic of viscous nanofluid flow over a stretching sheet by considering various categories of nanoparticles in the water-based fluid. They concluded that entropy generation has a direct relation with Hartmann number, magnetic parameter and dimensionless group parameter. Berrehal and Maougal [8] presented entropy generation analysis for nanofluid flow over a wedge with convective boundary condition and thermal radiation. The authors found that the entropy generation can be decreased with deducing convection through boundary and with increase radiation parameter. Malvandi et al. [9] discussed analytical solution of entropy generation for nanofluids over a flat plate. They demonstrated that the entropy generation depends on Prandtl number, Reynold number, Eckert number and volume fraction. Acharya et al. [10] presented the entropy generation and heat transfer in a regenerative cooling channel of a rocket engine. Govindaraju et al. [2] studied analytically the entropy generation of nanofluid flow over a stretching sheet with a uniform heat source-sink and inclined magnetic field. The results exposed that entropy generation has a direct relation with Ec and ϕ, while entropy generation decreases with the increase in A.

Magnetohydrodynamics attracted the attention of researchers due to natural phenomena, i.e., geophysics, astrophysics to several applications in engineering such as electromagnetic casting, liquid metal, plasma confinement and so on [11]. Nadeem et al. [12] discussed two-dimensional boundary layer flow over a stretching sheet with the effect of magnetohydrodynamics. The results exposed that when increasing the value of the Prandtl number and magnetic parameters, an opposite behavior is seen in the Sherwood number and Nusselt number. Rudraiah et al. [13] numerically examined the impact of magnetic field on natural convection. The authors noted that heat transfer rate decreases with the effect of the magnetic field. Abbas et al. [14] numerically discussed the effects of magnetohydrodynamics over a stretching continuous sheet in a rotating fluid.

The concept of a nanofluid was introduced by Choi [15] who examined notable results with several possibilities of usage. Nanofluids are the latest class of nanotechnology concerning nanoparticles dispersed in base fluid. The heat transfer fluid became a significant interest in research due to various applications. The heat transfer rate of a nanofluid is greater due to greater thermal conductivities as compared to the base fluids [16]. There are several papers on water-Ag nanofluid flow by many researchers. Atashafrooz [17] studied the impact of water-Ag nanofluid flow over an inclined step by using numerical technique. Upreti et al. [18] examined the water-Ag nanofluid over a flat porous plate with influences of injection/suction, heat absorption/generation and viscous-ohmic dissipation. Suleman et al. [19] discussed the water-Ag nanofluid flow in a stretching cylinder with the effects of homogenous-heterogeneous reaction and Newtonian heating.

The geometry effect of nanoparticles is very important to change the thermal conductivity of the nanofluid [20]. Several researchers have worked on Ag-water nanofluid flow. To the best of our knowledge, there is no such research on the shape effect of nanoparticles on magnetohydrodynamic Ag-water nanofluid flow over a stretching sheet with entropy generation. The motivation of the present paper is to investigate the shape effect of Ag nanoparticles on magnetohydrodynamic Ag-water nanofluid flow and heat transfer with entropy generation. Finally, graphs of velocity, temperature, Nusselt number and entropy generation are plotted and all their aspects are discussed.

This study is organized as follows: in the first section, we construct the mathematical formulation for the proposed model. The numerical solution of proposed model is obtained by homotopy analysis method in Section 3. The entropy generation is discussed in Section 4. The computational results are reported in Section 5. Finally, conclusions are described in Section 6.

2. Mathematical Formulation

This section considers entropy analysis for a two-dimensional steady nanofluid laminar flow over a linear semi-infinite stretching sheet. Furthermore, it is assumed that $B_0$ magnetic field is imposed normally to the stretching sheet. It is also assumed that the nanofluid is water based containing various shapes of Ag nanoparticles. The thermo-physical properties of Ag nanoparticles and water are presented in

Table 1. Thermo-physical properties.

Physical Properties	Ag	Water
$\rho$ (kg/m3)	10,500	997.1
$Cp$ (J/kg K)	235	4179
k (W/m K)	429	0.60

. The values of nanoparticles shape related parameter is presented in

Table 2. Shape factor (m) parameters.

Shapes	Blades	Sphere	Lamina
m	8.6	3	16.1576

. The governing equations are

$$\frac{\partial \left(u\right)}{\partial x} + \frac{\partial \left(v\right)}{\partial y} = 0,$$

(1)

$$u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = \frac{{\mu }_{n_{f }}}{{\rho }_{n_f}}\frac{{\partial }^{2}u}{\partial y^2} + \frac{{\sigma }_{n_f}{B}_0{}^2{u}^2}{{\rho }_{n_f}},$$

(2)

$$u\frac{\partial T}{\partial x} + v\frac{\partial T}{\partial y} = \frac{k_{n_f}}{{\left(\rho Cp\right)}_{nf}}\frac{{\partial }^{2}T}{\partial y^2}.$$

(3)

The boundary value conditions are

$$u = U_{\mathrm{w}}\left(x\right) = ax,v=0,T =T_w, \mathrm{at} y=0, u = v = 0, T = T_{\infty }, \mathrm{at} y \to \infty .$$

(4)

The nondimensionalize variables are defined as

$$x = \frac{x}{\sqrt{\frac{{\nu }_f}{a}}}, y = \frac{y}{\sqrt{\frac{{\nu }_f}{a}}}, U = \frac{u}{\sqrt{a{\nu }_f}}, V = \frac{v}{\sqrt{a{\nu }_f}}, \theta \left(\eta \right) = \frac{T-T_{\infty }}{T_f-T_{\infty }}$$

(5)

After substituting Equation (5) into Equations (1)–(3) we obtain

$$f^{‴}+B1\left\{B2\left(f{f}^{″}-f^{\prime \prime 2}\right)-M{f}^{\prime }\right\} = 0,$$

(6)

$${\theta }^{″}+Pr\frac{B3}{B4}f{\theta }^{\prime } = 0.$$

(7)

The corresponding boundary value conditions are

$$f\left(0\right) = f^{\prime }\left(0\right) = 1, f^{\prime }\left(\infty \right) = 0, \theta \left(0\right) = 1,\theta \left(\infty \right) = 0,$$

(8)

where,

$$\begin{gathered} B1={(1-\varphi )}^{2\cdot 5}, \\ B2=\left(1-\varphi \right)+\frac{{\rho }_s}{{\rho }_f}\varphi , \\ B3=\left(1-\varphi \right)+\frac{{\left(\rho Cp\right)}_s}{{\left(\rho Cp\right)}_f}\varphi , \\ B4=\frac{k_{n_f}}{k_f}=\frac{\left[k_s+\left(m-1\right)k_f\right]-\left(m-1\right)\varphi \left(k_f-k_s\right)}{\left[k_S+\left(m-1\right)k_f\right]+\varphi \left(k_f-k_s\right)}. \end{gathered}$$

The important physical quantity of interest, the Nusselt number (N_u) is defined as [1]

$$N_{u_x}, = \frac{x q_w\left(x\right)}{k_f\left[ T_f- T_{\infty }\right]},$$

(9)

where $q_w\left(x\right)$ is given by

$$q_w\left(x\right) = -k_{n_f}{\left(\frac{\partial T}{\partial y}\right)}_{y = 0}.$$

(10)

Substituting Equation (9) into Equation (10), we obtain

$$\frac{N_{u_x}}{\sqrt{R{e}_x}} = -\frac{k_{n_f}}{k_f}{\theta }^{\prime }\left(0\right).$$

(11)

3. Solution via Homotopy Analysis Method

For the solution of the homotopy analysis method, we choose the following auxiliary linear operators.

$${\mathcal{L}}_f=f^{‴}-f^{\prime }, {\mathcal{L}}_{\theta }=f^{″}-f,$$

(12)

which satisfied the following properties

$${\mathcal{L}}_f\left[Z_1+Z_2{\mathrm{e}}^{\eta }+Z_3{\mathrm{e}}^{-\eta }\right]=0, {\mathcal{L}}_{\theta }\left[Z_4{\mathrm{e}}^{\eta }+Z_5{\mathrm{e}}^{-\eta }\right]=0,$$

(13)

where $Z_1$, $Z_2$, $Z_3$, $Z_4$ and $Z_5$ are arbitrary constants.

It is assumed that p ∈ [0, 1] denotes an embedding parameter and $\hslash$f, $\hslash$θ ≠ 0 are convergence control parameters.

The zeroth order deformations are

$$\left(1-p\right){\mathcal{L}}_f\left[\widehat{f}\left(\eta ,p\right)-f_0\left(\eta \right)\right]=p{\hslash }_f{\mathcal{N}}_f\left[\widehat{f}\left(\eta ,p\right),{\widehat{\theta }}_0\left(\eta ,p\right)\right],$$

(14)

$$\begin{gathered} \widehat{f}\left(0,p\right)=0, {\widehat{f}}^{\prime }\left(0,p\right)=1, \widehat{f}\left(\infty ,p\right)=0. \\ \left(1-p\right){\mathcal{L}}_{\theta }\left[\widehat{\theta }\left(\eta ,p\right)-{\theta }_0\left(\eta \right)\right]=p{\hslash }_{\theta }{\mathcal{N}}_{\theta }\left[ {\widehat{\theta }}_0\left(\eta ,p\right),\widehat{f}\left(\eta ,p\right)\right], \\ \widehat{\theta }\left(0,p\right)=1, \widehat{\theta }\left(\infty ,p\right)=0, \end{gathered}$$

(15)

and

$$\begin{array}{ll}{N}_f\left[ \widehat{f}\left(\eta ,p\right),\widehat{\theta }\left(\eta ,p\right)\right]& \\ & =\frac{{\partial }^3 \widehat{f}\left(\eta ,p\right)}{\partial {\eta }^3}+B1\{B2\left(\widehat{f}\left(\eta ,p\right)\frac{{\partial }^2 \widehat{f}\left(\eta ,p\right)}{\partial {\eta }^2}-\frac{{\partial }^2 \widehat{f}\left(\eta ,p\right)}{\partial {\eta }^2}\frac{{\partial }^2 \widehat{f}\left(\eta ,p\right)}{\partial {\eta }^2}\right)\\ & -M\frac{\partial \widehat{f}\left(\eta ,p\right)}{\partial \eta }\}=0,\end{array}$$

(16)

$$N_{\theta }\left[ \widehat{\theta }\left(\eta ,p\right), \widehat{ f}\left(\eta ,p\right)\right]= \frac{{\partial }^2\widehat{\theta }\left(\eta ,p\right)}{\partial {\eta }^2}+Pr\frac{B3}{B4}\widehat{f}\left(\eta ,p\right)\frac{\partial \widehat{\theta }\left(\eta ,p\right)}{\partial \eta }=0.$$

(17)

Due to Taylor’s series with respect to p, we have

$$f\left(\eta ,p\right)=f_0\left(\eta \right)+{\displaystyle \sum }_{m=1}^{\infty }{f}_m(\eta )p^m, \theta \left(\eta ,p\right)={\theta }_0\left(\eta \right)+{\displaystyle \sum }_{m=1}^{\infty }{\theta }_m(\eta )p^m,$$

(18)

where,

$$f_m\left(\eta \right)={\left.\frac{1}{m!}\frac{{\partial }^{m}f\left(\eta ,p\right)}{\partial {\eta }^m}\right|}_{p=0}, \theta { }_m\left(\eta \right)={\left.\frac{1}{m!}\frac{{\partial }^m \theta \left(\eta ,p\right)}{\partial {\eta }^m}\right|}_{p=0},$$

(19)

and thus, higher order deformation problems are

$$\begin{gathered} {\mathcal{L}}_f \left[f_m\left(\eta \right)-{\chi }_m{f}_{m-1}\left(\eta \right)\right]={\hslash }_f{\mathcal{R}}_f^m\left(f_{m-1}\left(\eta \right),{\theta }_{m-1}\left(\eta \right)\right), \\ {f}_m\left(0\right)=0,f_m^{\prime }\left(0\right)=0,f_m^{\prime }\left(\infty \right)=0, \end{gathered}$$

(20)

$$\begin{gathered} \mathcal{L}{L}_f \left[{\theta }_m\left(\eta \right)-{\chi }_m{\theta }_{m-1}\left(\eta \right)\right]={\hslash }_{\theta }{\mathcal{R}}_{\theta }^m\left({\theta }_{m-1}\left(\eta \right),f_{m-1}\left(\eta \right)\right), \\ {\theta }_m\left(0\right)=0,{\theta }_m\left(\infty \right)=0. \end{gathered}$$

(21)

Here,

$${\chi }_m=\left\{\begin{array}{c}0 m \le 1 \\ 1 m>1\end{array} \right.$$

(22)

and

$${\mathcal{R}}_f^m{f}_m\left(\eta \right)=f_{m-1}^{\prime \prime \prime \prime }\left(\eta \right)+B1\{B2{\displaystyle \sum }_{z=0}^{m-1}{f}_z{f}_{m-1-z}^{\prime \prime }-\frac{1}{2}{\displaystyle \sum }_{z=0}^{m-1}{f}_z^{\prime }{f}_{m-1-z}^{\prime })+f_z{f}_{m-1}^{\prime }\left(\eta \right)\},$$

(23)

$${\mathcal{R}}_f^m{\theta }_m\left(\eta \right)={\theta }_{m-1}^{\prime \prime }+Pr\frac{B3}{B4}{\displaystyle \sum }_{z=0}^{m-1}{f}_z{\theta }_{m-1-z}^{\prime }.$$

(24)

4. Entropy Generation

The entropy generation in the influence of $B_0$ magnetic field is defined as:

$$S_a=\frac{k_{n_f}}{ T_{\infty }}\left[{\left(\frac{\partial T}{\partial x}\right)}^2+{\left(\frac{\partial T}{\partial y}\right)}^2\right]+\frac{{\mu }_{n_{f }}}{ T_{\infty }}{\left(\frac{\partial u}{\partial y}\right)}^2+\frac{{\sigma }_{n_f}{B}_0{ }^2{u}^2}{ T_{\infty }}.$$

(25)

The first term of Equation (25) expresses the entropy generation due to heat transfer, the middle term expresses the viscous dissipation and the last term denotes the entropy generation due to magnetic field effect. The entropy generation rate in dimensionless method is described as the ratio of volumetric entropy ($S_a$) to characteristic entropy generation rate ($S_b$), the characteristic entropy generation rate is defined as:

$$S_b=\frac{k_{n_f}{\left(\nabla T\right)}^2}{x^2 T_{\infty }^2}.$$

(26)

The entropy number is defined as [1]

$$N=\frac{ S_a}{ S_b},$$

(27)

$$N={\theta }^{\prime 2}\left(\eta \right)+\frac{B}{\Omega }{f}^{\prime \prime 2}\left(\eta \right)+\frac{B {\left( H_a\right)}^2}{ \Omega R_e}{f}^{\prime 2}\left(\eta \right).$$

(28)

5. Results and Discussion

In the present section, the effects of several parameters, such as magnetic parameter (M) and solid volume fraction ($\varphi )$, on velocity $f^{\prime }\left(\eta \right),$temperature $\theta \left(\eta \right)$ and Nusselt number ($Nu$) are displayed in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7. The effect of nanoparticles has also been plotted. The variation of velocity and temperature profiles with the influence of $\varphi$ is displayed in Figure 1 and Figure 2. In Figure 1 and Figure 2, a decrement in velocity and temperature of nanofluid increase with an increase in the$\varphi$can be seen. It is noticed that thermal conductivity of nanofluid has direct relation with solid volume fraction, thus the thermal boundary layer thickness is increased. It is also noted from Figure 2, performance of lamina-shaped nanoparticles on temperature distribution is better than blade and sphere shapes nanoparticles. The velocity and temperature profiles with influences (M) are depicted in Figure 3 and Figure 4. It is observed that velocity and temperature of the nanofluid have an inverse and direct relationship with M, respectively. Physically, the Lorentz force produced due to M opposes the motion of the nanofluid, also Lorentz force heating temperature equation will provide an extra heat source to the boundary layer of the nanofluid. It is observed from Figure 4 that the performance of lamina-shaped nanoparticles on temperature distribution is better than blade- and sphere-shaped nanoparticles. Figure 5 displays the effect of nanoparticles in temperature profile. Figure 5 shows that the performance of lamina-shaped nanoparticles on temperature distribution is better than blade- and sphere-shaped nanoparticles. This is because the nanoparticles in the form of a lamina have low viscosity and thermal conductivity compared to nanoparticles in the form of a blade and a sphere. The impacts of magnetic parameter (M) and solid volume fraction ($\varphi )$ on Nusselt number ($Nu)$ are plotted in Figure 6 and Figure 7. From Figure 6 and Figure 7, it is noted that lamina-shaped nanoparticles have a higher rate of heat transfer as compared to blade- and sphere-shaped nanoparticles.

Table 3. Results comparison for Nusselt number $\left(Nu\right)$ when $M = \varphi = 0.$

$\mathit{N}\mathit{u}$	Present Results	Govindaraju et al. [1]	Wang Results [21]
0.7	0.0654	0.065562	0.0656
0.2	0.1690	0.169089	0.1691
70	6.4682	6.462200	6.4622

shows that excellent agreement has been found in the comparison of results with Govindaraju et al. [1] and Wang Results [21].

6. Conclusions

This paper explores the heat transfer and entropy generation of magnetohydrodynamic (MHD) Ag-water nanofluid flow over a stretching sheet with the effect of nanoparticles shapes. The homotopy analysis method is used for analytical solutions of velocity profile, temperature profile, Nusselt number and entropy generation. The fixed value of Prandtl number $Pr$ = 6.2 is used. A agree ment has been obtained Lamina-shaped nanoparticles show dramatic role in the disturbance of temperature profile, heat transfer and increase in entropy generation, while the performance of sphere-shaped nanoparticle is lowest in disturbance of temperature profile, heat transfer and increase in entropy generation.