# Natural Convection of Viscoelastic Fluid from a Cone Embedded in a Porous Medium with Viscous Dissipation

Authors: Gilbert Makanda; O. D. Makinde; Precious Sibanda
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

We study natural convection from a downward pointing cone in a viscoelastic fluid
embedded in a porous medium. The fluid properties are numerically computed for different viscoelastic, porosity, Prandtl and Eckert numbers. The governing partial differential equations are converted to a system of fourth order ordinary
differential equations using the similarity transformations and then solved together by
using the successive linearization method (SLM). Many studies have been carried out on
natural convection from a cone but they did not consider a cone embedded in a porous
medium with linear surface temperature. The results in this work are validated by the
comparison with other authors.


---

## Body

## 1. Introduction

Natural convection of viscoelastic fluid in a porous medium with viscous dissipation is the transfer of heat due to density differences caused by temperature gradients through a permeable medium and heat generated due to the interaction of fluid molecules is considered. There are examples in practical application such as thermal insulation, extraction of petroleum resources and the so-called fracking, metal processing, performance of lubricants, application of paints, and extrusion of plastic sheets. The study of second grade fluids has been studied but there is no single constitutive equation that can fully describe non-Newtonian fluids [1]; due to this fact many authors did not consider the appropriate constitutive energy equation for second grade fluids.Natural convection on a cone geometry has been studied by among others Alim et al. [2], Awad et al. [3], Cheng [4, 5], and Kairi and Murthy [6]. Studies have been done on other geometries such as flow over a flat plate, cylinders, vertical surfaces, stretching sheets, and inclined surfaces by, among others, Abbas et al. [7] who considered unsteady second grade fluid flow on an unsteady stretching sheet; they did not consider the energy equation mainly due to difficulties in its characterization. Anwar et al. [8] studied mixed convection boundary layer flow of a viscoelastic fluid over a horizontal circular cylinder; they solved the fourth order ordinary differential equations by considering the insufficiency of the boundary conditions by taking the zeroth, first, and second order of the viscoelastic parameter and coming up with three systems of ordinary differential equations. Cortell [9] investigated flow and heat transfer of a viscoelastic fluid over a stretching sheet. Damseh et al. [10] studied the transient mixed convection flow of a second grade viscoelastic fluid over a vertical surface. They used McCormack’s method to solve their differential equations. Hayat et al. [11] studied mixed convection in a stagnation point flow adjacent to a vertical surface in a viscoelastic fluid.The model in this work has been originally developed from the work of Ece [5] who studied heat and mass transfer from a downward pointing cone in a Newtonian fluid. In this paper the work of Ece [5] is extended to take into account the flow of a second grade fluid in a porous medium and the effect of viscous dissipation is considered. Several other studies have been done in natural convection in a viscoelastic fluid by among others Hsiao [12] who studied mixed convection for viscoelastic fluid past a porous wedge. Kasim et al. [13] investigated free convection boundary layer flow of a viscoelastic fluid in the presence of heat generation. Massoudi et al. [14] studied natural convection flow of generalized second grade fluid between two vertical walls. Olajuwon [15] studied the convection heat and mass transfer in a hydromagnetic flow of a second grade fluid in the presence of thermal radiation and thermal diffusion; it was shown that increasing the second grade parameter causes reduction in the rate of the fluid flow and mass transfer, but heat transfer increases. Sajid et al. investigated fully developed mixed convection flow of a viscoelastic fluid between permeable parallel vertical plates [16].Studies for viscous dissipation in a second grade fluid have been done by many authors but some assumed that fluids are more viscous than elastic resulting in the energy equation without the elastic term. Viscous dissipation has been studied by among others Subhas Abel et al. [17] who studied viscoelastic MHD flow and heat transfer over a stretching sheet with viscous and ohmic dissipations and [18] in which a Newtonian fluid was considered. The viscous dissipation term which they used in [17] assumes that the fluid is more viscous in nature than elastic. Jha [19] investigated the effects of viscous dissipation on natural convection flow between parallel plates with time periodic boundary conditions. Chen [20] studied the analytic solution of MHD flow and heat transfer for two types of viscoelastic fluids over a stretching sheet with energy dissipation, internal heat source, and thermal radiation. Cortell [21] worked on viscous dissipation and thermal radiation effects on the flow and heat transfer of a power law fluid past an infinite porous plate. Hsiao [22] investigated multimedia physical feature for unsteady MHD mixed convection viscoelastic fluid over a vertical stretching sheet with viscous dissipation. Kameswaran et al. [23] studied hydromagnetic nanofluid flow due to a stretching sheet or shrinking sheet with viscous dissipation and chemical reaction effects.Studies have been done in porous media by among others Awad et al. [3, 4, 6, 24] and Singh and Agarwal [25] who studied heat transfer in a second grade fluid over an exponentially stretching sheet through porous medium with thermal radiation and elastic deformation under the effect of magnetic field.An investigation of available literature shows that, to the best of our knowledge, no analysis has been done on natural convection of a viscoelastic fluid embedded in a porous medium with viscous dissipation under the given boundary conditions. The study takes into consideration a temperature that changes linearly along the surface of the cone (see Ece [5]).


## 2. Mathematical Formulation

A cone in a viscoelastic fluid embedded in a porous medium is heated and maintained at a linearly changing temperatureT (>T∞), and the ambient conditions are maintained at T∞; the fluid has a constant viscosity ν. The vertex angle of the cone is 2ϕ. The velocity components u and v are in the directions of x and y, respectively, with the x-axis being inclined at an angle ϕ to the vertical. A sketch of the system and coordinate axis is illustrated in Figure 1.Figure 1
Physical model and coordinate system.The governing equations in this buoyant-driven flow are given by(1)∂∂x(ru)+∂∂y(rv)=0,u∂u∂x+v∂u∂y=ν∂2u∂y2-νKu-ko{u∂3u∂x∂y2+v∂3u∂y3ν∂2u∂y2-νKu-ko{l+∂u∂x∂2u∂y2-∂2u∂x∂y∂u∂y}+gβ(T-T∞)cosϕ,u∂T∂x+v∂T∂y=α∂2T∂y2+νCp(∂u∂y)2+k0ρCp(u∂2u∂x∂y∂u∂y+v∂2u∂y2∂u∂y),
where r=xsinϕ, g is the acceleration due to gravity, ν is the kinematic viscosity for the fluid, ko is the non-Newtonian parameter of the viscoelastic fluid, β is the coefficient of thermal expansion, α is the thermal diffusivity, Cp is the specific heat capacity for the fluid, ρ is the density of the fluid, and K is the permeability coefficient of the porous medium. The boundary conditions are given as
(2)u=v=0,T=Tw(x)=T∞+A(xL)aty=0,∂u∂y,u⟶0,T⟶T∞,asy⟶∞,
where A>0 is a constant, L>0 is the characteristic length, and the subscript ∞ refers to the ambient condition.We introduce the nondimensional variables:(3)X=xL,Y=Gr1/4yL,R=rL,U=uU0,V=Gr1/4vU0,T¯=T-T∞Tw-T∞,Gr=(U0Lν)2,
where U0=[gβcosϕL(Tw-T∞)]1/2. Using (3) in (1) gives the following equations:
(4)∂∂X(RU)+∂∂Y(RV)=0,U∂U∂X+V∂U∂Y=∂2U∂Y2-νUK-Λ{U∂3U∂X∂Y2+V∂3U∂Y3+∂U∂X∂2U∂Y2=∂2U∂Y2-νUK-Λ{l-∂2U∂X∂Y∂U∂Y}+T¯,U∂T¯∂X+V∂T¯∂Y=1Pr∂2T¯∂Y2+Ec(∂U∂Y)2+ΛEc(U∂2U∂X∂Y∂U∂Y+V∂2U∂Y2∂U∂Y),
where R=Xsinϕ, Λ=(k0U0/νL) is the viscoelastic parameter known as the Deborah number, Gr is the Grashof number, Pr=ν/α is the Prandtl number, and Ec=(U02/CpA) is the Eckert number. The corresponding boundary conditions are given as
(5)U=V=0,T¯=XatY=0,∂U∂Y,U⟶0,T¯⟶0asY⟶∞.
We now introduce the following similarity variables ψ=XRf(Y) and T¯=Xθ(Y) defined by
(6)U=1R∂ψ∂Y,V=-1R∂ψ∂X.Substituting (6) and the similarity variables in (4) gives the following ordinary differential equations:
(7)f′′′+2ff′′-(f′)2+θ-γf′-Λ(2f′f′′′-2ffiv-(f′′)2)=0,(8)θ′′+Pr(2fθ′-f′θ)+PrEcf′′2+ΛPrEc(f′f′′2-ff′′f′′′)=0.
With boundary conditions,
(9)f(0)=f′(0)=0,θ(0)=1,(10)f′(∞)⟶0,f′′(∞)⟶0,θ(∞)⟶0.
It is of interest to discuss the skin friction and the heat transfer coefficient in this context. The shear stress at the surface of the cone is defined as (see Olajuwon [15])
(11)τw=μ[∂u∂y]y=0+k0[u∂2u∂x∂y-2∂u∂x∂u∂y]y=0,
where μ is the coefficient of viscosity. The skin friction is defined as
(12)cf=τw(1/2)ρU∞2,cf=2XGr1/4f′′(0)(1+3Λf′(0)).
The skin friction coefficient can be expressed as
(13)CfGr1/42X=f′′(0).
The heat transfer rate at the surface of the cone is given by
(14)qw=-kX[∂T∂y]y=0.
The Nusselt number can be expressed as
(15)Nu=Lqwk(Tw-T∞).
Using the nondimensional variables (9)-(10), the dimensionless wall heat rate is given by
(16)NuGr-1/4=-θ′(0).


## 3. Method of Solution

In this study, (7)–(10) were solved using the successive linearization method. The inclusion of the non-Newtonian term brings about the fourth order ordinary differential equation for the momentum equation. The given boundary conditions are insufficient to obtain a unique solution. To overcome this problem the system is decomposed into the zeroth, first, and second order systems of the viscoelastic parameter. Subhas Abel et al. [17] showed that if this method is applied small values of the viscoelastic parameter can be used without difficulty in convergence. It is also noticed in this study that the direct application of the successive linearization method has difficulties in convergence for small values of the viscoelastic parameter. Anwar et al. [8] also confirmed the same observation and solved a system of differential equations simultaneously and obtained better convergence for small values of the viscoelastic parameter. In this work we solve the system using the successive linearization method. To solve the equations we seek the series solution of the form
(17)f(y)=f0(y)+Λf1(y)+Λ2f2(y)+⋯,θ(y)=θ0(y)+Λθ1(y)+Λ2θ2(y)+⋯.
The skin friction can be computed using
(18)f′′(0)=f0′′(0)+Λf1′′(0)+Λ2f2′′(0)+⋯.Then substituting (17) into the system (7)–(10). We then take the zeroth, first, and second order of the viscoelastic parameter Λ. We obtain the following system.Zeroth order:(19)f0′′′+2f0f0′′-f0′2+θ0-γf0′=0,(20)θ′′+2Prf0θ0′-Prf0′θ0+PrEcf0′′2=0,(21)f0(0)=0,f0′(0)=0,θ0(0)=1,(22)f0′(∞)=0,θ0(∞)=0.First order:(23)f1′′′+2f0f1′′+2f1f0′′-2f0′f1′+θ1-γf1′-2f0′f0′′′+2f0f0iv+f0′′2=0,(24)θ1′′+2Prf0θ1+2Prf1θ0′-Prf0′θ1-Prf1′θ0+2PrEcf0′′f1′′+PrEcf0′f0′′2-PrEcf0f0′′f0′′′=0,(25)f1(0)=0,f1′(0)=0,f1′(∞)=0,(26)θ1(0)=0,θ1(∞)=0.Second order:(27)f2′′′+2f0f2′′+2f1f1′′+2f2′f0′′-2f0′f2′-f1′2+θ2-γf2′-2f0′f1′′′-2f1′f0′′′+2f0f1iv+2f1f0iv+2f0′′f1′′=0,(28)θ2′′+2Prf0θ2′+2Prf1θ1′+2Prf2θ0′-Prf0′θ2-Prf1′θ1+PrEc(2f0′′f2′′+f1′′2+2f0′f0′′f1′′+f1′f0′′2+f2′f0′′2+PrEc(lllll+f0f0′′f1′′′+f0f1′′f0′′′+f1f0′′f0′′′)=0,(29)f2(0)=0,f2′(0)=0,f2′(∞)=0,(30)θ2(0)=0,θ2(∞)=0.
The functions in the system (19)–(30) may be expanded in series form as
(31)f0(y)=f0i(y)+∑m=0i-1‍f0m(y),θ0(y)=θ0i(y)+∑m=0i-1‍θ0m(y),f1(y)=f1i(y)+∑m=0i-1‍f1m(y),θ1(y)=θ1i(y)+∑m=0i-1‍θ1m(y),f2(y)=f2i(y)+∑m=0i-1‍f2m(y),θ2(y)=θ2i(y)+∑m=0i-1‍θ2m(y),
where f0i, f1i, and f2i and θ0i, θ1i, and θ2i (i=1,2,3,…) are unknown functions and f0m, f1m, and f2m and θ0m, θ1m, and θ2m are approximations that are found by successively solving the linear part of equations that are obtained after substituting (31) into system (19)–(30). These linear equations have the form
(32)f0i′′′+a01,i-1f0i′′+a02,i-1f0i′+a03,i-1f0i+a04,i-1θ0i=r01,i-1,(33)θ0i′′+b01,i-1θ0i′+b02,i-1θ0i+b03,i-1f0i′′+b04,i-1f0i′+b05,i-1f0i=r02,i-1,(34)f0i(0)=0,f0i′(0)=0,θ0i(0)=0,(35)f0i′(∞)=0,θ0i(∞)=0,f0i′′(∞)=0,(36)f1i′′′+a11,i-1f1i′′+a12,i-1f1i′+a13,i-1f1i+a14,i-1θ1i=r11,i-1,(37)θ1i′′+b11,i-1θ1i′+b12,i-1θ1i+b13,i-1f1i′′+b14,i-1f1i′+b15,i-1f1i=r12,i-1,(38)f1i(0)=0,f1i′(0)=0,θ1i(0)=0,(39)f1i′(∞)=0,θ1i(∞)=0,f1i′′(∞)=0,(40)f2i′′′+a21,i-1f2i′′+a22,i-1f2i′+a23,i-1f2i+a24,i-1θ2i=r21,i-1,(41)θ2i′′+b21,i-1θ2i′+b22,i-1θ2i+b23,i-1f2i′′+b24,i-1f2i′+b25,i-1f2i=r22,i-1,(42)f2i(0)=0,f2i′(0)=0,θ2i(0)=0,(43)f2i′(∞)=0,θ2i(∞)=0,f2i′′(∞)=0.
The coefficients ajk,i-1, bjk,i-1 (j=0,1,2,k=1,…,5), rj1,i-1, and rj2,i-1 are defined as
(44)a01,i-1=a11,i-1=a21,i-1=2∑m=0i-1‍f0m,a02,i-1=a12,i-1=a22,i-1=-(2∑m=0i-1‍f0m′+γ),a03,i-1=a13,i-1=a23,i-1=2∑m=0i-1‍f0m′′,a04,i-1=a14,i-1=a24,i-1=I,b01,i-1=b11,i-1=b21,i-1=2Pr∑m=0i-1‍f0m,b02,i-1=b12,i-1=b22,i-1=-Pr∑m=0i-1‍f0m′,b03,i-1=b13,i-1=b23,i-1=PrEc∑m=0i-1‍f0m′′,b04,i-1=b14,i-1=b24,i-1=-Pr∑m=0i-1‍θ0m,b05,i-1=b15,i-1=b25,i-1=2Pr∑m=0i-1‍θ0m′,r01,i-1=-[∑m=0i-1‍f0m′′′+2∑m=0i-1‍f0m∑m=0i-1‍f0m′′-(∑m=0i-1‍f0m′)2r01,i-1=-[l-∑m=0i-1‍θ0m-γ∑m=0i-1‍f0m′(∑m=0i-1‍f0m′)2],r02,i-1=-[(∑m=0i-1‍f0m′′)2∑m=0i-1‍θ0m′′+2Pr∑m=0i-1‍f0m∑m=0i-1‍θ0m′r02,i-1=-[l-Pr∑m=0i-1‍f0m′∑m=0i-1‍θ0m+PrEc(∑m=0i-1‍f0m′′)2],r11,i-1=-[∑m=0i-1‍f1m′′′+4∑m=0i-1‍f0m∑m=0i-1‍f1m′+4∑m=0i-1‍f1m∑m=0i-1‍f0m′′r11,i-1=-[l-6∑m=0i-1‍f0m′′∑m=0i-1‍-γ∑m=0i-1‍f0m′+3(∑m=0i-1‍f0m′′)2r11,i-1=-[l+6∑m=0i-1‍f0m∑m=0i-1‍f0miv-6∑m=0i-1‍f0m′∑m=0i-1‍f0m′′′+∑m=0i-1‍],r12,i-1=-[∑m=0i-1‍θ1m′′+4Pr∑m=0i-1‍f0m∑m=0i-1‍θ1m′r12,i-1=-[l-2Pr∑m=0i-1‍f0m′∑m=0i-1‍θ1m-2Pr∑m=0i-1‍f1m′∑m=0i-1‍θ0mr12,i-1=-[l+4Pr∑m=0i-1‍f1m∑m=0i-1‍θ0m′+4PrEc∑m=0i-1‍f0m′′∑m=0i-1‍f1m′′r12,i-1=-[l+4PrEc∑m=0i-1‍f0m′(∑m=0i-1‍f0m)2r12,i-1=-[l-4PrEc∑m=0i-1‍f0m∑m=0i-1‍f0m′′∑m=0i-1‍f0m′′′],r21,i-1=-[∑m=0i-1‍f2′′′+4∑m=0i-1‍f0m∑m=0i-1‍f2m′′+6∑m=0i-1‍f1m∑m=0i-1‍f1m′′r21,i-1=-[l-4∑m=0i-1‍f2m∑m=0i-1‍f0m′′-4∑m=0i-1‍f2m′∑m=0i-1‍f0m′r21,i-1=-[l-3(∑m=0i-1‍f1m′)2-γ∑m=0i-1‍f2m′+∑m=0i-1‍θ2mr21,i-1=-[l-6∑m=0i-1‍f1m′′′∑m=0i-1‍f0m′-6∑m=0i-1‍f0m′′′∑m=0i-1‍f1m′r21,i-1=-[l+6∑m=0i-1‍f1miv∑m=0i-1‍f0m+6∑m=0i-1‍f0miv∑m=0i-1‍f1mr21,i-1=-[l+6∑m=0i-1‍f1m′′∑m=0i-1‍f0m′′],r22,i-1=-[∑m=0i-1‍θ2m′′+4Pr∑m=0i-1‍f0m∑m=0i-1‍θ2m′r22,i-1=-[l-6Pr∑m=0i-1‍f1m′∑m=0i-1‍θ1m′-6Pr∑m=0i-1‍f1m∑m=0i-1‍θ1m′r22,i-1=-[l+4Pr∑m=0i-1‍f2m∑m=0i-1‍f0m′-2Pr∑m=0i-1‍θ2m∑m=0i-1‍f0m′r22,i-1=-[l-3Pr∑m=0i-1‍θ1m∑m=0i-1‍f1m′r22,i-1=-[l+10PrEc∑m=0i-1‍f0m′∑m=0i-1‍f0m′′∑m=0i-1‍f1m′r22,i-1=-[l+2(∑m=0i-1‍f0m′′)2∑m=0i-1‍f1m′+4(∑m=0i-1‍f0m′′)2∑m=0i-1‍f2m′r22,i-1=-[l+4∑m=0i-1f0m∑m=0i-1f0m′′∑m=0i-1f1m′′′+4∑m=0i-1‍f0m∑m=0i-1‍f1m′′∑m=0i-1‍f0m′′′r22,i-1=-[l+4∑m=0i-1‍f0m′′∑m=0i-1‍f0m′′′∑m=0i-1‍f1m].
Equations (32)–(43) must be solved simultaneously subject to certain initial approximations f0 and θ0. We choose these initial approximations so that they satisfy the given boundary conditions. In this case suitable initial approximations are
(45)f0(Y)=1-e-Y-Ye-Y,θ0(Y)=e-Y.
We note that when fi and θi (i>1) have been found, the approximate solutions f(Y) and θ(Y) are obtained as
(46)f(Y)≈∑n=0M‍fn(Y),θ(Y)≈∑n=0M‍θn(Y),
where M is the order of the SLM approximation. Equations (32) and (43) can be solved by any numerical method. In this work the equations have been solved by the Chebyshev spectral collocation method. The method of solution is fully described in Awad et al. [3]. The system of differential equations is solved simultaneously using the MATLAB SLM code.


### 3.1. Results and Discussion

The problem that is investigated in this study is the steady laminar flow and natural convection from a cone in a viscoelastic fluid in the presence of viscous dissipation in a porous medium. The coupled nonlinear differential equations (7)–(10) were solved numerically using the successive linearisation method (SLM). In this section we discuss the effects of the viscoelastic parameter (Λ), porosity parameter (γ), Prandtl number (Pr), and Eckert numbers (Ec) on both the velocity and temperature profiles.In Table1 the comparison between our results for the local skin friction and Nusselt numbers and those of Ece [5] who used the Thomas algorithm shows that our method gives satisfactory results, thus confirming that the method is accurate.Table 1
Comparison of the values off′′(0) obtained by SLM against the Thomas algorithm of Ece [5] when Λ=0.
Pr
Ece [5]
Present
f
′
′
(
0
)
-
θ
′
(
0
)
f
′
′
(
0
)
-
θ
′
(
0
)
1
0.681483
0.638855
0.68148334
0.63885473
10
0.433268
1.275499
0.43327820
1.27552877To get a clear understanding of natural convection effects on the physics of the problem of a flow from a cone in a viscoelastic fluid with viscous dissipation, the investigation has been carried out for different viscoelastic numbersΛ, porosity parameter γ, the Eckert number Ec, and the Prandtl number Pr. The results for the skin friction and heat transfer coefficients are depicted in Tables 1 and 2.Table 2
Effect of the viscoelastic and porosity parameters and Eckert numberΛ,γ, and Ec on the local skin friction and heat transfer for Pr = 1.
Λ
γ
Ec
Pr
f
′
′
(
0
)
-
θ
′
(
0
)
−0.1
1
0.1
1
0.51437649
0.64214087
−0.05
1
0.1
1
0.53489736
0.59964040
−0.01
1
0.1
1
0.55491411
0.56204002
0
1
0.1
1
0.56041829
0.55213993
0.01
1
0.1
1
0.56612247
0.54203983
0.05
1
0.1
1
0.59093920
0.49963495
0.1
1
0.1
1
0.62646010
0.44213898
0.01
0
0.1
1
0.68990728
0.61800108
0.01
1
0.1
1
0.56612247
0.54203983
0.01
2
0.1
1
0.49213165
0.48868800
0.01
3
0.1
1
0.44155007
0.44852699
0.01
1
0.1
1
0.56612247
0.54203983
0.01
1
0.2
1
0.56678919
0.53544933
0.01
1
0.3
1
0.56745956
0.52882111
0.01
1
0.4
1
0.56813364
0.52215479
0.01
1
0.1
0.7
0.59466242
0.47714847
0.01
1
0.1
1
0.56612247
0.54203983
0.01
1
0.1
2
0.50981746
0.68640396
0.01
1
0.1
10
0.38403617
1.13367723In Table2 the effect of increasing the viscoelastic parameter increases the skin friction coefficient and the opposite effect is noted on the Nusselt number in the presence of the porous medium and viscous dissipation. Cortell [9] noted the same result. A faster increase is noted in the absence of the porous medium and the Eckert number. Increasing the porosity parameter reduces local skin friction and the same trend is noted on the Nusselt number. Skin friction increases with increasing Eckert number and the opposite trend is noted on the Nusselt number. The skin friction decreases with increasing Prandtl number, and the opposite trend is noted on the Nusselt number.Figures2–9 show the effects of various fluid properties on the velocity and temperature profiles.Figure 2
Velocity profiles for different values of the viscoelastic parameterΛ at Pr=1, Ec=0.1, and γ=1.Figure 3
Velocity profiles for different values of the Prandtl numberPr at Ec=0.1, γ=1, and Λ=0.1.Figure 4
Velocity profiles for different values of the porosity parameterγ at Pr=1, Ec=0.1, and Λ=0.1.Figure 5
Velocity profiles for different values of the Eckert numberEc at Pr=1, γ=1, and Λ=0.1.Figure 6
Temperature profiles for different values of the viscoelastic parameterΛ at Pr=1, γ=1, and Ec=0.1.Figure 7
Temperature profiles for different values of the Prandtl numberPr at Ec=0.1, γ=0.1, and Λ=0.1.Figure 8
Temperature profiles for different values of the porosity parameterγ at Pr=1, Λ=0.1, and Ec=0.1.Figure 9
Temperature profiles for different values of the Eckert numberEc at Pr=1, γ=1, and Λ=0.1.Figure2 shows that increasing the viscoelastic parameter increases the velocity across the boundary layer (see Butt et al. [24]).Increasing the Prandtl number decreases the velocity profile in the boundary layer as shown in Figure3; This is because when the Prandtl number is increased the conduction process is more enhanced than convection suggesting lower molecular motion causing fluid velocity to decrease.Figure4 shows the variation of the porosity parameter with velocity profile for the linear surface temperature. Increasing porosity parameter reduces the velocity profile across the boundary layer. The fluid particles move slower as the medium becomes less porous (see Singh and Agarwal [25]).Figure5 shows the variation of the Eckert number with velocity profile across the boundary layer. Increasing the Eckert number increases the velocity profile; this is caused by the increase in the kinetic energy caused by viscous dissipation in the boundary layer which leads to a small temperature gradient.Figure6 shows the effect of increasing the viscoelastic parameter on the temperature profiles. Increasing the viscoelastic parameter increases the temperature profile.Figure7 depicts the variation of the Prandtl number with temperature profiles. Increasing the Prandtl number decreases the temperature profile; The thermal diffusivity becomes smaller than the viscous diffusion rate causing smaller temperature profiles.Figure8 shows the variation of the porosity parameter with the temperature profile. Increasing the porosity parameter increases the temperature profile; when the fluid moves much slower due to the reduction in porosity heat transfer becomes more rapid.In Figure9 increasing the Eckert number increases the temperature profile; the heat produced due to viscous dissipation increases the temperature across the boundary layer.Figure10 shows the variation of the skin friction with the viscoelastic parameter at different values of the porosity parameter. Skin friction increases with increasing viscoelastic parameter and increasing the porosity parameter reduces skin friction.Figure 10
Skin frictionf′′(0) versus viscoelastic parameter Λ for different values of porosity parameter.Figure11 shows the variation of the Nusselt number with the viscoelastic parameter; increasing the viscoelastic parameter reduces Nusselt number and increasing the porosity parameter reduces the Nusselt number.Figure 11
Nusselt number-θ′(0) versus viscoelastic parameter Λ for different values of porosity parameter.Figure12 shows the effect of increasing the Eckert number on the skin friction and viscoelastic parameter. Increasing viscoelastic parameter increases skin friction and increasing the Eckert number increases the skin friction.Figure 12
Skin frictionf′′(0) versus viscoelastic parameter Λ for different values of Eckert numbers.In Figure13 the increase of viscoelastic parameter reduces the Nusselt number and increasing the Eckert number reduces the Nusselt number.Figure 13
Nusselt number-θ′(0) versus viscoelastic parameter Λ for different values of Eckert numbers.Figure14 shows that generally increasing the viscoelastic parameter increases the skin friction and increasing the Prandtl number reduces skin friction.Figure 14
Skin frictionf′′(0) versus viscoelastic parameter Λ for different values of the Prandtl numbers.In Figure15 increasing the viscoelastic parameter reduces the Nusselt number and increasing the Prandtl number increases the Nusselt number.Figure 15
Nusselt number-θ′(0) versus viscoelastic parameter Λ for different values of Prandtl numbers.


## 3.1. Results and Discussion

The problem that is investigated in this study is the steady laminar flow and natural convection from a cone in a viscoelastic fluid in the presence of viscous dissipation in a porous medium. The coupled nonlinear differential equations (7)–(10) were solved numerically using the successive linearisation method (SLM). In this section we discuss the effects of the viscoelastic parameter (Λ), porosity parameter (γ), Prandtl number (Pr), and Eckert numbers (Ec) on both the velocity and temperature profiles.In Table1 the comparison between our results for the local skin friction and Nusselt numbers and those of Ece [5] who used the Thomas algorithm shows that our method gives satisfactory results, thus confirming that the method is accurate.Table 1
Comparison of the values off′′(0) obtained by SLM against the Thomas algorithm of Ece [5] when Λ=0.
Pr
Ece [5]
Present
f
′
′
(
0
)
-
θ
′
(
0
)
f
′
′
(
0
)
-
θ
′
(
0
)
1
0.681483
0.638855
0.68148334
0.63885473
10
0.433268
1.275499
0.43327820
1.27552877To get a clear understanding of natural convection effects on the physics of the problem of a flow from a cone in a viscoelastic fluid with viscous dissipation, the investigation has been carried out for different viscoelastic numbersΛ, porosity parameter γ, the Eckert number Ec, and the Prandtl number Pr. The results for the skin friction and heat transfer coefficients are depicted in Tables 1 and 2.Table 2
Effect of the viscoelastic and porosity parameters and Eckert numberΛ,γ, and Ec on the local skin friction and heat transfer for Pr = 1.
Λ
γ
Ec
Pr
f
′
′
(
0
)
-
θ
′
(
0
)
−0.1
1
0.1
1
0.51437649
0.64214087
−0.05
1
0.1
1
0.53489736
0.59964040
−0.01
1
0.1
1
0.55491411
0.56204002
0
1
0.1
1
0.56041829
0.55213993
0.01
1
0.1
1
0.56612247
0.54203983
0.05
1
0.1
1
0.59093920
0.49963495
0.1
1
0.1
1
0.62646010
0.44213898
0.01
0
0.1
1
0.68990728
0.61800108
0.01
1
0.1
1
0.56612247
0.54203983
0.01
2
0.1
1
0.49213165
0.48868800
0.01
3
0.1
1
0.44155007
0.44852699
0.01
1
0.1
1
0.56612247
0.54203983
0.01
1
0.2
1
0.56678919
0.53544933
0.01
1
0.3
1
0.56745956
0.52882111
0.01
1
0.4
1
0.56813364
0.52215479
0.01
1
0.1
0.7
0.59466242
0.47714847
0.01
1
0.1
1
0.56612247
0.54203983
0.01
1
0.1
2
0.50981746
0.68640396
0.01
1
0.1
10
0.38403617
1.13367723In Table2 the effect of increasing the viscoelastic parameter increases the skin friction coefficient and the opposite effect is noted on the Nusselt number in the presence of the porous medium and viscous dissipation. Cortell [9] noted the same result. A faster increase is noted in the absence of the porous medium and the Eckert number. Increasing the porosity parameter reduces local skin friction and the same trend is noted on the Nusselt number. Skin friction increases with increasing Eckert number and the opposite trend is noted on the Nusselt number. The skin friction decreases with increasing Prandtl number, and the opposite trend is noted on the Nusselt number.Figures2–9 show the effects of various fluid properties on the velocity and temperature profiles.Figure 2
Velocity profiles for different values of the viscoelastic parameterΛ at Pr=1, Ec=0.1, and γ=1.Figure 3
Velocity profiles for different values of the Prandtl numberPr at Ec=0.1, γ=1, and Λ=0.1.Figure 4
Velocity profiles for different values of the porosity parameterγ at Pr=1, Ec=0.1, and Λ=0.1.Figure 5
Velocity profiles for different values of the Eckert numberEc at Pr=1, γ=1, and Λ=0.1.Figure 6
Temperature profiles for different values of the viscoelastic parameterΛ at Pr=1, γ=1, and Ec=0.1.Figure 7
Temperature profiles for different values of the Prandtl numberPr at Ec=0.1, γ=0.1, and Λ=0.1.Figure 8
Temperature profiles for different values of the porosity parameterγ at Pr=1, Λ=0.1, and Ec=0.1.Figure 9
Temperature profiles for different values of the Eckert numberEc at Pr=1, γ=1, and Λ=0.1.Figure2 shows that increasing the viscoelastic parameter increases the velocity across the boundary layer (see Butt et al. [24]).Increasing the Prandtl number decreases the velocity profile in the boundary layer as shown in Figure3; This is because when the Prandtl number is increased the conduction process is more enhanced than convection suggesting lower molecular motion causing fluid velocity to decrease.Figure4 shows the variation of the porosity parameter with velocity profile for the linear surface temperature. Increasing porosity parameter reduces the velocity profile across the boundary layer. The fluid particles move slower as the medium becomes less porous (see Singh and Agarwal [25]).Figure5 shows the variation of the Eckert number with velocity profile across the boundary layer. Increasing the Eckert number increases the velocity profile; this is caused by the increase in the kinetic energy caused by viscous dissipation in the boundary layer which leads to a small temperature gradient.Figure6 shows the effect of increasing the viscoelastic parameter on the temperature profiles. Increasing the viscoelastic parameter increases the temperature profile.Figure7 depicts the variation of the Prandtl number with temperature profiles. Increasing the Prandtl number decreases the temperature profile; The thermal diffusivity becomes smaller than the viscous diffusion rate causing smaller temperature profiles.Figure8 shows the variation of the porosity parameter with the temperature profile. Increasing the porosity parameter increases the temperature profile; when the fluid moves much slower due to the reduction in porosity heat transfer becomes more rapid.In Figure9 increasing the Eckert number increases the temperature profile; the heat produced due to viscous dissipation increases the temperature across the boundary layer.Figure10 shows the variation of the skin friction with the viscoelastic parameter at different values of the porosity parameter. Skin friction increases with increasing viscoelastic parameter and increasing the porosity parameter reduces skin friction.Figure 10
Skin frictionf′′(0) versus viscoelastic parameter Λ for different values of porosity parameter.Figure11 shows the variation of the Nusselt number with the viscoelastic parameter; increasing the viscoelastic parameter reduces Nusselt number and increasing the porosity parameter reduces the Nusselt number.Figure 11
Nusselt number-θ′(0) versus viscoelastic parameter Λ for different values of porosity parameter.Figure12 shows the effect of increasing the Eckert number on the skin friction and viscoelastic parameter. Increasing viscoelastic parameter increases skin friction and increasing the Eckert number increases the skin friction.Figure 12
Skin frictionf′′(0) versus viscoelastic parameter Λ for different values of Eckert numbers.In Figure13 the increase of viscoelastic parameter reduces the Nusselt number and increasing the Eckert number reduces the Nusselt number.Figure 13
Nusselt number-θ′(0) versus viscoelastic parameter Λ for different values of Eckert numbers.Figure14 shows that generally increasing the viscoelastic parameter increases the skin friction and increasing the Prandtl number reduces skin friction.Figure 14
Skin frictionf′′(0) versus viscoelastic parameter Λ for different values of the Prandtl numbers.In Figure15 increasing the viscoelastic parameter reduces the Nusselt number and increasing the Prandtl number increases the Nusselt number.Figure 15
Nusselt number-θ′(0) versus viscoelastic parameter Λ for different values of Prandtl numbers.


## 4. Conclusion

This study presented an analysis of flow and heat transfer in natural convection of viscoelastic fluid from a cone embedded in a porous medium with viscous dissipation. The nonlinear coupled governing equations were solved using the successive linearization method (SLM). The equations were first split into the zeroth, first, and second order of the viscoelastic parameter and solved together under the linear surface boundary conditions. The velocity and temperature profiles together with local skin friction and local Nusselt numbers were presented and investigated. It was found that increasing the viscoelastic parameter increased the skin friction, reduced the Nusselt number, and increased the velocity and temperature profiles. Increasing the porosity parameter decreased the skin friction and Nusselt number and decreased the velocity profile and the opposite effect was noted in the temperature profile. Increasing the Eckert number increased both velocity and temperature profiles and decreased the Nusselt number and the opposite was noted on the skin friction. The results compared well with those of Ece [5] in case when γ=Λ=Ec=0.



---

*Source: 934712-2013-12-24.xml*