Research Article  Open Access
Analytic Approximate Solutions for Unsteady TwoDimensional and Axisymmetric Squeezing Flows between Parallel Plates
Abstract
The flow of a viscous incompressible fluid between two parallel plates due to the normal motion of the plates is investigated. The unsteady NavierStokes equations are reduced to a nonlinear fourthorder 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 fourthorder RungeKutta.
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, 4–14]. 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 NavierStokes 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) [16–21]. 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 threedimensional MHD flow over a porous stretching sheet [28], the wire coating analysis using MHD Oldroyd 8constant fluid [29], the axisymmetric flow and heat transfer of a secondgrade fluid past a stretching sheet [30], the MHD flow of a secondgrade fluid in a porous channel [31], the generalized Couette flow [32], the Glauertjet 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 twodimensional 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 twodimensional flow, Wang introduced the following transforms [36]: where Substituting (2.1) into the unsteady twodimensional NavierStokes 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 socalled zerothorder 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 zerothorder 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 mthorder deformation equation, we differentiate (3.5) m times with respect to p, divide by m!, and then set The resulting deformation equation at the mthorder 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 socalled 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 15thorder approximation for the axisymmetric and twodimensional 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 twodimensional 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 fourthorder RungeKutta 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 twodimensional 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.


6. Conclusions
In this paper, the unsteady axisymmetric and twodimensional 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 socalled 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
 M. J. Stefan, “Versuch Über die scheinbare adhesion,” Sitzungsberichte der Akademie der Wissenschaften in Wien. MathematikNaturwissen, vol. 69, pp. 713–721, 1874. View at: Google Scholar
 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
 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
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 W. A. Wolfe, “Squeeze film pressures,” Applied Scientific Research, vol. 14, no. 1, pp. 77–90, 1965. View at: Publisher Site  Google Scholar
 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
 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
 C. Y. Wang and L. T. Watson, “Squeezing of a viscous fluid between elliptic plates,” Applied Scientific Research, vol. 35, no. 23, pp. 195–207, 1979. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. M. Laun, M. Rady, and O. Hassager, “Analytical solutions for squeeze flow with partial wall slip,” Journal of NonNewtonian Fluid Mechanics, vol. 81, no. 12, pp. 1–15, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. T. Nhan, “Squeeze flow of a viscoelastic solid,” Journal of NonNewtonian Fluid Mechanics, vol. 95, no. 23, pp. 343–362, 2000. View at: Publisher Site  Google Scholar
 J. D. Jackson, “A study of squeezing flow,” Applied Science Research A, vol. 11, pp. 148–152, 1962. View at: Google Scholar  Zentralblatt MATH
 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
 S.J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems, Ph.D. thesis, Shanghai Jiao Tong University, Shanghai, China, 1992.
 S.J. Liao, “A uniformly valid analytic solution of twodimensional viscous flow over a semiinfinite flat plate,” Journal of Fluid Mechanics, vol. 385, pp. 101–128, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S.J. Liao, “An explicit, totally analytic approximate solution for Blasius' viscous flow problems,” International Journal of NonLinear Mechanics, vol. 34, no. 4, pp. 759–778, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 S.J. Liao, “On the analytic solution of magnetohydrodynamic flows of nonNewtonian fluids over a stretching sheet,” Journal of Fluid Mechanics, vol. 488, pp. 189–212, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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
 S.J. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Sajid, T. Hayat, and S. Asghar, “Comparison between the HAM and HPM solutions of thin film flows of nonNewtonian fluids on a moving belt,” Nonlinear Dynamics, vol. 50, no. 12, pp. 27–35, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. 23, pp. 144–150, 2008. View at: Publisher Site  Google Scholar
 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
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Hayat and T. Javed, “On analytic solution for generalized threedimensional MHD flow over a porous stretching sheet,” Physics Letters A, vol. 370, no. 34, pp. 243–250, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 M. Sajid, A. M. Siddiqui, and T. Hayat, “Wire coating analysis using MHD Oldroyd 8constant fluid,” International Journal of Engineering Science, vol. 45, no. 28, pp. 381–392, 2007. View at: Publisher Site  Google Scholar
 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. 12, pp. 75–84, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Bouremel, “Explicit series solution for the Glauertjet problem by means of the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 12, no. 5, pp. 714–724, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. Ziabakhsh and G. Domairry, “Solution of the laminar viscous flow in a semiporous 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
 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
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2008 Mohammad Mehdi Rashidi 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.