Journal of Computational Engineering

Volume 2014 (2014), Article ID 102197, 11 pages

http://dx.doi.org/10.1155/2014/102197

Research Article

A Modified Homotopy Perturbation Transform Method for Transient Flow of a Third Grade Fluid in a Channel with Oscillating Motion on the Upper Wall

Mohammed Abdulhameed, Rozaini Roslan, and Mahathir Bin Mohamad

Centre for Research in Computational Mathematics, Faculty of Science, Technology and Human Development, Universiti Tun Hussein Onn, 86400 Batu Pahat, Johor, Malaysia

Received 20 January 2014; Revised 19 June 2014; Accepted 19 June 2014; Published 17 July 2014

Academic Editor: Clement Kleinstreuer

Copyright © 2014 Mohammed Abdulhameed et al. 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.

A new analytical algorithm based on modified homotopy perturbation transform method is employed for solving the transient flow of third grade fluid in a porous channel generated by an oscillating upper wall. This method incorporates the He’s polynomial into the HPM, combined with Laplace transform. Comparison with HPM and OHAM analytical solutions reveals that the proposed algorithm is highly accurate. This proves the validity and great potential of the proposed algorithm as a new kind of powerful analytical tool for transient nonlinear problems. Graphs representing the solutions are discussed, and appropriate conclusions are drawn.

1. Introduction

The equations describing the motion of non-Newtonian fluids are strongly of nonlinear higher order than the Navier-Stokes equation for Newtonian fluids. These nonlinear equations form a very complex structure, with a small number of exact solutions. Mostly, numerical methods have largely been used to handle these equations. The class of problems with known exact solution is related to the problem for infinite flat plate. The related studies in the recent years are as follows: Fakhar et al. [1] examine the exact unsteady flow of an incompressible third grade fluid along an infinite plane porous plate. They obtained results by applying a translational type of symmetries combined with finite difference method. Danish and Kumar [2] analysed a steady flow of a third grade between two parallel plates using similarity transformation. Abdulhameed et al. [3] consider an unsteady viscoelastic fluid of second grade for an infinite plate. They applied Laplace transform together with the regular perturbation techniques to obtain the exact solution. Ayub et al. [4] analysed the problem of steady flow of a third grade fluid for an infinite plate porous plate using homotopy analysis method (HAM).

Homotopy perturbation method developed by He [5] for solving linear and nonlinear initial-boundary value problem merges two techniques, the perturbation and standard homotopy. Recently, the homotopy perturbation method has been modified by some scientists to obtain more accurate results and rapid convergence and also to reduce the amount of computation. Ghorbani [6] introduced He’s polynomials based on homotopy perturbation method for nonlinear differential equations. The homotopy perturbation transform method (HPTM) introduced by Khan and Wu [7] is a combination of the homotopy perturbation method and Laplace transform method that is used to solve various types of linear and nonlinear systems of partial differential equations. The modified homotopy perturbation transform method (MHPTM) by Khan and Smarda [8] is based on the application of Laplace transform to solve the third-order boundary layer equation on semi-infinite domain. Nazari-Golshan et al. [9] developed a modified homotopy perturbation Fourier transform method for nonlinear and singular Lane-Emden equations.

The goal of the present work is to present an algorithm base on MHPTM to handle the problem of transient flow of third grade fluid in a channel with oscillating upper wall, study the fluid behaviour in particular, and examine the effect of viscoelastic parameter on velocity field.

2. New Analytical Algorithm

Consider the following differential equation:

Usually the operators can be decomposed into two parts, a linear part and a nonlinear part :

We construct a homotopy as follows:

where is an embedding parameter. Taking the Laplace transform of both sides of (2) we obtain

Considering the linear operator in (3), the concept of the homotopy perturbation method with embedding parameter is used to generate a series expansion for as follows:

and for the nonlinear operator in (3), follow the concept of He’s polynomial, , as follows:

where He’s polynomials (Ghorbani [6]), , are defined as

Substituting (5) and (6) into (4) we obtain

The first few components of He’s polynomials, for example, are given by

However, (8) can be rewritten in the form

Using (10) we introduce the recursive relation

such that

The recursive equation deduced from (12) can be rewritten as

3. Model of the Problem

Consider a third grade viscoelastic fluid which is unsteady flows between two porous infinite vertical and parallel plane walls. The distance between the walls, that is, the channel width, is . The lower plate is stationary and the upper plate is oscillating with periodic velocity . The lower and the upper plates are, accordingly, located in the planes and of an orthogonal coordinate system with -axis in the direction of flow. The -axis is orthogonal to the channel walls, and the origin of the axes is such that the positions of the channel walls are and , respectively.

The fluid velocity vector is assumed to be parallel to the -axis, so that only the -component of the velocity vector does not vanish but the transpiration cross-flow velocity remains constant, where is the velocity of blowing and is the velocity of suction. Initially, both the channel walls and the fluid are at rest. The external pressure gradient is zero and the fluid velocity is described by the governing equation:

where is the fluid density, is the coefficient of viscosity, is the viscoelastic parameter for a second grade fluid, and is the viscoelastic parameter for a third grade fluid.

The initial and boundary conditions are

where is the amplitude of wall oscillations, is the frequency of the wall velocity, and is the imaginary unit. Using the wall velocity given in the expression (17), the cosine and sine oscillations can be obtained by taking the real and imaginary parts of the velocity field .

Consider the following set of nondimensional variables:

We obtain the nondimensional initial-boundary values problem (dropping the notation)

4. Solution Technique

4.1. Application of New Algorithm

To solve the problem formulated in the previous section, we apply the new algorithm formulated in Section 2.

By applying the Laplace transform with respect to time of (19)–(22) we get the following problem:

where is the Laplace transform of the function . Substituting the recursive (12) into (23) leads to the following equation:

The recursive equation deduced from (26) can be written as follows:

The solutions of the recursive (27) can be compactly written as

Using the Maple symbolic code, the inverse Laplace transform of (28) is

Consequently, the first-order approximate analytical solution of (19) is given by

4.2. Application of Homotopy Perturbation Method (HPM)

Rewrite (19) as

Integrate (31) with respect to over the interval  as follows:

According to the HPM technique by He [5], we construct the following homotopy which satisfies the following relation:

where , , and is an initial approximation to the transient solution .

Taking as small parameter we assume a power series solution of (33) in the form

where are unknown functions of , . Now letting , (34) yields the approximate solution of in the following form:

We now substitute (34) into (33) and the initial and boundary conditions (20)–(22) and equate the coefficients of like powers of to obtain first-order problem:

We can now solve these problems to find and :

The first-order approximate solution of (19) by HPM method is

4.3. Application of Optimal Homotopy Asymptotic Method (OHAM)

By means of the OHAM proposed by Marinca and Herişanu [10], we construct an optimal homotopy which satisfies the following relation:

We have a great freedom to choose the auxiliary function as

where , , and are functions depending on the variable , .

Let us consider the solutions of (39) in the form

Substituting (41) into (39) and equating the coefficients of like powers of we obtain the governing equations of ; that is,

The first-order approximate solution of the problem is

where the zeroth-order and first-order problems from (42) are

The solution of (44) can be solved using widely symbolic computational Maple software. The values of the constants , , and are obtained using collocation method.

5. Analysis of Result

In Figure 1 and Table 1, we show that the transient approximate solutions obtained using newly technique are in better agreement with the OHAM method, as compared to the HPM solutions, for small values of time   and different values of all the non-Newtonian parameters and constants. This proves that the accuracy of the solution obtained by the new method is significant and more accurate for small values of the time . However, in Figure 2 and Table 2 for large values of time , the transient approximate solutions obtained using newly method become divergent with the approximate solutions obtained using HPM and OHAM which proves that the accuracy of the new method fails for growing values of time .

Table 1: Numerical results show the comparison of new algorithm results with the results obtained from HPM and OHAM for the case of small value of time  , when ,  , , , , , and .

New algorithm
HPM method
OHAM method
Absolute error
Absolute error

HPM method
OHAM method

Table 2: Numerical results show the comparison of new algorithm results with the results obtained from HPM and OHAM for the case of large value of time  , when ,  , , , , , and .

New method
HPM method
OHAM method
Absolute error
Absolute error

HPM method
OHAM method

0.1622669
0.1584984

0.3224548
0.2768760

0.4771476
0.3532050

0.6200831
0.4017600

0.7403246
0.4384540

0.8200037
0.4726490

0.8316530
0.4998180

0.7355870
0.4908880

0.4788530
0.3737000

0
0

Figure 1: Comparison of the present results with HPM and OHAM for small value of time .

Figure 2: Comparison of the present results with HPM and OHAM for large value of time .

To see the physical impact of the oscillating walls on the third grade fluid on the flow field, we have plotted the graphs for velocity profiles given by (30). Figure 3 is plotted to show the variation for different amplitude of wall oscillations in cases of cosine and sine oscillation; it is noted that the results are physically satisfied with different values for both cosine and sine excitation of the upper wall. Figure 4 shows the effect of third grade viscoelastic parameter with blowing case on the velocity profile for small values of time ; it is clearly seen that, by increasing , the velocity increases across the channel, and this increase is rapid near the upper wall. This abrupt change in velocity near the upper wall is due to oscillatory nature of the wall boundary which generates depressive harmonic waves into the velocity field. From this figure we can also compare the velocity field of third grade fluid with corresponding velocity field of second grade fluid . For both cosine and sine oscillation, the third grade fluid flows faster than second grade fluid across the channel. It is noted from Figure 5 that, as we increase the value of the time , the influence of the third grade parameter on the fluid motion is not significant as compared in Figure 4 for small values of time .

(a)

(b)

Figure 3: Plot of velocity field for varying values of for small time and , , and .

(a)

(b)

Figure 4: Plot of velocity field for varying values of for small value of time and , and .

(a)

(b)

Figure 5: Plot of velocity field for varying values of for large value of time and , and .

In Figure 6 the wall stress is plotted against the dimensionless space coordinate . For cosine oscillations, it is observed that the wall stress increases near the lower plate, whereas it decreases near the upper plate. This is because of the simultaneous suction and blowing phenomena at the lower plate and upper plate, respectively. The effects get reversed in the case of sine oscillations.

(a)

(b)

Figure 6: Effect of blowing on wall stress when , , , and are fixed.

6. Conclusion Remark

In this paper, a new technique is proposed to obtain analytical solutions of the transient flow of viscoelastic third grade fluid. We applied a new technique to our problem, which is also solved using HPM and OHAM. We obtained an explicit analytic solution of the 2D laminar transient third grade flow over a vertical channel with oscillation motions on a transpiration wall. This explicit analytic solution is valid in the whole transient region . The results obtained using HPM and OHAM are compared with the solutions of the new technique. We observed that we get accurate results using new method, by comparing with OHAM results.

The result has shown that the influence of the third grade parameter on the fluid motion is significant only for small values of time .

This approach seems to be useful and can be used to obtain other analytical solutions for other moving boundary layer transient equations in fluid mechanics Fan et al. [11] and Noghrehabadi et al. [12].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Associate Professor Dr. Rozaini Roslan and Dr. Mahathir Bin Mohamad have made significant contributions in the revised version of the paper. They assisted immensely in the comparison between HPM and OHAM. Their contributions have helped in making positive change in the paper.

References

K. Fakhar, Z. Xu, and C. Yi, “Exact solutions of a third grade fluid flow on a porous plate,” Applied Mathematics and Computation, vol. 202, no. 1, pp. 376–382, 2008.

M. Danish and S. Kumar, “Exact analytical solutions for the Poiseuille and Couette-Poiseuille flow of third grade fluid between parallel plates,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 3, pp. 1089–1097, 2012.

M. Abdulhameed, I. Khan, D. Vieru, and S. Shafie, “Exact solutions for unsteady flow of a second grade fluid generated by an oscillating wall with transpiration,” Applied Mathematics and Mechanics, vol. 35, no. 7, pp. 821–830, 2014.

M. Ayub, A. Rasheed, and T. Hayat, “Exact flow of a third grade fluid past a porous plate using homotopy analysis method,” International Journal of Engineering Science, vol. 41, no. 18, pp. 2091–2103, 2003.

J. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005.

A. Ghorbani, “Beyond Adomian polynomials: he polynomials,” Chaos, Solitons and Fractals, vol. 39, no. 3, pp. 1486–1492, 2009.

Y. Khan and Q. Wu, “Homotopy perturbation transform method for nonlinear equations using He's polynomials,” Computers & Mathematics with Applications, vol. 61, no. 8, pp. 1963–1967, 2011.

Y. Khan and Z. Smarda, “A novel computing approach for third order boundary layer equation,” Sains Malaysiana, vol. 41, no. 11, pp. 1489–1493, 2012.

A. Nazari-Golshan, S. S. Nourazar, H. Ghafoori-Fard, A. Yildirim, and A. Campo, “A modified homotopy perturbation method coupled with the Fourier transform for nonlinear and singular Lane-Emden equations,” Applied Mathematics Letters, vol. 26, no. 10, pp. 1018–1025, 2013.

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.

T. Fan, H. Xu, and I. Pop, “Mixed convection heat transfer in horizontal channel filled with nanofluids,” Applied Mathematics and Mechanics, vol. 34, no. 3, pp. 339–350, 2013.

A. Noghrehabadi, A. Behseresht, and M. Ghalambaz, “Natural convection of nanofluid over vertical plate embedded in porous medium: prescribed surface heat flux,” Applied Mathematics and Mechanics, vol. 34, no. 6, pp. 669–686, 2013.