Thermally Radiative Rotating Magneto-Nanofluid Flow over an Exponential Sheet with Heat Generation and Viscous Dissipation: A Comparative Study
Sagheer M., Bilal M. *, Hussain S., Ahmed R. N.
Department of Mathematics, Capital University of Science and Technology, Islamabad, Pakistan

 

† Corresponding author. E-mail: m.bilal@cust.edu.pk

Abstract
Abstract

This article examines a mathematical model to analyze the rotating flow of three-dimensional water based nanofluid over a convectively heated exponentially stretching sheet in the presence of transverse magnetic field with additional effects of thermal radiation, Joule heating and viscous dissipation. Silver (Ag), copper (Cu), copper oxide (CuO), aluminum oxide (Al2O3) and titanium dioxide (TiO2) have been taken under consideration as the nanoparticles and water (H2O) as the base fluid. Using suitable similarity transformations, the governing partial differential equations (PDEs) of the modeled problem are transformed to the ordinary differential equations (ODEs). These ODEs are then solved numerically by applying the shooting method. For the particular situation, the results are compared with the available literature. The effects of different nanoparticles on the temperature distribution are also discussed graphically and numerically. It is witnessed that the skin friction coefficient is maximum for silver based nanofluid. Also, the velocity profile is found to diminish for the increasing values of the magnetic parameter.

PACS: ;47.15.Cb;;47.50.Cd;
Nomenclature
1 Introduction

In the present fast growing and developing computer age, the transportation, communication, heavy mechanical industries, electronics industries and house hold appliances, all are running by some mechanical and electronic devices. Almost in all such devices, according to the requirements of devices, a system of cooling or heating is built-in, by which a fluid flows through or around the device to prevent these devices from overheating or cooling down from certain temperature threshold. To meet the human requirements and demand of the market, it is essential that these devices work round the clock. To keep the devices at a constant temperature, the heat dissipated must be equal to the heat generated. The conventional fluids with low thermal conductivity do not meet the temperature requirements of many mechanical and electronics devices, which results in poor performance of these devices and reduces their efficiency and working age. Therefore, it is imperative to improve the thermal conductivity of the conventional fluids.

The conventional fluids used for the transfer of heat energy were first time replaced by the nanofluids by Choi[1] followed by many researchers. A nanofluid is a mixture of nanoparticles in a conventional heat transfer fluid. The nanoparticles (1–100 nm) in size are usually metals, metallic oxides, nanofibers, etc. Choi[1] experimentally found that the nanoparticles when added to the base fluids, considerably improve the thermal conductivity of the base fluid. Magyari and Keller[2] focused on the heat and mass transfer analysis in the boundary layer flow due to an exponentially continuous stretching sheet. Eastman et al.[3] observed that the thermal conductivity of pure ethylene glycol is much increased when copper nanoparticles are added to it. Li et al.[4] investigated the MHD nanofluid flow in a thin film through unsteady stretching sheet with additional effects of thermal radiation, heat generation, Brownian motion, and thermophoresis. They used the MATLAB built-in bvp4c solver to solve their ODEs. It is found form their investigation that temperature and nanoparticle concentration have opposite behavior for the thermophoresis parameter. Nadeem et al.[5] analyzed the influence of nanoparticles on the two-dimensional flow of Maxwell nanofluid over a stretching sheet for the heat and mass transfer effects. By applying the boundary-layer approximation, they also incorporated the effects of MHD and elasticity parameter. Sheikholeslami et al.[6] presented an analysis focusing on the unsteady squeezing flow of electrically conducting nanofluid using the homotopy perturbation method. Two phase simulation model for the nanofluid is considered along with the magnetohydrodynamics effects. They concluded that the Nusselt number is a decreasing function of the squeezing parameter. Another useful contribution of Sheikholeslami and Ganji,[7] is a review work addressing both the single and the double phase models for the nanofluids. They describe briefly the various attempts of different scientists on heat transfer of convective nanofluids. It is further analyzed that while increasing the Reynolds number and Rayleigh number, the rate of heat transfer is increased. Sheikholeslami et al.[8] discussed the thermal radiation on MHD free convection of Al2O3-water nanofluid. Chopkar et al.[9] studied the effect of the size of the nanoparticles on the thermal conductivity of nanofluid and found that the thermal conductivity decreases by increasing the size of the particles. In last half decade, many articles related to the nanofluid dynamics are published in literature.[1024]

A reasonable number of applications emphasizing the role of steady and unsteady rotating flows may be found in chemical and geophysical fluid mechanics. These all are of applied nature like in the thermal power generating systems, food processing, the skins of high speed air crafts and in rotor stator systems. The pioneering work highlighting the rotating flow was done by Wang.[25] Takhar et al.[26] discussed the effects of magnetohydrodynamic in a rotating flow. They concluded that the skin friction along the x-axis increases for the higher values of the magnetic parameter and has a reverse relation for the y-axis skin friction coefficient. Zaimi et al.[27] applied the numerical technique to examine the rotating flow of viscoelastic fluid. Turkyilmazoglu[28] applied the spectral numerical integration method for the problem related to the shrinking rotating disk with the effect of magnetohydrodynamic. Some recent attempts emphasizing the rotating flow can be found in Refs. [2931].

During the study of nanofluid, the thermal radiative properties of Newtonian and non-Newtonian fluids for the heat transfer phenomenon have got much attention. Because of the insertion of the nanoparticles in the base fluid, the thermal properties are enhanced which resultantly rises the temperature of the nanofluids and for the higher temperature differences, the effects of the thermal radiation cannot be neglected. The operating systems performing the energy conversion at high temperature show a comparable effect of the thermal radiation. In other engineering and chemical processes such as solar water technology, fossil fuel combustion, astrophysical flows, hypersonic flights, gas turbines, space vehicles, nuclear reactors etc., the effects of thermal radiations are quite phenomenal. Many researchers have considered the influence of thermal radiation on the boundary layer flow of Newtonian and non-Newtonian fluids. Mushtaq et al.[32] considered the nonlinear thermal radiation in the two-dimensional stagnation point flow with additional effects of Joule heating and viscous dissipation over a convectively heated surface. They deduced that both the temperature and its gradient are increasing functions of thermal radiation parameter. Pourmehran et al.[33] numerically investigated the MHD boundary layer flow of nanofluid through convectively heated vertical stretching sheet. During the study, the influence of thermal radiation and buoyancy effects got special attention. Pourmehran et al. considered three different types of base fluid i.e., pure water, ethylene glycol 30% and ethylene glycol 50% while the four types of nanoparticles i.e., copper, silver, alumina, and titanium oxide.

Motivated by the above mentioned literature, the primary objective of the present study is to examine the effects of thermal radiation and viscous dissipation on heat transfer flow over a bi-directional convectively heated exponentially stretching sheet in the presence of transverse magnetic field and volumetric rate of heat generation. Five different nanoparticles (silver, copper, copper oxide, titanium oxide, alumina) are assumed to be suspended in the pure water. A detailed comparative study of these nanofluids for the flow and heat transfer is presented and discussed graphically and numerically. Boundary layer approximations are used to govern the partial differential equations, which are then transformed to the ordinary differential equations with the help of transformations. The modeled problem is solved numerically by the shooting method using Runga-Kutta integration scheme of order 4. Effects of emerging parameters on velocity and temperature profiles are discussed in detail. Nusselt number and skin friction coefficient are also calculated. In the limiting case, the results are verified by reproducing the results of previously published article[3435]

2 Mathematical Formulation

A laminar, incompressible and steady water-based electrically conducted nanofluid flow over an exponentially bidirectional stretching sheet is considered. Sheet temperature Tf is controlled via convection by considering hot fluid below it. The temperature faraway from the surface where its difference is negligible is known as the ambient temperature and is denoted by T. The fluid is assumed to rotate with angular velocity Ω = Ω0 ex/L along z-axis having the coriolis effect. A transverse variable magnetic field B(x) = B0 ex/2L is applied along z-axis with the assumption of small Reynolds number, ignoring the induced magnetic field. The fluid has internal volumetric rate of heat generation Q0. In the x direction, the velocity of the sheet is taken as Uw = U ex/L as shown in Fig. 1.

Fig. 1 Geometry of the Problem.

Further the effects of the thermal radiation, Joule heating and viscous dissipation are considered in the formulation of energy equation. For the nanofluid model, the Tiwari and Das model.[36] has been utilized. Applying the boundary layer by incorporating the Boussinesq approximations, the conservation equations of mass, momentum and energy in the mathematical form can be expressed as

The corresponding boundary conditions for the velocity and temperature expressions are

where, hf is the convective heat transfer coefficient and μnf is the dynamic viscosity of nanofluid proposed by Brinkman[37] and these are given by the relations

The effective heat capacity (ρcp)nf, the effective density ρnf and the thermal diffusivity αnf of the nanofluid are formulated as[3840]

The effective thermal conductivity modeled by Hamilton-Crosser (H-C)[41] is given by

In Table 1, the thermo-physical properties of different nanoparticles and pure water are shown. The following dimensionless variables are used to convert the system of the non-linear PDEs to the system of ODEs.

In Eq. (4), qr is a radiative heat flux which is defined as

Furthermore, the temperature difference within the flow is assumed such that T4 may be expanded in a Taylor series. Hence, expanding T4 about T and neglecting the higher order terms, we get

Using the similarity transformation defined in Eq. (9), Eq. (1) is identically satisfied while Eqs. (2)–(5) are converted into the following nonlinear ordinary differential equations:

The transformed boundary conditions are:

Different dimensionless parameters appearing in Eqs. (11)–(15) are defined as

The important quantities of interest are the skin friction coefficient and the local Nusselt number which are defined as

where τwx, τwy and qw denote the wall shear stress along x-axis, y-axis and wall heat flux respectively and are given as under:

The skin friction coefficient and the local Nusselt number, in the dimensionless form, are as follows

where Rex = U0z2/νl.

Table 1

Thermo-physical properties of H2O and nanoparticles.

.
3 Solution Methodology

An efficient numerical technique, namely the shooting method has been employed to solve the transformed ordinary differential equations along with the boundary conditions for different values of the emerging parameters. While applying the shooting method,[42] first the higher order boundary value problem is converted to a system of first order initial value problem (IVP). During the conversion, f is denoted by y1, g by y4 and θ by y6. The missing initial conditions are supposed to be ς1, ς2 and ς3. The converted first order IVP takes the following form

Fourth order Runge-Kutta method is utilized to solve this IVP. The refinement of initial guesses is carried out by the Newton’s method. Because the numerical solution cannot be computed on the unbounded domain [0, ∞), a bounded domain [0, η] has been considered, where η is an appropriate real number. After performing a number of computational experiments, η is set to 4, because there is no significant variation in the results for η > 4. The stoping criteria set for the Newton’s iterative process is

Throughout this article, ϵ is chosen as 10−6. For the validation of the MATLAB code of the shooting method, it is affectively applied to reproduce the numerical results of Javed et al.[35] and Ahmad and Mustafa.[34] The successful comparison has been presented in Table 2.

Table 2

Comparison of present results with those of Javed et al.[35] and Ahmad and Mustafa.[34]

.
4 Results and Discussions

In this section, we discuss the influence of different parameters such as nanoparticles volume fraction ϕ, rotational parameter λ, magnetic parameter M, thermal radiation parameter R, Eckert number Ec, heat generation/absorption parameter Qh on the velocity, temperature, skin-friction and Nusselt number, both graphically and numerically in the tabular form.

In Table 3, the influence of the nanoparticle volume fraction ϕ, rotational parameter λ and magnetic parameter M on the skin friction coefficient along x-axis is presented for different nanoparticles. It is observed that due to the addition of more nanoparticles in the base fluid, the skin-friction is enhanced. This enhancement is more rapid in Ag-H2O nanofluid whereas in case of Al2O3-H2O, the increase in the skin-friction is less as compared to the other nanolfuids. Quite similar behavior is noticed for the rotational parameter λ. When the magnetic field is intensified along the z-axis, the skin-friction escalates along the x-axis due to the presence of the Lorentz force. Again, Ag-H2O nanofluid has more frictional force as compared to the other nanofluids.

Table 3

Effect of ϕ, λ and M on the skin friction coefficient along x-axis when Qh = 0.1, Bi = 0.6, Pr = 6.2, n = 3.0, R = 0.2, Ec = 0.01.

.
Table 4

Numerical values of skin friction coefficient along y-axis for different values of parameters when Qh = 0.1, R = 0.2, Ec = 0.01.

.

Table 3 shows the effect of variation of nanoparticles volume fraction ϕ, rotational parameter λ, and magnetic parameter M on the skin friction coefficient along y-axis. From the table, it is highlighted that fractional force between the fluid and the solid surface in the y-direction is enhanced when the nanoparticles volume fraction is increased. However this increase in skin friction is very small. Again it is observed that Al2O3 water base nanoflluid has least increase when compared with the other nanofluids. By escalating the angular velocity of the nanoparticles, the skin friction rises along the y-axis. The effect of magnetic field M on the skin friction along the y-axis is quite similar as already shown in Table 3. By enhancing the magnetic parameter, the surface fractional force also increased.

The effect of thermal radiation parameter R, Eckert number Ec, heat generation/absorption parameter Qh and Biot number on Nusselt number is shown in Table 5. From this table, it is observed that these parameters have increasing effect on the Nusselt number. The comparison among different nanoparticles exhibits that Al2O3-H2O nanofluid possesses the highest value of Nusselt number for thermal radiation parameter. It is also perceived that Ag-H2O nanofluid has more heat transfer rate as compared to the other nanofluids for the increasing values of Eckert number and heat generation/absorption parameter Qh. Nusselt number is also enhanced for the higher values of the Biot number Bi. It happens because for the higher values of the Biot number, a stronger convection is produced which causes higher rate of change in the temperature.

Table 5

Numerical values of local Nusselt for different values of parameters when ϕ = 0.01, λ = 0.2, M = 0.3.

.

To visualize the effect of different physical parameters on the velocity f′(η) and the temperature profile θ(η), Figs. 26 are plotted. In Figs. 26, the effects of nanoparticle volume fraction ϕ, rotational parameter λ and magnetic parameter M on the velocity profile for alumina, copper oxide, titanium oxide, copper and silver based nanofluids are displayed. For all the nanofluids, it is observed that the velocity as well as the boundary layer thickness of the nanofluid decreases when the quantity of the nanoparticles in the base fluid is increased. Velocity distribution is dominant at the surface of the sheet. The effect of rotational parameter λ which is the associated with the angular velocity of the fluid, on velocity profile is displayed in Figs. 26. From these figures, it is noticed that the velocity profile and its momentum boundary layer thickness is reduced for the increasing values of λ. Hence, the rotational effects resist the fluid flow in the x-direction. For higher values of rotational parameter, the velocity becomes negative in some part of the boundary layer thickness and an interesting phenomenon of oscillatory decaying profile are also observed. The Lorentz forces, which are resistive in nature are produced when the magnetic field is applied across the fluid flow. These forces are responsible for the reduction in the fluid particle’s motion for the higher values of magnetic parameter M. Hence for all the nanofluids, the speed of the fluid decreases for the increasing values of magnetic parameter.

Fig. 2 Influence of ϕ, λ, and M on f′(η) for Al2O3-H2O nanofluid.
Fig. 3 Influence of ϕ, λ, and M on f′(η) for CuO-H2O nanofluid.
Fig. 4 Influence of ϕ, λ and M on f′(η) for TiO2-H2O nanofluid.
Fig. 5 Influence of ϕ, λ and M on f′(η) for Cu-H2O nanofluid.
Fig. 6 Influence of ϕ, λ and M on f′(η) for Ag-H2O nanofluid.

To observe the effect of the variation in the Biot number Bi, Eckert number Ec, and heat generation parameter Qh on the temperature distribution Figs. 711 are plotted. It is observed that the higher values of the Biot number escalate the temperature distribution and the thermal boundary layer thickness. The same observation is preserved for all the nanofluids. The strength of the convected heating is signified for the higher values of the Biot number which resultantly rise the temperature distribution. The temperature is enhanced when the Eckert number is increased. Eckert number appears in the energy equation because of the consideration of the viscous dissipation effects in the fluid motion. It is inversely proportional to the difference between the fluid temperature on the surface and the ambient temperature. An increase in the Eckert number means there is a slight temperature difference between the surface and the thermal boundary layer and hence the rate of heat transfer is reduced. This reduction in heat transfer rate leads to escalate the temperature of the nanofluid as shown in Figs. 711. The effect of heat generation parameter Qh is also displayed in the same figures. It is quite obvious that if heat is generated from any external or internal source the temperature of the fluid is increased. In Figs. 1216, the influence of the thermal radiation parameter R, magnetic parameter M and nanoparticle volume fraction ϕ is displayed for the temperature distribution. Higher values of thermal radiation produces more heat in the working fluid which rises the temperature and the thermal boundary layer thickness of the nanofluid as shown in these figures. By increasing the magnetic field across the fluid, the resistive forces are enhanced. Temperature is increased due to these resistive forces. Lastly, by inserting the more quantity of nanopartices in the base fluid, the thermal properties of the fluid go up and hence the temperature of the fluid is increased.

Fig. 7 Influence of Bi, Ec and Qh on θ(η) for Al2O3-H2O nanofluid.
Fig. 8 Influence of Bi, Ec and Qh on θ(η) for CuO-H2O nanofluid.
Fig. 9 Influence of Bi, Ec and Qh on θ(η) for TiO2-H2O nanofluid.
Fig. 10 Influence of Bi, Ec and Qh on θ(η) for Cu-H2O nanofluid.
Fig. 11 Influence of Bi, Ec and Qh on θ(η) for Ag-H2O nanofluid.
Fig. 12 Influence of R, M and ϕ on θ(η) for Al2O3-H2O nanofluid.
Fig. 13 Influence of R, M and ϕ on θ(η) for CuO-H2O nanofluid.
Fig. 14 Influence of R, M and ϕ on θ(η) for TiO2-H2O nanofluid.
Fig. 15 Influence of R, M and ϕ on θ(η) for Cu-H2O nanofluid.
Fig. 16 Influence of R, M and ϕ on θ(η) for Ag-H2O nanofluid.
5 Concluding Remarks

This article encompasses the three-dimensional MHD rotating flow of electrically conducting nanofluid over an exponentially stretching sheet. The effect of heat generation, viscous dissipation and thermal radiation for five different nanoparticles is analyzed graphically and numerically. The main findings of the investigation are as follows.

Al2O3-H2O nanofluid has more capacity to transfer heat as compared to the other discussed nanofluids when the thermal radiation is enhanced.

The skin friction coefficient is maximum for Ag-H2O nanofluid.

An increase in the Eckert number Ec and the heat generation parameter Qh reduces the Nusselt number. This reduction in the heat transfer rate is much lower for Ag-H2O nanofluid.

The velocity profile diminishes for increasing values of the magnetic parameter M.

Ag-H2O and Cu-H2O nanofluids have greater values of the Nusselt number as compared to Al2O3-H2O and TiO2-H2O

Reference
[1] Choi S. U. S. ASME Int. Mech. Engr. Cong. Exp. 78 1995 99
[2] Magyari E. Keller B. J. Phys. D 32 1999 577
[3] Eastman J. A. Choi S. U. S. Li S. et al. Appl. Phys. Let. 78 2001 718
[4] Li J. Liu L. Zheng L. Mohsin B. B. J. Taiwan Inst. Chem. Engr. 67 2016 226
[5] Nadeem S. Haq R. U. Khan Z. H. J. Taiwan Inst. Chem. Engr. 45 2014 121
[6] Sheikholeslami M. Hatami M. Domairry G. J. Taiwan Inst. Chem. Engr. 46 2015 43
[7] Sheikholeslami M. Ganji D. D. J. Taiwan Inst. Chem. Engr. 65 2016 43
[8] Sheikholeslami M. Hayat T. Alsaedi A. Int. J. Heat Mass Transfer 96 2016 513
[9] Chopkar M. Sudarshan S. Das P. K. Manna I. Metall. Matter. Trans. A 39 2009 1535
[10] Ramzan M. Bilal M. PLoS ONE 10 2015 e0124929
[11] Ramzan M. Bilal M. J. Mol. Liq. 215 2016 212
[12] Bilal M. Sagheer M. Hussain S. Alex. Engr. J. org/10.1016/j.aej.2017.03.039 2017
[13] Sheikholeslami M. Shehzad S. A. Int. J. Heat Mass Transfer 109 2017 82
[14] Bilal M. Hussain S. Sagheer M. Bull. Po. Acad. Sci. Tech. Sci. 65 2017 383
[15] Sheikholeslami M. Shehzad S. A. Int. J. Heat Mass Transfer 113 2017 796
[16] Sheikholeslami M. Hayat T. Alsaedi A. Int. J. Heat Mass Transfer 108 2017 1870
[17] Sheikholeslami M. Bhatti M. M. Int. J. Heat Mass Transfer 109 2017 115
[18] Sheikholeslami M. Rokni H. B. Int. J. Heat Mass Transfer 107 2017 288
[19] Sheikholeslami M. Bhatti M. M. Int. J. Heat Mass Transfer 111 2017 1039
[20] Sheikholeslami M. Sadoughi M. K. Int. J. Heat Mass Transfer 116 2018 909
[21] Sheikholeslami M. Bhatti M. M. Int. J. Heat Mass Transfer 111 2017 1039
[22] Sheikholeslami M. Phys. B 516 2017 55
[23] Sheikholeslami M. Eur. Phys. J. Plus. 132 2017 55
[24] Sheikholeslami M. Rokni H. B. Int. J. Heat Mass Transfer 115 2017 1203
[25] Wang C. Y. Zeitschrift für angewandte Math-ematik und Physik ZAMP 39 1988 177
[26] Takhar H. S. Chamkha A. J. Nath G. Int. J. Therm. Sci. 42 2003 23
[27] Zaimi K. Ishak A. Pop I. Appl. Math. Mech. 34 2013 945
[28] Turkyilmazoglu M. Comp. & Fluids 90 2014 51
[29] Mushtaq A. Mustafa M. Hayat T. Alsaedi A. Comp. & Fluids 27 2016 2223
[30] Mushtaq A. Mustafa M. Hayat T. Alsaedi A. Adv. Powder Technol. 27 2017 2223
[31] Rehman A. U. Mehmood R. Nadeem S. Appl. Therm. Engr. 112 2017 832
[32] Mushtaq A. Mustafa M. Hayat T. Alsaedi A. J. Taiwan Inst. Chem. Engr. 45 2014 1176
[33] Pourmehran O. Rahimi-Gorji M. Ganji D. D. J. Taiwan Inst. Chem. Engr. 65 2016 162
[34] Ahmed R. Mustafa M. J. Mol. Liq. 220 2016 635
[35] Javed T. Sajid M. Abbas Z. Ali N. Int. J. Num. Meth. Heat & Fluid Flow 21 2011 903
[36] Tiwari R. Das S. Int. J. Heat Mass Trans. 50 2007 2002
[37] Brinkman H. C. J. Chem. Phys. 20 1952 571
[38] Sheikholeslami M. J. Mol. Liq. 234 Supplement C 2017 364
[39] Sheikholeslami Mohsen Magnetic Int. J. Hydrogen Energy 42 2017 19611
[40] Sheikholeslami M. Zeeshan A. Comput. Methods Appl. Mech. Eng. 320 Supplement C 2017 68
[41] Hamilton R. L. Crosser O. K. Ind. Engr. Chem. Fund. 1 1962 187
[42] Na T. Y. Computational Methods in Engineering Boundary Value Problem Acad. Press New York 1979 71 76