Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2009 (2009), Article ID 965603, 13 pages
http://dx.doi.org/10.1155/2009/965603
Research Article

Thermal Radiation Effects on Heat and Mass Transfer over an Unsteady Stretching Surface

1Department of Mathematics, University of Venda, Private Bag X5050, Thohoyandou 0950, South Africa
2Department of Mathematics, University of Swaziland, Private Bag 4, Kwaluseni M201, Swaziland

Received 9 June 2009; Revised 25 September 2009; Accepted 5 November 2009

Academic Editor: Mehrdad Massoudi

Copyright © 2009 Stanford Shateyi and Sandile Sydney Motsa. 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 unsteady heat, mass, and fluid transfer over a horizontal stretching sheet has been numerically investigated. Using a similarity transformation the governing time-dependent boundary layer equations for the momentum, heat, and mass transfer were reduced to a sets of ordinary differential equations. These set of ordinary differential equations were then solved using the Chebyshev pseudo-spectral collocation method, and a parametric analysis was carried out. The study observed, among other observations that the local Sherwood number increases as the values of the stretching parameter ?? and the Schmidt number ???? increase. Also the fluid temperature was found to be significantly reduced by increases in the values of the Prandtl number ????, the unsteadiness parameter ??, and the radiation parameter ??. The velocity and concentration profiles were found to be reduced by increasing values of the unsteadiness parameter ??.

1. Introduction

The continuous surface heat and mass transfer problem has many practical applications in industrial manufacturing processes. The knowledge of flow and heat and mass transfer within a thin liquid film is crucial in understanding the coating process and design of heat exchangers and chemical processing equipments. This phenomenon is also applied in wire and fibre coatings, food stuff processing reactor fluidization, and transpiration cooling. The prime aim in almost every extrusion is to maintain the surface quality of the extrudate. The problem of extrusion of thin surface layers needs special attention to gain some knowledge for controlling the coating efficiently. In the pioneering work of Crane [1], the flow of Newtonian fluid over a linearly stretching surface was studied. Subsequently, the pioneering works of Crane are extended by many authors to explore various aspects of the flow and heat transfer occurring in an infinite domain of the fluid surrounding the stretching sheet (see [28]). However, these studies dealt with a steady flow only. In some cases the flow field and heat transfer can be unsteady due to a sudden stretching of the flat sheet or by a steep change of the temperature of the sheet.

Hossain et al. [9] determined the effect of radiation on natural convection flow of an optically thick viscous incompressible flow past a heated vertical porous plate with a uniform surface temperature and a uniform rate of suction where radiation is included by assuming the Rosseland discussion approximation.

Andersson et al. [10] investigated using a similarity transformation the flow of a thin liquid film of a power-law fluid by unsteady stretching of a surface. Later on, Andersson et al. [11] analyzed the momentum and heat transfer in a laminar liquid film on a horizontal stretching sheet governed by time-dependent boundary layer equations. Dandapat et al. [12] explored how the hydrodynamics and heat transfer in a liquid film on unsteady stretching surface are affected by thermo-capillarity. Tsai et al. [13] studied the non-uniform heat source/sink effect on the flow and heat transfer from an unsteady stretching sheet through a quiescent fluid medium extending to infinity.

Liu and Andersson [14] generalized the analysis by Andersson et al. [11] of the thermal characteristics of a liquid film driven by an unsteady stretching surface. Here, they considered a more general form of the prescribed temperature variation of the stretching sheet than that considered in Andersson et al. [11].

Elbashbeshy and Bazid [15] presented similarity solutions of the boundary layer equations, which describe the unsteady flow and heat transfer over an unsteady stretching sheet. However, as was to be later discovered by Abd El-Aziz [16], some physically unrealistic phenomena are encountered for specific values of the unsteadiness parameter. To that end, Abd El-Aziz [16] was concerned with correcting the similarity transformation introduced by Elbashbeshy and Bazid [15] and extended their analysis to include the effect of thermal radiation. Thermal radiation plays a very significant role in controlling heat transfer in polymer processing industry. The quality of the final product depends to a great extend on the heat controlling factors, and the knowledge of radiative heat transfer in the system can perhaps lead to a desired product with sought qualities.

In this paper, we extend the work of Abd El-Aziz [16] to include mass transfer. This problem arises in a number of industrial manufacturing processes such as polymer extrusion, wire drawing, metal and plastic extrusion, continuous casting, glass fibre production, crystal growing and paper production. The physical and thermal characteristics of such unsteady process are investigated in the boundary layer approximation, assuming variation of the surface temperature and concentration with the horizontal coordinate ?? and time, ??. The main objective of this work is to try to solve this type of model using Chebyshev pseudospectral collocation method. In previous studies finite-difference, Runge-Kutta integration schemes were used. However, it is now known that Chebyshev pseudospectral collocation method gives better approximations than most of these numerical methods. Our main focus is on the comparison of the results obtained by Chebyshev pseudospectral collocation method and those obtained by Abd El-Aziz [16]. We also look at the effects of the unsteadiness parameter on the concentration characteristics of the overlying above the stretching surface.

2. Problem Formulation

We consider the flow of a viscous and incompressible fluid on a horizontal sheet, which comes through a slot at the origin. We are considering a gray, absorbing-emitting radiation fluid but nonscattering medium, and the Rosseland approximation is used to describe the radiative heat flux in the energy equation. The radiative heat flux in the ??-direction is negligible in comparison with that in the ??-direction. The fluid flow over the unsteady stretching sheet is composed of a reacting chemical species. The fluid motion arises due to the stretching of the elastic sheet. The continuous sheet aligned with the ??-axis at ??=0 moves in its own plane with a velocity ????(??,??), the temperature distribution ????(??,??) and the concentration distribution ????(??,??) vary both along the sheet and with time. The velocity, temperature, and concentration fields in the boundary layer are governed by the two-dimensional boundary layer equations for mass, momentum, thermal energy, and chemical species given by:

????+????????????=0,(2.1)????????+??????????+????????????=??2??????2,(2.2)????????+??????????+??????????=??0??2??????2-1????????????????,(2.3)????????+??????????+????????????=??2??????2,(2.4) where ?? and ?? are the velocity components along the ??- and ??-axes, respectively, ?? is the fluid, temperature, ?? is the concentration, ?? is the kinematic coefficient of viscosity, ??0=??/?????? is the thermal diffusivity with ?? as the fluid thermal conductivity, ?? the fluid density, ???? the specific heat at constant pressure, ?? the mass diffusivity, and ???? as the radiative heat flux.

By using the Rosseland diffusion approximation (Hossain et al. [9]) and following Raptis [17] among other researchers, the radiative heat flux, ???? is given by

????=-4??*3??1????4,????(2.5) where ??* and ??1 are the Stefan-Boltzman constant and the Roseland mean absorption coefficient, respectively. We assume that the temperature differences within the flow are sufficiently small such that ??4 may be expressed as a linear function of temperature ??4˜4??38??-3??48.(2.6) Using (2.5) and (2.6) in the last term of (2.3), we obtain

??????????=-16??*??383??1????????2??????2.(2.7) The associated boundary conditions for the model are

??(??,0)=????(??,??),??(??,0)=0,??(??,0)=????(??,??),??(??,0)=????(??,??),(2.8)??(??,8)=0,??(??,8)=??8,??(??,8)=??8.(2.9) Following Andersson et al. [10], the stretching velocity ????(??,??) is assumed ????=????/(1-????), where both ?? and ?? are positive constants with dimension reciprocal time. We have ?? as the initial stretching rate ??/(1-????) is increasing with time. In the context of polymer extrusion, the material properties and in particular the elasticity of the extruded sheet may vary with time even though the sheet is being pulled by a constant force.

With unsteady stretching, however, ??-1 becomes the representative time scale of the resulting unsteady boundary layer problem. We assume both the surface temperature ????(??,??) and the surface concentration ????(??,??) of the stretching sheet to vary with the distance ?? along the sheet and time in the following form:

????(??,??)=??8+??0?????2???(1-????)-3/2??,(2.10)??(??,??)=??8+??0?????2???(1-????)-3/2,(2.11) where ??0 is a heating or cooling reference temperature and ??0 is a positive concentration reference. The equation for the temperature increases (reduces) if ??0 is positive (negative) from ??0 at the leading edge in proportion to ??2 and such that the amount of temperature increase (reduction) along the sheet increases with time. The equation for the concentration ??(??,??) of the sheet also represents a situation in which the sheet species increases if ??0 is positive from ??0 at the leading edge in proportion to ??2 and such that the amount of concentration increase along the sheet increases with time. In order to reduce (2.1)–(2.4) into a set of convenient ordinary differential equations, we introduce the similarity variable ?? and the dimensionless variables ??, ??, and ?? (Abd El-Aziz [16]) as: ??=(????)1/2(1-????)-1/2???????(??),??=???1/2(1-????)-1/2??,??=??8+??0?????2?2??(1-????)-3/2??(??),??=??8+??0?????2?2??(1-????)-3/2??(??),(2.12) where ??(??,??) is the physical stream function which automatically satisfies the continuity equation. The velocity components are then derived from the stream function expression and obtained as

??=????????=???????(??),??=-????????=-(????)1/2(1-????)-1/2??(??).(2.13) Governing equations are then transformed into a set of differential equations and associated boundary conditions as given below:

?????+??????-?????2???-???+??2?????=0,(2.14)(3??+4)?????+3??Pr?????-2???????-2?3??+?????????=0,(2.15)???+Sc?????-2???????-2?3??+???????=0,(2.16) where the prime indicates differentiation with respect to ??, ??=??/?? is the unsteady parameter, ??=(??-??8)/(????-??8) is the non-dimensional temperature, ??=(??-??8)/(????-??8) is the non-dimensional concentration, ??=????1/3??*??38 is the radiation parameter, Pr=??/??0 is the Prandtl number and Sc=??/?? is the Schmidt number. We note that equations (2.15) and (2.16) are equivalent if Sc in equation (2.16) is replaced by 3??Pr/(3??+4). In view of equations (2.12), the boundary conditions from (2.8) to (2.10), transform into

??(0)=0,?????(0)=1,??(0)=??(0)=1,(2.17)?(8)=0,??(8)=0,??(8)=0.(2.18)

3. Method of Solution

To solve the governing equations (2.14)–(2.17), we use the Chebyshev pseudospectral collocation method. The domain in the ?? direction is approximated by [0,??], where ?? is the edge of the boundary layer. We use the algebraic mapping ??=(??+1)??/2 to map the physical region into the Chebyshev spectral method computation domain of [-1,1]. The unknown variables ??,??, and ?? in (2.14)–(2.17) are approximated by the following truncated series of polynomials:

??(??)=?????=0?????????(??),??(??)=?????=0????????(??),??(??)=?????=0?????????(??),(3.1) where the values ?????, ????, and ????? are the coefficients of the series to be determined, ???? is the Chebyshev polynomial of degree ?? which is defined on the interval ???[-1,1] as

?????(??)=cos??cos-1???,??=0,1,.(3.2) We use the Gauss-Lobatto collocation points to define the Chebyshev nodes in [-1,1], namely;

????=cos??????-1=??=1,??=0,1,,??.(3.3) Derivatives of the functions ??,??, and ?? at the collocation points are represented as

????=2???????????=0???????????,????=2???????????=0??????????,????=2???????????=0???????????,(3.4) where ?? is the Chebyshev spectral method differentiation matrix (see, e.g., Canuto et. al [18] for details). Higher order derivatives are computed as multiple powers of ??, that is, (with ?? as an example)

????????????=?2??????????=0?????????????,??=0,1,2,??,(3.5) where ?? is the order of the derivative.

Substituting (3.1)–(3.5) into the governing equations and boundary conditions (2.14)–(2.17), we obtain the following equations:

??3??+????2???-????????????-????????+2??2???=0,(3.6)(3??+4)??2???T+3??Pr??T-2T??????-2?3T+????T??=0,(3.7)??2???F+Sc??F-2F??????-2?3F+????F??=0,(3.8) where ??=(2/??)??, ?????={0,???1???,,??}??, ??T={0,??1??,,??}??, ???F={0,???1???,,??}?? and the superscript ?? denotes Transpose. The boundary conditions become

?????=0,?????=0????????????=1,?????=0???0?????????=0,(3.9)??=???????=1,0=???0=0.(3.10) We note that (3.6) can be solved separately because it depends on ?? only. To solve the system of equations we start by solving (3.6) subject to the boundary conditions (3.9) is a MATLAB nonlinear equation solver fzero which is based on the quasi-Newton method. Once the solution for ?? is obtained, it is then substituted into equation (3.7) to get:

??T=??,thenT=??-1??,(3.11) where ??=(3??+4)??2???+3??Pr??-2??????-2?3??+????[]??,??=00?01??,(3.12) and ?? is an identity matrix of size ??+1. The boundary conditions are imposed by setting the first and last row of B to be [10?00] and [00?01], respectively. The solution for (3.8) is obtained in a similar manner.

4. Results and Discussion

The nonlinear ordinary equations (2.14)–(2.16) (in the case when ??=0, Sc=2, and Pr set to be Pr=2(3??+4)/3??) together with the associated boundary conditions have exact solutions in the form ??(??)=1-??-??,??(??)=??-2??,??(??)=??-2??.(4.1) To validate out numerical method of solution, we compared our numerical results to exact solution (4.1). Table 1 presents the comparison between the exact solution and the Chebyshev spectral method of the skin friction -????(0) and the Nusselt number -??'(0). The table also gives the maximum errors (??????,?? and (??????,??)) which refer to the maximum difference between the approximated and exact values of ?? and ?? at the Gauss-Lobatto collocation points. In computing the numerical results, we used ??=20 and ??=100.

tab1
Table 1: Comparison of the values of -????(0), -???(0) and the maximum error between the exact solution and the Chebyshev spectral method.

From Table 1 it can be seen that the Chebyshev spectral method is in good agreement with the exact solution for the skin friction, Nusselt number, and the values of ?? and ??. Figure 1 shows the numerical solutions for the velocity and temperature profiles compared against the exact solutions (4.1). As depicted in this figure, Chebyshev spectral method gives excellent approximations to the model under consideration.

965603.fig.001
Figure 1: Graph of the velocity ??? and temperature ?? profiles plotted against the exact solutions with ??=1, ??=1, Pr = 0.7, Sc = 0.2.

The present numerical results were also compared with results of El-Aziz [16] (see Table 2).

tab2
Table 2: Comparison of the Chebyshev solution for -???(0) with the results of El-Aziz [16].

In Table 2, we have the data comparisons with the previous published paper and the results are found to be in excellent agreement. It is clearly seen in this table that for given values of ??, the local Nusselt number is increased as the Prandtl number increases. Also, when all other parameters are kept constant, higher local Nusselt number are observed at larger values of Pr. This can be elucidated from the fact that as the Prandtl number increases, the thermal boundary layer thickness decreases and the wall gradient increases. We can also observe in this table that for fixed values of Pr, the local Nusselt number increases as the unsteadiness parameter ?? increases.

Figure 2 illustrates the concentration profile for Sc=0.2 and varying unsteadiness parameters ?? = 0, 0.5, and 2.0. The results show that the concentration ??(??) decreases with increases in the values of the unsteadiness parameter.

965603.fig.002
Figure 2: Effect of various values of unsteadiness parameter ?? on ?? at Sc=0.2, ??=1, Pr=0.7.

In Figure 3, we have the local Sherwood number in terms of -???(0) as a function of the unsteadiness parameter ?? for some typical Schmidt number values. It can be clearly seen from this figure that the Sherwood number is increased as the values of the Schmidt number increase. We can also see that, for a given value of ??, higher Sherwood number values are observed at large values of Sc. This can be explained from the fact that as the Schmidt increases, the concentration boundary layer thickness decreases and the wall concentration gradient increases.

965603.fig.003
Figure 3: Local Sherwood number -??'(0) versus ?? for various values of Sc with ??=1, Pr = 0.7.

In Figure 4, we show the effect of the Schmidt number on the concentration distributions. The Schmidt number represents the relative ease of molecular momentum and mass transfer and is very important in calculations of binary mass transfer in multiphase flows. The effect of an increase in the Schmidt number values is to reduce the momentum boundary layer and this leads to the thinning of the diffusion layer. Diffusing chemical species of most interest in air have Schmidt numbers in the range from 0.1 to 10 (Hussain and Hossain [19]).

965603.fig.004
Figure 4: Effect of varying the Schmidt number Sc on the concentration profiles with Pr=0.7, ??=1, ??=1.

In Figure 5, we elucidate the effect of the unsteadiness parameter ?? on the typical temperature profiles. At fixed values of other parameters, increasing ?? values, decreases the temperature with an accompanying decreases in the boundary layer thickness.

965603.fig.005
Figure 5: Effect of various values of unsteadiness parameter ?? on on the temperature profiles with Pr=0.7, Sc=0.2, ??=1.

The effects of the Prandtl number on the temperature distributions in the boundary layer are illustrated in Figure 6. We observe that the thermal boundary layer is reduced by increasing the Prandtl number values thereby reducing the fluid temperature. The effects of thermal radiation parameter ?? on the temperature profiles in the boundary layer are illustrated in Figure 7. We observe in this figure that increasing the thermal radiation parameter produces significant decrease in the thermal condition of fluid and its boundary layer. This can be explained by the fact that a decrease in the values of ??(=????1/3????38) for the given values of ?? and ??8 means a decrease in the Rosseland radiation absorptivity ??1. We can deduce in equations (2.3) and (2.5) that the divergence of the radiative heat flux ??????/???? increases as ??1 decreases the rate of radiative heat transferred to the fluid and consequently the fluid temperature increases. It is shown in Figure 7 the effect of thermal radiation becomes more significant as ???0 though not equals to zero and ??=8 corresponds to non-radiating case so the effects of radiation can be neglected when ???8.

965603.fig.006
Figure 6: The variation of the temperature distribution profiles with increasing Prandtl number.
965603.fig.007
Figure 7: The variation of the temperature distribution profiles with increasing thermal radiation parameter.

In Figure 8, we present velocity profiles when all other parameters are kept constant for different values of the unsteadiness parameter ??. We observe in this figure that the velocity decreases with the distance from the stretching sheet for all ??. Also, increasing the values of ?? decreases the velocity in the boundary layer. When ??=0, we have a steady state flow and for ??>0, we have an unsteady flow. As can be seen in Figures 2, 5, and 8, increasing the unsteadiness parameter ?? reduces the flow properties such as velocity, temperature and concentration. When ?? values are increased in the system, the boundary layer thicknesses are reduced and this inhibits the development of transition of laminar to turbulent flow. This shows that stretching of surfaces can be used as a flow stabilizing mechanism.

965603.fig.008
Figure 8: Effect of various values of unsteadiness parameter ?? on on the velocity profiles.

5. Conclusion

The unsteady heat, mass, and fluid transfer over a horizontal stretching sheet has been investigated in this study. The study put much emphasis on the numerical method used to solve the set of non-linear ordinary differential equations and on how the concentration distribution and mass transfer change due to stretching of the sheet. The study found out that the Nusselt number increases with the increase of the unsteadiness parameter and Prandtl number. The fluid temperature was found to be reduced by increases in the values of the unsteadiness parameter, Prandtl number and the radiation parameter. We also observed that the velocity profiles and the concentration distributions were reduced as the values of the unsteadiness parameter ?? were increased. The local Sherwood number was found to be increased by the increased values of ?? and the Schmidt number.

References

  1. L. J. Crane, “Flow past a stretching plate,” Zeitschrift für angewandte Mathematik und Physik, vol. 12, pp. 645–647, 1970. View at Google Scholar
  2. M. K. Laha, P. S. Gupta, and A. S. Gupta, “Heat transfer characteristics of the flow of an incompressible viscous fluid over a stretching sheet,” Warme und Stoffubertrag, vol. 24, no. 3, pp. 151–153, 1989. View at Publisher · View at Google Scholar
  3. 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 Publisher · View at Google Scholar · View at Zentralblatt MATH
  4. K. V. Prasad, A. S. Abel, and P. S. Datti, “Diffusion of chemically reactive species of non-Newtonian fluid immersed in a porous medium over a stretching sheet,” The International Journal of Non-Linear Mechanics, vol. 38, no. 5, pp. 651–657, 2003. View at Publisher · View at Google Scholar
  5. M. S. Abel and N. Mahesha, “Heat transfer in MHD viscoelastic fluid flow over a stretching sheet with variable thermal conductivity, non-uniform heat source and radiation,” Applied Mathematical Modelling, vol. 32, no. 10, pp. 1965–1983, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  6. P. G. Abel, Siddheshwar, and M. M. Nandeppanava, “Heat transfer in a viscoelastic boundary layer flow over a stretching sheet with viscous dissipation and non-uniform heat source,” International Journal of Heat and Mass Transfer, vol. 50, no. 5-6, pp. 960–966, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  7. R. Cortell, “Viscoelastic fluid flow and heat transfer over a stretching sheet under the effects of a non-uniform heat source, viscous dissipation and thermal radiation,” International Journal of Heat and Mass Transfer, vol. 50, no. 15-16, pp. 3152–3162, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  8. C. Wang, “Analytic solutions for a liquid film on an unsteady stretching surface,” Heat and Mass Transfer, vol. 42, pp. 43–48, 2006. View at Google Scholar
  9. M. A. Hossain, M. A. Alim, and D. A. S. Rees, “The effect of radiation on free convection from a porous vertical plate,” International Journal of Heat and Mass Transfer, vol. 42, no. 1, pp. 181–191, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  10. H. I. Andersson, J. B. Aarseth, N. Braud, and B. S. Dandapat, “Flow of a power-law fluid film on an unsteady stretching surface,” Journal of Non-Newtonian Fluid Mechanics, vol. 62, no. 1, pp. 1–8, 1996. View at Publisher · View at Google Scholar
  11. H. I. Andersson, J. B. Aarseth, and B. S. Dandapat, “Heat transfer in a liquid film on unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 43, no. 1, pp. 69–74, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  12. B. S. Dandapat, B. Santra, and H. I. Andersson, “Thermocapillary in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 46, pp. 3009–3015, 2003. View at Google Scholar
  13. R. Tsai, K. H. Huang, and J. S. Huang, “Flow and heat transfer over an unsteady stretching surface with non-uniform heat source,” International Communications in Heat and Mass Transfer, vol. 35, no. 10, pp. 1340–1343, 2008. View at Publisher · View at Google Scholar
  14. I.-C. Liu and H. I. Andersson, “Heat transfer in a liquid film on an unsteady stretching sheet,” International Journal of Thermal Sciences, vol. 47, no. 6, pp. 766–772, 2008. View at Publisher · View at Google Scholar
  15. E. M. A. Elbashbeshy and M. A. A. Bazid, “Heat transfer over an unsteady stretching surface with internal heat generation,” Applied Mathematics and Computation, vol. 138, no. 2-3, pp. 239–245, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  16. M. Abd El-Aziz, “Radiation effect on the flow and heat transfer over an unsteady stretching sheet,” International Communications in Heat and Mass Transfer, vol. 36, no. 5, pp. 521–524, 2009. View at Publisher · View at Google Scholar
  17. A. Raptis, “Flow of a micropolar fluid past a continuously moving plate by the presence of radiation,” International Journal of Heat and Mass Transfer, vol. 41, no. 18, pp. 2865–2866, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  18. C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. Zang, Spectral Method in Fluid Dynamics, Springer, New York, NY, USA, 1998.
  19. S. Hussain and M. A. Hossain, “Natural convection flow from a vertical permeable flat plate with variable surface temperature and species concentration,” Engineering Computations, vol. 17, no. 7, pp. 789–813, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH