# Non-Newtonian Effects on the Squeeze Film Characteristics between a Sphere and a Flat Plate: Rabinowitsch Model

Authors: Udaya P. Singh; Ram S. Gupta
Journal: Advances in Tribology
 (2012)

---

## Abstract

The use of additives (polyisobutylene, ethylene-propylene, lithium hydroxy stearate, hydrophobic silica, etc.) changes lubricants’ rheology due to which they show pseudoplastic and dilatant nature, which can be modelled as cubic stress fluid model (Rabinowitsch fluid model). The present theoretical analysis investigates the effects of non-Newtonian pseudoplastic and dilatant lubricants on the squeezing characteristics of a sphere and a flat plate. The modified Reynolds equation has been derived and an asymptotic solution for film pressure is obtained. The results for the film pressure distribution, load carrying capacity, and squeezing time characteristics have been calculated for various values of pseudoplastic parameter and compared with the Newtonian results. These characteristics show a significant variation with the non-Newtonian pseudoplastic and dilatant behavior of the fluids.


---

## Body

## 1. Introduction

Squeeze film between a sphere and plate is observed in various machine elements such as ball bearings, cam and followers, and gears. The mechanical action (squeezing, shearing, etc.) leading to the generation of high pressure at the contacts [1–4] changes the rheology of the lubricants such as viscosity and density which account for the performance characteristics of machine elements. Dowson [5], Wada and Hayashi [6], and Yadav and Kapur [7] emphasized the variation of viscosity and density with temperature and pressure and reported significant changes in bearing characteristics. Denn [8], Rajagopal [9], and Renardy [10] indicated that in high pressure lubrication applications, the variation of viscosity becomes more important than the density. Variation of the viscosity also causes the instability of the lubricants’ nature by changing its shearing stress-strain rate relation due to which the estimated characteristics of lubricated contacts such as sphere-plate contacts (point contacts) may deviate from the desired value. This situation is avoided by enhancing the efficiency of stabilizing properties of lubricants by the addition of additives (polyisobutylene, ethylene propylene, etc.). The use of additives minimizes the sensitivity of the lubricant to the change in the shearing strain rate and the lubricants behave like non-Newtonian pseudoplastic, dilatant, and viscoplastic fluids depending on the nature and quantity of the additives. To account for the effects of lubricant additives on the performance characteristics of lubricated point contacts, various non-Newtonian fluid models like power law, micropolar, and couple stress fluid models have been studied by researchers from time to time [11–14]. Among these fluid models, Rabinowitsch fluid model [6] is an established model to predict the effects of additives on the performance characteristics of the lubricated bearings. The shearing stress-strain relation in this model for one-dimensional fluid flow is given by
(1)τ-rz+κ-τ-rz3=μ-∂u-∂z-,
where μ- is the initial viscosity of lubricant and κ- is the nonlinear factor responsible for the non-Newtonian effects of the fluid which will be referred to as the coefficient of pseudoplasticity. This model can be applied to Newtonian, dilatant, and pseudoplastic lubricants for κ-=0, κ-<0, and κ->0, respectively. The advantage of this model lies in the fact that the theoretical analysis for the present model was verified with the experimental justification by Wada and Hayashi [6]. They used spindle oil as the base lubricant and concluded a decrease of dimensionless film pressure with the increase of additive (polyisobutylene). Afterwards, the theoretical study of bearing performance with non-Newtonian lubricants using this models was done by Bourgin and Gay [15] on journal bearing, Hashimoto and Wada [16] on circular plates bearing, and Lin [17] on parallel circular plates. Recently Singh et al. [18–21] used this model to study the performance characteristics of hydrostatic thrust bearings and slider bearings.The objective of this paper is to extend the results [18–21] to squeeze film characteristics between a sphere and a plate by introducing a quantitative analysis using Rabinowitsch fluid model, which accounting for the effect of additives in the lubricant. The modified Reynolds equation governing the squeeze film pressure is derived. Squeeze film characteristics such as film pressure, load-carrying capacity, and squeezing time are presented. The importance of the present analysis lies in that the earlier theoretical investigations on sphere-plate squeezing [11–14] are based on couple stress or micropolar fluid models, which suffer the scarcity of experimental verification.


## 2. Constitutive Equations and BoundaryConditions

The physical configuration of a sphere-plate system is shown in Figure1. The sphere is approaching towards the plate with a normal velocity (-dh-/dt-), separated by a lubricant thin film. The lubricant in the system is taken as a non-Newtonian Rabinowitsch fluid. The body forces and body couples are assumed to be absent.Figure 1
Schematic diagram of squeeze film between a sphere and a plate.Under the assumptions of hydrodynamic lubrication applicable to thin film as considered by Dowson [5], the field equations governing the one-dimensional motion of an incompressible non-Newtonian Rabinowitsch fluid model in polar coordinates (r,θ,z) system [18] are
(2)1r-∂∂r-(r-u-)+∂w-∂z-=0,(3)∂p-∂r-=∂τ-rz∂z-,(4)∂p-∂z-=0
which are solved under the no-slip boundary conditions:
(5)u-=0,w-=0,atz-=0,u-=0,w-=-∂h-∂t-,atz-=h-,
where u- and w- are the velocity components in r- and z- directions, respectively, and h- is the film thickness between the sphere and plate.


## 3. Analysis

Integrating (3) with respect to z- under boundary conditions (5) and using (1), the expression for velocity u- is obtained as
(6)u-=12μ-[∂p-∂r-z-(z--h-)+κ-(∂p-∂r-)3×(12z-4-z-3h-+34z-2h-2-14z-h-3){(∂p-∂r-)3}].
Integrating (2) with respect to z- under the relevant boundary conditions (5) for w- and using (6), the modified Reynolds equation is obtained as
(7)1r-∂∂r-[r-{h-3∂p-∂r-+3κ-h-520(∂p-∂r-)3}]=-12μ-∂h-∂t-.
In the limiting case of κ→0, (7) reduces to the Newtonian form of Reynolds equation obtained by Conway and Lee [22]:
(8)1r-∂∂r-[r-h-3∂p-∂r-]=-12μ-∂h-∂t-.
The expression for film thickness between the sphere and plate at a time t- is taken of the form [1]
(9)h-=h-m+r-22R,
where R denote the radius of the sphere. In case of squeezing between two spheres, the value of the radius R can be taken as (i) R-1=R1-1+R2-1 for external contact and (ii) R-1=R1-1-R2-1 for internal contact, where R1 and R2 are the radii of the spheres.The modified Reynolds equation (7) takes the dimensionless form:
(10)1rddr[r{h3dpdr+3αh520(dpdr)3}]=-12β,
where p=h-mo2p¯/μ-R(dh-m/dt-) is the dimensionless pressure, α=κ-(μ-R(dh-m/dt-)/h-mo2)2 is the parameter of pseudoplasticity, β=h-mo/R is the sphere parameter, hm=h-m/h-mo is the minimum film thickness at a time t, and h-mo is the minimum film thickness at t-=0. The value of the pseudoplastic coefficient κ- depends on the type of and the quantity of additives which can be determined experimentally [6]. Thus, the values of R, h-mo, and μ- being known for a particular bearing and lubricant and the  vales of α can be calculated with the appropriate value of κ-. However, for the validity of the present analysis, the value of α is restricted to |α|<0.01.As (10) is a nonlinear equation in p, it is not easy to solve it using analytical methods. Therefore, the classical perturbation method is used to solve it. The perturbation series for p can be expressed in the form:
(11)p=po+αp1+α2p2+⋯.
For α≪1, it is sufficient, for analysis, to consider the first order term in α as follows:
(12)p=po+αp1.
For the higher values of α, second and higher order terms can be considered to increase the accuracy of the results. However, for the higher values of α, it will be more appropriate to adopt a numerical solution procedure such as the finite element method to solve the Reynolds equation.Substituting (12) in (10), the perturbation equations are obtained as
(13)1rddr[rh3dpodr]=-12β,1rddr[r{h3dp1dr+3h520(dpodr)3}]=0.
Solving (13) under the boundary conditions
(14)dpdr=0atr=0,p=0atr=1,
the dimensionless pressure developed in the film region is:
(15)p=12β2[1(1+2hmβ)2-1(r2+2hmβ)2]-345625αβ4[3+hmβ(1+2hmβ)6-3r2+hmβ(r2+2hmβ)6].


### 3.1. Load Carrying Capacity

The load carrying capacity can be obtained by integrating the film pressure over the squeezing film area as follows:(16)w-=2π∫0Rp-r-dr-
which takes the dimensionless form:
(17)w=∫01prdr,
where
(18)w=w-[h-mo2(2πμ-R3(dh-m/dt-))].


### 3.2. Squeezing Time

The squeezing time can be calculated by integrating (18) with respect to t under the condition that hm=1 at t=0 as follows:
(19)t=∫hf1wdh=∫hf1∫01rpdrdhm,
where
(20)t=w-h-mo2πμ-R3t-.


## 3.1. Load Carrying Capacity

The load carrying capacity can be obtained by integrating the film pressure over the squeezing film area as follows:(16)w-=2π∫0Rp-r-dr-
which takes the dimensionless form:
(17)w=∫01prdr,
where
(18)w=w-[h-mo2(2πμ-R3(dh-m/dt-))].


## 3.2. Squeezing Time

The squeezing time can be calculated by integrating (18) with respect to t under the condition that hm=1 at t=0 as follows:
(19)t=∫hf1wdh=∫hf1∫01rpdrdhm,
where
(20)t=w-h-mo2πμ-R3t-.


## 4. Results and Discussions

Based on the Rabinowitsch fluids model, the effects of non-Newtonian rheology on the squeeze-film characteristics between a sphere and a plate are investigated using a dimensionless parameterα which accounts for the non-Newtonian nature of the lubricant, that is, for the induced nature due to the use of additives. The parameters α=0, α<0, and α>0 describe the Newtonian, dilatant, and pseudoplastic lubricants, respectively. For the validity of the analysis, the numerical results for non-Newtonian lubricants are compared with the Newtonian results [12].In order to analyze the non-Newtonian effects of fluids on the squeeze-film performance of sphere-plate system, various squeeze-film characteristics are presented with the following values:(i)
pseudoplastic parameterα=-0.01 to 0.01, [2, 18];(ii)
sphere parameterβ=0.03,0.05, [12, 14].Figure 2 shows the variation of dimensionless film pressure (p) with respect to the dimensionless coordinate (r). It is clear from the figure that the pressure is the maximum at r=0, that is, at the minimum film thickness and decreases towards the outer of the sphere and hence, the analysis obeys the basic theory of film pressure in sphere-plate system. Again, the pressure for dilatant lubricants (α<0) is higher than the pressure for Newtonian lubricants whereas the pressure for pseudoplastic lubricants (α>0) is lower than the pressure for Newtonian lubricants. Further, the pressure increases as α decrease from 0.01 to -0.01. It is observed that the effect of non-Newtonian pseudoplastic and dilatant lubricants produces a remarkable change in the film pressure near r=0, and it decreases towards periphery r=1. For design parameters β=0.05 and hm=0.7, a small value of pseudoplastic parameter α=-0.01 (dilatant fluids) increases the film pressure by nearly 25% at r=0.1 and 3% at r=0.5. For the same design parameters, the value of pseudoplastic parameter α=0.01 (pseudoplastic fluids) decreases the film pressure by nearly 30% at r=0.1 and 2% at r=0.5. It shows that the pseudoplastic and dilatant lubricants produce larger effects with the higher film pressure.Variation of dimensionless pressure with respect to the dimensionless radiusr for different values of pseudoplastic parameter α.
(a)
β
=
0.03
,
h
m
=
1
(b)
β
=
0.05
,
h
m
=
0.70Figure3 shows the variation of dimensionless maximum film pressure (pmax) with respect to the dimensionless minimum film thickness (hm). The effect of dilatant lubricants is observed to increase the maximum film pressure from its value in Newtonian case whereas the effect of pseudoplastic lubricant decreases the maximum film pressure from its value in the Newtonian case.  Furthermore, the maximum pressure increases with the decrease of pseudoplastic parameter α from 0.01 to -0.01. There is also a relative change in the maximum pressure for pseudoplastic and dilatant lubricants, which increases with the decrease of film thickness. For the sphere parameter β=0.05, a small value of pseudoplastic parameter α=-0.01 increases the maximum film pressure by nearly 12% to 25% as the minimum film thickness hm decreases from 1 to 0.7. For the same value of β and hm, the value of pseudoplastic parameter α=0.01 (pseudoplastic fluids) decreases the film pressure nearly 15% to 30%. Therefore, it can be safely said that the higher the film pressure, the greater the change produced by pseudoplastic and dilatant lubricants is.Variation of the dimensionless maximum film pressure with respect to the minimum film thicknesshm for different values of parameter α.
(a)
β
=
0.03
(b)
β
=
0.05Figure4 shows the variation of dimensionless load carrying capacity (W) of the system with respect to the dimensionless minimum film thickness (hm). It is observed that the load capacity obtained with dilatant lubricants is higher than that with Newtonian lubricants, and the load capacity obtained with pseudoplastic lubricants is lower than its value obtained with Newtonian lubricants.  Furthermore, the load capacity increases with the decrease of pseudoplastic parameter α from 0.01 to -0.01. It is also observed that there is a relative change in load capacity obtained with the different values of pseudoplastic parameter, which increases with the decrease of film thickness. For the sphere parameter β=0.05, the effect of dilatant lubricant α=-0.01 increases the load capacity by 10% to 15% as the minimum film thickness hm decreases from 1 to 0.7 and for the same value of β and hm, the effect of pseudoplastic lubricants  α=0.01 decreases the film pressure from 15% to 20%.Variation of the dimensionless  load capacity  with  respect  to  the  dimensionless  minimum film thicknesshm  for values of  parameter  α.
(a)
β
=
0.03
(b)
β
=
0.05Figure5 shows the time (t) elapsed in squeezing the film from its initial thickness hm=1 to a final thickness hm=hf. It is observed that for each value of hf and β, the squeeze time for dilatant lubricants is longer than that with the Newtonian lubricants, whereas, the squeeze time for pseudoplastic lubricants is less than its value in the Newtonian case. Further, the squeeze time increases with the decrease of pseudoplastic parameter from 0.01 to -0.01. This phenomenon can be interpreted as a result of the increase in the film pressure from pseudoplastic to dilatant lubricants. For the sphere parameter β=0.03, the time to squeeze the film to hf=0.7 is increased by nearly 23% with dilatant lubricant α=-0.01 in comparison with the Newtonian case, whereas the same is decreased by nearly 18% for pseudoplastic lubricants. For the sphere parameter β=0.05, the time to squeeze the film from h=1 to h=0.7 is increased by nearly 14% with dilatant lubricant α=-0.01 in comparison with the Newtonian case. For the same value of β and hf, the squeeze time is reduced by nearly 16% with pseudoplastic lubricants  α=0.01. Thus, dilatant lubricants increase and pseudoplastic lubricants reduce the squeeze time of the bearing.Variation of dimensionless squeeze time with   respect to the squeezed film thicknesshf  for different values of parameter α.
(a)
β
=
0.03
(b)
β
=
0.05


## 5. Conclusions

Based on the Rabinowitsch fluid model (cubic stress model) for non-Newtonian pseudoplastic and dilatant fluids, the effects of lubricant additives on the performance characteristics of squeezing film between a sphere and a plate are presented avoiding the inertia and cavitation effects. The analytical solution for pressure distribution is obtained using a classical perturbation technique. Based on the present theoretical analysis, the following results have been drawn.(1)
Dilatant lubricants increase the pressure and load carrying capacity significantly, whereas the case is reversed with the pseudoplastic lubricants.(2)
On comparing with the Newtonian case, dilatant lubricants increase the squeeze time, whereas the pseudoplastic lubricants decrease it.(3)
As the squeezing time of the sphere-plate system is significantly increased with the dilatant lubricants, it is expected that the use of additives can reduce the vibration in the sphere-plate systems.Thus, the present analysis can also provide a guideline to control the vibration in the system. Hence, the results are expected to be more helpful for better bearing performance and stability. However, an experimental validation of these results is required at laboratory level.



---

*Source: 571036-2012-12-18.xml*