# Heat Transfer and Mass Diffusion in Nanofluids over a Moving Permeable Convective Surface

Authors: Muhammad Qasim; Ilyas Khan; Sharidan Shafie
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Heat transfer and mass diffusion in nanofluid over a permeable moving surface are investigated. The surface exhibits convective boundary conditions and constant mass diffusion. Effects of Brownian motion and thermophoresis are considered. The resulting partial differential equations are reduced into coupled nonlinear ordinary differential equations using suitable transformations. Shooting technique is implemented for the numerical solution. Velocity, temperature, and concentration profiles are analyzed for different key parameters entering into the problem. Performed comparative study shows an excellent agreement with the previous analysis.


---

## Body

## 1. Introduction

Usually, the conventional heat transfer fluids such as oil, water, and ethylene glycol mixtures are poor heat transfer fluids because of their poor thermal conductivity. Therefore, several attempts have been made by many researchers to enhance the thermal conductivity of these fluids by suspending nano/microparticles in liquids. Amongst them Choi [1] was the first who introduced a fluid with enhanced thermal conductivity known as nanofluid. At present, the flow problem involving nanofluids has attracted the investigators to the field. These fluids are engineered colloidal suspensions of nanoparticles (nanometer-sized particles of metals, oxides, nitrides, carbides, or nanotubes) in the ordinary base fluid. Thermal conductivity of nanofluids is higher than the base fluids. Such fluids over a moving surface with heat transfer seem to be very important in microelectronics, fuel cells, hybrid-powered engines, and pharmaceutical processes. It should be pointed out that several metallurgical processes involve the cooling of continuous strips or filaments by drawing them through a nanofluid [2]. Having that in mind, the boundary layer flows of nanofluids have been studied through different approaches, in the recent attempts [3–15] there has been an increasing interest of the recent researchers in the studies of boundary layer flows over a moving surface with heat transfer. This interest in fact stems from several applications of such flows in aerodynamic extrusion process, paper production, food process, glass fiber production, metallurgical process, and so forth. Sakiadis [16] made an initial attempt for boundary layer flow induced by a continuous solid surface moving with constant speed. Tsou et al. [17] extended the work of Sakiadis [16] to heat transfer concept in the boundary layer flow on a continuous moving surface. Abdelhafez [18] studied the boundary layer flow over a continuous moving flat surface in a parallel free stream. He studied the case when surface and free stream move in the same direction. Afzal et al. [19] revisited similar problem and reported the case when the surface and free stream move in the opposite directions. Ishak et al. [20] extended the work of Afzal et al. [19] by considering viscous dissipation and mass transfer effects. They investigated both the cases when surface and free stream move in the same direction and opposite directions. They obtained the numerical solution of the problem by a finite difference scheme known as Keller box method. Very recently, Aziz [21] obtained a similarity solution for Blasius flow of a viscous fluid employing convective boundary conditions. Hayat et al. [22] investigated the flow of Maxwell fluid over a stretching sheet with convective boundary conditions. Some more interesting problems with convective boundary conditions have been reported [23–29].The present paper concentrates on the numerical study of the boundary layer flow of a nanofluid over a permeable moving surface. Mathematical model is constructed in the presence of Brownian and thermophoresis effects. Governing nonlinear analysis is computed by shooting method. Results are compared and analyzed in detail.


## 2. Problem Statement

We study flow of nanofluid over a moving permeable surface with constant velocityuw in the parallel direction of the uniform free stream velocity u∞. The constant temperature and concentration of wall are, Tw and Cw, respectively. The ambient values of temperature and concentration are respectively, T∞ and C∞. In Cartesian coordinate system, x- and y-axes are chosen parallel and perpendicular to the moving surface (see Figure 1).Figure 1
Physical flow model.The boundary layer equations (in absence of viscous dissipation) for the physical problems under examination are(1)∂u∂x+∂v∂y=0,(2)u∂u∂x+v∂u∂y=-1ρf∂p∂x+ν∂2u∂y2,(3)u∂T∂x+v∂T∂z=α∂2T∂y2+τ[DB(∂C∂y∂T∂y)+DTT∞(∂T∂y)2],(4)u∂C∂x+v∂C∂z=DB∂2C∂y2+DTT∞∂2T∂y2,
where u and v are the components of the velocity along the x- and y-directions, respectively, ρf is the density of the base fluid, ν(=μ/ρf) is the kinematic viscosity, α is the thermal diffusivity, DB is the Brownian motion coefficient, DT is the thermophoretic diffusion coefficient, and τ=((ρc)p/(ρc)f) is the ratio of effective heat capacity of the nanoparticle material to the heat capacity of the fluid.The boundary conditions in view of physics of the present problem are(5)u=uw,v=vw,-k∂T∂y=hf(Tf-T),C=Cwaty=0,u⟶u∞,T⟶T∞,C⟶C∞asy⟶∞.
It seems worth mentioning to point out that uw=0 corresponds to the Blasius problem and for u∞=0 we have the Sakiadis problem. Here h is the heat transfer coefficient and  Tf is the convective fluid temperature below the moving surface. A stream of cold fluid at temperature T∞ moving over the right surface of the plate with uniform velocity u∞ while the left surface of the plate is heated below by the convection from the hot fluid of temperature Tf which provides a heat transfer coefficient hf. As a result, convective boundary conditions arise.We define the dimensionless quantities given by(6)ψ=xνUf(η),θ(η)=T-T∞Tf-T∞,ϕ(η)=C-C∞Cw-C∞,η=Uxνy,
where U=uw+u∞, and the free stream function ψ satisfies
(7)u=∂ψ∂y,v=-∂ψ∂x.
The above expression also satisfies the continuity equation (1) and (2)–(5) are reduced to the following forms:
(8)f′′′+12ff′′=0,θ′′+Pr2fθ′+PrNbθ′ϕ′+PrNtθ′2=0,ϕ′′+Le2fϕ′+NtNbθ′′=0,f(0)=S,f′(0)=ϵ,f′(∞)=1-ϵ,θ′(0)=-γ(1-θ(0)),θ(∞)⟶0,ϕ(0)=1,ϕ(∞)=0.
Here primes denote differentiation with respect to η,f(0)=S with S>0 corresponding to suction case and S<0 implying injection, Pr(=ν/α) is the Prandtl number, Le(=ν/DB) is the Lewis number, Nb(=(ρc)pDB(Cw-C∞)/(ρc)fν) is the Brownian motion parameter, Nt(=(ρc)pDT(Tf-T∞)/(ρc)fT∞ν) is the thermophoresis parameter, γ=(hf/kνx/U) is the Biot number, and ϵ(=uw/U) is the velocity ratio parameter. Further we noticed that ϵ=0 corresponds to the flow over a stationary surface caused by the free stream velocity while ϵ=1 is subjected to a moving plate in an ambient fluid, respectively. The case 0<ϵ<1 holds when the plate and fluid are moving in the same direction. If ϵ<0, the free stream is directed towards the positive x-direction, while the plate moves towards negative x-direction. On the other hand, if ϵ>1, the free stream is directed towards negative x-direction, while the plate moves towards the positive x-direction. Here we only discussed the case when 0≤ϵ≤1.Expressions for the local Nusselt numberNux and the local Sherwood number Shx are
(9)Nux=xqwk(Tw-T∞),Shx=xjwDB(Cw-C∞),
where the wall heat flux qw and the mass flux jw are given by
(10)qw=-k(∂T∂y)y=0,jw=-De(∂C∂y)y=0.
Dimensionless form of (10) is given by
(11)(Rew-Re∞)-1/2Nux=-θ′(0),(Rew-Re∞)-1/2Shx=-ϕ′(0),
where the Reynolds numbers are defined as
(12)Rew=uwxν,Re∞=u∞xν.


## 3. Results and Discussion

Here the velocity, temperature, and concentration profiles are analyzed for the velocity ratioϵ, suction parameter S, Brownian motion parameter Nb, thermophoresis parameter Nt, and Lewis number Le. Such theme is achieved through the plots of Figures 2–12 which are sketched. Figure 2 describes the effect of ϵ on f′. It is found that initially f′ decreases, but after η=1.0, it increases when ϵ decreases. Figures 3 and 4 study the influence of S on f′ when ϵ=0.0 and ϵ=0.3, respectively. The boundary layer thickness is found to decrease with the increasing values of S. Sucking fluid particles through porous wall reduces the growth of the boundary layer. This is quite reliable as the fact that suction causes reduction in the boundary layer thickness. Hence a porous character of wall provides a powerful mechanism for controlling the momentum boundary layer thickness. Influence of parameter γ on dimensionless temperature θ is seen in Figure 5. A gradual increase in γ increases the thermal boundary layer thickness. This is expected because the thermal resistance on the hot fluid side is proportional to hf. Hence when γ increases, the hot fluid side convection decreases and consequently the surface temperature increases. Also for α→∞, the result approaches the classical solution for the constant surface temperature. For fixed values of cold fluid properties and free stream velocity, α at any location x is directly proportional to heat transfer coefficient associated with the hot fluid, namely, hf. The thermal resistance on the hot fluid side is inversely proportional to hf. Thus when α increases, then hot fluid side convection resistance decreases and, consequently, the surface temperature increases [21]. Figure 6 elucidates the effects of S on θ. Temperature field θ decreases when S increases. The thermal boundary layer thickness also decreases by increasing S. Effects of thermophoresis parameter Nt and Brownian motion parameter Nb on the temperature θ are shown in Figures 7 and 8. An appreciable increase in the temperature and thermal boundary layer thickness is noticed with an increase in Nt and Nb. The Brownian motion of nanoparticles contributes to thermal conduction enhancement and hence both the temperature and thermal boundary layer thickness increase. It is also noticed that such increase is larger in the case of Nb when compared with Nt. Figure 9 illustrates the effect of Lewis number Le on mass fraction field ϕ. An increase in Le leads to a decrease in molecular diffusivity. Thus increasing the values of Le gradually decreases the concentration of boundary layer. Figures 10 and 11 are plotted to show the effects of thermophoresis parameter Nt and Brownian motion parameter Nb, respectively, on concentration of field ϕ. It is observed that an increase in Nb increases the concentration boundary layer whereas an increase in Nt causes a decrease in ϕ. Figure 12 shows the effect of γ on the mass fraction field ϕ. It is also observed that Nt=0=Nb corresponds to the case when there is no transport driven by the moment of nanoparticles from the surface to the fluid. Further for ϵ=0.0, our results are in excellent agreement with those presented in [21] (see Table 1).Table 1
Values ofθ(0) and -θ′(0) for various values of γ when Nb=0=Nt=S=ϵ=0.0.
γ
θ
(
0
)
-
θ
′
(
0
)
θ
(
0
)
-
θ
′
(
0
)
Pr
=
0.1
Pr
=
10
[21]
Present
[21]
Present
[21]
Present
[21]
Present
0.05
0.2536
0.25362
0.0373
0.03733
0.0643
0.06427
0.0468
0.04683
0.10
0.4046
0.40463
0.0594
0.05941
0.1208
0.12075
0.0879
0.08778
0.20
0.5761
0.57613
0.848
0.84821
0.2155
0.21552
0.1569
0.15685
0.40
0.7310
0.73104
0.1076
0.10762
0.3546
0.35458
0.2582
0.25823
0.60
0.8030
0.80306
0.1182
0.11824
0.4518
0.45177
0.3289
0.32877
0.80
0.8446
0.84463
0.1243
0.12434
0.5235
0.52354
0.3812
0.38134
1
0.8717
0.87172
0.1283
0.12833
0.5787
0.57868
0.4213
0.42125
5
0.9714
0.97144
0.1430
0.14395
0.8729
0.87289
0.6356
0.63558
10
0.9855
0.98554
0.1450
0.14597
0.9321
0.93207
0.6787
0.67856
20
0.9927
0.99275
0.1461
0.14607
0.9649
0.96491
0.7026
0.70246Figure 2
Velocity profilef′(η) for various values of ϵ when S=1.0.Figure 3
Velocity profilef′(η) for various values of S when ϵ=0.0.Figure 4
Velocity profilef′(η) for various values of S when ϵ=0.0.Figure 5
Temperature profileθ(η) for various values of γ when Pr=0.7,  S=1.0,  Nt=Nb=0.1, and ϵ=0.3.Figure 6
Temperature profileθ(η) for various values of  S  when Pr=0.7,  γ=1.0,  Nb=0.2,  Le=5.0, and ϵ=0.3.Figure 7
Temperature profileθ(η) for various values of Nt when Pr=0.7,  γ=1.0,  Nb=0.2,  Le=5.0, and ϵ=0.3.Figure 8
Temperature profileθ(η) for various values of Nb when Pr=0.7,  γ=1.0,Nt=0.3, Le=5.0, and ϵ=0.3.Figure 9
Concentration profileϕ(η) for various values of Le when Pr=0.7,  γ=1.0,  Nt=0.3,  Nt=0.1, and ϵ=0.3.Figure 10
Concentration profileϕ(η) for various values of  Nt when Pr=0.7,  γ=1.0,  Nb=0.3, and ϵ=0.3.Figure 11
Concentration profileϕ(η) for various values of Nb when Pr=0.7,  γ=1.0,Nt=0.3, and ϵ=0.3.Figure 12
Concentration profileϕ(η) for various values of γ when Pr=0.7,  Nb=0.1,  Nt=0.3, and ϵ=0.3.


### 3.1. Conclusions

An incompressible two-dimensional boundary layer flow of nanofluids past a permeable moving surface with convective boundary conditions is studied numerically. The governing boundary layer equations are converted into highly nonlinear coupled ordinary differential equations using some suitable transformations. The resulting equations are solved numerically using shooting technique. Based on the results, the following conclusions are drawn.(i)
f′ decreases initially, but after η=1.0, it increases when ϵ decreases.(ii)
Boundary layer thickness decreases with the increasing values ofS.(iii)
Porous character of wall provides a powerful mechanism for controlling the momentum of boundary layer thickness.(iv)
Asγ increases, it increases the thermal boundary layer thickness.(v)
Temperature decreases whenS increases.(vi)
Temperature and thermal boundary layer thickness increase with increasingNt and Nb.(vii)
An increase inNb increases the concentration of boundary layer whereas an increase in Nt causes a decrease in ϕ.(viii)
Results in [21] are found to be special cases of the present work.


## 3.1. Conclusions

An incompressible two-dimensional boundary layer flow of nanofluids past a permeable moving surface with convective boundary conditions is studied numerically. The governing boundary layer equations are converted into highly nonlinear coupled ordinary differential equations using some suitable transformations. The resulting equations are solved numerically using shooting technique. Based on the results, the following conclusions are drawn.(i)
f′ decreases initially, but after η=1.0, it increases when ϵ decreases.(ii)
Boundary layer thickness decreases with the increasing values ofS.(iii)
Porous character of wall provides a powerful mechanism for controlling the momentum of boundary layer thickness.(iv)
Asγ increases, it increases the thermal boundary layer thickness.(v)
Temperature decreases whenS increases.(vi)
Temperature and thermal boundary layer thickness increase with increasingNt and Nb.(vii)
An increase inNb increases the concentration of boundary layer whereas an increase in Nt causes a decrease in ϕ.(viii)
Results in [21] are found to be special cases of the present work.



---

*Source: 254973-2013-12-09.xml*