# Nanofluid Flow over a Permeable Surface with  Convective Boundary Conditions and Radiative Heat Transfer

Authors: P. K. Kameswaran; P. Sibanda; A. S. N. Murti
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

We investigate the effects of thermal radiation and convective boundary conditions on heat and mass transfer in nanofluid flow over a permeable flat plate. The mathematical model for the nanofluid incorporates variations in the nanoparticle volume fraction of up to 20%. The performance of two water-based nanofluids, namely, stable suspensions of copper and gold nanoparticles in water was investigated. The governing partial differential equations were transformed into ordinary ones using a similarity transformation and solved numerically. The numerical results were validated by comparison with previously published results in the literature. The main focus of this paper is to study the fluid and surface parameters such as the radiation parameter, and suction/injection parameter, solute concentration profiles, as well as the skin friction coefficient and heat and mass transfer rates were conducted.


---

## Body

## 1. Introduction

Heat and mass transfer characteristics of a Newtonian fluid flow past a vertical plate have received considerable attention in recent years because of its extensive applications in engineering processes. Forced convection over a flat plate has been widely studied from both theoretical and experimental standpoints over the past few decades. Early studies were mainly concerned with investigating the similarity solutions of boundary layer flows, particularly flow past a flat plate placed in a uniform free stream first investigated by Blasius [1]. This work was subsequently extended by many authors to explore various aspects of flow and heat transfer on flat surfaces (Pohlhausen [2], Howarth [3], Abu-Sitta [4], Wang [5], Cortell [6], and Sakiadis [7]). Mukhopadhyay et al. [8] presented a study of boundary layer flow and heat transfer over a porous moving plate in the presence of thermal radiation. They found that dual solutions existed when the plate and the fluid moved in opposite directions. In the case of a porous plate, the fluid velocity increases, whereas the temperature decreased for increasing values of the suction parameter.Radiative heat transfer has important applications in physics and engineering including in space technology and other high-temperature processes. Thermal radiation effects may also play an important role in controlling heat transfer in manufacturing processes where the quality of the final product may depend on heat control factors. High-temperature plasmas, the cooling of nuclear reactors, liquid metal fluids, and power generation systems are some important applications of radiative heat transfer from a vertical wall to conductive gray fluids. In high-temperature chemical operations such as in combustion and fire science, it may be necessary to simulate thermal radiation heat transfer effects in combination with conduction, convection, and mass transfer. The effect of radiation on heat transfer problems was studied by Hossain and Takhar [9]. Hamad et al. [10] studied radiation effects on heat and mass transfer in MHD stagnation-point flow over a permeable flat plate with convective surface boundary conditions and temperature-dependent viscosity. They observed that the convective heat transfer parameter lowers the fluid velocity and the wall heat transfer rate. Radiation, however, increases the fluid velocity, temperature, and the heat transfer rate. Singh et al. [11] studied the effects of thermophoresis on hydromagnetic mixed convection and mass transfer flow past a vertical permeable plate with variable suction and thermal radiation. They observed that thermophoresis has a dominant effect on mass transfer in particle deposition processes. Salleh et al. [12] studied free convection over a permeable horizontal flat plate embedded in a porous medium with radiation effects and mixed thermal boundary conditions. They found that the velocity and temperature decreased with an increase in the radiation parameter. A similar solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition was studied by Magyari [13]. He found that the wall temperature decreased monotonically with increasing values of Pr for all values of Bi. Merkin and Pop [14] studied the forced convection flow of a uniform stream over a flat surface with a convective surface boundary condition. They found that the heat transfer near the leading edge is to be dominated by the surface heat flux. The effect of transpiration on self-similar boundary layer flow over moving surfaces was studied by Weidman et al. [15].In recent years tremendous effort has been given to the study of nanofluids. The word nanofluid describes a liquid suspension containing ultra-fine particles (diameter less than 50 nm), Choi [16]. Experimental studies (Masuda et al. [17], Das et al. [18], Pak and Cho [19], and  Xuan and Li [20]) showed that even with a small volumetric fraction of nanoparticles (usually less than 5%), the thermal conductivity of the base liquid is enhanced by 10–50% with a remarkable improvement in the convective heat transfer coefficient.The thermophysical properties of nanofluids such as thermal conductivity, diffusivity, and viscosity have been studied by Kang et al. [21], Velagapudi et al. [22], and Rudyak et al. [23]. The flow of a nanofluid with suction or injection on a moving surface was studied by Bachok et al. [24]. They found that dual solutions exist when the surface and the free stream move in opposite directions. Their results further indicated that suction delays boundary layer separation, while injection accelerated it. Merkin [25] studied the boundary layer flow of a Newtonian fluid past a permeable static or moving flat plate. The literature on nanofluids has been reviewed by Trisaksri and Wongwises [26], Wang and Mujumdar [27], and Eastman et al. [28] among several others. The unsteady stagnation-point flow in a nanofluid was investigated by Bachok et al. [29]. They found that the inclusion of nanoparticles in the water-based fluid produced an increase in the skin friction and heat transfer coefficients. Bachok et al. [30] studied the stagnation-point flow and heat transfer over an exponentially stretching/shrinking sheet in a nanofluid. They observed that the skin friction coefficient and the local Nusselt number are higher for a Cu-water nanofluid compared to those for Al2O3 and TiO2 nanoparticles. Bachok et al. [31] investigated flow and heat transfer characteristics on a moving plate in a nanofluid. They observed that reduced thermal diffusivity leads to higher temperature gradients and, therefore, enhanced heat transfer. These studies were subsequently extended by the authors ( Bachok et al. [32], Bachok et al. [33], Bachok et al. [34], Bachok et al. [35], and  Bachok et al. [36]) to various aspects of heat transfer mechanisms on various flow geometries.The objective of this paper is to discuss the effects of various material and physical parameters such as thermal radiation and nanoparticle volume fraction on heat and mass transfer characteristics in two water-based nanofluids with suction/injection and radiation and solute concentration. The momentum and energy equations are coupled and nonlinear. By using suitable similarity variables these equations are converted into coupled ordinary differential equations and solved numerically. This work has not been considered previously in the literature.


## 2. Mathematical Formulation

Consider steady, incompressible two-dimensional nanofluid flow over a vertical flat plate. The coordinate system is chosen such thatx measures the distance along the plate and y is the distance normal to it. The ambient temperature T and solute concentration  C are denoted by T∞ and C∞, respectively. The fluid is a water-based nanofluid containing either copper Cu or gold Au nanoparticles. The base fluid and the nanoparticles are in thermal equilibrium. Under these assumptions, the boundary layer equations governing the flow, heat, and concentration fields are (Tiwari and Das [37])
(1)∂u∂x+∂v∂y=0,(2)u∂u∂x+v∂u∂y=μnfρnf∂2u∂y2,(3)u∂T∂x+v∂T∂y=αnf∂2T∂y2-1(ρCp)nf∂qr∂y,(4)u∂C∂x+v∂C∂y=D∂2C∂y2.The boundary conditions for (1)–(4) are assumed in the form:
(5)u=0,v=Vw(x),-kf∂T∂y=hf(Tf-T),C=Cwaty=0,u→U∞,T→T∞,C→C∞asy→∞,
where u,v are the velocity components in the x and y directions, respectively, ρ is the density, Cp is the specific heat at constant pressure, qr is the radiative heat flux, D is the species diffusivity, Vw(x) is the velocity at the surface of the plate with Vw(x)>0 for injection, Vw(x)<0 for suction, and Vw(x)=0 corresponds to an impermeable plate. The injection/suction velocity is defined as
(6)Vw(x)=-12U∞νfxfw,
where fw is a constant that determines the transpiration rate at the surface with fw>0 for suction, fw<0 for injection, and fw=0 corresponding to an impermeable surface.At the bottom surface of the plate the fluid is heated via convection from a hot fluid of temperatureTf and heat transfer coefficient hf. Here hf is proportional to x-1/2, that is, hf=cx-1/2, where c is a constant. Following Rosseland's approximation the radiative heat flux qr is modeled as
(7)qr=-4σ*3k*∂T4∂y,
where σ* is the Stefan-Boltzman constant and k* is the mean absorption coefficient. Assuming that the temperature differences within the flow are sufficiently small, then T4 may be expressed as a linear function T4≡4T∞3T-3T∞4 and
(8)∂qr∂y=-16σ*T∞33k*∂2T∂y2.The effective dynamic viscosity of the nanofluid was given by Brinkman [38] as
(9)μnf=μf(1-ϕ)2.5,
where ϕ is the solid volume fraction of nanoparticles. The effective density of the nanofluids is given as
(10)ρnf=(1-ϕ)ρf+ϕρs.The thermal diffusivity of the nanofluid is(11)αnf=knf(ρCp)nf,
where the heat capacitance of the nanofluid is given by
(12)(ρCp)nf=(1-ϕ)(ρCp)f+ϕ(ρCp)s.The thermal conductivity of nanofluids restricted to spherical nanoparticles is approximated by the Maxwell-Garnetss (MG) model (see Maxwell-Garnett [39] and Guérin et al. [40]):
(13)knf=kf[ks+2kf-2ϕ(kf-ks)ks+2kf+ϕ(kf-ks)].Here, the subscriptsnf, f, and s represent the thermophysical properties of the nanofluid, base fluid, and nanosolid particles, respectively. The continuity equation (1) is satisfied by introducing a stream function ψ(x,y) such that
(14)u=∂ψ∂y,v=-∂ψ∂x,
where ψ=(U∞νfx)1/2f(η), f(η) is the dimensionless stream function, and η=(U∞/νfx)1/2y.The velocity components are then given by(15)u=U∞f′(η),v=12(U∞νfx)1/2[ηf′(η)-f(η)].The temperature and solute concentrations are represented as(16)T=T∞+(Tf-T∞)g(η),C=C∞+(Cw-C∞)h(η),
where g(η) is the dimensionless temperature and h(η) is the dimensionless solute concentration. On using (7)–(16), (2)–(5) transform into the following two-point boundary value problem:
(17)f′′′+12ϕ1ff′′=0,(18)(1+43NR)g′′+12Pr(kfknf)ϕ2fg′=0,(19)h′′+12Scfh′=0,(20)f(0)=fw,f′(0)=0,f′(∞)→1,(21)g′(0)=-Bi(1-g(0)),g(∞)→0,(22)h(0)=1,h(∞)→0.The nondimensional constants in (17)–(21) are the radiation parameter NR, the Prandtl number Pr, the Schmidt number Sc, and the Biot number Bi which is used in heat transfer calculations. They are respectively defined as
(23)NR=knfk*4σ*T∞3,Pr=νf(ρCp)fkf,Sc=νfD,Bi=ckfνfU∞,
where
(24)ϕ1=(1-ϕ)2.5[1-ϕ+ϕ(ρsρf)],ϕ2=1-ϕ+ϕ(ρCp)s(ρCp)f.The thermophysical properties of the nanofluid are given in Table1 (see Oztop and Abu-Nada [41], Alloui et al. [42]).Table 1
Thermophysical properties of water and nanoparticles.
Properties→
ρ(kg/m3)
C
p(J/kgK)
k(W/mK)
β
×
1
0
5 (K−1)
Pure water
997.1
4179
0.613
21
Cu
8933
385
401
1.67
Au
19282
129
310
—


## 3. Skin Friction, Heat and Mass Transfer Coefficients

In addition to the velocity and temperature, the quantities of engineering interest in heat and mass transport problems are the skin friction coefficientCf, the local Nusselt number Nux, and the local Sherwood number Shx. These parameters, respectively, characterize the surface drag and wall heat transfer rates. The shearing stress at the surface τw is given by
(25)τw=-μnf[∂u∂y]y=0=-μf(1-ϕ)2.5U∞U∞νfxf′′(0),
where μnf is the coefficient of viscosity. The skin friction coefficient is defined as
(26)Cf=2τwρfu∞2,
and using (25) in (26) we obtain
(27)Cf(1-ϕ)2.5Rex=-2f′′(0).The heat transfer rate at the surface flux at the wall is given by(28)qw=-knf[∂T∂y]y=0=-knf(Tf-T∞)U∞νfxg′(0),
where knf is the thermal conductivity of the nanofluid. The Nusselt number is defined as
(29)Nux=xkfqwTf-T∞.Using (28) in (29) the dimensionless wall heat transfer rate is obtained as mentioned later
(30)NuxRex(kfknf)=-g′(0).The mass flux at the surface of the wall is given by(31)Jw=-D[∂C∂y]y=0=-D(Cw-C∞)U∞νfxh′(0),
and the Sherwood number is defined as
(32)Shx=xDJwCw-C∞.Using (31) in (32) the dimensionless wall mass transfer rate is obtained as
(33)ShxRex=-h′(0).In (27), (30), and (33), Rex represents the local Reynolds number defined by Rex=xU∞/νf.


## 4. Results and Discussion

This study investigated the effects of radiation and nanoparticle volume fraction on heat and mass transfer on nanofluid flow over a permeable flat surface. We considered two different nanoparticles, copperCu and gold Au, with water as the base fluid. The Prandtl number of the base fluid (water) was kept at constant Pr=6.7850.The nonlinear differential equations (17)–(19) with boundary conditions (20)–(22) do not have a closed form solution. The equations were solved numerically using the Matlab bvp4c function. To check the accuracy of the solutions, the local Nusselt number g′(0) for a clear fluid is compared with results reported by Aziz [43] and Ishak [44] in Table 2. From Table 2 we observe that the values of the wall temperature gradient for different values of Bi and Pr when fw=0, ϕ=0, NR=1010, and Sc=0 are in good agreement with results in the literature.Table 2
Comparison of wall temperature gradient-g′(0) for different values of Bi, Pr when fw=0,ϕ=0,NR=1010, and Sc=0 (impermeable surface).
Bi
Aziz [43]
Ishak [44]
Present Results
Pr = 0.1
Pr = 10
Pr = 0.1
Pr = 10
Pr = 0.1
Pr = 10
0.05
0.0373
0.0468
0.036844
0.046787
0.036844
0.046787
0.10
0.0594
0.0879
0.058338
0.087925
0.058338
0.087925
0.20
0.0848
0.1569
0.082363
0.156903
0.082363
0.156903
0.40
0.1076
0.2582
0.103720
0.258174
0.103720
0.258174
0.60
0.1182
0.3289
0.113533
0.328945
0.113533
0.328945
0.80
0.1243
0.3812
0.119170
0.381191
0.119170
0.381191
1
0.1283
0.4213
0.122830
0.421344
0.122830
0.421344
5
0.1430
0.6356
0.136215
0.635583
0.136215
0.635583
10
0.1450
0.6787
0.138096
0.678721
0.138096
0.678721
20
0.1461
0.7026
0.139056
0.702563
0.139056
0.702563Figures1 and 2 show the effect of nanoparticle volume fraction on the velocity, temperature, and concentration profiles, respectively, in the case of Cu-water nanofluid. From Figure 1(a), it is clear that as the nanoparticle volume fraction ϕ increases, the nanofluid velocity profiles width decreases (became narrower) and vice versa. Increasing the volume fraction of nanoparticles increases the thermal conductivity of the nanofluid,  and we observe a thickening of the thermal boundary layer. We also observe that in the dynamic boundary layer, the gradient of the velocity from wall to the free stream is more rapid when ϕ is increased. The same trend is observed in the case of a Au-water nanofluid. We also found that the axial velocity in a Au-water nanofluid is relatively smaller than that of a Cu-water nanofluid.Effects of the nanoparticle volume fractionϕ on (a) velocity and (b) temperature profiles, when fw=0, Bi=10, NR=1, and Sc=1.5.
(a)
(b)Figure 2
Effects of the nanoparticle volume fractionϕ on concentration profiles when fw=0, Bi=10, NR=1, and Sc=1.5.Figures1(b) and 2 illustrate the effect of nanoparticle volume fraction on the temperature and concentration profiles, respectively. It is clear that as the nanoparticle volume fraction increases the nanofluid temperature increases. The conductivity of the nanofluid increases as the nanoparticle volume fraction increases leading to a broadening of the thermal boundary layer profiles, see Hady et al. [45]. Since the thermal conductivity of Cu is higher than that of Au, we observe that the profiles for a Cu-water nanofluid are proportionally broader than those of a Ag-water nanofluid. In addition, we note that the solute concentration decreases with an increase in the nanoparticle volume fraction for both Ag-water and Cu-water nanofluids.The variation of the skin friction coefficient, wall heat, and mass transfer rates with the suction/injection parameterfw for different nanoparticle volume fractions ϕ is shown in Figures 3 and 4. We note that the inclusion of nanoparticles into the base fluid increases the skin friction coefficient. We also observe that skin friction coefficient increases as the nanoparticle volume fraction increases and that the increment is more pronounced in the case of a Au-water nanofluid as compared to a Cu-water nanoparticle. The same trend is observed in the case of the mass transfer coefficient. Figure 3(b) shows that the lowest heat transfer rate is obtained for Au nanoparticles due to the domination of the conduction mode of heat transfer. This is due to a lower thermal conductivity of Au nanoparticles compared with Cu nanoparticles.Effects of suction/injection parameterfw on (a) skin friction and (b) heat transfer coefficient, when NR=1, Bi=10, and Sc=1.5.
(a)
(b)Figure 4
Effects of suction/injection parameterfw on mass transfer coefficient, when NR=1, Bi=10, and Sc=1.5.Figures5 and 6 show the velocity, temperature, and concentration distributions for different values of the suction or injection parameter fw. The parameter Bi has no influence on the flow field, which is also clear from (17)–(22). We note that the velocity gradient at the surface, which is proportional to the skin friction coefficient, increases with increasing fw for both nanofluids. Thus the surface shear stress is higher for suction (fw>0) compared to injection (fw<0). We also observe that for all values of fw the radial velocity f′(η) attains a maximum value close to the surface of the wall and decreases monotonically along the boundary layer for both nanofluids. Figures 5(b) and 6 show that the temperature and concentration profiles at the surface are higher for suction compared to those for injection in the case of a Cu-water nanofluid. This is due to the fact that the surface shear stress increases when suction is introduced, which in turn increases the local Nusselt number. We note that the lowest heat transfer rate is obtained in the case of Au nanoparticles due to the domination of the conduction mode of heat transfer. This is due to the fact that Au nanoparticles have a lower value of thermal conductivity compared to Cu nanoparticles. These findings are similar to the results obtained by Ishak [44] in the case of newtonian fluids (ϕ=0).Effects of suction/injection parameterfw on (a) velocity and (b) temperature profiles, when ϕ=0.1, Bi=10, NR=1, and Sc=1.5.
(a)
(b)Figure 6
Effects of suction/injection parameterfw on concentration, when ϕ=0.1, Bi=10, NR=1 and Sc=1.5.The effects of thermal radiationNR and the Biot number Bi on the temperature profiles are shown in Figure 7. The temperature decreases as the thermal radiation increases. This in turn leads to an increase in the heat transfer rate, Hady et al. [45]. The drop in the temperature is higher for a Au-water nanofluid compared to a Cu-water nanofluid, and Cu nanoparticles have higher temperature distribution values than the Au nanoparticles.Effects of (a) thermal radiationNR, when Bi=10 and (b) effects of Bi,  when NR=1 on temperature for various values of ϕ=0.1, fw=1, and Sc=1.5.
(a)
(b)The temperature profiles for different values of the Biot number are shown in Figure7(b). The temperature profiles for a clear fluid (when Pr=0.72 and fw=1) match with the results obtained by Aziz [43]. The surface temperature increases with Bi for both nanofluids. The parameter Bi is directly proportional to the heat transfer coefficient associated with the hot fluid hf. The thermal resistance on the hot fluid side is inversely proportional to hf. Thus as Bi increases, the hot fluid side convection resistance decreases, and consequently, the surface temperature increases.Figures8, 9, and 10 show the behaviour of the heat and mass transfer rates under the influence of NR, Bi, and Sc, respectively, for the two different nanofluids when Pr=6.7850 and ϕ=0,0.1, and 0.2. The change in the heat transfer rate is important in cooling and heating processes. Figures 8 and 9 show that an increase in the thermal radiation parameter and Biot number reduces the heat transfer rate. Figure 10 shows that  by increasing the Schmidt numbers and solid volume fraction, the mass transfer rate increases for both nanofluids, and this increment is more in the case of Au nanoparticles compared to Cu nanoparticles.Figure 8
Effects of thermal radiation parameterNR on heat transfer rate, when fw=1, Bi=1, and Sc=1.5.Figure 9
Effects ofBi on heat transfer rate, when fw=1, NR=1, and Sc=1.5.Figure 10
Effects ofSc on mass transfer rate, when fw=1, NR=1, and Bi=10.The temperature profiles for various physical parameters in the caseBi→∞ (i.e., g(0)=1) are displayed in Figures 11–13. We note, from the convective boundary conditions (21), that g(η)→1 as Bi→∞. The convective parameter Bi has no influence on the flow field. We observe further that, as the nanoparticle volume fraction increases, the nanofluid temperature increases which is consistent with expected physical behaviour and may be attributed to the fact that as the thermal conductivity of the nanofluid increases, thermal boundary layer also increases. The temperature profiles for various suction/injection parameter values are shown in Figure 12. We observe that heat transfer rate at the surface is higher in the case of suction compared to injection due to an increase in the shear stress when suction is introduced in the case of a Cu-water nanofluid. The same trend is observed in the case of Au nanoparticles, but the change in the temperature is less significant than for Cu nanoparticles due to differences in thermal conductivity.Figure 11
Effects ofϕ on temperature, when fw=0, NR=1, Sc=1.5, and Bi→∞.Figure 12
Effects of suction/injection parameterfw on temperature, when NR=1, Sc=1.5, ϕ=0.1, and Bi→∞.Figure 13
Effects of thermal radiationNR on temperature, when fw=1, Sc=1.5, ϕ=0.1, and Bi→∞.Figure13 shows the effect of the radiation parameter on the temperature profile. An increase in the radiation parameter reduces the thermal boundary layer thickness across the plate. The temperature g(η) decreases as the radiation parameter NR increases. This may be explained by the fact that a decrease in the radiation parameter NR=knfk*/4σ*T∞3 for given k* and T∞ implies a decrease in the Rosseland radiation absorptivity. The same trend is observed for both Cu and Au nanoparticles. It was also observed that the change in temperature profiles in the case of Au-water nanofluid was slightly less significant than for a Cu-water nanofluid.


## 5. Conclusions

The effects of thermal radiation on heat and mass transfer in a nanofluid flow over a permeable flat plate were studied. The governing equations were transformed into a set of coupled nonlinear differential equations and solved by using the Matlab bvp4c ODE Solver. The results were compared for selected values of the governing physical parameters with results available in the literature. We found that the velocity gradient at the surface increases with increasingfw for both nanofluids considered. The surface temperature increases as Biot number increases for both nanofluids.



---

*Source: 201219-2013-03-21.xml*