Couple Stress Squeeze Films with VPD in a Curved Circular Geometry

The problem of couple stress squeeze films considering viscosity pressure dependence (VPD) has been analysed in a curved circular geometry. Using Stokes microcontinuum theory and the Barus formula, the Reynolds type equation has been derived. The approximate analytical expressions for the squeeze film pressure and load carrying capacity are obtained using a perturbation technique.The numerical solutions for the squeeze film pressure and load carrying capacity are presented for the sinusoidal motion of the upper curved disk, assuming an exponential form for the curvature. The effects of curvature, the non-Newtonian couple stresses, and VPD and their combined effects are investigated through the squeeze film pressure and the load carrying capacity of the squeeze film.


Introduction
Many areas of Engineering and Science such as lubricated joints, gears, viscous dampers, bearings, and automotive and aircraft engines witness the major role of squeeze film phenomena.Squeeze film is a thin layer of viscous fluid between approaching surfaces.As the viscous lubricant is contained between the surfaces, it takes certain time for these surfaces to come into contact with each other and the lubricant cannot be squeezed out instantaneously.Since the viscous lubricant has a resistance to extrusion, a pressure is built up during this interval and then the load is supported by the lubricant film.
In most of the studies of classical hydrodynamic lubrication, it is assumed that the lubricant behaves as a Newtonian fluid as studied by Pinkus and Sternlicht [1] and Jones and Wilson [2].In recent years, experimental results show that the addition of a small amount of long-chained polymers to a Newtonian fluid gives the most desirable lubricant.It has been experimentally showed by Spikes [3] that base oil blended with viscosity index improvers can enhance the behaviour of the lubricants in hydrodynamic contacts.Many microcontinuum theories have been developed by Ariman and Sylvester [4,5] and Stokes [6] to describe the rheological behaviour of such non-Newtonian lubricants.Among these, the Stokes microcontinuum theory is considered as the simplest generalisation of the classical continuum theory of fluids, which allows the polar effects such as the presence of couple stresses and body couples.Several investigators such as Ramanaiah [7], Bujurke and Naduvinamani [8], and Alyaqout and Elsharkawy [9] considered Stokes microcontinuum theory for the study of various bearing systems because of its relative mathematical simplicity.In literature, investigations on the couple stress squeeze films between flat circular disks between a sphere and a flat plate and between porous surfaces have been carried out by Lin and Hung [10], Lin [11], and Naduvinamani et al. [12], respectively.However, the curved plate geometry is representative of many practical machineries involving hydrodynamic lubrication.It is observed that in a squeeze flow problem, the gap width between the approaching surfaces is small, when compared with the dimensions of the surfaces.Thus even if the order of waviness in the curved plates is small, the effects of nonflatness on the squeeze film flow is not necessarily negligible 2 Advances in Tribology as shown by Gupta and Kapur [13], Usha and Vimala [14,15].Therefore, the effects of curvature on the squeeze film behaviour are also significant in this study.
The study of variation of viscosity with pressure on the isothermal, incompressible squeeze films between parallel disks by Butler [16] shows that the predicted fluid pressure is higher when VPD is taken into consideration than in the isoviscous case.The effects of variation of viscosity with pressure and temperature have been studied by Gould [17] and it has been shown that the above effects play a major role in squeeze film behaviour.Kalogirou et al. [18] have studied VPD in incompressible Poiseuille flow of Newtonian liquids.Lu and Lin [19] have studied the combined effects of couple stresses and VPD for the sphere-plate squeeze film system and observed that these effects enhance the load carrying capacity of the squeeze film.
The circular squeeze film system can be used in biolubricated joints such as synovial joints and clutch plates in automotive transmission.In this paper, a squeeze film flow of couple stress fluids between a flat circular fixed disk and a curved circular moving disk is considered taking VPD into account.Using Stokes microcontinuum theory and the Barus formula, the nonlinear Reynolds type equation is derived.The approximate analytical expressions for the squeeze film pressure and the load carrying capacity are obtained using a perturbation technique.Numerical results are obtained for the sinusoidal motion of the upper curved disk with the curvature described by an exponential form.The effects of couple stresses, curvature, and variation of viscosity with pressure and their combined effects on the squeeze film characteristics are investigated.

Mathematical Formulation
The unsteady, axi-symmetric laminar flow of an incompressible couple stress fluid between a flat circular lower disk at  = 0 and a curved circular upper disk at  = ℎ(, ) of radius   is shown in Figure 1.The upper curved disk approaches the lower flat disk with a squeezing velocity   = −ℎ  = −ℎ/.
Assuming that the body forces and body couples are absent, the governing equations of motion of the incompressible couple stress fluid in the axi-symmetric case using Stokes microcontinuum theory take the form where  and  are the velocity components in  and  directions, respectively,  is the film pressure,  is the fluid density,  is the shear viscosity, and  represents the new material constant with the couple stress property.To consider the variation of viscosity with pressure, the Barus formula for isothermal VPD is where  0 is the viscosity at ambient pressure and at constant temperature and  is the pressure-viscosity coefficient.The boundary conditions for the velocity components are the noslip conditions given by and the no-couple stress conditions given by The boundary conditions for the squeeze film pressure are  = 0, on  =   ,   = 0, at  = 0. (5)

Nondimensional Form of the Problem
Introduce the following nondimensional quantities Advances in Tribology 3 Here  = √/ is the couple stress parameter of dimension length and this length may be regarded as the chain-length of polar additives in the lubricant.Using ( 2) and ( 6), the nondimensional form of the governing equations given by ( 1), take the form where the nondimensional form of squeeze film pressure  is given by  = ℎ 3 0 / 2    and the viscosity parameter is  =  0  2    /ℎ 3 0 .On nondimensionalization, the boundary conditions of the velocity components and pressure become Solving (8) with the boundary conditions ( 10) and ( 11), the radial velocity component is obtained as Integrating the continuity equation ( 7) with respect to  and using the boundary conditions of axial velocity component from ( 10) and ( 11), the non-Newtonian Reynolds type equation is obtained as where (, , , ) =  3 exp(−) − 12 2  exp(−2) + 24 3 exp(−2.5)tanh( exp(0.5)/2).
To find the approximate analytical solution for the highly nonlinear equation ( 14), the squeeze film pressure is expanded, using a small perturbation, in the form with 0 ≤  ≪ 1. Substituting (15) in the Reynolds type equation ( 14) and in the pressure conditions the zero order problem becomes where  0 (, ) =  3 − 12 2  + 24 3 tanh(/2).The pressure conditions are given by By solving (16) along with the boundary conditions given by ( 17), the zero order solution is obtained as Similarly, the first order problem becomes where The first order solution is obtained as where  1 and  2 are dummy variables of integration representing the radial coordinate .Substituting  0 and  1 in (15), the approximate squeeze film pressure  is obtained as The dimensionless load carrying capacity is now obtained by integrating the film pressure given in (22) as where  1 ,  2 , and  3 are dummy variables of integration representing  and the nondimensionalization for the load carrying capacity  sq is given by  sq =  sq ℎ 3 0 / 0  4    .Equations ( 22) and (23) give the squeeze film pressure and the load carrying capacity for an arbitrary shape of the upper curved disk.

Solution for Sinusoidal Motion
Assuming a specific form of ℎ(, ) given by ℎ(, ) = ℎ 1 ()() the problem is analysed.Here ℎ 1 () represents the central film thickness and () is the curvature profile of the curved disk.The solution is obtained for a sinusoidal motion ℎ 1 () = ℎ 10 +  sin  of the upper curved disk.Here ℎ 10 is the initial central film thickness and  and  are the amplitude and the angular frequency of the sinusoidal motion, respectively.Here the curvature profile is assumed to take an exponential form () =  −ĉ 2 where ĉ is the curvature parameter with ĉ > 0 for concave disk, ĉ < 0 for convex disk, and ĉ = 0 for a flat disk.Using the nondimensional quantities from (6) along with the dimensionless form of ℎ(, ), ℎ 1 (), and () is given by where  = ĉ 2  .The configuration of the curved disk for various values of the curvature parameter  is shown in Figure 2. The dimensionless squeeze film pressure and force can be obtained from ( 22) and (23) using (25).Considering the flat disk case, the gap-width becomes independent of the radial coordinate; that is, ℎ(, ) = ℎ().In this case, the squeeze film pressure and the load carrying capacity are obtained from ( 22) and (23) as

Results and Discussions
The combined effects of curvature, couple stresses, and variation of viscosity with pressure on the squeeze film behaviour are analysed using Stokes microcontinuum theory and Barus formula for VPD.The effects of couple stresses on the squeeze film characteristics are observed through the nondimensional couple stress parameter , the curvature effects through the curvature parameter , and the variation of viscosity with pressure through the viscosity parameter .
Considering the sinusoidal motion of the upper curved disk, the squeeze film characteristics for the flow in the flat disk case are obtained, respectively, from ( 26) and ( 27).Here, a small amplitude  = 0.2 is taken for the sinusoidal motion and time  = 0.8 is chosen for a squeezing state.Figures 3-5 present the squeeze film characteristics for the flat disk case.Figure 3 shows the squeeze film pressure  as a function of radial coordinate  for a small amplitude of sinusoidal motion  = 0.2 at time  = 0.8 and for the viscosity parameter  = 0.02 and for various values of the couple stress parameter .Higher pressure is obtained for larger values of  due to the rotations of higher number of microstructures present in the fluid.Figure 4 shows  as a function of  for  = 0.2,  = 0.8, and  = 0.2 and for different values of .Higher film pressure is obtained by considering VPD ( = 0.02), compared to the isoviscous case ( = 0), and an increase in the squeeze film pressure is observed with an increase in the viscosity parameter .show the pressure graphs for  = 0.8,  = 0.02 and various values of  in the concave ( > 0) and convex ( < 0) disk cases, respectively.It is observed that in both the cases, the squeeze film pressure increases with an increase in the couple stress parameter .Further, it is seen that the squeeze film pressure is more significant for larger values of  for the concave nature of the upper curved disk.Figures 7(a) and 7(b) display the squeeze film pressure  as a function of radial coordinate  with  = 0.8,  = 0.2 and for various values of  in the concave and convex disk cases, respectively.In both concave and convex disk cases, higher squeeze film pressure is observed by considering VPD with  = 0.02, compared to the isoviscous case with  = 0 and a further increase in the squeeze film pressure is observed with an increase in the viscosity parameter , particularly in the concave disk case.Figures 8(a) and 8(b) describe the variation of load carrying capacity  sq as a function of response time  for  = 0.2 and various values of  and  in the concave ( > 0) and convex ( < 0) disk cases, respectively.An increase in the load carrying capacity is obtained with an increase in the couple stress parameter .Further increase in the load carrying capacity is observed when considering VPD for both concave and convex nature of the upper curved disk.Figure 9 shows the variation of load carrying capacity as a function of response time at  = 0.2,  = 0.02 and various values of curvature parameter .It is observed that the load carrying capacity is higher in a concave disk than in a flat disk whereas it is lower in a convex disk than in a flat disk.

Conclusion
The combined effects of couple stresses and VPD on the curved squeeze film behaviour are studied through the squeeze film pressure and load carrying capacity.Using Stokes microcontinuum theory and the Barus formula, the Reynolds type equation is derived.The expressions for squeeze film pressure and load carrying capacity are derived using a perturbation technique.It is observed that the effects of couple stresses, variation of viscosity with pressure, and curvature and their combined effects strongly influence the squeeze film behaviour.That is, the larger the value of the couple stress parameter, the greater the load carrying capacity in the squeeze film.This effect is more significant when the viscosity pressure dependence (VPD) is considered.The magnitude of the squeeze film pressure and load carrying capacity on a concave disk is always more than that on a flat disk whereas it is always less for a convex disk which suggests that the load carrying capacity in the squeeze film is strongly influenced by the curvature of the upper moving disk in addition to VPD and couple stress parameters.On the whole the combined effects of curvature VPD and couple stresses on the squeeze film pressure and load carrying capacity are seen to be more significant in a concave disk.Thus, a concave disk is more suitable for practical applications of squeeze film flows with couple stress fluids.

Figure 5 :Figure 6 :
Figure 5: Effects of couple stresses and VPD on the squeeze film force-flat disk case.

Figure 5
Figure 5 describes the variation of the load carrying capacity  sq as a function of response time  for  = 0.2 and for various values of  and  in the flat disk case.Compared to the isoviscous Newtonian case with  = 0 and  = 0, an increase in the load carrying capacity is observed with an increase in the couple stress parameter with  = 0.1 and  = 0. Further increase is seen when considering VPD with  = 0.02 and  = 0.1.Combined effects of couple stresses and VPD, with  = 0.2 and  = 0.04, further enhance the load carrying capacity.Figures 6-9 present the squeeze film characteristics of the flow in the curved disk case with the same amplitude  = 0.2 of the sinusoidal motion.Figures 6(a) and 6(b)

Figure 7 :
Figure 7: Effects of VPD on the squeeze film pressure.(a) Concave disk and (b) convex disk.

Figure 8 :Figure 9 :
Figure 8: Effects of couple stresses and VPD on the squeeze film force.(a) Concave disk and (b) convex disk.