Research Article  Open Access
Effects of Thermocapillarity and Thermal Radiation on Flow and Heat Transfer in a Thin Liquid Film on an Unsteady Stretching Sheet
Abstract
This paper examines the effects of thermocapillarity and thermal radiation on the boundary layer flow and heat transfer in a thin film on an unsteady stretching sheet with nonuniform heat source/sink. The governing partial differential equations are converted into ordinary differential equations by a similarity transformation and then are solved by using the homotopy analysis method (HAM). The effects of the radiation parameter, the thermocapillarity number, and the temperaturedependent parameter in this study are discussed and presented graphically via velocity and temperature profiles.
1. Introduction
The analysis of heat transfer of boundary layer flow with thermal radiation is important in electrical power generation, astrophysical flows, solar power technology, space vehicle reentry, and other industrial engineering processes.
Wang [1] was the pioneer in investigating the hydrodynamics of a flow in a thin liquid film on an unsteady stretching surface. Later, Andersson et al. [2] studied the heat transfer characteristics of the hydrodynamical problem solved by Wang [1]. Liu and Andersson [3] examined the problem with a more general form of prescribed temperature variation of a stretching sheet. Wang [4] investigated the same problem of Andersson et al. [2], presenting analytic solutions. Several researchers have extended Wang’s [1] classical problem, taking into consideration nonNewtonian [5–9], thermocapillarity [10–13], and magnetic effects [14, 15].
Wang [4] was probably the first to analyze the flow and heat transfer in a thin liquid film on an unsteady stretching surface using the homotopy analysis method (HAM). The HAM is a general analytic method for obtaining series solutions for various types of nonlinear equations [16]. Several studies have successfully applied HAM to various nonlinear problems in science and engineering [17–22]. Studies by Dandapat et al. [10, 11] reported that thermocapillarity generates surfacetension gradients along the horizontal interface between the passive gas and the liquid film. Both Abd ElAziz [23] and Mahmoud [24] have studied the influence of thermal radiation on the flow and heat transfer on an unsteady stretching sheet.
Motivated by these studies, in this paper we extend the results of Wang [4] and Mahmoud [24] for thin film flow on an unsteady stretching sheet with combined effects of thermocapillarity, thermal radiation, and internal heating. Furthermore, we employ the similarity transformation introduced by Wang [4] to transform the extent of the independent variable into a finite range 01. The solutions reached using HAM are presented and implications discussed.
2. Problem Formulation
The fluid flow, modeled as an unsteady, twodimensional, incompressible viscous laminar flow on a horizontal thin elastic sheet, emerges from a narrow slot at the origin of a Cartesian coordinate system. Fluid motion and heat transfer arise in the stretching of the horizontal elastic sheet. It is assumed that the elastic sheet has internal heat generation/absorption and that the sheet temperature varies with the coordinate and time . Under these assumptions, the governing conservation equations of mass, momentum, and energy at unsteady state can be expressed as subject to where and are the velocity components of the fluid in the  and directions, is the time, is the temperature, is the kinematic viscosity, is the density, is the thermal diffusivity, is the specific heat at constant pressure, is the viscosity, and is the uniform thickness of the liquid film. Note that the surface tension varies linearly with temperature, defined as , where is a positive fluid property. The velocity of the stretching surface is defined as , with and as positive constants. The term is the heat generated (>0) or absorbed (<0) per unit volume and is defined as where is the temperaturedependent heat generation/absorption, that is, is positive in the case of the elastic sheet’s generation of heat and negative in the case of the sheet’s absorption of heat from the fluid flow. The radiative heat flux under Rosseland approximation (Brewster, 1992, cited in [23]) is given in the form where is the Stefan Boltzmann constant and is the mean absorption coefficient. We assume that the temperature difference within the flow is sufficiently small for to be expressed as a linear function of temperature. This is accomplished by expanding in a Taylor series about and neglecting higherorder terms, thus According to Wang [4], the surface temperature of the stretching sheet varies with the distance from the slot and time in the form where is the temperature at the slit, is the constant reference temperature for all . The surface of the planar liquid film is assumed to be smooth and free of surface waves while viscous shear stress and heat flux are assumed to vanish at the adiabatic free surface. Liu and Andersson [3] developed similarity transformations, which are given as where is the dimensionless film thickness and is the stream function defined by where a prime denotes differentiation with respect to . Apparently, (2.10) have already satisfied (2.1). Consequently, (2.1)–(2.4) are transformed to the following nonlinear boundary value problem: subject to where a prime denoting differentiation with respect to , is the dimensionless measure of unsteadiness, is the Prandtl number, is the radiation parameter defined as , and is an unknown constant which must be determined as a part of the present problem. The thermocapillarity number is given as It is worth noting that the momentum boundary layer problem defined by (2.11) subject to the relevant boundary conditions (2.13) is decoupled from the thermal boundary layer problem, while the temperature field is on the other hand coupled to the velocity field. The most important characteristics of flow and heat transfer are shear stress and heat flux of the stretching sheet, which are defined as where is the fluid dynamic viscosity. The local skinfriction coefficient and the local Nusselt number can be defined as Thus, the skin friction and the rate of heat transfer for fluid flow in a thin film can be expressed as where is the local Reynolds number.
3. Solution Approach
In this section we apply HAM to solve system (2.11)–(2.14). We assume that the solutions of and can be expressed by a set of base functions with the following forms: where and are constants. Under the rule of solution expression given by (2.11) and (2.12), subject to the boundary conditions (2.13) and (2.14), it is straightforward to choose as the initial guesses of and . The auxiliary linear operators and are chosen with the following properties: where , , , , and are constants of integration. From (2.11) and (2.12), the nonlinear operators are defined as where and are both unknown functions of and , while is a function dependent on . The prime denotes differentiation with respect to . Let and denote the nonzero auxiliary parameters, whereas and denote nonzero auxiliary functions. The zeroorder deformation equation can be constructed as subject to the boundary conditions where is an embedding parameter. From (3.2), it is straightforward to show that when , the solutions of (3.5) and (3.6) are Since and when , (3.5) and (3.6) are equivalent to (2.11)–(2.14), respectively, provided that Thus, as increases from 0 to 1, and vary from the initial guesses and to the solutions and in (2.11)–(2.14). also varies from the initial guess to the timescale parameter . Applying the Taylor series to (3.7) and (3.9), , , and can be expanded as series of , where Thus, using (3.8), we have By differentiating (3.5) times with respect to , then setting , and finally dividing by , the thorder deformation equations are obtained: subject to the boundary conditions for , where Let and denote the particular solutions of (3.13) as The general solutions for (3.3) are It should be pointed out that and contain the unknown parameter , which should be determined along with ,, , , and by the boundary conditions (3.14). We solve (3.13) for successfully, and the thorder approximations of , and are given, respectively, by
4. Convergence of the HAM Solution
We note the presence of the auxiliary parameters and in (3.13). Liao [16, 25] indicates that the convergence and rate of approximation of such series depend on the values of and . For the purpose of choosing the proper values for and , we set , , , , and . Figure 1(a) shows the variation of with using the tenthorder approximation. From the figure, it is seen that convergent result can be obtained when . Thus, we can choose an appropriate value for in this range to get convergent solution for . In Figure 1(b), we observe that convergent results of can be obtained by choosing a value of in the range of under the tenthorder approximation. From (2.13), we can see that depends on . Figure 1(c) shows the appropriate value for in the range of when under the tenthorder approximation. To see the accuracy of the solution, we define and as the residual errors of the system. Figure 2 shows the residual errors of the HAM solution in the case and .
(a)
(b)
(c)
(a)
(b)
5. Results and Discussion
We solved (3.13) and (3.14) with the aid of Maple, a symbolic computation software. The auxiliary functions and in (3.13) were set to be equal to 1 in all calculations reported in this paper.
In order to assess the accuracy of the numerical method, Table 1 compares our present results with those of previous investigations. We can claim that our numerical results are in excellent agreement with those of [4, 14, 24] under some limiting cases.

In Table 2, a good agreement was found in comparisons of values of free surface temperature between [24] and present work.

The effects of the radiation parameter on the film thickness , surface shear stress , free surface temperature , and heat flux at a specific thermocapillarity number are shown in Table 3. It is concluded that the film thickness and the free surface temperature increase as the radiation parameter increases and consequently the surface shear stress and dimensionless heat flux decrease.

Table 4 shows the effects of thermocapillarity number on the film thickness , surface shear stress , free surface temperature , and heat flux . We agree with [24] that the thermocapillarity parameter has the effect of enhancing the velocity and the local Nusselt number while decreasing the temperature and the surface shear stress.

The heat absorption sink leads to a decrease in the thermal boundary layer whereas the boundary layer thickness increases with increase in . These effects can be observed in Table 5. As the temperaturedependent parameter increases for both cases (heat source or heat sink), the free surface temperature increases while the surface shear stress and wall heat flux decrease.

The effects of the different values of radiation parameter on the temperature profiles are depicted in Figure 3. By fixing , , , and , the thin film flow swings from slight deceleration to higher velocity while the temperature increases as increases.
(a)
(b)
Thermocapillarity produced an outward flow along the free surface. Figure 4 shows that by increasing the value of the thermocapillarity number, the velocity decreases until it arrives at a minimum and then increases to its free surface value while the temperature consistently cools down.
(a)
(b)
Figure 5 shows that when the temperaturedependent parameter increases, the velocity decreases while the temperature increases. This is due to the fact that the internal heat source enhances or damps heat transport.
(a)
(b)
The unsteadiness parameter has the effect of increasing both velocity profiles and temperature distribution at a specific , , and , as shown in Figure 6. The effects of Prandtl number on temperature distribution is presented in Figure 7. As increases, the velocity profile increases and the temperature distribution decreases. The higher the Prandtl number, the cooler the temperature of the fluid flow.
(a)
(b)
(a)
(b)
6. Concluding Remarks
The effects of thermal radiation and thermocapillarity in a thin liquid film on an unsteady stretching sheet with nonuniform heat source/sink was analyzed successfully by means of the homotopy analysis method (HAM). With the presence of internal heat generation/absorption, the radiation parameter plays a significant role in controlling the temperature of the fluid flow by enhancing the temperature of fluid flow and decreasing the dimensionless heat flux. The thermocapillarity enhances the velocity and the dimensionless heat flux while reducing the temperature and the surface shear stress.
References
 C. Y. Wang, “Liquid film on an unsteady stretching surface,” Quarterly of Applied Mathematics, vol. 48, no. 4, pp. 601–610, 1990. View at: Google Scholar  Zentralblatt MATH
 H. I. Andersson, J. B. Aarseth, and B. S. Dandapat, “Heat transfer in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 43, no. 1, pp. 69–74, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 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 Site  Google Scholar
 C. Wang, “Analytic solutions for a liquid film on an unsteady stretching surface,” Heat and Mass Transfer, vol. 42, no. 8, pp. 759–766, 2006. View at: Publisher Site  Google Scholar
 H. I. Andersson, J. B. Aarseth, N. Braud, and B. S. Dandapat, “Flow of a powerlaw fluid film on an unsteady stretching surface,” Journal of NonNewtonian Fluid Mechanics, vol. 62, no. 1, pp. 1–8, 1996. View at: Publisher Site  Google Scholar
 C. H. Chen, “Heat transfer in a powerlaw fluid film over a unsteady stretching sheet,” Heat and Mass Transfer, vol. 39, no. 89, pp. 791–796, 2003. View at: Publisher Site  Google Scholar
 C. H. Chen, “Effect of viscous dissipation on heat transfer in a nonNewtonian liquid film over an unsteady stretching sheet,” Journal of NonNewtonian Fluid Mechanics, vol. 135, no. 23, pp. 128–135, 2006. View at: Publisher Site  Google Scholar
 C. Wang and I. Pop, “Analysis of the flow of a powerlaw fluid film on an unsteady stretching surface by means of homotopy analysis method,” Journal of NonNewtonian Fluid Mechanics, vol. 138, no. 23, pp. 161–172, 2006. View at: Publisher Site  Google Scholar
 M. A. A. Mahmoud and A. M. Megahed, “MHD flow and heat transfer in a nonNewtonian liquid film over an unsteady stretching sheet with variable fluid properties,” Canadian Journal of Physics, vol. 87, no. 10, pp. 1065–1071, 2009. View at: Publisher Site  Google Scholar
 B. S. Dandapat, B. Santra, and H. I. Andersson, “Thermocapillarity in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 46, no. 16, pp. 3009–3015, 2003. View at: Publisher Site  Google Scholar
 B. S. Dandapat, B. Santra, and K. Vajravelu, “The effects of variable fluid properties and thermocapillarity on the flow of a thin film on an unsteady stretching sheet,” International Journal of Heat and Mass Transfer, vol. 50, no. 56, pp. 991–996, 2007. View at: Publisher Site  Google Scholar
 C. H. Chen, “Marangoni effects on forced convection of powerlaw liquids in a thin film over a stretching surface,” Physics Letters Section A, vol. 370, no. 1, pp. 51–57, 2007. View at: Publisher Site  Google Scholar
 N. F. M. Noor and I. Hashim, “Thermocapillarity and magnetic field effects in a thin liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 53, no. 910, pp. 2044–2051, 2010. View at: Publisher Site  Google Scholar
 M. S. Abel, N. Mahesha, and J. Tawade, “Heat transfer in a liquid film over an unsteady stretching surface with viscous dissipation in presence of external magnetic field,” Applied Mathematical Modelling, vol. 33, no. 8, pp. 3430–3441, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 N. F. M. Noor, O. Abdulaziz, and I. Hashim, “MHD flow and heat transfer in a thin liquid film on an unsteady stretching sheet by the homotopy analysis method,” International Journal for Numerical Methods in Fluids, vol. 63, no. 3, pp. 357–373, 2010. View at: Publisher Site  Google Scholar
 S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, vol. 2 of CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2004.
 S. Abbasbandy, “Homotopy analysis method for the Kawahara equation,” Nonlinear Analysis. Real World Applications, vol. 11, no. 1, pp. 307–312, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 O. Abdulaziz, N. F. M. Noor, and I. Hashim, “Homotopy analysis method for fully developed MHD micropolar fluid flow between vertical porous plates,” International Journal for Numerical Methods in Engineering, vol. 78, no. 7, pp. 817–827, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 O. Abdulaziz and I. Hashim, “Fully developed free convection heat and mass transfer of a micropolar fluid between porous vertical plates,” Numerical Heat Transfer Part A, vol. 55, no. 3, pp. 270–288, 2009. View at: Publisher Site  Google Scholar
 T. Hayat, T. Javed, and Z. Abbas, “MHD flow of a micropolar fluid near a stagnationpoint towards a nonlinear stretching surface,” Nonlinear Analysis. Real World Applications, vol. 10, no. 3, pp. 1514–1526, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 T. Hayat, M. Mustafa, and S. Asghar, “Unsteady flow with heat and mass transfer of a third grade fluid over a stretching surface in the presence of chemical reaction,” Nonlinear Analysis. Real World Applications, vol. 11, no. 4, pp. 3186–3197, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 N. Kousar and S. Liao, “Unsteady nonsimilarity boundarylayer flows caused by an impulsively stretching flat sheet,” Nonlinear Analysis. Real World Applications, vol. 12, no. 1, pp. 333–342, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Abd ElAziz, “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 Site  Google Scholar
 M. A. A. Mahmoud, “Thermal radiation effects on the flow and heat transfer in a liquid film on an unsteady stretching sheet,” International Journal for Numerical Methods in Fluids, vol. 67, no. 11, pp. 1692–1701, 2011. View at: Publisher Site  Google Scholar
 S. Liao, “Notes on the homotopy analysis method: some definitions and theorems,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 983–997, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2012 R. C. Aziz et al. 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.