Mathematical Problems in Engineering

Volume 2010 (2010), Article ID 706840, 18 pages

http://dx.doi.org/10.1155/2010/706840

## Application of Homotopy Analysis Method to the Unsteady Squeezing Flow of a Second-Grade Fluid between Circular Plates

^{1}Department of Mechanical Engineering, Faculty of Engineering, Bu Ali Sina University, P.O. Box 65175-4161, Hamedan, Iran^{2}Department of Mathematics, Pennsylvania State University, York Campus, York, PA 17403, USA

Received 17 September 2010; Accepted 23 November 2010

Academic Editor: Victoria Vampa

Copyright © 2010 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.

#### Abstract

We investigated an axisymmetric unsteady two-dimensional flow of nonconducting, incompressible second grade fluid between two circular plates. The similarity transformation is applied to reduce governing partial differential equation (PDE) to a nonlinear ordinary differential equation (ODE) in dimensionless form. The resulting nonlinear boundary value problem is solved using homotopy analysis method and numerical method. The effects of appropriate dimensionless parameters on the velocity profiles are studied. The total resistance to the upper plate has been calculated.

#### 1. Introduction

Squeezing flows are induced by externally normal stresses or vertical velocities by means of moving boundary. Squeezing flows have many applications in food industry, especially in chemical engineering. 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. The study of squeezing flows has its origins in the 19th century and continues to receive considerable attention due to the practical applications in physical and biophysical areas. Stefan [1] published a classical paper on squeezing flow by using lubrication approximation. Such types of flow exist in lubrication when there is squeezing flow between two parallel plates. The tackiness of liquid adhesives also reflects squeeze film effects [2]. The squeeze film geometry has been studied extensively since 1947. Other applications in the biomechanics area relate to squeezing flow between parallel plates and the alternation between contraction and expansion of the blood vessels. In addition, polymer extrusion processes are modeled using squeezing flow of viscous fluids [3]. The squeezing flow between parallel plates when the confining walls have a transverse motion has great importance in hydrodynamic lubrication theory Langlois [4] and Salbu [5] have analyzed isothermal compressible squeeze films neglecting inertial effects. Thorpe [6] presented an explicit solution of the squeeze flow problem taking inertial terms into account. Later, P. S. Gupta and A. S. Gupta [7] showed that the solution given by [6] fails to satisfy the boundary conditions. Gupta solution, however, is restricted only to small Reynolds number. Squeeze film between two plane annuli with fluid inertia effects has been studied by Elkouh [8]. Some numerical solutions of squeezing flow between parallel plates has been conducted by Verma [9] and later by Singh et al. [10]. In addition, Hamza [11] has considered suction and injection effects on flow between parallel plates reflecting squeezing flow. Further work on second grade fluid flow between parallel plates has been carried out by Rajagopal and Gupta [12] and by Dandapat and Gupta [13]. Both papers consider the problem of flow between rotating parallel plates. Barron and Wiley [14] have also considered this fluid model to extend the Newtonian flow theory for slender bodies in a dusty gas. This fluid model has also been considered even in flow of a dusty liquid. For example, Debnath and Ghosh [15, 16] used this fluid model to study the flow of dusty fluids in rotating and oscillating cases, in the presence of magnetic fields. Hamdan and Barron [17] used this fluid model to analyse the steady squeezing flow of incompressible dusty fluids between two cylindrical plates. The theoretical and experimental studies of squeezing flows have been conducted by many researchers [18–28]. In this work the unsteady 2-dimensional flow of nonconducting, incompressible second grade fluid between circular plates that are moving symmetrically about the line of axial symmetry and giving rise to the squeezing flow is studied. The unsteady equations of motion are reduced to a single nonlinear fifth-order ordinary differential after employing a similarity transformation. The resulting nonlinear boundary value problem is solved with homotopy analysis method (HAM) and numerical method.

Most scientific problems and phenomena are modeled by nonlinear ordinary or partial differential equations. In most cases, these problems do not have precise analytic solution. Recently, much attention has been devoted to the newly developed methods to construct approximate analytic solutions of nonlinear equations. The HAM is developed in 1992 by Liao [28–35] and Rashidi et al. [36–38]. Liao applied this method to solve many types of nonlinear equations in science and engineering and then, this method has been successfully applied to solve many other nonlinear evolution equations. In recent years, the application of HAM in nonlinear problems has been developed by scientists and engineers. We know all perturbation methods require small parameter in nonlinear equation and the approximate solutions of equation containing this parameter are expressed as series expansions in the small parameter. Selection of small parameter requires a special skill. A proper choice of small parameter gives acceptable results, while an improper choice may result in incorrect solutions. The HAM, does not require a small parameter in equation modeling phenomena. Briefly speaking, Liao’s method is a universal one which can solve various kinds of nonlinear equations, thus many researchers applied this method to various linear and nonlinear problems. A substantial amount of research work has been invested in the study of linear and nonlinear systems of PDEs.

In Sections 2 and 3 description and mathematical formulation of the problem is presented. Sections 4 and 5 discussed the methods of solution of the problem. The results are discussed in Section 6 and Section 7 concludes finding of the present work.

#### 2. Problem Description

The unsteady two-dimensional squeezing flow of nonconducting, incompressible second grade fluid between two circular plates is studied. The distance between the plates at any time is . The central axis of the channel is taken as the -axis and -axis is normal to it. It is assumed that the circular plates move symmetrically with respect to the central axis and flow is axisymmetric about . The velocity components along the radial and axial directions are and , respectively (Figure 1). The fluid is assumed to have density and kinematics viscosity .

#### 3. Mathematical Formulation

The basic equations governing the motion of a homogeneous incompressible second grade fluid neglecting the thermal effects are [37]
where
where is the velocity vector, is the pressure, is the stress tensor, *μ* is the coefficient of viscosity and and are material constants. The Rivlin-Ericksen tensors are defined as
Now we formulate the equations of motion for unsteady two dimensional flow. Assuming that
and introducing the vorticity function and generalized pressure gradient as
where
We find that the equations of motion (3.1) and (3.2) reduce to
where is Laplacian operator.

Eliminating the generalized pressure between (3.10) and (3.11) to obtain The boundary conditions on and are where denotes the velocity of the plates. The first two conditions are due to the no-slip condition at the upper plate and the remaining two follow from the symmetry of the flow at .

If the dimensionless variable is introduced, (3.6), (3.9), and (3.12) become The boundary conditions on and are If we define the velocity components as follows Singh et al. [10] and Birkhoff [39] It follows that (3.14) takes the form substituting (3.18) and (3.19) in (3.15) and (3.16), we find that the continuity equation is identically satisfied and (3.16) becomes where prime denotes differentiation with respect to . The boundary conditions are determined from (3.17), (3.18), and (3.19) to be Thus for a similarity solution we define where , and and are functions of but for a similarity solution , , , and are considered to be constants. After integrating the first equation of (3.23), we have where and are constants. When and , the plates move apart symmetrically with respect to . In addition, when and , the plates approach each other and squeezing flow exists with similar velocity profile as long as . From (3.23) and (3.24), it follows that , and (3.21) becomes where This nonlinear boundary value problem is solved by the homotopy analysis method as well as numerically in the following sections.

#### 4. Basic Idea of the HAM

Let us consider the following differential equation where is a nonlinear operator, denotes independent variable, is an unknown function, respectively. For simplicity, we ignore all boundary or initial conditions, which can be treated in the similar way. By means of generalizing the traditional HAM, Liao [30] constructs the so-called zero-order deformation equation where is the embedding parameter, is a nonzero auxiliary parameter, is an auxiliary function, is an auxiliary linear operator, is an initial guess of is a unknown function, respectively. It is important, that one has great freedom to choose auxiliary things in the HAM. Obviously, when and , it holds respectively. Thus as increases from 0 to 1, the solution varies from the initial guess to the solution . Expanding in Taylor series with respect to , we have where If the auxiliary linear operator, the initial guess, the auxiliary parameter , and the auxiliary function are so properly chosen, the series (4.4) converges at , then we have which must be one of solutions of original nonlinear equation, as proved by Liao [30]. As and , (4.2) becomes which is used mostly in the homotopy analysis method, where as the solution obtained directly, without using Taylor series [38].

According to the definition (4.5), the governing equation can be deduced from the zero-order deformation equation (4.2). Define the vector Differentiating (4.2) times with respect to the embedding parameter and then setting and finally dividing them by , we have the so-called th-order deformation equation where It should be emphasized that for is governed by the linear equation (4.9) with the linear boundary conditions that come from original problem, which can be easily solved by symbolic computation software such as Maple and Mathematica.

#### 5. Application

Consider (3.22) and (3.25) and let us solve them through the HAM. Following the HAM, we start with initial approximation and the linear operator and nonlinear operator is defined as Using the above definition, with assumption , we construct the zero-order deformation equation Obviously, when and , Differentiating the zero-order deformation equation (5.3) times with respect to , and finally dividing by , we have the th-order deformation equation where Then We choose auxiliary functions as follow: Following the HAM and trying higher iterations with the unique and proper assignment of the results converge to the exact solution: using the symbolic software Mathematica to solve the system of linear equations, (5.5), with the boundary conditions (3.22), and successively obtain

#### 6. Results and Discussion

In this work the unsteady 2-dimensional flow of nonconducting, incompressible second grade fluid between two circular plates is considered. The squeezing flow is generated by moving the plates symmetrically with respect to the central region . The unsteady equations of motion are reduced to a single nonlinear 5th-order ordinary differential by using a similarity transformation. The resulting nonlinear boundary value problem is solved using the homotopy analysis method and numerical technique. The effect of the Reynolds number on velocity profiles is depicted. In these profiles we fixed the nonNewtonian parameters , and varied as (Figures 11 and 12). It is noted that the normal velocity increases as the Reynolds number increases. It is also observed that at a given time and for a fixed positive value of Reynolds number the normal velocity increases monotonically from to presents effect of Reynolds number on the longitudinal velocity. It is observed that this component of velocity increases near the central axis of the channel but deceases near the walls. The influence of nonNewtonian parameter on the normal and longitudinal velocity components is depicted. we varied fixing , . This observation on the velocity components is similar to that of Reynolds number presents the effects of nonNewtonian parameter on the squeezing flow of the second grade fluid. It is noted that effect of on the velocity components is opposite to that of the nonNewtonian parameter .

The HAM was applied successfully to find an explicit, totally analytic, uniformly valid solution for unsteady squeezing flow between circular parallel plates. The validity of our analytic solutions is verified by numerical results. The results show that HAM is very accurate (Figures 3, 4, 6, 7, 9, and 10). Also, the HAM gives rapidly convergent series with specific significant features for each scheme (Figures 2, 5, and 8). The accuracy of the method is very good.

#### 7. Conclusions

A similarity solution of the problem for the unsteady flow of nonconducting, incompressible second grade fluid between two circular plates approaching or receding from each other has been investigated. It is noted that a similarity solution exists only when the distance between the plates varies as , and squeezing flow takes place for and as for as , Approximate solutions for the fluid velocity have been found for the flow of unsteady 2-dimensional second grade fluid between two circular plates by the homotopy analysis method and numerical method. The major finding of the present paper can be summarized as follows:The Reynolds number and nonNewtonian parameter have similar effect on the normal and longitudinal velocity components.The nonNewtonian parameter has opposite trend on the velocity profiles to that of Reynolds number and . It is found that at a given time and for a fixed positive value of Reynolds number, nonNewtonian parameters , , the normal velocity increases monotonically from to . The longitudinal component of velocity increases near the central axis of the channel but deceases near the walls when increases from 1.0 to 2.0. Similar effect of on this component is observed when varies from 0.0 to 0.1. The results are recovered for .The obtained results are valid for all values of Reynolds number and nonNewtonian parameters.

#### Nomenclature

: | Density |

: | Kinematic viscosity |

: | Velocity vector |

: | Pressure |

: | Stress tensor |

μ: | Viscosity coefficient |

: | Material constant |

: | Material constant |

: | Vorticity function |

: | Pressure gradient |

: | Laplacian operator |

: | Dimensionless variable |

: | Auxiliary linear parameter |

: | Rivlin-Ericksen tensor |

: | Velocity component along the radial |

: | Velocity component along the axial |

: | Distance between the plates at any time |

: | Non-Newton parameter |

: | Non-Newton parameter |

: | Reynolds number. |

#### References

- M. J. Stefan, “Versuch Uber die scheinbare adhesion,”
*Akademie der Wissenschaften in Wien. Mathematik-Naturwissen*, vol. 69, p. 713, 1874. View at Google Scholar - J. J. Bickerman,
*Development in Theoretical and Applied Mathematics*, vol. 3, Academic Press, New York, NY, USA, 1958. - D. Yao, V. L. Virupaksha, and B. Kim, “Study on squeezing flow during nonisothermal embossing of polymer microstructures,”
*Polymer Engineering and Science*, vol. 45, no. 5, pp. 652–660, 2005. View at Publisher · View at Google Scholar · View at Scopus - W. E. Langlois, “Isothermal squeeze films,”
*Applied Mathematics*, vol. 20, p. 131, 1962. View at Google Scholar - E. O. Salbu, “Compressible squeeze films and squeeze bearings,”
*Journal of Basic Engineering*, vol. 86, p. 355, 1964. View at Google Scholar - J. F. Thorpe, in
*Development in Theoretical and Applied Mathematics*, W. A. Shah, Ed., vol. 3, Pergamon Press, Oxford, UK, 1967. - P. S. Gupta and A. S. Gupta, “Squeezing flow between parallel plates,”
*Wear*, vol. 45, no. 2, pp. 177–185, 1977. View at Publisher · View at Google Scholar · View at Scopus - A. F. Elkouh, “Fluid inertia effects in a squeeze film between two plane annuli,”
*Journal of Tribology*, vol. 106, no. 2, pp. 223–227, 1984. View at Publisher · View at Google Scholar · View at Scopus - R. L. Verma, “A numerical solution for squeezing flow between parallel channels,”
*Wear*, vol. 72, no. 1, pp. 89–95, 1981. View at Publisher · View at Google Scholar · View at Scopus - P. Singh, V. Radhakrishnan, and K. A. Narayan, “Squeezing flow between parallel plates,”
*Ingenieur-Archiv*, vol. 60, no. 4, pp. 274–281, 1990. View at Publisher · View at Google Scholar · View at Scopus - E. A. Hamza, “Suction and injection effects on a similar flow between parallel plates,”
*Journal of Physics D: Applied Physics*, vol. 32, no. 6, pp. 656–663, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - K. R. Rajagopal and A. S. Gupta, “On a class of exact solutions to the equations of motion of a second grade fluid,”
*International Journal of Engineering Science*, vol. 19, no. 7, pp. 1009–1014, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. S. Dandapat and A. S. Gupta, “Stability of a thin layer of a second-grade fluid on a rotating disk,”
*International Journal of Non-Linear Mechanics*, vol. 26, no. 3-4, pp. 409–417, 1991. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - R. M. Barron and J. T. Wiley, “Newtonian flow theory for slender bodies in a dusty gas,”
*Journal of Fluid Mechanics*, vol. 108, pp. 147–157, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - A. K. Ghosh and L. Debnath, “Hydromagnetic Stokes flow in a rotating fluid with suspended small particles,”
*Applied Scientific Research*, vol. 43, no. 3, pp. 165–192, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - Lokenath Debnath and A. K. Ghosh, “On unsteady hydromagnetic flows of a dusty fluid between two oscillating plates,”
*Applied Scientific Research*, vol. 45, no. 4, pp. 353–365, 1988. View at Google Scholar · View at Zentralblatt MATH - M. H. Hamdan and R. M. Barron, “Analysis of the squeezing flow of dusty fluids,”
*Applied Scientific Research*, vol. 49, no. 4, pp. 345–354, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - W. A. Wolfe, “Squeeze film pressures,”
*Applied Scientific Research*, vol. 14, no. 1, pp. 77–90, 1965. View at Publisher · View at Google Scholar · View at Scopus - D. C. Kuzma, “Fluid inertia effects in squeeze films,”
*Applied Scientific Research*, vol. 18, no. 1, pp. 15–20, 1968. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - . TICHY JA and . WINER WO, “Inertial considerations in parallel circular squeeze film bearings,”
*Transactions of the ASME: Journal of Lubrication Technology*, vol. 92, pp. 588–592, 1970. View at Google Scholar · View at Scopus - 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. View at Google Scholar - 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 · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. H. Hamdan and R. M. Barron, “Analysis of the squeezing flow of dusty fluids,”
*Applied Scientific Research*, vol. 49, no. 4, pp. 345–354, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - P. T. Nhan, “Squeeze flow of a viscoelastic solid,”
*Journal of Non-Newtonian Fluid Mechanics*, vol. 95, pp. 43–362, 2000. View at Google Scholar · View at Zentralblatt MATH - S. J. Liao,
*The proposed homotopy analysis technique for the solution of nonlinear problems*, Ph.D. thesis, Shanghai Jiao Tong University, 1992. - 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 · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. J. Liao,
*Beyond Perturbation*, vol. 2 of*Introduction to the Homotopy Analysis Method*, Chapman & Hall/CRC, Boca Raton, FL, 2004. View at Zentralblatt MATH - 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. Liao, “A new branch of solutions of boundary-layer flows over an impermeable stretched plate,”
*International Journal of Heat and Mass Transfer*, vol. 48, no. 12, pp. 2529–2539, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. Liao, J. Su, and A. T. Chwang, “Series solutions for a nonlinear model of combined convective and radiative cooling of a spherical body,”
*International Journal of Heat and Mass Transfer*, vol. 49, no. 15-16, pp. 2437–2445, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. M. Rashidi and S. A. Mohimanian Pour, “Analytic approximate solutions for unsteady boundary-layer flow and heat transfer due to a stretching sheet by homotopy analysis method,”
*Nonlinear Analysis. Modelling and Control*, vol. 15, no. 1, pp. 83–95, 2010. View at Google Scholar · View at Zentralblatt MATH - M. M. Rashidi, D. D. Ganji, and S. Dinarvand, “Approximate traveling wave solutions of coupled Whitham-Broer-Kaup shallow water equations by homotopy analysis method,”
*Differential Equations & Nonlinear Mechanics*, vol. 2008, Article ID 243459, 8 pages, 2008. View at Google Scholar - 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 · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - G. Birkhoff,
*Hydrodynamics*, Princeton University Press, Princeton, NJ, USA, 1960.