`Mathematical Problems in EngineeringVolumeΒ 2011, Article IDΒ 349803, 10 pageshttp://dx.doi.org/10.1155/2011/349803`
Research Article

## An Axisymmetric Squeezing Fluid Flow between the Two Infinite Parallel Plates in a Porous Medium Channel

2Islamia College Peshawar (Chartered University), Khyber Pakhtunkhwa, Peshawar 25120, Pakistan
3Malakand University, Khyber Pukhtonnkhawa, Chakdara 1800, Pakistan

Received 5 October 2010; Revised 14 December 2010; Accepted 22 February 2011

Copyright Β© 2011 S. Islam 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

The flow between two large parallel plates approaching each other symmetrically in a porous medium is studied. The Navier-Stokes equations have been transformed into an ordinary nonlinear differential equation using a transformation . Solution to the problem is obtained by using differential transform method (DTM) by varying different Newtonian fluid parameters and permeability of the porous medium. Result for the stream function is presented. Validity of the solutions is confirmed by evaluating the residual in each case, and the proposed scheme gives excellent and reliable results. The influence of different parameters on the flow has been discussed and presented through graphs.

#### 1. Introduction

The study through a porous medium is an interesting and hot issue in these days, especially, with the introduction of the modified Darcyβs Law [1], in contrast to the classical Darcyβs Law [2] and the Brinkman model [3]. The flow through a porous media has wide spread applications in engineering and science, such as ground water hydrology, petroleum engineering, reservoir engineering, chemical engineering, chemical reactors to agriculture irrigation and drainage and the recovery of crude oil from the pores of the reservoir rocks [4β9].

Squeezing flows are common in moulding, food industry, and chemical engineering, and they have, therefore, been studied for a long time as researchers and scientists have sought to optimize processing operations to produce improved components. These flows are also interesting from rheological perspective [10]. The fundamental analysis of the squeezing flow between two large parallel plates was carried out by Stefan [11] in which he developed an ad hoc asymptotic solution for a thin film of Newtonian fluid. Naduvinamani et al. in [12, 13] studied static and dynamic behavior of squeeze-film lubrication of narrow porous journal bearings with coupled stress fluid and a squeeze film lubrication of a short porous journal bearing with couple stress fluids, respectively.

The determination of squeeze flow characteristics has attracted the attention of several investigators due to its importance in the practical problems of improving the performance of hydraulic machine elements, food industry, chemical engineering, polymer processing, compression, and injection molding [14β20].

Differential transform method (DTM) is one of the numerical methods for differential equations. The basic idea of differential transformation was initially introduced by Zhou [21] in 1986. Its main application therein is to solve both linear and nonlinear initial value problems in electrical circuit analysis. This method constructs the solution in the form of a polynomial. It is different from the traditional higher-order Taylor series method. The DTM is an alternative procedure for getting Taylor series solution of the differential equations. This method reduces the size of computational domain and is easily applicable to a variety of problems. The method was applied to differential algebraic equations (DAEs) of index-1 by Ayaz [22]. Liu and Song [23] analyzed higher index differential algebraic equations using this technique where they showed that the method is effective in case of index-2 DAEs but not suitable for DAEs of index-3. Application of the two-dimensional differential transform method was studied by Ayaz [24] for partial differential equations. Comparison of this method with ADM was done by Hassan [25] to solve PDEs. The same author has used this method to solve higher-order initial value problems where he studied second- and third-order initial value problems to show the efficiency of the method. Islam et al. [26] used this method for the solution of special twelfth-order boundary value problems.

The study reported in this paper considers the axisymmetric squeezing flow of a viscous fluid between two large parallel plates in a porous medium separated by a small distance and the plates approaching each other with a low constant velocity , and the flow can be assumed to be quasisteady. We use DTM to find the approximate solutions of the modeled problem.

#### 2. Basic Equations and Problem Formulation

Employing the modified Darcyβs Law, the two-dimensional flow of an incompressible axisymmetric homogenous fluid in a porous medium in the absence of body force is governed by the following equations [27, 28]: where is the velocity vector, denotes the material time derivative, is the constant density, is the Cauchy stress tensor where in which is the Rivlin-Ericksen tensor, is the pressure, is the viscosity, and is the Darcyβs resistance given by the relation, here, is the permeability and is the effective viscosity of the porous medium and is the Darcian velocity which is related to the fluid velocity by , () being the porosity of the medium. This assumes the same velocity in each pore, and the fluid in the pores is averaged over the volume. In general, the effective viscosity and the fluid viscosity are different. However, at the macrolevel we may take them equal, though this assumption does not hold at the microlevel. Following Naduvinamani et al. [12] and Breugem [29], (2.3) can now be written as Continuity equation (2.1) implies the existence of a stream function such that Substitution of (2.5) together with (2.4) in (2.2) and elimination of pressure from the resulting equations using yield where.

#### 3. Basic Idea of Differential Transform Method (DTM)

If is a given function, its differential transform is defined as The inverse transform of is defined by In actual application, the function is expressed by a finite series Equation (3.2) implies that is negligibly small.

The fundamental operations of the DTM are given in Table 1.

Table 1

#### 4. Analysis of the Method

Consider a fourth-order boundary value problem with the following boundary conditions: where , are given values.

The differential transform of (4.1) is where is the differential transform of .

The transformed boundary conditions (4.2) are given by Using (4.3) and (4.4), values of , are obtained which give the following series solution up to ,

#### 5. Our Problem

Here, we consider an incompressible Newtonian fluid, squeezed between two large planar, parallel plates which are separated by a small distance and moving towards each other with velocity . For small values of the velocity , the gap distance between the plates changes slowly with time , so that it can be taken as constant, the flow is steady [30, 31]: with the following boundary conditions: where , .

The transformed boundary conditions (5.2) are where and will be determined later.

The differential transform of (5.1) is Using (5.3) and (5.4), we obtain the following values of, : We fix and . The permeability and the fluid parameters , are taken without units as they appear in the value of and which have been nondimensionalized. For these values, the unknowns and are determined by the following system: We find that , .

Now, using the inverse differential transform, the following approximate solution of is obtained: Residual of the solution is In Table 2, numerical results of the DTM solution are compared with the results of Mathematica NDSolve solution. Residuals of both the solutions are given for comparison. List Interpolation is used for the construction of approximating polynomial of the numerical solution by NDSolve.

Table 2

In Table 3, residuals of the DTM solutions have been evaluated for various values of , , , keeping and fixed.

Table 3

In Table 4, we fix , , and vary , . The reliability of the solution can be seen by looking into the residual at different mesh points.

Table 4

#### 6. Conclusion

In this paper, we studied squeezing flow in a porous medium between two parallel plates. Flow pattern for various parameters of Newtonian fluids are derived. When the parameters , , and are fixed as 1, and velocity of the plates is varied, different velocity profiles are obtained. It is clearly visible that when the velocity of the plates is increased the fluid velocity is increased (Figure 1). For , the variation in flow pattern is negligible but for higher values of it does matter as seen in Table 2, and Figure 2. It is also noticed that at lower velocity of the plates the flow rate is also lower. The results obtained by the application of DTM are reliable with high accuracy. Applicability of the method is simple and needs no restrictions of large and small parameters. It avoids the difficulties and massive computational work encountered in other numerical techniques such linearization, discretization and perturbation. Moreover this method has superiority over the Adomian decomposition method as it does not require calculation of Adomian polynomials in case of nonlinearity. Furthermore, the solution obtained by this method converges rapidly to analytical solution in case of integrable system.

Figure 1: Velocity profiles of the squeezing flow in a porous medium when the plates velocity is varied.
Figure 2: Velocity profiles of fluid flow when permeability , is increased for a fixed.

#### Acknowledgment

The authors would like to thank the anonymous reviewers for their helpful and encouraging comments on earlier draft of this paper.

#### References

1. W. Tan and T. Masuoka, βStokes' first problem for a second grade fluid in a porous half-space with heated boundary,β International Journal of Non-Linear Mechanics, vol. 40, no. 4, pp. 515β522, 2005.
2. J. Bear, Dynamics of Fluids in Porous Media, Elsevier, New York, NY, USA, 1972.
3. H. C. Brinkman, βA calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles,β Applied Scientific Research, vol. 1, no. 1, pp. 27β34, 1947.
4. M. H. Hamdan, βAn alternative approach to exact solutions of a special class of Navier-Stokes flows,β Applied Mathematics and Computation, vol. 93, no. 1, pp. 83β90, 1998.
5. M. H. Hamdan and F. M. Allan, βA note on the generalized Beltrami flow through porous media,β International Journal of Pure and Applied Mathematics, vol. 27, no. 4, pp. 491β500, 2006.
6. K. Vafai and C. L. Tien, βBoundary and inertia effects on flow and heat transfer in porous media,β International Journal of Heat and Mass Transfer, vol. 24, no. 2, pp. 195β203, 1981.
7. S. Islam, M. R. Mohyuddin, and C. Y. Zhou, βFew exact solutions of non-Newtonian fluid in porous medium with hall effect,β Journal of Porous Media, vol. 11, no. 7, pp. 669β680, 2008.
8. S. Islam and C. Y. Zhou, βCertain inverse solutions of a second-grade magnetohydrodynamic aligned fluid flow in a porous medium,β Journal of Porous Media, vol. 10, no. 4, pp. 401β408, 2007.
9. K. Vafai, Hand Book of Porous Media, CRC Press, Taylor & Francis Group, Boca Raton, Fla, USA, 2nd edition, 2005.
10. A. S. Burbidge and C. Servais, βSqueeze flows of apparently lubricated thin films,β Journal of Non-Newtonian Fluid Mechanics, vol. 124, no. 1β3, pp. 115β127, 2004.
11. M. J. Stefan, βVersuch Uber die scheinbare adhesion,β Akademie der Wissenschaften in Wien. Mathematisch-Naturwissenschaftliche, vol. 69, p. 713, 1874.
12. N. B. Naduvinamani, P. S. Hiremath, and G. Gurubasavaraj, βStatic and dynamic behaviour of squeeze-film lubrication of narrow porous journal bearings with coupled stress fluid,β Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, vol. 215, no. 1, pp. 45β62, 2001.
13. N. B. Naduvinamani, P. S. Hiremath, and G. Gurubasavaraj, βSqueeze film lubrication of a short porous journal bearing with couple stress fluids,β Tribology International, vol. 34, no. 11, pp. 739β747, 2001.
14. Q. K. Ghori, M. Ahmed, and A. M. Siddiqui, βApplication of homotopy perturbation method to squeezing flow of a newtonian fluid,β International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 179β184, 2007.
15. X. J. Ran, Q. Y. Zhu, and Y. Li, βAn explicit series solution of the squeezing flow between two infinite plates by means of the homotopy analysis method,β Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 1, pp. 119β132, 2009.
16. R. J. Grimm, βSqueezing flows of Newtonian liquid films an analysis including fluid inertia,β Applied Scientific Research, vol. 32, no. 2, pp. 149β166, 1976.
17. W. F. Hughes and R. A. Elco, βMagnetohydrodynamic lubrication flow between parallel rotating disks,β Journal of Fluid Mechanics, vol. 13, pp. 21β32, 1962.
18. S. Kamiyama, βInertia Effects in MHD hydrostatic thrust bearing,β ASME Transactions Journal of Lubrication Technology, vol. 91, pp. 589β596, 1969.
19. E. A. Hamza, βThe Magnetohydrodynamic squeeze film,β Journal of Tribology, vol. 110, no. 2, pp. 375β377, 1988.
20. S. Bhattacharyya and A. Pal, βUnsteady MHD squeezing flow between two parallel rotating discs,β Mechanics Research Communications, vol. 24, no. 6, pp. 615β623, 1997.
21. J. K. Zhou, Deferential Transformation and Its Application for Electrical Circuits, Huazhong University Press, Wuhan, China, 1986.
22. F. Ayaz, βApplications of differential transform method to differential-algebraic equations,β Applied Mathematics and Computation, vol. 152, no. 3, pp. 649β657, 2004.
23. H. Liu and Y. Song, βDifferential transform method applied to high index differential-algebraic equations,β Applied Mathematics and Computation, vol. 184, no. 2, pp. 748β753, 2007.
24. F. Ayaz, βOn the two-dimensional differential transform method,β Applied Mathematics and Computation, vol. 143, pp. 361β374, 2003.
25. I. H. A. Hassan, βComparison deferential transformation technique with Adomian decomposition method for linear and nonlinear initial value problems,β Chaos, Solution and Fractals, vol. 36, no. 1, pp. 53β65, 2007.
26. S. U. Islam, S. Haq, and J. Ali, βNumerical solution of special 12th-order boundary value problems using differential transform method,β Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1132β1138, 2009.
27. M. T. Kamel and H. M. Hamdan, βRiabounchisky flow through porous media,β International Journal of Pure and Applied Mathematics, vol. 27, no. 1, pp. 113β126, 2006.
28. N. Merabet, H. Siyyam, and M. H. Hamdan, βAnalytical approach to the Darcy-Lapwood-Brinkman equation,β Applied Mathematics and Computation, vol. 196, no. 2, pp. 679β685, 2008.
29. W. -P. Breugem, βThe effective viscosity of a channel-type porous medium,β Physics of Fluids, vol. 19, no. 10, 2007.
30. T. C. Papanastasiou, G. C. Georgio, and A. N. Alexandrou, Viscous Fluid Flow, CRC Press, New York, NY, USA, 1994.
31. M. Idrees, S. Islam, S. Haq, and S. Islam, βApplication of the optimal homotopy asymptotic method to squeezing flow,β Computers and Mathematics with Applications, vol. 59, no. 12, pp. 3858β3866, 2010.