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

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.

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 temperature 𝑇 (>𝑇_∞), and the ambient conditions are maintained at 𝑇_∞; the fluid has a constant viscosity 𝜈. The vertex angle of the cone is 2𝜙. The velocity components 𝑢 and 𝑣 are in the directions of 𝑥 and 𝑦, respectively, with the 𝑥-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 𝜕()+𝜕𝜕𝑥𝑟𝑢()𝑢𝜕𝑦𝑟𝑣=0,𝜕𝑢𝜕𝑥+𝑣𝜕𝑢𝜕𝜕𝑦=𝜈²𝑢𝜕𝑦²−𝜈𝐾𝑢−𝑘_𝑜𝑢𝜕³𝑢𝜕𝑥𝜕𝑦²𝜕+𝑣³𝑢𝜕𝑦³+𝜕𝑢𝜕𝜕𝑥²𝑢𝜕𝑦²−𝜕²𝑢𝜕𝑥𝜕𝑦𝜕𝑢𝜕𝑦+𝑔𝛽𝑇−𝑇_∞𝑢cos𝜙,𝜕𝑇𝜕𝑥+𝑣𝜕𝑇𝜕𝜕𝑦=𝛼²𝑇𝜕𝑦²+𝜈𝐶_𝑝𝜕𝑢𝜕𝑦²+𝑘₀𝜌𝐶_𝑝𝑢𝜕²𝑢𝜕𝑥𝜕𝑦𝜕𝑢𝜕𝜕𝑦+𝑣²𝑢𝜕𝑦²𝜕𝑢,𝜕𝑦 (1) where 𝑟=𝑥sin𝜙, 𝑔 is the acceleration due to gravity, 𝜈 is the kinematic viscosity for the fluid, 𝑘_𝑜 is the non-Newtonian parameter of the viscoelastic fluid, 𝛽 is the coefficient of thermal expansion, 𝛼 is the thermal diffusivity, 𝐶_𝑝 is the specific heat capacity for the fluid, 𝜌 is the density of the fluid, and 𝐾 is the permeability coefficient of the porous medium. The boundary conditions are given as 𝑢=𝑣=0,𝑇=𝑇_𝑤(𝑥)=𝑇_∞𝑥+𝐴𝐿at𝑦=0,𝜕𝑢𝜕𝑦,𝑢⟶0,𝑇⟶𝑇_∞,as𝑦⟶∞, (2) where 𝐴>0 is a constant, 𝐿>0 is the characteristic length, and the subscript ∞ refers to the ambient condition.

We introduce the nondimensional variables: 𝑥𝑋=𝐿,𝑌=Gr^1/4𝑦𝐿𝑟,𝑅=𝐿,𝑢𝑈=𝑈₀,𝑉=Gr^1/4𝑣𝑈₀,𝑇=𝑇−𝑇_∞𝑇_𝑤−𝑇_∞𝑈,Gr=⁰𝐿𝜈², (3) where 𝑈₀=[𝑔𝛽cos𝜙𝐿(𝑇_𝑤−𝑇_∞)]^1/2. Using (3) in (1) gives the following equations: 𝜕()+𝜕𝜕𝑋𝑅𝑈()𝑈𝜕𝑌𝑅𝑉=0,𝜕𝑈𝜕𝑋+𝑉𝜕𝑈=𝜕𝜕𝑌₂𝑈𝜕𝑌²−𝜈𝑈𝐾𝑈𝜕−Λ³𝑈𝜕𝑋𝜕𝑌²𝜕+𝑉³𝑈𝜕𝑌³+𝜕𝑈𝜕𝜕𝑋²𝑈𝜕𝑌²−𝜕²𝑈𝜕𝑋𝜕𝑌𝜕𝑈+𝜕𝑌𝑈𝜕𝑇,𝑇𝜕𝜕𝑋+𝑉𝑇=1𝜕𝑌𝜕Pr²𝑇𝜕𝑌²+Ec𝜕𝑈𝜕𝑌²𝑈𝜕+ΛEc²𝑈𝜕𝑋𝜕𝑌𝜕𝑈𝜕𝜕𝑌+𝑉²𝑈𝜕𝑌²𝜕𝑈,𝜕𝑌 (4) where 𝑅=𝑋sin𝜙, Λ=(𝑘₀𝑈₀/𝜈𝐿) is the viscoelastic parameter known as the Deborah number, Gr is the Grashof number, Pr=𝜈/𝛼 is the Prandtl number, and Ec=(𝑈²₀/𝐶_𝑝𝐴) is the Eckert number. The corresponding boundary conditions are given as 𝑈=𝑉=0,𝑇=𝑋at𝑌=0,𝜕𝑈𝜕𝑌,𝑈⟶0,𝑇⟶0as𝑌⟶∞. (5) We now introduce the following similarity variables 𝜓=𝑋𝑅𝑓(𝑌) and 𝑇=𝑋𝜃(𝑌) defined by 1𝑈=𝑅𝜕𝜓1𝜕𝑌,𝑉=−𝑅𝜕𝜓.𝜕𝑋 (6)

Substituting (6) and the similarity variables in (4) gives the following ordinary differential equations: 𝑓^′′′+2𝑓𝑓^′′−𝑓^′2+𝜃−𝛾𝑓^′−Λ2𝑓^′𝑓^′′′−2𝑓𝑓^𝑖𝑣−𝑓^′′2𝜃=0,_′′+Pr2𝑓𝜃^′−𝑓^′𝜃+PrEc𝑓^′′2+ΛPrEc𝑓^′𝑓^′′2−𝑓𝑓^′′𝑓^′′′=0. (7)(8) With boundary conditions, 𝑓(0)=𝑓^′(0)(0)𝑓=0,𝜃=1,_′(∞)⟶0,𝑓^′′(∞)⟶0,𝜃(∞)⟶0. (9)(10) 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]) 𝜏_𝑤=𝜇𝜕𝑢𝜕𝑦_𝑦=0+𝑘₀𝑢𝜕²𝑢𝜕𝑥𝜕𝑦−2𝜕𝑢𝜕𝑥𝜕𝑢𝜕𝑦_𝑦=0, (11) where 𝜇 is the coefficient of viscosity. The skin friction is defined as 𝑐_𝑓=𝜏_𝑤(1/2)𝜌𝑈²_∞,𝑐_𝑓=2𝑋Gr^1/4𝑓^′′(0)1+3Λ𝑓^′(0). (12) The skin friction coefficient can be expressed as 𝐶_𝑓Gr^1/42𝑋=𝑓^′′(0). (13) The heat transfer rate at the surface of the cone is given by 𝑞_𝑤𝑘=−𝑋𝜕𝑇𝜕𝑦_𝑦=0. (14) The Nusselt number can be expressed as Nu=𝐿𝑞_𝑤𝑘𝑇_𝑤−𝑇_∞. (15) Using the nondimensional variables (9)-(10), the dimensionless wall heat rate is given by NuGr^−1/4=−𝜃^′(0). (16)

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) The skin friction can be computed using 𝑓^′′(0)=𝑓^′′₀(0)+Λ𝑓^′′₁(0)+Λ²𝑓^′′₂(0)+⋯. (18)

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: 𝑓^′′′₀+2𝑓₀𝑓^′′₀−𝑓^′2₀+𝜃₀−𝛾𝑓^′₀𝜃=0,_′′+2Pr𝑓₀𝜃^′₀−Pr𝑓^′₀𝜃₀+PrEc𝑓^′′2₀𝑓=0,₀(0)=0,𝑓^′₀(0)=0,𝜃₀(0)𝑓=1,_′0(∞)=0,𝜃₀(∞)=0. (19)(20)(21)(22)

First order: 𝑓^′′′₁+2𝑓₀𝑓^′′₁+2𝑓₁𝑓^′′₀−2𝑓^′₀𝑓^′₁+𝜃₁−𝛾𝑓^′₁−2𝑓^′₀𝑓^′′′₀+2𝑓₀𝑓₀^𝑖𝑣+𝑓^′′2₀𝜃=0,_′′1+2Pr𝑓₀𝜃₁+2Pr𝑓₁𝜃^′₀−Pr𝑓^′₀𝜃₁−Pr𝑓^′₁𝜃₀+2PrEc𝑓^′′₀𝑓^′′₁+PrEc𝑓^′₀𝑓^′′2₀−PrEc𝑓₀𝑓^′′₀𝑓^′′′₀𝑓=0,₁(0)=0,𝑓^′₁(0)=0,𝑓^′₁(∞)𝜃=0,₁(0)=0,𝜃₁(∞)=0. (23)(24)(25)(26)

Second order: 𝑓^′′′₂+2𝑓₀𝑓^′′₂+2𝑓₁𝑓^′′₁+2𝑓^′₂𝑓^′′₀−2𝑓^′₀𝑓^′₂−𝑓^′2₁+𝜃₂−𝛾𝑓^′₂−2𝑓^′₀𝑓^′′′₁−2𝑓^′₁𝑓^′′′₀+2𝑓₀𝑓₁^𝑖𝑣+2𝑓₁𝑓₀^𝑖𝑣+2𝑓^′′₀𝑓^′′₁𝜃=0,_′′2+2Pr𝑓₀𝜃^′₂+2Pr𝑓₁𝜃^′₁+2Pr𝑓₂𝜃^′₀−Pr𝑓^′₀𝜃₂−Pr𝑓^′₁𝜃₁+PrEc2𝑓^′′₀𝑓^′′₂+𝑓^′′2₁+2𝑓^′₀𝑓^′′₀𝑓^′′₁+𝑓^′₁𝑓^′′2₀+𝑓^′₂𝑓^′′2₀+𝑓₀𝑓^′′₀𝑓^′′′₁+𝑓₀𝑓^′′₁𝑓^′′′₀+𝑓₁𝑓^′′₀𝑓^′′′₀𝑓=0,₂(0)=0,𝑓^′₂(0)=0,𝑓^′₂(∞)𝜃=0,₂(0)=0,𝜃₂(∞)=0. (27)(28)(29)(30) The functions in the system (19)–(30) may be expanded in series form as 𝑓₀𝑦=𝑓_0𝑖𝑦+^𝑖−1_𝑚=0𝑓_0𝑚𝑦,𝜃₀𝑦=𝜃_0𝑖𝑦+^𝑖−1_𝑚=0𝜃_0𝑚𝑦,𝑓₁𝑦=𝑓_1𝑖𝑦+^𝑖−1_𝑚=0𝑓_1𝑚𝑦,𝜃₁𝑦=𝜃_1𝑖𝑦+^𝑖−1_𝑚=0𝜃_1𝑚𝑦,𝑓₂𝑦=𝑓_2𝑖𝑦+^𝑖−1_𝑚=0𝑓_2𝑚𝑦,𝜃₂𝑦=𝜃_2𝑖𝑦+^𝑖−1_𝑚=0𝜃_2𝑚𝑦, (31) where 𝑓_0𝑖, 𝑓_1𝑖, and 𝑓_2𝑖 and 𝜃_0𝑖, 𝜃_1𝑖, and 𝜃_2𝑖 (𝑖=1,2,3,…) are unknown functions and 𝑓_0𝑚, 𝑓_1𝑚, and 𝑓_2𝑚 and 𝜃_0𝑚, 𝜃_1𝑚, and 𝜃_2𝑚 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 𝑓^′′′_0𝑖+𝑎_01,𝑖−1𝑓^′′_0𝑖+𝑎_02,𝑖−1𝑓^′_0𝑖+𝑎_03,𝑖−1𝑓_0𝑖+𝑎_04,𝑖−1𝜃_0𝑖=𝑟_01,𝑖−1,𝜃^′′_0𝑖+𝑏_01,𝑖−1𝜃^′_0𝑖+𝑏_02,𝑖−1𝜃_0𝑖+𝑏_03,𝑖−1𝑓^′′_0𝑖+𝑏_04,𝑖−1𝑓^′_0𝑖+𝑏_05,𝑖−1𝑓_0𝑖=𝑟_02,𝑖−1,𝑓_0𝑖(0)=0,𝑓^′_0𝑖(0)=0,𝜃_0𝑖(0)𝑓=0,_′0𝑖(∞)=0,𝜃_0𝑖(∞)=0,𝑓^′′_0𝑖(∞)𝑓=0,_{′′′1𝑖}+𝑎_11,𝑖−1𝑓^′′_1𝑖+𝑎_12,𝑖−1𝑓^′_1𝑖+𝑎_13,𝑖−1𝑓_1𝑖+𝑎_14,𝑖−1𝜃_1𝑖=𝑟_11,𝑖−1,𝜃^′′_1𝑖+𝑏_11,𝑖−1𝜃^′_1𝑖+𝑏_12,𝑖−1𝜃_1𝑖+𝑏_13,𝑖−1𝑓^′′_1𝑖+𝑏_14,𝑖−1𝑓^′_1𝑖+𝑏_15,𝑖−1𝑓_1𝑖=𝑟_12,𝑖−1,𝑓_1𝑖(0)=0,𝑓^′_1𝑖(0)=0,𝜃_1𝑖(0)𝑓=0,_′1𝑖(∞)=0,𝜃_1𝑖(∞)=0,𝑓^′′_1𝑖(∞)𝑓=0,_{′′′2𝑖}+𝑎_21,𝑖−1𝑓^′′_2𝑖+𝑎_22,𝑖−1𝑓^′_2𝑖+𝑎_23,𝑖−1𝑓_2𝑖+𝑎_24,𝑖−1𝜃_2𝑖=𝑟_21,𝑖−1,𝜃^′′_2𝑖+𝑏_21,𝑖−1𝜃^′_2𝑖+𝑏_22,𝑖−1𝜃_2𝑖+𝑏_23,𝑖−1𝑓^′′_2𝑖+𝑏_24,𝑖−1𝑓^′_2𝑖+𝑏_25,𝑖−1𝑓_2𝑖=𝑟_22,𝑖−1,𝑓_2𝑖(0)=0,𝑓^′_2𝑖(0)=0,𝜃_2𝑖(0)𝑓=0,_′2𝑖(∞)=0,𝜃_2𝑖(∞)=0,𝑓^′′_2𝑖(∞)=0. (32)(33)(34)(35)(36)(37)(38)(39)(40)(41)(42)(43) The coefficients 𝑎_{𝑗𝑘,𝑖−1}, 𝑏_{𝑗𝑘,𝑖−1} (𝑗=0,1,2,𝑘=1,…,5), 𝑟_{𝑗1,𝑖−1}, and 𝑟_{𝑗2,𝑖−1} are defined as 𝑎_01,𝑖−1=𝑎_11,𝑖−1=𝑎_21,𝑖−1=2^𝑖−1_𝑚=0𝑓_0𝑚,𝑎_02,𝑖−1=𝑎_12,𝑖−1=𝑎_22,𝑖−1=−2^𝑖−1_𝑚=0𝑓^′_0𝑚,𝑎+𝛾_03,𝑖−1=𝑎_13,𝑖−1=𝑎_23,𝑖−1=2^𝑖−1_𝑚=0𝑓^′′_0𝑚,𝑎_04,𝑖−1=𝑎_14,𝑖−1=𝑎_24,𝑖−1𝑏=𝐈,_01,𝑖−1=𝑏_11,𝑖−1=𝑏_21,𝑖−1=2Pr^𝑖−1_𝑚=0𝑓_0𝑚,𝑏_02,𝑖−1=𝑏_12,𝑖−1=𝑏_22,𝑖−1=−Pr^𝑖−1_𝑚=0𝑓^′_0𝑚,𝑏_03,𝑖−1=𝑏_13,𝑖−1=𝑏_23,𝑖−1=PrEc^𝑖−1_𝑚=0𝑓^′′_0𝑚,𝑏_04,𝑖−1=𝑏_14,𝑖−1=𝑏_24,𝑖−1=−Pr^𝑖−1_𝑚=0𝜃_0𝑚,𝑏_05,𝑖−1=𝑏_15,𝑖−1=𝑏_25,𝑖−1=2Pr^𝑖−1_𝑚=0𝜃^′_0𝑚,𝑟_01,𝑖−1=−^𝑖−1_𝑚=0𝑓^′′′_0𝑚+2^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝑓^′′_0𝑚−^𝑖−1_𝑚=0𝑓^′_0𝑚²−^𝑖−1_𝑚=0𝜃_0𝑚−𝛾^𝑖−1_𝑚=0𝑓^′_0𝑚,𝑟_02,𝑖−1=−^𝑖−1_𝑚=0𝜃^′′_0𝑚+2Pr^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝜃^′_0𝑚−Pr^𝑖−1_𝑚=0𝑓^′_0𝑚^𝑖−1_𝑚=0𝜃_0𝑚+PrEc^𝑖−1_𝑚=0𝑓^′′_0𝑚²,𝑟_11,𝑖−1=−^𝑖−1_𝑚=0𝑓^′′′_1𝑚+4^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝑓^′_1𝑚+4^𝑖−1_𝑚=0𝑓_1𝑚^𝑖−1_𝑚=0𝑓^′′_0𝑚−6^𝑖−1_𝑚=0𝑓^′′_0𝑚^𝑖−1_𝑚=0−𝛾^𝑖−1_𝑚=0𝑓^′_0𝑚+3^𝑖−1_𝑚=0𝑓^′′_0𝑚²+6^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝑓^𝑖𝑣_0𝑚−6^𝑖−1_𝑚=0𝑓^′_0𝑚^𝑖−1_𝑚=0𝑓^′′′_0𝑚+^𝑖−1_𝑚=0,𝑟_12,𝑖−1=−^𝑖−1_𝑚=0𝜃^′′_1𝑚+4Pr^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝜃^′_1𝑚−2Pr^𝑖−1_𝑚=0𝑓^′_0𝑚^𝑖−1_𝑚=0𝜃_1𝑚−2Pr^𝑖−1_𝑚=0𝑓^′_1𝑚^𝑖−1_𝑚=0𝜃_0𝑚+4Pr^𝑖−1_𝑚=0𝑓_1𝑚^𝑖−1_𝑚=0𝜃^′_0𝑚+4PrEc^𝑖−1_𝑚=0𝑓^′′_0𝑚^𝑖−1_𝑚=0𝑓^′′_1𝑚+4PrEc^𝑖−1_𝑚=0𝑓^′_0𝑚^𝑖−1_𝑚=0𝑓_0𝑚²−4PrEc^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝑓^′′_0𝑚^𝑖−1_𝑚=0𝑓^′′′_0𝑚,𝑟_21,𝑖−1=−^𝑖−1_𝑚=0𝑓^′′′₂+4^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝑓^′′_2𝑚+6^𝑖−1_𝑚=0𝑓_1𝑚^𝑖−1_𝑚=0𝑓^′′_1𝑚−4^𝑖−1_𝑚=0𝑓_2𝑚^𝑖−1_𝑚=0𝑓^′′_0𝑚−4^𝑖−1_𝑚=0𝑓^′_2𝑚^𝑖−1_𝑚=0𝑓^′_0𝑚−3^𝑖−1_𝑚=0𝑓^′_1𝑚²−𝛾^𝑖−1_𝑚=0𝑓^′_2𝑚+^𝑖−1_𝑚=0𝜃_2𝑚−6^𝑖−1_𝑚=0𝑓^′′′_1𝑚^𝑖−1_𝑚=0𝑓^′_0𝑚−6^𝑖−1_𝑚=0𝑓^′′′_0𝑚^𝑖−1_𝑚=0𝑓^′_1𝑚+6^𝑖−1_𝑚=0𝑓^𝑖𝑣_1𝑚^𝑖−1_𝑚=0𝑓_0𝑚+6^𝑖−1_𝑚=0𝑓^𝑖𝑣_0𝑚^𝑖−1_𝑚=0𝑓_1𝑚+6^𝑖−1_𝑚=0𝑓^′′_1𝑚^𝑖−1_𝑚=0𝑓^′′_0𝑚,𝑟_22,𝑖−1=−^𝑖−1_𝑚=0𝜃^′′_2𝑚+4Pr^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝜃^′_2𝑚−6Pr^𝑖−1_𝑚=0𝑓^′_1𝑚^𝑖−1_𝑚=0𝜃^′_1𝑚−6Pr^𝑖−1_𝑚=0𝑓_1𝑚^𝑖−1_𝑚=0𝜃^′_1𝑚+4Pr^𝑖−1_𝑚=0𝑓_2𝑚^𝑖−1_𝑚=0𝑓^′_0𝑚−2Pr^𝑖−1_𝑚=0𝜃_2𝑚^𝑖−1_𝑚=0𝑓^′_0𝑚−3Pr^𝑖−1_𝑚=0𝜃_1𝑚^𝑖−1_𝑚=0𝑓^′_1𝑚+10PrEc^𝑖−1_𝑚=0𝑓^′_0𝑚^𝑖−1_𝑚=0𝑓^′′_0𝑚^𝑖−1_𝑚=0𝑓^′_1𝑚+2^𝑖−1_𝑚=0𝑓^′′_0𝑚^2𝑖−1_𝑚=0𝑓^′_1𝑚+4^𝑖−1_𝑚=0𝑓^′′_0𝑚^2𝑖−1_𝑚=0𝑓^′_2𝑚+4^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝑓^′′_0𝑚^𝑖−1_𝑚=0𝑓^′′′_1𝑚+4^𝑖−1_𝑚=0𝑓_0𝑚^𝑖−1_𝑚=0𝑓^′′_1𝑚^𝑖−1_𝑚=0𝑓^′′′_0𝑚+4^𝑖−1_𝑚=0𝑓^′′_0𝑚^𝑖−1_𝑚=0𝑓^′′′_0𝑚^𝑖−1_𝑚=0𝑓_1𝑚. (44) Equations (32)–(43) must be solved simultaneously subject to certain initial approximations 𝑓₀ and 𝜃₀. We choose these initial approximations so that they satisfy the given boundary conditions. In this case suitable initial approximations are 𝑓₀(𝑌)=1−𝑒^−𝑌−𝑌𝑒^−𝑌,𝜃₀(𝑌)=𝑒^−𝑌. (45) We note that when 𝑓_𝑖 and 𝜃_𝑖 (𝑖>1) have been found, the approximate solutions 𝑓(𝑌) and 𝜃(𝑌) are obtained as 𝑓(𝑌)≈^𝑀_𝑛=0𝑓_𝑛(𝑌),𝜃(𝑌)≈^𝑀_𝑛=0𝜃_𝑛(𝑌), (46) where 𝑀 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 Table 1 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 of 𝑓^′′(0) obtained by SLM against the Thomas algorithm of Ece [5] when Λ=0.

To 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.

In Table 2 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.

Figures 2–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 number Pr 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 number Ec 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 number Pr 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 number Ec at Pr=1, 𝛾=1, and Λ=0.1.

Figure 2 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 Figure 3; 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.

Figure 4 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]).

Figure 5 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.

Figure 6 shows the effect of increasing the viscoelastic parameter on the temperature profiles. Increasing the viscoelastic parameter increases the temperature profile.

Figure 7 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.

Figure 8 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 Figure 9 increasing the Eckert number increases the temperature profile; the heat produced due to viscous dissipation increases the temperature across the boundary layer.

Figure 10 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 friction 𝑓′′(0) versus viscoelastic parameter Λ for different values of porosity parameter.

Figure 11 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.

Figure 12 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 friction 𝑓′′(0) versus viscoelastic parameter Λ for different values of Eckert numbers.

In Figure 13 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.

Figure 14 shows that generally increasing the viscoelastic parameter increases the skin friction and increasing the Prandtl number reduces skin friction.

Figure 14 

Skin friction 𝑓′′(0) versus viscoelastic parameter Λ for different values of the Prandtl numbers.

In Figure 15 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.