# A Novel Numerical Technique for Two-Dimensional Laminar Flow between Two Moving Porous Walls

Authors: Zodwa G. Makukula; Precious Sibanda; Sandile S. Motsa
Journal: Mathematical Problems in Engineering
 (2010)

---

## Abstract

We investigate the steady two-dimensional flow of a viscous incompressible fluid in a rectangular domain that is bounded by two permeable surfaces. The governing fourth-order nonlinear differential equation is solved by applying the spectral-homotopy analysis method  and a novel successive linearisation method. Semianalytical results are obtained and the convergence rate of the solution series was compared with numerical approximations and with earlier results where the homotopy analysis and homotopy perturbation methods were used. We show that both the spectral-homotopy analysis method and successive linearisation method are computationally efficient and accurate in finding
solutions of nonlinear boundary value problems.


---

## Body

## 1. Introduction

Laminar viscous flow in tubes that allow seepage across contracting or expanding permeable walls is encountered in the transport of biological fluids such as blood and filtration in kidneys and lungs. Such flows have many other practical applications such as in binary gas diffusion, chromatography, ion exchange, and ground water movement [1–6]. In addition, flow in channels with permeable walls provides a good starting point for the study of flow in multichannel filtration systems such as the wall flow monolith filter used to reduce emissions from diesel engines introduced by Oxarango et al. in [7]. Consequently, in the past four decades a considerable amount of research effort has been expended in the study of laminar flows in rectangular domains that are bounded by permeable walls [8–15].The equations governing such flows are generally nonlinear and in the past asymptotic techniques, and numerical methods have been used to analyze such flows and to solve the equations; for example, in the pioneering study by Berman [8] asymptotic methods were used to solve the steady flow problem for small suction. In the study by Uchida and Aoki [16], numerical methods were used to solve the governing nonlinear equations and to explain the flow characteristics. Majdalani and Roh [4] and Majdalani [3] studied the oscillatory channel flow with wall injection, and the resulting governing equations were solved using asymptotic formulations (WKB and multiple-scale techniques). The multiple-scale solution was found to be advantageous over the others in that its leading-order term is simpler and more accurate than other formulations, and it displayed clearly the relationship between the physical parameters that control the final motion. It also provided means of quantifying important flow features such as corresponding vortical wave amplitude, rotational depth of penetration, and near wall velocity overshoot to mention a few. Jankowski and Majdalani [12] used the same approach and drew similar conclusions about the multiple-scale solution for oscillatory channel flow with arbitrary suction.An analytical solution by means of the Liouville-Green transformation was developed for laminar flow in a porous channel with large wall suction and a weakly oscillatory pressure by Jankowski and Majdalani [13]. The scope of the problem had many limitations, for example, the study did not consider variations in thermostatic properties and the oscillatory pressure amplitude was taken to be small in comparison with the stagnation pressure. Zhou and Majdalani [17] investigated the mean-flow for slab rocket motors with regressing walls. The transformed governing equation was solved numerically, using finite differences, and asymptotically, using variation of parameters and small parameter perturbations in the blowing Reynolds number. The results from the two methods were compared for different Reynolds numbers Re and the wall regression rate α, and it was observed that accuracy of the analytical solution deteriorates for small Re and large α. A good agreement between the solutions was observed for large values of Re. A similar analysis was done by Majdalani and Zhou [6] for moderate-to-large injection and suction driven channel flows with expanding or contracting walls.In recent years, the use of nonperturbation techniques such as the Adomian decomposition method [18, 19]. He's homotopy perturbation method [20, 21], and the homotopy analysis method [22, 23] has been increasingly preferred to solve nonlinear differential equations that arise in science and engineering. Dinarvand et al. [2] solved Berman’s model of two-dimensional viscous flow in porous channels with wall suction or injection using both the HAM and the homotopy perturbation method (HPM). They concluded that the HPM solution is not valid for large Reynolds numbers, a weakness earlier observed in the case of other perturbation techniques. Using the homotopy analysis method, Xu et al. [24] developed highly accurate series approximations for two-dimensional viscous flow between two moving porous walls and obtained multiple solutions associated with this problem. The multiple solutions associated with this problem were also reported by Zarturska et al. [25]. Although the homotopy analysis method is a reliable and efficient semi-analytical technique, it however suffers from a number of limiting assumptions such as the requirements that the solution ought to conform to the so-called rule of solution expression and the rule of coefficient ergodicity. A modification of the homotopy analysis method, see Motsa et al. [26, 27], seeks to produce a more efficient method while also addressing the limitations of the HAM. In this paper, we use the spectral homotopy analysis method to solve the nonlinear differential equation that governs the flow of a viscous incompressible fluid in a rectangular domain bounded by two permeable walls. The problem is also solved using a new and highly efficient technique, the successive linearisation method (see [28, 29]) so as to independently corroborate and validate the SHAM results. The results are also compared with numerical approximations and the recent results reported in Xu et al. [24] and Dinarvand and Rashidi [30].


## 2. Governing Equations

Consider two-dimensional laminar, isothermal, and incompressible viscous fluid flow in a rectangular domain bounded by two permeable surfaces that enable the fluid to enter or exit during successive expansions or contractions. The walls are placed at a separation2a and contract or expand uniformly at a time-dependent rate ȧ(t). The governing Navier-Stokes equations are given in Majdalani et al. [31] as
(2.1)∂û∂x̂+∂v̂∂ŷ=0,(2.2)∂û∂t+û∂û∂x̂+v̂∂û∂ŷ=-1ρ∂p̂∂x̂+ν∇2û,(2.3)∂v̂∂t+û∂v̂∂x̂+v̂∂v̂∂ŷ=-1ρ∂p̂∂ŷ+ν∇2v̂,
where û and v̂ are the velocity components in the x̂ and ŷ directions, respectively, p̂, ρ, ν and t are the dimensional pressure, density, kinematic viscosity, and time, respectively. Assuming that inflow or outflow velocity is vw, then the boundary conditions are(2.4)û(x̂,a)=0,v̂(a)=-vw=-ȧ/c,∂û∂ŷ(x̂,0)=0,v̂(0)=0,û(0,ŷ)=0,
where c(=ȧ/vw) is the injection or suction coefficient. Introducing the stream function ψ̂=νx̂F̂(y,t)/a and the transformations (2.5)ψ=ψ̂aȧ,u=ûȧ,v=v̂ȧ,x=x̂a,y=ŷa,F=F̂Re,
Majdalani et al. [31] and Dinarvand and Rashidi [30] showed that (2.1)–(2.3) reduce to the normalized nonlinear differential equation(2.6)FIV+α(yF′′′+3F′′)+Re(FF′′′-F′F′′)=0,
subject to the boundary conditions(2.7)F=0,F′′=0,aty=0,(2.8)F=1,F′=0,aty=1,
where α(t)=ȧa/ν is the nondimensional wall dilation rate defined to be positive for expansion and negative for contraction, and Re=avw/ν is the filtration Reynolds number defined positive for injection and negative for suction through the walls. Equation (2.6) is strongly nonlinear and not easy to solve analytically, and most researchers have studied the classic Berman formula [8]; that is, when α=0. In this paper, we seek to solve (2.6) subject to the boundary conditions (2.7) and (2.8) using a novel spectral modification of the homotopy analysis method and the successive linearisation method. By comparison with the numerical approximations and previously obtained results, we show that these new techniques are accurate and more efficient than the standard homotopy analysis method.


## 3. Spectral Homotopy Analysis Method Solution

In applying the spectral-homotopy analysis method, it is convenient to first transform the domain of the problem from [0,1] to [-1,1] and make the governing boundary conditions homogeneous by using the transformations(3.1)y=ξ+12,U(ξ)=F(y)-F0(y),F0(y)=32y-12y3.
Substituting (3.1) in the governing equation and boundary conditions (2.6)–(2.8) gives(3.2)16U1V+8a1U′′′+4a2U′′+2a3U′-3ReU+8Re(UU′′′-U′U′′)=ϕ(y),
subject to(3.3)U=0,U′′=0,ξ=-1,U=0,U′=0,ξ=1,
where the primes denote differentiation with respect to ξ and(3.4)a1=αy+Re(32y-12y3),a2=3α-32Re(1-y2),a3=3yRe,ϕ(y)=12αy+3Rey3.The initial approximation is taken to be the solution of the nonhomogeneous linear part of the governing equations (3.2) given by(3.5)16U01V+8a1U0′′′+4a2U0′′+2a3U0′-3ReU0=ϕ(y),
subject to(3.6)U0=0,U0′′=0,ξ=-1,U0=0,U0′=0,ξ=1.
We use the Chebyshev pseudospectral method to solve (3.5)–(3.6). The unknown function U0(ξ) is approximated as a truncated series of Chebyshev polynomials of the form(3.7)U0(ξ)≈U0N(ξj)=∑k=0NÛkTk(ξj),j=0,1,…,N,
where Tk is the kth Chebyshev polynomial, Ûk, are coefficients and ξ0,ξ1,…,ξN are Gauss-Lobatto collocation points (see [32]) defined by(3.8)ξj=cosπjN,j=0,1,…,N.
Derivatives of the functions U0(ξ) at the collocation points are represented as(3.9)drU0dξr=∑k=0N𝒟kjrU0(ξj),
where r is the order of differentiation and 𝒟 is the Chebyshev spectral differentiation matrix ([32, 33]). Substituting (3.7)–(3.9) in (3.5)–(3.6) yields(3.10)AU0=Φ,
subject to theboundary conditions
(3.11)U0(ξ0)=0,U0(ξN)=0,(3.12)∑k=0N𝒟Nk2U0(ξk)=0,∑k=0N𝒟0kU0(ξk)=0,
where(3.13)A=16𝒟4+8a1𝒟3+4a2𝒟2+2a3𝒟-3ReI,U0=[U0(ξ0),U0(ξ1),…,U0(ξN)]T,Φ=[ϕ(y0),ϕ(y1),…,ϕ(yN)]T,(3.14)as=diag([as(y0),as(y1),…,as(yN-1),as(yN)]),s=1,2,3.
In the above definitions, the superscript T denotes transpose, diag is a diagonal matrix and I is an identity matrix of size (N+1)×(N+1).To implement the boundary conditions (3.11), we delete the first and the last rows and columns of A and delete the first and last rows of U0 and Φ. The boundary conditions (3.12) are imposed on the resulting first and last rows of the modified matrix A and setting the resulting first and last rows of the modified matrix Φ to be zero. The values of [U0(ξ1),U0(ξ2),…,U0(ξN-1)] are then determined from (3.15)U0=A-1Φ.
To find the SHAM solutions of (3.2) we begin by defining the following linear operator:(3.16)ℒ[Ũ(ξ;q)]=16∂4Ũ∂ξ4+8a1∂3Ũ∂ξ3+4a2∂2Ũ∂ξ2+2a3∂Ũ∂ξ-3ReŨ,
where q∈[0,1] is the embedding parameter, and Ũ(ξ;q) is an unknown function.The zeroth order deformation equation is given by(3.17)(1-q)ℒ[Ũ(ξ;q)-U0(ξ)]=qℏ{𝒩[Ũ(ξ;q)]-Φ},
where ℏ is the nonzero convergence controlling auxiliary parameter and 𝒩 is a nonlinear operator given by(3.18)𝒩[Ũ(ξ;q)]=16∂4Ũ∂ξ4+8a1∂3Ũ∂ξ3+4a2∂2Ũ∂ξ2+2a3∂Ũ∂ξ-3ReŨ+8Re(U∂3Ũ∂ξ3-∂Ũ∂ξ∂2Ũ∂ξ2).Differentiating (3.17) 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.19)ℒ[Um(ξ)-χmUm-1(ξ)]=ℏRm,
subject to the boundary conditions(3.20)Um(-1)=Um(1)=Um′′(-1)=Um′(1)=0,
where (3.21)Rm(ξ)=16Um-11V+8a1Um-1′′′+4a2Um-1′′+2a3Um-1′-3ReUm-1+8Re∑n=0m-1(UnUm-1-n′′′-Un′Um-1-n′′)-ϕ(y)(1-χm),(3.22)χm={0,m≤11,m>1.
Applying the Chebyshev pseudospectral transformation on (3.19)–(3.21) gives(3.23)AUm=(χm+ℏ)AUm-1-ℏ(1-χm)Φ+ℏPm-1
subject to the boundary conditions
(3.24)Um(ξ0)=0,Um(ξN)=0,(3.25)∑k=0N𝒟Nk2Um(ξk)=0,∑k=0N𝒟0kUm(ξk)=0,
where A and Φ, are as defined in (3.13) and(3.26)Um=[Um(ξ0),Um(ξ1),…,Um(ξN)]T,Pm-1=8Re∑n=0m-1[Un(𝒟3Um-1-n)-(𝒟Un)(𝒟2Um-1-n)].
To implement the boundary conditions (3.24) we delete the first and last rows of Pm-1 and Φ and delete the first and last rows and first and last columns of A in (3.23). The boundary conditions (3.25) are imposed on the first and last row of the modified A matrix on the left side of the equal sign in (3.23). The first and last rows of the modified A matrix on the right of the equal sign in (3.23) are the set to be zero. This results in the following recursive formula for m≥1:(3.27)Um=(χm+ℏ)A-1ÃUm-1+ℏA-1[Pm-1-(1-χm)Φ].
Thus, starting from the initial approximation, which is obtained from (3.15), higher-order approximations Um(ξ) for m≥1 can be obtained through the recursive formula (3.27).


## 4. Successive Linearisation Method

In this section, we solve (2.6) using the successive linearisation method. The main assumptions underpinning the use of the successive linearisation method are the following.  (i)
The unknown functionF(y) maybe expanded as
(4.1)F(y)=Fi(y)+∑m=0i-1Fm(y),i=1,2,3,…,
where Fi are unknown functions and Fm (m≥1) are approximations which are obtained by recursively solving the linear part of the equation that results from substituting (4.1) in the governing equation (2.6).(ii)
Fi becomes progressively smaller as i becomes large, that is,(4.2)limi→∞Fi=0.
Substituting (4.1) in the governing equation gives(4.3)Fi(iv)+a1,i-1Fi′′′+a2,i-1Fi′′+a3,i-1Fi′+a4,i-1Fi+Re(FiFi′′′-Fi′Fi′′)=ri-1,
where the coefficient parameters ak,i-1, (k=1,2,3,4), and ri-1 are defined as(4.4)a1,i-1=Re∑m=0i-1Fm+αy,a2,i-1=-Re∑m=0i-1Fm′+3α,a3,i-1=-Re∑m=0i-1Fm′′,a4,i-1=Re∑m=0i-1Fm′′′,ri-1=-(∑m=0i-1Fm(iv)+αy∑m=0i-1Fm′′′+3α∑m=0i-1Fm′′)-Re(∑m=0i-1Fm∑m=0i-1Fm′′′-∑m=0i-1Fm′∑m=0i-1Fm′′).The SLM algorithm starts from the initial approximation(4.5)F0(y)=12(3+β)y-12(1+3β)y3+βy4,
which is chosen to satisfy the boundary conditions (2.7)–(2.8). The parameter β in (4.5) is an arbitrary constant which when varied results in multiple solutions. The subsequent solutions for Fm, m≥1 are obtained by successively solving the linearized form of (4.3) and which is given as(4.6)Fi(iv)+a1,i-1Fi′′′+a2,i-1Fi′′+a3,i-1Fi′+a4,i-1Fi=ri-1,
subject to the boundary conditions(4.7)Fi(0)=0,Fi′′(0)=0,Fi(1)=1,Fi′(1)=0.
Once each solution for Fi, (i≥1) has been found from iteratively solving (4.6) for each i, the approximate solution for F(y) is obtained as(4.8)F(y)≈∑m=0MFm(y),
where M is the order of SLM approximation. Since the coefficient parameters and the right hand side of (4.6), for i=1,2,3,…, are known from previous iterations, (4.6) can easily be solved using analytical means or any numerical methods such as finite differences, finite elements, Runge-Kutta-based shooting methods, or collocation methods. In this paper, (4.6) is integrated using the Chebyshev spectral collocation method [32, 33] as described in the previous section.Applying the Chebyshev spectral method to (4.6) leads to the matrix equation(4.9)Ai-1Fi=Ri-1,
in which A is an (N+1)×(N+1) square matrix and Y and R are (N+1)×1 column vectors defined by(4.10)A=D4+a1,i-1D3+a2,i-1D2+a3,i-1D+a4,i-1,Ri-1=ri-1,
with(4.11)Fi=[Fi(x0),Fi(x1),…,Fi(xN-1),Fi(xN)]T,ri-1=[ri-1(x0),ri-1(x1),…,ri-1(xN-1),ri-1(xN)]T.
In the above definitions, N is the number of collocation points, x=2y-1, ak,i-1, (k=1,2,3,4) are diagonal matrices of size (N+1)×(N+1), and D=2𝒟. After modifying the matrix system (4.9) to incorporate boundary conditions, the solution is obtained as(4.12)Ỹi=Ai-1-1Ri-1.


## 5. Results and Discussion

In this section, we compare the results obtained using the various methods: the SHAM, the SLM, and the numerical approximations with those obtained using the HAM in Dinarvand and Rashidi [30] and the homotopy-Páde method in Xu et al. [24]. The solution obtained using most numerical solutions depends on the initial approximation. Using different initial guesses can give rise to multiple solutions. Multiple solutions were obtained if the initial guess in (4.5) is used in the SHAM method in place of F0(y) in (3.1). In this paper, it was observed that using different values of β results in multiple solutions. For the multiple solutions comparison was made against the HAM results of [24].An optimalℏ value can easily be sought that can considerably improve the convergence rate of the results. However, for comparison purposes we used ℏ=-1. It is however worth noting, as pointed out in Dinarvand et al. [2], that when ℏ=-1, the solution series obtained by the HAM is the same solution series obtained by means of the homotopy perturbation method.In Table1 we compare the values of F(y) when α=-1 and Re=-2,0 and 2 with the numerical and the HAM results reported in Dinarvand and Rashidi [30]. In [30], convergence up to six decimal places was achieved at the sixth order of the HAM approximation for Re=0 and 2. In this study, the same level of convergence and accuracy was achieved at the first order approximation for the same values of Re. For Re=-2 the convergence of the homotopy analysis method series solution was achieved at the eighth order of approximation while with the spectral homotopy analysis method series solution gives the same level of convergence at the second order.Table 1
Comparison of the numerical results against the SHAM approximate solutions forF(y) when α=-1 with N=60 and ℏ=-1.
Rey0th order1st order2nd order3rd orderNumericalRef. [30]-20.20.2738280.2738310.2738320.2738320.2738320.2738320.40.5328270.5328390.5328390.5328390.5328390.5328390.60.7594420.7594670.7594680.7594680.7594680.7594680.80.9289670.9289900.9289900.9289900.9289900.92899000.20.2794490.2794490.2794490.2794490.2794490.2794490.40.5422430.5422430.5422430.5422430.5422430.5422430.60.7689500.7689500.7689500.7689500.7689500.7689500.80.9338890.9338890.9338890.9338890.9338890.93388920.20.2839960.2839830.2839830.2839830.2839830.2839830.40.5497590.5497380.5497380.5497380.5497380.5497380.60.7763280.7763060.7763060.7763060.7763060.7763060.80.9375180.9375070.9375070.9375070.9375070.937507In Table2, we demonstrate the computational efficiency of the SHAM solution for large values of Re. As has been noted in the introduction, some semi-analytical methods fail to converge at large values of Re, for example, Dinarvand et al. [2] have shown that for |Re|>9.5 the HPM fails to converge. However, in using the SHAM convergence up seven decimal places is achieved at the third order of approximation for values of Re as large as Re=150. For 5≤Re<100 convergence up six decimal places is achieved at the second order. In Table 3, Re=2 is fixed and F′′(1) evaluated for -2.5≤α≤2.5. Convergence up to seven digits is achieved at the first order for α=0, at the second order for -2≤α≤-0.5, and at the third-order approximation for α=-2.5. Clearly the SHAM gives faster convergence than the HAM under the same conditions.Table 2
Comparison of the numerical results against the SHAM approximate solutions forF′′(1) when α=-1 with N=60 and ℏ=-1 for different values of Re.
Re0th order1st order2nd order3rd order4th orderNumerical0−3.8213722−3.8213723−3.8213723−3.8213723−3.8213723−3.82137235−3.1725373−3.1731774−3.1731800−3.1731800−3.1731800−3.173180010−2.9069253−2.9069653−2.9069654−2.9069654−2.9069654−2.906965415−2.7783086−2.7784366−2.7784369−2.7784369−2.7784369−2.778436920−2.7056150−2.7060556−2.7060557−2.7060557−2.7060557−2.706055725−2.6596223−2.6603920−2.6603907−2.6603907−2.6603907−2.660390730−2.6281102−2.6291718−2.6291682−2.6291682−2.6291682−2.629168240−2.5879160−2.5894333−2.5894247−2.5894247−2.5894247−2.589424750−2.5634446−2.5652868−2.5652732−2.5652733−2.5652733−2.5652733100−2.5139141−2.5165255−2.5164964−2.5164967−2.5164967−2.5164967150−2.4972813−2.5001845−2.5001479−2.5001484−2.5001484−2.5001484Table 3
Comparison of the numerical results against the SHAM approximate solutions forF′′(1) when α=-1 with Re=10 and ℏ=-1 for different values of N.
N0th order1st order2nd order3rd order4th orderNumerical10−2.7059921−2.7064302−2.7064302−2.7064302−2.7064302−2.706055715−2.7056141−2.7060549−2.7060550−2.7060549−2.7060549−2.706055720−2.7056150−2.7060556−2.7060557−2.7060557−2.7060557−2.706055730−2.7056150−2.7060556−2.7060557−2.7060557−2.7060557−2.706055740−2.7056150−2.7060556−2.7060557−2.7060557−2.7060557−2.706055760−2.7056150−2.7060556−2.7060557−2.7060557−2.7060557−2.7060557Table3 also gives a comparison of the SHAM approximate results with the numerical results generated using different values of N. It is evident however that for small values of N (say less than N=20) the SHAM results are not very accurate. Accuracy however improves with an increase in N.Table4 gives a comparison of the numerical and the SHAM results of F′′(1) when -2.5≤α≤2.5. Convergence is generally achieved at either the third- or the fourth-order of the SHAM approximation. Figures 1 and 2 give a comparison of the numerical and the SHAM solutions for the characteristic mean-flow function F(y)=-ν/c and F'(y)=μc/x at different Reynolds numbers and α. The mean-flow function F(y), increases with increasing (positive) values of Re and α while F'(y) decreases, which makes sense since F(y) is inversely proportional to the injection or suction coefficient c=α/Re while F'(y) is directly proportional. Figure 1 further illustrates the efficiency of the solution method with excellent agreement for Re as large as 200.Table 4
Comparison of the numerical results against the SHAM approximate solutions forF′′(1) when Re=2 with N=60 and ℏ=-1 for different values of α.
α1st order2nd order3rd order4th orderNumerical−2.5−4.5258487−4.5258506−4.5258505−4.5258505−4.5258505−2.0−4.1673848−4.1673892−4.1673892−4.1673892−4.1673892−1.5−3.8209704−3.8209740−3.8209740−3.8209740−3.8209740−1.0−3.4873966−3.4873982−3.4873982−3.4873982−3.4873982−0.5−3.1674332−3.1674334−3.1674334−3.1674334−3.16743340.0−2.8618116−2.8618116−2.8618116−2.8618116−2.86181160.5−2.5712067−2.5712055−2.5712055−2.5712055−2.57120561.0−2.2962176−2.2962075−2.2962076−2.2962076−2.29620761.5−2.0373492−2.0373088−2.0373092−2.0373092−2.03730922.0−1.7949940−1.7948795−1.7948810−1.7948810−1.79488102.5−1.5694172−1.5691503−1.5691557−1.5691556−1.5691556Comparison between numerical and SHAM approximate solution ofF(y) and F′(y) for different values of Re when α=-1 when ℏ=-1.14 (for Re=-10) and ℏ=-1 (for Re=0,200).
(a)(b)Comparison between numerical and SHAM approximate solution ofF(y) and F′(y) for different values of α when Re=2 when ℏ=-1 (for α=-5,-2.5,0,2.5) and ℏ=-0.94 (for α=5).
(a)(b)Using the initial approximationF0(y) in (4.5) with different values of β in place of (3.1) in the SHAM implementation leads to multiple solutions. When β=0 and -5, the SHAM gives the multiple solutions observed in Xu et al. [24]. A comparison of the SHAM results against the HAM results reported in [24] is presented in Table 5. It is evident that the SHAM results converge much more rapidly than the HAM results of [24] for both branches of the solution.Table 5
Comparison between the multiple solutions of HAM results (see [24]) and the present SHAM results in the case of Re=-10 and α=4.
HAM solution [24]SHAM solutionOrder of approx.F′(0)F′′′(0)Order of approx.F′(0)F′′′(0)First solution100.624 1618.267 5650.6255498.24662200.624 9678.256 03100.6250168.25532300.625 0058.255 48150.6250088.25544400.625 0078.255 45200.6250078.25545500.625 0078.255 45250.6250078.25545Second solution10-1.189 9535.89845-1.19052935.8564720-1.190 0335.847410-1.19032235.8541430-1.190 4135.855515-1.19032335.8541640-1.190 3135.853920-1.19032335.8541650-1.190 3335.854225-1.19032335.85416Tables6 and 7 give, first, the analytical approximations of F′(0)/Re and F′′′(0)/Re for the two solutions obtained using the successive linearisation method. Secondly, the tables give a direct comparison of the convergence rates of the SLM and the [m,m] homotopy-Padé method used by Xu et al. [24]. The fourth-order SLM approximation gives the same level of accuracy as the twenty-fourth-order of the [m,m] homotopy-Páde approximation, which suggests that the successive linearisation method is much more computationally efficient and accurate compared to the [m,m] homotopy-Páde approximation (although it is not clear at this stage whether this could be attributed to the use of a more efficient initial guess).Table 6
Comparison ofF′(0)/Re and F′′′(0)/Re obtained at different orders for the SLM, and the [m,m] homotopy-Páde approximation when Re=-10 and α=4. For the SLM first solution, β=-1, and N=25, and for the SLM second solution, β=-5, and N=25.
[m,m] homotopy-Padé [24]SLMorderF′(0)/ReF′′′(0)/ReorderF'(0)F′′′(0)First solution40.6244787328.26544422220.6244858958.2628376580.6250053368.25547742230.6250075168.25544430160.6250073958.25544612740.6250073968.25544612200.6250073968.25544612560.6250073968.25544612240.6250073968.25544612480.6250073968.25544612Second solution4-1.21989136.010912-1.19093435.860428-1.17846535.178783-1.19032335.8541616-1.19031935.854094-1.19032335.8541620-1.19032335.854156-1.19032335.8541624-1.19032335.854168-1.19032335.85416Table 7
Comparison ofF′(0)/Re and F′′′(0)/Re obtained at different orders for the SLM, and the [m,m] homotopy-Páde approximation when Re=-11 and α=3/4. For the SLM first solution, β=-5, and N=25, β=1, and N=25 for the SLM second solution, and for the SLM third solution β=3, and N=25.
[m,m] homotopy-Páde [24]SLMorderF′(0)/ReF′′′(0)/ReorderF′(0)F′′′(0)First solution4-1.023162124.29258512-1.325016827.86404868-1.023770024.28637973-1.076577724.909500616-1.023771224.28630884-1.025952724.311998720-1.023771224.28630886-1.023771224.286308824-1.023771224.28630888-1.023771224.2863088Second solution40.166859010.28286020.21349509.68258180.167998010.23945130.171842010.216566160.169351810.24502640.169357310.245102200.169353210.24515060.169353110.245151240.169353210.24515180.169353110.245151Third solution4……22.76262-15.526682.81591-15.895032.76111-15.5123162.76154-15.515742.76111-15.5122202.76113-15.512362.76111-15.5122242.76111-15.512382.76111-15.5122


## 6. Conclusion

In this paper, we have used the spectral homotopy analysis method (SHAM) and the successive linearisation method (SLM) to solve a fourth-order nonlinear boundary value problem that governs the two-dimensional Laminar flow between two moving porous walls. Multiple solutions recently reported in Xu et al. [24] have been obtained, depending on the initial approximation used. Comparison of the computational efficiency and accuracy of the results between the current methods and the previous homotopy analysis method results described in Dinarvand and Rashidi [30] and Xu et al. [24] has been made. Our simulations show that the convergence of the SHAM solution series to the numerical solution (up to six decimal place accuracy) is achieved at the second order (for Re=-2) and first order for R=0,2. In contrast to the standard homotopy analysis method (c.f. Dinarvand and Rashidi [30]) convergence was achieved at the eighth order (for Re=-2) and sixth order for Re=0,2. The SHAM is apparently more efficient because it offers more flexibility in choosing linear operators compared to the standard HAM. It is however important to note that if the same initial guess and linear operators were to be used, the two methods would give the same solution.We have further shown that notwithstanding the acceleration of the convergence ratio of the homotopy series by means of the homotopy-Padé technique, the successive linearisation techniques is more computationally efficient (although this again could be due, in part, to the use of a different initial guess) and gives accurate results.



---

*Source: 528956-2010-09-19.xml*