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

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 2𝑎(𝑡). 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 𝑧=0 and flow is axisymmetric about 𝑟=0. The velocity components along the radial and axial directions are 𝑢(𝑟,𝑧,𝑡) and 𝑤(𝑟,𝑧,𝑡), respectively (Figure 1). The fluid is assumed to have density 𝜌 and kinematics viscosity 𝜐.

Figure 1 

Geometry of the problem.

3. Mathematical Formulation

The basic equations governing the motion of a homogeneous incompressible second grade fluid neglecting the thermal effects are [37]1𝑟𝜕()+𝜕𝜕𝑟𝑟𝑢(𝑤)𝜌𝜕𝑧=0,𝜕𝑢𝜕𝑡+𝑢𝜕𝑢𝜕𝑟+𝑤𝜕𝑢𝜕𝑧=−𝜕𝑃𝜕𝜕𝑟+𝜇²𝑢𝜕𝑟²+1𝑟𝜕𝑢−𝑢𝜕𝑟𝑟²+𝜕²𝑢𝜕𝑧²+𝐹_𝑟,𝜌𝜕𝑤𝜕𝑡+𝑢𝜕𝑤𝜕𝑟+𝑤𝜕𝑤𝜕𝑧=−𝜕𝑃𝜕𝜕𝑧+𝜇²𝑤𝜕𝑟²+1𝑟𝜕𝑤+𝜕𝜕𝑟²𝑤𝜕𝑧²+𝐹_𝑧, (3.1)(3.2) where 𝑇=−𝑝𝐼+𝜇𝐴₁+𝛼₁𝐴₂+𝛼₂𝐴²₁, (3.3) 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𝐴₁=∇𝑉+(∇𝑉)^𝑇,𝐴₂=𝑑𝐴₁𝑑𝑡+𝐴₁(∇𝑉)+(∇𝑉)^𝑇𝐴₁. (3.4) Now we formulate the equations of motion for unsteady two dimensional flow. Assuming that 𝑉=[𝑢(𝑟,𝑧,𝑡),0,𝑤(𝑟,𝑧,𝑡)], (3.5) and introducing the vorticity function Ω(𝑟,𝑧,𝑡) and generalized pressure gradient ℎ(𝑟,𝑧,𝑡) as Ω=𝜕𝑤−𝜕𝑟𝜕𝑢,𝜌𝜕𝑧ℎ=2𝑢²+𝑤²+𝑝−𝛼₁𝑤𝜕+1𝜕𝑟𝑟𝜕−𝑢1𝜕𝑧Ω−43𝛼₁+2𝛼₂𝐴²₁, (3.6)(3.7) where 𝐴²₁=4𝜕𝑢𝜕𝑟²+4𝜕𝑤𝜕𝑧²𝑢+4𝑟²+2𝜕𝑢+𝜕𝑧𝜕𝑤𝜕𝑟². (3.8) We find that the equations of motion (3.1) and (3.2) reduce to𝜕𝑢+𝑢𝜕𝑟𝑟+𝜕𝑤𝜕𝑧=0,𝜕ℎ𝜕𝑟+𝜌𝜕𝑢𝜕𝑡−𝑤Ω=−𝜇+𝛼₁𝜕𝜕𝑡𝜕Ω𝜕𝑡−𝛼₁𝑤𝜕²𝜕𝑟²+1𝑟𝜕+𝜕𝜕𝑟𝜕𝑧²²ΩΩ−𝑟²+𝛼₁+𝛼₂2𝑟𝜕()+Ω𝜕𝑧𝑢Ω²𝑟,𝜕ℎ𝜕𝑧+𝜌𝜕𝑤=𝜕𝑡+𝑢Ω𝜇+𝛼₁𝜕𝜕𝑡𝜕Ω+Ω𝜕𝑟𝑟+𝛼₁𝑢𝜕²𝜕𝑟²+1𝑟𝜕+𝜕𝜕𝑟𝜕𝑧²²ΩΩ−𝑟²−𝛼₁+𝛼₂2𝑟𝜕(),𝜕𝑧𝑢Ω (3.9)(3.10)(3.11) where ∇²=𝜕²/𝜕𝑟²+(1/𝑟)(𝜕/𝜕𝑟)+𝜕²/𝜕𝑧² is Laplacian operator.

Eliminating the generalized pressure ℎ between (3.10) and (3.11) to obtain𝜌𝜕Ω+𝑢𝜕𝜕𝑡𝜕𝜕𝑟+𝑤−𝑢𝜕𝑧𝑟Ω=𝜇+𝛼₁𝜕𝜕𝜕𝑡𝜕𝑟²²+1𝑟𝜕+𝜕𝜕𝑟𝜕𝑧²²ΩΩ−𝑟²+𝛼₁𝑢𝜕𝜕𝜕𝑟+𝑤−𝑢𝜕𝑧𝑟𝜕²𝜕𝑟²+1𝑟𝜕+𝜕𝜕𝑟𝜕𝑧²²ΩΩ−𝑟²𝛼−2¹+𝛼₂𝑟𝜕²𝜕𝑟²−1𝑟𝜕+𝜕𝜕𝑟𝜕𝑧²²1𝑢Ω+2𝜕Ω𝜕𝑧². (3.12) The boundary conditions on 𝑢(𝑟,𝑧,𝑡) and 𝑤(𝑟,𝑧,𝑡) are𝑢(𝑟,𝑧,𝑡)=0,𝑤(𝑟,𝑧,𝑡)=𝑣_𝑤(𝑡)𝑤()(),𝑧=𝑎,𝑟,𝑧,𝑡=0,𝜕𝑢𝑟,𝑧,𝑡𝜕𝑧=0,𝑧=0, (3.13) 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 𝑧=0.

If the dimensionless variable 𝜂=𝑧/𝑎(𝑡) is introduced, (3.6), (3.9), and (3.12) become𝜕𝑤−𝜕𝑟𝜕𝑢𝑎𝜕𝜂=Ω,𝜕𝑢+𝑢𝜕𝑟𝑟+𝜕𝑤𝜌𝑎𝜕𝜂=0,𝜕Ω+𝑢𝜕𝜕𝑡𝜕𝜕𝑟+𝑤−𝑢𝑎𝜕𝜂𝑟Ω=𝜇+𝛼₁𝜕𝜕𝜕𝑡𝜕𝑟²²+1𝑟𝜕+𝜕𝜕𝑟𝜕𝑧²²ΩΩ−𝑟²+𝛼₁𝑢𝜕𝜕𝜕𝑟+𝑤−𝑢𝑎𝜕𝜂𝑟𝜕²𝜕𝑟²+1𝑟𝜕+𝜕𝜕𝑟𝜕𝑧²²ΩΩ−𝑟²𝛼−2¹+𝛼₂𝑟𝜕²𝜕𝑟²−1𝑟𝜕+𝜕𝜕𝑟𝜕𝜂²²1𝑢Ω+𝜕2𝑎Ω𝜕𝜂². (3.14)(3.15)(3.16) The boundary conditions on 𝑢(𝑟,𝑧,𝑡) and 𝑤(𝑟,𝑧,𝑡) are𝑢𝑟,𝜂,𝑡=0,𝑤𝑟,𝜂,𝑡=𝑣_𝑤(𝑡)𝑤,𝜂=1,𝑟,𝜂,𝑡=0,𝜕𝑢𝑟,𝜂,𝑡𝜕𝜂=0,𝜂=0. (3.17)(3.18) If we define the velocity components as follows Singh et al. [10] and Birkhoff [39]𝑢=−𝑟(𝑡)𝑣2𝑎^𝑤(𝑡)𝑓^′𝜂,𝑤=𝑣_𝑤(𝑡)𝑓𝜂. (3.19) It follows that (3.14) takes the form 𝑟Ω=2𝑎(𝑡)²𝑣_𝑤(𝑡)𝑓^′′𝜂, (3.20) 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𝑎𝑣_𝑤𝜐𝑓𝑓^′′′−2𝑓^′′−𝜂𝑓^′′′+𝑎𝑑𝑣_𝑤/𝑑𝑡𝜐𝑣_𝑤𝑓^′′=𝑓^𝑖𝑣+𝛼₁+𝛼₂𝑣_𝑤𝜇𝑎2𝑓^′′𝑓^′′′+𝑓^′𝑓^𝑖𝑣+𝛼₁𝑣_𝑤𝜇𝑎𝑓𝑓^𝑣−𝜂𝑓^𝑣−4𝑓^𝑖𝑣+𝑎𝑑𝑣_𝑤/𝑑𝑡𝜐𝑣²_𝑤𝑓^𝑖𝑣, (3.21) where prime denotes differentiation with respect to 𝜂. The boundary conditions are determined from (3.17), (3.18), and (3.19) to be𝑓(1)=1,𝑓^′(1)𝑓(0)=0,=0,𝑓^′′(0)=0. (3.22) Thus for a similarity solution we define𝑎𝑣_𝑤𝜐=𝑅_𝑒,𝑎²𝑑𝑣_𝑤/𝑑𝑡𝜐𝑣_𝑤=𝑅_𝑒𝛼𝑄,₁𝑣_𝑤𝜇𝑎=𝑊_𝑒1,𝛼₂𝑣_𝑤𝜇𝑎=𝑊_𝑒2, (3.23) where 𝑅_𝑒, 𝑄 and 𝑊_𝑒1 and 𝑊_𝑒2 are functions of 𝑡 but for a similarity solution 𝑅_𝑒, 𝑄, 𝑊_𝑒1, and 𝑊_𝑒2 are considered to be constants. After integrating the first equation of (3.23), we have𝑎(𝑡)=𝐾𝑡+𝑎₀^1/2, (3.24) where 𝐾 and 𝑎₀ are constants. When 𝐾>0 and 𝑎₀>0, the plates move apart symmetrically with respect to 𝜂=0. In addition, when 𝐾<0 and 𝑎₀>0, the plates approach each other and squeezing flow exists with similar velocity profile as long as (𝐾𝑡+𝑎₀)>0. From (3.23) and (3.24), it follows that 𝑄=−1, and (3.21) becomes𝑓^𝑖𝑣+𝐾²𝑓^′′=𝐿𝑓−𝜂𝑓^′′′+𝑀𝜂−𝑓𝑓^𝑣−𝑁2𝑓^′′𝑓^′′′+𝑓^′𝑓^𝑖𝑣, (3.25) where𝐾=3𝑅_𝑒1−5𝑊_𝑒1,𝑅𝐿=^𝑒1−5𝑊_𝑒1,𝑊𝑀=^𝑒11−5𝑊_𝑒1,𝑊𝑁=^𝑒1+𝑊_𝑒21−5𝑊_𝑒1. (3.26) 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𝒩𝑓𝜂=0, (4.1) 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 1−𝑝ℒ𝜑𝜂;𝑝−𝑓₀𝜂=𝑝ℏ𝐻𝜂𝒩𝜑𝜂;𝑝, (4.2) where 𝑝∈[0,1] is the embedding parameter, ℏ≠0 is a nonzero auxiliary parameter, 𝐻(𝜂)≠0 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 𝑝=0 and 𝑝=1, it holds 𝜑𝜂;0=𝑓₀𝜂,𝜑𝜂;1=𝑢𝜂, (4.3) 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𝜑𝜂;𝑝=𝑓₀𝜂+^+∞_𝑚=1𝑓_𝑚𝜂𝑝^𝑚, (4.4) where 𝑓_𝑚𝜂=1𝜕𝑚!^𝑚𝜑𝜂;𝑝𝜕𝑝^𝑚_𝑝=0. (4.5) 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 𝑝=1, then we have𝑓𝜂=𝑓₀𝜂+^+∞_𝑚=1𝑓𝜂, (4.6) which must be one of solutions of original nonlinear equation, as proved by Liao [30]. As ℏ=−1 and 𝐻(𝜂)=1, (4.2) becomes1−𝑝ℒ𝜑𝜂;𝑝−𝑓₀𝜂+𝑝𝒩𝜑𝜂;𝑝=0, (4.7) 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𝑓_𝑛𝜂=𝑓₀𝜂,𝑓₁𝜂,…,𝑓_𝑛𝜂. (4.8) Differentiating (4.2) 𝑚 times with respect to the embedding parameter 𝑝 and then setting 𝑝=0 and finally dividing them by 𝑚!, we have the so-called 𝑚th-order deformation equationℒ𝑓_𝑚𝜂−𝜒_𝑚𝑓_𝑚−1𝜂=ℏ𝐻𝜂𝑅_𝑚𝑓_𝑚−1, (4.9) where 𝑅_𝑚𝑓_𝑚−1=1()!𝜕𝑚−1^𝑚−1𝒩𝜑𝜂;𝑝𝜕𝑝^𝑚−1_𝑝=0,𝜒_𝑚=0,𝑚≤1,1,𝑚>1. (4.10) It should be emphasized that 𝑓_𝑚(𝜂) for 𝑚≥1 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 𝑓₀(𝜂)=−0.5𝜂³+1.5𝜂 and the linear operator ℒ𝑓𝜂=𝜕⁴𝜑𝜂;𝑝𝜕𝜂⁴, (5.1) and nonlinear operator is defined as𝒩𝑓𝜂𝜕=−⁴𝑓𝜕𝜂⁴−𝐾²𝜕²𝑓𝜕𝜂²𝜕+𝐿𝑓−𝜂³𝑓𝜕𝜂³𝜕+𝑀𝜂−𝑓⁵𝑓𝜕𝜂⁵2𝜕−𝑁²𝑓𝜕𝜂²𝜕³𝑓𝜕𝜂³+𝜕𝑓𝜕𝜕𝜂⁴𝑓𝜕𝜂⁴. (5.2) Using the above definition, with assumption 𝐻(𝜂)=1, we construct the zero-order deformation equation1−𝑝ℒ𝜑𝜂;𝑝−𝑓₀𝜂=𝑝ℏ𝒩𝜑𝜂;𝑝. (5.3) Obviously, when 𝑝=0 and 𝑝=1,𝜑𝜂;0=𝑓₀𝜂,𝜑𝜂;1=𝑓𝜂. (5.4) Differentiating the zero-order deformation equation (5.3) 𝑚 times with respect to 𝑝, and finally dividing by 𝑚!, we have the 𝑚th-order deformation equationℒ𝑓_𝑚𝜂−𝜒_𝑚𝑓_𝑚−1𝜂=ℏ𝑅_𝑚𝑓_𝑚−1, (5.5) where 𝑅𝑓_𝑚−1𝜕=−⁴𝑓_𝑚−1𝜂𝜕𝜂⁴−𝐾²𝜕²𝑓_𝑚−1𝜂𝜕𝜂²𝜕−𝐿𝜂³𝑓_𝑚−1𝜂𝜕𝜂³𝜕+𝑀𝜂⁵𝑓_𝑚−1𝜂𝜕𝜂⁵+𝐿^𝑚−1_𝑛=0𝑓_𝑛𝜂𝜕³𝑓_{𝑚−1−𝑛}𝜂𝜕𝜂³−𝑀^𝑚−1_𝑛=0𝑓_𝑛𝜂𝜕⁵𝑓_{𝑚−1−𝑛}𝜂𝜕𝜂⁵−2𝑁^𝑚−1_𝑛=0𝜕²𝑓_𝑛𝜂𝜕𝜂²𝜕³𝑓_{𝑚−1−𝑛}𝜂𝜕𝜂³−𝑁^𝑚−1_𝑛=0𝜕𝑓_𝑛𝜂𝜕𝜕𝜂⁴𝑓_{𝑚−1−𝑛}𝜂𝜕𝜂⁴. (5.6) Then ℒ𝑓_𝑚𝜂−𝜒_𝑚𝑓_𝑚−1𝜂=ℏ𝑅_𝑚𝑓_𝑚−1,𝜒_𝑚=0,𝑚≤1,1,𝑚>1. (5.7) We choose auxiliary functions as follow:ℋ𝜂=1. (5.8) Following the HAM and trying higher iterations with the unique and proper assignment of the results converge to the exact solution:𝑓𝜂≈𝑓₀𝜂+𝑓₁𝜂+𝑓₂𝜂+⋯+𝑓_𝑚𝜂, (5.9) using the symbolic software Mathematica to solve the system of linear equations, (5.5), with the boundary conditions (3.22), and successively obtain𝑓₁𝜂=−0.058567ℏ𝜂+0.115348ℏ𝜂³−0.054996ℏ𝜂⁵−0.001786ℏ𝜂⁷ (5.10)

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 𝑧=0. 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 𝑊_𝑒1=0.0, 𝑊_𝑒2=0.05 and varied 𝑅_𝑒 as  𝑅_𝑒=1.0,1.5,2.0 (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 𝜂=0 to 𝜂=1 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 𝑊_𝑒1 on the normal and longitudinal velocity components is depicted. we varied 𝑊_𝑒1=0.0,0.05,0.75 fixing 𝑊_𝑒2=0.05, 𝑅_𝑒=1.0. This observation on the velocity components is similar to that of Reynolds number 𝑅_𝑒 presents the effects of nonNewtonian parameter 𝑊_𝑒2 on the squeezing flow of the second grade fluid. It is noted that effect of 𝑊_𝑒2 on the velocity components is opposite to that of the nonNewtonian parameter 𝑊_𝑒1.

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.

Figure 2 

ℏ curve by the HAM with 𝑥=0.00001, 𝐾=1.732, 𝐿=1, 𝑀=0, and 𝑁=0.05, by 20th-order approximation solution.

Figure 3 

Obtained solution by the HAM with ℏ=−1, 𝐾=1.732, 𝐿=1, 𝑀=0, and 𝑁=0.05, by 20th-order approximation solution.

Figure 4 

Obtained solution by the HAM with ℏ=−1, 𝐾=1.732, 𝐿=1, 𝑀=0, and 𝑁=0.05, by 20th-order approximation solution.

Figure 5 

ℏ curve by the HAM with 𝑥=0.00001, 𝐾=2.121, 𝐿=1.5, 𝑀=0, and 𝑁=0.05, by 20th-order approximation solution.

Figure 6 

Obtained solution by the HAM with ℏ=−1, 𝐾=2.121, 𝐿=1.5, 𝑀=0, and 𝑁=0.05, by 20th-order approximation solution.

Figure 7 

Obtained solution by the HAM with ℏ=−1, 𝐾=2.121, 𝐿=1.5, 𝑀=0, and 𝑁=0.05, by 20th-order approximation solution.

Figure 8 

ℏ curve by the HAM with 𝑥=0.00001, 𝐾=2.449, 𝐿=2, 𝑀=0, and 𝑁=0.05, by 10th-order approximation solution.

Figure 9 

Obtained solution by the HAM with ℏ=−1, 𝐾=2.449, 𝐿=2, 𝑀=0, and 𝑁=0.05, by 20th-order approximation solution.

Figure 10 

Obtained solution by the HAM with ℏ=−1, 𝐾=2.449, 𝐿=2, 𝑀=0, and 𝑁=0.05, by 20th-order approximation solution.

Figure 11 

The effect of 𝑅𝑒 on the axial velocity profiles.

Figure 12 

The effect of 𝑅𝑒 on the radial velocity profiles.

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 (𝐾𝑡+𝑎₀)^1/2, and squeezing flow takes place for 𝐾<0 and 𝑎₀>0 as for as 𝐾𝑡+𝑎₀>0, 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 𝑊_𝑒1 have similar effect on the normal and longitudinal velocity components.The nonNewtonian parameter 𝑊_𝑒2 has opposite trend on the velocity profiles to that of Reynolds number and 𝑊_𝑒1. It is found that at a given time and for a fixed positive value of Reynolds number, nonNewtonian parameters 𝑊_𝑒1, 𝑊_𝑒2, the normal velocity increases monotonically from 𝜂=0 to 𝜂=1. 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 𝑊_𝑒2 on this component is observed when 𝑊_𝑒2 varies from 0.0 to 0.1. The results are recovered for 𝑊_𝑒→0.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
𝛼1:	Material constant
𝛼2:	Material constant
Ω:	Vorticity function
ℎ:	Pressure gradient
∇2:	Laplacian operator
𝜂:	Dimensionless variable
ℏ:	Auxiliary linear parameter
𝐴𝑖:	Rivlin-Ericksen tensor
𝑢:	Velocity component along the radial
𝑤:	Velocity component along the axial
2𝑎(𝑡):	Distance between the plates at any time 𝑡
𝑊𝑒1:	Non-Newton parameter
𝑊𝑒2:	Non-Newton parameter
𝑅𝑒:	Reynolds number.