# Mixed Convection Flow along a Stretching Cylinder in a Thermally Stratified Medium

Authors: Swati Mukhopadhyay; Anuar Ishak
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

An analysis for the axisymmetric laminar boundary layer mixed convection flow of a viscous and incompressible fluid towards a stretching cylinder immersed in a thermally stratified medium is presented in this paper. Similarity transformation is employed to convert the governing partial differential equations into highly nonlinear ordinary differential equations. Numerical solutions of these equations are obtained by a shooting method. It is found that the heat transfer rate at the surface is lower for flow in a thermally stratified medium compared to that of an unstratified medium. Moreover, both the skin friction coefficient and the heat transfer rate at the surface are larger for a cylinder compared to that for a flat plate.


---

## Body

## 1. Introduction

Many convection processes occur in environments with stratification. Good examples are closed containers and environmental chambers with heated walls. Also the convection flow associated with heat-rejection systems for long-duration deep ocean power modules where the ocean environment is stratified. Thermally stratified flows are also of great interest in various buoyant flow systems including geothermal systems, geological transport, power plant condensation systems, lake thermohydraulics, and volcanic flows and also in industrial thermal treatment processes. Stratification of the medium may arise due to a temperature variation, which gives rise to a density variation in the medium. This is known as thermal stratification and usually arises due to thermal energy input into the medium from heated bodies and thermal sources. Another situation of interest is the one in which stratification arises due to concentration differences. This is relevant in many natural processes such as transport processes in the sea where stratification exists due to salinity variation. Stratification may arise due to the presence of different fluids so that a stable situation arises when the lighter fluid overlies the denser one. The temperature difference varies from layer to layer and these types of flows have wider applications in oceanography, industry, and agriculture [1].The flow due to a heated surface immersed in a stable stratified viscous fluid has been investigated experimentally and analytically in several studies such as Yang et al. [2], Jaluria and Gebhart [3], Chen and Eichhorn [4], and Ishak et al. [5]. The study of the hydrodynamic flow and heat transfer over a stretching cylinder or a flat plate has gained considerable attention due to its wide applications in industries and important bearings on several technological processes. Crane [6] investigated the flow caused by a stretching plate. Other researchers such as P. S. Gupta and A. S. Gupta [7], Dutta et al. [8], and Chen and Char [9] extended the work of Crane [6] by including the heat and mass transfer analysis under different physical situations. Recently, various aspects of similar problem have been investigated by many authors such as Xu and Liao [10], Cortell [11, 12], Hayat et al. [13], Hayat and Sajid [14], and Ishak et al. [15, 16].Lin and Shih [17, 18] considered the laminar boundary layer and heat transfer along cylinders moving horizontally and vertically with constant velocity and found no similarity solutions due to the curvature effect of the cylinder. Ishak and Nazar [19] showed that the similarity solutions could be obtained by assuming that the cylinder is stretched with a linear velocity in the axial direction. In fact, the study by Ishak and Nazar [19] is an extension of the problem considered by Grubka and Bobba [20] and Ali [21], that is, from a stretching sheet to a stretching cylinder. By considering the effects of mixed convection and thermal stratification parameters, a new dimension is added to the above-mentioned study of Ishak and Nazar [19].Since no attempt has been made to analyse the effects of thermal stratification on boundary layer axisymmetric mixed convection flow along a stretching cylinder, it is considered in this paper. Using a similarity transformation, a third-order ordinary differential equation corresponding to the momentum equation and a second-order ordinary differential equation corresponding to the heat equation are derived. Using a shooting method, up to a desired level of accuracy, the numerical computations are carried out for different values of the dimensionless parameters. The analysis of the results obtained shows that the flow field is influenced appreciably by the mixed convection parameter and the thermal stratification parameter. Estimations of the skin friction and heat transfer coefficients which are very important due to their applications in industries 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. Problem Formulation

Let us consider the steady axisymmetric mixed convection flow of an incompressible viscous fluid along a stretching cylinder embedded in a thermally stratified fluid-saturated medium of variable ambient temperatureT∞(x), where Tw(x)>T∞(x) (heated surface). The continuity, momentum, and energy equations governing such type of flow are
(2.1)∂(ru)∂x+∂(rv)∂r=0,(2.2)u∂u∂x+v∂u∂r=νr∂∂r(r∂u∂r)+gβ(T-T∞),(2.3)u∂T∂x+v∂T∂r=κr∂∂r(r∂T∂r),
where u and v are the components of velocity in the x and r directions, respectively, ν=μ/ρ is the kinematic viscosity, ρ is the fluid density, μ is the coefficient of fluid viscosity, κ is the thermal diffusivity of the fluid, and T is the fluid temperature. It is assumed that the convecting fluid and the medium are in local thermodynamic equilibrium.The boundary conditions for the problem are given by(2.4)u=U(x),v=0,T=Tw(x),atr=R,u→0,T→T∞(x),asr→∞,
where R is the radius of the cylinder, U(x)=U0x/L is the stretching velocity, Tw(x)=T0+b(x/L) is the prescribed surface temperature, and T∞(x)=T0+c(x/L) is the variable ambient temperature. Further, U0, is the reference velocity, T0, the reference temperature, L the characteristic length, and b and c are positive constants.To get similarity solutions of (2.1)–(2.3) subject to the boundary conditions (2.4), we introduce the following similarity transformation:
(2.5)η=r2-R22R(Uνx)1/2,ψ=(Uνx)1/2Rf(η),θ(η)=T-T∞Tw-T0,
where ψ is the stream function defined as u=r-1(∂ψ/∂r) and v=-r-1(∂ψ/∂x) which identically satisfies the continuity equation (2.1). Substituting (2.5) into (2.2) and (2.3) gives
(2.6)(1+2Mη)f′′′+2Mf′′+ff′′-f'2+λθ=0,(1+2Mη)θ′′+2Mθ′+Pr(fθ′-f′θ-f′S)=0,
subject to the boundary conditions
(2.7)f′=1,f=0,θ=1-S,atη=0,f′→0,θ→0,asη→∞,
where prime denotes differentiation with respect to η, S=c/b is the stratification parameter, λ=gβLb/U02 is the mixed convection parameter, and M=(νL/U0R2)1/2 is the curvature parameter. We note that S=0 is for an unstratified environment and M=0 is for a flat plate.One can note that ifM=0 (i.e., R→∞), the problem under consideration (with λ=0, S=0) reduces to the boundary layer flow along a stretching flat plate considered by Ali [21], with m=1 in that paper. Moreover, when M=0, λ=0, and S=0, the analytical solution for the flow field was given by Crane [6] and that of the thermal field was given by Grubka and Bobba [20].


## 3. Results and Discussion

In order to get a clear insight of the physical problem, numerical results are displayed graphically. The results are given through a parametric study showing the influence of several nondimensional parameters, namely, curvature parameter (M), mixed convection parameter (λ), thermal stratification parameter (S), and Prandtl’s number (Pr).Let us first concentrate on the effects of the curvature parameterM on the velocity and temperature distributions. In Figure 1, velocity profiles are shown for different values of M. The velocity curves show that the rate of transport decreases with increasing distance (η) from the surface and vanishes asymptotically at some large distance. The velocity gradient at the surface is larger for larger values of M, which produces larger skin friction coefficient. On the other hand, the temperature is found to decrease with increasing the values of the curvature parameter M (Figure 2). The thermal boundary layer thickness decreases as M increases, which implies an increase in the wall temperature gradient and in turn increases the heat transfer rate at the surface -θ′(0).Figure 1
Variation of velocityf′(η) with η for several values of the curvature parameter M when λ=0.1, S=0, and Pr = 0.7.Figure 2
Variation of temperatureθ(η) with η for several values of the curvature parameter M when λ=0.1, S = 0, and Pr = 0.7.The effects of the stratification parameterS on the temperature and the temperature gradient are exhibited in Figure 3. Due to stratification, the temperature in the boundary layer decreases, which results in a decreasing manner of the temperature gradient (in absolute sense). The thermal boundary layer thickness also decreases with an increase in the stratification parameter S. With the increase in the stratification parameter, the buoyancy factor (Tw-T∞) reduces within the boundary layer. Ambient thermal stratification  causes a significant decrease the local buoyancy levels, which reduces the velocities in the boundary layer. All temperature profiles decay from the maximum value at the wall to zero in the free stream, that is, converge at the outer edge of the boundary layer.Figure 3
Variation of temperatureθ(η) and temperature gradient θ′(η) with η for several values of the stratification parameter S when M = 0.25, λ=0.1, and Pr = 0.7.The effect of the Prandtl number (Pr) on the temperature profiles is exhibited in Figure4. Temperature is found to decrease with increasing Pr. An increase in the Prandtl number reduces the thermal boundary layer thickness. The Prandtl number signifies the ratio of momentum diffusivity to thermal diffusivity. Fluids with lower Prandtl number possess higher thermal conductivities (and thicker thermal boundary layer structures) so that heat can diffuse from the wall faster than higher Pr fluids (thinner boundary layers). Hence, the Prandtl number can control the rate of cooling in conducting flows.Figure 4
Variation of temperatureθ(η) with η for several values of the Prandtl number Pr when M=1, λ=0.1, and S=0.The variation of the velocity gradient at the surfacef′′(0) which is proportional to the skin friction coefficient is presented in Figure 5. It is seen in this figure that the skin friction coefficient is higher (in absolute sense) for a cylinder (with M=0.25) compared to a flat plate (M=0). The negative value of f′′(0) presented in Figure 5 means the solid surface exerts a drag force on the ambient fluid. In the presence of a buoyancy force (λ≠0), the magnitude of the skin friction coefficient |f′′(0)| increases as the stratification parameter S increases.Figure 5
Skin friction coefficientf′′(0) against the stratification parameter S for two values of the curvature parameter M when λ=0.1 and Pr = 0.7.Finally, Figure6 shows that the quantity θ′(0) increases with increasing values of S for both M=0 and M=0.25. Thus, the heat transfer rate at the surface -θ′(0) decreases as S increases. Moreover, the heat transfer rate at the surface is higher for a cylinder compared to a flat plate.Figure 6
Heat transfer coefficientθ′(0) against the stratification parameter S for two values of the curvature parameter M when λ=0.1 and Pr = 0.7.


## 4. Conclusions

The present study gives the numerical solutions for steady mixed convection boundary layer flow and heat transfer along a stretching cylinder in a thermally stratified medium. In the presence of the buoyancy force, the magnitude of the skin friction coefficient increases, but the heat transfer rate at the surface decreases as the stratification parameter increases. In the absence of the curvature parameter, the present problem reduces to that of a flat plate. Both the magnitude of the skin friction coefficient and the heat transfer rate at the surface are higher for a cylinder compared to a flat plate.



---

*Source: 491695-2012-12-18.xml*