A study is performed on two-dimensional flow and
diffusion of chemically reactive species of Casson fluid from an unsteady stretching surface in porous medium in the presence of a magnetic field. The boundary layer velocity, temperature, and concentration profiles are numerically computed for different governing parameters. The paper intends to show unique results of a combination of heat transfer and chemical reaction in Casson fluid flow. The resulting partial differential equations are converted to a system of ordinary differential equations using the appropriate similarity transformation, which are solved by using the Runge-Kutta-Fehlberg numerical scheme. The results in this work are validated by the comparison with other authors.
1. Introduction
The study of Casson fluid has attracted attention to many researchers due to its application in the field of metallurgy, food processing, drilling operations, and bioengineering operations. Its application extends to the manufacturing of pharmaceutical products, coal in water, china clay, paints, synthetic lubricants, and biological fluids such as synovial fluids, sewage sludge, jelly, tomato sauce, honey, soup, and blood due to its contents such as plasma, fibrinogen, and protein, making the study of Casson fluid important in fluid dynamics. Casson fluid is classified as a non-Newtonian fluid due to its rheological characteristics. These characteristics show shear stress-strain relationships that are significantly different from Newtonian fluids. The study of non-Newtonian fluids has not been thoroughly covered due to the complex representation of their constitutive equations. It is therefore important to undertake this study of Casson fluid. Most studies have concentrated on viscoelastic fluids in which different constitutive equations have been suggested. This work can be applied to chemical processing equipment in which some fluids react chemically with some ingredients present in them.
The driving force for mass transfer is a combination of temperature and concentration gradients. In this study the effect of chemical reaction on the fluid is considered as in Mukhopadhyay and Vajravelu [1]. The study of boundary layer flow over a stretching sheet has been studied by Mukhopadhyay et al. [2], among others, who investigated Casson fluid flow over an unsteady stretching surface. In their work they did not consider mass transfer and they considered a different wall temperature expression. Abd El-Aziz [3] studied mixed convection flow of a micropolar fluid from an unsteady stretching surface with viscous dissipation. In this work he considered a similar stretching velocity, wall temperature, and wall concentration distribution. We extended the work of Grubka and Bobba [4] who investigated heat transfer characteristics of a continuous stretching surface with variable temperature in which we introduced the MHD and porous medium source terms and chemical reaction effects. Sharidan et al. [5] studied similarity solutions for the unsteady boundary layer flow and heat transfer due to a stretching sheet; Nadeem et al. [6, 7] and Ahmed and Nazar [8] also studied Casson fluid over a stretching sheet and in their paper they assumed that the velocity of the stretching surface is linearly proportional to the distance from fixed origin.
In the study of non-Newtonian fluids many authors have studied the flow of blood as Casson fluid. The studies were carried out by Rohlf and Tenti [9], among others, who investigated the role of Womersley number in pulsatile blood flow, a theoretical study of the Casson model; Sankar and Lee [10, 11] investigated two-fluid nonlinear model for flow in catheterized blood vessels and two-fluid Casson model for pulsatile blood flow through stenosed arteries, respectively. Shaw et al. [12] studied pulsatile Casson fluid flow through stenosed bifurcated artery. In relation to blood flow there are other research works that were done in different geometries such as flows in microslit channels, slightly curved channels, and peristaltic transport as in [13–15].
The study of Casson fluid in porous media was also studied by Nadeem et al. [16] who considered MHD three-dimensional Casson fluid flow past a porous linearly stretching sheet. Dash et al. [17] studied Casson fluid flow in a pipe filled with homogeneous porous medium. Tripathi [18] investigated the transient peristaltic heat flow through a finite porous channel. Pramanik [19] studied Casson fluid flow and heat transfer past an exponentially porous stretching surface in the presence of thermal radiation. Ramachandra et al. [20] investigated flow and heat transfer of Casson fluid from a horizontal circular cylinder with partial slip in a non-Darcy porous medium. In their work they considered slip conditions on the wall.
There are many studies that investigated fluid flow with chemical reactions. Kameswaran et al. [21] investigated homogeneous-heterogeneous reactions in a nanofluid flow due to a porous stretching sheet, Shaw et al. [22] studied homogeneous-heterogeneous reactions in a nanofluid flow due to a porous stretching sheet, and Chamkha et al. [23] investigated similarity solutions for unsteady heat and mass transfer from a stretching surface embedded in a porous medium with suction/injection and chemical reaction effects.
Although there are many applications and use of non-Newtonian fluids in industry, the study of Casson fluid has not been thoroughly investigated for heat and mass transfer past a stretching surface. In this work we extend the work of Mukhopadhyay and Vajravelu [1] and Grubka and Bobba [4] in which the energy equation, the source terms for porous medium, and magnetic field are introduced. Similarity transformations are used to convert the partial differential equations into ordinary differential equations which are then solved by using Runge-Kutta-Fehlberg integration scheme and the successive linearization method described by Makanda et al. [24]. In this work we investigate the effect of varying unsteadiness parameter, Casson, Schmidt, and Prandtl numbers, and the reaction rate parameter on the velocity, temperature, and concentration profiles with the depiction of graphical illustrations.
2. Mathematical Formulation
Consider two-dimensional laminar boundary layer flow, temperature, and mass transfer of an incompressible Casson fluid flow over an unsteady stretching sheet. The flow of heat and mass transfer starts at t=0. The sheet is pulled out of the slit at the origin (x=0,y=0) and moves with velocity U(x,t)=ax/(1-αt), a>0, α≥0 are constants, and a is the initial stretching rate. The rheological equation of state for an isotropic and incompressible flow of a Casson fluid is given as in [1, 2, 20] by(1)τij=2μB+Py2πeij,π>πc2μB+Py2πceij,π<πc,where π=eijeij and eij is the (i,j)th component of the deformation rate, π is the product of the deformation rate with itself, πc is a critical value of this product based on the non-Newtonian model, μB is the plastic dynamic viscosity of the non-Newtonian fluid, and Py is the yield stress of the fluid. Given that Tw and Cw are, respectively, the temperature and concentration at the sheet and T∞ and C∞ are, respectively, the ambient conditions, the positive x coordinate is measured along the stretching sheet and the positive y coordinate is measured perpendicular to the sheet. It is assumed that both temperature and concentration at the surface vary with distance from the origin and time. The temperature Tw and concentration Cw at the surface are therefore given by(2)Twx,t=T∞+bx1-αt2,Cwx,t=C∞+cx1-αt2,where b and c are constants. The surface temperature and surface concentration increase if b and c are positive and reduce if they are negative from T∞ and C∞ at the origin to x and the temperature and concentration increase/decrease along the sheet. It is assumed that radiation effects and viscous dissipation are negligible. The expressions Uw(x,t), Tw(x,t), and Cw(x,t) are only valid for t<α-1 but not when α=0. Under these assumptions the governing equations in this flow are given as(3)∂∂xu+∂∂yv=0,(4)∂u∂t+u∂u∂x+v∂u∂y=ν1+1β∂2u∂y2-νKu-σB02ρu,∂T∂t+u∂T∂x+v∂T∂y=α0∂2T∂y2,∂C∂t+u∂C∂x+v∂C∂y=D∂2C∂y2-kC-C∞,where ν is kinematic viscosity of Casson fluid, β=μB2πc/Py is the non-Newtonian Casson parameter, K is the permeability of the porous medium, σ is the electrical conductivity, B0 is the strength of the magnetic field, ρ is the density of the Casson fluid, D is the diffusion coefficient of species in the fluid, α0 is the thermal diffusivity, and k(t)=k0/(1-αt) is the time dependent reaction rate, where k>0 represents destructive reaction, k<0 represents constructive reaction, and k0 is a constant. The boundary conditions are given as(5)u=Ux,t,v=0,T=Twx,t,C=Cwx,t,y=0,u⟶0,T⟶T∞,C⟶C∞,asy⟶∞,where the subscript ∞ refers to the ambient condition.
We introduce the nondimensional variables(6)u=∂ψ∂y,v=-∂ψ∂x,η=aν1-αty,ψ=νa1-αtxfη,Twx,t=T∞+bx1-αt2θη,Cwx,t=C∞+cx1-αt2ϕη,where ψ(x,y,t) is the stream function which satisfies the continuity equation (3). The velocity components are defined as(7)u=∂ψ∂y=Uwf′η,v=-∂ψ∂x=-νa1-αt.The governing equations reduce to(8)1+1βf′′′+ff′′-f′2-A22f′+ηf′′-Λ+M2f′=0,1Prθ′′+fθ′-f′θ-A24θ+ηθ′=0,1Scϕ′′+fϕ′-f′ϕ-A24ϕ+ηϕ′-Rϕ=0,with boundary conditions(9)f0=0,f′0=1,θ0=1,ϕ0=1,f′∞⟶0,θ∞⟶0,ϕ∞⟶0,where A=α/a is the unsteadiness parameter, Pr=ν/α0 is the Prandtl number, Sc=ν/D is the Schmidt number, and R=k0/a is the reaction parameter. Λ=ν(1-αt)/aK is the permeability parameter; M2=σB02(1-αt)/ρa is the magnetic parameter. The nondimensional temperature and concentration are, respectively, given by θ=(T-T∞)/(Tw-T∞) and ϕ=(C-C∞)/(Cw-C∞).
The parameters of engineering interests are the local skin friction and the Nusselt and Sherwood numbers which are defined as(10)Cfx=21+1βRex-1/2f′′0.The local Nusselt and Sherwood numbers are defined as(11)Nux=xα0qwTw-T∞,Shx=xDJwCw-C∞,(12)qw=-α0∂T∂yy=0,Jw=-D∂C∂yy=0.Using expressions (11) and (6), (13)Nux=-Rex1/2θ′0,Shx=-Rex1/2ϕ′0,where Rex is the Reynolds number defined as Rex=Uwx/ν. It is important at this stage to mention that the system of equations (8)-(9) reduce to those of Grubka and Bobba [4] when (1/β→0), A=λ=M=Sc=0. The present problem reduces to that of Grubka and Bobba [4], A=0 denote steady flow, and in their paper they obtained an exact solution in terms of Kummer’s functions written in terms of the confluent hypergeometric functions.
The solution of the boundary value problem (8)-(9) was solved using the Runge-Kutta-Fehlberg integration scheme. In the method we choose finite values of η→∞. This value is the boundary layer thickness given by η∞. We begin by choosing an initial guess of η∞ to obtain the values f′′0, -θ′0, and -ϕ′0. The solution is repeated with new values until two consecutive values differ by 10-6.
3. Results and Discussion
To obtain a clear understanding of the flow of Casson fluid, we discuss the physics of the problem by studying the effects of the unsteadiness (A), permeability (Λ), magnetic (M), Prandtl (Pr), Schmidt (Sc), and reaction rate (R) numbers on velocity, temperature, and concentration profiles. We also study the variation of skin friction, the local Nusselt number, and the Sherwood number with unsteadiness parameter. For validation of the numerical method used in this study, results for the Nusselt number -θ(0) for the Newtonian fluid were compared to those of Abd El-Aziz [3] and Grubka and Bobba [4] for the unsteadiness parameter A=Λ=M=Sc=0. The comparison is shown in Table 1 and it is found to be in agreement with at least four decimal places. To further verify the accuracy of the numerical scheme used, the successive linearization method (SLM) was used and there was agreement with Runge-Kutta-Fehlberg integration scheme.
Comparison of the values of -θ(0) for Λ=M=Sc=0 and various values of A and Pr with previously published data.
A
Pr
Grubka and Bobba [4]
Abd El-Aziz [3]
SLM
Present results
0
0.72
0.8086
0.80873135
0.80873007
0.80863761
1
1.0000
1.00000000
1.00000000
1.00000006
3
1.9237
1.92368255
1.92367361
1.92367736
10
3.7207
3.72067395
3.72066225
3.72066701
To get a clear understanding of the behavior of velocity, temperature, and concentration profiles of Casson fluid, a detailed numerical calculation is done for different parameter values that describe the nature of flow and the results are depicted through Figures 1–5.
Effects of different parameters on velocity profiles.
Velocity profiles for different values of the Casson parameter β and unsteadiness parameter A at Λ=0.5, M=1, Pr=1, Sc=1, and R=0.5
Velocity profiles for different values of the permeability parameter Λ and unsteadiness parameter A at β=2, M=1, Pr=1, Sc=1, and R=0.5
Velocity profiles for different values of the Casson parameter β and magnetic parameter M at A=0.5, Λ=0.5, Pr=1, Sc=1, and R=0.5
Velocity profiles for different values of the permeability parameter Λ and magnetic parameter M at β=2, A=0.5, Pr=1, Sc=1, and R=0.5
Effects of different parameters on temperature profiles.
Temperature profiles for different values of the Casson parameter β and unsteadiness parameter A at Λ=0.5, M=1, Pr=1, Sc=1, and R=0.5
Temperature profiles for different values of the permeability parameter Λ and unsteadiness parameter A at β=2, M=1, Pr=1, Sc=1, and R=0.5
Temperature profiles for different values of the Prandtl numbers Pr and unsteadiness parameter A at β=2, Λ=0.5, M=1, Sc=1, and R=0.5
Temperature profiles for different values of the Casson parameter β and magnetic parameter M at A=0.5, Λ=0.5, Pr=1, Sc=1, and R=0.5
Effects of different parameters on temperature profiles.
Temperature profiles for different values of the permeability parameter Λ and magnetic parameter M at β=2, A=0.5, Pr=1, Sc=1, and R=0.5
Temperature profiles for different values of the Prandtl number Pr and magnetic parameter M at β=2, A=0.5, Λ=0.5, Sc=1, and R=0.5
Effects of different parameters on concentration profiles.
Concentration profiles for different values of the Casson parameter β and unsteadiness parameter A at Λ=0.5, M=1, Pr=1, Sc=1, and R=0.5
Concentration profiles for different values of the permeability parameter Λ and unsteadiness parameter A at β=2, M=1, Pr=1, Sc=1, and R=0.5
Concentration profiles for different values of the Casson parameter β and magnetic parameter M at A=0.5, Λ=0.5, Pr=1, Sc=1, and R=0.5
Concentration profiles for different values of the permeability parameter Λ and magnetic parameter M at A=0.5, β=2, Pr=1, Sc=1, and R=0.5
Effects of different parameters on concentration profiles.
Concentration profiles for different values of the Schmidt number Sc and unsteadiness parameter A at β=2, Λ=0.5, Pr=1, M=1, and R=0.5
Concentration profiles for different values of the Schmidt number Sc and magnetic parameter M at A=0.5, β=2, Λ=0.5, Pr=1, and R=0.5
Concentration profiles for different values of the Casson parameter β and reaction rate R at A=0.5, M=1, Λ=0.5, Pr=1, and Sc=1
Concentration profiles for different values of the permeability parameter Λ and reaction rate R at A=0.5, β=2, M=1, Pr=1, and Sc=1
Figure 1 shows that increasing the Casson parameter β and the permeability parameter Λ decreases velocity profiles. Both situations retard fluid motion across the boundary layer as shown in Figures 1(a) to 1(d). The velocity profiles for unsteady flow become similar to steady flow far away from the boundary at η>2. This indicates that for A=0.5 the flow away from the boundary is steady. Increasing A retards velocity profiles across the boundary layer as shown in Figures 1(a) and 1(b). In Figure 1(b) it is noted that increasing A reduces velocity profiles and the reverse effect is noted around η=2.5 indicating slightly steady flow far away from the boundary. A slightly different observation is noted if the magnetic parameter is increased from M=1 to M=3; the velocity profiles decrease significantly reducing the momentum boundary layer thickness as shown in Figures 1(c) and 1(d).
Figure 2 shows the variation of different parameters with temperature profiles. Increasing the Casson parameter increases temperature profiles; this decreases the yield stress. The fluid behaves as Newtonian as β increases; this retards fluid motion. The effect of increasing β leads to enhancing the temperature profiles for both steady and unsteady flows as shown in Figure 2(a). The same observation is noted in Figure 2(d). Increasing the permeability parameter Λ enhances temperature profiles. High values of Λ lead to the reduction of fluid velocity enhancing temperature profiles. This effect is more pronounced in steady flow and the thickening of the thermal boundary layer increases as Λ increases. In Figure 2(c) increasing the Prandtl number Pr reduces temperature profiles; Pr is inversely proportional to the thermal diffusivity which is low for Casson fluid. In Figures 2(d), 3(a), and 3(b), increasing M increases temperature profiles; the slowing down of fluid flow leads to the buildup of heat. It is more pronounced in low values of M. Increasing A leads to the reduction of temperature profiles; less heat is transferred from the sheet to the fluid; hence, temperature decreases as depicted in Figures 2(a), 2(b), and 2(c).
In Figure 4, increasing the Casson parameter β results in the increase of the solute concentration; as the fluid becomes Newtonian, solute transfer is enhanced as shown in Figures 4(a), 4(d), and 5(c). This is more pronounced in steady flow. The opposite effect is noted in increasing A which reduces solute concentration. Increasing Λ increases concentration profiles; since this reduces fluid motion, it leads to the buildup of solute concentration as shown in Figures 4(b), 4(d), and 5(d). Increasing Sc decreases concentration profiles; the diffusion rate is smaller in Casson fluid and hence reducing solute concentration as depicted in Figures 5(a) and 5(b). In Figures 4(c), 4(d), and 5(b), increasing M increases profiles; the slowing down of fluid flow leads to the buildup of mass. It is more pronounced in low values of M.
Figures 5(c) and 5(d) show the effects of reaction rate R on the concentration profiles in Casson fluid in porous medium. Increasing the absolute value of R decreases the concentration profiles. The reaction rate is in general a destructive agent which leads to the reduction of the solute boundary layer. It shows the increase of solute boundary layer in the case of constructive reaction R<0 and the decrease of solute boundary layer in the case of destructive reaction R>0. The effect of increasing β and Λ is more pronounced in unsteady flow.
Figures 6(a) and 6(b) show effects of Casson parameter β and permeability parameter Λ on velocity gradient at the wall with unsteadiness parameter. The magnitude of f′′(0) related to skin friction decreases with increasing unsteadiness parameter A and also with Casson parameter β and permeability parameter Λ. The magnitude of the heat transfer rate at the surface -θ′(0) decreases with increasing β and Λ and increases with A as shown in Figure 7. The mass transfer coefficient at the surface -ϕ′(0) decreases with increasing β and Λ and increases with A as shown in Figure 8.
Variation of (a) skin friction coefficient with unsteadiness parameter for different β values and (b) with different Λ values.
Skin friction versus unsteadiness parameter at different values of the Casson parameter β at M=1, Λ=0.5, Pr=1, Sc=1, and R=0.5
Skin friction versus unsteadiness parameter at different values of the permeability parameter Λ at M=1, β=2, Pr=1, Sc=1, and R=0.5
Variation of (a) heat transfer coefficient with unsteadiness parameter for different β values and (b) with different Λ values.
Heat transfer coefficient versus unsteadiness parameter at different values of the Casson parameter β at M=1, Λ=0.5, Pr=1, Sc=1, and R=0.5
Heat transfer coefficient versus unsteadiness parameter at different values of the permeability parameter Λ at M=1, β=0.5, Pr=1, Sc=1, and R=0.5
Variation of (a) mass transfer coefficient with unsteadiness parameter for different β values and (b) with different Λ values.
Mass transfer coefficient versus unsteadiness parameter at different values of the Casson parameter β at M=1, Λ=0.5, Pr=1, Sc=1, and R=0.5
Heat transfer coefficient versus unsteadiness parameter at different values of the permeability parameter Λ at M=1, β=2, Pr=1, Sc=1, and R=0.5
4. Conclusion
The study presented in this analysis of diffusion of chemically reactive species in Casson fluid flow over an unsteady stretching surface in a porous medium in the presence of a magnetic field provides numerical solutions for the boundary layer flow and heat and mass transfer. The coupled nonlinear governing equations were solved using the Runge-Kutta-Fehlberg integration scheme. Increasing the unsteadiness parameter decreases velocity profiles. Increasing the Casson parameter decreases the velocity profiles. The presence of heat transfer does not seem to significantly change the results obtained by Mukhopadhyay and Vajravelu [1]; the main difference is noted in the presence of magnetic and permeability parameters which they did not consider. In this study it is generally noted that increasing the magnetic parameter and permeability parameter decreases the velocity profiles but increases the skin friction; it then decreases the coefficient of heat transfer and the concentration profiles.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
MukhopadhyayS.VajraveluK.Diffusion of chemically reactive species in Casson fluid flow over an unsteady permeable stretching surface201325459159810.1016/S1001-6058(11)60400-X2-s2.0-84884168675MukhopadhyayS.deP. R.BhattacharyyaK.LayekG. C.Casson fluid flow over an unsteady stretching surface20134493393810.1016/j.asej.2013.04.0042-s2.0-84887815558Abd El-AzizM.Mixed convection flow of a micropolar fluid from an unsteady stretching surface with viscous dissipation201321338539410.1016/j.joems.2013.02.010MR3116467ZBL06244814GrubkaL. J.BobbaK. M.Heat transfer characteristics of a continuous stretching surface with variable temperature1985107124825010.1115/1.32473872-s2.0-0022013926SharidanS.MahmoodT.PopI.Similarity solutions for the unsteady boundary layer flow and heat transfer due to a stretching sheet200611647654NadeemS.Ul HaqR.LeeC.MHD flow of a Casson fluid over an exponentially shrinking sheet20121961550155310.1016/j.scient.2012.10.0212-s2.0-84871963705NadeemS.HaqR. U.KhanZ. H.Numerical study of MHD boundary layer flow of a Maxwell fluid past a stretching sheet in the presence of nanoparticles201445112112610.1016/j.jtice.2013.04.0062-s2.0-84894247694AhmedK.NazarR.Magnetohydrodynamic three dimensional flow and heat transfer over a stretching surface in a viscoelastic fluid2011313346RohlfK.TentiG.The role of the Womersley number in pulsatile blood flow a theoretical study of the Casson model200134114114810.1016/S0021-9290(00)00103-22-s2.0-0035217964SankarD. S.LeeU.Two-fluid non-linear model for flow in catheterized blood vessels200843762263110.1016/j.ijnonlinmec.2008.02.0072-s2.0-46849104712SankarD. S.LeeU.Two-fluid Casson model for pulsatile blood flow through stenosed arteries: a theoretical model20101582086209710.1016/j.cnsns.2009.08.021MR25926222-s2.0-76049124163ShawS.GorlaR. S. R.MurthyP. V. S. N.NgC.-O.Pulsatile casson fluid flow through a stenosed bifurcated artery2009361436310.1615/InterJFluidMechRes.v36.i1.302-s2.0-64749102780MernoneA. V.MazumdarJ. N.LucasS. K.A mathematical study of peristaltic transport of a Casson fluid2002357-889591210.1016/S0895-7177(02)00058-4MR1901295ZBL1022.760582-s2.0-0037177273DasB.BatraR. L.Secondary flow of a Casson fluid in a slightly curved tube199328556757710.1016/0020-7462(93)90048-P2-s2.0-0027659584NgC.-O.Combined pressure-driven and electroosmotic flow of Casson fluid through a slit microchannel20131981910.1016/j.jnnfm.2013.03.0032-s2.0-84878122440NadeemS.Ul HaqR.AkbarN. S.KhanZ. H.MHD three-dimensional Casson fluid flow past a porous linearly stretching sheet201352457758210.1016/j.aej.2013.08.0052-s2.0-84889684506DashR. K.MehtaK. N.JayaramanG.Casson fluid flow in a pipe filled with a homogeneous porous medium199634101145115610.1016/0020-7225(96)00012-2ZBL0900.760222-s2.0-0001195756TripathiD.Study of transient peristaltic heat flow through a finite porous channel2013575-61270128310.1016/j.mcm.2012.10.030MR30341312-s2.0-84873088108PramanikS.Casson fluid flow and heat transfer past an exponentially porous stretching surface in presence of thermal radiation20145120521210.1016/j.asej.2013.05.0032-s2.0-84897603369RamachandraP. V.SubbaR. A.AnwaB. O.Flow and heat transfer of casson fluid from a horizontal cylinder with partial slip in non-darcy porous medium20132, article 2KameswaranP. K.ShawS.SibandaP.MurthyP. V. S. N.Homogeneous-heterogeneous reactions in a nanofluid flow due to a porous stretching sheet201357246547210.1016/j.ijheatmasstransfer.2012.10.0472-s2.0-84869139569ShawS.KameswaranP. K.SibandaP.Homogeneous-heterogeneous reactions in micropolar fluid flow from a permeable stretching or shrinking sheet in a porous medium20132013, article 7710.1186/1687-2770-2013-772-s2.0-84887979143ChamkhaA. J.AlyA. M.MansourM. A.Similarity solution for unsteady heat and mass transfer from a stretching surface embedded in a porous medium with suction/injection and chemical reaction effects2010197684685810.1080/009864409033590872-s2.0-77949335296MakandaG.MakindeO. D.SibandaP.Natural convection of viscoelastic fluid from a cone embedded in a porous medium with viscous dissipation2013201311934712MR314841910.1155/2013/934712