Advances in Mathematical Physics
Volume 2015 (2015), Article ID 642835, 4 pages
http://dx.doi.org/10.1155/2015/642835

Research Article

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

Fazle Mabood¹ and Nopparat Pochai2,3

¹Department of Mathematics, University of Peshawar, Peshawar, Pakistan
²Department of Mathematics, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand
³Centre of Excellence in Mathematics, CHE, Si Ayutthaya Road, Bangkok 10400, Thailand

Received 24 March 2015; Accepted 11 May 2015

Academic Editor: John D. Clayton

Copyright © 2015 Fazle Mabood and Nopparat Pochai. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 by Runge-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.

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 Section 2, 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]:

where is the fluid velocity, is the density, is the dynamic viscosity, and are the material constants of the third grade fluid, is acceleration due to gravity, is the uniform thickness of the fluid film, and is the speed of the belt.

Here, we introduce the following dimensionless variables:

From (1)-(2), we obtain the dimensionless form as

 (for simplicity we removed )

 subject to the boundary conditions

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:

where is problem domain, , where , are linear and nonlinear operator, is an unknown function, and is a known function.

(ii) Construct an optimal homotopy equation as

where is an embedding parameter and is auxiliary function on which the convergence of the solution greatly depends. The auxiliary function 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 , where , are auxiliary functions depending upon and unknown parameters . This means that we could have more convergence-control parameters even at the first order of approximation.

(iii) Expand in Taylor’s series about ; one has an approximate solution:

Many researchers have observed that the convergence of the series (7) depends upon , if it is convergent then, we obtain

(iv) Substituting (8) in (6), we have the following residual:

If , then will be the exact solution. For nonlinear problems, generally this will not be the case. For determining , 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)

where primes denote differentiation with respect to .

We consider and as the following:

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

with boundary conditions

Its solution is

First order problem is

with boundary conditions

having solution

Second order problem is

with boundary conditions

Its solution becomes

We obtain the three terms solution using OHAM for

From least squares method we obtain the unknown convergent constants , in (21).

For the particular case, if and , we have , .

5. Results and Discussion

Table 1 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 is preset. Whilst Figure 3 depicts an increase in 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 , .

OHAM
ADM
NM
Error (ADM)
Error (OHAM)

0.0
1
1
1
0
0

0.1
0.959206
0.962523
0.959268

0.2
0.921830
0.926065
0.9219695

0.3
0.888047
0.892574
0.8882986

0.4
0.858081
0.862771
0.8584654

0.5
0.832177
0.836998
0.8326892

0.6
0.810581
0.815502
0.8111930

0.7
0.793520
0.798502
0.7941932

0.8
0.781186
0.786196
0.7818867

0.9
0.773724
0.778742
0.7744327

1.0
0.771227
0.776245
0.7719358

Figure 1: Comparison of velocity profile using OHAM (three terms) and ADM [3] (five terms) for , .

Figure 2: Effects on velocity profile for various values of at .

Figure 3: Effects on velocity profile for various values of at .

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 . 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.

Conflict of Interests

The authors declare no conflict of interests.

Acknowledgments

The second author appreciates the Center of Excellence in Mathematics, the Commission on Higher Education, Thailand. The authors greatly appreciate the valuable comments received from the referees.

References

A. M. Siddiqui, R. Mahmood, and Q. K. Ghori, “Homotopy perturbation method for thin film flow of a third grade fluid down an inclined plane,” Chaos, Solitons & Fractals, vol. 35, no. 1, pp. 140–147, 2008.

A. M. Siddiqui, R. Mahmood, and Q. K. Ghori, “Thin film flow of a third grade fluid on a moving belt by He's homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 7–14, 2006.

A. M. Siddiqui, A. A. Farooq, T. Haroon, and B. S. Babcock, “A comparison of variational iteration and Adomian decomposition methods in solving nonlinear thin film flow problems,” Applied Mathematical Sciences, vol. 6, no. 97-100, pp. 4911–4919, 2012.

N. Herişanu, V. Marinca, T. Dordea, and G. Madescu, “A new analytical approach to nonlinear vibration of an electrical machine,” Proceedings of the Romanian Academy Series A: Mathematics Physics Technical Sciences Information Science, vol. 9, no. 3, pp. 229–236, 2008.

V. Marinca and N. Herisanu, “Optimal homotopy perturbation method for strongly nonlinear differential equations,” Nonlinear Science Letters A, vol. 1, no. 3, pp. 273–280, 2010.

V. Marinca and N. Herişanu, “Application of optimal homotopy asymptotic method for solving nonlinear equations arising in heat transfer,” International Communications in Heat and Mass Transfer, vol. 35, no. 6, pp. 710–715, 2008.

F. Mabood, W. A. Khan, and A. I. M. Ismail, “Optimal homotopy asymptotic method for heat transfer in hollow sphere with robin boundary conditions,” Heat Transfer—Asian Research, vol. 43, no. 2, pp. 124–133, 2014.

V. Marinca and N. Herişanu, “Determination of periodic solutions for the motion of a particle on a rotating parabola by means of the optimal homotopy asymptotic method,” Journal of Sound and Vibration, vol. 329, no. 9, pp. 1450–1459, 2010.

N. Herisanu, V. Marinca, and G. Madescu, “An analytical approach to non-linear dynamical model of a permanent magnet synchronous generator,” Wind Energy, 2014.