Abstract
The paper investigates the radiation effect on the magnetohydrodynamic Newtonian fluid flow over an exponentially stretching sheet. The effects of frictional heating and viscous dissipation on the heat transport are taken into account. The governing partial differential equations are transformed into ordinary differential equations using a suitable similarity transformation. Zeroorder analytical solutions of the momentum equation and confluent hypergeometric solutions of heat and mass transport equations are obtained. The accuracy of analytical solutions is verified by numerical solutions obtained using a shooting technique that uses a RungeKuttaFelhberg integration scheme and a NewtonRaphson correction scheme. The effects of the radiation parameter, the magnetic parameter, Gebhart and Schmidt numbers on the momentum, heat and mass transports are discussed. The skin friction and heat and mass transfer coefficients for various physical parameters are discussed.
1 Introduction
The study of laminar boundary layer flow over a stretching sheet has received considerable attention in the recent past due to its immense application in industry, for example, in extrusion processes such as the polymer extrusion from a dye and wire drawing. Other engineering applications of the stretching sheet problem include polymer sheet extrusion from a dye, drawing, tinning and annealing of copper wires, glass fiber and paper production, the cooling of a metallic plate in a cooling bath and so on. There has been tremendous amount of work on the stretching sheet problem in the past several decades (see Crane [1], Gupta and Gupta [2], Grubka and Bobba [3], Dutta and Gupta [4], Siddappa and Abel [5], Chen and Char [6], Laha et al.[7], Chakrabarti and Gupta [8], Anderson et al.[9], Siddheshwar and Mahabaleswar [10], Abel and Mahesha [11], Abel et al.[12] and the references therein).
The above studies concern the linear stretching sheet problem but most of the practical situations involve a nonlinear stretching sheet such as an exponential one. With this in mind, several authors have considered the velocity of the sheet to vary exponentially with the distance from the slit. Elbashbeshy [13] was among the first to study the exponentially stretching sheet problem. He considered a perforated sheet and examined the effect of wall mass suction on the flow and heat transfer over an exponentially stretching surface. Using a suitable similarity transformation, he transformed the momentum equation into a nonlinear Riccati type equation and solved it iteratively. Ishak [14] studied the MHD boundary layer flow due to an exponentially stretching sheet with radiation effect. He found that the local heat transfer rate at the surface decreased with increasing values of the magnetic and radiation parameters. The flow and heat transfer from an exponentially stretching surface was considered by Magyari and Keller [15]. They examined the heat and mass transfer characteristics and compared with the wellknown results of the powerlaw models. Sanjayanand and Khan [16] studied the heat and mass transfer in a viscoelastic boundary layer flow over an exponentially stretching sheet. They found that the viscoelastic parameter enhances the thermal boundary layer thickness. The effect of viscous dissipation on the mixed convection heat transfer from an exponentially stretching surface was studied by Partha et al.[17]. They observed a rapid growth in the nondimensional skin friction coefficient with the mixed convection parameter. The influence of thermal radiation on the boundary layer flow due to an exponentially stretching sheet is studied by Sajid and Hayat [18]. Khan [19] presented an elegant solution of the viscoelastic boundary layer flow over an exponentially stretching sheet in terms of Whittaker’s function.
The characteristics desired of the final product in an extrusion process depend on the rate of stretching and cooling. Hence, it is very important to have a controlled cooling environment where the flow over the stretching sheet can be regulated by external agencies like a magnetic field. An exponential variation of a magnetic field is used, among other applications, to determine the diamagnetic susceptibility of plasma. Steenbeck [20] determined the diamagnetic susceptibility of a cylindrical plasma for axial magnetic fields with various gas pressure and magnetic field strengths. Tonks [21] studied the effects of a magnetic field in the plasma of an arc. Pavlov [22] considered the magnetohydrodynamic flow of an incompressible viscous fluid over a linearly stretching surface. Sarpakaya [23] extended Pavlov’s work to nonNewtonian fluids. Subsequent studies by Andersson [24], Lawrence and Rao [25], Abel et al.[26], Cortell [27] concerned the magnetohydrodynamic flow of viscoelastic liquids over a stretching sheet. Radiation effects on MHD flow past an exponentially accelerated isothermal vertical plate with uniform mass diffusion in the presence of a heat source was studied by Reddy et al.[28]. They observed that the velocity decreases with an increase in the magnetic parameter due to a resistive drag force which tends to resist the fluid flow and thus reduces the velocity. The boundary layer thickness was also found to decrease with an increase in the magnetic parameter.
Most of the earlier work neglected radiation effects. If the polymer extrusion process is placed in a thermally controlled environment, radiation could become important. As with magnetohydrodynamics, careful control of thermal radiative heat transfer has an effect on the characteristics of the final product. Many researchers have considered the effect of thermal radiation on flows over stretching sheets. Studies by Raptis [29], Raptis and Perdikis [30] address the effect of radiation in various situations. Siddheshwar and Mahabaleswar [10] studied the effects of radiation and heat source on MHD flow of a viscoelastic liquid and heat transfer over a stretching sheet. Bidin and Nazar [31] studied the effects of numerical solution of the boundary layer flow over an exponentially stretching sheet with thermal radiation. They observed that the temperature profiles and the thermal boundary layer thickness increase slightly with an increase in the Eckert number. They also showed that an increase in Pr causes a decrease in temperature profiles and the thermal boundary layer thickness. Physically, if Pr increases, the thermal diffusivity decreases, and these phenomena lead to the decreasing of energy ability that reduces the thermal boundary layer. Elbashbeshy and Dimian [32] analyzed boundary layer flow in the presence of radiation effect and heat transfer over the wedge with a viscous coefficient. Thermal radiation effects on hydromagnetic flow due to an exponentially stretching sheet were studied by Reddy and Reddy [33]. They found that as radiation increases, the temperature profiles and thermal boundary layer thickness also increase. They also observed that the temperature profiles and thermal boundary layer thickness increase slightly with an increase in the Eckert number. Raptis et al.[34] studied the effect of thermal radiation on the magnetohydrodynamic flow of a viscous fluid past semiinfinite stationary plate and Hayat et al.[35] extended the analysis for the second grade fluid.
In addition to a magnetic field and thermal radiation, one has to consider the viscous dissipation effects due to frictional heating between fluid layers. The effect of viscous dissipation in natural convection processes has been studied by Gebhart [36] and Gebhart and Mollendorf [37]. They observed that the effect of viscous dissipation is predominant in vigorous natural convection and mixed convection processes. They also showed the existence of a similarity solution for the external flow over an infinite vertical surface with an exponential variation of surface temperature. Vajravelu and Hadjinicalaou [38] studied the heat transfer characteristics over a stretching surface with viscous dissipation in the presence of internal heat generation or absorption.
In this paper, we investigate the effects of various physical and fluid parameters such as the magnetic parameter, radiation parameter and viscous dissipation parameter on the flow and heat transfer characteristics of an exponentially stretching sheet. The momentum, energy and concentration equations are coupled and nonlinear. By using suitable similarity variables, these equations are converted into coupled ordinary differential equations and are solved analytically and numerically by using the RungeKuttaFehlberg and NewtonRaphson schemes.
2 Mathematical formulation
Consider the twodimensional magnetohydrodynamic flow of a Newtonian fluid over a stretching sheet. The origin of the system is located at the slit from which the sheet is drawn. The xaxis is taken along the continuous stretching surface and points in the direction of motion. The yaxis is perpendicular to the plate. The sheet velocity is assumed to vary as an exponential function of the distance x from the slit. The temperature and concentration far away from the fluid are assumed to be and respectively as shown in Figure 1. The sheetambient temperature and concentration differences are also assumed to be exponential functions of the distance x from the slit. A variable magnetic field of strength is applied normally to the sheet. Under the usual boundary layer approximation, subject to radiation and viscous dissipation effects, the equations governing the momentum, heat and mass transports can be written as
where u, v are the velocity components in the x, y directions respectively, ν is the kinematic viscosity, ρ is the density, σ is the electrical conductivity of the fluid, T is the temperature, C is the concentration, is the thermal diffusivity, k is the thermal conductivity, is the specific heat at constant pressure, is the radiative heat flux, and D is the species diffusivity.
Figure 1 . Schematics of the problem.
The boundary conditions for Equations (1)(4) have the form
Here the subscripts w, ∞ refer to the surface and ambient conditions respectively, , are positive constants, is the characteristic velocity, and L is the characteristic length.
To facilitate a similarity solution, the magnetic field is assumed to be of the form
where is a constant. It is also assumed that the fluid is weakly electrically conducting so that the induced magnetic field is negligible. Following Rosseland’s approximation, the radiative heat flux is modeled as
where is the StefanBoltzman constant, is the mean absorption coefficient. Assuming that the temperature differences within the flow are sufficiently small such that may be expressed as a linear function of temperature , we have
Continuity Equation (1) is satisfied by introducing a stream function ψ such that
The following similarity variables are used:
where η is the similarity variable, is the dimensionless stream function, is the dimensionless temperature, and is the dimensionless concentration.
On using Equations (6), (8) and (10), Equations (2)(5) transform into the following twopoint boundary value problem:
The nondimensional constants appearing in Equations (11)(13) are the magnetic parameter M, the radiation parameter K, the Prandtl number Pr, the Gebhart number Gb, and the Schmidt number Sc respectively defined as
3 Skin friction, heat and mass transfer coefficients
The parameters of engineering interest in heat and mass transport problems are the skin friction coefficient , the local Nusselt number , and the local Sherwood number . These parameters respectively characterize the surface drag, wall heat and mass transfer rates.
The shearing stress at the surface of the wall is given by
where μ is the coefficient of viscosity and is the Reynolds number. The skin friction coefficient is defined as
and using Equation (17) in Equation (18), we obtain
The heat transfer rate at the surface flux at the wall is given by
where k is thermal conductivity of the fluid. The Nusselt number is defined as
Using Equation (20) in Equation (21), the dimensionless wall heat transfer rate is obtained as follows:
The mass flux at the surface of the wall is given by
and the Sherwood is defined as
Using (23) in (24), the dimensionless wall mass transfer rate is obtained as
In Equations (19), (22) and (25), represents the local Reynolds number and it is defined as .
4 Analytical solution
4.1 Solution of momentum equation
The momentum boundary layer equation is partially decoupled from the energy and species equations. Integrating Equation (11) with η once over to the interval , we obtain
Integrating Equation (26) once again, we obtain
The solution procedure of Equation (28) can be reduced to the sequential solutions of the Riccati type equation of the form
This iteration algorithm has to be solved by substituting suitable zeroorder approximations for into the righthand side of Equation (28). We assume a zeroorder approximation as
which satisfies the condition at infinity. Integrating Equation (30) with respect to η and using the condition , we get
Using the above solution in Equation (27), the approximate value of s can be obtained as
Now substituting all the derivatives of zeroorder approximation into the righthand side of Equation (28), we obtain the equation for firstorder iteration as follows:
Further, we assume that the firstorder iterate of f satisfies the boundary conditions on f as given in (14). The above nonlinear Riccati type equation can be solved in terms of a confluent hypergeometric Whittaker function as discussed by Khan [19]. However, we restrict ourselves to the zeroorder solution, and similarly, to heat and mass transport equations.
4.2 Solution of heat transfer equation
Using the zeroorder approximations of f and and further introducing a new variable
Equation (12) and the thermal boundary conditions (15) take the form
where is the modified Prandtl number. The solution of Equation (35) is assumed in the form of
where is the complementary solution and is the particular solution. The complementary solution of Equation (35) is obtained in terms of confluent hypergeometric function in the following form:
where
is Kummer’s function (see Abramowitz and Stegun [39]) and
The particular solution is obtained as
where
Now, the complete solution can be written as
Making use of the boundary conditions (36) and rewriting the solution in terms of the variable η, we get
where
4.3 Solution of mass transfer equation
Using the zeroorder approximation of f and and further introducing a new variable
Equation (13) and the thermal boundary conditions in (16) take the form
where is the modified Schmidt number. Following the solution procedure discussed in the case of the energy equation, the solution of Equation (42) is obtained in terms of confluent hypergeometric function as
5 Solution procedure
The set of nonlinear differential Equations (11)(13) subject to the boundary conditions (14)(16) were solved numerically using an efficient RungeKuttaFehlberg method with a shooting technique, which is described in Pal and Shivakumara [40]. The most important step in this method is to choose an appropriate finite value of . In order to determine for the boundary value problem described by Equations (11)(13), we start with initial guess values for a particular set of physical parameters to obtain , and . The solution procedure is repeated with another large value of until two successive values of , and differ only by a specified significant digit. The value of η may change for a different set of physical parameters. Once the appropriate value of η is determined, the coupled boundary value problem given by Equations (11)(13) is solved numerically using the method of superposition. In this method, thirdorder nonlinear Equation (11), secondorder Equations (12) and (13) have been reduced to five ordinary differential equations as follows:
where
and a prime denotes differentiation with respect to η. The boundary conditions now become
where , and are determined such that , and . Thus, to solve this system, we require six initial conditions. However, since we have only three initial conditions for f and two initial conditions for θ and ϕ, the conditions , , are to be determined by the shooting method using the initial guess values , and until the conditions , and are satisfied. In this paper, we employed the shooting technique with the RungeKuttaFehlberg scheme to determine two more unknowns in order to convert the boundary value problem to an initial value problem. Once all the six initial conditions were determined, the resulting differential equations were integrated using an initial value solver. For this purpose, the fifthorder RungeKuttaFehlberg integration scheme was used.
6 Results and discussion
Analytical and numerical solutions were obtained for the effects of radiation and viscous dissipation for the MHD flow over an exponentially stretching sheet. Similarity transformations were used to transform the governing partial differential equations of flow, heat and mass transfer into a system of nonlinear ordinary differential equations. The zeroorder approximate solution for the dimensionless stream function f has been obtained analytically. Solutions of the energy and species equations were obtained in terms of confluent hypergeometric functions. The accuracy of the method was established by comparing the analytical solution with the numerical solution obtained by a shooting method together with RungeKuttaFehlberg and NewtonRaphson schemes. The skin friction, heat and mass transfer coefficients are tabulated in Tables 13. The effects of magnetic, radiation and viscous dissipation parameters on the velocity , temperature and concentration profiles are shown in Figures 26.
Figure 2 . Effect of the magnetic parameter (M) on velocity profile for,,,.
Table 1 . A comparison ofobtained by the analytical method with the shooting technique for different values ofM
Table 1 provides values of the skin friction coefficient for different values of the magnetic parameter M. Increasing values of M result in considerable opposition to the flow in the form of a Lorenz drag which enhances the values of the skin friction coefficient. Table 2 highlights the effect of the magnetic field, radiation and dissipation on the dimensionless wall temperature gradient. It is evident that all the three parameters reduce the values of the wall temperature gradient. Table 3 shows that the increase in Schmidt numbers leads to the increase in the dimensionless wall concentration gradient, while the opposite trend is observed in the case of the magnetic parameter. The results confirm a good agreement between analytical and numerical results.
Table 2 . A comparison ofobtained by the analytical method with the shooting technique for different values ofM,GbandKfor fixed values of
Table 3 . A comparison ofobtained by the analytical method with the shooting technique for different values ofM,Sc
The skin friction coefficients are shown in Table 4 for different values of the magnetic parameter in the absence of the physical parameters (i.e., ). We observe that skin friction coefficient increases with an increase in the magnetic parameter. It is interesting to note that the value of the wall skinfriction coefficient in the nonmagnetic () and magnetic () cases are in good agreement with the results presented by Reddy and Reddy [33].
Table 4 . A comparison offor different values ofMfor fixed values of
Figure 2 shows the variation of the velocity profile against the magnetic parameter. We notice that the effect of the magnetic parameter is to reduce the velocity of the fluid in the boundary layer region. This is due to an increase in the Lorentz force, similar to Darcy’s drag observed in the case of flow through a porous medium. This adverse force is responsible for slowing down the motion of the fluid in the boundary layer region. These results are similar to the results obtained by Reddy and Reddy [33].
The variation of the temperature distribution with the magnetic parameter is shown in Figure 3. The thermal boundary layer thickness increases with increasing values of the magnetic parameter. The opposing force introduced in the form of the Lorentz drag contributes in increasing the frictional heating between the fluid layers, and hence energy is released in the form of heat. This results in thickening of the thermal boundary layer.
Figure 3 . Effect ofMon temperature profile for,,,.
The effect of the magnetic parameter on the concentration profile is shown in Figure 4. It is observed that increases in the values in M result in thickening of the species boundary layer.
Figure 4 . Effect ofMon concentration profile for,,,.
The influence of the thermal radiation parameter K on temperature is shown in Figure 5. It is clear that thermal radiation enhances the temperature in the boundary layer region. Thus radiation should be kept at its minimum in order to facilitate better cooling environment. The radiation parameter K defines the relative contribution of conduction heat transfer to thermal radiation transfer.
Figure 5 . Effect of the radiation parameter on temperature profile for,,,.
The effect of the Gebhart number Gb on the heat transfer is shown in Figure 6. It is clear that the temperature in the boundary layer region increases with an increase in the viscous dissipation parameter. We also note that since the energy equation is partially decoupled from the momentum and species conservation equations, the parameters affecting the energy equation, namely, the Prandtl number, the radiation parameter and the Gebhart number, do not alter velocity and concentration profiles. We also observe a good agreement between the analytical and numerical solutions through Figures 26.
Figure 6 . Effect of viscous dissipation on temperature profile for,,,.
7 Conclusions
The problem of hydromagnetic Newtonian liquid flow due to an exponentially stretching sheet in the presence of radiation and viscous dissipation effects has been analyzed. Exact solutions were found in terms of hypergeometric functions, and a comparison of analytical and numerical results was shown. We found that the effect of the magnetic parameter is to reduce the velocity of the fluid in the boundary layer region. It was also observed that the increase in values of M results in thickening of the species boundary layer. The combined and individual effects of the magnetic parameter M, the radiation parameter K, and the viscous dissipation parameter Gb are to increase the heat transfer rates. Under some limiting conditions when the parameters Pr, Sc, K, Gb are zero, the current results agree well with available results in the literature.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
The problem was conceived in discussions by the authors. PKK and MN carried out the analytical and numerical computations, while GM and PS participated in the design of the study and drafted the manuscript.
Acknowledgements
The authors are grateful to the National Research Foundation (NRF) and the University of KwaZuluNatal for financial support.
References

Crane, LJ: Flow past a stretching plate . Z. Angew. Math. Phys.. 21, 645–647 (1970). Publisher Full Text

Gupta, PS, Gupta, AS: Heat and mass transfer on a stretching sheet with suction or blowing . Can. J. Chem. Eng.. 55, 744–746 (1977). Publisher Full Text

Grubka, LG, Bobba, KM: Heat transfer characteristics of a continuous stretching surface with variable temperature . J. Heat Transf.. 107, 248–250 (1985). Publisher Full Text

Dutta, BK, Gupta, AS: Cooling of a stretching sheet in a various flow . Ind. Eng. Chem. Res.. 26, 333–336 (1987). Publisher Full Text

Siddappa, B, Abel, MS: NonNewtonian flow past a stretching plate . Z. Angew. Math. Phys.. 36, 47–54 (1985). Publisher Full Text

Chen, CK, Char, MI: Heat transfer of a continuous stretching surface with suction or blowing . J. Math. Anal. Appl.. 135, 568–580 (1988). Publisher Full Text

Laha, MK, Gupta, PS, Gupta, AS: Heat transfer characteristics of the flow of an incompressible viscous fluid over a stretching sheet . WärmeStoffübertrag.. 24, 151–153 (1989). PubMed Abstract  PubMed Central Full Text

Chakarabarti, A, Gupta, AS: Hydromagentic flow and heat transfer over a stretching sheet . Q. Appl. Math.. 37, 73–78 (1979)

Andersson, HI, Bech, KH, Dandapat, BS: Magnetohydrodynamic flow of a powerlaw fluid over a stretching sheet . Int. J. NonLinear Mech.. 27, 929–936 (1992). Publisher Full Text

Siddheshwar, PG, Mahabaleswar, US: Effects of radiation and heat source on MHD flow of a viscoelastic liquid and heat transfer over a stretching sheet . Int. J. NonLinear Mech.. 40, 807–820 (2005). Publisher Full Text

Abel, MS, Mahesha, N: Heat transfer in MHD viscoelastic fluid flow over a stretching sheet with variable thermal conductivity, nonuniform heat source and radiation . Appl. Math. Model.. 32, 1965–1983 (2008). Publisher Full Text

Abel, MS, Siddheshwar, PG, Mahesha, N: Effects of thermal buoyancy and variable thermal conductivity on the MHD flow and heat transfer in a powerlaw fluid past a vertical stretching sheet in the presence of a nonuniform heat source . Int. J. NonLinear Mech.. 44, 1–12 (2009). Publisher Full Text

Elbashbeshy, EMA: Heat transfer over an exponentially stretching continuous surface with suction . Arch. Mech.. 53(6), 643–651 (2001)

Ishak, A: MHD boundary layer flow due to an exponentially stretching sheet with radiation effect . Sains Malays.. 40(4), 391–395 (2011)

Magyari, E, Keller, B: Heat and mass transfer in the boundary layers on an exponentially stretching continuous surface . J. Phys. D, Appl. Phys.. 32, 577–585 (1999). Publisher Full Text

Sanjayanand, E, Khan, SK: On heat and mass transfer in a viscoelastic boundary layer flow over an exponentially stretching sheet . Int. J. Therm. Sci.. 45, 819–828 (2006). Publisher Full Text

Partha, MK, Murthy, PVSN, Rajasekhar, GP: Effect of viscous dissipation on the mixed convection heat transfer from an exponentially stretching surface . Heat Mass Transf.. 41, 360–366 (2005). Publisher Full Text

Sajid, M, Hayat, T: Influence of thermal radiation on the boundary layer flow due to an exponentially stretching sheet . Int. Commun. Heat Mass Transf.. 35, 347–356 (2008). Publisher Full Text

Khan, SK: Boundary layer viscoelastic fluid flow over an exponentially stretching sheet . Int. J. Appl. Mech. Eng.. 11(2), 321–335 (2006)

Steenbeck, M: Wissen. Veroff . Siemens. 15(2), Article ID 1 (1936)

Tonks, L: Theory of magnetic effects in the plasma of an arc . Phys. Rev.. 56, 360–373 (1939). Publisher Full Text

Pavlov, KB: Magnetohydromagnetic flow of an incompressible viscous fluid caused by deformation of a surface . Magn. Gidrodin.. 4, 146–147 (1974)

Sarpakaya, T: Flow of nonNewtonian fluids in a magnetic field . AIChE J.. 7, 324–328 (1961). Publisher Full Text

Andersson, HI: MHD flow of a viscoelastic fluid past a stretching surface . Acta Mech.. 95, 227–230 (1992). Publisher Full Text

Lawrence, PS, Rao, BN: The nonuniqueness of the MHD flow of a viscoelastic fluid past a stretching sheet . Acta Mech.. 112, 223–228 (1995). Publisher Full Text

Abel, MS, Joshi, A, Sonth, RM: Heat transfer in MHD viscoelastic fluid flow over a stretching surface . Z. Angew. Math. Mech.. 81, 691–698 (2001)

Cortell, R: Flow and heat transfer of an electrically conducing fluid of second grade over a stretching sheet subject to suction and to a transverse magnetic field . Int. J. Heat Mass Transf.. 49, 1851–1856 (2006). Publisher Full Text

Reddy, PBA, Reddy, NB, Suneetha, S: Radiation effects on MHD flow past an exponentially accelerated isothermal vertical plate with uniform mass diffusion in the presence of heat source . J. Appl. Fluid Mech.. 5(3), 119–126 (2012)

Raptis, A: Radiation and viscoelastic flow . Int. Commun. Heat Mass Transf.. 26(6), 889–895 (1999). Publisher Full Text

Raptis, A, Perdikis, C: Viscoelastic flow by the presence of radiation . Z. Angew. Math. Mech.. 78, 277–279 (1998). Publisher Full Text

Bidin, B, Nazar, R: Numerical solution of the boundary layer flow over an exponentially stretching sheet with thermal radiation . Eur. J. Sci. Res.. 33(4), 710–717 (2009)

Elbashbeshy, EMA, Dimian, MF: Effect of radiation on the flow and heat transfer over a wedge with variable viscosity . Appl. Math. Comput.. 132, 445–454 (2002). Publisher Full Text

Reddy, PBA, Reddy, NB: Thermal radiation effects on hydromagnetic flow due to an exponentially stretching sheet . Int. J. Appl. Math. Comput.. 3(4), 300–306 (2011)

Raptis, A, Perdikis, C, Takhar, HS: Effect of thermal radiation on MHD flow . Int. J. Heat Mass Transf.. 153, 645–649 (2004)

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 Transf.. 50, 931–941 (2007). Publisher Full Text

Gebhart, B: Effect of viscous dissipation in natural convection . J. Fluid Mech.. 14, 225–232 (1962). Publisher Full Text

Gebhart, B, Mollendorf, J: Viscous dissipation in external natural convection flows . J. Fluid Mech.. 38, 97–107 (1969). Publisher Full Text

Vajravelu, K, Hadjinicolaou, A: Heat transfer in a viscous fluid over a stretching sheet with viscous dissipation and internal heat generation . Int. Commun. Heat Mass Transf.. 20, 417–430 (1993). Publisher Full Text

Abramowitz, M, Stegun, LA: Handbook of Mathematical Functions (1972)

Pal, D, Shivakumara, IS: Mixed convection heat transfer from a vertical plate embedded in a sparsely packed porous medium . Int. J. Appl. Mech. Eng.. 11(4), 929–939 (2006)