Generalized Hybrid Nanofluid Model with the Application of Fully Developed Mixed Convection Flow in a Vertical Microchannel?

Hang Xu(徐航)and Qiang Sun(孫強)

1State Key Laboratory of Ocean Engineering,School of Naval Architecture,Ocean and Civil Engineering,Shanghai Jiao Tong University,Shanghai 200240,China

2Particulate Fluids Processing Centre,Department of Chemical Engineering,The University of Melbourne,Parkville,VIC 3010,Australia

AbstractThe fully developed mixed convection hybrid nanofluid flow in a vertical microchannel is examined in detail.The simplified hybrid model that omits the nonlinear terms due to the interaction of different nanoparticle volumetric fractions is derived and compared with the existing one.The generalized model describing hybrid nanofluid suspended with multiple kinds of solid particles is suggested.The argument that the corresponding nanofluid solutions obtained by the homogenous model can be recovered from the results of the regular problems through simple arithmetic operations is checked.Solutions in similarity form for this flow problem are formulated by means of a set of similarity variables.The effects of various parameters on important physical quantities are analyzed and discussed.

Key words:hybrid nanofluid,vertical channel,mixed convection,slippery effect

Nomenclature

?

1 Introduction

Many recent studies revealed that nanofluids have better heat transfer capability than regular fluids.Therefore it is possible to replace traditional heat transfer fluids by nanofluids in the design of various heat transfer systems such as cooling systems,heat regenerators,and heat exchangers.Choi[1]noticed that,by suspending nanometer-sized metallic particles in conventional heat transfer fluids,the resulting nanofluids hold higher thermal conductivities than those of currently used ones.Xuan and Li[2]attributed the heat transfer enhancement of nanofluids to the increase of thermal conductivity of the nanofluid.Eastman et al.[3]found that the particle shape has stronger effects on effective nanofluid thermal conductivity than particle size or particle thermal conductivity.Wen and Ding[4]speculated that possible reasons for the heat transfer enhancement of nanofluids are due to the migration of nanoparticles and the resulting disturbance of the boundary layer.Buongiorno[5]concluded that the Brownian diffusion and thermophoresis are dominant factors for heat enhancement within the boundary layer owing to the effect of the temperature gradient and thermophoresis.Other classic researches on nanofluids have been experimentally done by Pak and Cho,[6]Xie et al.,[7]Williams et al.[8]and numerically done by Tiwari and Das,[9]Oztop and Abu-Nada,[10]Raza et al.,[11]Khan,[12]Sajid et al.,[13]Sheikholeslami and Sadoughi,[14]Sheremet and Pop,[15]Kefayati[16]and Alsabery et al.[17]

Some researchers made attempts to investigate the characters of nanofluids containing different kinds of nanoparticles.Suresh et al.[18]found that both thermal conductivity and viscosity of hybrid nanofluids increase with the nanoparticle volume concentration while the viscosity increase is substantially higher than the increase in thermal conductivity for an Al2O3-Cu hybrid nanofluid.The behaviours of hybrid nanofluids then were examined in detail by different researchers such as Esfe et al.,[19]Rostamian et al.[20]and Ebrahimi and Saghravani.[21]From modelling point of view,Devi and Devi[22]proposed a mathematical model to investigate the effects of Lorentz force over a three-dimensional stretching surface subject to Newtonian heating.Tayebi and Chamkha[23?24]considered natural convection in an annulus between two confocal elliptic cylinders and eccentric horizontal cylinders filled with a Cu-Al2O3/water hybrid nanofluid.Huminic and Huminic[25]examined the influence of hybrid nanofluids on the performances of elliptical tube.Rostami et al.[26]considered mixed convective stagnation-point flow of an aqueous silica Calumina hybrid nanofluid.

This paper intends to analyze a fully developed mixed convection hybrid nanofluid flow in a vertical microchannel by means of a generalized hybrid nanofluid model.We are to simplify the Devi and Devi’s model[22]by omitting the nonlinear terms due to the interaction of different nanoparticle volumetric fractions.Then we extend Devi and Devi’s model[22]to the case that the hybrid nanofluids contain various kinds of nanoparticles.The argument by Magayari[27]that corresponding nanofluid results can be recovered from the solutions of already solved regular problems by simple arithmetic operations is then checked.The similarity solutions for this microchannel fl ow and heat transfer of a hybrid nanofluid are formulated by means of a set of similarity variables.It should be mentioned at this end,that the studies on hybrid nanofluids are still very new at this stage.There is no conclusive idea on how nanoparticles act on fluid flow and heat transfer.Complementary studies are urgently needed to understand the heat transfer characteristics of hybrid nanofluids,especially for those in suspension of multiple kinds of small particles.

2 Generalized Hybrid Nanofluid Model

In experimental and numerical studies on nanofluids’behaviours,it is a common practice to model their physical quantities by using simplified mathematical relations between the corresponding ones of base fluid and solid particles,as presented by many researchers such as Vajravelu et al.and[28]Devi and Devi.[22]Several experiments have been carried out to confirm the validity of such expressions for dilute nanofluids in suspension of one single kind of solid particles[6]and two types of mixed solid particles.[18]Devi and Devi[22]suggested a group of correlations for physical quantities of hybrid nanofluids.In their approach,they took the fluid containing one kind of nanoparitcles as the base fluid and the other kind of nanoparticles as the individual particles.The correlations of viscosity and thermal conductivity matched the experimental results given by Suresh et al.[18]

In Devi and Devi’s approach,[22]there are nonlinear terms due to the interaction of two kinds of different nanoparticles.However,in dilute solutions in which the nanoparticle volumetric fractions are usually small,the effects of these nonlinear terms may not be significant. Therefore,we reasonably neglect the nonlinear terms in Devi and Devi’s model.[22]Our simplified model of hybrid nanofluid,as well as the classic nanofluid model and Devi and Devi’s model[22]are listed in Table 1 in which Type I denotes the traditional nanofluid model(nanofluid in suspension one kind of small particles),Types II and III,respectively,denote Devi and Devi’s hybrid nanofluid model[22]and our simplified hybrid nanofluid model(nanofluid in suspension two different kind of small particles).

Devi and Devi’s approach[22]used the recurrence formulae to represent the viscosity,density,specific heat and thermal conductivity of the hybrid nanofluid corresponding to the n-th kinds of nanoparticles as

and

Table 1 Models of nanofluid and hybrid nanofluid.

By neglecting the nonlinear terms in above correlations,we obtain

Note that we keep the recurrence formula(5)for the thermal conductivity khnfsince the interactions between different particles can hardly be expressed using the Maxwell equation.Also,M=3 is chosen throughout this work that means that the particle shape is spherical.

3 Mathematical Description

Consider a mixed convection flow of a hybrid nanofluid in a constant porosity medium between two parallel vertical infinite walls separated by a distance of 2H.As shown in Fig.1,the Cartesian coordinate system(x,y)is chosen with the x-axis being along the walls and the y-axis being perpendicular to the walls.The temperatures on both walls are assumed to vary linearly along the height that are prescribed as T1+ax and T2+ax on the left and the right walls,respectively.Since the hydrodynamically fl ow is fully-developed,the velocity along the wall is only a function of y.Invoking the Boussinesq approximation,the governing equations are written as

subject to the boundary conditions

Fig.1 Physical sketch.

It is easy to see from Eq.(12)that ?2p/?x?y=0.This indicates that?p/?x is a constant and that all terms on the right-hand side of Eq.(11)are only dependent on y.Based on this fact,we define the following variables

where Ur=gβfKaH/νfis a reference velocity.

Substituting Eq.(16)into Eqs.(10),(11),and(13),the continuity equation(10)is automatically satisfied,and the rest of equations are reduced to

with the boundary conditions

where

4 Comparison Analysis Between Models

We calculate the coefficients ε1,ε2,and ε3in Eqs.(17)and(18)by Type II and Type III hybrid nanofluid models listed in Table 1.The data regarding to the basic thermophysical properties of the base fluid and nanoparticles are given in Table 2 where the thermophysical properties of water is chosen at 25C.

Table 2 Thermophysical properties of fluid and nanoparticles.

Substituting the quantities in Table 2 into different hybrid nanofluid models,the values of ε1,ε2,and ε3can be obtained,as shown in Table 3.It can be seen from the table,when the solution is dilute,namely,the nanoparticle volumetric fractions are small,the difference between the values of ε1,ε2and ε3obtained by both models is imperceptible.Recalling the experimental results and modelling tests by Pak et al.[6]and Suresh et al.,[18]it is clear that the computational range of nanoparticle volumetric fractions listed in Table 3 is acceptable and widely used.As a result,we can infer that our simplified hybrid nanofluid model is meaningful for prediction of nanofluids’behaviours.Further verification shows that,for the cases that two kinds of nanoparticles coexist,each nanoparticle volumetric fraction needs to be small to keep the solution remaining dilute.

Magyari[27]once found that,without consideration of velocity-slip effects,the governing equations of homogeneous nanofluid models can be reduced via elementary scaling transformations to the corresponding equations of the regular fluids.Thus he concluded that the corresponding nanofluid results can be recovered from the solutions of already solved problems with regular Newtonian fluids by simple arithmetic operations.

Here,we would like to check if this applicability is valid on the hybrid nanofluid flow problems.We introduce the following scaling transformations:

Substituting Eq.(21)into Eqs.(17)and(18),we obtain

To keep Eqs.(22)and(23)invariant in forms,the two relationships below must hold:

which leads to

Equation(26)clearly indicates that,for the problem considered in this work,there is no alternative scaling transformation that can be used to obtain solutions from the existing results.We therefore are able to conclude that Mayari’s conclusion[27]on that nanofluid results can be recovered from the solutions of already solved regular Newtonian fluid problems by simple arithmetic operations is only valid for several special cases in nanofluid researches.

Table 3 Computation of nanoparticles related parameters.?

5 Results

It is known that Eq.(17)contains an unknown constant σ,which requires an additional boundary condition.In the studies on channel flow problems,it is a common practice to sprecify the mass flow rate as a prescribed quantity.We thus obtain

which can be simplified,by using the similarity variables(16),to

where Umis the constant average flow velocity across the channel,and λ=2Um/Ur.For convenience,we let Um=Urwhich leads to λ=2.

The homotopy analysis method(HAM)is used to solve this flow problem.Since the similar HAM procedures are available in Refs.[29–30],we omit the detailed process but just give the core information as shown in Table 4.

To check the accuracy of our solutions,we define the following functions to evaluate errors:

where

When all physical parameters are prescribed,the cor-responding errors can be obtained.For example,if we set Ra=10,γ=1/100,K1=1,N1=1/10,N2= ?1/10,and θw=1/10,and prescribe ?1and ?2for a range of values,at a certain HAM computational order,the errors can be determined by Eq.(29)as shown in Table 5.

Table 4 HAM computation related quantities.

?

Fig.2 Comparisons of U(η)and θ(η)in the case of Ra=10,γ =1/100,K1=1,N1=1/10,N2= ?1/10,θw=1/10.Line with gradients:solutions by Devi and Devi’s model,[22]Line with circles:solutions by our model.(a)?1= ?2=1/10.(b)?1= ?2=1/10.(c)?1= ?2= ?3=3/100.(d)?1= ?2= ?3=3/100.

Further to check the validity and accuracy of our simplified model,we compare our results of velocity and temperature profiles for a hybrid nanofluid in suspension of two types of nanoparticles with those given by Devi and Devi.[22]It can be seen in Figs.2(a)and 2(b)that very good agreement is found.Note that here Al2O3and Cu nanoparticles are chosen for comparison.As shown in Figs.2(a)and 2(b),we also notice that the results by our simplified hybrid nanofluid model match to those given by the generalized hybrid nanofluid model when three types of nanoparticles,namely,Al2O3,Cu and TiO2are employed.

In our computation,it is found that the variation of nanoparticle volumetric fraction plays limited influence on velocity profiles while it has significant effect on the temperature profiles,as shown in Fig.3.This indicates that the increase of nanoparticle volumetric fraction can enhance heat transfer significantly.In another words,this also verifies the fact that the nanofluids have better thermal transport capability than traditional ones.

Physically,the skin friction and the Nusselt number are important quantities to measure the fluid behaviours.Since the flow and heat transfer exhibit similar characters on both walls,we therefore only consider those quantities on the left wall.In this situation,they are defined by

where

Substituting similarity variables in Eq.(16)into Eq.(32),we obtain

where Re=UmH/νnhfis the Reynolds number.

Fig.3 Variation of(a)U(η)and(b) θ(η)with ?2for Ra=10, γ =1/100,K1=1,N1=1/10,N2= ?1/10,θw=1/10,?1=25/1000.

To test the effects of the nanoparticle volumetric fractions ?1and ?2on various physical quantities,we select Al2O3and Cu nanoparticles in following analysis.As shown in Fig.4(a),for a certain value of ?1,the absolute value of the skin friction coefficient CfLreduces as ?2enlarges.Similarly,when ?2is prescribed,the absolute value of CfLdecreases as ?1evolves.This clearly shows that the nanofluids can effectively diminish the skin friction.The slip effects between the velocities of nanoparticles and the base fluid is the key factor to affect the Nusselt number.The trend of NuLvaries with ?2is similar to that of CfL,namely,when ?1is given,the increase of ?2causes the decrease of the absolute value of NuL,or verse visa,as shown in Fig.4(b).As concluded by Buongiorno,[5]the temperature difference between the walls and the fluid can alter the temperature gradient and thermophoresis,which could result in a significant decrease of viscosity within the boundary layer,thus leading to heat transfer enhancement.

The effect of nanoparticle volumetric fractions on the pressure constant is shown in Fig.4(c).For a given value of ?1,it is found that the pressure decreases gradually as ?2grows.Same trend is found for the variation of the pressure constant σ with ?1at a fixed value of ?2.This reflects another aspect that the flow velocity reduces owing to the reduction of skin friction caused by the increase of nanoparticle volumetric fractions,either for nanofluids or hybrid nanofluids.

In microchannel studies,the slip of the channel wall is of great importance to alter flow and heat transfer behaviours.Take the hybrid nanofluid containing Al2O3,Cu and TiO2nanopaticles as an example.As shown in Fig.5(a),the absolute value of the skin friction coefficient CfLdecreases monotonously as N1grows.However,the absolute value of the Nusselt number NuLincreases continuously as N1increases,as shown in Fig.5(b),while the pressure constant σ reduces gradually as N1enlarges,as shown in Fig.5(c).It is seen from Fig.6(a)that the increase of N1leads to the increase of the velocity near the left wall.This velocity variation leads to the enhancement of temperature in the channel,as presented in Fig.6(b).Physically,the increase of the slip length indicates the decrease of the skin friction,which leads to the increase of fl ow velocity near that wall.Nevertheless,due to the conservation of flow flux,the flow velocity far from the left wall decreases with N1increasing.

Fig.4 Variation of(a)reduced CfL,(b)reduced NuLand(c) σ with ?2for some values of ?1in the case of Ra=10,γ=1/100,K1=1,N1=1/10,N2=?1/10,θw=1/10.

Fig.5 Variation of(a)reduced CfL,(b)reduced NuLand(c)σ with N1in the case of Ra=10,γ=1/100,K1=1,N2=0,θw=1/10,and ?1= ?2= ?3=3/100.

Fig.6 Variation of(a)U(η)and(b)θ(η)with N1in the case of Ra=10,γ =1/100,K1=1,N2=0,θw=1/10,and ?1= ?2= ?3=3/100.

6 Conclusion

The generalized hybrid nanofluid model and its simplifi ed form have been proposed to study the flow and heat transfer behaviours of a hybrid nanofluid convection in a vertical microchannel.It has been found that when the solution is dilute,our simplified model can well predict the flow and heat transfer behaviours of hybrid nanofluids.The argument by Magyari[27]with regards to the homogenous model for expressions of nanofluid solutions by the results of already solved regular Newtonian fluid problems via simple arithmetic operations has been found problematic when it is applied to hybrid nanofluid flows.The effects of various parameters on important physical quantities are analysed and discussed with the following conclusions can be reached:

(i) The variations of nanoparticle volumetric fractions have more obvious effects on temperature distribution than on velocity distribution.

(ii)The nanoparticle volumetric fractions play a significant role on altering flow and heat transfer behaviours.

(iii)The slip effect of the channel wall are of great importance to affect the flow and heat transfer behaviours.