About this Journal Submit a Manuscript Table of Contents 
ISRN Thermodynamics
Volume 2012 (2012), Article ID 795814, 6 pages
doi:10.5402/2012/795814
Research Article

Heat Generation and Thermal Radiation Effects over a Stretching Sheet in a Micropolar Fluid

M. Gnaneswara Reddy

Department of Mathematics, Acharya Nagarjuna University, Ongole Campus, Ongole 523 001, India

Received 6 July 2012; Accepted 26 July 2012

Academic Editors: C. D. Daub and S. Hashimoto

Copyright © 2012 M. Gnaneswara Reddy. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The effects of radiation and heat generation on steady thermal boundary layer flow induced by a linearly stretching sheet immersed in an incompressible micropolar fluid with constant surface temperature are investigated. Similarity transformation is employed to transform the governing partial differential equations into ordinary ones, which are then solved numerically using the Runge-Kutta fourth order along shooting method. Results for the local Nusselt number as well as the temperature profiles are presented for different values of the governing parameters. It is observed that the velocity increases with an increase in the material parameter. It is seen that the temperature profile is influenced considerably and increases when the value of heat generation parameter increases along the boundary layer. Also, the temperature distribution of the fluid increases with an increase in the radiation parameter. Comparisons with previously published work are performed and the results are found to be in very good agreement.

1. Introduction

Flow of a viscous fluid past a stretching sheet is a classical problem in fluid dynamics. The development of boundary layer flow induced solely by a stretching sheet was first studied by Crane [1] who first obtained an elegant analytical solution to the boundary layer equations for the problem of steady two dimensional flow due to a stretching surface in a quiescent incompressible fluid. Flow and heat transfer characteristics due to a stretching sheet in a stationary fluid occur in a number of industrial manufacturing processes and include both metal and polymer sheets, for example, the cooling of an infinite metallic plate in a cooling bath, the boundary layer along material handling conveyers, the aerodynamic extrusion of plastic sheets, paper production, metal spinning, and drawing plastic films. The quality of the final product depends on the rate of heat transfer at the stretching surface. This problem was then extended by P. S. Gupta and A. S. Gupta [2] to a permeable surface. The flow problem due to a linearly stretching sheet belongs to a class of exact solutions of the Navier-Stokes equations. Thus, the exact solutions reported by Crane [1] and P. S. Gupta and A. S. Gupta [2] are also the exact solutions to the Navier-Stokes equations. The heat transfer aspects of similar problems were studied by Grubka and Bobba [3], Chen and Char [4], Dutta et al. [5], Ali [6, 7], Afzal and Varshney [8], Afzal [9], and many others. On the other hand, the effects of buoyancy force on the development of velocity and thermal boundary layer flows over a stretching sheet have been investigated by Chen [10], Ali and Al-Yousef [11], Daskalakis [12], Partha et al. [13], Abd El-Aziz [14], Mahapatra et al. [15], and Ishak et al. [16–20], among others.

The study of flow and heat transfer past a stretching sheet has gained tremendous interest among researchers due to its industrial and engineering applications. This includes extrusion of plastic sheets, annealing and tinning of copper wire, paper production, crystal growing, and glass blowing. The final products depend mainly on the stretching and cooling rates at the surface. Their studies are not only restricted to the Newtonian fluids but also include the non-Newtonian fluids such as micropolar fluids. Such studies have been carried out by Chiam [21], Heruska et al. [22], Agarwal et al. [23], Hassanian and Gorla [24], Kelson and Desseaux [25], Kelson and Farrell [26], Nazar et al. [27], and very recently by Hayat et al. [28].

In the present paper we study the development of thermal boundary layer flow induced by a stretching sheet immersed in a micropolar fluid with the effect of heat generation and thermal radiation is taken into consideration. The governing partial differential equations are transformed into ordinary ones using similarity transformation, before being solved numerically by the Runge-Kutta fourth order along shooting method.

2. Mathematical Model

A steady laminar boundary layer flow over a stretching sheet immersed in a quiescent and incompressible micropolar fluid with uniform surface temperature 𝑇 𝑤 is considered. It is assumed that the sheet is stretched with a linear velocity 𝑈 𝑤 = 𝑎 𝑥 , where 𝑎 is a positive constant and 𝑥 is the distance from the slit where the sheet is issued. The simplified two-dimensional equations governing the flow may be written as

Continuity equation 𝜕 𝑢 + 𝜕 𝑥 𝜕 𝑣 𝜕 𝑦 = 0 . ( 1 )

Momentum equation 𝑢 𝜕 𝑢 𝜕 𝑥 + 𝑣 𝜕 𝑢 = 𝜅 𝜕 𝑦 𝜈 + 𝜌 𝜕 2 𝑢 𝜕 𝑦 2 + 𝜅 𝜌 𝜕 𝑁 𝜕 𝑦 . ( 2 )

Angular Momentum equation, 𝑢 𝜕 𝑁 𝜕 𝑥 + 𝑣 𝜕 𝑁 = 𝛾 𝜕 𝑦 𝜕 𝜌 𝑗 2 𝑁 𝜕 𝑦 2 − 𝜅 𝜌 𝑗 2 𝑁 + 𝜕 𝑢 𝜕 𝑦 . ( 3 )

Energy equation, 𝑢 𝜕 𝑇 𝜕 𝑥 + 𝑣 𝜕 𝑇 = 𝑘 𝜕 𝑦 𝜌 𝑐 𝑝 𝜕 2 𝑇 𝜕 𝑦 2 − 1 𝜌 𝑐 𝑝 𝜕 𝑞 𝑟 + 𝑄 𝜕 𝑦 0 𝜌 𝑐 𝑝 𝑇 − 𝑇 ∞ . ( 4 )

The spin-gradient viscosity 𝛾 (Ahmadi [29]) can be defined as 𝜅 𝛾 = 𝜇 + 2 𝐾 𝑗 = 𝜇 1 + 2 𝑗 , ( 5 ) where 𝜇 is the dynamic viscosity, 𝐾 = 𝜅 / 𝜇 is the dimensionless viscosity ratio and is called the material parameter, and we take 𝑗 = 𝜈 / 𝑎 as a reference length. The Equation (5) is invoked to allow the field of equations predicts the correct behavior in the limiting case when the microstructure effects become negligible and the total spin 𝑁 reduces to the angular velocity.

The boundary conditions for the velocity, angular velocity, and fields are 𝑢 = 𝑈 𝑤 , 𝑣 = 0 , 𝑁 = − 𝑚 𝜕 𝑢 , 𝜕 𝑦 𝑇 = 𝑇 𝑤 a t 𝑦 = 0 𝑢 ⟶ 0 , 𝑁 ⟶ 0 , 𝑇 ⟶ 𝑇 ∞ a s 𝑦 ⟶ ∞ , ( 6 ) where 𝑇 ∞ is the ambient fluid temperature and 𝑚 is the boundary parameter with 0 ≤ 𝑚 ≤ 1 .

By using the Rosseland approximation, Brewster [30], the radiative heat flux 𝑞 𝑟 is given by 𝑞 𝑟 = − 4 𝜎 𝑠 3 𝑘 𝑒 𝜕 𝑇 4 , 𝜕 𝑦 ( 7 ) where 𝜎 𝑠 is the Stefan-Boltzmann constant and 𝑘 𝑒 is the mean absorption coefficient. It should be noted that by using the Rosseland approximation, the present analysis is limited to optically thick fluids. If the temperature differences within the flow are sufficiently small, then (7) can be linearized by expanding 𝑇 4 into the Taylor series about 𝑇 ∞ , which after neglecting higher order terms takes the form 𝑇 4 ≅ 4 𝑇 3 ∞ 𝑇 − 3 𝑇 4 ∞ . ( 8 ) In view of (7) and (8), (4) reduces to 𝑢 𝜕 𝑇 𝜕 𝑥 + 𝑣 𝜕 𝑇 𝜕 𝜕 𝑦 = 𝛼 ( 1 + 𝑅 ) 2 𝑇 𝜕 𝑦 2 + 𝑄 0 𝜌 𝑐 𝑝 𝑇 − 𝑇 ∞ , ( 9 ) where 𝛼 = 𝑘 / 𝜌 𝑐 𝑝 is the thermal diffusivity and 𝑅 = 1 6 𝜎 ∗ 𝑇 3 ∞ / 3 𝑘 𝑘 ∗ is the radiation parameter.

We introduce now the following similarity transformation: 𝑈 𝜂 = 𝑤 𝜈 𝑥 1 / 2 𝑦 , 𝜓 = 𝜈 𝑥 𝑈 𝑤 1 / 2 𝑓 ( 𝜂 ) , 𝑁 = 𝑈 𝑤 𝑈 𝑤 𝜈 𝑥 1 / 2 ℎ ( 𝜂 ) , 𝜃 ( 𝜂 ) = 𝑇 − 𝑇 ∞ 𝑇 𝑤 − 𝑇 ∞ , 𝑄 𝑄 = 0 𝜈 𝜌 𝑐 𝑝 𝑈 2 𝑤 , ( 1 0 ) where 𝜂 is the similarity variable and 𝜓 is the stream function defined as 𝑢 = 𝜕 𝜓 / 𝜕 𝑦 , 𝑣 = − 𝜕 𝜓 / 𝜕 𝑥 which identically satisfies the mass conservation equation (1). Substituting (10) into (2), (3), and (9) we obtain the following ordinary differential equations: ( 1 + 𝐾 ) 𝑓 + 𝑓 𝑓 − 𝑓 2 𝐾 + 𝐾 ℎ ′ = 0 , ( 1 1 ) 1 + 2 ℎ + 𝑓 ℎ − 𝑓 ℎ − 𝐾 2 ℎ + 𝑓 = 0 , ( 1 2 ) ( 1 + 𝑅 ) 𝜃 + P r 𝑓 𝜃 + P r 𝑄 𝜃 = 0 , ( 1 3 ) where primes denote differentiation with respect to 𝜂 and P r = 𝜇 / 𝛼 is the Prandtl number.

The corresponding dimensionless boundary conditions are 𝑓 = 0 , 𝑓 = 1 , ℎ = − 𝑚 𝑓 𝑓 ( 0 ) , 𝜃 = 1 a t 𝜂 = 0 ⟶ 0 , ℎ ⟶ 0 , 𝜃 ⟶ 0 a s 𝜂 ⟶ ∞ . ( 1 4 )

It is note that K = 0 corresponds to viscous fluid.

3. Results and Discussion

The set of nonlinear ordinary differential equations (11)–(13) with boundary conditions (14) have been solved by using the Runge-Kutta fourth order along with Shooting method. First of all, higher order nonlinear differential Equations (11)–(13) are converted into simultaneous linear differential equations of first order and they are further transformed into initial value problem by applying the shooting technique (Jain et al. [31]). The resultant initial value problem is solved by employing Runge-Kutta fourth order technique. The step size Δ 𝜂 = 0 . 0 1 is used to obtain the numerical solution with five decimal place accuracy as the criterion of convergence. To analyze the results, numerical computation has been carried out using the method described in the previous section for variations in the governing parameters, namely, material parameter 𝐾 , radiation parameter 𝑅 , heat generation parameter 𝑄 , the Prandtl number P r , and boundary parameter 𝑚 . In the present study following default parameter values are adopted for computations: 𝐾 = 1 . 0 , 𝑅 = 1 . 0 , 𝑄 = 1 . 0 , P r = 0 . 7 1 , and 𝑚 = 0 . 5 . All graphs therefore correspond to these values unless specifically indicated on the appropriate graph.

In order to assess the accuracy of our computed results, the present result has been compared with Ishak [32] for different values of 𝐾 as shown in Figure 1 with 𝑄 = 0 . 0 . It is observed that the agreements with the solution of velocity profiles are excellent.

795814.fig.001
Figure 1: Comparison of temperature profiles.

Figure 2 presents the velocity profiles for various values of 𝐾 when 𝑚 = 0 . 5 . We note that the parameters 𝑅 and P r have no influence on the flow field, which is clear from (11)–(13). It is evident from this figure that the boundary layer thickness increases with 𝐾 . The velocity gradient at the surface 𝑓 ( 0 ) decreases (in absolute sense) as 𝐾 increases. Thus, micropolar fluids show drag reduction compared to viscous fluids. The negative velocity gradient at the surface 𝑓 ( 0 ) < 0 as shown in Figure 2 means the stretching sheet exerts a drag force on the fluid. This is not surprising since the development of the boundary layer is solely induced by it. The effect of material parameter 𝐾 on the temperature is shown in Figure 3. The temperature decreases with an increase in the material parameter.

795814.fig.002
Figure 2: Velocity profiles for different values of 𝐾 .
795814.fig.003
Figure 3: Temperature profiles for different values of 𝐾 .

Figure 4 shows the temperature profiles for different values of 𝑅 . The radiation parameter 𝑅 defines the relative contribution of conduction heat transfer to thermal radiation transfer. It is obvious that an increase in the radiation parameter results in increasing temperature within the boundary layer. The effect of heat generation parameter 𝑄 on the temperature is shown in Figure 5. From this figure, we observe that when the value of heat generation parameter increases, the temperature distribution also increases along the boundary layer.

795814.fig.004
Figure 4: Temperature profiles for different values of 𝑅 .
795814.fig.005
Figure 5: Temperature profiles for different values of 𝑄 .

The influence of the Prandtl number P r on temperature field is shown in Figure 6. The numerical results show that the effect of increasing values of the Prandtl number results in a decreasing temperature. It is observed that an increase in the Prandtl number results in a decrease of the thermal boundary layer thickness and in general lower average temperature within the boundary layer. The reason is that smaller values of Pr are equivalent to increasing the thermal conductivities, and therefore heat is able to diffuse away from the heated surface more rapidly than for higher values of Pr. Hence in the case of the smaller Prandtl numbers as the boundary layer is thicker and the rate of heat transfer is reduced.

795814.fig.006
Figure 6: Temperature profiles for different values of P r .

The effects of 𝑚 on velocity, angular velocity, and temperature profiles are depicted in Figures 7, 8, and 9, respectively Figure 7 shows that the velocity gradient at the surface is larger for larger values of 𝑚 . Different behaviors are observed for the effect of 𝑚 on the heat transfer rate at the surface as presented in Figure 9. As expected, the couple stress ℎ ( 0 ) is more dominant for larger values of m, as shown in Figure 8.

795814.fig.007
Figure 7: Velocity profiles for different values of 𝑚 .
795814.fig.008
Figure 8: Angular velocity profiles for different values of 𝑚 .
795814.fig.009
Figure 9: Temperature profiles for different values of 𝑚 .

4. Conclusions

A steady two-dimensional laminar and heat transfer due to a stretching sheet immersed in an incompressible micropolar fluid has been investigated. The effect of thermal radiation and heat generation on the development of the thermal boundary layer flow has been taken into consideration. The effects of the governing parameters 𝐾 , 𝑅 , 𝑄 , P r , and 𝑚 on the fluid flow and heat transfer characteristics are discussed. It is found that the temperature distribution increases as the radiation parameter 𝑅 increases. It is observed that the velocity gradient at the surface is larger for larger values of m.

Nomenclature

𝑎 , 𝑏 : Constants
𝑐 𝑝 :Specific heat at constant pressure
𝑓 :Dimensionless stream function
ℎ :Dimensionless microrotation
𝑗 :Microinertia density
𝑘 :Thermal conductivity
𝑘 ∗ :Mean absorption coefficient
𝐾 :Material parameter
𝑚 :Boundary parameter
𝑁 :Microrotation or angular velocity
P r :Prandtl number
𝑞 𝑟 :Radiative heat flux
𝑄 0 :Heat generation constant
𝑄 :Heat generation parameter
𝑅 :Radiation parameter
𝑇 :Fluid temperature
𝑇 𝑤 :Surface temperature
𝑇 ∞ :Ambient temperature
𝑢 , 𝑣 :Velocity components in the x- and y-directions, respectively
𝑈 𝑤 :Velocity of the stretching sheet
𝑥 , 𝑦 : Cartesian coordinates along the sheet andnormal to it, respectively.
Greek Letters
𝛼 Thermal diffusivity
𝛽 :Thermal expansion coefficient
𝛾 :Spin gradient viscosity
𝜂 :Similarity variable
𝜃 :Dimensionless temperature
𝜅 :Vortex viscosity
𝜈 :Kinematic viscosity
𝜇 :Dynamic viscosity
𝜌 :Fluid density
𝜎 ∗ : Stefan-Boltzmann constant
𝜓 :Stream function.
Subscripts
𝑤 :Condition at the solid surface
∞ :Ambient condition.
Superscript
( ) :Differentiation with respect to 𝜂 .

Acknowledgments

The author is very thankful to the editor and reviewers for their encouraging comments and constructive suggestions to improve the presentation of this paper.

References

  1. L. J. Crane, “Flow past a stretching plate,” Zeitschrift für Angewandte Mathematik und Physik, vol. 21, no. 4, pp. 645–647, 1970. View at Publisher · View at Google Scholar · View at Scopus
  2. P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction or blowing,” The Canadian Journal of Chemical Engineering, vol. 55, no. 6, pp. 744–746, 1977. View at Publisher · View at Google Scholar
  3. L. J. Grubka and K. M. Bobba, “Heat transfer characteristics of a continuous, stretching surface with variable temperature,” Journal of Heat Transfer, vol. 107, no. 1, pp. 248–250, 1985. View at Scopus
  4. C. K. Chen and M. I. Char, “Heat transfer of a continuous, stretching surface with suction or blowing,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 568–580, 1988. View at Scopus
  5. B. K. Dutta, P. Roy, and A. S. Gupta, “Temperature field in flow over a stretching sheet with uniform heat flux,” International Communications in Heat and Mass Transfer, vol. 12, no. 1, pp. 89–94, 1985. View at Scopus
  6. M. E. Ali, “Heat transfer characteristics of a continuous stretching surface,” Heat and Mass Transfer, vol. 29, no. 4, pp. 227–234, 1994. View at Publisher · View at Google Scholar · View at Scopus
  7. M. E. Ali, “On thermal boundary layer on a power-law stretched surface with suction or injection,” International Journal of Heat and Fluid Flow, vol. 16, no. 4, pp. 280–290, 1995. View at Scopus
  8. N. Afzal and I. S. Varshney, “The cooling of a low heat resistance stretching sheet moving through a fluid,” Wärme- und Stoffübertragung, vol. 14, no. 4, pp. 289–293, 1980. View at Publisher · View at Google Scholar · View at Scopus
  9. N. Afzal, “Heat transfer from a stretching surface,” International Journal of Heat and Mass Transfer, vol. 36, no. 4, pp. 1128–1131, 1993. View at Scopus
  10. C. H. Chen, “Laminar mixed convection adjacent to vertical, continuously stretching sheets,” Heat and Mass Transfer, vol. 33, no. 5-6, pp. 471–476, 1998. View at Scopus
  11. M. Ali and F. Al-Yousef, “Laminar mixed convection from a continuously moving vertical surface with suction or injection,” Heat and Mass Transfer, vol. 33, no. 4, pp. 301–306, 1998. View at Scopus
  12. J. E. Daskalakis, “Free convection effects in the boundary layer along a vertically stretching flat surface,” Canadian Journal of Physics, vol. 70, pp. 1253–1260, 1993.
  13. M. K. Partha, P. V. S. N. Murthy, and G. P. Rajasekhar, “Effect of viscous dissipation on the mixed convection heat transfer from an exponentially stretching surface,” Heat and Mass Transfer, vol. 41, no. 4, pp. 360–366, 2005. View at Publisher · View at Google Scholar · View at Scopus
  14. M. Abd El-Aziz, “Temperature dependent viscosity and thermal conductivity effects on combined heat and mass transfer in MHD three-dimensional flow over a stretching surface with Ohmic heating,” Meccanica, vol. 42, no. 4, pp. 375–386, 2007. View at Publisher · View at Google Scholar · View at Scopus
  15. T. R. Mahapatra, S. Dholey, and A. S. Gupta, “Momentum and heat transfer in the magnetohydrodynamic stagnation-point flow of a viscoelastic fluid toward a stretching surface,” Meccanica, vol. 42, no. 3, pp. 263–272, 2006. View at Publisher · View at Google Scholar · View at Scopus
  16. A. Ishak, R. Nazar, and I. Pop, “Unsteady mixed convection boundary layer flow due to a stretching vertical surface,” Arabian Journal for Science and Engineering, vol. 31, no. 2, pp. 165–182, 2006. View at Scopus
  17. A. Ishak, R. Nazar, and I. Pop, “Magnetohydrodynamic stagnation point flow towards a stretching vertical sheet,” Magnetohydrodynamics, vol. 43, no. 1, pp. 83–98, 2006. View at Publisher · View at Google Scholar · View at Scopus
  18. A. Ishak, R. Nazar, and I. Pop, “Mixed convection boundary layers in the stagnation-point flow toward a stretching vertical sheet,” Meccanica, vol. 41, no. 5, pp. 509–518, 2006. View at Publisher · View at Google Scholar · View at Scopus
  19. A. Ishak, R. Nazar, and I. Pop, “Mixed convection on the stagnation point flow toward a vertical, continuously stretching sheet,” Journal of Heat Transfer, vol. 129, no. 8, pp. 1087–1090, 2007. View at Publisher · View at Google Scholar · View at Scopus
  20. A. Ishak, R. Nazar, and I. Pop, “Hydromagnetic flow and heat transfer adjacent to a stretching vertical sheet,” Heat and Mass Transfer, vol. 44, no. 8, pp. 921–927, 2008. View at Publisher · View at Google Scholar · View at Scopus
  21. T. C. Chiam, “Micropolar fluid flow over a stretching sheet,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 62, pp. 565–568, 1982.
  22. M. W. Heruska, L. T. Watson, and K. K. Sankara, “Micropolar flow past a porous stretching sheet,” Computers and Fluids, vol. 14, no. 2, pp. 117–129, 1986. View at Scopus
  23. R. S. Agarwal, R. Bhargava, and A. V. S. Balaji, “Finite element solution of flow and heat transfer of a micropolar fluid over a stretching sheet,” International Journal of Engineering Science, vol. 27, no. 11, pp. 1421–1428, 1989. View at Scopus
  24. I. A. Hassanien and R. S. R. Gorla, “Heat transfer to a micropolar fluid from a non-isothermal stretching sheet with suction and blowing,” Acta Mechanica, vol. 84, no. 1–4, pp. 191–199, 1990. View at Publisher · View at Google Scholar · View at Scopus
  25. N. A. Kelson and A. Desseaux, “Effect of surface conditions on flow of a micropolar fluid driven by a porous stretching sheet,” International Journal of Engineering Science, vol. 39, no. 16, pp. 1881–1897, 2001. View at Publisher · View at Google Scholar · View at Scopus
  26. N. A. Kelson and T. W. Farrell, “Micropolar flow over a porous stretching sheet with strong suction or injection,” International Communications in Heat and Mass Transfer, vol. 28, no. 4, pp. 479–488, 2001. View at Publisher · View at Google Scholar · View at Scopus
  27. R. Nazar, N. Amin, D. Filip, and I. Pop, “Stagnation point flow of a micropolar fluid towards a stretching sheet,” International Journal of Non-Linear Mechanics, vol. 39, no. 7, pp. 1227–1235, 2004. View at Publisher · View at Google Scholar · View at Scopus
  28. T. Hayat, Z. Abbas, and T. Javed, “Mixed convection flow of a micropolar fluid over a non-linearly stretching sheet,” Physics Letters A, vol. 372, no. 5, pp. 637–647, 2008. View at Publisher · View at Google Scholar · View at Scopus
  29. G. Ahmadi, “Self-similar solution of imcompressible micropolar boundary layer flow over a semi-infinite plate,” International Journal of Engineering Science, vol. 14, no. 7, pp. 639–646, 1976. View at Scopus
  30. M. Q. Brewster, Thermal Radiative Transfer Properties, Wiley, New York, NY, USA, 1992.
  31. M. K. Jain, S. R. K. Iyengar, and R. K. Jain, Numerical Methods for Scientific and Engineering Computation, Wiley Eastern, New Delhi, India, 1985.
  32. A. Ishak, “Thermal boundary layer flow over a stretching sheet in a micropolar fluid with radiation effect,” Meccanica, vol. 45, no. 3, pp. 367–373, 2010. View at Publisher · View at Google Scholar · View at Scopus