# A Note on Improved Homotopy Analysis Method for Solving the Jeffery-Hamel Flow

Authors: Sandile Sydney Motsa; Precious Sibanda; Gerald T. Marewo; Stanford Shateyi
Journal: Mathematical Problems in Engineering
 (2010)

---

## Abstract

This paper presents the solution of the nonlinear equation that governs the flow of a viscous, incompressible fluid between two converging-diverging rigid walls using an improved homotopy
analysis method. The results obtained by this new technique show that the improved homotopy
analysis method converges much faster than both the homotopy analysis method and the optimal
homotopy asymptotic method. This improved technique is observed to be much more accurate
than these traditional homotopy methods.


---

## Body

## 1. Introduction

The mathematical study of the flow of a viscous incompressible two-dimensional fluid in a wedge-shaped channel with a sink or source at the vertex was pioneered by Jeffery [1] and Hamel [2]. The problem has since been studied extensively by,  among others, Axford [3] who included the effects of an externally applied magnetic field and Rosenhead [4] who obtained a general solution containing elliptic functions.Instability and bifurcation are other aspects of the Jeffery-Hamel problem that have attracted widespread interest; see, for example, Akulenko and Kumakshev [5, 6]. Three-dimensional extensions to and bifurcations of the Jeffery-Hamel flow have been made by Stow et al. [7] while McAlpine and Drazin [8] presented a normal mode analysis of two-dimensional perturbations of a viscous incompressible fluid driven between inclined plane walls by a line source at the intersection of the walls. Banks et al. [9] investigated various perturbations and the linear temporal stability of such flows and found evidence of a strong interaction between the disturbances up- and down-stream if the angle between the planes exceeds a certain Reynolds-number-dependent critical value. Makinde and Mhone [10] investigated the temporal stability of MHD Jeffery-Hamel flows. They showed that an increase in the magnetic field intensity has a strong stabilizing effect on both diverging and converging channel geometries. A review of the theory of instabilities and bifurcations in channels is given by Drazin [11].As with most problems in science and engineering, the equations governing the Jeffery-Hamel problem are highly nonlinear and so generally do not have closed form analytical solutions. Nonlinear equations can, in principle, be solved by any one of a wide variety of numerical methods. However, numerical solutions rarely give intuitive insights into the effects of various parameters associated with a problem. Consequently, most recent studies of flows in diverging and converging channels have centred on the use of the Jeffery-Hamel flow equations as a testing and proving tool for the accuracy, reliability, and robustness of new techniques for solving nonlinear equations. Examples of such techniques include the summation series technique [12], the homotopy analysis method [13, 14], the decomposition method [15], the homotopy perturbation method [16], Hermite-Padé approximation [17], and the spectral-homotopy analysis method [18–21]. The study by Joneidi et al. [14] used three methods: the differential transform method (DTM), the homotopy analysis method (HAM), and the homotopy perturbation method (HPM) to solve the Jeffery-Hamel problem. The study confirmed that although both the DTM and the HPM give acceptable accuracy, the HAM is by far the superior method delivering faster convergence and better accuracy. Nonetheless, important improvements have been made to the HAM by Motsa et al. [18, 19]. The spectral modification of the HAM or SHAM proposed by Motsa et al. [18, 19] removes some of the prescriptive assumptions associated with the HAM and further accelerates the convergence rate of the method.A recent study by Esmaeilpour and Ganji [22] reported on the solution of the Jeffery-Hamel problem using an optimal homotopy asymptotic method (OHAM). This method was developed by Marinca and Herişanu [23] for the approximate solution of nonlinear problems of thin film flow of a fourth-grade fluid. The OHAM was used by these authors and others to solve the equations for the steady flow of a fourth-grade fluid in a porous medium [24, 25] and has since been applied to many other nonlinear problems including the squeezing flow problem by Idrees et al. [26]. Studies by Islam et al. [27, 28] further suggest that the OHAM is more general than both the HAM and the HPM with the latter methods being special cases of the OHAM.In this paper, we report on a new and improved method known as improved homotopy analysis method (IHAM) for solving general boundary value problems. The IHAM is an algorithm that seeks to improve the initial approximation that is later used in the HAM to solve the governing nonlinear equation resulting in significant improvement in the accuracy and convergence rate of the solutions. We seek to demonstrate the application of this method by solving the Jeffery-Hamel problem and to show its accuracy and rapid convergence by comparing the present solutions with those in the literature, including the HAM, SHAM, and the OHAM. The IHAM solutions are further compared with numerical solutions. Finally, we believe that this work will motivate further improvements to the homotopy-based semi-analytical methods.


## 2. Governing Equations

We consider the steady two-dimensional flow of a viscous incompressible fluid between two nonparallel rigid planes with angle2α. Assuming radial flow with velocity v=u(r,θ), the motion of such a fluid is described by the equations (see [14, 19, 22])(2.1)ρr∂∂r(ru)=0,(2.2)u∂u∂r=-1ρ∂p∂r+ν[∂2u∂r2+1r∂u∂r+1r2∂2u∂θ2-ur2],(2.3)-1ρr∂p∂θ+2νr2∂u∂θ=0,
where ρ is the fluid density, p the pressure, and ν the kinematic viscosity. The continuity equation implies that (2.4)h(θ)=ru(r,θ),
and by using the nondimensional parameters (2.5)f(y)=f(θ)fmax⁡,y=θα,
Equations (2.2)-(2.3) reduce to(2.6)f′′′+2αReff′+4α2f′=0,
where (2.7)Re=αfmax⁡ν=Umax⁡rαν
is the Reynolds number and Umax is the maximum velocity at the centre of the channel. The appropriate boundary conditions are(2.8)f(0)=1,f′(0)=0,f(1)=0.


## 3. Improved Homotopy Analysis Method (IHAM) Solution

In this section, we describe the use of the improved homotopy analysis method (IHAM) in the governing equation (2.6). To apply the IHAM, we assume that the solution f(y) can be expanded as(3.1)f(y)=fi(y)+∑n=0i-1fn(y),i=1,2,3,…,
where fi are unknown functions whose solutions are obtained using the HAM approach at the ith iteration and fn, (1≤n≤i-1) are known from previous iterations. The algorithm starts with the initial approximation f0(y) which is chosen to satisfy the boundary conditions (2.8). An appropriate initial guess is(3.2)f0(y)=1-y2.
Substituting (3.1) in the governing equation (2.6)-(2.8) gives(3.3)fi′′′+a1,i-1fi′+a2,i-1fi+2αRefifi′=ri-1,
subject to the boundary conditions(3.4)fi(0)=0,fi′(0)=0,fi(1)=0,
where the coefficient parameters ak,i-1, (k=1,2) and ri-1 are defined as(3.5)a1,i-1=2αRe∑n=0i-1fn+4α2,a2,i-1=2αRe∑n=0i-1fn′,ri-1=-[∑n=0i-1fn′′′+2αRe∑n=0i-1fn′∑n=0i-1fn+4α2∑n=0i-1fn′].
Starting from the initial approximation (3.2), the subsequent solutions fi (i≥1) are obtained by recursively solving equation (3.3) using the HAM approach [29, 30]. To find the HAM solutions of (3.3), we begin by rewriting (3.3) as(3.6)N[fi(y)]=ri-1,
where𝒩is a nonlinear operator defined by
(3.7)N[fi(y)]=fi′′′+a1,i-1fi′+a2,i-1fi+2αRefifi′.
Let fi,0(y) denote the initial guess for the unknown function fi(y), and let ℏ≠0 be an auxiliary parameter. Using an embedding parameter q∈[0,1] we construct a homotopy (3.8)H[Fi(y;q);fi,0(y),ℏ,q]=(1-q)L[Fi(y;q)-fi,0(y)]-qℏ{N[Fi(y;q)]-ri-1},
whereLis an auxiliary linear operator with the property that ℒ[0]=0 and Fi(y;q) is an unknown function. Upon equating ℋ to 0, we obtain the zero-order deformation equation
(3.9)(1-q)L[Fi(y;q)-fi,0(y)]=qℏ{N[Fi(y;q)]-ri-1}.
When, q=0 (3.9) becomes
(3.10)L[Fi(y;0)-fi,0(y)]=0.
This equation holds provided
(3.11)Fi(y;0)=fi,0(y)
as ℒ[0]=0. Whenq=1, (3.9) is simplified to
(3.12)N[Fi(y;1)]=ri-1
asℏ≠0. Equation (3.12) is the same as equation (3.6) provided
(3.13)Fi(y;1)=fi(y).
It follows from (3.11) and (3.13) that as q increases from 0 to 1, the unknown function F(y;q) varies continuously from initial guess fi,0(y) to exact solution fi(y) of (3.6).Differentiating (3.9) m times with respect to q and then setting q=0 and finally dividing the resulting equations by m! yields the mth-order deformation equations(3.14)L[fi,m(η)-χmfi,m-1]=ℏ(fi,m-1′′′+a1,i-1fi,m-1′+a2,i-1fi,m-1+2αRe∑j=0m-1fi,jfi,m-1-j′′-(1-χm)ri-1),
subject to the boundary conditions(3.15)fi,m(0)=fi,m′(0)=fi,m(1)=0,
where(3.16)χm={0,m≤1,1,m>1..
If we define the linear operator ℒ by(3.17)L[Fi(y;q)]=∂3Fi∂y3,
then the initial approximation fi,0 that is used in the higher-order equations (3.14) is obtained by solving the differential equation(3.18)L[fi,0]=fi,0′′′=ri-1,
subject to the boundary conditions(3.19)fi,0(0)=fi,0′(0)=fi,0(1)=0.
Thus, starting from the initial approximation, which is obtained-from (3.18), higher order approximations fi,m(y) for m≥1 can be obtained through the recursive formula (3.14). We note that (3.14) forms a set of linear ordinary differential equations and can be easily solved analytically, especially by means of symbolic computation software such as Maple, Mathematica, Matlab, and others.ExpandingFi(y;q) in Taylor series about q=0 gives(3.20)Fi(y;q)=Fi(y;0)+∑k=1+∞qkk!∂kFi(y;q)∂yk|q=0.
If we set q=1 in (3.20) and define
(3.21)fi,k(y):=1k!∂kFi(y;q)∂yk|q=0
for each k=1,2,…, then making use of (3.11) and (3.13) transforms (3.20) to
(3.22)fi(y)=fi,0(y)+∑k=1+∞fi,k(y).
Upon truncating the infinite series in (3.22), the solutions for fi are then generated using the solutions for fi,m as follows:(3.23)fi=fi,0+fi,1+fi,2+fi,3+fi,4+⋯+fi,m.
The [i,m] approximate solution for f(y) is then obtained by substituting fi (obtained from (3.23)) in (3.1).


## 4. Results and Discussion

The accuracy and reliability of the IHAM is determined by comparing the current results with the optimal homotopy analysis method (OHAM) of Esmaeilpour and Ganji [22], the homotopy analysis results (HAM) of Joneidi et al. [14], and the numerical results obtained using the MATLAB bvp4c routine which is a boundary value solver based on the adaptive Lobatto quadrature scheme [31, 32]. In Table 1 we show a comparison of the current velocity results at various dimensionless angles y for different orders [m,n] of the IHAM. The Reynolds number is fixed at Re=50, and we consider a diverging channel with angle α=5∘. Convergence of the current results to the numerical results to six decimal places is achieved at order [2,2] while convergence to ten decimal places is achieved at order [2,3]. On the other hand, convergence of the OHAM is evidently slow, and up to the twentieth order of approximation, the method mostly fails to converge to the numerical results.Table 1
Comparison between the numerical results and the order[m,n] IHAM approximate results (using ℏ=-1) for f(y) against the OHAM results reported in [22] when Re=50, α=5.
yOHAM [22][2,2][2,3][2,4]Numerical01.000000001.000000001.000000001.000000001.000000000.10.982518090.982431210.982431240.982431240.982431240.20.931565890.931225830.931225960.931225960.931225960.30.851381550.850610330.850610610.850610620.850610620.40.748260400.746790350.746790800.746790800.746790800.50.629538650.626947610.626948170.626948170.626948170.60.502428940.498233890.498234450.498234450.498234450.70.372933830.366965890.366966340.366966340.366966340.80.245081980.238123470.238123750.238123750.238123750.90.120715620.115151810.115151930.115151930.1151519310.000000000.000000000.000000000.000000000.00000000The slow convergence of the OHAM is confirmed in Table2 where the absolute errors of the two methods in relation to the numerical solution are given. At any angle y, the absolute error using the OHAM is much larger than that obtained using the IHAM. The largest absolute error at any angle y using the IHAM is 5.6×10-7 while, by comparison, the largest absolute error when using the OHAM is 6.958×10-3.Table 2
Comparison between the errors of the order[m,n] IHAM results against the OHAM results reported in [22] for Re=50 and α=5.
OHAM[2,2][2,3][2,4]0.000000000.000000000.000000000.000000000.000086850.000000030.000000000.000000000.000339930.000000130.000000000.000000000.000770930.000000290.000000010.000000000.001469600.000000450.000000000.000000000.002590480.000000560.000000000.000000000.004194490.000000560.000000000.000000000.005967490.000000450.000000000.000000000.006958230.000000280.000000000.000000000.005563690.000000120.000000000.000000000.000000000.000000000.000000000.00000000Table3 gives a comparison of the convergence rate of the current method with the HAM [14]. Convergence of the IHAM to the numerical results is found to be rapid, with agreement to ten decimal places being archived at order [3,3] of the IHAM algorithm. The sixteenth-order HAM on the other hand fails to achieve the accuracy of both the IHAM and the numerical results. The poor performance of the HAM is further confirmed in Table 4 where the largest absolute error incurred by the IHAM at order [3,2] is 7.2×10-9 while the largest absolute error incurred by using the sixteenth-order HAM at any angle y is approximately a hundred times larger at 4.7×10-7.Table 3
Comparison between the numerical results and the order[m,n] IHAM approximate results (using ℏ=-1) for f(y) against the HAM results reported in [14] when Re=110, α=3.
yHAM [14][3,2][3,3][3,4]Numerical01.00000000001.00000000001.00000000001.00000000001.00000000000.10.97923570620.97923570590.97923570650.97923570650.97923570650.20.91926588420.91926588390.91926588560.91926588560.91926588560.30.82653361020.82653361050.82653361230.82653361230.82653361230.40.71022118380.71022118310.71022118330.71022118320.71022118320.50.58049947000.58049946170.58049945880.58049945880.58049945880.60.44693509410.44693507280.44693506710.44693506700.44693506700.70.31740845450.31740843480.31740842760.31740842760.31740842760.80.19764106610.19764110120.19764109450.19764109450.19764109450.90.09123022880.09123042520.09123042110.09123042110.091230421110.00000047000.00000000000.00000000000.00000000000.0000000000Table 4
Comparison between the errors of the order[m,n] IHAM results against the HAM results reported in [14] for Re=110 and α=3.
yHAM [14][3,2][3,3][3,4]00.00000000000.00000000000.00000000000.00000000000.10.00000000030.00000000060.00000000000.00000000000.20.00000000140.00000000170.00000000000.00000000000.30.00000000210.00000000180.00000000000.00000000000.40.00000000060.00000000010.00000000010.00000000000.50.00000001120.00000000290.00000000000.00000000000.60.00000002710.00000000580.00000000010.00000000000.70.00000002690.00000000720.00000000000.00000000000.80.00000002840.00000000670.00000000000.00000000000.90.00000019230.00000000410.00000000000.000000000010.00000047000.00000000000.00000000000.0000000000Figure1 shows the effect of the Reynolds number on the fluid velocity for different values of α. It can be seen from the figure that the fluid velocity increases with Reynolds numbers in the case of convergent channels (α<0) but decreases with Re in the case of divergent channels (α>0). This observation is consistent with the observations made in [13, 14, 19].Velocity profilef(y) for different values of Re when α=5 and α=-5, respectively.
(a)(b)Figure2 shows the effect of the Reynolds number on f′′(0), which is related to wall shear stress, for different values of α. It can be seen from the figure that f′′(0) decreases monotonically for α>0. In the case of α<0, f′′(0) increases for a range of Re until it reaches a peak then decreases.Variation off′′(0) against the Reynolds number when α=5 and α=-5, respectively.
(a)(b)


## 5. Conclusion

In this brief note we have proposed an improved homotopy analysis method (IHAM) for the solution of general nonlinear differential equations. We have compared the performance of the new algorithm against the optimal homotopy analysis method (OHAM), the standard homotopy analysis method (HAM), and numerical approximations by solving the Jeffery-Hamel problem for Newtonian flow in converging/diverging channels. Numerical computations show that the IHAM is accurate and converges to the numerical approximations at lower orders compared to both the HAM and the OHAM. However, at this juncture, we cannot say with certainty that this method is better than these other existing methods. We recommend that we need to use this improved homotopy analysis method (IHAM) to solve more nonlinear differential equations.



---

*Source: 359297-2010-12-30.xml*