# Effect of Internal Heat Generation/Absorption on Dusty Fluid Flow over an Exponentially Stretching Sheet with Viscous Dissipation

Authors: G. M. Pavithra; B. J. Gireesha
Journal: Journal of Mathematics
 (2013)

---

## Abstract

A numerical analysis has been carried out to describe the boundary layer flow and heat transfer of a dusty fluid over an exponentially stretching surface in the presence of viscous dissipation and internal heat generation/absorption. The governing partial differential equations are reduced to nonlinear ordinary differential equations by a similarity transformation, before being solved numerically by Runge-Kutta-Fehlberg 45 method. The heat transfer analysis has been carried out for both PEST and PEHF cases. The numerical results are compared with the earlier study and found to be in excellent agreement. Some important features of the flow and heat transfer in terms of velocities and temperature distributions for different values of the governing parameters like fluid-particle interaction parameter, Prandtl number, Eckert number, Number density, heat source/sink parameter, and suction parameter which are of physical and engineering interests are analyzed, discussed, and presented through tables and graphs.


---

## Body

## 1. Introduction

An investigation on boundary layer flow and heat transfer of viscous fluids over a moving continuous stretching surface has considerable practical applications in industries and engineering, since the study of heat transfer has become important industrially for determining the quality of final product which greatly depends on the rate of cooling. Many researchers inspired by Sakiadis [1, 2] who initiated the boundary layer behavior studied the stretching flow problem in various aspects. Extension to that, an exact solution was given by Crane [3] for a boundary layer flow caused by stretching surface.  A new dimension to the boundary layer flow was given by Magyari and Keller [4] by considering the nonstandard stretching flow known as exponentially stretching surface. They described the mass and heat transfer characteristics of the boundary layer. After that, Elbashbeshy [5] was the first who considered the heat transfer problem over an exponentially stretching sheet with suction parameter. The effect of viscous dissipation on the boundary layer flow along vertical exponential stretching sheet was explained by Partha et al. [6].Sanjayanand and Khan [7] discussed heat and mass transfer in a viscoelastic boundary layer flow over an exponentially stretching sheet. The numerical solution for the boundary layer flow with thermal radiation over an exponentially stretching sheet was given by Bidin and Nazar [8]. Samad and Mohebujjaman [9] worked on both heat and mass transfer of free convective flow along a vertical stretching sheet in presence of magnetic field with heat generation. Again Elbashbeshy and Aldawody [10] investigated the heat transfer over an unsteady stretching surface embedded in a porous medium with variable heat flux in the presence of heat source or sink. Ishak [11] and Bala and Bhaskar [12] discussed the boundary layer flow over an exponentially stretching sheet with magnetic field and thermal radiation effect. On the other hand, Sahoo and Poncet [13] carried out an exponentially stretching flow problem by considering partial slip boundary condition for third grade fluid. The flow and heat transfer of a Jeffrey fluid past an exponentially stretching sheet with thermal radiation is analyzed by Nadeem et al. [14]. Kumar [15] discussed the heat transfer over a stretching porous sheet subjected to power law heat flux in presence of heat source. Very recently Bhattacharyya [16] analyzed the effects of radiation and heat source/sink on unsteady MHD boundary layer flow and heat transfer over a shrinking sheet with suction/injection.All the works as mentioned above are restricted only for fluids induced by stretching sheet. In fact, a fluid flow with dust particles is a significant type of flow due to its wide range of its applications in the boundary layer that include the use of dust in gas cooling systems, centrifugal separation of matter from fluid, polymer technology, and fluid droplets sprays, that is, powder technology and paint spraying. The fluid-particle system was initially described by Saffman [17] and derived the motion of gas equations carrying the dust particles. Further, Datta and Mishra [18] have investigated dusty fluid in boundary layer flow over a semi-infinite flat plate. Gireesha et al. [19] have obtained the boundary layer and heat transfer of a dusty fluid flow over a stretching sheet with nonuniform heat source/sink and concluded that the temperature dependent heat sinks are better suited for cooling purposes. Also Gireesha et al. [20, 21] have discussed the boundary layer flow and heat transfer of a dusty fluid over a stretching sheet by considering the viscous dissipation for both steady and unsteady flows.On basis of previous observations, we have considered the flow over an exponentially stretching sheet. The main aim of the present investigation is to study the effect of viscous dissipation and internal heat generation/absorption on flow and heat transfer of a dusty fluid over an exponentially stretching sheet with suction parameter by taking PEST and PEHF cases. The study of heat generation or absorption in moving fluids is important in problems dealing with chemical reactions and these concerned with dissociating fluids. Flows in porous media have several applications in geothermal, oil reservoir engineering, and astrophysics. The governing boundary layer equations have been simplified using suitable similarity transformations and then have been solved numerically using Runge-Kutta-Fehlberg 45 method with the help of Maple.


## 2. Mathematical Formulation and Solution of the Problem

A steady two-dimensional laminar boundary layer flow and heat transfer of an incompressible viscous dusty fluid near a permeable plane wall stretching with velocityUw=U0e(x/L) is considered (Figure 1). The x-axis is chosen along the sheet and y-axis normal to it. Two equal and opposite forces are applied along the sheet, so that the wall is stretched exponentially.Figure 1
Schematic representation of boundary layer flow.Under these assumptions, the two-dimensional boundary layer equations can be written as(1)∂u∂x+∂v∂y=0,(2)u∂u∂x+v∂u∂y=υ∂2u∂y2+KNρ(up-u),(3)∂up∂x+∂vp∂y=0,(4)up∂up∂x+vp∂up∂y=Km(u-up),
where x and y represent coordinate axes along the continuous surface in the direction of motion and perpendicular to it, respectively. (u,v) and (up,vp) denote the velocity components of the fluid and particle phase along the x and y directions, respectively, υ is the coefficient of viscosity of fluid, ρ is the density of the fluid phase, K is Stoke’s resistance, N is the number density of dust particles, m is the mass concentration of dust particles, and τv=m/K is the relaxation time of particle phase.In order to solve the governing boundary layer equations, consider the following appropriate boundary conditions on velocity:(5)u=Uw(x),v=Vw(x)aty=0,u⟶0,up⟶0,vp⟶v,asy⟶∞,
where Uw(x)=U0e(x/L) is the sheet velocity, Vw(x)=-SU0υ/2Le(x/2L) is the suction velocity, U0 is reference velocity, and L is the reference length.Equations (1)–(4) that are subject to boundary condition (5) admit self-similar solutions in terms of the similarity function f and the similarity variable η as
(6)u=U0ex/Lf′(η),v=-u0υ2Lex/2L[f(η)+ηf′(η)],up=U0ex/LF′(η),vp=-u0υ2Lex/2L[F(η)+ηF′(η)],η=u02υLex/2Ly.These equations identically satisfy the governing equations (1) and (3). Substitute (6) into (2)–(4) and on equating the coefficient of (x∖L)0 on both sides, and then one can get
(7)f′′′(η)+f(η)f′′(η)-2f′(η)2+2lβ[F′(η)-f′(η)]=0,F(η)F′′(η)-2F′(η)2+2β[f′(η)-F′(η)]=0,
where prime denotes the differentiation with respect to η, l=mN/ρ is the mass concentration, and β=L/τvU0 is the fluid-particle interaction parameter for velocity.Similarity boundary conditions (5) will become
(8)f′(η)=1,f(η)=Satη=0,f′(η)=0,F′(η)=0,F(η)=f(η)+ηf′(η)-ηF′(η)asη⟶∞,
where S>0 is a suction parameter.The important physical parameter for the boundary layer flow is the skin-friction coefficient which is defined as(9)Cf=τwρUw2,
where the skin friction τw is given by
(10)τw=μ(∂u∂y)y=0.Using the nondimensional variables, one obtains(11)2ReCf=f′′(0),
where Re=U0L/υ is the Reynolds number.


## 3. Heat Transfer Analysis

The governing steady, boundary layer heat transport equations for both fluid and dust phases with viscous dissipation and heat generation/absorption are given by(12)ρcp[u∂T∂x+v∂T∂y]=k∂2T∂y2+NcpτT(Tp-T)+Nτv(up-u)2+μ(∂u∂y)2+Q(T-T∞),Ncm[up∂Tp∂x+vp∂Tp∂y]=-NcpτT(Tp-T),
where T and Tp are the temperatures of the fluid and dust particles inside the boundary layer, cp and cm are the specific heat of fluid and dust particles, τT is the thermal equilibrium time; that is, it is time required by a dust cloud to adjust its temperature to the fluid, k is the thermal conductivity, τv is the relaxation time of the of dust particle, that is, the time required by a dust particle to adjust its velocity relative to the fluid, and Q represents the heat source when Q>0 and the sink when Q<0.We solved the heat transfer phenomenon for two types of heating process, namely,(1)
prescribed exponential order surface temperature (PEST);(2)
prescribed exponential order heat flux (PEHF).Case 1 (prescribed exponential order surface temperature).
For this heating process, we employ the following boundary conditions:(13)T=Tw(x)aty=0,T⟶T∞,Tp⟶T∞asy⟶∞,
where Tw=T∞+T0ec1x/2L is the temperature distribution in the stretching surface, T0 is a reference temperature, and c1 is constant.
Introduce the dimensionless variables for the temperaturesθ(η) and θp(η) as follows:
(14)θ(η)=T-T∞Tw-T∞,θp(η)=Tp-T∞Tw-T∞,
where T-T∞=T0ec1x/2Lθ(η).
Using the similarity variableη and (14) into (12) and on equating the coefficient of (x∖L)0 on both sides, one can arrive at the following system of equations:
(15)θ′′(η)+Pr[f(η)θ′(η)-c1f′(η)θ(η)]+2NρβTPr[θp(η)-θ(η)]+2NρβPrEc[F′(η)-f′(η)]2+PrEc[f′′(η)]2+2Prλθ(η)=0,c1F′(η)θp(η)-F(η)θp′(η)+2βTγ[θp(η)-θ(η)]=0,
where Pr=μcp/k is the Prandtl number, Ec=U02/cpT is the Eckert number, β=L/τvU0 and βτ=L/τTU0 are the fluid-particle interaction parameter for velocity and temperature, and γ=cp/cm is the ratio of specific heat.
Corresponding thermal boundary conditions becomes(16)θ(η)=1atη=0,θ(η)⟶0,θp(η)⟶0asη⟶∞.Case 2 (prescribed exponential order heat flux).
For this heating process, consider the boundary conditions as follows:(17)∂T∂y=-qw(x)katy=0,T⟶T∞,Tp⟶T∞asη⟶∞,
where qw(x)=T1e(c2+1)x/2L, T1 is reference temperature, and c2 is constant.
Using the similarity variableη and (14) into (12) and on equating the coefficient of (x∖L)0 on both sides, one can arrive at the following system of equations:
(18)θ′′(η)+Pr[f(η)θ′(η)-c1f′(η)θ(η)]+2NρβTPr[θp(η)-θ(η)]+2NρβPrEc[F′(η)-f′(η)]2+PrEc[f′′(η)]2+2Prλθ(η)=0,c1F′(η)θp(η)-F(η)θp′(η)+2βTγ[θp(η)-θ(η)]=0,
where  Ec=kU02/cpT1U0/2υL is the Eckert number.Corresponding thermal boundary conditions becomes(19)θ′(η)=-1atη=0,θ(η)⟶0,θp(η)⟶0asη⟶∞.


## 4. Numerical Solution

Two-dimensional, boundary layer flow and heat transfer of a dusty fluid over an exponential stretching sheet is considered. The system of highly nonlinear ordinary differential equations (7), (15) for PEST case and (18) for PEHF case are solved numerically using RKF-45 method with the help of an algebraic software Maple. In this method, we choose suitable finite values of η→∞  as η=5.Here we have given the comparison of our results of-θ′(0) with Nadeem et al. [14] and Bidin and Nazar [8] as in Table 1 for various values of Pr and Ec. From this table, one can notice that there is a close agreement with these approaches and thus verifies the accuracy of the method used. The thermal characteristics at the wall that are examined for the values of temperature gradient θ′(0) in PEST case and the temperature θ(0) in PEHF case are tabulated in Table 2. The computed values of skin friction coefficient f′′(0) is tabulated in Table 3 for different values of fluid-particle interaction parameter (β) and suction parameter (S). Further, we studied the effects of viscous dissipation and heat generation/absorption on velocity and temperature profiles that are depicted graphically for different values of fluid-particle interaction parameter (β), Suction parameter (S), heat source/sink parameter (λ), Number density (N), Prandtl number (Pr), and Eckert number (Ec).Table 1
Comparison of the results for the dimensionless wall temperature gradientθ′(0) (PEST case) by varying the Pr and Ec with β=N=λ=S=0.
Pr
Ec
=
0
Ec
=
0.2
Ec
=
0.9
Nadeem et al. [14]
Bidin and Nazar [8]
Present studyθ′(0)
Nadeem et al. [14]
Bidin and Nazar [8]
Present studyθ′(0)
Nadeem et al. [14]
Bidin and Nazar [8]
Present studyθ′(0)
1
−0.955
−0.955
−0.955
−0.863
−0.862
−0.863
−0.539
−0.539
−0.539
2
−1.471
−1.471
−1.471
−1.306
−1.306
−1.305
−0.725
−0.725
−0.724
3
−1.870
−1.869
−1.869
−1.639
−1.688
−1.638
−0.830
−0.830
−0.830Table 2
Values of wall temperature gradientθ′(0) (for PEST case) and wall temperature θ(0) (for PEHF case).
β
Pr
Ec
λ
N
θ
′
(
0
)
θ
(
0
)
0.20.61.0
0.72
0.5
0.5
0.5
−1.75024−1.66464−1.64077
0.633480.674080.68511
0.6
0.721.01.5
0.5
0.5
0.5
−1.66464−2.01014−2.53037
0.674080.596340.52418
0.6
0.72
00.51.0
0.5
0.5
−2.03930−1.66464−1.28998
0.490360.674080.85780
0.6
0.72
0.5
0.51.01.5
0.5
−1.66464−1.39884−1.05000
0.674080.778240.96894
0.6
0.72
0.5
0.5
0.51.01.5
−1.66464−2.39962−2.94768
0.674080.514310.44797Table 3
Values of skin friction coefficientf′′(0).
β
S
f
′
′
(
0
)
0.20.61.0
0.5
−1.55365−1.56928−1.57897
0.6
0.51.01.5
−1.56928−1.88171−2.23651


## 5. Results and Discussion

Numerical calculations are performed for velocity and temperature profiles for various values of physical parameters such as fluid-particle interaction parameter (β), Suction parameter (S), heat source/sink parameter (λ), Number density (N), Prandtl number (Pr), and Eckert number (Ec) and are depicted graphically (from Figure 2 to Figure 11). Comparison values of wall-temperature gradient are tabulated in Table 1. From Table 3 we observed that the skin friction coefficient decreases with increasing values of fluid-particle interaction parameter (β) as well as for the suction parameter (S). Physically negative values of f′′(0) mean that the surface exerts a drag force on the fluid, so that the stretching surface will induce the flow.Figure 2
Effect ofβ on velocity profiles with N=λ=Ec=S=0.5 and Pr=0.72.From Figure2, the observation shows that an increase in fluid-particle interaction parameter β decreases the fluid velocity f′(η) and increases the particle velocity F′(η), and also as β increases, the velocity of fluid will be equal to the velocity of the dust particle. Figure 3 represents the velocity profiles for different values of suction parameter S. It shows that the velocity of fluid and dust decreases when suction parameter increases. So that the momentum boundary layer thickness becomes thinner.Figure 3
Effect ofS on velocity profiles with N=λ=Ec=0.5, Pr=0.72, and β=0.6.Figure4 depicts the temperature profiles θ(η) and θp(η) versus η for different values of fluid-particle interaction parameter β. We infer from this figure that the temperature increase with increase in the fluid-particle interaction parameter β, and it indicates that the fluid-particle temperature is parallel to that of dust phase. Also one can observ that fluid phase temperature is higher than that of dust phase. From Figure 5, the observation shows the effect of heat source/sink parameter λ on temperature profiles. As λ increases, temperature profiles for both fluid and dust phases increase in both PEST and PEHF cases. It is clear that the temperature in the case of heat source is higher than in the case of sink. This is very much significant in which the heat transfer is given prime importance.Effect ofβ on temperature profiles for both PEST and PEHF cases with  N=λ=Ec=S=0.5 and Pr=0.72.
(a)
(b)Effect ofλ on temperature profiles for both PEST and PEHF cases with N=Ec=S=0.5, Pr=0.72, and β=0.6.
(a)
(b)The effects of suction parameterS on the temperature profiles are depicted as in Figure 6. This figure explains that the temperature will decrease as S increases. This is due to the fact that the fluid at ambient conditions is brought closer to the surface which results in thinning of thermal boundary layer thickness. This causes an increase in the rate of heat transfer. So, suction can be used as a means for cooling the surface. Figure 7 shows the temperature distributions θ(η) and θp(η) versus η for different values of number density N. We infer from this figure that the temperature decreases with increases in N for both cases.Effect ofS on temperature profiles for both PEST and PEHF cases with N=Ec=λ=0.5, Pr=0.72, and β=0.6.
(a)
(b)Effect ofN on temperature profiles for both PEST and PEHF cases with S=Ec=λ=0.5, Pr=0.72, and β=0.6.
(a)
(b)The temperature field for various values of the Prandtl number (Pr) is represent in Figure 8. The relative thickening of momentum and thermal boundary layers is controlled by Prandtl number (Pr). Since small values of Pr will possess higher thermal conductivities, so that heat can diffuse from the sheet very quickly compared to the velocity. From this figure, it reveals that the temperature decreases with increase in the value of Pr. Hence Prandtl number can be used to increase the rate of cooling. Figure 9 explains the effect of viscous dissipation on temperature profiles. Viscous dissipation changes the temperature distribution by playing a role like an energy source, which leads to affect heat transfer rates. Here the temperature increases with increase in the value of Ec, due to the heat energy that is stored in the liquid and frictional heating, and this is true in both cases.Effect ofPr on temperature profiles for both PEST and PEHF cases with S=Ec=N=λ=0.5 and β=0.6.
(a)
(b)Effect ofEc on temperature profiles for both PEST and PEHF cases with S=N=λ=0.5, Pr=0.72, and β=0.6.
(a)
(b)Figures10(a) and 10(b) show the effect of internal heat source/sink parameter in the presence and in the absence of suction parameter S. This shows that there is an increase of temperature as λ increases, which results in the reduction of the thermal boundary layer thickness. It is also noticed that the temperature is less in presence of suction parameter than in absence of suction parameter.(a) Effect ofλ on temperature profiles for PEST case; (b) effect of λ on temperature profiles for PEHF case.
(a)
(b)Heat transfer characteristics for different values ofPr and λ for both PEST and PEHF cases.
(a)
(b)The rate of heat transfer from the sheet that is evaluated by the variation of wall temperature gradientθ′(0) at sheet is present in Figure 11 for various values of Pr and λ. It is seen from this figure that the rate of heat transfer decreases with increases in Pr. It is also evident that θ′(0)  which is negative means heat transfer and θ(0) which is positive means heat absorption, and also it is clear from Table 2. We have used the values of βT=0.6, c1=1, l=0.1, and ρ=1 throughout our thermal analysis.


## 6. Conclusions

The present work deals with the boundary layer flow and heat transfer of a steady dusty fluid over an exponential stretching sheet with viscous dissipation and heat generation or absorption for both PEST and PEHF cases. The set of nonlinear ordinary differential equations (7), (15) for PEST case or (18) for PEHF case is solved numerically by applying RKF-45 order method using the software Maple. The results of the thermal characteristics at the wall are examined for the values of temperature gradient function θ′(0) in PEST case and the temperature function θ(0) in PEHF case, which are tabulated in Table 2. Also the results of skin friction coefficient f′′(0) are tabulated in Table 3 for various values of fluid particle interaction parameter (β) and suction parameter (S). The velocity and temperature profiles are obtained for various values of physical parameters like fluid-particle interaction parameter (β), Suction parameter (S), heat source/sink parameter (λ), Number density (N), Prandtl number (Pr), and Eckert number (Ec). The numerical results obtained are agrees with previously reported cases available in the literature [8, 14].The major findings from the present study can be summarized as follows(i)
Suction parameter reduces the velocity and temperature profiles for both PEST and PEHF cases.(ii)
Heat source/sink effect is less in permeable than in impermeable stretching sheet.(iii)
The PEHF boundary condition is better suited for effective cooling of the stretching sheet.(iv)
The effect of increasing the values ofβ, Ec, and λ is to increase the wall temperature gradient function θ′(0) and wall temperature function θ(0) and decrease the increasing values of Pr and N both PEST and PEHF cases.(v)
The effect of heat source/sink on temperature is quite opposite to that of suction parameter.(vi)
Ifβ→0, N→0, and S→0, then our results coincide with the results of Nadeem et al. [14] and Bidin and Nazar [8] for different values of Prandtl and Eckert numbers.(vii)
The effect ofEc increases, while Pr decreases the thermal boundary layer thickness.(viii)
Fluid phase temperature is higher than that of dust phase.



---

*Source: 583615-2013-06-25.xml*