This paper investigates the entropy generation in a third-grade fluid flow with variable properties through a channel. Approximate solutions to the nonlinear boundary-value problem are obtained using Adomian decomposition method (ADM). Variation of important parameters on the fluid velocity, temperature distribution, entropy generation and irreversibility ratio are presented graphically and discussed.
1. Introduction
Over the past few decades, there has been a tremendous increase in the study of heat transfer to viscous fluids due to its numerous applications in many industrial and engineering processes. As postulated in the second law of thermodynamics, the study could give insight into the thermal performance of the system by improving energy that is available for work [1, 2].
In the class of third-grade fluid, quite a lot has been done in recent times on entropy generation by assuming a constant thermal conductivity. For instance, Pakdemirli and Yilbas investigated the entropy generation in the flow of third-grade fluid through a pipe with constant viscosity in [3] while the temperature dependent viscosity has been investigated in [4] by using Vogel model. We refer interested readers to [5–14] for more interesting results on third-grade fluid flow. Moreover, Makinde and Aziz [15] presented the second law analysis of a pressure-driven temperature dependent fluid flow with asymmetry at the walls. Kahraman and Yürüsoy [16] examined the entropy generation due to non-Newtonian fluid flow in an annular pipe with relative rotation using a third-grade fluid model while Chauhan and Kumar [17] presented a non-Newtonian third-grade fluid flow in an annulus partially filled by a porous medium of very small permeability and many more results on entropy generation in literature.
Surprisingly, in spite of the enormous amount of work done on the entropy generation, it is observed that not much has been done on the exergy analysis of third-grade fluid flow with variable thermal conductivity and internal heat generation. As shown by Hayat et al. [18], these variations can significantly affect the flow field. Hence a more accurate result could be obtained by taking these variations into consideration.
Therefore, the specific objective of this paper is to investigate the entropy generation in the flow of third-grade fluid channel flow with temperature dependent properties. The problem under discussion is strongly nonlinear boundary valued problem. Approximate solution will be obtained using Adomian decomposition method as presented by Siddiqui et al. [19] for the third-grade fluid flow. Therefore, in the absence of internal heat generation and neglecting variations with temperature, the result obtained in [19] will be fully recovered. The seminumerical method has also been used to obtain approximate solution to different fluid flow problems in [20–22]. To the best of the author’s knowledge, the study reported here on entropy generation in a heat generating third-grade fluid flow with variable properties has not been undertaken in literature.
The rest of the paper is organized as follows. In the following section, the mathematical analysis of the flow is presented. Section 3 describes the method of obtaining the approximate solution of the nonlinear problem. In Section 4, results are presented and discussed before the paper is concluded in Section 5.
2. Mathematical Analysis
Consider the steady flow of a third-grade fluid with variable viscosity and thermal conductivity through infinite parallel plates of distance 2a apart as shown in Figure 1.
Channel geometry.
The fluid is assumed to be reactive and internal heat generation is assumed to be a linear function of temperature. Then the governing equations for the fully developed flow can be written as [15, 18]
(1)0=-dPdx+ddy′(μ′du′dy′)+6β3d2u′dy′2(du′dy′)2,0=ddy′(k′dTdy′)+μ′(du′dy′)2+2β3(du′dy′)4+Q0(T-T0),
and the entropy generation is given as
(2)EG=k′T02(dTdy′)2+1T0(du′dy′)2(μ′+2β3(du′dy′)2)
with appropriate boundary conditions
(3)u′=0,T=T0,ony′=a,du′dy′=0=dTdy′,ony′=0;
the temperature dependent viscosity and thermal conductivity take the form
(4)μ′=μ-η0(T-T0),k′=k+η1(T-T0),
where u′ is the fluid velocity, T is the fluid temperature, P is the fluid pressure, a is the channel half width, β3 is the material coefficient, k is the thermal conductivity, μ′ is the dynamic viscosity, η0,1 are the viscosity and thermal conductivity variation parameters, Q0 is the heat generated internally, T0 and T1 are referenced fluid temperature and EG is the entropy generation.
Introducing the following dimensionless parameters and variables:
(5)y=y′a,u=u′U,γ=β3U2a2μ,θ=T-T0T1-T0,α=η0(T1-T0)μ,Ω=T1-T0T0,G=-a2μUdPdx,δ=Q0a2k,Br=μU2k0(T1-T0),λ=η1(T1-T0)k,Ns=T02a2EGk(T1-T0)2,
we get the following dimensionless equations together with appropriate boundary conditions:
(6)d2udy2=α[θd2udy2+dudydθdy]-6γd2udy2(dudy)2-G;kkkkkkkkkkkkkkkkkkkkkkkkkgikkiu(±1)=0,d2θdy2=-λ[(dθdy)2+θd2θdy2]d2θdy2=-Br(dudy)2{(1-αθ)+2γ(dudy)2}+δθ;kkkkkkkkkkkkkkkkkkkkkkkkkkkkkθ(±1)=0,(7)Ns=(1+λθ)(dθdy)2+BrΩ(dudy)2{(1-αθ)+2γ(dudy)2}.
If we define
(8)N1=(1+λθ)(dθdy)2,N2=BrΩ(dudy)2{(1-αθ)+2γ(dudy)2},
then the irreversibility ratio denoted by Bejan number (Be) can be written as
(9)Be=N1Ns=N1N1+N2=11+Φ,Φ=N2N1,
here u is the dimensionless velocity, α is the viscosity variation parameter, γ is the third-grade material effect, G is the pressure gradient, θ is the dimensionless fluid temperature, λ is the thermal conductivity variation parameter, Br is Brinkman number, δ is the internal heat generation parameter, Ω is the parameter that measures the temperature difference between the two heat reservoirs, Ns is the dimensionless entropy generation rate, and U represents the characteristic velocity.
3. Adomian Decomposition Method of Solution
The integral form of (6) can be written as
(10)u(y)=a0+∬0y{α[θd2udY2+dudYdθdY](dudY)2kkkkkkkkkkkkkkki-6γd2udY2(dudY)2-G}dYdY,θ(y)=b0-∬0y{λ[(dθdY)2+θd2θdY2]+Br(dudY)2kkkkkkkikkkkkk×{(1-αθ)+2γ(dudY)2}+δθ}dYdY,
where the constants a0 and b0 are to be determined later; the standard Adomian decomposition method [19–21] assumes series solutions in the form
(11)u=∑n=0∞un,θ=∑n=0∞θn;
substituting (11) into the integral equations (10) yields the following recursive algorithm:
(12)θ0(y)=b0,u0(y)=a0-∬0y(G)dYdY,un+1(y)=∬0y{αAn-6γBn}dYdY,n≥1,θn+1(y)=-∬0y[λCn+BrFn{(1-αθn)+2γFn}+δθn]dYdY;kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkn≥1;
the approximate solution is given by the partial sum
(13)θ=∑n=0mθn,u=∑n=0mun,
where m represent the truncation point. The nonlinear terms in (12) represented by
(14)An=dundYdθndY+θnd2undY2,Bn=d2undY2(dundY)2,Cn=(dθndY)2+θnd2θndY2,Fn=(dundY)2
are decomposed into Adomian polynomials as follows:
(15)A0=dθ0dYdu0dY+θ0d2u0dY2,A1=dθ0dYdu1dY+dθ1dYdu0dY+θ1d2u0dY2+θ0d2u1dY2,A2=dθ0dYdu2dY+dθ1dYdu1dY+dθ2dYdu0dY+θ2d2u0dY2A2=+θ1d2u1dY2+θ0d2u2dY2,
similarly,
(16)B0=d2u0dY2(du0dY)2,B1=2d2u0dY2(du0dY)(du1dY)+d2u1dY2(du0dY)2,B2=d2u0dY2(du1dY)2+2d2u0dY2(du0dY)(du2dY)+2d2u1dY2(du0dY)(du1dY)+d2u2dY2(du0dY)2,
with
(17)F0=(du0dY)2,F1=2du0dYdu1dY,F2=2du0dYdu2dY+(du1dY)2,…,(18)C0=(dθ0dY)2+θd2θ0dY2,C1=2dθ0dYdθ1dY+θ0d2θ1dY2+θ1d2θ0dY2,C2=2dθ0dYdθ2dY+(dθ1dY)2+θ2d2θ0dY2+θ1d2θ1dY2+θ0d2θ2dY2,….
4. Results and Discussion
The major objective of this work is to investigate the effects of temperature dependent viscosity and thermal conductivity on heat generating third-grade fluid flow through a channel with isothermal temperature. Approximate solutions to the nonlinear boundary-value problem are obtained in the form of a rapidly convergent Adomian series solution. Table 1 confirmed that series solution is convergent and reliable. Figure 2 shows the velocity profile for variations in the viscosity variation parameter. As observed from the graph, an increase in the viscosity variation parameter enhances the flow. This is physically true due to the melting effect on the fluid viscosity. Similarly, Figure 3 shows that an increase in the viscosity variation parameter increases the temperature profile due to rise in the temperature difference within the channel. Figures 4 and 5 show the velocity and temperature profiles, respectively, for variation in Brinkman number. From the plots it is observed that an increase in the Brinkman number leads to an increase in both the fluid velocity and temperature due to rise in the kinetic energy within the moving fluid. Moreover, Figures 6 and 7 represent the effect of third-grade material effect on the velocity and temperature profiles; the result shows that an increase in the non-Newtonian material effect decreases both the velocity and temperature profiles due to fluid thickening. Figure 8 depicts the velocity profile for variations in internal heat generation parameter. The result shows that an increase in heat generation parameter increases the flow velocity due to rise in the concentration of the reacting fluid. However, as observed in Figure 9, an increase in the internal heat generation increases the temperature distribution within the channel. Figures 10 and 11 display the effect of variations in thermal conductivity on both the velocity and temperature distributions within the channel. As observed from Figure 10, an increase in the thermal conductivity variation parameter increases the velocity distribution within the channel. On the other hand, an increase in the thermal conductivity variation parameter implies a decrease in the thermal conductivity of the fluid. This ultimately decreases the fluid temperature as observed in Figure 11. As observed from the plots, Figures 12, 13, and 14 show that, in an increase in thermal conductivity variation parameter, Brinkman number and third-grade material effect enhance the entropy generation rate within the channel, while increase in viscosity variation parameter decreases the entropy generation rate within the channel as seen in Figure 15. In Figure 16, rise in the internal heat generation parameter is observed to increase entropy generation in the center line of the channel while entropy generation decreases at the walls. Moreover, Figures 17 and 18 shows that as the internal heat generation and Brinkman number increases heat transfer dominates over the fluid viscosity within the channel. Figure 19 represents the effect of viscosity variation parameter on the irreversibility ratio; from the plot it is observed that as the viscosity variation parameter increases, heat transfer dominates over the fluid viscosity at the walls whereas at the centerline of the channel fluid viscosity is observed to dominate over heat transfer. Finally, Figures 20 and 21 show that rise in increase in both third-grade material effect and thermal conductivity variation parameter shows that fluid viscosity dominates over the heat transfer within the channel.
Convergence result for γ=0.2,δ=1,Br=0.4, α=0.5=λ,y=0.1.
n
un
∑n=0mun
θn
∑n=0mθn
0
0.445252
0.445252
0.0507153
0.0507153
1
−0.000116788
0.445135
−0.000256831
0.0504585
2
-1.96106×10-6
0.445134
6.56661×10-6
0.0504650
3
-2.58151×10-8
0.445133
-1.93224×10-7
0.0504648
4
1.74273×10-10
0.445133
5.79716×10-9
0.0504648
Effect of viscosity variation parameter on velocity profile.
Effect of viscosity variation parameter on temperature profile.
Effect of Brinkman number on velocity profile.
Effect of Brinkman number on temperature profile.
Effect of non-Newtonian material on velocity profile.
Effect of non-Newtonian material on temperature profile.
Effect of internal heat generation parameter on velocity profile.
Effect of internal heat generation parameter on temperature profile.
Effect of thermal conductivity variation parameter on velocity profile.
Effect of thermal conductivity variation parameter on temperature profile.
Effect of thermal conductivity variation parameter on entropy generation rate.
Effect of Brinkman number on entropy generation rate.
Effect of third-grade material effect on entropy generation rate.
Effect of viscosity variation parameter on entropy generation rate.
Effect of internal heat generation parameter on entropy generation rate.
Effect of internal heat generation on the irreversibility ratio.
Effect of Brinkman number on the irreversibility ratio.
Effect of viscosity variation parameter on the irreversibility ratio.
Effect of third-grade material effect on irreversibility ratio.
Effect of thermal conductivity variation parameter on irreversibility ratio.
5. Conclusion
In this paper, the entropy generation analysis for third-grade fluid flow with variable properties through a channel with uniform wall temperature is studied. In the limiting case, there is perfect agreement between the present result and that obtained in [19] when λ=δ=α=0. Summarily, the main contributions to knowledge from the computation are as follows.
Rise in the viscosity variation parameter is observed to decrease the entropy generation rate within the channel. Moreover, in the middle of the channel, fluid viscosity is observed to dominate over heat transfer. A reverse behaviour is observed at the channel walls.
An increase in the fluid thermal conductivity enhances entropy generation rate. Interestingly, it supports the dominance of fluid viscosity over the heat transfer within the channel.
An increase in the internal heat generation parameter increases the entropy generation rate in the centreline of the channel only.
One major area of future research is the flow behavior for large values of the thermal conductivity parameter. To achieve this further analysis will be required to determine the inherent thermal criticality.Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
AdesanyaS. O.MakindeO. D.Entropy generation in couple stress fluid flow through porous channel with fluid slippageInternational Journal of Exergy. In pressAdesanyaS. O.MakindeO. D.Effects of couple stresses on entropy generation rate in a porous channel with convective heatingPakdemirliM.YilbasB. S.Entropy generation in a pipe due to non-Newtonian fluid flow: constant viscosity casePakdemirliM.YilbasB. S.Entropy generation for pipe flow of a third grade fluid with Vogel model viscosityEllahiR.A study on the convergence of series solution of non-Newtonian third grade fluid with variable viscosity: by means of homotopy analysis methodEllahiR.HayatT.MahomedF. M.Generalized couette flow of a third-grade fluid with slip: the exact solutionsEllahiR.RiazA.Analytical solution for MHD flow in a third-grade fluid with variable viscosityEllahiR.HayatT.MahomedF. M.AsgharS.Effects of slip on the non-linear flows of a third grade fluidAbbasbandyS.HayatT.MahomedF. M.EllahiR.On comparison of exact and series solutions for thin film flow of a third-grade fluidEllahiR.ArielP. D.HayatT.AsgharS.Effect of heat transfer on a third grade fluid in a flat channelHayatT.EllahiR.MahomedF. M.The analytical solutions for magnetohydrodynamic flow of a third order fluid in a porous mediumHayatT.EllahiR.MahomedF. M.Exact solutions for thin film flow of a third grade fluid down an inclined planeEllahiR.AfzalS.Effects of variable viscosity in a third grade fluid with porous medium: an analytical solutionEllahiR.The thermodynamics, stability, applications and techniques of differential type: a reviewMakindeO. D.AzizA.Second law analysis for a variable viscosity plane Poiseuille flow with asymmetric convective coolingKahramanA.YürüsoyM.Entropy generation due to non-newtonian fluid flow in annular pipe with relative rotation: constant viscosity caseChauhanD. S.KumarV.Entropy analysis for third-grade fluid flow with temperature-dependent viscosity in annulus partially filled with porous mediumHayatT.ShehzadS. A.AsgharS.MHD flow, thixotropic fluid, thermal radiation, variable thermal conductivity and thermal radiationSiddiquiA. M.HameedM.SiddiquiB. M.GhoriQ. K.Use of Adomian decomposition method in the study of parallel plate flow of a third grade fluidAdesanyaS. O.Linear stability analysis of a plane-Poiseuille hydromagnetic flow using Adomian decomposition methodAdesanyaS. O.Steady magnetohydrodynamic visco-elastic heat generating /absorbing slip flow through a porous medium with radiation effectAdesanyaS. O.MakindeO. D.Heat transfer to magnetohydrodynamic non-newtonian couple stress pulsatile flow between two parallel porous plates