# Perturbation Solutions for Hagen-Poiseuille Flow and Heat Transfer of Third-Grade Fluid with Temperature-Dependent Viscosities and Internal Heat Generation

Authors: B. Y. Ogunmola; A. T. Akinshilo; M. G. Sobamowo
Journal: International Journal of Engineering Mathematics
 (2016)

---

## Abstract

Regular perturbation technique is applied to analyze the fluid flow and heat transfer in a pipe containing third-grade fluid with temperature-dependent viscosities and heat generation under slip and no slip conditions. The obtained approximate solutions were used to investigate the effects of slip on the heat transfer characteristics of the laminar flow in a pipe under Reynolds’s and Vogel’s temperature-dependent viscosities. Also, the effects of parameters such as variable viscosity, non-Newtonian parameter, viscous dissipation, and pressure gradient at various values were established. The results of this work were compared with the numerical results found in literature and good agreements were established. The results can be used to advance the analysis and study of the behavior of third-grade fluid flow and steady state heat transfer processes such as those found in coal slurries, polymer solutions, textiles, ceramics, catalytic reactors, and oil recovery applications.


---

## Body

## 1. Introduction

The third-grade fluid being a subclass of non-Newtonian fluid whose viscosity varies based on the applied force is a favored fluid due to the exciting phenomena it captures such as the shear thinning and thickening effects. The mathematical model of a third-grade fluid represents a more realistic description of the behavior of non-Newtonian fluids. The model also represents a further attempt towards the study of the flow structure of non-Newtonian fluids. Consequently, considerable interest has been shown in the third-grade fluid over the past few decades by various researchers due to its potential application in industry and technology. It should be stated that the governing equations for the third-grade fluid model are nonlinear and much more complicated than those of Newtonian fluids. They require additional boundary conditions to obtain a physically meaningful solution. This issue has been discussed in detail by researchers [1–23] who made a complete thermodynamical analysis of a third-grade fluid and derived the restriction on the stress constitutive equation. They investigated some stability characteristics of third-grade fluids and showed that they exhibit features different from those of Newtonian and second-grade fluids. In an attempt to improve the characteristics of the viscoelastic properties of the fluid, the stability of the third-grade fluid model was studied by Fosdick and Rajagopal [1] while Majhi and Nair [2] investigated the effects of stenotic geometry and the non-Newtonian parameter of the third-grade fluid on the resistive impedance and wall shearing stress. Their results were compared with similar study submitted by Massoudi and Christie [3] who presented the numerical solutions on the effects of variable viscosity and viscous dissipation on the flow of third-grade fluid in a pipe using the finite difference method. Yürüsoy and Pakdemirli [4] developed approximate analytical solution for the flow of the third-grade fluid in a pipe using constant viscosity model and variable viscosity models under no slip condition while Vajravelu et al. [5] presented a numerical solution for the third-grade fluid flows between rotating cylinder using Schauder theory and perturbation technique. The fluctuating flow of a magnetohydrodynamic rotational flow of the third-grade fluid on a porous plate was studied by Hayat et al. [6]. The similarity solutions to boundary layer equation for the third-grade fluid were developed by Muhammet [7] using the special coordinate system generated by the potential flow. The steady flow analysis of the third-grade fluid between circular concentric cylinders with heat transfer was presented by Yürüsoy [8]. In the study, the pipe temperature is assumed to be higher than fluid temperature. Also, Pakdemirli and Yilbas [9] developed approximate analytical solution of non-Newtonian fluid using Vogel’s viscosity model and entropy generation in a pipe. The steady flow of a third-grade fluid past a porous horizontal plate with partial slip was investigated by Sajid et al. [10]. Ellahi et al. [11] used implicit finite difference method to analyze the unsteady free convective flow of a third-grade fluid past an infinite vertical plate when uniform suction is applied at the plate. Jayeoba and Okoya [12] presented the analytical approximate solution to determine the temperature fields for steady flow of a third-grade fluid in a pipe with models of viscosities including a heat generation term for the no slip boundary condition. Most of the studies previously carried out on the third-grade fluid are limited to no slip flow condition which is a simplified method of predicting the actual behavior of the fluid in various applications. In practice, most problems of fluid flow exist as either partial slip or slip condition, Hence, in this work, regular perturbation method is used to develop approximate analytical solutions that are used to investigate the effects of slip on the heat transfer analysis in pipe flow of third-grade fluid with temperature-dependent viscosities and internal heat generation.


## 2. Model Development and Analytical Solution

Consider the flow of a third-grade fluid in a pipe as shown in Figure1. The momentum and temperature equations with incorporated source term are given by [3, 4, 13](1)1rddrr2∝1+∝2dwdr2=∂p∂r,(2)0=∂p∂∅,(3)1rddrrμdwdr+1rddr2rβ3dwdr3=∂p∂z,(4)k1rddrrdTdr+μdwdr2+2β3dwdr4+QC0T-T0=0.The slip condition at the pipe wall can be introduced in terms of shear stress (5)wR=γdwdr+2β3μdwdr3r=R,hTR-Ts=-kdθdr+μdwdr2+2β3dwdr4r=R,which makes the boundary condition nonlinear. Linearizing the boundary condition at the wall gives a corresponding linear boundary condition at the wall under slip condition as (6)w1=γldwdrr=R,hTR-Ts=-kldθdrr=R.At the center of the pipe, the boundary condition is given as(7)dwdrr=0=0,dTdrr=0=0.The boundary conditions under slip condition are given as (8)wR=APdwdrr=R,hTR-Ts=-kldTdrr=R.With the dimensionless parameters as stated in the nomenclature, the dimensionless equations (leaving out the bars on the equations) for (3) and (4) are given as (9)dμdrdwdr+μrdwdr+rd2wdr2+Λrdwdr2dωdr+3rd2wdr2=C,(10)d2θdr2+1rdθdr+Γdwdr2μ+Λdwdr2+δθ=0.And the dimensionless boundary conditions under slip condition are given as (11)dwdr0=0,dθdr0=0,w1=APdwdrr=1/D,θ1=-1Bidθdrr=1/D,where D=1+8Kn and Ap=-2Kn.Figure 1
Physical model and coordinate system [12].In this work, two equations of temperature-dependent viscosity are used, Reynolds and Vogel’s models.(12a)μT=μ0e-MT-T0Reynolds  model(12b)μT=μ0ea/b-T0Vogel's  model.Reynolds model is a model of viscosity varying constantly with temperature while Vogel’s model of viscosity is varying exponentially with temperature. The nondimensionalized form of both models yields(13a)μ=e-YθReynolds  model(13b)μ=eA/B-θ-T0Vogel's  model.


### 2.1. Reynolds Model of Viscosity

The regular perturbation technique is used to determine the approximate analytical solution for the Reynolds viscosity model.Taking Maclaurin’s series, (13a) takes the form(14)μ=exp⁡-Yθ=1-Yθ+OY2.Series solutions of the equation of the velocity and temperature may be obtained by using perturbation method taking ϵ as the perturbation parameter. Solutions are obtained in the following form:(15)w=w0+ϵw1+Oϵ2,θ=θ0+ϵθ1+Oϵ2.Substituting (15) into (9) and changing parameters into terms yield(16)Oϵ0:rd2w0dr2+dw0dr=Cr,(17)Oϵ1:dw1dr+rd2w1dr2-Pθ0dw0dr-Pθ0d2w0dr2-Prdw0drdθ0dr=-Mdw0dr2dw0dr+3rd2w0dr2+Cr.Also, substituting (14) and (15) into (10) and changing parameters into terms yield(18)Oϵ0:rd2θ0dr2+dθ0dr+dθ0dr+Γrdw0dr2=0,(19)θϵ1:rd2θ1dr2+dθ0dr+dθ1dr+Γrdw0dr2·Mdw0dr2-Pθ0+2Γrdw0drdw1dr+rNθ0=0.The boundary conditions for the leading order equation are(20)dw0dr0=dθ0dr0=0,w01=APdw0drr=1/D,θ01=-1Bidθ0drr=1/D,where Λ=ϵM, δ=ϵN, γ=ϵP.With the boundary conditions, it could be easily shown that (16) and (18) are given as(21)w0=C4r2+2ApD-1D2,θ0=ΓC2641-r4+4BiD3.Also, for the first-order equation, the boundary conditions are given as (22)dw1dr0=dθ1dr0=0,w11=APdw1drr=1/D,θ11=-1Bidθ1dr,where D=1+8Kn and Ap=-2Kn.Substituting the corresponding terms from the solutions in (21) into (17) and (19) and integrating the resulting equations two times successively, with the applications of the boundary conditions in (22) yields(23)w1=APPΓC36412D-16D5+2BiD4-APPΓC3192D2-APMC58D3-PΓC31152+PΓC364r24-r636+r2BiD3-14D2+136D6-1BiD5-MC332r4-1D4,θ1=MΓC496BiD5-PΓ2C4256Bi14D3-18D7+1BiD6+PΓ2C446Bi18D3-148D7+12BiD6-PΓ2C4153BiD7-MC4Γ48Bid5+ΓC2N64Bi12D-16D5+2BiD4-MΓC4576r4-1D6+PΓ2C4256r416-r864+r44BiD3-116D4+164D8-114BiD7-PΓ2C464r432-r8384+r36BiD3-132D4+1384D8-16BiD6+PΓ2C412288r8-1D8+MC4Γ288r6-1D6-ΓC2N64r24-r636+r2BiD3-14D+136D6-1BiD5.Substituting (21) and (23) into (15) the series solution at the order of expansion and then changing terms back to original parameters finally give(24)w=C4r2+2APD-1D2+γAPΓC36412D-16D5+2BiD4-γAPΓC3192D5-ΛAPC38D3+γΓC364r24-r636+r2BiD3-14D2+136D6-1BiD5-γΓC31152r6-1D6-ΛC332r6-1D6,(25)θ=ΓC2641-r4+4BiD3+ΛΓC496BiD5-γΓ2C4256Bi14D3-18D7+1BiD6+γΓ2C464Bi18D3-148D7+12BiD5-γΓ2C41536BiD7-ΛC4Γ48BiD5+∂ΓC264Bi12D-16D5+2BiD4-ΛΓC4576r6-1D6+YΓ2C4256r416-r864+r44BiD3-116D4+164D8-114BiD7-YΓ2C464r432-r8384+r36BiD3-132D4+1384D8-16BiD6+PΓ2C412288r8-1D8+ΛC4Γ288r6-1D6-∂ΓC264r24-r636+r2BiD3-14D+136D6-1BiD5.


### 2.2. Vogel’s Model of Viscosity

With the aid of Maclaurin series, Vogel’s viscosity model can be written as(26)μ=exp⁡AB-T01-ϵAθB2+Oϵ2=exp⁡AB-T01-ϵAθB2.Taking the series solution of the velocity and temperature fields yields the expansion (27)w=w+ϵw1+Oϵ2,(28)θ=ϵθ0+ϵ2θ1+Oϵ3.Therefore substituting (26) and (27) into (9) and changing parameters into term yield(29)Oϵ0:rd2w0dr2+dw0dr=C∗r,where C∗=C/exp⁡A/B-T0,(30)Oϵ0:rd2θ0dr2+dw1dr=AB2θ0rd2w0dr2+dw0dr+AB2rdθ0drdw0dr-MC∗Cdw0dr23rd2w0dr2+dω0dr,where Λ=ϵM,δ=ϵN,Γ=ϵQ.Substituting (26), (27), and (28) into (10) and changing parameters into term yield(31)Oϵ1:rd2θ0dr2+dθ0dr+Qrdw0dr2CC∗=0,where C∗=C/exp⁡A/B-T0,  Λ=ϵM,  δ=ϵN,  Γ=ϵQ, (32)θϵ2:rd2θ1dr2+dθ1dr=QrCC∗AB2θ0dw0dr2-2rQCC∗dw0drdw1dr-rQMdw0dr4-Nrθ0=0.The boundary conditions for the leading order are as follows:(33)dw0dr0=dθ0dr0=0,w01=APdw0drr=1/D,θ01=-1Bidθ0drr=1/D.Using the boundary condition, it could be easily shown that (29) and (30) are given as(34)w0=C∗4r2+2ApD-1D2,(35)θ0=QCC∗644D3Bi-r4+1D4.Also, the first-order boundary condition is as follows:(36)dw1dr0=dθ1dr0=0,w11=APdw1drr=1/D,θ11=-1Bidθ1dr,where D=1+8Kn and Ap=-2Kn.Substituting the corresponding terms from the solutions in (33) and (34) into (30) and (32) and integrating the resulting equations two times successively, with the applications of the boundary conditions in (36) yields(37)w1=APAQCC∗264B22D4Bi-16D5+12D5+APAQCC∗192B2D5+ApMC∗48CD3+AQCC∗264B2r2D3Bi-r636+r24D4-1D5Bi+136D6-14D6+AQCC∗1152B2r6-1D6+MC∗432Cr4-1D4,(38)θ1=Q2C2C∗A256B2·Bi1D6Bi-18D7+14D4+AQ2C2C∗264B2·Bi12D6Bi-148D7+18D7+AQ2C2C∗1536B2·BiD7+MCC∗4q48CBi·D5+C∗4QM96Bi·D5+QCC∗N64Bi2D4Bi-16D5+12D5+Q2CC∗2A256B2r44D3Bi-r864+r416D4-14D7Bi+164D8-116D8-AQ2C2C∗264B2r48D3Bi-r8384+r432D4-18D7Bi+138D8-132D8-AQ2C2C∗12288B2r8-1D8-MCC∗4Q288Cr6-1D6-C∗4QM576r6-1D6-QCC∗N64r2D3Bi-r636+r24D4-1D5Bi+136D6-14D6.Substituting (34) and (37), (35), and (38) into (27) and (28), the series solution at the order of expansion, and then changing terms back to original parameters finally give(39)w=C∗4r2+2APD-1D2+APΓACC∗264B22D4Bi-16D5+12D5+APΓACC∗192B2D5+APΛC∗48CD3+APΓACC∗264B2r2D3Bi-r636+r24D4-1D5Bi+136D6-14D6+AΓCC∗1152B2r6-1D6+ΛC∗432Cr4-1D4,θ=ΓCC∗644D3Bi-r4+1D4+Γ2C2C∗2A256B2·Bi1D6Bi+18D2+14D4+Γ2AC2C∗264B2·Bi12D6Bi-148D7+18D7+Γ2AC2C∗1536B2·BiD7+ΛC∗4Γ48Bi·D5+C∗4ΓΛ96Bi·D5+ΓCC∗∂64Bi2D4·Bi-16D5+12D5+Γ2C2C∗2A256B2r44D3Bi-r864+r416D4-14D7Bi+164D8-116D8Γ2AC2C∗21228B2r8-1D8-ΛC∗4288r6-1D6-ΓC∗4Λ576r6-1D6-ΓCC∗∂64r2D3Bi-r636+r24D4-1D5Bi-136D6-14D6,where(40)C∗=cexp⁡A/B-T0.


## 2.1. Reynolds Model of Viscosity

The regular perturbation technique is used to determine the approximate analytical solution for the Reynolds viscosity model.Taking Maclaurin’s series, (13a) takes the form(14)μ=exp⁡-Yθ=1-Yθ+OY2.Series solutions of the equation of the velocity and temperature may be obtained by using perturbation method taking ϵ as the perturbation parameter. Solutions are obtained in the following form:(15)w=w0+ϵw1+Oϵ2,θ=θ0+ϵθ1+Oϵ2.Substituting (15) into (9) and changing parameters into terms yield(16)Oϵ0:rd2w0dr2+dw0dr=Cr,(17)Oϵ1:dw1dr+rd2w1dr2-Pθ0dw0dr-Pθ0d2w0dr2-Prdw0drdθ0dr=-Mdw0dr2dw0dr+3rd2w0dr2+Cr.Also, substituting (14) and (15) into (10) and changing parameters into terms yield(18)Oϵ0:rd2θ0dr2+dθ0dr+dθ0dr+Γrdw0dr2=0,(19)θϵ1:rd2θ1dr2+dθ0dr+dθ1dr+Γrdw0dr2·Mdw0dr2-Pθ0+2Γrdw0drdw1dr+rNθ0=0.The boundary conditions for the leading order equation are(20)dw0dr0=dθ0dr0=0,w01=APdw0drr=1/D,θ01=-1Bidθ0drr=1/D,where Λ=ϵM, δ=ϵN, γ=ϵP.With the boundary conditions, it could be easily shown that (16) and (18) are given as(21)w0=C4r2+2ApD-1D2,θ0=ΓC2641-r4+4BiD3.Also, for the first-order equation, the boundary conditions are given as (22)dw1dr0=dθ1dr0=0,w11=APdw1drr=1/D,θ11=-1Bidθ1dr,where D=1+8Kn and Ap=-2Kn.Substituting the corresponding terms from the solutions in (21) into (17) and (19) and integrating the resulting equations two times successively, with the applications of the boundary conditions in (22) yields(23)w1=APPΓC36412D-16D5+2BiD4-APPΓC3192D2-APMC58D3-PΓC31152+PΓC364r24-r636+r2BiD3-14D2+136D6-1BiD5-MC332r4-1D4,θ1=MΓC496BiD5-PΓ2C4256Bi14D3-18D7+1BiD6+PΓ2C446Bi18D3-148D7+12BiD6-PΓ2C4153BiD7-MC4Γ48Bid5+ΓC2N64Bi12D-16D5+2BiD4-MΓC4576r4-1D6+PΓ2C4256r416-r864+r44BiD3-116D4+164D8-114BiD7-PΓ2C464r432-r8384+r36BiD3-132D4+1384D8-16BiD6+PΓ2C412288r8-1D8+MC4Γ288r6-1D6-ΓC2N64r24-r636+r2BiD3-14D+136D6-1BiD5.Substituting (21) and (23) into (15) the series solution at the order of expansion and then changing terms back to original parameters finally give(24)w=C4r2+2APD-1D2+γAPΓC36412D-16D5+2BiD4-γAPΓC3192D5-ΛAPC38D3+γΓC364r24-r636+r2BiD3-14D2+136D6-1BiD5-γΓC31152r6-1D6-ΛC332r6-1D6,(25)θ=ΓC2641-r4+4BiD3+ΛΓC496BiD5-γΓ2C4256Bi14D3-18D7+1BiD6+γΓ2C464Bi18D3-148D7+12BiD5-γΓ2C41536BiD7-ΛC4Γ48BiD5+∂ΓC264Bi12D-16D5+2BiD4-ΛΓC4576r6-1D6+YΓ2C4256r416-r864+r44BiD3-116D4+164D8-114BiD7-YΓ2C464r432-r8384+r36BiD3-132D4+1384D8-16BiD6+PΓ2C412288r8-1D8+ΛC4Γ288r6-1D6-∂ΓC264r24-r636+r2BiD3-14D+136D6-1BiD5.


## 2.2. Vogel’s Model of Viscosity

With the aid of Maclaurin series, Vogel’s viscosity model can be written as(26)μ=exp⁡AB-T01-ϵAθB2+Oϵ2=exp⁡AB-T01-ϵAθB2.Taking the series solution of the velocity and temperature fields yields the expansion (27)w=w+ϵw1+Oϵ2,(28)θ=ϵθ0+ϵ2θ1+Oϵ3.Therefore substituting (26) and (27) into (9) and changing parameters into term yield(29)Oϵ0:rd2w0dr2+dw0dr=C∗r,where C∗=C/exp⁡A/B-T0,(30)Oϵ0:rd2θ0dr2+dw1dr=AB2θ0rd2w0dr2+dw0dr+AB2rdθ0drdw0dr-MC∗Cdw0dr23rd2w0dr2+dω0dr,where Λ=ϵM,δ=ϵN,Γ=ϵQ.Substituting (26), (27), and (28) into (10) and changing parameters into term yield(31)Oϵ1:rd2θ0dr2+dθ0dr+Qrdw0dr2CC∗=0,where C∗=C/exp⁡A/B-T0,  Λ=ϵM,  δ=ϵN,  Γ=ϵQ, (32)θϵ2:rd2θ1dr2+dθ1dr=QrCC∗AB2θ0dw0dr2-2rQCC∗dw0drdw1dr-rQMdw0dr4-Nrθ0=0.The boundary conditions for the leading order are as follows:(33)dw0dr0=dθ0dr0=0,w01=APdw0drr=1/D,θ01=-1Bidθ0drr=1/D.Using the boundary condition, it could be easily shown that (29) and (30) are given as(34)w0=C∗4r2+2ApD-1D2,(35)θ0=QCC∗644D3Bi-r4+1D4.Also, the first-order boundary condition is as follows:(36)dw1dr0=dθ1dr0=0,w11=APdw1drr=1/D,θ11=-1Bidθ1dr,where D=1+8Kn and Ap=-2Kn.Substituting the corresponding terms from the solutions in (33) and (34) into (30) and (32) and integrating the resulting equations two times successively, with the applications of the boundary conditions in (36) yields(37)w1=APAQCC∗264B22D4Bi-16D5+12D5+APAQCC∗192B2D5+ApMC∗48CD3+AQCC∗264B2r2D3Bi-r636+r24D4-1D5Bi+136D6-14D6+AQCC∗1152B2r6-1D6+MC∗432Cr4-1D4,(38)θ1=Q2C2C∗A256B2·Bi1D6Bi-18D7+14D4+AQ2C2C∗264B2·Bi12D6Bi-148D7+18D7+AQ2C2C∗1536B2·BiD7+MCC∗4q48CBi·D5+C∗4QM96Bi·D5+QCC∗N64Bi2D4Bi-16D5+12D5+Q2CC∗2A256B2r44D3Bi-r864+r416D4-14D7Bi+164D8-116D8-AQ2C2C∗264B2r48D3Bi-r8384+r432D4-18D7Bi+138D8-132D8-AQ2C2C∗12288B2r8-1D8-MCC∗4Q288Cr6-1D6-C∗4QM576r6-1D6-QCC∗N64r2D3Bi-r636+r24D4-1D5Bi+136D6-14D6.Substituting (34) and (37), (35), and (38) into (27) and (28), the series solution at the order of expansion, and then changing terms back to original parameters finally give(39)w=C∗4r2+2APD-1D2+APΓACC∗264B22D4Bi-16D5+12D5+APΓACC∗192B2D5+APΛC∗48CD3+APΓACC∗264B2r2D3Bi-r636+r24D4-1D5Bi+136D6-14D6+AΓCC∗1152B2r6-1D6+ΛC∗432Cr4-1D4,θ=ΓCC∗644D3Bi-r4+1D4+Γ2C2C∗2A256B2·Bi1D6Bi+18D2+14D4+Γ2AC2C∗264B2·Bi12D6Bi-148D7+18D7+Γ2AC2C∗1536B2·BiD7+ΛC∗4Γ48Bi·D5+C∗4ΓΛ96Bi·D5+ΓCC∗∂64Bi2D4·Bi-16D5+12D5+Γ2C2C∗2A256B2r44D3Bi-r864+r416D4-14D7Bi+164D8-116D8Γ2AC2C∗21228B2r8-1D8-ΛC∗4288r6-1D6-ΓC∗4Λ576r6-1D6-ΓCC∗∂64r2D3Bi-r636+r24D4-1D5Bi-136D6-14D6,where(40)C∗=cexp⁡A/B-T0.


## 3. Results and Discussion

Plotting the results for the temperature distribution for both Reynolds and Vogel viscosity models of Jayeoba and Okoya [12] who studied the flow and heat transfer under no slip and no temperature jump conditions and that of the temperature distribution obtained in this work under slip and temperature jump conditions are reported graphically below at varying parameters.Figure2 shows the effect of pressure drop parameter (C) for the Reynolds viscosity model on the temperature distribution. From the result, it shows that as C becomes more negative for no slip condition, the maximum temperature of the fluid which is at the center of the pipe increases, while for the slip condition the maximum temperature increases but the fully developed temperature profile does not begin at the origin which was shown in Figure 3.Figure 2
Effect of pressure gradient parameter on the temperature distribution of Reynolds viscosity model under no slip condition whenγ=Γ=Λ=δ=1.Figure 3
Effect of varying parameterC for the temperature distribution of Reynolds viscosity model under slip condition when Λ=δ=Γ=γ=1.When heat generation term is excluded, the effect is seen asδ=0. For the Reynolds viscosity model as seen in Figures 4 and 5, as δ increases the temperature distribution of the fluid increases and the effect of the heat generation on the fluid is maximum at the center of the pipe for no slip condition, while for slip condition the effect of increasing δ causes an increase in temperature distribution but the temperature profile does not begin at the origin. Increasing values of viscous dissipation (Γ) gives a corresponding increase in temperature distribution for the Reynolds viscosity model because mechanical energy is converted into thermal energy. When the condition is no slip, Γ is maximum at the pipe center as depicted in Figure 6. However, when the flow is under the slip, temperature profile is not from the origin as shown in Figure 7. The third-grade fluid can be seen to exhibit Newtonian character for the Reynolds viscosity model when Λ=0 and also the effect when the fluid behaves in non-Newtonian manner at increasing values of Λ. It can be seen from Figure 8 that, at increasing values of Λ, the temperature distribution of the fluid decreases under the no slip condition and the maximum effect of Λ occurs at the center of the pipe, while for slip condition at increasing values of Λ the temperature distribution decreases and the fully developed temperature profile does not begin at the origin which can be seen in Figure 9. Viscosity variation parameter (γ) for the Reynolds viscosity model can be viewed from Figures 10 and 11. As viscosity variation parameter (γ) increases, the temperature distribution increases for no slip condition and the effect of γ is maximum at the center of the pipe, in Figure 10. As γ increases for the slip condition, the temperature distribution increases and the fully developed profile does not begin at the origin as seen from Figure 11.Figure 4
Effect of varying parameterδ for the temperature distribution of Reynolds viscosity model under no slip condition when Λ=-C=Γ=γ=1.Figure 5
Effect of varying parameterδ for the temperature distribution of Reynolds viscosity model under slip condition when Λ=-C=Γ=γ=1.Figure 6
Effect of varying parameterΓ for the temperature distribution of Reynolds viscosity model under no slip condition when γ=δ=Λ=-C=1.Figure 7
Effect of varying parameterΓ for the temperature distribution of Reynolds viscosity model under slip condition when Λ=-C=δ=γ=1.Figure 8
Effect of varying parameterΛ for the temperature distribution of Reynolds viscosity model under no slip condition when γ=δ=Γ=-C=1.Figure 9
Effect of varying parameterΛ for the temperature distribution of Reynolds viscosity model under slip condition when Γ=-C=δ=γ=1.Figure 10
Effect of varying parameterγ for the temperature distribution of Reynolds viscosity model under no slip condition for Λ=δ=Γ=-C=1.Figure 11
Effect of varying parameterγ for the temperature distribution of Reynolds viscosity model under slip condition when Γ=-C=δ=Λ=1.Increasing values of parameterA in Vogel’s viscosity model gives decreasing temperature distribution and the effect of Vogel’s parameter A is maximum at the center of the pipe for the no slip condition shown in Figure 12. But the temperature distribution profile does not begin at the origin for slip condition which can be seen in Figure 13. Also, for the Vogel viscosity parameter B the influence of increasing value of B shows a corresponding increase in temperature distribution, where the maximum temperature occurs at the center of the pipe. Fully developed temperature profile begins from the origin which can be glanced in Figure 14. And for slip condition, the temperature distribution profile is not from the origin as in Figure 15. As pressure drop parameter (C) becomes more negative for Vogel’s viscosity model, the temperature distribution increases with the maximum temperature at the center of the pipe and temperature profile begins at the origin for no slip condition as shown in Figure 16, while effects of pressure drop parameter on the the temperature distribution for slip condition are shown in Figure 17. The fully developed temperature profile is not from the origin. When the heat generation term is excluded, the effect is seen as δ=0 for Vogel’s viscosity model. When the condition is no slip, temperature profile begins at the origin and effect of heat generation term is maximum at the center of the pipe as shown in Figure 18. And when the condition is slip, the fully developed profile does not begin at the origin, the effect of heat generation is also maximum at the center, and the temperature distribution increases as shown in Figure 19.Figure 12
Effect of varying parameterA for the temperature distribution of Vogel’s viscosity model under no slip condition when B=-C=T0=Γ=Λ=δ=1.Figure 13
Effect of varying parameterA for the temperature distribution of Vogel’s viscosity model under slip condition when B=-C=T0=Γ=Λ=δ=1.Figure 14
Effect of varying parameterB for the temperature distribution of Vogel’s viscosity model under no slip condition when A=-C=T0=Γ=Λ=δ=1.Figure 15
Effect of varying parameterB for the temperature distribution of Vogel’s viscosity model under slip condition when A=-C=T0=Γ=Λ=δ=1.Figure 16
Effect of varying parameterC for the temperature distribution of Vogel’s viscosity model under no slip condition when A=B=T0=Γ=Λ=δ=1.Figure 17
Effect of varying parameterC for the temperature distribution of Vogel’s viscosity model under slip condition when A=B=T0=Γ=Λ=δ=1.Figure 18
Effect of varying parameterδ for the temperature distribution of Vogel’s viscosity model under no slip condition when A=B=-C=T0=Γ=Λ=1.Figure 19
Effect of varying parameterδ for the temperature distribution of Vogel’s viscosity model under slip condition when A=B=-C=T0=Γ=Λ=1.The viscous dissipation parameter (Γ) for Vogel’s viscosity model shows that, for increasing values of Γ, the temperature distribution increases. The effect of Γ is maximum at the pipe center and the temperature profile begins at the origin for no slip condition which can be seen from Figure 20, while for slip condition the effect of Γ is maximum at the center of the pipe and the fully developed temperature profile does not begin from the origin as shown in Figure 21. Increasing values of non-Newtonian parameter (Λ) for Vogel’s model gives decreasing values of temperature distribution where maximum effect occurs at the center of the pipe. The fluid behaves in a Newtonian manner when Λ=0, as can be seen from Figure 22. When condition is slip, the effect of increasing values of Λ gives a corresponding decrease in temperature distribution and the fully developed temperature profile does not begin at the origin which can be seen in Figure 23. When the initial temperature term T0 for Vogel’s viscosity model is increased, there is a corresponding increase in temperature distribution. The effect of T0 is maximum at the center of the pipe. When condition is no slip, the fully developed temperature profile begins from the origin as shown in Figure 24 and for slip condition increase in T0 gives a corresponding increase in temperature distribution with the fully developed temperature profile not beginning from the origin, as depicted in Figure 25. From the velocity distribution in Figure 26, as the Knudsen number increases from zero, it can be seen that the velocity profile reduces, with the arrow indicating the direction of increasing Knudsen number. Thus, as the Knudsen number increases, the profile of the velocity decreases.Figure 20
Effect of varying parameterΓ for the temperature distribution of Vogel’s viscosity model under no slip condition when A=B=-C=T0=δ=Λ=1.Figure 21
Effect of varying parameterΓ for the temperature distribution of Vogel’s viscosity model under slip condition when for A=B=-C=T0=δ=Λ=1.Figure 22
Effect of varying parameterΛ for the temperature distribution of Vogel’s viscosity model under no slip condition when A=B=-C=T0=δ=Γ=1.Figure 23
Effect of varying parameterΛ for the temperature distribution of Vogel’s viscosity model under slip condition when for A=B=-C=T0=δ=Γ=1.Figure 24
Effect of varying parameterT0 for the temperature distribution of Vogel’s viscosity model under no slip condition when B=-C=Λ=δ=Γ=1.Figure 25
Effect of varying parameterT0 for the temperature distribution of Vogel’s viscosity model under slip condition when A=B=-C=Λ=δ=Γ=1.Figure 26
Velocity distribution profile as the Knudsen number increases from 0 to 0.08.


## 4. Conclusion

In this work, approximate analytical solutions for the nonlinear analysis of heat transfer in a pipe flow of a third-grade fluid with temperature-dependent viscosities and heat generation under slip and no slip conditions have been developed using regular perturbation technique. The obtained analytical solutions were used to investigate the effects of slip boundary conditions on the heat transfer in the pipe flow. Also, the effects of parameters such as variable viscosity, non-Newtonian parameter, viscous dissipation, and pressure gradient at various values (which indicate different fluid flow situations) are established. The results can be used to advance the analysis and study of the behavior of third-grade fluid flow and steady state heat transfer processes such as those found in coal slurries, polymer solutions, textiles, ceramics, catalytic reactors, and oil recovery applications.



---

*Source: 8915745-2016-10-10.xml*