# Comparison of Optimal Homotopy Asymptotic and Adomian Decomposition Methods for a Thin Film Flow of a Third Grade Fluid on a Moving Belt

Authors: Fazle Mabood; Nopparat Pochai
Journal: Advances in Mathematical Physics
 (2015)

---

## Abstract

We have investigated a thin film flow of a third grade fluid on a moving belt using a powerful and relatively new approximate analytical technique known as optimal homotopy asymptotic method (OHAM). The variation of velocity profile for different parameters is compared with the numerical values obtained byRunge-Kutta Fehlberg fourth-fifth order method and with Adomian Decomposition Method (ADM). An interesting result of the analysis is that the three terms OHAM solution is more accurate than five terms of the ADM solution and this thus confirms the feasibility of the proposed method.


---

## Body

## 1. Introduction

Many physical systems in fluid mechanics generally lead to nonlinear ordinary or partial differential equations. Due to complexity of Non-Newtonian fluid, it is difficult to solve nonlinear differential equation. A second grade fluid is one of the most acceptable fluids in this sub clam of Non-Newtonian fluids, because of its mathematical simplicity in comparison to third grade and fourth grade fluids. In related literature many authors have effectively treated the complicated nonlinear equations governing the flow of a third grade fluid [1, 2].Since the Non-Newtonian fluids are of the great challenges in the solution of governing nonlinear differential equations, many numerical and analytical techniques have been proposed by many researchers. But an efficient approximate analytical solution still finds enormous appreciations. Keeping this fact in mind, we have solved the governing nonlinear equation of the present problem using the two techniques. It is important to mention here that the analytical and numerical solutions are in a good agreement but better than the results of Siddiqui et al. [3].In this study, it is also observed that the optimal homotopy asymptotic method is a powerful approximate analytical tool that is simple and straightforward and does not require the existence of any small or large parameter as does traditional perturbation method. Optimal homotopy asymptotic method has successfully been applied to a number of nonlinear problems arising in the science and engineering by various researchers [4–9]. This proves the validity and acceptability of OHAM as a useful solution technique.This paper is organized as follows. First in Section2, we formulate the problem. In Section 3 we present basic principles of OHAM. The OHAM solution is given in Section 4. In Section 5, we analyze the comparison of the solution using OHAM with existing solution of ADM. Section 6 is devoted for the conclusion.


## 2. Governing Equation

The thin film flow of third grade fluid on a moving belt is governed by the following nonlinear boundary value problem [3]:(1)
d
2
v
d
x
2
+
6
β
2
+
β
3
μ
d
v
d
x
2
d
2
v
d
x
2
-
ρ
g
μ
=
0
,
v
0
=
U
0
,
d
v
d
x
=
0
at
x
=
δ
,where v is the fluid velocity, ρ is the density, μ is the dynamic viscosity, β
2 and β
3 are the material constants of the third grade fluid, g is acceleration due to gravity, δ is the uniform thickness of the fluid film, and U
0 is the speed of the belt.Here, we introduce the following dimensionless variables:(2)
x
∗
=
x
δ
,
v
∗
=
v
U
0
,
β
=
β
2
+
β
3
μ
δ
2
U
0
,
m
=
ρ
g
μ
U
0
δ
2
.From (1)-(2), we obtain the dimensionless form as(for simplicity we removed∗)(3)
d
2
v
d
x
2
+
6
β
d
v
d
x
2
d
2
v
d
x
2
-
m
=
0
,
subject to the boundary conditions(4)
v
0
=
1
,
v
′
1
=
0
.


## 3. Optimal Homotopy Asymptotic Method

We review the basic principles of OHAM as expounded in [4–8] in five steps.(i) Let us consider the following differential equation:(5)
A
v
τ
+
a
τ
=
0
,
τ
∈
Ω
,where Ω is problem domain, A
(
v
)
=
L
(
v
)
+
N
(
v
), where L, N are linear and nonlinear operator, v
(
τ
) is an unknown function, and a
(
τ
) is a known function.(ii) Construct an optimal homotopy equation as(6)
1
-
p
L
ϕ
τ
;
p
+
a
τ
-
H
p
A
ϕ
τ
;
p
+
a
τ
=
0
,where 0
≤
p
≤
1 is an embedding parameter and H
(
p
)
=
∑
k
=
1
m
p
k
C
k is auxiliary function on which the convergence of the solution greatly depends. The auxiliary function H
(
p
) also adjust the convergence domain and control the convergence region. According to the new development in OHAM [9] the more generalized form of the auxiliary function is H
(
τ
;
p
,
C
i
)
=
p
H
1
(
τ
,
C
i
)
+
p
2
H
2
(
τ
,
C
i
)
+
…, where H
i
(
τ
,
C
i
), i
=
1,2
,
… are auxiliary functions depending upon τ and unknown parameters C
j. This means that we could have more convergence-control parameters even at the first order of approximation.(iii) Expandϕ
(
τ
;
p
,
C
j
) in Taylor’s series about p; one has an approximate solution:(7)
ϕ
τ
;
p
,
C
j
=
v
0
τ
+
∑
k
=
1
∞
v
k
τ
,
C
j
p
k
,
j
=
1,2
,
3
,
…
.Many researchers have observed that the convergence of the series (7) depends upon C
j
,
(
j
=
1,2
,
…
,
m
), if it is convergent then, we obtain(8)
v
~
=
v
0
τ
+
∑
k
=
1
m
v
k
τ
;
C
j
.(iv) Substituting (8) in (6), we have the following residual:(9)
R
τ
;
C
j
=
L
v
~
τ
;
C
j
+
a
τ
+
N
v
~
τ
;
C
j
.If R
(
τ
;
C
j
)
=
0, then v
~ will be the exact solution. For nonlinear problems, generally this will not be the case. For determining C
j
,
(
j
=
1,2
,
…
,
m
), Galerkin’s Method, Ritz Method, or the method of least squares can be used.(v) Finally, substitute these constants in (8) and one can get the approximate solution.


## 4. Solution of the Problem via OHAM

According to the OHAM, applying (6) to (3)(10)
1
-
p
v
′
′
-
H
p
v
′
′
+
6
β
v
′
2
v
′
′
-
m
=
0
,where primes denote differentiation with respect to x.We considerv and H
(
p
) as the following:(11)
v
=
v
0
+
p
v
1
+
p
2
v
2
,
H
p
=
p
C
1
+
p
2
C
2
.Put (11) in (10) and arrange the terms according to the powers of p to get zeroth, first, and second order problems as follows.Zeroth order problem is(12)
v
0
′
′
x
=
mwith boundary conditions(13)
v
0
0
=
1
,
v
0
′
1
=
0
.Its solution is(14)
v
0
=
1
2
2
-
2
m
x
+
m
x
2
.First order problem is(15)
v
1
′
′
x
,
C
1
=
-
m
-
m
C
1
+
6
β
C
1
v
0
′
2
v
0
′
′
+
1
+
C
1
v
0
′
′
xwith boundary conditions(16)
v
1
0
=
0
,
v
1
′
1
=
0
,having solution(17)
v
1
x
,
C
1
=
1
2
-
4
m
3
β
x
C
1
+
6
m
3
β
x
2
C
1
-
4
m
3
β
x
3
C
1
+
m
3
β
x
4
C
1
.Second order problem is(18)
v
2
′
′
x
,
C
1
,
C
2
=
-
m
C
2
+
C
2
v
0
′
′
+
6
β
C
2
v
0
′
2
v
0
′
′
+
12
β
C
1
v
0
′
v
1
′
v
0
′
′
+
6
β
C
1
v
0
′
2
v
1
′
′
+
1
+
C
1
v
1
′
′with boundary conditions(19)
v
2
0
=
0
,
v
2
′
1
=
0
.Its solution becomes(20)
v
2
x
,
C
1
,
C
2
=
1
2
-
4
m
3
β
x
C
1
+
6
m
3
β
x
2
C
1
-
4
m
3
β
x
3
C
1
+
m
3
β
x
4
C
1
-
4
m
3
β
x
C
1
2
-
24
m
5
β
2
x
C
1
2
+
6
m
3
β
x
2
C
1
2
+
60
m
5
β
2
x
2
C
1
2
-
4
m
3
β
x
3
C
1
2
-
80
m
5
β
2
x
3
C
1
2
+
m
3
β
x
4
C
1
2
+
60
m
5
β
2
x
4
C
1
2
-
24
m
5
β
2
x
5
C
1
2
+
4
m
5
β
2
x
6
C
1
2
-
4
m
3
β
x
C
2
+
6
m
3
β
x
2
C
2
-
4
m
3
β
x
3
C
2
+
m
3
β
x
4
C
2
.We obtain the three terms solution using OHAM for p
=
1
(21)
v
~
x
,
C
1
,
C
2
=
v
0
x
+
v
1
x
,
C
1
+
v
2
x
,
C
1
,
C
2
.From least squares method we obtain the unknown convergent constants C
1, C
2 in (21).For the particular case, ifβ
=
0.5 and m
=
0.2, we have C
1
=
-
0.877411, C
2
=
-
0.003097.


## 5. Results and Discussion

Table1 shows the comparison of absolute error between OHAM (three terms) and ADM (five terms) [3]. It is noteworthy to mention here that OHAM low error is remarkable, while the effectiveness of the proposed method (OHAM) can be seen from Figure 1. The effect of fluid parameter β is displayed in Figure 2. From Figure 2 it is found that the boundary layer thickness is increased with an increase in fluid parameter β whereas the value of m is preset. Whilst Figure 3 depicts an increase in m for the fixed value of fluid parameter β, a decrease can be seen in boundary layer thickness. However, opposite observation is made by comparing of Figures 2 and 3.Table 1
Comparison of absolute error using OHAM (three terms) and ADM (five terms) [3] for β
=
0.5, m
=
0.5.
x
OHAM
ADM
NM
Error (ADM)
Error (OHAM)
0.0
1
1
1
0
0
0.1
0.959206
0.962523
0.959268
3.2
×
10
-
3
6.2
×
10
-
5
0.2
0.921830
0.926065
0.9219695
4.1
×
10
-
3
1.3
×
10
-
4
0.3
0.888047
0.892574
0.8882986
4.3
×
10
-
3
2.5
×
10
-
4
0.4
0.858081
0.862771
0.8584654
4.3
×
10
-
3
3.8
×
10
-
4
0.5
0.832177
0.836998
0.8326892
4.3
×
10
-
3
5.1
×
10
-
4
0.6
0.810581
0.815502
0.8111930
4.3
×
10
-
3
6.1
×
10
-
4
0.7
0.793520
0.798502
0.7941932
4.3
×
10
-
3
6.7
×
10
-
4
0.8
0.781186
0.786196
0.7818867
4.3
×
10
-
3
7.1
×
10
-
4
0.9
0.773724
0.778742
0.7744327
4.3
×
10
-
3
7.1
×
10
-
4
1.0
0.771227
0.776245
0.7719358
4.3
×
10
-
3
7.1
×
10
-
4Figure 1
Comparison of velocity profile using OHAM (three terms) and ADM [3] (five terms) for β
=
0.5, m
=
0.3.Figure 2
Effects on velocity profile for various values ofβ at m
=
0.5.Figure 3
Effects on velocity profile for various values ofm at β
=
0.5.


## 6. Conclusion

Optimal homotopy asymptotic method is employed to investigate the approximate solution for a thin film flow of third grade fluid on a moving belt. Both numeric and analytic results are obtained for the problem. The results are sketched and discussed for the fluid parameterβ and for constant m. From the study it is revealed that the solution using OHAM is better than ADM results. Finally, we conclude that OHAM provide a simple and easy way to control and adjust the convergence region for strong nonlinearity and is applicable to highly nonlinear fluid problems.



---

*Source: 642835-2015-08-31.xml*