This work is concerned with the series solutions for the flow of third-grade non-Newtonian fluid with variable viscosity. Due to the nonlinear, coupled, and highly complicated nature of partial differential equations, finding an analytical solution is not an easy task. The homotopy analysis method (HAM) is employed for the presentation of series solutions. The HAM is accepted as an elegant tool for effective solutions for complicated nonlinear problems. The solutions of (Hayat et al., 2007) are developed, and their convergence has been discussed explicitly for two different models, namely, constant and variable viscosity. An error analysis is also described. In addition, the obtained results are illustrated graphically to depict the convergence region. The physical features of the pertinent parameters are presented in the form of numerical tables.
1. Introduction
During the last few years, there has been substantial progress in the steady and unsteady flows of non-Newtonian fluids. A huge amount of literature is now available on the topic (see some studies [1–6]). All real fluids are diverse in nature. Hence in view of rheological characteristics, all non-Newtonian fluids cannot be explained by employing one constitutive equation. This is the striking difference between viscous and the non-Newtonian fluids. The rheological parameters appearing in the constitutive equations lead to a higher-order and complicated governing equations than the Navier-Stokes equations. The simplest subclass of differential-type fluids is called the second grade. In steady flow such fluids can predict the normal stress and does not show shear thinning and shear thickening behaviors. The third-grade fluid puts forward the explanation of shear thinning and shear thickening properties. Therefore, the present paper aims to study the pipe flow of a third-grade fluid. Some progress on the topic is mentioned in the studies [7, 8] and many references therein. In all these studies, variable viscosity is used. Massoudi and Christie [9] numerically examined the pipe flow of a third-grade fluid when viscosity depends upon temperature. Hayat et al. [10] presented the homotopy solution of the problem considered in [10] up to second-order deformation.
In this paper, the motivation comes from a desire to understand the convergence of the problem discussed in [10]. The relevant equations for flow and temperature have been solved analytically by using homotopy analysis method [11–15]. Here the convergence of the obtained solutions is explicitly shown,and that was not previously given in [10].
2. Problem
From [10], we have the equations (2.1) to (3.4) in nondimensional and nonlinear coupled partial differential equations of the form 1rddr(rμ(r)(dvdr))+Λrddr(r(dvdr)3)=c,d2θdr2+1r(dθdr)+Γ(dvdr)2(μ(r)+Λ(dvdr)2)=0,
subject to boundary conditionsv(1)=θ(1)=0,dv(0)dr=dθ(0)dr=0,v(1)=θ(1)=0,dvdr(0)=dθdr(0)=0.
3. Solution of the Problem
Our interest is to carry out the analysis for the homotopy solutions for two cases of viscosity, namely, constant and space-dependent viscous dissipation.
Case I.
For constant viscosity model, we choose
μ=1.
For HAM solution, we select
v0(r)=c4(r2-1),θ0=c2Γ(1-r4)64,
as initial approximations of v and θ, respectively, which satisfy the linear operator and corresponding boundary conditions. We use the method of higher-order differential mapping [16] to choose the linear operator ℒ which is defined by
L1=d2dr2+1rddr,
such that
L1(C1+C2lnr)=0,
where C1 and C2 are the arbitrary constants.
If the convergence parameter is ħ and 0≤p≤1 is an embedding parameter, then the zeroth-order problems become
(1-p)L1[v*(r,p)-v0(r)]=pħN1[v*(r,p),θ*(r,p)],(1-p)L1[θ*(r,p)-θ0(r)]=pħN2[v*(r,p),θ*(r,p)],v*(1,p)=θ*(1,p)=0,∂v*(r,p)∂r|r=0=∂θ*(r,p)∂r|r=0=0,
where the nonlinear parameters 𝒩1 and 𝒩2 are defined by
N1[v*(r,p),θ*(r,p)]=1rdv*dr+d2v*dr2+Λr(dv*dr)3+3Λ(dv*dr)2d2v*dr2-c,N2[v*(r,p),θ*(r,p)]=1rdθ*dr+d2θ*dr2+Γ(dv*dr)2+ΓΛ(dv*dr)4.
For p=0 and p=1, we have
v*(r,0)=v0(r),θ*(r,0)=θ0(r),v*(r,1)=v(r),θ*(r,1)=θ(r).
When p increases from 0 to 1, v*(r,p),θ*(r,p) vary from v0(r),θ0(r) to v(r),θ(r), respectively. By Taylor’s theorem and (3.7), one can get
v*(r,p)=v0(r)+∑m=1∞vm(r)pm,θ*(r,p)=θ0(r)+∑m=1∞θm(r)pm,
where
vm(r)=1m!∂mv*(r,p)∂pm|p=0,θm(r)=1m!∂mθ*(r,p)∂pm|p=0.
The convergence of the series (3.8) depends upon ħ. We choose ħ in such a way that the series (3.8) is convergent at p=1; then, due to (3.7), we get
v(r)=v0(r)+∑m=1∞vm(r),θ(r)=θ0(r)+∑m=1∞θm(r).
The mth-order deformation problems are
L1[vm(r)-χmvm-1(r)]=ħR1m(r),L1[θm(r)-χmθm-1(r)]=ħR2m(r),vm(1)=θm(1)=0,vm′(0)=θm′(0)=0,
where the recurrence formulae ℜ1 and ℜ2 are given by
R1m(r)=1rdvm-1dr+d2vm-1dr2+Λr∑k=0m-1∑i=0k(dvm-1-kdr)dvk-idrdvidr+3Λ∑k=0m-1∑i=0k(dvm-1-kdr)dvk-1drd2vidr2-(1-χm)c,R2m(r)=1rdθm-1dr+d2θm-1dr2+Γ∑k=0m-1(dvm-1-kdr)dvkdr+ΛΓ∑k=0m-1∑j=0k∑i=0j(dvm-1-kdr)dvk-jdrdvj-idrdvidr
in which
χm={0,m≤1,1,m>1.
For constant viscosity, the velocity and temperature expressions up to second-order deformation are
v(r)=c4(r2-1)+hc3Λ(2h+3)(r4-1)16+h2c5Λ2(r6-1)32,θ(r)=[M1(r4-1)+M2(r6-1)+M3(r8-1)+M4(r10-1)+M5(r12-1)+M6(r14-1)+M7(r16-1)+M8(r18-1)+M9(r20-1)+M10(r22-1)].
Case II.
For space-dependent viscosity, we take
μ=r.
For HAM solution, we select
v0(r)=c6(r2-1),θ0=c4ħ4Γ(1-r2)64.
As the initial approximation of v and θ. We select
L2=d2dr2+2rddr,
such that
L2(C3+C4r)=0,
where C3 and C4 are arbitrary constants. The zeroth- and mth-order deformation problems are
(1-p)L2[v*(r,p)-v0(r)]=pħN3[v*(r,p),θ*(r,p)],(1-p)L2[θ*(r,p)-θ0(r)]=pħN4[v*(r,p),θ*(r,p)],v*(1,p)=θ*(1,p)=0,∂v*(r,p)∂r|r=0=∂θ*(r,p)∂r|r=0=0,L2[vm(r)-χmvm-1(r)]=ħR3m(r),L2[θm(r)-χmθm-1(r)]=ħR4m(r),vm(1)=θm(1)=0,vm′(0)=θm′(0)=0,
where
N3[v*(r,p),θ*(r,p)]=2rdv*dr+d2v*dr2+Λr2(dv*dr)3+3Λr(dv*dr)2d2v*dr2-cr,N4[v*(r,p),θ*(r,p)]=1rdθ*dr+d2θ*dr2+Γ(dv*dr)2+ΓΛ(dv*dr)4+Γr(dv*dr)2,R3m(r)=2rdvm-1dr+r2d2vm-1dr2+Λ∑k=0m-1∑i=0k(dvm-1-kdr)dvk-idrdvidr+3Λr∑k=0m-1∑i=0k(dvm-1-kdr)dvk-idrd2vidr2-(1-χm)cr,R4m(r)=1rdθm-1dr+d2θm-1dr2+Γr∑k=0m-1(dvm-1-kdr)dvkdr+ΛΓ∑k=0m-1∑j=0k∑i=0j(dvm-1-kdr)dvk-jdrdvj-idrdvidr.
For variable viscosity, the velocity and temperature expressions up to second-order deformation are
v(r)=hc2(r-1)+c(2h+3)(r2-1)18+c3hΛ(r3-1)81,θ(r)=[M11(r2-1)+M12(r3-1)+M13(r4-1)+M14(r5-1)+M15(r6-1)],
where the constant coefficients M1–M15 can be easily obtained through the routine calculation.
mth-order solutions
In both cases, for p=0 and p=1, we have
v*(r;0)=v0(r),θ*(r;0)=θ0(y),v*(r;1)=v(r),θ*(r;1)=θ(r).
When p increases from 0 to 1, v*(r,p), θ*(r,p)ϕ*(r,p) varies from v0(r),θ0(r)ϕ0(r) to v(r),θ(r) and ϕ(r), respectively. By Taylor’s theorem and (3.24) the general solutions can be written as
v*(r,p)=v0(r)+∑m=1∞vm(r)pm,θ*(r,p)=θ0(r)+∑m=1∞θm(r)pm,
where
vm(r)=1m!∂mv*(r,p)∂pm|p=0,θm(r)=1m!∂mθ*(r,p)∂pm|p=0.
The convergence of (3.25) depends upon ħ; therefore, we choose ħ in such a way that it should be convergent at p=1. In view of (3.24), finally the general form of mth-order solutions is
v(r)=v0(r)+∑m=1∞vm(r),θ(r)=θ0(r)+∑m=1∞θm(r).
4. Discussion
It is noticed that the explicit, analytical expressions (3.11), (19), (3.19), and (3.20) contain the auxiliary parameter ħ. As pointed out by Liao [17], the convergence region and rate of approximations given by the HAM are strongly dependent upon ħ. Figures 1 and 2 show the ħ-curves of velocity and temperature profiles, respectively, just to find the range of ħ for the case of constant viscosity. The range for admissible values of ħ for velocity is -2.4≤ħ≤0.4 and for temperature is -2.2≤ħ≤0.5. Figures 4 and 5 represent the ħ-curves for variable viscosity. The admissible ranges for both velocity and temperature profiles are -3≤ħ≤0.4 and -2.8≤ħ≤0.8, respectively. In Figures 3 and 6, the graphs of residual error are plotted for constant and variable viscosity, respectively. The error of norm 2 of two successive approximations over [0,1] with HAM by 10th-order approximations is calculated byE2=111∑i=010(v10(i10))2=f.(say)
It is seen that the error is minimum at ħ=-0.01. These values of ħ also lie in the admissible range of ħ.
ħ-curve for velocity in case of constant viscosity at 10th-order approximation.
ħ-curve for temperature in case of constant viscosity at 10th-order approximation.
Residual error curve for constant viscosity.
ħ-curve for velocity in case of variable viscosity at 10th-order approximation.
ħ-curve for temperature in case of variable viscosity at 10th-order approximation.
Residual error curve for variable viscosity.
We use the widely applied symbolic computation software MATHEMATICA to see the effects of sundry parameters by Tables 1, 2, and 3.
Illustrating the variation of the velocity and temperature with c.
h
Λ
c
V
θ
−0.01
1
−1
1.673
0.006
−2
3.191
0.068
−3
4.4331
0.270
−4
5.339
0.661
−5
5.924
1.205
Illustrating the variation of the velocity and temperature with Λ.
h
c
Γ
Λ
V
θ
−0.01
−1
1
0
1.700
0.243
5
1.571
2.002
10
1.455
3.209
15
1.353
4.011
20
1.263
4.520
Illustrating the variation of temperature with Γ.
h
c
Λ
Γ
V
θ
−0.01
1
1
0
0
0
5
0.075
3.242
10
0.158
6.484
15
0.249
9.726
20
0.351
12.969
5. Conclusion
In this paper, the convergence of series solution for constant and variable viscosity in a third-grade fluid is presented. The steady pipe flow is considered. Convergence values and residual error are also examined in Figures 1 to 6. To see the effects of emerging parameters for constant and variable viscosity, Tables 1 to 3 have been displayed. In Tables 1 and 2, it is found that the velocity and temperature increase with the decrease in pressure gradient and third-grade parameter, respectively, whereas Table 3 explains the variation of viscous dissipation parameter on velocity and temperature distributions. Here, it is revealed that the velocity and temperature decrease by increasing the viscous dissipation. It is observed that the results and figures [10] for important parameters c,Λ and Γ are correct and remain unchanged.
Acknowledgments
R. Ellahi thanks the United State Education Foundation Pakistan and CIES USA for honoring him by the Fulbright Scholar Award for the year 2011-2012. R. Ellahi is also grateful to the Higher Education Commission and PCST of Pakistan to award him the awards of NRPU and Productive Scientist, respectively.
TanW.MasuokaT.Stability analysis of a Maxwell fluid in a porous medium heated from below200736034544602-s2.0-3375073748410.1016/j.physleta.2006.08.054HayatT.MahomedF. M.Note on an exact solution for the pipe flow of a third-grade fluid20071901–42332362-s2.0-3424736894510.1007/s00707-006-0428-2ZBL1117.76005TanW.MasuokaT.Stokes' first problem for a second grade fluid in a porous half-space with heated boundary20054045155222-s2.0-954424674910.1016/j.ijnonlinmec.2004.07.016MalikM. Y.drmymalik@hotmail.comHussainA.NadeemS.Flow of a Jeffery-six constant fluid between coaxial cylinders with heat transfer analysis201156234535110.1088/0253-6102/56/2/27MalikM. Y.HussainA.NadeemS.HayatT.Flow of a third grade fluid between coaxial cylinders with variable viscosity2009649-105885962-s2.0-73249142874HameedM.NadeemS.Unsteady MHD flow of a non-Newtonian fluid on a porous plate20073251724733227356010.1016/j.jmaa.2006.02.002ZBL1102.76077EllahiR.2009Saarbrücken, GermanyVDMEllahiR.AfzalS.Effects of variable viscosity in a third grade fluid with porous medium: an analytic solution200914520562072247446510.1016/j.cnsns.2008.05.006ZBL1221.76015MassoudiM.ChristieI.Effects of variable viscosity and viscous dissipation on the flow of a third grade fluid in a pipe19953056876992-s2.0-000252328210.1016/0020-7462(95)00031-IZBL0865.76005HayatT.EllahiR.AsgharS.The influence of variable viscosity and viscous dissipation on the non-Newtonian flow: an analytical solution200712330031310.1016/j.cnsns.2005.04.0082276973ZBL1102.76003LiaoS. J.1992Shanghai, ChinaShanghai Jiao Tong UniversityLiaoS.20042Boca Raton, Fla, USAChapman & Hall2058313AbbasbandyS.The application of homotopy analysis method to nonlinear equations arising in heat transfer20063601109113228811810.1016/j.physleta.2006.07.065EllahiR.rahmatellahi@yahoo.comEffects of the slip boundary condition on non-Newtonian flows in a channel20091441377138410.1016/j.cnsns.2008.04.002HayatT.EllahiR.ArielP. D.AsgharS.Homotopy solution for the channel flow of a third grade fluid2006451-25564225116610.1007/s11071-005-9015-7ZBL1100.76005Van GorderR. A.VajraveluK.On the selection of auxiliary functions, operators, and convergence control parameters in the application of the homotopy analysis method to nonlinear differential equations: a general approach2009141240784089253756910.1016/j.cnsns.2009.03.008ZBL1221.65208LiaoS.-J.An analytic approximate technique for free oscillations of positively damped systems with algebraically decaying amplitude200338811731183195518210.1016/S0020-7462(02)00062-8