AOSIS OpenJournals
http://www.openjournals.net
info@openjournals.net
http://www.sajs.co.za
0038-2353
0038-2353
<article_title>Effects of internal heat generation, thermal radiation and buoyancy force on a boundary layer over a vertical plate with a convective surface boundary condition
</article_title> In this paper we analyse the effects of internal heat generation, thermal radiation and buoyancy force on the laminar boundary layer about a vertical plate in a uniform stream of fluid under a convective surface boundary condition. In the analysis, we assumed that the left surface of the plate is in contact with a hot fluid whilst a stream of cold fluid flows steadily over the right surface; the heat source decays exponentially outwards from the surface of the plate. The similarity variable method was applied to the steady state governing non-linear partial differential equations, which were transformed into a set of coupled non-linear ordinary differential equations and were solved numerically by applying a shooting iteration technique together with a sixth-order Runge–Kutta integration scheme for better accuracy. The effects of the Prandtl number, the local Biot number, the internal heat generation parameter, thermal radiation and the local Grashof number on the velocity and temperature profiles are illustrated and interpreted in physical terms. A comparison with previously published results on similar special cases showed excellent agreement.
Philip O. Olanrewaju1
1Department of Mathematics, Covenant University, Ota, Ogun State, Nigeria
Jacob A. Gbadeyan1
1Department of Mathematics, Covenant University, Ota, Ogun State, Nigeria
Tasawar Hayat2
2Department of Mathematics, Quaid-I-Azam University, Islamabad, Pakistan
Awatif A. Hendi3
3Department of Physics, Faculty of Science, King Saud University, Riyadh, Saudi Arabia
Philip Olanrewaju
oladapo_anu@yahoo.ie
Department of Mathematics, Covenant University, Ota, Ogun State, PMB 1023, Nigeria
476
107
9/10
10.4102/sajs.v107i9/10.476
14 Oct. 2010
03 May 2011
07 Sept. 2011 Olanrewaju PO, Gbadeyan JA, Hayat T, Hendi AA. Effects of internal heat generation, thermal radiation and buoyancy force on a boundary layer over a vertical plate with a convective surface boundary condition. S Afr J Sci. 2011;107(9/10), Art. #476, 6 pages. doi:10.4102/sajs.v107i9/10.476
© 2011. The Authors.
Licensee: AOSIS OpenJournals. This work
is licensed under the
Creative Commons
Attribution License. IntroductionBoundary-layer flows over a moving or stretching plate are of great importance in view of their relevance to a wide variety of technical applications, particularly in the manufacture of fibres in glass and polymer industries. The first and foremost work regarding boundary-layer behaviour in moving surfaces in a quiescent fluid was performed by Sakiadis.1
Subsequently, many researchers
2,3,4,5,6,7,8,9
worked on the problem of moving or stretching plates under different situations. In the boundary-layer theory, similarity solutions were found to be useful in the interpretation of certain fluid motions at large Reynolds numbers. Similarity solutions often exist for the flow over semi-infinite plates and stagnation point flow for two-dimensional, axisymmetrical and three-dimensional bodies. In special cases, when there is no similarity solution, one has to solve a system of non-linear partial differential equations. For similarity boundary-layer flows, velocity profiles are similar. But this kind of similarity is lost for non-similarity flows.
10,11,12,13,14
Obviously, the non-similarity boundary-layer flows are more general in nature and are more important, not only in theory but also in application.
The heat-transfer analysis of boundary-layer flows with radiation is also important in electrical power generation, astrophysical flows, solar power technology, space vehicle re-entry and other industrial areas. Extensive literature that deals with flows in the presence of radiation effects is now available. Raptis et al.
15 studied the effect of thermal radiation on the magnetohydrodynamic flow of a viscous fluid past a semi-infinite stationary plate. Hayat et al.
16 extended the analysis of reference
15 for a second-grade fluid.
Convective heat transfer studies are very important in processes involving high temperatures, such as gas turbines, nuclear plants and thermal energy storage. Recently, Ishak
17 examined the similarity solutions for flow and heat transfer over a permeable surface with convective boundary condition. Moreover, Aziz
18,19
studied a similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition and also studied hydrodynamic and thermal slip flow boundary layers over a flat plate with a constant heat flux boundary condition. Very recently, Makinde and Olanrewaju
20 investigated the buoyancy effects on a thermal boundary layer over a vertical plate with a convective surface boundary condition.
In this study, the recent work of Ishak
17 , Aziz
18 and Makinde and Olanrewaju
20 was extended to include the effect of thermal radiation and internal heat generation. The numerical solutions of the resulting momentum and the thermal similarity equations are reported for representative values of the thermophysical parameters embedded in the fluid-convection process. The objective of this paper was to explore the effects of thermal radiation and internal heat generation on the fluid under a convective surface boundary condition. The non-linear equations governing the flow were solved numerically using a shooting technique together with a sixth-order Runge–Kutta integration scheme, which gives better accuracy than a fourth-order Runge–Kutta method. Graphical results are first reported for emerging parameters and then discussed.
Mathematical formulation
We considered a two-dimensional steady incompressible fluid flow coupled with heat transfer by convection over a vertical plate. A stream of cold fluid at temperature T∞ moved over the right surface of the plate with a uniform velocity U∞ whilst the left surface of the plate was heated by convection from a hot fluid at temperature Tf , which provided a heat transfer coefficient hf . The density variation as a result of buoyancy force effects was taken into account in the momentum equation and the thermal radiation and the internal heat generation effects were taken into account in the energy equation (the Boussinesq approximation). The continuity, momentum and energy equations describing the flow are, respectively:
1 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2568
2 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2569
3 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2570
where u and v are the x (along the plate) and the y (normal to the plate) components of the velocities, respectively. T is the temperature, υ is the kinematics viscosity of the fluid, α is the thermal diffusivity of the fluid, β is the thermal expansion coefficient, Q is the heat released per unit per mass, g is the gravitational acceleration, qr is the radiative heat flux and k is the thermal conductivity. The velocity boundary conditions can be expressed as
4 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2571
and 5 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2572
The boundary conditions at the plate surface and far into the cold fluid may be written as
6 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2573
and 7 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2574
The radiative heat flux qr is described by the Rosseland approximation such that
8 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2575
where σ* and K are the Stefan-Boltzmann constant and the mean absorption coefficient, respectively. Following Chamkha21, we assume that the temperature differences within the flow are sufficiently small such that T 4 can be expressed as a linear function after using the Taylor series to expand T 4 about the free stream temperature T∞ and neglecting higher-order terms. This result is the following approximation:
9 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2576
Using [Eqn 8] and [Eqn 9] in [Eqn 3], we obtain
10 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2577
We then introduce a similarity variable η and a dimensionless stream function f(η) and temperature θ(η) as
11 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2578
where the prime symbol denotes differentiation with respect to η and Rex = U∞x/υ is the local Reynolds number. [Eqn 1] to [Eqn 7] reduce to:
12
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2579
13 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2580
14 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2581
and
15 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2582
where 16 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2583
Bix is the local Biot number, Pr is the Prandtl number, Grx is the local Grashof number, Ra is the radiation parameter and λx is the internal heat generation parameter. For the momentum and energy equations to have a similarity solution, the parameters Grx, λx and Bix must be constants and not functions of x as in [Eqn 16]. This condition can be met if the heat-transfer coefficient hf is proportional to x–½ , the thermal expansion coefficient β is proportional to x–1 and the heat-release coefficient Q is proportional to x–1. We therefore assume
17 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2584
where c, d and m are constants. Substituting [Eqn 17] into [Eqn 16], we get
18 http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2585
Bi, λ and Gr defined by [Eqn 18] are the Biot number, the internal heat generation parameter and the Grashof number, respectively. The solutions of [Eqn 12] to [Eqn 15] yield the similarity solutions. However, the solutions generated are the local similarity solutions whenever Bix, λx and Grx are defined as in [Eqn 13].
Results and discussion
The ordinary differential equations ([Eqn 9] and [Eqn 10]) subject to the boundary conditions ([Eqn 11] and [Eqn 12]) were solved numerically using the symbolic algebra software Maple.
22 Table 1 presents a comparison of the values of –θ' (0) and θ(0) with those reported by Aziz
18 , Ishak
17 and Makinde and Olanrewaju
20 , which show an excellent agreement for
Pr = 0.72. Table 2 shows the values of the skin-friction coefficient f ''(0) and the local Nusselt number –θ'(0), for various values of embedded parameters. From Table 2, it can be seen that the skin friction and the rate of heat transfer at the plate surface increased with an increase in the local Grashof number, the convective surface heat transfer parameter, the internal heat generation parameter and the radiation absorption parameter. However, an increase in the fluid Prandtl number decreased the skin friction but increased the rate of heat transfer at the plate surface. Figures 1–6 depict the fluid velocity profiles. Generally, the fluid velocity is zero at the plate surface and increases gradually away from the plate towards the free stream value satisfying the boundary conditions. It is clearly seen from Figure 1 that the Grashof number had profuse effects on the velocity boundary layer thickness. It is interesting to note that an increase in the intensity of the convective surface heat transfer (Bix) produced a slight increase in the fluid velocity within the boundary layer (Figure 2). The local internal heat generation parameter, the Prandtl number and the local Biot number had little or no influence on the velocity profiles (Figures 3–4 and Figure 6), which could be justified by [Eqn 12], which had a small value for the Grashof number. When Ra = 0.1, the Prandtl number effected the velocity profile (Figure 5) and the negative part could have been caused by the value of the Grashof number used (see [Eqn 12]) and the value of the radiation parameter used (see [Eqn 13]). Figures 7–12 illustrate the fluid temperature profiles within the boundary layer. The fluid temperature was at a maximum at the plate surface and decreased exponentially to zero away from the plate, thus satisfying the boundary conditions. From these figures, it is noteworthy that the thermal boundary layer thickness increased with an increase in Bix, λx and Ra and decreased with increasing values of Grx and Pr. Hence, the convective surface heat transfer, the internal heat generation parameter and the radiation parameter enhanced thermal diffusion whilst an increase in the Prandtl number and the intensity of buoyancy force slowed down the rate of thermal diffusion within the boundary layer. Figure 13 shows the influence of Prandtl numbers on the thermal boundary layer, as obtained by Aziz
18
1
A comparison of the values obtained for the Nusselt number [θ'(0)] and surface temperature [θ(0)] with an increase in the local Biot number (Bix) in this study and by Aziz18, Ishak17 and Makinde and Olanrewaju20 in previous studies.
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2586
2 A comparison of the values of the skin-friction coefficient [ f "(0)], temperature at the wall surface [θ(0)] and the Nusselt number [θ'(0)] for different parameter values embedded in the flow model.
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2587 1
The fluid velocity profile with increasing distance from the plate surface (where f ’(η) = 0) and increasing Grashof number (Grx).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2588
2
The fluid velocity profile with increasing distance from the plate surface (where f ’(η) = 0) and increasing Biot number (Bix).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2589
3 The fluid velocity profile with increasing distance from the plate surface (where f ’(η) = 0) and increasing internal heat generation (λx).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2590
4 The fluid velocity profile with increasing distance from the plate surface (where f ’(η) = 0) and increasing Prandtl number (Pr) when the radiation parameter, Ra = 0.5.
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2591
5 The fluid velocity profile with increasing distance from the plate surface (where f ’(η) = 0) and increasing Prandtl number (Pr) when the radiation parameter, Ra = 0.1.
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2592
6 The fluid velocity profile with increasing distance from the plate surface (where f ’(η) = 0) and increasing radiation (Ra).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2593
7 The fluid temperature profile with increasing distance from the plate surface (where θ(η) = maximum) and increasing Grashof number (Grx).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2594
8 The fluid temperature profile with increasing distance from the plate surface (where θ(η) = maximum) and increasing Biot number (Bix).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2595
9 The fluid temperature profile with increasing distance from the plate surface (where θ(η) = maximum) and increasing internal heat generation (λx).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2596
10 The fluid temperature profile with increasing distance from the plate surface (where θ(η) = maximum) and increasing Prandtl number (Pr) when the radiation parameter, Ra = 0.5.
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2597
11 The fluid temperature profile with increasing distance from the plate surface (where θ(η) = maximum) and increasing Prandtl number (Pr) when the radiation parameter, Ra = 0.1.
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2598
12 The fluid temperature profile with increasing distance from the plate surface (where θ(η) = maximum) and increasing radiation (Ra).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2599
13 The fluid temperature profile obtained by Aziz18 with increasing Prandtl number (Pr).
http://sajs.co.za/index.php/SAJS/article/downloadSuppFile/476/2600
ConclusionsWe analysed the effects of internal heat generation, thermal radiation and buoyancy force on the laminar boundary layer about a vertical plate in a uniform stream of fluid under a convective surface boundary. A similarity solution for the momentum and the thermal boundary layer equations is possible if the convective heat transfer of the fluid heating the plate on its left surface is proportional to x–½ and if the thermal expansion coefficient β and the heat released per unit per mass Q are proportional to x–1. Numerical solutions of the similarity equations were reported for the various parameters embedded in the problem. The combined effect of increasing the Prandtl number and the Grashof number tended to reduce the thermal boundary layer thickness along the plate whilst the effects of increasing the Biot number, the internal heat generation parameter and the radiation absorption parameter enhanced thermal diffusion.
AcknowledgementO.P.O. wishes to thank Covenant University, Ogun State, Nigeria for its generous financial support. A.A.H., as a VP, appreciates the support of King Saud University, Saudi Arabia (KSU-VPP-117).
1.10.1002/aic.690070108
2.Crane IJ Flow past a stretching plate Z Angew Math Phys. 1970 21 56 1 37 3.10.1002/cjce.5450550619 4.10.1002/zamm.19820621009 5.Danberg JE Fansler KS A nonsimilar moving wall boundary-layer problem Quart Appl Math1976 34 305 309 6.Chakrabarti A Gupta AS Hydromagnetic flow and heat transfer over a stretching sheet Quart Appl Math.1979 37 73 78 7.0.1007/BF01450393 8.10.1007/BF01179221 9.10.1016/0017-9310(90)90181-S 10. Wanous KJ Sparrow EM Heat transfer for flow longitudinal to a cylinder with surface mass transfer J Heat Transf ASME Ser C1965 87 1 317 319 11. 10.1098/rspa.1965.0069 12.Sparrow EM Quack H Boerner CJ Local non-similarity boundary layer solutions J AIAA1970 8 11 1936 1942 13.Sparrow EM Yu HS Local non-similarity thermal boundary-layer solutions J Heat Transf ASME1971 93 328 334 14.10.1016/S0020-7462(00)00061-5 15.10.1016/S0096-3003(03)00657-X 16.10.1016/j.ijheatmasstransfer.2006.08.014 17. 10.1016/j.amc.2010.06.026 18.Aziz A A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition Commun Nonlinear Sci Numer Simulat2009 14 1064 1068 19.10.1016/j.cnsns.2009.04.026 20. 10.1115/1.4001386 21.10.1016/S0020-7225(96)00122-X 22.Heck A Introduction to Maple. 3rd ed. New York: Springer-Verlag2003