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

A new analytical algorithm based onmodified homotopy perturbation transformmethod 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.


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

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 ∈ [0,1] 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

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 2ℎ. 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 V = [ ( , ), −V ] 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 V remains constant, where V < 0 is the velocity of blowing and V > 0 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 u( , ) = ( , )i is described by the governing equation: where is the fluid density, is the coefficient of viscosity, 1 is the viscoelastic parameter for a second grade fluid, and 3 is the viscoelastic parameter for a third grade fluid.
Journal of Computational Engineering 3 The initial and boundary conditions are where 0 is the amplitude of wall oscillations, > 0 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)

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 ( , ) = ∫ ∞ 0 ( , ) − 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 Journal of Computational Engineering Using the Maple symbolic code, the inverse Laplace transform of (28) is Consequently, the first-order approximate analytical solution of (19) is given by Integrate (31) with respect to over the interval [0, ] as follows: According to the HPM technique by He [5], we construct the following homotopy V( , ; ) : R × [0, 1] → R which satisfies the following relation: where , and 0 is an initial approximation to the transient solution ( , ).
Taking as small parameter we assume a power series solution of (33) in the form where V ( , ) are unknown functions of , . Now letting → 1, (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 V 0 and V 1 : The first-order approximate solution of (19) by HPM method is (38)

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

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 0 < < 1.5. 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].