# Analysis of Boundary Layer Flow and Heat Transfer along a Stretching Cylinder in a Porous Medium

Authors: Swati Mukhopadhyay
Journal: ISRN Thermodynamics
 (2012)

---

## Abstract

This paper presents an axi-symmetric laminar boundary layer flow of a viscous incompressible fluid and heat transfer towards a stretching cylinder embedded in a porous medium. The partial differential equations corresponding to the momentum and heat equations are converted into highly nonlinear ordinary differential equations with the help of similarity transformations. Numerical solutions of these equations are obtained by shooting method. It is found that the velocity decreases with increasing permeability parameter. The skin friction as well as the heat transfer rate at the surface is larger for a cylinder compared to a flat plate.


---

## Body

## 1. Introduction

The study of hydrodynamic flow and heat transfer over a stretching cylinders or flat plates has gained considerable attention due to its applications in industries and important bearings on several technological processes. Crane [1] investigated the flow caused by the stretching of a sheet. Many researchers such as  Gupta and Gupta [2], Dutta et al. [3], Chen and Char [4] extended the work of Crane [1] by including the effect of heat and mass transfer analysis under different physical situations. Recently, various aspects of such problem have been investigated by many authors such as Xu and Liao [5], Cortell [6, 7], Hayat et al. [8] and Hayat and Sajid [9].Flow over cylinders is considered to be two-dimensional if the body radius is large compared to the boundary layer thickness. On the other hand for a thin or slender cylinder, the radius of the cylinder may be of the same order as that of the boundary layer thickness. Therefore, the flow may be considered as axi-symmetric instead of two-dimensional [10–12]. The study of steady flow in a viscous and incompressible fluid outside of a stretching hollow cylinder in an ambient fluid at rest has been done by Wang [13]. The effect of slot suction/injection over a thin cylinder as studied by Datta et al. [10] and Kumari and Nath [12] may be useful in the cooling of nuclear reactors during emergency shutdown, where a part of the surface can be cooled by injecting a coolant (Ishak et al. [14]). Lin and Shih [15, 16] considered the laminar boundary layer and heat transfer along horizontally and vertically moving cylinders with constant velocity and found that the similarity solutions could not be obtained due to the curvature effect of the cylinder. Ishak and Nazar [17] showed that the similarity solutions may be obtained by assuming that the cylinder is stretched with linear velocity in the axial direction and claimed that their study may be regarded as the extension of the papers by Grubka and Bobba [18] and Ali [19], from a stretching sheet to a stretching cylinder.The study of hydrodynamic flow and heat transfer in porous medium becomes much more interesting due to its vast applications on the boundary layer flow control. Heat removal from nuclear fuel debris, underground disposal of radioactive waste material, storage of food stuffs and exothermic and/or endothermic chemical reactions and dissociating fluids in the packed-bed reactors, and so forth are some porous media applications. It is well known that Darcy’s law is an empirical formula relating the pressure gradient, the bulk viscous fluid resistance and the gravitational force for a forced convective flow in a porous medium. Deviations from Darcy’s law occur when the Reynolds number based on the pore diameter is within the range of 1 to 10 (Ishak et al. [20]).No attempt has been made yet to analyze the flow and thermal characteristics of boundary layer axi-symmetric flow and heat transfer along a stretching cylinder in a porous medium. Therefore, an attempt is made to study the steady forced convection flow and heat transfer past a stretching cylinder placed in a fluid-saturated porous medium using the Darcy model. Using similarity transformation, a third order ordinary differential equation corresponding to the momentum equation and a second-order ordinary differential equation corresponding to heat equation are derived. Using shooting method, numerical calculations up to desired level of accuracy were carried out for different values of dimensionless parameters of the problem under consideration for the purpose of illustrating the results graphically. The results obtained are then compared with those of Grubka and Bobba [18], Ali [19], and Ishak and Nazar [17] who reported the results for some special case of the presented study. The analysis of the results obtained shows that the flow field is influenced appreciably by the permeability parameter. Estimation of skin friction and heat transfer coefficient which are very important from the industrial application point of view are also presented in the analysis. It is hoped that the results obtained will not only provide useful information for applications but also serve as a complement to the previous studies.


## 2. Equations of Motion

Consider the steady axi-symmetric flow of an incompressible viscous fluid along a stretching cylinder embedded in a porous medium. The continuity, momentum, and energy equations governing such type of flow are written as(1)∂(ru)∂x+∂(rv)∂r=0,(2)u∂u∂x+v∂u∂r=vr∂∂r(r∂u∂r)-vku,u∂T∂x+v∂T∂r=κr∂∂r(r∂T∂r),
where u  and v are the components of velocity, respectively, in the x and r directions, v=μ/ρ is the kinematic viscosity, ρ is the fluid density (assumed constant), μ is the coefficient of fluid viscosity, k is the porosity of the medium, κ is the thermal diffusivity of the fluid, and T is the fluid temperature.


### 2.1. Boundary Conditions

The appropriate boundary conditions for the problem are given by(3)u=U(x),v=0,T=Tw(x),atr=R,u⟶0,T⟶T∞,asr⟶∞.
Here, U(x)=U0(x/L)  is the stretching velocity, Tw(x)=T∞+T0(x/L)N is the prescribed surface temperature, U0, T0 are the reference velocity and temperature, respectively, T∞ is the ambient temperature, L is the characteristic length, and N is the temperature exponent.


### 2.2. Method of Solution

The continuity equation is automatically satisfied by the introduction of stream functionψ as u=(1/r)(∂ψ/∂r),v=-(1/r)(∂ψ/∂x). Introducing the similarity variables as
(4)η=r2-R22R(Uvx)1/2,ψ=(Uvx)1/2Rf(η),θ(η)=T-T∞Tw-T∞,
and upon substitution of (4) in (2) and (3), the governing equations and the boundary conditions reduce to(5)(1+2Mη)f///+2Mf//+ff//-f/2-k1f/=0,(1+2Mη)θ//+2Mθ/+Pr(fθ/-Nf/θ)=0,(6)f/=1,f=0,θ=1atη=0,(7)f/⟶0,θ/⟶0asη⟶∞,
where the prime denotes differentiation with respect toη, k1=vL/U0k is the permeability parameter of the porous medium and  M=(vL/U0R2)1/2 is the curvature parameter. The case of nonporous medium is recovered for k1=0. k1-1 will reflect the effect of Darcian flow on the present problem.One can note that ifM=0 (i.e., R→∞), the problem under consideration (with k1=0) reduces to the boundary layer flow along a stretching flat plate considered by Ali [19], with m=1 in that paper. Moreover, when M=0 (stretching flat plate) subjected to (6) with k1=0 (i.e., for non-porous medium), the analytical solutions of (5) are given by Crane [1] and Grubka and Bobba [18], respectively.


## 2.1. Boundary Conditions

The appropriate boundary conditions for the problem are given by(3)u=U(x),v=0,T=Tw(x),atr=R,u⟶0,T⟶T∞,asr⟶∞.
Here, U(x)=U0(x/L)  is the stretching velocity, Tw(x)=T∞+T0(x/L)N is the prescribed surface temperature, U0, T0 are the reference velocity and temperature, respectively, T∞ is the ambient temperature, L is the characteristic length, and N is the temperature exponent.


## 2.2. Method of Solution

The continuity equation is automatically satisfied by the introduction of stream functionψ as u=(1/r)(∂ψ/∂r),v=-(1/r)(∂ψ/∂x). Introducing the similarity variables as
(4)η=r2-R22R(Uvx)1/2,ψ=(Uvx)1/2Rf(η),θ(η)=T-T∞Tw-T∞,
and upon substitution of (4) in (2) and (3), the governing equations and the boundary conditions reduce to(5)(1+2Mη)f///+2Mf//+ff//-f/2-k1f/=0,(1+2Mη)θ//+2Mθ/+Pr(fθ/-Nf/θ)=0,(6)f/=1,f=0,θ=1atη=0,(7)f/⟶0,θ/⟶0asη⟶∞,
where the prime denotes differentiation with respect toη, k1=vL/U0k is the permeability parameter of the porous medium and  M=(vL/U0R2)1/2 is the curvature parameter. The case of nonporous medium is recovered for k1=0. k1-1 will reflect the effect of Darcian flow on the present problem.One can note that ifM=0 (i.e., R→∞), the problem under consideration (with k1=0) reduces to the boundary layer flow along a stretching flat plate considered by Ali [19], with m=1 in that paper. Moreover, when M=0 (stretching flat plate) subjected to (6) with k1=0 (i.e., for non-porous medium), the analytical solutions of (5) are given by Crane [1] and Grubka and Bobba [18], respectively.


## 3. Numerical Method for Solution

The above equations (5) along with boundary conditions are solved by converting them to an initial value problem. We set(8)f/=z,z/=p,p/=[z2+k1z-fp-2Mp](1+2Mη),θ/=q,q/=-[Pr⁡(fq-Nzθ)+2Mq](1+2Mη),
with the boundary conditions
(9)f(0)=0,f/(0)=1,θ(0)=1.
In order to integrate (8) as initial value problems, one requires a value for p(0), that is, f//(0) and a value for q(0) that is, θ/(0) but no such values are given at the boundary. The suitable guess values for f//(0) and θ/(0) are chosen and then the integration is carried out. Comparing the calculated values for f/ and θ at η=10 (say) with the given boundary conditions f/(10)=0 and θ(10)=0 and adjusting the estimated values, f//(0) and θ/(0), a better approximation for the solution is given.Taking the series of values forf//(0)  and θ/(0) and applying the fourth-order classical Runge-Kutta method with step size h=0.01, the above procedure is repeated until the results up to the desired degree of accuracy (10-5) are obtained.


## 4. Results and Discussions

In order to analyse the results, numerical computation has been carried out using the method described in the previous section for various values of the curvature parameter (M), permeability parameter (k1), temperature exponent (N) and Prandtl number (Pr). For illustrations of the results, numerical values are plotted in Figures 1(a) to 5(b).(a) Variation of velocityf/(η) with η for several values of  curvature parameter M of the stretching cylinder. (b) Variation of temperature θ(η) with η for several values of  curvature parameter M of the stretching cylinder.
(a)(b)For the verification of accuracy of the applied numerical scheme, a comparison of the present results corresponding to the heat transfer coefficient[-θ/(0)] for k1=0 (i.e., in case of nonporous medium) and M=0 (i.e., for stretching flat plate) with the available published results of Ishak and Nazar [17], Grubka and Bobba [18], Ali [19] is made and presented in Table 1. The results are found in excellent agreement.Table 1
Values ofNuxRex-1/2=-θ/(0) for several values of temperature exponent N for flat plate (M=0) in nonporous medium (k1=0) and Pr=1.
NIshak and Nazar [17]Grubka and Bobba [18]Ali [19]Present study00.58200.58200.58010.582111.00001.00000.99611.000021.33331.33331.32691.3332Let us first concentrate on the effects of curvature parameterM on velocity distribution in presence of porous medium. In Figure 1(a), horizontal velocity profiles are shown for different values of M. The horizontal velocity curves show that the rate of transport decreases with the increasing distance (η) of the sheet. In all cases, the velocity vanishes at some large distance from the sheet (at η=10). The velocity increases with increasing values of M. The velocity gradient at the surface is larger for larger values of M, which produces larger skin friction coefficient.Effects of curvature parameter on the temperature distribution are presented in Figure1(b). Temperature is found to decrease with the increasing curvature parameter M. The thermal boundary layer thickness decreases as M increases, which implies increasing in the wall temperature gradient and in turn increase the surface heat transfer rate. Hence, the Nusselt number increases as M increases.Now, the velocity profiles are presented for the variation of permeability parameter for flat plate and stretching cylinder, respectively. Figures2(a) and 2(c) demonstrate the effects of permeability parameter (k1) for a stretching flat plate (i.e., for M=0) on velocity and temperature respectively. With the increasing k1, the horizontal velocity is found to decrease (Figure 2(a)) but the temperature increases in this case. This feature prevails up to certain heights and then the process is slowed down. It is noted that k1 has a substantial effect on the solutions. It is obvious that the presence of porous medium causes higher restriction to the fluid, which reduces the fluid velocity (Figure 2(a)) and enhances the temperature (Figure 2(b)). The porous medium presents resistance to the flow, thus, the flow becomes slower. Therefore, as the inverse Darcy number (k1-1) increases, the resistance due to the porous medium increases and the surface velocity gradient increases. It is seen from the figure that the skin friction increases monotonically.k1=0 corresponds to the case of non-porous medium. Figures 2(b) and 2(d) present the effects of permeability parameter k1 on velocity and temperature for the stretching cylinder (M>0). Here also velocity decreases with k1 whereas the temperature increases with increasing k1 (Figure 2(d)).(a) Variation of velocityf/(η) with η for several values of  Permeability parameter k1 for flat plate with M=0. (b) Variation of velocity f/(η) with η for several values of permeability parameter k1 for stretching cylinder with M=1. (c) Variation of temperature θ(η) with η for several values of  permeability parameter k1 for flat plate with M=0. (d) Variation of temperature θ(η) with η  for several values of permeability parameter k1 for stretching cylinder with M=0.25.
(a)(b)(c)(d)Figures3(a) and 3(b) represent the temperature and temperature gradient profiles for variable temperature exponent N for stretching flat plate (M=0) in the presence of porous medium. It is seen that temperature increases with increasing values of N and temperature overshoot is noted near the plate for N>0. The temperature gradient initially increases with N but after a certain distance η from the wall it decreases. Far away from the wall, such feature is smeared out. It is observed that, for higher positive values of N(=2,3), initially θ/(η) is positive but finally it becomes negative. Temperature increases with N but no temperature overshoot is noted for stretching cylinder (M=1) (Figure 3(c)). For stretching cylinder, the temperature gradient is negative for all values of N (Figure 3(d)).(a) Variation of temperatureθ(η) with η for several values of temperature exponent N for flat plate with M=0. (b) Variation of temperature gradient θ/(η) with η for several values of temperature exponent, N for flat plate with M=0. (c) Variation of temperature θ(η) with η for several values of temperature exponent N for stretching cylinder with M=1. (d) Variation of temperature gradient θ/(η) with η for several values of temperature exponent N for stretching cylinder with M=1.
(a)(b)(c)(d)Figures4(a) and 4(b) present the behaviour of skin friction and heat transfer coefficients with the permeability parameter k1 of the porous medium for two values of curvature parameter. Magnitude of the skin friction coefficient increases with increasing permeability parameter k1 and also with the curvature parameter M. This feature also supports the earlier findings in Figures 1(a), 2(a) and 2(b). From the figure, it is very clear that shear stress at the wall is negative here. Physically, negative sign of f//(0) implies that surface exerts a dragging force on the fluid and positive sign implies the opposite. This is consistent with the present case as a stretching cylinder which induces the flow is considered here. From Figure 4(b) it is very clear that the magnitude of heat transfer coefficient decreases with permeability parameter k1 as well as with the curvature parameter M. This is consistent with the findings in Figures 1(b), 2(c) and 2(d).(a) Variation of skin friction coefficientf//(0) with permeability parameter k1 for two values of curvature parameter M. (b) Variation of heat transfer coefficient θ/(0) with permeability parameter k1 for two values of curvature parameter M.
(a)(b)(a) Variation of temperatureθ(η) with η for several values of  Prandtl number Pr for flat plate with M=0. (b) Variation of temperature θ(η) with η for several values of Prandtl number Pr for stretching cylinder with M=1.
(a)(b)It is noted that temperature decreases with increasing Pr. An increase in Prandtl number reduces the thermal boundary layer thickness. It is also observed that the effects of Pr is much more prominent for flat plate (Figure5(a)) compared to stretching cylinder (Figure 5(b)). Prandtl number signifies the ratio of momentum diffusivity to thermal diffusivity. Fluids with lower Prandtl number will possess higher thermal conductivities (and thicker thermal boundary layer structures) so that heat can diffuse from the wall, faster than for higher Pr fluids (thinner boundary layers). Hence, Prandtl number can be used to increase the rate of cooling in conducting flows.



---

*Source: 704984-2012-06-13.xml*