# A New Numerical Approach for the Laminar Boundary Layer Flow and Heat Transfer along a Stretching Cylinder Embedded in a

Authors: S. Shateyi; G. T. Marewo
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

The study presents an axisymmetric laminar boundary layer flow of a viscous
incompressible fluid and heat transfer over a stretching cylinder embedded in a
porous medium. A suitable similarity transformation is  employed to transform
the partial differential equations corresponding to the momentum and heat
equations into nonlinear ordinary differential equations. The resultant
ordinary differential equations are then solved using a successive relaxation
method (SRM). The effects of significant parameters on the velocity and temperature
profiles have been analyzed graphically. The obtained results are also compared
with previously published results in some special cases and were found to be
in excellent agreement. The skin friction as well as the heat transfer rate at
the surface are increased as the values of the curvature parameter increase.


---

## Body

## 1. Introduction

The study of boundary layer flow and heat transfer due to stretching flat plates or cylinders has gained considerable attention due to its applications in fibre technology and extrusion processes, as well as theoretical interest. Such applications include the cooling of metallic plates, the boundary layer along a liquid film in condensation processes, boundary layer along material handling conveyers, among others. The rate of heat transfer at the stretching surfaces determines the quality of the final product. Sakiadis [1, 2] pioneered the study of boundary layer flow on a moving continuous solid surface. Thereafter, Crane [3] extended this concept to a stretching sheet with linear surface speed. The study presented an exact solution for the steady two-dimensional flow over a stretching surface in a quiescent fluid. Since then, considerable work has been done by many authors who considered various aspects of this important field (e.g., Laha et al. [4]; Afzal [5]; Prasad et al. [6]; Abel and Mahesha [7]; Abel et al. [8]; Bataller [9]).Flow over cylinders are of two types. They may be considered to be two-dimensional if the body radius is large compared to the boundary layer thickness. On the other hand, if the cylinder is thin/slender, the radius of the cylinder may be of the same order as that of the boundary layer thickness. In such a scenario, the flow may be considered as axisymmetric instead of two-dimensional (Elbarbary and Elgazery [10], Datta et al. [11], among others).Lin and Shih [12, 13] studied laminar boundary layer and heat transfer along horizontally and vertically moving cylinders with constant velocity. The studies found that similarity solutions could not be obtained due to the curvature effect of the cylinder. However, Ishak and Nazar [14] showed that the similarity solutions may be obtained by assuming that the cylinder is stretched with linear velocity in the axial direction and ascertained that their study may be regarded as the extension of the papers by Grubka and Bobba [15] and Ali [16] from a stretching sheet to a stretching cylinder.Mukhopadhyay [17] presented an axisymmetric laminar boundary layer mixed convection flow of a viscous incompressible fluid and heat transfer towards a stretching cylinder embedded in a porous medium. Mukhopadhyay [18] analyzed the flow and thermal characteristics of boundary layer axisymmetric flow and heat transfer along a stretching cylinder in a porous medium. Rangi and Ahmad [19] numerically studied the boundary layer flow of viscous incompressible fluid over a stretching cylinder with variable conductivity. Elbashbeshy et al. [20] discussed laminar boundary layer flow of an incompressible viscous fluid along a stretching horizontal surface embedded in a porous medium in the presence of a heat source or sink with suction/injection. Varjravelu et al. [21] presented numerical solutions to the unsteady convective boundary layer flow of a viscous fluid at a vertical stretching surface with variable transport properties and thermal radiation.Sarkar et al. [22] investigated the buoyancy driven mixed convective flow and heat transfer characteristics of water-based nanofluid past a cylinder. Soid et al. [23] studied a steady two-dimensional magnetohydrodynamics flow and heat transfer over a stretching vertical sheet influenced by radiation and porosity. Elbashbeshy et al. [24] discussed laminar boundary layer flow of an incompressible viscous fluid along a stretching horizontal cylinder in a porous medium. More recently, Chatterjee and Raja [25] investigated numerically the problem of fluid flow and heat transfer subjected to superimposed thermal buoyancy around tandem square cylinders. More recently, Mukhopadhyay and Gorla [26] presented an axisymmetric laminar boundary layer flow of a viscous incompressible fluid and heat transfer towards a stretching cylinder.The present study seeks to extend the work of Rangi and Ahmad [19] to include porosity while also extending the work of Mukhopadhyay [18] to include variable thermal conductivity. The study will be carried out using a newly developed numerical scheme known as the successive relaxation method. This method is based on simple iteration schemes which are formed by reducing the order of the momentum equation which is then followed by linearizing the nonlinear equations. For more details, please see [27]. The rest of the paper is outlined as follows: in Section 2 we give the model formulation of the problem; in Section 3, the successive relaxation method is given in detail; in Section 4, the results and discussion are given; and Section 5 gives the conclusions based on the findings.


## 2. Mathematical Formulation

We consider the steady, axisymmetric boundary layer flow of a viscous and incompressible fluid along a continuously stretching horizontal cylinder of radiusR embedded in a porous medium. We assume that the stretching velocity U(x) and the surface temperature are of the forms U(x)=U0(x/l) and Tw(x)=T∞+T∞(x/l), respectively, where U0 and T0 are constants and T∞ and l are the ambient temperature and characteristic length, respectively. Under these assumptions together with the boundary layer approximations, the equations which model  the current problem under consideration are given as
(1)∂u∂x(ru)+∂v∂r(rv)=0,(2)u∂u∂x+v∂u∂r=νr∂∂r(r∂u∂r)-νKpu,(3)u∂T∂x+v∂T∂r=1r∂∂r(αr∂T∂r),
where u and v are velocity components in the x- and r-directions, respectively, T is the fluid temperature, α is the thermal diffusivity, ν is the kinematic viscosity, and Kp is the permeability parameter. The appropriate boundary conditions are
(4)u=Uw(x),v=0,T=Tw(x),atr=R,u⟶0,T⟶T∞asr⟶∞.


### 2.1. Similarity Transformation

We then transform the momentum and energy equations into the corresponding ordinary differential equations using the following (Ishak and Nazar [14]):
(5)η=r2-R22R(Uwνx)1/2,ψ(x,r)=RνxUwf(η),θ(η)=T-T∞Tw-T∞,
where f(η) is the dimensionless stream function, η is the similarity variable, and θ is the dimensionless temperature. The continuity equation is automatically satisfied through the variables. For liquid metals, it had been observed that the thermal conductivity α varies with temperature in an approximately linear relationship in the range from 0°F to 400°F. To that end, we therefore assume α=α∞(1+ϵθ). Upon introducing the relations (5) into (2)-(3), we obtain the following nonlinear system of ordinary differential equations:
(6)(1+2γη)f′′′+2γf′′+ff′′-f′2-Kf′=0,(7)(1+2γη)θ′′+2γθ′+Prfθ′+ϵ[(1+2γη)(θ′2+θθ′′)+2γθθ′]=0,
subject to the boundary conditions
(8)f(0)=0,f′(0)=1,θ(0)=1,(9)f′⟶0,θ⟶0asη⟶∞,
where prime denotes the differentiation with respect to η,γ=R-1(νl/U0)1/2 is the curvature parameter, K=νl/KpU0 is the permeability parameter, and Pr=μcp/k is the Prandtl number.The physical quantities of interest in engineering are the skin friction coefficientCf and the local Nusselt number Nux which are defined as
(10)Cf=2τwρUw2,Nux=xqwκ(Tw-T∞),
respectively, where τw=μ(∂u/∂r)r=R is the surface shear stress and qw=-κ(∂T/∂r)r=R is the surface heat flux. Upon substituting the similarity transformations (5) into (10) we obtain
(11)12Rex=f′′(0),NuxRex=-θ′(0),
where Rex=Ux/ν is the local Reynolds number.


## 2.1. Similarity Transformation

We then transform the momentum and energy equations into the corresponding ordinary differential equations using the following (Ishak and Nazar [14]):
(5)η=r2-R22R(Uwνx)1/2,ψ(x,r)=RνxUwf(η),θ(η)=T-T∞Tw-T∞,
where f(η) is the dimensionless stream function, η is the similarity variable, and θ is the dimensionless temperature. The continuity equation is automatically satisfied through the variables. For liquid metals, it had been observed that the thermal conductivity α varies with temperature in an approximately linear relationship in the range from 0°F to 400°F. To that end, we therefore assume α=α∞(1+ϵθ). Upon introducing the relations (5) into (2)-(3), we obtain the following nonlinear system of ordinary differential equations:
(6)(1+2γη)f′′′+2γf′′+ff′′-f′2-Kf′=0,(7)(1+2γη)θ′′+2γθ′+Prfθ′+ϵ[(1+2γη)(θ′2+θθ′′)+2γθθ′]=0,
subject to the boundary conditions
(8)f(0)=0,f′(0)=1,θ(0)=1,(9)f′⟶0,θ⟶0asη⟶∞,
where prime denotes the differentiation with respect to η,γ=R-1(νl/U0)1/2 is the curvature parameter, K=νl/KpU0 is the permeability parameter, and Pr=μcp/k is the Prandtl number.The physical quantities of interest in engineering are the skin friction coefficientCf and the local Nusselt number Nux which are defined as
(10)Cf=2τwρUw2,Nux=xqwκ(Tw-T∞),
respectively, where τw=μ(∂u/∂r)r=R is the surface shear stress and qw=-κ(∂T/∂r)r=R is the surface heat flux. Upon substituting the similarity transformations (5) into (10) we obtain
(11)12Rex=f′′(0),NuxRex=-θ′(0),
where Rex=Ux/ν is the local Reynolds number.


## 3. Solution Using the Successive Relaxation Method

The Successive Relaxation Method has been used with great success to solve boundary problems; see, for example, [27–29]. For the current problem, the method begins by letting p=f′ so that p′=f′′, p′′=f′′′, and (6) becomes
(12)(1+2γη)p′′+2γp′+fp′-p2-Kp=0,
but (7) remains unchanged. Proceeding in a manner similar to the Gauss-Seidel method, we replace both (6) and (12) with the iterative scheme defined by
(13)fr+1′=pr(14)(1+2γη)pr+1′′+2γpr+1′+fr+1pr+1′-pr2-Kpr+1=0,(15)(1+2γη)θr+1′′+2γθr+1′+Prfθr+1′+ϵ[(1+2γη)(θr′2+θrθr′′)+2γθrθr′]=0,
where it follows from boundary conditions (8) and (9) that
(16)fr+1(0)=0,(17)pr+1(0)=1,pr+1(∞)=0,(18)θr+1(0)=1,θr+1(∞)=0.
Equations (13)–(15) subject to conditions (16)–(18) will be solved using the Chebyshev Spectral Collocation method. Before we do this, we start by using the domain truncation method to replace the semi-infinite interval [0,∞) with the finite interval [0,L] on the η axis, where L is sufficiently large. For the sake of convenience, we use the transformation
(19)η(ξ)=L2(ξ+1);-1≤ξ≤1
to map interval [0,L] on the η axis to interval [-1,1] on the ξ axis. On the latter, we form a computational grid using the Chebyshev collocation points generated by the formula
(20)ηj=cos(πpN);j=0,1,…,N.
When Chebyshev differentiation [30] is applied on iterative scheme defined by (13)–(18), we end up with the discrete form
(21)A1fr+1=B1,fr+1(ξN)=0,(22)A2pr+1=B2,pr+1(ξN)=1,pr+1(ξ0)=0,(23)A3Θr+1=B3,θr+1(ξN)=1,θr+1(ξ0)=0,
where
(24)A1=D,B1=pr,(25)A2=(1+2γη)D2+diag(2γ+fr+1)D-KI,B2=pr2,(26)A3=(1+2γη)D2+diag(2γ+Prfr+1)D,B3=ϵ[(1+2γη)(Θr′2+ΘrΘr′′)+2γΘΘ′].
As a consequence of the transformation (19), we have D=(2/L)D^ where D^ is the Chebyshev differentiation matrix. Since (21)–(23) are decoupled, they may be solved seperately. This is preceded by applying boundary conditions as illustrated below:
(27)(A10⋯01)(fr+1(ξ0)⋮fr+1(ξN))=(B20),(10⋯0A30⋯01)(θr+1(ξ0)⋮θr+1(ξN))=(0B21),(10⋯0A20⋯01)(pr+1(ξ0)⋮pr+1(ξN))=(0B11).
In order to satisfy the boundary conditions (16), (17), and (18), we choose initial approximations
(28)f0(η)=1-e-η,p0(η)=θ0(η)=e-η,
which if used with iterative scheme defined by (13)–(18) generates subsequent approximations fr, pr, θr for each r=1,2,3,….


## 4. Results and Discussion

In this section, we give the SRM results for the main parameters that have significant effects on the fluid flow velocity and temperature. We remark that all the SRM results presented in this work were obtained usingN=50 collocation points, and also convergence was achieved after as few as five iterations. Also the infinity value (η∞) was taken to be 50. Unless otherwise stated, the default values for the parameters are taken as Pr=0.71; γ=0.1; ϵ=0.2; K=1. In order to validate the numerical method, it was compared with the MATLAB routine bvp4c which is an adaptive Lobatto quadrature iterative scheme. Table 1 presents a comparison between SRM approximate results, the bvp4c results and previous results for Pr=1, and all other parameters being equal to zero. From this table, we clearly observe that the SRM results excellently agree with the bvp4c results up to eight significant figures; whereas previously obtained results are correct only to four significant figures. It can also be remarked that the cpu time for SRM is significantly less than that of bvp4c method. Convergence is achieved by SRM with as few as ten iterations.Table 1
Comparison of the SRM results of-θ′(0) with those obtained by bvp4c as well as previously obtained results.
Pr
-
θ
′
(
0
)
[14]
[15]
[16]
[18]
SRM
b
v
p
4
c
1.0
0.5820
0.5820
0.5801
0.5821
0.58197671
0.58197671Table2 displays the effect of the curvature parameter on the skin friction and Nusselt number when all other parameters are kept constant. Both the temperature and velocity gradients at the surface are larger for larger values of curvature parameter (γ), which produces larger skin friction coefficient and Nusselt number. The effect of the permeability parameter K and thermal conductivity parameter ϵ, on the skin friction and Nusselt number, are, respectively, depicted on Table 3. We observe clearly from this table that the absolute values of f′′(0) are increased by increasing the permeability parameter. Hence, in order to minimize the skin friction value which is usually looked for in an industrial application, one needs to reduce the permeability of the medium. The Nusselt number is reduced as the thermal conductivity of the material increases.Table 2
Comparison of the SRM results of-f′′(0),-θ′(0) with those obtained by bvp4c for different values of the curvature parameter.
γ
-
f
′
′
(
0
)
-
θ
′
(
0
)
b
v
p
4
c
SRM
b
v
p
4
c
SRM
0.0
1.41421356
1.41421356
0.33010615
0.33010615
0.25
1.52316810
1.52316810
0.37955209
0.37955209
0.5
1.62648794
1.62648794
0.42630480
0.42630480Table 3
Comparison of the SRM results of-f′′(0),-θ′(0) with those obtained by bvp4c for different values of the permeability parameter and thermal conductivity parameter, respectively.
K
-
f
′
′
(
0
)
ϵ
-
θ
′
(
0
)
b
v
p
4
c
SRM
b
v
p
4
c
SRM
1.0
1.50188161
1.50188161
0.0
0.43791391
0.4379139
2.0
1.82359290
1.82359290
0.2
0.37955209
0.37955209
3.0
2.09341552
2.09341552
0.4
0.32962783
0.32962783The effects of the pertinent physical parameters in this study on the velocity components and temperature distributions are depicted from Figures1 to 7. Figure 1 displays the effect of the curvature parameter γ on the horizontal velocity profiles. We clearly observe in this figure that the effect of this parameter on the horizontal component is very insignificant within the dynamic region [0,0.7] near the surface. Outside this region, we observe that this velocity component increases as γ increases.Figure 1
Graph of the SRM solutions for the horizontal velocity for different values ofγ.The influence of the permeability parameter on the horizontal velocity is depicted in Figure2. From this figure, we observe that when the permeability parameter is increased, the velocity boundary layer is decreased. Thus, the fluid velocity in the horizontal direction decreases. Figure 3 shows the effect of the curvature of the cylinder on the transverse velocity profiles. We observe that the transverse velocity component is insignificantly affected by curvature parameter within the dynamic region [0,1.3]. Outside this dynamic region, we see that as the curvature of the stretching cylinder increases, this velocity component increases. Physically, as γ→0, the outer surface of the cylinder behaves like a flat surface. Thus, as γ increases, the viscosity effect is reduced due to that the contact area of the surface with fluid tends to the tangential position.Figure 2
Effect of theK on the horizontal velocity profiles.Figure 3
Graph of the SRM solutions of the transverse velocity for different values ofγ.In Figure4, we display the effect of the permeability parameter K on the transverse velocity profiles. The transverse velocity profiles are greatly reduced as the permeability parameter increases. Figure 5 shows the effect of the curvature parameter on the temperature profiles. It can be clearly observed in this figure that there is no significant effect within the dynamic region less than a unit. However, after this region the temperature is enhanced by increasing the values of the curvature parameter. Figure 6 displays the effect of varying thermal conductivity on the temperature distribution. As expected, increasing thermal conductivity (or a higher viscosity) results in the thinning of the thermal boundary layer and hence higher heat transfer rate at the surface. This is clearly depicted in Figure 7.Figure 4
Effect of theK on the v-velocity profiles.Figure 5
Graph of the SRM solutions of the temperature profiles for different values ofγ.Figure 6
Effect of theϵ on the temperature profiles.Figure 7
Effect of thePr on the temperature profiles.


## 5. Conclusion

We have numerically analyzed the problem of laminar boundary layer flow and heat transfer along a stretching cylinder embedded in a porous medium with variable conductivity applying a recently developed numerical method. The system of partial differential equations governing the current problem was transformed into a set of ordinary differential equations by using appropriate similarity transformations. The investigation observed that the curvature of the stretching cylinder has a very significant effect on both the velocity and temperature fields. Both the skin friction coefficient and local Nusselt number increase as the curvature increases. The flow properties were found to be significantly influenced by the permeability parameter. The present study also shows that the SRM is a very reliable, easy, and accurate method which we trust that it can be used to solve even more complex and complicated systems. It is hoped that the current findings can be used as the basis for many scientific and engineering applications.



---

*Source: 576453-2013-11-21.xml*