Research Article | Open Access

Volume 2008 |Article ID 935095 | 13 pages | https://doi.org/10.1155/2008/935095

# Analytic Approximate Solutions for Unsteady Two-Dimensional and Axisymmetric Squeezing Flows between Parallel Plates

Accepted01 Dec 2008
Published25 Feb 2009

#### Abstract

The flow of a viscous incompressible fluid between two parallel plates due to the normal motion of the plates is investigated. The unsteady Navier-Stokes equations are reduced to a nonlinear fourth-order differential equation by using similarity solutions. Homotopy analysis method (HAM) is used to solve this nonlinear equation analytically. The convergence of the obtained series solution is carefully analyzed. The validity of our solutions is verified by the numerical results obtained by fourth-order Runge-Kutta.

#### 1. Introduction

The problem of unsteady squeezing of a viscous incompressible fluid between two parallel plates in motion normal to their own surfaces independent of each other and arbitrary with respect to time is a fundamental type of unsteady flow which is met frequently in many hydrodynamical machines and apparatuses. Some practical examples of squeezing flow include polymer processing, compression, and injection molding. In addition, the lubrication system can also be modeled by squeezing flows. Stefan [1] published a classical paper on squeezing flow by using lubrication approximation. In 1886, Reynolds [2] obtained a solution for elliptic plates, and Archibald [3] studied this problem for rectangular plates. The theoretical and experimental studies of squeezing flows have been conducted by many researchers [4, 414]. Earlier studies of squeezing flow are based on Reynolds equation. The inadequacy of Reynolds equation in the analysis of porous thrust bearings and squeeze films involving high velocity has been demonstrated by Jackson [13], Ishizawa [14]. The general study of the problem with full Navier-Stokes equations involves extensive numerical study requiring more computer time and larger memory. However, many of the important features of this problem can be grasped by prescribing the relative velocity of the plates suitably. If the relative normal velocity is proportional to where t is the time and a constant of dimension which characterizes unsteadiness, then the unsteady Navier–Stokes equations admit similarity solution.

In 1992, Liao [15] employed the basic ideas of the homotopy in topology to propose a general analytic method for nonlinear problems, namely, homotopy analysis method (HAM) [1621]. Based on homotopy of topology, the validity of the HAM is independent of whether or not there exist small parameters in the considered equation. Therefore, the HAM can overcome the foregoing restrictions and limitations of perturbation methods [22]. Furthermore, the HAM always provides us with a family of solution expressions in the auxiliary parameter the convergence region, and the rate of each solution might be determined conveniently by the auxiliary parameter The HAM also avoids discretization and provides an efficient numerical solution with high accuracy, minimal calculation, and avoidance of physically unrealistic assumptions. Besides, the HAM is rather general and contains the homotopy perturbation method (HPM) [21], the Adomian decomposition method (ADM) [23], and δ-expansion method. In fact, HPM and ADM are always special cases of HAM when The convergence of HAM solution series is dependent upon three factors, that is, the initial guess, the auxiliary linear operator, and the auxiliary parameter However, as a special case of homotopy analysis method when the convergence of HPM solution series is only dependent upon two factors: the auxiliary linear operator and the initial guess. So, given the initial guess and the auxiliary linear operator, HPM cannot provide other ways to ensure that the solution is convergent. HAM provides us with a family of solution expression in the auxiliary parameter and the solution given by ADM is only one of them.

In recent years, the HAM has been successfully employed to solve many types of nonlinear problems such as the nonlinear equations arising in heat transfer [24], the nonlinear model of diffusion and reaction in porous catalysts [25], the chaotic dynamical systems [26], the nonhomogeneous Blasius problem [27], the generalized three-dimensional MHD flow over a porous stretching sheet [28], the wire coating analysis using MHD Oldroyd 8-constant fluid [29], the axisymmetric flow and heat transfer of a second-grade fluid past a stretching sheet [30], the MHD flow of a second-grade fluid in a porous channel [31], the generalized Couette flow [32], the Glauert-jet problem [33], the Burger and regularized long wave equations [34], the laminar viscous flow in a semiporous channel in the presence of a uniform magnetic field [35], and other problems. All of these successful applications verified the validity, effectiveness, and flexibility of the HAM.

In this paper, we use homotopy analysis method to investigate the problem of unsteady squeezing of a viscous incompressible fluid between two parallel plates. The paper is organized as follows. In Section 2, the mathematical formulation is presented. In Section 3, we extend the application of the HAM to construct the approximate solutions for the governing equation. The convergence of the obtained series solution is carefully analyzed in Section 4. Section 5 contains the results and discussion. The conclusions are summarized in Section 6.

#### 2. Mathematical formulation

Let the position of the two plates be at where is the position at time as shown in Figure 1. We assume that the length 1 (in the two-dimensional case) or the diameter D (in the axisymmetric case) is much larger than the gap width at any time that the end effects can be neglected. When is positive the two plates are squeezed until they touch at When is negative the two plates are separated. Let u, v, and w be the velocity components in the x, y, and z directions, respectively. For two-dimensional flow, Wang introduced the following transforms [36]: where Substituting (2.1) into the unsteady two-dimensional Navier-Stokes equations yields a nonlinear ordinary differential equation in form where (squeeze number) is the nondimensional parameter. The flow is characterized by this parameter. It should be mentioned that is the kinematic viscosity. The boundary conditions are such that on the plates the lateral velocities are zero and the normal velocity is equal to the velocity of the plate, that is, Similarly, Wang’s transforms [36] for axisymmetric flow are Using transforms (2.5), unsteady axisymmetric Navier–Stokes equations reduce to subject to the boundary conditions (2.4).

Consequently, we should solve the nonlinear ordinary differential equation where and subject to boundary conditions (2.4).

#### 3. HAM solution

To investigate the explicit and totally analytic solutions of (2.7) by using HAM, we choose as initial approximation of which satisfies the boundary conditions (2.4). Besides, we select the auxiliary linear operator as It is easy to check that this operator satisfies the following equation: where are arbitrary constants. Based on (2.7), we are led to define the following nonlinear operator: Using these operators, we can construct the so-called zeroth-order deformation equation as where is an embedding parameter and is an auxiliary nonzero parameter. It should be emphasized that one has great freedom to choose the initial guess, the auxiliary linear operator and the auxiliary parameter However, (3.5), the original equation of HAM, is the origin of the mathematical term “homotopy" (parameter is the head letter of “homotopy"). In addition, if (3.5) will always change to original equation of HPM. The boundary conditions for (3.5) are Obviously, when and the above zeroth-order deformation equation has the following solutions: As p increases from 0 to 1, varies from to Now expanding by its Taylor series in terms of p, one would obtain where As pointed out by Liao [19], the convergence of series (3.8) strongly depends upon the auxiliary parameter Assume that is selected such that series (3.8) is convergent at then due to (3.7), the final series solution becomes For the mth-order deformation equation, we differentiate (3.5) m times with respect to p, divide by m!, and then set The resulting deformation equation at the mth-order is with the following boundary conditions where We use the symbolic software MATHEMATICA to solve the system of linear equations (3.11) with the boundary conditions (3.12), and successively obtain Therefore, like (3.10), the analytical solution of the problem can be expressed as an infinite series of the form (see [37])

#### 4. Convergence of HAM solution

The series solution contains the auxiliary parameter The validity of the method is based on such an assumption that series (3.8) converges at It is the auxiliary parameter which ensures that this assumption can be satisfied. In general, by means of the so-called it is straightforward to choose a proper value of which ensures that the series solution is convergent. For the different values of the squeeze number , the obtained by the 15th-order approximation for the axisymmetric and two-dimensional cases are shown in Figures 2 and 3, respectively. From these figures, the valid regions of correspond to the line segments nearly parallel to the horizontal axis. Figures 2 and 3 elucidate that the size of the valid region strongly depends on S. In fact, the interval for admissible values of shrinks toward zero by increasing the squeeze number. As mentioned above, the homotopy analysis method is rather general and always contains the homotopy perturbation method (HPM) [21] and the Adomian decomposition method (ADM) [23] when From Figures 2 and 3, is not valid for the large values of S.

#### 5. Results and discussion

Our main concern is the various values of and These quantities describe the flow behavior. For several values of , the function obtained by the different orders of approximation for the axisymmetric and two-dimensional cases are compared with the numerical results in Tables 1 and 2, respectively. It is worth mentioning that the numerical results have been obtained using the fourth-order Runge-Kutta in C++ program. We can see a very good agreement between the purely analytic results of the HAM and numerical results. The variation of with the change in the positive values of S for the two-dimensional case is plotted in Figure 4. Figure 5 shows the influence of negative S on for the axisymmetric case. Note that for the large negative values of , the results of similarity analysis are not reliable. gives skin friction, and represents pressure gradient. and as functions of S are illustrated in Figures 6 and 7, respectively.

 2nd-order 4th-order 6th-order 7th-order 8th-order Numerical 0.2 0.319474 0.319526 0.319526 0.319526 0.319526 0.319526 0.4 0.603652 0.603825 0.603830 0.603830 0.603830 0.603830 0.6 0.822574 0.822863 0.822875 0.822876 0.822876 0.822876 0.8 0.956580 0.956789 0.956800 0.956801 0.956801 0.956801 0.2 0.302545 0.302582 0.302582 0.302582 0.302582 0.302582 0.4 0.578028 0.578082 0.578082 0.578082 0.578082 0.578082 0.6 0.800737 0.800780 0.800780 0.800780 0.800780 0.800780 0.8 0.947686 0.947702 0.947702 0.947702 0.947702 0.947702 0.2 0.290353 0.290322 0.290322 0.290322 0.290322 0.290322 0.4 0.559299 0.559253 0.559252 0.559252 0.559252 0.559252 0.6 0.784341 0.784304 0.784303 0.784303 0.784303 0.784303 0.8 0.940717 0.940704 0.940703 0.940703 0.940703 0.940703 0.2 0.281032 0.281010 0.281010 0.281010 0.281010 0.281010 0.4 0.544851 0.544780 0.544779 0.544779 0.544779 0.544779 0.6 0.771493 0.771374 0.771371 0.771371 0.771371 0.771371 0.8 0.935127 0.935038 0.935036 0.935036 0.935036 0.935036 0.2 0.273767 0.273683 0.273682 0.273682 0.273682 0.273682 0.4 0.533516 0.533255 0.533247 0.533246 0.533246 0.533246 0.6 0.761299 0.760868 0.760848 0.760847 0.760847 0.760847 0.8 0.930617 0.930299 0.930281 0.930280 0.930280 0.930280
 2nd-order 4th-order 6th-order 7th-order 8th-order Numerical 0.2 0.332883 0.333591 0.333617 0.333618 0.333618 0.333618 0.4 0.623190 0.624315 0.624358 0.624358 0.624358 0.624358 0.6 0.838219 0.839284 0.839324 0.839325 0.839325 0.839325 0.8 0.962441 0.962961 0.962983 0.962984 0.962984 0.962984 0.2 0.305436 0.305543 0.305545 0.305545 0.305545 0.305545 0.4 0.582314 0.582468 0.582470 0.582470 0.582470 0.582470 0.6 0.804271 0.804390 0.804392 0.804392 0.804392 0.804392 0.8 0.949065 0.949107 0.949108 0.949108 0.949108 0.949108 0.2 0.288347 0.288261 0.288260 0.288260 0.288260 0.288260 0.4 0.556268 0.556145 0.556143 0.556143 0.556143 0.556143 0.6 0.781768 0.781670 0.781671 0.781671 0.781671 0.781671 0.8 0.939674 0.939641 0.939640 0.939640 0.939640 0.939640 0.2 0.276526 0.276433 0.276432 0.276432 0.276432 0.276432 0.4 0.537929 0.537754 0.537752 0.537752 0.537752 0.537752 0.6 0.765463 0.765252 0.765249 0.765249 0.765249 0.765249 0.8 0.932607 0.932474 0.932471 0.932471 0.932471 0.932471 0.2 0.268041 0.267797 0.267791 0.267791 0.267791 0.267791 0.4 0.524558 0.524057 0.524045 0.524045 0.524045 0.524045 0.6 0.753285 0.752627 0.752605 0.752605 0.752605 0.752605 0.8 0.927166 0.926724 0.926704 0.926703 0.926703 0.926703

#### 6. Conclusions

In this paper, the unsteady axisymmetric and two-dimensional squeezing flows between two parallel plates are studied analytically using the HAM. The convergence of the results is explicitly shown. Graphical results and tables are presented to investigate the influence of the squeeze number S on the velocity, skin friction, and pressure gradient. The solution obtained by means of the HAM is an infinite power series for appropriate initial approximation, which can be, in turn, expressed in a closed form. Unlike perturbation methods, the HAM does not depend on any small physical parameters. Thus, it is valid for both weakly and strongly nonlinear problems. Besides, different from all other analytic methods, the HAM provides us with a simple way to adjust and control the convergence region of the series solution by means of auxiliary parameter Thus the auxiliary parameter plays an important role within the frame of the HAM which can be determined by the so-called Consequently, the present success of the homotopy analysis method for the highly nonlinear problem of squeezing flows verifies that the method is a useful tool for nonlinear problems in science and engineering.

#### References

1. M. J. Stefan, “Versuch Über die scheinbare adhesion,” Sitzungsberichte der Akademie der Wissenschaften in Wien. Mathematik-Naturwissen, vol. 69, pp. 713–721, 1874. View at: Google Scholar
2. O. Reynolds, “On the theory of lubrication and its application to Mr. Beauchamp Tower's experiments, including an experimental determination of the viscosity of olive oil,” Philosophical Transactions of the Royal Society of London, vol. 177, pp. 157–234, 1886. View at: Publisher Site | Google Scholar
3. F. R. Archibald, “Load capacity and time relations for squeeze films,” Journal of Lubrication Technology, vol. 78, pp. A231–A245, 1956. View at: Google Scholar
4. R. J. Grimm, “Squeezing flows of Newtonian liquid films an analysis include the fluid inertia,” Applied Scientific Research, vol. 32, no. 2, pp. 149–166, 1976.
5. W. A. Wolfe, “Squeeze film pressures,” Applied Scientific Research, vol. 14, no. 1, pp. 77–90, 1965. View at: Publisher Site | Google Scholar
6. D. C. Kuzma, “Fluid inertia effects in squeeze films,” Applied Scientific Research, vol. 18, no. 1, pp. 15–20, 1968. View at: Publisher Site | Google Scholar
7. J. A. Tichy and W. O. Winer, “Inertial considerations in parallel circular squeeze film bearings,” Journal of Lubrication Technology, vol. 92, pp. 588–592, 1970. View at: Google Scholar
8. C. Y. Wang and L. T. Watson, “Squeezing of a viscous fluid between elliptic plates,” Applied Scientific Research, vol. 35, no. 2-3, pp. 195–207, 1979.
9. R. Usha and R. Sridharan, “Arbitrary squeezing of a viscous fluid between elliptic plates,” Fluid Dynamics Research, vol. 18, no. 1, pp. 35–51, 1996.
10. H. M. Laun, M. Rady, and O. Hassager, “Analytical solutions for squeeze flow with partial wall slip,” Journal of Non-Newtonian Fluid Mechanics, vol. 81, no. 1-2, pp. 1–15, 1999.
11. M. H. Hamdan and R. M. Baron, “Analysis of the squeezing flow of dusty fluids,” Applied Scientific Research, vol. 49, no. 4, pp. 345–354, 1992.
12. P. T. Nhan, “Squeeze flow of a viscoelastic solid,” Journal of Non-Newtonian Fluid Mechanics, vol. 95, no. 2-3, pp. 343–362, 2000. View at: Publisher Site | Google Scholar
13. J. D. Jackson, “A study of squeezing flow,” Applied Science Research A, vol. 11, pp. 148–152, 1962. View at: Google Scholar | Zentralblatt MATH
14. S. Ishizawa, “The unsteady flow between two parallel discs with arbitary varying gap width,” Bulletin of the Japan Society of Mechanical Engineers, vol. 9, no. 35, pp. 533–550, 1966. View at: Google Scholar
15. S.-J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems, Ph.D. thesis, Shanghai Jiao Tong University, Shanghai, China, 1992.
16. S.-J. Liao, “A uniformly valid analytic solution of two-dimensional viscous flow over a semi-infinite flat plate,” Journal of Fluid Mechanics, vol. 385, pp. 101–128, 1999.
17. S.-J. Liao, “An explicit, totally analytic approximate solution for Blasius' viscous flow problems,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 759–778, 1999. View at: Publisher Site | Google Scholar | MathSciNet
18. S.-J. Liao, “On the analytic solution of magnetohydrodynamic flows of non-Newtonian fluids over a stretching sheet,” Journal of Fluid Mechanics, vol. 488, pp. 189–212, 2003.
19. S.-J. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, vol. 2 of CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2004. View at: Zentralblatt MATH | MathSciNet
20. S.-J. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.
21. S.-J. Liao, “Comparison between the homotopy analysis method and homotopy perturbation method,” Applied Mathematics and Computation, vol. 169, no. 2, pp. 1186–1194, 2005.
22. M. Sajid, T. Hayat, and S. Asghar, “Comparison between the HAM and HPM solutions of thin film flows of non-Newtonian fluids on a moving belt,” Nonlinear Dynamics, vol. 50, no. 1-2, pp. 27–35, 2007. View at: Publisher Site | Google Scholar | MathSciNet
23. F. M. Allan, “Derivation of the Adomian decomposition method using the homotopy analysis method,” Applied Mathematics and Computation, vol. 190, no. 1, pp. 6–14, 2007.
24. S. Abbasbandy, “The application of homotopy analysis method to nonlinear equations arising in heat transfer,” Physics Letters A, vol. 360, no. 1, pp. 109–113, 2006.
25. S. Abbasbandy, “Approximate solution for the nonlinear model of diffusion and reaction in porous catalysts by means of the homotopy analysis method,” Chemical Engineering Journal, vol. 136, no. 2-3, pp. 144–150, 2008. View at: Publisher Site | Google Scholar
26. F. M. Allan, “Construction of analytic solution to chaotic dynamical systems using the Homotopy analysis method,” Chaos, Solitons & Fractals. In press. View at: Publisher Site | Google Scholar
27. F. M. Allan and M. I. Syam, “On the analytic solutions of the nonhomogeneous Blasius problem,” Journal of Computational and Applied Mathematics, vol. 182, no. 2, pp. 362–371, 2005.
28. T. Hayat and T. Javed, “On analytic solution for generalized three-dimensional MHD flow over a porous stretching sheet,” Physics Letters A, vol. 370, no. 3-4, pp. 243–250, 2007. View at: Publisher Site | Google Scholar | MathSciNet
29. M. Sajid, A. M. Siddiqui, and T. Hayat, “Wire coating analysis using MHD Oldroyd 8-constant fluid,” International Journal of Engineering Science, vol. 45, no. 2-8, pp. 381–392, 2007. View at: Publisher Site | Google Scholar
30. T. Hayat and M. Sajid, “Analytic solution for axisymmetric flow and heat transfer of a second grade fluid past a stretching sheet,” International Journal of Heat and Mass Transfer, vol. 50, no. 1-2, pp. 75–84, 2007.
31. T. Hayat, N. Ahmed, M. Sajid, and S. Asghar, “On the MHD flow of a second grade fluid in a porous channel,” Computers & Mathematics with Applications, vol. 54, no. 3, pp. 407–414, 2007.
32. T. Hayat, M. Sajid, and M. Ayub, “A note on series solution for generalized Couette flow,” Communications in Nonlinear Science and Numerical Simulation, vol. 12, no. 8, pp. 1481–1487, 2007.
33. Y. Bouremel, “Explicit series solution for the Glauert-jet problem by means of the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 12, no. 5, pp. 714–724, 2007.
34. M. M. Rashidi, G. Domairry, and S. Dinarvand, “Approximate solutions for the Burger and regularized long wave equations by means of the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 3, pp. 708–717, 2009.
35. Z. Ziabakhsh and G. Domairry, “Solution of the laminar viscous flow in a semi-porous channel in the presence of a uniform magnetic field by using the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1284–1294, 2009. View at: Publisher Site | Google Scholar
36. C.-Y. Wang, “The squeezing of fluid between two plates,” Journal of Applied Mechanics, vol. 43, no. 4, pp. 579–583, 1976. View at: Google Scholar | Zentralblatt MATH
37. T. Hayat, F. Shahzad, and M. Ayub, “Analytical solution for the steady flow of the third grade fluid in a porous half space,” Applied Mathematical Modelling, vol. 31, no. 11, pp. 2424–2432, 2007.