Flow behavior of unsteady incompressible Newtonian fluid flow between two parallel plates via homotopy analysis method

Abstract

An analysis has been performed to study the problem of the flow of incompressible Newtonian fluid between two parallel plates where the upper plate is impermeable and can move up or down and the lower one is fixed and has a porous surface. The governing equations for this problem are reduced to an ordinary form and is solved using Homotopy Analysis Method (HAM) and numerically by fourth order Runge-Kutta technique. Also, Velocity fields have been computed and shown graphically for various values of physical parameters. As an important outcome, HAM is able to solve a large class of nonlinear problems effectively, more easily and accurately; and thus it has been widely applicable in engineering and physics.

Keywords:
Homotopy Analysis Method; parallel porous plates; unsteady flow

1 INTRODUCTION

The problem of unsteady time-dependent flow between parallel plates has many crucial applications in science and technology. Among them are hydrodynamic lubrication, aerodynamic heating, polymer technology, petroleum industry and biomechanics. Many researchers have investigated such flows with different geometries and different flow conditions (Sharma and Singh, 2008;Ishizawa, 1966;Hamza, 1999;Rashidi et al., 2008;Hasanzadeh et al., 2013).

In recent decades many attempts have been made to develop analytical methods for solving such nonlinear equations. One of them is the perturbation method (Nayfeh, 2000), which is strongly dependent on a so called small parameter to be defined according to the physics of the problem. Since these equations cannot be solved via the conventional analytical techniques, recent attempts have been focused on constructing an analytical solution for these equations using the advanced developed methods such as Adomian's Decomposition Method (ADM) (Sheikholeslami, 2012), Homotopy Perturbation Method (HPM) (He, 2006;Rahimi et al., 2012), Variational x'Iteration Method (VIM) (He, 2007), Differential Transformation Method (Ganji and Azimi, 2013), Homotopy Aanalysis Method (HAM) (Liao, 2003;Liao, 2012;Abbasi et al., 2014;Ganji et al., 2014) and Least Square Method (LSM) (Hatami et al., 2013;2014).

Liao introduced the basic idea of Homotopy in topology to propose a general analytical method for nonlinear problems, namely the Homotopy Analysis Method (Hatami and Ganji, 2013;2014), that does not need any small parameter. This method has been successfully applied to solve many types of nonlinear problems (Hatami et al., 2014;Hasanzadeh et al., 2013;Aziz, 2006).

The aim of this study is to investigate, the effect of physical parameters on Flow behavior of unsteady incompressible Newtonian fluid flow between two parallel plates. In addition, the convergence of the series solution is also explicitly discussed. Obtaining the analytical solution of the models and comparing with numerical result reveal the capability, effectiveness and convenience of HAM.

2 PROBLEM STATEMENT AND MATHEMATICAL FORMULATION

We consider an incompressible two dimensional flow of Newtonian fluid between two parallel infinite rectangular plates in Cartesian coordinates. These two plates are placed a distance a(t) at apart from each other and t denotes time.

We also consider that the upper plate which is at y = a(t) at = is moving toward the lower plate with velocity=(t)and the lower porous plate which is at y = o is fixed. The equations of motion for this flow are (Hasanzadeh et al., 2013):

where ρ denotes density, P denotes pressure and ν denotes kinematic viscosity. For boundary conditions we have:

As we can see for the axial velocity we have no slip boundary condition. Also as the distance between two plates varies with(t) and the lower plate is stationary we can write for the upper plate that v =(t). For the lower plate we insert a constant parameter like K which is a candidate for the strength of the suction or injection. Here K > 0corresponds to suction and K < 0 corresponds to injection. An assumption made by (Singh et al., 1990) proposes that,

Substituting these changes in Eq. (1) to Eq. (3), Equation of continuity will be satisfied automatically and the Navier-Stokes equations of motion reduce to:

Where

Now differentiating Eq. (7) with respect to η gives,

Similarity solution exists only when R and S are constants. Now integrating the first equation of (9) we can find the distance between two plates as,

From this equation, we can say that when K > 0 the upper plate moves away from the lower plate and when K < 0, it moves toward it. Combining Eq. (9) with Eq. (11) we see that S =-1. Equation (10) will then be written as,

With the boundary conditions,

Eq. (12) with boundary conditions (13) will be solved by analytical method.

3 IMPLEMENTATION OF THE HOMOTOPY ANALYSIS METHOD

For HAM solutions, we choose the initial guess and auxiliary linear operator in the following form:

Where c_i (i = 1,2,3,4) are constants. Let P ∈ [0,1] denotes the embedding parameter and h indicates non-zero auxiliary parameters. We then construct the following equations:

Zeroth -order deformation equations

For p = 0 and p = 1 we have

When p increases from 0 to 1 then F(η;p) varies from f₀(η) to f(η). By Taylor's theorem and using equation (20), F(η;p) can be expanded in a power series of p as follows:

In which h is chosen in such a way that this series is convergent at p = 1, therefore we have through equation (26) that

mth -order deformation equations

Now we determine the convergency of the result, the differential equation, and the auxiliary function according to the solution expression. So let us assume:

We have found the answer by maple analytic solution device. The first deformation is presented below

The solutions f(η) were too long to be mentioned here, therefore, they are shown graphically.

4 CONVERGENCE OF THE HAM SOLUTION

As pointed out by (Liao, 2003,2012), the convergence region and rate of solution series can be adjusted and controlled by means of the auxiliary parameter h. To influence of h on the convergence of solution, we plot the so-called h curve of f'''(0), as shown inFigures 2and3. The solutions converge for h values which are corresponding to the horizontal line segment in h curve.

In order to investigate the range of admissible values of the auxiliary parameter h for various quantities of R and K, the curves of h were derived 6th-order approximations.Figures 2and3shows a typical h curve for f'''(0) which presentation admissible values for auxiliary parameter h. In our case study, it is easy to discover that h = 1 is suitable value which is used for values of -2.0 < K < 2.0 (M ≥ 0) and -10 < R < 10 (0 ≤ T ≤ 6).

5 RESULTS

In the present study HAM method is applied to obtain an explicit analytic solution of the flow of incompressible Newtonian fluid between two parallel plates (Figure 1). First, a comparison between the applied methods, obtained by the numerical method and HAM for different values of active parameters is shown inFigures 4till 6. The numerical solution is performed using the algebra package Maple 16.0, to solve the present case. The package uses a fourth order Runge-Kutta procedure for solving nonlinear boundary value (B-V) problem by (Aziz, 2006). Validity of HAM is shown inTable 1. In these tables, the %Error is defined as:

Figure 1
Schematic diagram and the coordinate system for the considered flow.

Figure 2
The h validity for R = 0.1and different value of K.

Figure 3
The h validity for K = 0.9 and different value of R.

Table 1
The results of HAM and Numerical methods for f(η) and f'(η) for K = -0.3 and Re = 4.

η	f(η)		f'(η)
	HAM	NUM	Error	HAM	NUM	Error
0.00	-0.300000000	-0.300000000	0.00E+00	0.000000000	0.000000000	0.00E+00
0.05	-0.291354676	-0.291354676	5.00E-10	0.345993417	0.345993396	2.13E-08
0.10	-0.265475699	-0.265475701	1.90E-09	0.687476197	0.687476163	3.44E-08
0.15	-0.222869707	-0.222869711	3.80E-09	1.013210723	1.013210683	4.00E-08
0.20	-0.164599515	-0.164599521	5.50E-09	1.312281007	1.312280969	3.80E-08
0.25	-0.092254286	-0.092254293	7.35E-09	1.574636006	1.574635976	3.00E-08
0.30	-0.007893760	-0.007893768	8.82E-09	1.791579558	1.791579538	2.00E-08
0.35	0.086030154	0.086030144	9.96E-09	1.956181991	1.956181985	6.00E-09
0.40	0.186770002	0.186769991	1.07E-08	2.063594851	2.063594855	4.00E-09
0.45	0.291392456	0.291392445	1.10E-08	2.111256426	2.111256446	2.00E-08
0.50	0.396895398	0.396895387	1.07E-08	2.098982036	2.098982067	3.10E-08
0.55	0.500326605	0.500326595	9.70E-09	2.028938543	2.028938584	4.10E-08
0.60	0.598898367	0.598898359	8.00E-09	1.905508152	1.905508196	4.40E-08
0.65	0.690092719	0.690092714	4.70E-09	1.735051163	1.735051202	3.90E-08
0.70	0.771752572	0.771752569	3.60E-09	1.525581730	1.525581765	3.50E-08
0.75	0.842154803	0.842154802	1.40E-09	1.286374302	1.286374324	2.20E-08
0.80	0.900062378	0.900062379	1.00E-09	1.027521158	1.027521170	1.20E-08
0.85	0.944753585	0.944753585	2.00E-10	0.759463805	0.759463805	1.00E-10
0.90	0.976027672	0.976027670	1.60E-09	0.492521856	0.492521847	8.80E-09
0.95	0.994187374	0.994187376	2.10E-09	0.236443433	0.236443428	5.00E-09
1.00	1.000000000	1.000000000	1.00E-09	4.00000E-09	0.000000000	4.00E-09

The results are proved to be precise and accurate in solving a wide range of mathematical and engineering problems especially Fluid mechanic cases. This accuracy gives high confidence to us about validity of this problem and reveals an excellent agreement of engineering accuracy. This investigation is completed by depicting the effects of some important parameters to evaluate how these parameters influence on this fluid.

From a physical point of view,Figures 4to6are prepared in order to see the effects of the K and R flow parameters on the velocity distribution.Figure 4is the graphical representation on the velocity distribution in both cases of approaching (R < 0) and receding (R > 0) plates for different values of the K number.

Figure 4
Dimensionless velocities predicted by HAM and numerical method (NUM) in different R number when K = 0.0, h = -1.0.

Figure 5
Dimensionless velocities predicted by HAM and numerical method (NUM) in different R number when K = -0.9, h = -1.0.

Figure 6
Dimensionless velocities predicted by HAM and numerical method (NUM) in different R number when K = 0.9, h = -1.0.

According to the obtained solution, It can be seen that when K = 0, in the case of approaching plates R < 0, increasing the magnitude of R will damp the maximum value of f'. By contrast, when the plates recede each other everything is reversed e.g. we have an increase in the middle region.

Moreover,Figures 5and6have been prepared for the variations of flow number R on the distribution of velocity for K < 0 and K > 0 respectively. As we seen that, in the case of K < 0, near the fixed plate the increase in R will increase the flow but near the moving plate this reverses. In addition,Figure 6shows that we have the same trend like the injection case but in the opposite direction. The increase of flow near the wall region is completely obvious but the flow near the upper plate will not change sensitively with the increase of R.

6 CONCLUSIONS

In this investigation, the analytical approach called Homotopy Analysis Method (HAM) has been successfully applied to find the most accurate analytical solution for the velocity distributions of unsteady incompressible Newtonian fluid flow between two parallel plates. Furthermore, the obtained solutions by proposed methods have been compared with the direct numerical solutions generated by the symbolic algebra package Maple 16. In addition, the effects of different physical parameters, such as K, and R on the velocity profiles of the problem have been investigated. The following main points can be concluded from the present study:

• •

  The comparison shows that the HAM solutions is highly accurate and provide the rapid achievement to compute the flow velocities. Also according to the previous publications this methods is a powerful technique for finding analytical solutions in science and engineering problems.

• •

  The results show that the reversal flow will take place near the wall region above a critical value of R in both case of injection and suction.

References

• Abbasi, M., Ganji, D.D., Rahimi Petroudi, I., Khaki, M., (2014). Comparative analysis of MHD Boundary-Layer flow of viscoelastic fluid in permeable channel with slip Boundaries by using HAM, VIM, HPM. Walailak Journal for science and Technology 11(7): 551-567.
• Aziz, A., (2006). Heat Conduction with Maple. R.T. Edwards, Philadelphia (PA).
• Ganji, D.D., Azimi, M., (2013). Application of DTM on MHD Jeffery Hamel problem with nanoparticles. U.P.B. Sci. Bull., Series A, 75(1).
• Ganji, D.D., Abbasi, M., Rahimi, J., Gholami, M., Rahimipetroudi, I., (2014). On the MHD squeezeflow between two parallel disks with suction or injection via HAM and HPM. Frontiers of Mechanical Engineering. DOI 10.1007/s11465-014-0303-0.
  » https://doi.org/10.1007/s11465-014-0303-0
• Hamza, E.A., (1999). Suction and injection effects on a similar flow between parallel plates. J. of Phys. D: Appl. Phys. 32: 656-663.
• Hasanzadeh, A., Ebrahimi, S., Ganji, D.D., (2013). Investigation on the problem of unsteady flow between parallel porous plates using Homotopy perturbation method. Int. J. of Appl. Math and Mech. 9(18): 96-115.
• Hatami, M., Ganji, D.D., (2013). Thermal performance of circular convective-radiative porous fins with different section shapes and materials. Energy Conversion and Management 76: 185-193.
• Hatami, M., Hasanpour, A., Ganji, D.D., (2013). Heat transfer study through porous fins (Si3N4and AL) with temperature-dependent heat generation. Energy Conversion and Management 74: 9-16.
• Hatami, M., Ganji, D.D., (2014). Thermal and flow analysis of microchannel heat sink (MCHS) cooled by Cu-water nanofluid using porous media approach and least square method. Energy Conversion and Management 78: 347-358.
• Hatami, M., Mehdizadeh Ahangar, G.H.R., Ganji, D.D., Boubaker, K., (2014). Refrigeration efficiency analysis for fully wet semi-spherical porous fins. Energy Conversion and Management 84: 533-540.
• He, J.H., (2006). New Interpretation of homotopy-perturbation method. International Journal of Modern Physics B 20: 2561-2568.
• He, J.H., (2007). Variational iteration method- some recent results and new interpretations. Journal of Computational and Applied Mathematics 207 (1): 3-17.
• Ishizawa, S., (1966). The unsteady flow between two parallel disks with arbitrary varying gap width. Bulletin of JSME 35: 533-550.
• Liao, S.J., (2003). Beyond Perturbation: Introduction to the homotopy analysis method. Chapman & Hall/CRC Press, Boca Raton.
• Liao, S.J., (2012). Homotopy analysis method in nonlinear differential equations. Springer & Higher Education Press, Heidelberg.
• Nayfeh, A.H., (2000). Perturbation methods. Wiley, New York.
• Rahimi Petroudi, I., Ganji, D.D., Shotorban, A.B., Khazayi Nejad, M., Rahimi, E., Rohollahtabar, R., Taherinia, F., (2012). Semi analytical method for solving nonlinear equation arising of natural convection porous fin. Thermal Science 16(5): 1303-1308.
• Rashidi, M.M., Shahmohamadi, H., Dinarvand, S., (2008). Analytic approximate solutions for unsteady two-dimensional and axisymmetric squeezing flowsbetween parallel plates. Math Probl Eng, article ID 935095. DOI:10.1155/2008/935095.
  » https://doi.org/10.1155/2008/935095
• Sharma, R.P., Singh, G., (2008). Unsteady MHD free convection flow and heat transfer along a vertical porous plate with variable suction and internal heat generation. International Journal of Applied Mathematics and Mechanics 4(5): 1-8.
• Sheikholeslami, M., Ganji, D.D., Ashorynejad, H.R., Rokni, H.B., (2012). Analytical investigation of Jeffery Hamel flow with High Magnetic Field and nanoparticle by Adomian Decomposition Method. Applied Mathematics and Mechanics 33 (11): 25-36.
• Singh, P., Radhakrishnan, V., Narayan, K.A., (1990). Squeezing flow between parallel plates. Ingenieur-Archiv 60(4): 274-281.

Publication Dates

• Publication in this collection
  Oct 2015

History

• Received
  28 Dec 2014
• Accepted
  04 Feb 2015