Page of 6 Effects of internal heat generation, thermal radiation buoyancy force on a boundary layer over a vertical plate with a convective surface boundary condition Authors: Philip O. Olanrewaju Jacob A. Gbadeyan Tasawar Hayat 2 Awatif A. Hendi 3 Affiliations: Department of Mathematics, Covenant University, Ota, Ogun State, Nigeria 2 Department of Mathematics, Quaid-I-Azam University, Islamabad, Pakistan 3 Department of Physics, Faculty of Science, King Saud University, Riyadh, Saudi Arabia Correspondence to: Philip Olanrewaju Email: oladapo_anu@yahoo.ie Postal address: Department of Mathematics, Covenant University, Ota, Ogun State, PMB 23, Nigeria Dates: Received: 4 Oct. 2 Accepted: 3 May 2 Published: 7 Sept. 2 How to cite this article: Olanrewaju PO, Gbadeyan JA, Hayat T, Hendi AA. Effects of internal heat generation, thermal radiation buoyancy force on a boundary layer over a vertical plate with a convective surface boundary condition. S Afr J Sci. 2;7(9/), Art. #476, 6 pages. doi:.42/sajs. v7i9/.476 2. The Authors. Licensee: AOSIS OpenJournals. This work is licensed under the Creative Commons Attribution License. In this paper we analyse the effects of internal heat generation, thermal radiation 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 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 Prtl number, the local Biot number, the internal heat generation parameter, thermal radiation the local Grashof number on the velocity temperature profiles are illustrated interpreted in physical terms. A comparison with previously published results on similar special cases showed excellent agreement. Introduction Boundary-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 polymer industries. The first foremost work regarding boundary-layer behaviour in moving surfaces in a quiescent fluid was performed by Sakiadis. 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 stagnation point flow for twodimensional, axisymmetrical 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.,,2,3,4 Obviously, the non-similarity boundary-layer flows are more general in nature 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 other industrial areas. Extensive literature that deals with flows in the presence of radiation effects is now available. Raptis et al. 5 studied the effect of thermal radiation on the magnetohydrodynamic flow of a viscous fluid past a semi-infinite stationary plate. Hayat et al. 6 extended the analysis of reference 5 for a second-grade fluid. Convective heat transfer studies are very important in processes involving high temperatures, such as gas turbines, nuclear plants thermal energy storage. Recently, Ishak 7 examined the similarity solutions for flow heat transfer over a permeable surface with convective boundary condition. Moreover, Aziz 8,9 studied a similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition also studied hydrodynamic thermal slip flow boundary layers over a flat plate with a constant heat flux boundary condition. Very recently, Makinde Olanrewaju 2 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 7, Aziz 8 Makinde Olanrewaju 2 was extended to include the effect of thermal radiation internal heat generation. The numerical solutions of the resulting momentum the thermal similarity equations are reported for representative values of the thermophysical parameters embedded in the fluid-convection process. The S Afr J Sci 2; 7(9/)
Page 2 of 6 objective of this paper was to explore the effects of thermal radiation 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 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 T f, which provided a heat transfer coefficient h f. The density variation as a result of buoyancy force effects was taken into account in the momentum equation the thermal radiation the internal heat generation effects were taken into account in the energy equation (the Boussinesq approximation The continuity, momentum energy equations describing the flow are, respectively: u v + =, x y u u + v u = υ 2 u + gβ(t T ), x y y 2 u T + v T = α 2 T + Q (T T ) α q r, x y y 2 ρc p k y [Eqn ] [Eqn 2] [Eqn 3] where u v are the x (along the plate) 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, q r is the radiative heat flux k is the thermal conductivity. The velocity boundary conditions can be expressed as q r = 4σ* T 4, 3K y [Eqn 8] where σ* K are the Stefan-Boltzmann constant the mean absorption coefficient, respectively. Following Chamkha 2, 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 exp T 4 about the free stream temperature T neglecting higher-order terms. This result is the following approximation: T 4 4T 3 T 3T 4. Using [Eqn 8] [Eqn 9] in [Eqn 3], we obtain q r 6σ* T 4 =. y 3K y [Eqn 9] [Eqn ] We then introduce a similarity variable a dimensionless stream function f() temperature θ() as U = y y = Re, u = f, v = U υ x (f f ), υx x U 2 x θ = T T, T f T [Eqn ] where the prime symbol denotes differentiation with respect to Re x = U x/υ is the local Reynolds number. [Eqn ] to [Eqn 7] reduce to: f + ff + Grx θ =, 2 θ [ + 4 Ra] + Pr fθ + Pr λx θ =, 3 2 f() = f () =, θ () = [ θ()] [Eqn 2] [Eqn 3] [Eqn 4] u(x, ) = v(x, ) = [Eqn 4] f ( ) =, θ( ) =, [Eqn 5] where u(x, ) = U. [Eqn 5] = h f υx υ, Pr =, Grx = υxgβ(t f T ), k U α U 2 The boundary conditions at the plate surface far into the cold fluid may be written as Ra = 4ασT 3, = xq. kk ρc p U [Eqn 6] k T (x, ) = h f [T f T(x, )] y T(x, ) = T. [Eqn 6] [Eqn 7] The radiative heat flux q r is described by the Rossel approximation such that is the local Biot number, Pr is the Prtl number, is the local Grashof number, Ra is the radiation parameter is the internal heat generation parameter. For the momentum energy equations to have a similarity solution, the parameters, must be constants not functions of x as in [Eqn 6]. This condition can be met if the heat-transfer coefficient h f is proportional to x ½, the thermal expansion coefficient β is proportional to x the heat-release coefficient Q is proportional to x. We therefore assume S Afr J Sci 2; 7(9/)
Page 3 of 6 h f = cx 2, β = mx, Q = dx, [Eqn 7] where c, d m are constants. Substituting [Eqn 7] into [Eqn 6], we get Bi = c υ υmg(t, Gr = f T ) dq, λ =, k U U 2 ρc U p [Eqn 8] Bi, λ Gr defined by [Eqn 8] are the Biot number, the internal heat generation parameter the Grashof number, respectively. The solutions of [Eqn 2] to [Eqn 5] yield the similarity solutions. However, the solutions generated are the local similarity solutions whenever, are defined as in [Eqn 3]. Results discussion The ordinary differential equations ([Eqn 9] [Eqn ]) subject to the boundary conditions ([Eqn ] [Eqn 2]) were solved numerically using the symbolic algebra software Maple. 22 Table presents a comparison of the values of θ' () θ() with those reported by Aziz 8, Ishak 7 Makinde Olanrewaju 2, which show an excellent agreement for Pr =.72. Table 2 shows the values of the skin-friction coefficient f ''() the local Nusselt number θ'(), for various values of embedded parameters. From Table 2, it can be seen that the skin friction 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 the radiation absorption parameter. However, an increase in the fluid Prtl number decreased the skin friction but increased the rate of heat transfer at the plate surface. Figures 6 depict the fluid velocity profiles. Generally, the fluid velocity is zero at the plate surface increases gradually away from the plate towards the free stream value satisfying the boundary conditions. It is clearly seen from Figure 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 ( ) produced a slight increase in the fluid velocity within the boundary layer (Figure 2 The local internal heat generation parameter, the Prtl number the local Biot number had little or no influence on the velocity profiles (Figures 3 4 Figure 6), which could be justified by [Eqn 2], which had a small value for the Grashof number. When Ra =., the Prtl number effected the velocity profile (Figure 5) the negative part could have been caused by the value of the Grashof number used (see [Eqn 2]) the value of the radiation parameter used (see [Eqn 3] Figures 7 2 illustrate the fluid temperature profiles within the boundary layer. The fluid temperature was at a maximum at the plate surface 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, Ra decreased with increasing values of Pr. Hence, the TABLE : A comparison of the values obtained for the Nusselt number [θ'()] surface temperature [θ()] with an increase in the local Biot number ( ) in this study by Aziz 8, Ishak 7 Makinde Olanrewaju 2 in previous studies. Present study Aziz 8 Ishak 7 Makinde Olanrewaju 2 θ'() θ() θ'() θ() θ'() θ'().5.42767.4466.428.447.42767.428..74724.25275.747.2528.74724.747.2.9295.4352.93.435.9295.93.4.69994.575.7.575.69994.7.6.985.6699.98.6699.985.98.8.25864.736.259.732.25864.259..22878.778.2282.778.22878.2282 5..2793.9447.279.944.2793.279..28746.9728.287.973.28746.287 2..29329.98543.293.9854.29329.293 3..292754.9924 - - -.2928 In these computations, the radiation parameter (Ra), the local Grashof number ( ) the internal heat generation parameter ( ) were zero the Prtl number was.72. TABLE 2: A comparison of the values of the skin-friction coefficient [ f "()], temperature at the wall surface [θ()] the Nusselt number [θ'()] for different parameter values embedded in the flow model. Pr Ra f "() θ'() θ()...72...38636.668.338...72...46825.7679.8232...72...48326.2388.9786..5.72...55724.6973.3269...72...7233.7736.28263.. 3... -.7454.2332 -.332.. 7... -.586.26733 -.6733...72.5..287.63 -.63...72.6..298365.252 -.252...72..5.392337.6535.34694...72...398724.63698.3639...72. 2..48879.677.388227, Biot number;, Grashof number; Pr, Prtl number;, internal heat generation parameter; Ra, radiation parameter. S Afr J Sci 2; 7(9/)
Page 4 of 6.2 =. =.5 + =.8 Pr =.72 + Pr Pr = 7..8.6 f ().6 f ().4.4.2.2 2 4 6 8 2 4 6 8 Prtl number, Pr =.72; internal heat generation parameter, =.; Biot number, =.; radiation parameter, Ra =.. FIGURE : The fluid velocity profile with increasing distance from the plate surface (where f () = ) increasing Grashof number ( Biot number, =.; internal heat generation parameter, =.; Grashof number, =.; radiation parameter, Ra =.5. FIGURE 4: The fluid velocity profile with increasing distance from the plate surface (where f () = ) increasing Prtl number (Pr) when the radiation parameter, Ra =.5..8 =. = 2 + =.7.8 Pr =.72 Pr = 7. + Pr.6.6 f ().4 f ().4.2.2 2 4 6 8 2 4 6 8 Prtl number, Pr =.72; internal heat generation parameter, =.; Grashof number, =.; radiation parameter, Ra =.. FIGURE 2: The fluid velocity profile with increasing distance from the plate surface (where f () = ) increasing Biot number ( Biot number, =.; internal heat generation parameter, =.; Grashof number, =.; radiation parameter, Ra =.. FIGURE 5: The fluid velocity profile with increasing distance from the plate surface (where f () = ) increasing Prtl number (Pr) when the radiation parameter, Ra =...8 =.5 + =.5 =.56 =..8 Ra =. Ra = 2 + Ra = Ra =.5.6.6 f ().4 f ().4.2.2 2 4 6 8 Prtl number, Pr =.72; Biot number, =.; Grashof number, =.; radiation parameter, Ra =.. FIGURE 3: The fluid velocity profile with increasing distance from the plate surface (where f () = ) increasing internal heat generation ( 2 4 6 8 Prtl number, Pr =.72; internal heat generation parameter, =.; Grashof number, =.; Biot number, =.. FIGURE 6: The fluid velocity profile with increasing distance from the plate surface (where f () = ) increasing radiation (Ra S Afr J Sci 2; 7(9/)
Page 5 of 6.3 =. + =.5 =.3 Pr =.72 Pr = 7. + Pr.2.2 θ() θ().. 2 4 6 8 2 4 6 8 Prtl number, Pr =.72; internal heat generation parameter, =.; Biot number, =.; radiation parameter, Ra =.. FIGURE 7: The fluid temperature profile with increasing distance from the plate surface (where θ() = maximum) increasing Grashof number ( Internal heat generation parameter, =.; Grashof number, =.; Biot number, =.; radiation parameter, Ra =.5. FIGURE : The fluid temperature profile with increasing distance from the plate surface (where θ() = maximum) increasing Prtl number (Pr) when the radiation parameter, Ra =.5. θ().9.8.7.6.5.4.3.2. 2 4 6 8 =. + =.7 = 2 θ().2 -.2 -.4 -.6 -.8 - -.2 -.4 -.6 2 6 8 Pr =.72 Pr = 7. + Pr Prtl number, Pr =.72; internal heat generation parameter, =.; local Biot number, Bi =.; radiation parameter, Ra =.. FIGURE 8: The fluid temperature profile with increasing distance from the plate surface (where θ() = maximum) increasing Biot number ( Internal heat generation parameter, =.; Grashof number, =.; Biot number, =.; radiation parameter, Ra =.. FIGURE : The fluid temperature profile with increasing distance from the plate surface (where θ() = maximum) increasing Prtl number (Pr) when the radiation parameter, Ra =...4.3 =.5 =.56 + =.5 =..4.3 Ra =. + Ra = Ra = 2 Ra =.5 θ().2 θ().2.. 2 4 6 8 Prtl number, Pr =.72; Grashof number, =.; Biot number, =.; radiation parameter, Ra =.. FIGURE 9: The fluid temperature profile with increasing distance from the plate surface (where θ() = maximum) increasing internal heat generation ( 2 4 6 8 S Afr J Sci 2; 7(9/) Prtl number, Pr =.72; internal heat generation parameter, =.; Grashof number, =.; Biot number, =.. FIGURE 2: The fluid temperature profile with increasing distance from the plate surface (where θ() = maximum) increasing radiation (Ra
Page 6 of 6 θ().25.2.5..5 convective surface heat transfer, the internal heat generation parameter the radiation parameter enhanced thermal diffusion whilst an increase in the Prtl number the intensity of buoyancy force slowed down the rate of thermal diffusion within the boundary layer. Figure 3 shows the influence of Prtl numbers on the thermal boundary layer, as obtained by Aziz 8. Conclusions 2 4 6 8 Pr =.72 + Pr = 7. Pr = Internal heat generation parameter, = ; Grashof number, = ; Biot number, =.; radiation parameter, Ra =. FIGURE 3: The fluid temperature profile obtained by Aziz 8 with increasing Prtl number (Pr We analysed the effects of internal heat generation, thermal radiation 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 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 ½ if the thermal expansion coefficient β the heat released per unit per mass Q are proportional to x. Numerical solutions of the similarity equations were reported for the various parameters embedded in the problem. The combined effect of increasing the Prtl number 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 the radiation absorption parameter enhanced thermal diffusion. Acknowledgement O.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-7 References. Sakiadis BC. Boundary-layer behavior on continuous solid surfaces. I. Boundary-layer equations for two-dimensional axisymmetric flow. AIChE J. 96;7:26 28. doi:.2/aic.6978 2. Crane IJ. Flow past a stretching plate. Z Angew Math Phys. 97;2(56): 37. 3. Gupta PS, Gupta AS. Heat mass transfer on a stretching sheet with suction blowing. Can J Chem Eng. 977;55:744 746. doi:.2/ cjce.5455569 4. Carragher P, Crane IJ. Heat transfer on a continuous stretching sheet. Z Angew Math Mech. 982;62:564 565. doi:.2/zamm.982629 5. Danberg JE, Fansler KS. A nonsimilar moving wall boundary-layer problem. Quart Appl Math. 976;34:35 39. 6. Chakrabarti A, Gupta AS. Hydromagnetic flow heat transfer over a stretching sheet. Quart Appl Math. 979;37:73 78. 7. Vajravelu K. Hydromagnetic flow heat transfer over a continuous moving porous, flat surface. Acta Mech. 986;64:79 85. doi:.7/ BF45393 8. Dutta BK. Heat transfer from a stretching sheet with uniform suction blowing. Acta Mech. 986;78:255 262. doi:.7/bf7922 9. Lee SL, Tsai JS. Cooling of a continuous moving sheet of finite thickness in the presence of natural convection. Int J Heat Mass Transfer. 99;33:457 464. doi:.6/7-93(9)98-s. Wanous KJ, Sparrow EM. Heat transfer for flow longitudinal to a cylinder with surface mass transfer. J Heat Transf ASME Ser C. 965;87():37 39.. Catherall D, Stewartson K, Williams PG. Viscous flow past a flat plate with uniform injection. Proc R Soc A. 965;284:37 396. doi:.98/ rspa.965.69 2. Sparrow EM, Quack H, Boerner CJ. Local non-similarity boundary layer solutions. J AIAA 97;8():936 942. 3. Sparrow EM, Yu HS. Local non-similarity thermal boundary-layer solutions. J Heat Transf ASME. 97;93:328 334. 4. Massoudi M. Local non-similarity solutions for the flow of a non- Newtonian fluid over a wedge. Int J Non-Linear Mech. 2;36:96 976. doi:.6/s2-7462()6-5 5. Raptis A, Perdikis C, Takhar HS. Effects of thermal radiation on MHD flow. Appl Math Comput. 24;53:645 649. doi:.6/s96-33(3)657-x 6. Hayat T, Abbas Z, Sajid M, Asghar S, The influence of thermal radiation on MHD flow of a second grade fluid. Int J Heat Mass Transfer. 27;5:93 94. doi:.6/j.ijheatmasstransfer.26.8.4 7. Ishak A. Similarity solutions for flow heat transfer over a permeable surface with convective boundary condition. Appl Math Comput. 2;27:837 842. doi:.6/j.amc.2.6.26 8. Aziz A. A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition. Commun Nonlinear Sci Numer Simulat. 29;4:64 68. 9. Aziz A. Hydrodynamic thermal slip flow boundary layers over a flat plate with constant heat flux boundary condition. Commun Nonlinear Sci Numer Simulat. 2;5:573 58. doi:.6/j.cnsns.29.4.26 2. Makinde OD, Olanrewaju PO. Buoyancy effects on thermal boundary layer over a vertical plate with a convective surface boundary condition. J Fluids Eng. 2;32:4452 4455. doi:.5/.4386 2. Chamkha AJ. Hydromagnetic natural convection from an isothermal inclined surface adjacent to a thermally stratified porous medium. Int J Eng Sci. 997;37:975 986. doi:.6/s2-7225(96)22-x 22. Heck A. Introduction to Maple. 3rd ed. New York: Springer-Verlag; 23. S Afr J Sci 2; 7(9/)