Advances in Mathematical Physics

Volume 2015, Article ID 642835, 4 pages

http://dx.doi.org/10.1155/2015/642835

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

^{1}Department of Mathematics, University of Peshawar, Peshawar, Pakistan^{2}Department of Mathematics, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand^{3}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 aswhere 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 iswith boundary conditionsIts solution isFirst order problem iswith boundary conditionshaving solutionSecond order problem iswith boundary conditionsIts solution becomesWe 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.