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

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 linearizationmethod (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.


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

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 V 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.The governing equations in this buoyant-driven flow are given by 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 where  > 0 is a constant,  > 0 is the characteristic length, and the subscript ∞ refers to the ambient condition.
We introduce the nondimensional variables: where gives the following equations: where  =  sin , Λ = ( 0  0 /]) is the viscoelastic parameter known as the Deborah number, Gr is the Grashof number, Pr = ]/ is the Prandtl number, and Ec = ( 2 0 /  ) is the Eckert number.The corresponding boundary conditions are given as We now introduce the stream functions  = () and  = () defined by Substituting ( 6) and the similarity variables in (4) gives the following ordinary differential equations: With boundary conditions, 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]) where  is the coefficient of viscosity.The skin friction is defined as The skin friction coefficient can be expressed as The heat transfer rate at the surface of the cone is given by The Nusselt number can be expressed as Using the nondimensional variables ( 9)-( 10), the dimensionless wall heat rate is given by Nu Gr −1/4 = −  (0) .( 16)

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 The skin friction can be computed using 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: First order: Second order: The functions in the system ( 19)-( 30) may be expanded in series form as 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 1 + The coefficients  ,−1 ,  ,−1 ( = 0, 1, 2,  = 1, . . ., 5),  1,−1 , and  2,−1 are defined as Equations ( 32)-( 43) must be solved simultaneously subject to certain initial approximations  0 and  0 .We choose these initial approximations so that they satisfy the given boundary conditions.In this case suitable initial approximations are We note that when   and   ( > 1) have been found, the approximate solutions () and () are obtained as 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.

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.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.
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 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 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 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 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.
In Figure 13 the increase of viscoelastic parameter reduces the Nusselt number and increasing the Eckert number reduces the Nusselt number.
Figure 14 shows that generally increasing the viscoelastic parameter increases the skin friction and increasing the Prandtl number reduces skin friction.In Figure 15 increasing the viscoelastic parameter reduces the Nusselt number and increasing the Prandtl number increases the Nusselt number.

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

Figure 1 :
Figure 1: Physical model and coordinate system.

Figure 7
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.Figure8shows 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 Figure9increasing the Eckert number increases the temperature profile; the heat produced due to viscous

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.