# Thermal Boundary Layer in Flow due to an Exponentially Stretching Surface with an Exponentially Moving Free Stream

Authors: Krishnendu Bhattacharyya; G. C. Layek
Journal: Modelling and Simulation in Engineering
 (2014)

---

## Abstract

A numerical investigation is made to study the thermal boundary layer for flow of incompressible Newtonian fluid over an exponentially stretching sheet with an exponentially moving free stream. The governing partial differential equations are transformed into self-similar ordinary differential equations using similarity transformations in exponential forms. Then those are solved numerically by shooting technique using Runge-Kutta method. The study reveals that the momentum boundary layer thickness for this flow is considerably smaller than the linear stagnation point flow past a linearly stretching sheet. The momentum and thermal boundary layer thicknesses reduce when the velocity ratio parameter increases. For the temperature distribution, in addition to the heat transfer from the sheet, the heat absorption at the sheet also occurs in certain situations and both heat transfer and absorption increase with the velocity ratio parameter and the Prandtl number. The temperature inside the boundary layer significantly decreases with higher values of velocity ratio parameter and the Prandtl number.


---

## Body

## 1. Introduction

The viscous fluid flow due to a stretching sheet is very significant problem in fluid dynamics due to its huge applications in many manufacturing processes, 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, the boundary layer along a liquid film in the condensation processes, hot rolling, paper production, metal spinning, glass-fiber production, and drawing of plastic films. The heat transfer from a stretching surface is of interest in polymer extrusion processes where the object, after passing through a die, enters the fluid for cooling below a certain temperature. During the processes, mechanical properties are greatly dependent upon the rate of cooling and the rate at which such objects are cooled has an important bearing on the properties of the final product. So, the quality of the final product depends on the rate of heat transfer from the stretching surface. The suitable choice of cooling liquid is very much crucial as it has a direct impact on rate of cooling of the surface.Crane [1] first investigated the steady boundary layer flow of an incompressible viscous fluid over a linearly stretching plate and gave an exact similarity solution in closed analytical form. Crane’s work was extended by many researchers such as P. S. Gupta and A. S. Gupta [2], Chen and Char [3], Pavlov [4], and Ali [5] considering the effects of heat and mass transfer and magnetic field under various physical conditions. Later, Ali [6] discussed the thermal boundary layer on a surface stretching with power law velocity. The mixed convection for the flow due to linearly stretching sheet with mass suction/injection was studied by Ali and Al-Yousef [7]. Ali [8] also investigated the buoyancy effects on the boundary layers induced by a rapidly stretching surface. Tsai et al. [9] discussed the effect of nonuniform heat source/sink on the flow and heat transfer from an unsteady stretching sheet through a quiescent fluid medium extending to infinity. In two important papers, Bhattacharyya and Layek [10, 11] explained the behaviour of chemically reactive solute distribution in MHD flow over a stretching sheet and in slip flow towards a vertical stretching sheet. Recently, Bhattacharyya [12, 13] presented effect of heat source/sink on the MHD boundary layer flow over steady and unsteady shrinking/stretching sheet and Bhattacharyya et al. [14] found the analytic solution of MHD flow of non-Newtonian Casson fluid flow induced due to stretching/shrinking sheet with wall mass transfer.On the other hand, Hiemenz [15] first studied the steady flow in the neighbourhood of a stagnation point. Chiam [16] investigated a problem which is a combination of the works of Hiemenz [15] and Crane [1], that is, the linear stagnation point flow towards a linear stretching sheet when the stretching rate of the plate is equal to the strain rate of the stagnation point flow and he found no boundary layer structure near the plate. After few years, Mahapatra and Gupta [17] reinvestigated the same stagnation point flow towards a stretching sheet with different stretching and straining velocities and they found two kinds of boundary layers near the sheet depending on the ratio of the stretching and straining velocity rates. Wang [18] demonstrated the stagnation point flow towards a shrinking sheet. Furthermore, the behaviour of stagnation point flow over stretching/shrinking sheet under different physical aspects was discussed by many researchers [19–35].In last few decades in almost all investigations of the flow over a stretching sheet, the flow occurs because of linear variation of stretching velocity of the flat sheet with the distance from the origin. So, the boundary layer flow induced by an exponentially stretching sheet is not studied much though it is very important and practical flow frequently appeared in many engineering processes. In 1999, Magyari and Keller [36] first consider the boundary layer flow due to an exponentially stretching sheet and he also studied the heat transfer in the flow taking exponentially varied wall temperature. After that, Elbashbeshy [37] numerically examined the flow and heat transfer over an exponentially stretching surface considering wall mass suction. Khan and Sanjayanand [38] investigated the flow of viscoelastic fluid and heat transfer over an exponentially stretching sheet with viscous dissipation effects. Partha et al. [39] presented a similarity solution for mixed convection flow past an exponentially stretching surface by taking into account the influence of viscous dissipation on the convective transport. Sanjayanand and Khan [40] discussed the effects of heat and mass transfer on the boundary layer flow of viscoelastic fluid taking nth order chemical reaction. Al-Odat et al. [41] investigated the effect magnetic field on thermal boundary layer on an exponentially stretching continuous surface with an exponential temperature distribution. Later, Sajid and Hayat [42] showed the influence of thermal radiation on the boundary layer flow past an exponentially stretching sheet and they reported series solutions for velocity and temperature using homotopy analysis method (HAM) and Bidin and Nazar [43] obtained a numerical solution of the same problem. Ishak [44] studied the MHD boundary layer flow over an exponentially stretching sheet in presence of thermal radiation. Mandal and Mukhopadhyay [45] analyzed the heat transfer in the flow due to an exponentially stretching porous sheet with surface heat flux embedded in a porous medium. Bhattacharyya [46] discussed the boundary layer flow and heat transfer past an exponentially shrinking sheet and Bhattacharyya and Pop [47] presented the effect of magnetic field on that flow. Recently, The MHD boundary layer nanofluid flow due to an exponentially stretching permeable sheet was demonstrated by Bhattacharyya and Layek [48].Inspired by the above investigations, in the present paper, the thermal boundary layer in incompressible flow over an exponentially stretching sheet in an exponentially moving free stream is studied. The wall temperature distribution is taken variable in exponential form. Using similarity transformations, the governing partial differential equations are transformed into self-similar ordinary differential equations. Then the transformed equations are solved by standard shooting technique using fourth order Runge-Kutta method. The numerical computations are presented through some figures and the various characteristics of flow and heat transfer are thoroughly described. This type of flow is possible in reality in some special situations. The motive of the investigation is to study new type of flow dynamics and the heat transfer in the flow field.


## 2. Analysis of Motion and Heat Transfer

Consider the steady two-dimensional laminar boundary layer flow of a viscous incompressible fluid and heat transfer over an exponentially stretching sheet (of velocityUw) with an exponentially moving free stream (of velocity U∞). The sheet coincides with the plane y=0 and the flow confined to y>0. The governing equations for the flow and the temperature are written in usual notation as
(1)∂u∂x+∂v∂y=0,(2)u∂u∂x+v∂u∂y=U∞dU∞dx+υ∂2u∂y2,(3)u∂T∂x+v∂T∂y=κρcp∂2T∂y2,
where u and v are the velocity components in x and y directions, respectively, U∞ is the fluid velocity in free stream, υ(=μ/ρ) is the kinematic fluid viscosity, ρ is the fluid density, μ is the coefficient of fluid viscosity, T is the temperature, κ is the fluid thermal conductivity, and cp is the specific heat. The diagram of the physical problem is given in Figure 1.Figure 1
A diagram of the physical problem.The boundary conditions corresponding to the velocity components and the temperature are given by(4)u=Uw(x),v=0aty=0;u⟶U∞(x)asy⟶∞,(5)T=Tw=T∞+T0exp⁡(λx2L)aty=0;T⟶T∞asy⟶∞,
where Tw is the variable temperature of the sheet, T∞ is the free stream temperature assumed to be constant, T0 is a constant which measures the rate of temperature increase along the sheet, L denotes the reference length, and λ is a parameter which is physically very important in controlling the exponential increment of temperature along the sheet and it may have both positive and negative values. The stretching and free stream velocities Uw and U∞ are, respectively, given by
(6)Uw(x)=cexp⁡(xL),U∞(x)=aexp⁡(xL),
where c and a are constants with c>0 and a≥0 (for without any free stream a=0).Now, the stream functionψ(x,y) is given by
(7)u=∂ψ∂y,v=-∂ψ∂x.
For relations in (7), the continuity equation (1) is identically satisfied and the momentum equation (2) and the temperature equation (3) are reduced to the following forms:
(8)∂ψ∂y∂2ψ∂x∂y-∂ψ∂x∂2ψ∂y2=U∞dU∞dx+υ∂3ψ∂y3,∂ψ∂y∂T∂x-∂ψ∂x∂T∂y=κρcp∂2T∂y2.
The boundary conditions in (4) for the flow become
(9)∂ψ∂y=Uw(x),∂ψ∂x=0aty=0;∂ψ∂y⟶U∞(x)asy⟶∞.
Next, the dimensionless variables for ψ and T are introduced as
(10)ψ=2υLcf(η)exp⁡(x2L),T=T∞+(Tw-T∞)θ(η),
where η is the similarity variable and is defined as η=yc/2υLexp⁡(x/2L).Using (10), finally the following nonlinear self-similar equations are obtained:
(11)f′′′+ff′′-2f′2+2(ac)2=0,(12)θ′′+Pr(fθ′-λf′θ)=0,
where a/c is the velocity ratio parameter and Pr=μcp/κ is the Prandtl number.The boundary conditions in (9) and (5) reduce to the following forms:
(13)f(η)=0,f′(η)=1atη=0;f′(η)⟶acasη⟶∞,θ(η)=1atη=0;θ(η)⟶0asη⟶∞.
When c=a, (11) gives closed form analytical solution f(η)=η which is the same as the linear stagnation-point flow over a stretching sheet with velocity varied linearly to the distance from the origin.The physical quantities of interest in this problem are the local skin-friction coefficient and the local Nusselt number. The dimensionless local skin-friction coefficient is expressed as(14)Cfx=τw|y=0ρUw2,
where τw is the shear stress at the wall and is given by τw=μ[∂u/∂y]y=0. So,
(15)2ReCfx=f′′(0),
where Re=UwL/υ is the Reynolds number.The local Nusselt number can be written as(16)Nux=-xTw-T∞∂T∂y|y=0.So,(17)NuxRex=-x2Lθ′(0),
where Rex=Uwx/υ is the local Reynolds number.


## 3. Numerical Method for Solution

The nonlinear coupled differential equations (11) and (12) along with the boundary conditions (13) form a two point boundary value problem (BVP) and is solved using shooting method, by converting it into an initial value problem (IVP). In this method, it is necessary to choose a suitable finite value of η→∞, say η∞. The following first-order system is set:
(18)f′=p,p′=q,q′=2p2-fq-2(ac)2,θ′=z,z′=-Pr(fz-λpθ)
with the boundary conditions
(19)f(0)=0,p(0)=1,θ(0)=1.
To solve (18) with (19) as an IVP, the values for (0), that is, f′′(0), and z(0), that is, θ′(0), are needed but no such values are given. The initial guess values for f′′(0) and θ′(0) are chosen and the fourth-order Runge-Kutta method is applied to obtain the solution. The calculated values of f′(η) and θ(η) at η∞(=25) are compared with the given boundary conditions f′(η∞)=a/c and θ(η∞)=0 and adjust values of f′′(0) and θ′(0) using “secant method” to give better approximation for the solution. The step size is taken as Δη=0.01. The process is repeated until we get the results correct up to the desired accuracy of 10-6 level.


## 4. Results and Discussion

The numerical computations have been carried out for various values of the physical parameters involved in the equations, namely, the velocity ratio parametera/c, the Prandtl number Pr, and the parameter λ. For illustration of the results, computed values are plotted through graphs and the physical explanations are rendered for all cases.To ensure the accuracy of the numerical scheme, we compare the values off′′(0) and f(∞) with the published data by Magyari and Keller [36] in Table 1 without any exponential free stream, U∞=0 that is a/c=0(a=0) and those are found in excellent agreement. Thus, we are feeling confident that the presented results are accurate.Table 1
The values off′′(0) and f(∞) for a/c=0.
Magyari and Keller [36]
Present study
−f′′(0)
1.281808
1.2818084
f
(
∞
)
0.905639
0.9056433As similar to that of stagnation point flow with linear straining velocity towards a linearly stretching sheet, in the boundary layer flow with exponentially varied free stream velocity over exponentially stretching sheet two different kinds of boundary layer structures have formed near the sheet depending upon the ratio of two constants relating to the stretching and free stream velocities, that is, on the velocity ratio parametera/c, for a/c>1 and a/c<1. Also, it is important to note that for a/c=1 no velocity boundary layer is formed near the sheet. The velocity profiles for various values of a/c are depicted in Figure 2. Also, the variations of f(η) and the dimensionless temperature for different values of a/c are plotted in Figures 3 and 4. It is observed that f(η) increases with increasing values of a/c and the temperature at a point decreases with a/c.Figure 2
Velocity profilesf′(η) for various values of a/c.Figure 3
The variations off(η) for various values of a/c.Figure 4
Temperature profilesθ(η) for various values of a/c.The momentum and thermal boundary layer thicknesses are denoted byδ and δT, respectively, and are described by the equations δ=ηδ2υL/cexp⁡(-x/2L) and δT=ηδT2υL/cexp⁡(-x/2L). The dimensionless boundary layer thicknesses ηδ and ηδT are defined as the values of η (nondimensional distance from the surface) at which the difference of dimensionless velocity f′(η) and the parameter a/c has been reduced to 0.001 and the dimensionless temperature θ(η) has been decayed to 0.001, respectively. Both boundary layer thicknesses are very vital in practical and theoretical standpoints. The values ηδ and ηδT are given in Table 2 for several values of a/c. From the table, it is seen that the momentum boundary layer thickness decreases when a/c increases (both for a/c>1 and a/c<1). But, it is worth noting that the momentum boundary layer thickness is considerably thinner than that of the boundary layer of linear stagnation point flow over a linearly stretching sheet which is obtained by Mahapatra and Gupta [17]. In exponential flow past an exponentially stretching sheet though the flow dynamics is similar to that of linear stagnation point flow, the momentum boundary layer thickness is smaller and it is physically important. Due to the fact that larger velocities of the free stream and the surface stretching make the nonzero viscosity layer thinner, it makes boundary layer thickness thinner. Similar to the velocity field, the thermal boundary layer thickness also decreases with increase of a/c.Table 2
Values ofηδ and ηδT for several values of a/c with Pr=0.5 and λ=1.
a
/
c
0.1
0.2
0.5
1.5
2.0
2.5
3.0
Mahapatra and Gupta [17]
η
δ
6.96
5.91
4.36
—
2.62
—
2.30
Present study
η
δ
5.04
4.17
2.91
1.87
1.79
1.68
1.59
η
δ
T
8.97
7.78
5.63
3.54
3.11
2.80
2.57The effect of the parameterλ on the temperature distribution is very significant because the increase of wall temperature along the sheet depends upon λ. For various values of λ, the dimensionless temperature profiles θ(η) are presented in Figures 5 and 6 for a/c=2 and a/c=0.1, respectively. In both the cases, the figures exhibit the overshoot of temperature profiles for some negative values of λ. Hence, when the wall temperature in inverse exponential forms, the thermal overshoot occurs. Also, it is noticed that the temperatures at particular point and the thermal boundary layer thickness decrease with increase of λ. Furthermore, it is very important to note that the overshoot pick is of higher length when free stream velocity is large compared to the stretching velocity (a/c=2). The dimensionless temperature gradient profiles θ′(η) for different λ are depicted in Figures 7 and 8 for a/c=2 and a/c=0.1, respectively. From these two figures, it is clearly seen that the temperature gradient at the sheet is positive for some values of λ and negative for some other values of λ. So, heat transfer from the sheet and heat absorption at the sheet occur depending on the values of λ for both values of a/c. Because of inverse exponential distribution of surface temperature (λ>0), the heat absorption is found in the thermal boundary layer. Thus the parameter λ related to wall temperature plays a vital role in controlling the cooling process of the surface and consequently it has a great impact on the quality of the manufactured products.Figure 5
Temperature profilesθ(η) for various values of λ with a/c=2.Figure 6
Temperature profilesθ(η) for various values of λ with a/c=0.1.Figure 7
Temperature gradient profilesθ′(η) for various values of λ with a/c=2.Figure 8
Temperature gradient profilesθ′(η) for various values of λ with a/c=0.1.The temperature profiles for various values of Prandtl number Pr are demonstrated in Figures9–12 with a/c=2, 0.1 and λ=1, −1.5. In Figures 9 and 10, for λ=1, it is observed that the temperature at a fixed point decreases with Pr for both values of a/c. Accordingly, the thermal boundary layer thickness decreases with increasing values of Pr and the heat transfers from the sheet to the adjacent fluid. But, when λ=-1.5 (Figures 11 and 12), the temperature overshoot is observed; that is, the sheet absorbs the heat and interestingly for both values of a/c the height of the overshoot pick increases with Pr, so the heat absorption at the sheet also increases. However, similar to the linear case, in this situation the thermal boundary layer thickness decreases with increasing values of Pr. For increase of Prandtl number, the fluid thermal conductivity reduces and consequently the thermal boundary layer thickness becomes thinner.Figure 9
Temperature profilesθ(η) for various values of Pr with a/c=2 and λ=1.Figure 10
Temperature profilesθ(η) for various values of Pr with a/c=0.1 and λ=1.Figure 11
Temperature profilesθ(η) for various values of Pr with a/c=2 and λ=-1.5.Figure 12
Temperature profilesθ(η) for various values of Pr with a/c=0.1 and λ=-1.5.The values off′′(0) related to local skin-friction coefficient and -θ′(0) related to local Nusselt number are plotted in Figures 13–15 for various values of a/c, λ and Pr. The value of f′′(0) increases with velocity ratio parameter a/c. Also, the value of -θ′(0) increases with λ. Moreover, for the larger negative values of λ, the heat absorption [-θ′(0)<0] at the surface is found. With the increase of a/c and Pr, the heat transfer enhances and also heat absorption becomes larger (Figures 14 and 15).Figure 13
Values off′′(0) and -θ′(0) for various values of a/c.Figure 14
Values of-θ′(0) versus λ for various values of a/c with Pr=0.5.Figure 15
Values of-θ′(0) versus λ for various values of a/c and Pr.


## 5. Concluding Remarks

The momentum and heat transfer characteristics of the laminar boundary layer flow induced by an exponentially stretching sheet in an exponentially moving free stream are investigated. The transformed governing equations are solved by shooting technique using Runge-Kutta method. The findings of this study can be summarized as follows.(a)
The momentum boundary layer thickness for this flow is significantly smaller than that of the linear stagnation point flow past a linearly stretching sheet.(b)
The momentum and thermal boundary layer thicknesses decrease with increase of velocity ratio parametera/c.(c)
For positive values ofλ and smaller negative values heat transfers from surface to the ambient fluid and importantly for greater negative value of λ heat transfers from the ambient fluid to the surface, the heat absorption occurs.(d)
Temperature and thermal boundary layer thicknesses decrease withλ, a parameter associated with wall temperature.(e)
The heat transfer from the sheet and the heat absorption at the sheet are enhanced with the increase in both velocity ratio parameter and Prandtl number.



---

*Source: 785049-2014-05-08.xml*