Medical, Pharma, Engineering, Science, Technology and Business

**Anuradha Parinam ^{*}, Ram Karan**

Department of Mechanical Engineering, University Institute of Engineering and Technology,Kurukshetra University, Kurukshetra, INDIA

- *Corresponding Author:
- Anuradha Parinam

Department of Mechanical Engineering

University Institute of Engineering and Technology

Kurukshetra University, Kurukshetra, INDIA

**E-mail:**[email protected]

**Visit for more related articles at** International Journal of Advancements in Technology

Most of the available works in elastohydrodynamic lubrication (EHL) have been carried out considering exponential pressure-viscosity relationship. However, recently, some lubricants such as 2, 3-dimethylpentane, pure hydrocarbon, water/glycol solutions etc. are reported to show linear piezo-viscous response. The investigations carried out using exponential pressure-viscosity relation for such lubricants yield to large errors as they highly overestimate the film thickness values. Also, a major aspect of practical EHL behavior i.e. starvation effect still remains largely unexplored. Starvation refers to the inadequate filling of the conjunction and has a considerable effect on EHL characteristics such as film thickness and coefficient of friction. Therefore, a numerical analysis is carried out for lubricants exhibiting linear piezo-viscosity to analyze the influence of starvation in EHL line contacts. Starvation effect is modeled by shifting the inlet meniscus towards the contact zone. The Reynolds equation is discretized using finite differences and solved along with the load balance equation using Newton-Raphson technique.

*EHL, Linear piezo-viscosity, Starvation, Line contact.*

EHL is a form of hydrodynamic lubrication in which elastic deformation of the surfaces and pressure dependence of viscosity play a fundamental role in the generation of load carrying film, which is thick enough to separate the surfaces completely. The accurate prediction of lubricant film thickness and friction coefficient in EHL contacts involved in mechanical components such as gears, cams, rolling element bearings etc. has long been the focus of attention. Conventionally, EHL studies employed piezo-viscous relationships involving an exponential function of lubricant pressure. Dowson and Higginson [1] and Hamrock and Dowson [2] presented film thickness equations using exponential piezo-viscous response and Newtonian fluid model for EHL line and point contacts, respectively. However, recently, Kumar et al. [3] presented examples of some practical situations involving high operating temperatures, low viscosity lubricants and water/glycol solutions where the lubricant’s piezo-viscous response is more closely described by a linear relationship. Low molecular weight liquids like a pure hydrocarbon, 2, 3-dimethylepentane [4], octane and toluene [5] exhibit a linear piezo-viscous response. In addition, concentrated contact lubrication using the process liquid of low molecular weight exhibits a linear pressure viscosity relationship and also, same trend is observed for ordinary lubricants at very high temperatures [3]. Further, it was reported that the viscosity pressure relationship of an aircraft engine lubricants (eg. L7808 and L23699) is linear. In fact, Barus’s coefficient (α_{B} ) decreases linearly as temperature increases for water/glycol solution. However, water/glycol solution forms a thinner film than conventional lubricants under the same operating conditions [6]. Furthermore, a substantial difference between the piezo-viscosities (and hence, EHL behaviors) represented by a simple exponential and a linear relationship was observed [7]. The use of conventional film thickness formulas for such lubricants may lead to large errors in the estimation of film thickness values. Therefore, using the linear piezo-viscous relationship of Barus [8-9], Kumar et al. offered central film thickness formulas which were validated experimentally. Anuradha and Kumar [7] extended their work by developing central and minimum film thickness formulas for EHL line contacts using linear piezo-viscosity.

Most of the above works were accomplished by considering fully flooded conditions. However, sometimes contact zone suffers from shortage of lubricant supply which leads to thinning of lubricant film. This is known as “starvation”. Moreover, due to inadequate filling of the conjunction, the EHL film generated is much thinner as compared to the corresponding fullyflooded film and therefore, traction coefficient is also very high. Starvation effect is found to be higher at high speeds, using high viscosity lubricants and limited lubricant supply [10]. The starved steady state EHL circular contact problem was studied by Chevalier, Lubrecht and Cann [10, 11] for grease lubricated contacts. The transient effect in EHL circular contacts was investigated by Dumont et al. [12] by introducing the surface feature effects such as micro-pitted surface. Venner et al. [13] performed the numerical simulations to study the deformation of transverse and isotropic harmonic waviness in EHL circular contacts under pure rolling in relation to lubricant supply to the contact. Kumar and Khonsari [14] studied the effect of starvation on traction and film thickness in Thermo-EHL line contacts with shear-thinning lubricants. Therefore, the present work focuses on the effect of starvation on film thickness and traction coefficient using lubricant with linear piezo-viscosity with Newtonian fluid model.

*Reynolds Equation*

The Reynolds equation which governs the generation of pressure in the EHL line contacts with compressible fluid in dimensionless form is:

---- (1)

*Film Thickness Equation*

The film thickness equation in dimensionless form is:

---- (2)

where D_{ij} are the influence coefficients for a uniform mesh size ΔX :

---- (3)

*Boundary Conditions*

The boundary conditions applicable are:

P = 0 at X= *X _{in}*---- (4)

---- (5)

*Density-Pressure Relationship*

The following Tait’s equation of state is used to describe density-pressure relationship:

---- (6)

The constants K_{o} and K '_{o} for the present analysis are taken as K_{o}= 10 and K_{o} =1.0 GPa (same as used in Ref [7]).

*Viscosity-Pressure Relationship*

The present analysis uses linear form of piezo-viscous relation as given by Barus:

---- (7)

where, α_{B}=linear pressure viscosity coefficient.

**Load Equilibrium Equation**

The load equilibrium equation in dimensionless form is:

---- (8)

**Percentage of Starvation**

Percentage of starvation is the ratio of difference between film thickness in fully flooded condition and starved condition to the film thickness in fully flooded condition. Percentage of starvation is given as:

---- (9)

*Coefficient of Friction*

Coefficient of friction is the ratio of shear force of lubricant over whole surface to the applied normal load. This is expressed in non-dimensional form as:

---- (10)

where is the non-dimensional shear stress.

The solution procedure consists of following steps as given below.

1. The solution begins by assuming an initial guess for the pressure distribution[ P_{i} ] , offset film thickness H_{o} and outlet boundary coordinate[X_{o} ]

2. The film profile [H_{i} ] is calculated using the pressure distribution [ P_{i} ] in the film thickness Eq. (2).

3. The fluid properties ( ) are computed using Eqs. (6) and (7).

4. The Reynolds Equation (1) is discretized using a mixed second order central and first order backward differencing scheme to obtain the equations as follows:

---- (11)

5. The load equilibrium Equation (8) is written in the following discrete form using Simpson’s coefficients:

---- (12)

---- (13)

6. The simultaneous system of N equations represented by discretized Reynolds Equation (11) and load equilibrium Equation (12) are solved using the Newton-Raphson technique similar to that illustrated by Kumar and Khonsari [15]. The N system unknowns are P_{2} , P_{3} , P_{4} ,……, P_{N} and H_{0} . The matrix equation of this system is

---- (14)

7. The pressure distribution [ P _{i}] and offset film thickness are updated using the correction calculated above.

8. The steps 2-7 are repeated until the following convergence criteria are satisfied:

When inlet conjunction does not fill adequately with lubricant then there is shortage of supply of lubricant in contact area which leads to film thinning and even breakdown takes place. This is called starvation effect. In the present analysis, the starvation effect is investigated by shifting the position of inlet meniscus towards contact zone. The effect of starvation is observed at different value of loads, velocities, pressure-viscosity coefficients and equivalent radii.

Fully-flooded condition is considered to be at the position of inlet meniscus where negligible change in film thickness is observed. In **Table 1**, central film thickness (nm) at three different loads is shown with change in inlet meniscus position from *X _{in}*=-7 to

It is observed that *X _{in}* =-5 ensures the fully-flooded condition for different operating conditions. At this position, it is assumed that the effect of starvation is zero and inlet meniscus occupies enough lubricating oil. So, effect of starvation is studied by shifting the inlet meniscus position from

* Percentage of Starvation*

Percentage of starvation is defined as the ratio of difference of film thickness in fullyflooded and starved condition to the fully-flooded conditions. In this section, starvation effect is investigated considering different loads, rolling speeds, pressure viscosity coefficients and equivalent radii. The range of input parameters used in the simulations is shown in **Table 3**.

The variation of percentage of starvation is shown in **fig. 1** for three different loads (*W*=3x10^{-5}, 1x10^{-4}, 5.2x10^{-4}). Percentage of starvation varies with the position of inlet meniscus (*X _{in}*). As apparent from

**Fig. 2** shows the variation of percentage of starvation with the position of inlet meniscus at different rolling speeds (U=10^{-12}, U=10^{-11}, U=10^{-10}). The starvation is zero at *X _{in}* =-5 due to the fully-flooded condition. Unlike the load, as inlet meniscus position is shifts towards contact zone a rapid increase in starvation is found for higher speed (U=10

*Pressure Profile*

Pressure distribution for two different inlet meniscus positions is shown in **figures 3** and 4 for two different loads W=3x10^{-5},and W=5.2x10^{-4}, respectively. One of inlet meniscus position is at fully-flooded condition (*X _{in}*=-5) and second (

It can be seen from **fig. 3** that significant pressure builds up in the inlet zone for fullyflooded inlet meniscus position. However, only a marginal difference is observed in the central pressure values. On the other hand, the pressure profiles almost coincide at the higher value of load as seen in **fig. 4**.

**Figures 5** and **6** show the effect of rolling speed on pressure distribution. It can be seen that at the lower speed (U=10^{-12}), the pressure profiles are almost coincident while at higher speed (U=10^{-10}), the pressure build up starts early for the fully flooded condition.

* Film Profile*

**Figures 7** and **8** show the film profiles at two different positions of inlet meniscus keeping other variables constant, at two different loads W=3x10-5 and W=5.2x10^{-4} respectively. It can be observed from both these figures that the film thickness obtained at *X _{in}* =-5 is greater than the film thickness at

This is due to the starvation effect which is much more near the contact zone. However, it can be seen from **fig. 10** that at higher value of load, the film thickness profile at two different position of inlet meniscus has a little difference. This shows that the effect of starvation is less pronounced at higher load at both inlet meniscus positions. This is due to the fact that the pressure distribution in contact zone approaches to Hertzian pressure and shrinkage of inlet zone takes place which minimize the effect of starvation.

The effect of rolling speed on the film profile is shown in **fig. 9** and **10** for two different positions of inlet meniscus at U=10^{-12} and U=10^{-10} respectively. It can be observed that unlike the load, the effect of starvation is more pronounced at higher rolling speed. Film thickness is much thinner at *X _{in}* =-1.06 in comparison with the film thickness at

As the lubricating film needs to be thicker at such high speeds, excess lubricant needs to be supplied which lead to shortage of lubricant at inlet zone.

*Coefficient of Friction*

**Fig. 11** shows the variation in COF_{Starved}/COF_{Fully-Flooded} with the percentage of starvation (φ) at three different loads, W=3x10-5, W=1x10^{-4}, W=5.2x10^{-4}. Coefficient of friction rises with rise in percentage of starvation. It can be seen from the figure that there is a steep increase in the friction value at lower load (W=3x10-5) in comparison with other two loads. This is because high starvation at lower load, the film gets thinner and offers much friction to the moving surfaces which increase the coefficient of friction.

**Fig. 12** compares the ratio of coefficient of friction with the percentage of starvation at three different speeds, U=10^{-12}, U=10^{-11}, U=10^{-10}. There is continuous rise in coefficient of friction with the increase in percentage of starvation at all speeds. However, a huge rise in coefficient of friction is noticed for higher speed, U=10^{-10}. This is attributed to high starvation effect at high rolling speed as discussed in previous sections. Due to high starvation, contact zone suffers from shortage of lubricant supply and film gets thinner to large extent and rubbing of surface takes place which increase the coefficient of friction.

A comprehensive EHL analysis is carried out to study the starvation effect on pressure distribution, film profile and coefficient of friction by changing the position of inlet meniscus for a lubricant exhibiting linear piezo-viscosity. It is found that starvation has a great influence on the film thickness and coefficient of friction. Starvation is largely influenced by rolling speed and viscosity of lubricant. Following are some of the salient conclusions:

1. Starvation is negligible in fully-flooded conditions and increases continuously while shifting inlet meniscus toward contact zone. It has been observed that it has maximum value near contact zone.

2. Starvation increases rapidly with increase in rolling speed. However, starvation effect decreases with increase in load. It is due to the fact that the pressure distribution in contact zone approaches to Hertzian pressure and shrinkage of inlet zone takes place which minimize the effect of starvation.

3. It has been observed that both pressure distribution and film profile are not affected by starvation at high load and low velocity. However, pressure profile has a significant difference for fully-flooded conditions and starved conditions at low load and high velocity, which is attributed to lower starvation at high load and low speed.

4. Coefficient of friction increases rapidly with increase in rolling speed. This is due to the fact that starvation is higher at high rolling speeds. Due to high value of starvation the contact zone suffers from lack of lubricant supply and thus, film gets thinner to a large extent and rubbing of surfaces takes place which lead to high value of coefficient of friction.

**Nomenclature**

**Dimensional Parameters**

b : half width of Hertzian contact zone,

E': effective elastic modulus of rollers,

*h* : film thickness (m)

*h _{c}* : central film thickness (m)

h_{min} : minimum film thickness (m)

p : pressure (Pa)

p_{H} : maximum Hertzian pressure,

R : equivalent radius of contact (m).

: average rolling speed,

:velocities of lower and upper surfaces, respectively (m/s)

*w* : applied load per unit length (N/m)

*Greek Symbols*

α : Bridgman (exponential) piezo-viscous coefficient (Pa^{-1})

α_{B} : Barus (linear) piezo-viscous coefficient (Pa^{-1})

: inlet density of the lubricant (kg/m^{3})

: lubricant density at the local pressure and temperature (kg/m^{3})

: low-shear viscosity at p=0 (Pa.s)

μ : limiting low-shear viscosity at local pressure (Pa.s)

**Dimensionless Parameters**

G : classical dimensionless load parameter,

G : dimensionless load parameter,

H : dimensionless film thickness,

H_{o} : dimensionless offset film thickness

K_{o} : isothermal bulk modulus at *p* = 0 (Pa)

: pressure rate of change of isothermal bulk modulus at *p *= 0

*P* : dimensionless pressure,

*U* : dimensionless speed parameter,

*W* : dimensionless load parameter,

*X* : coordinate in the direction of surface velocity (m),

*Greek Symbols*

: dimensionless fluid density,

: dimensionless viscosity, o

[1] Dowson, D. and Higginson G.R., “Elastohydrodynamic Lubrication”, Oxford: Pergamon Press, 1966.

[2] Hamrock B.J. and Dowson D., “Ball Bearing Lubrication–The Elastohydrodynamics of Elliptical Contacts”, New York: John Wiley & Sons, 1981.

[3] Kumar, P., Bair, S., Krupka, I. and Hartl, M., “Newtonian Quantitative Elastohydrodynamic Film Thickness with Linear Piezo-Viscosity,” Tribology International, 43, 2010, pp. 2159-2165.

[4] Pensado, A., Comunas M.J.P., Lugo L., and Fernandez J., “Experimental dynamic viscosities of 2, 3- Dimethylpentane up to 60 MPa and from (303.15 to 353.15) K using a rolling-ball viscometer”. Journal of Chemical Engineering Data, 50, 2005, pp. 849-855.

[5] Harris, K. R., Malhotra, R. and Woolf, L. A., “Temprature and Density Dependence of the Viscosity of Octane and Toluene,” Journal of Chemical Engineering Data, 42, 1997, pp. 1254-60.

[6] Spikes, H.A., “Wear and Fatigue Problems in Connection with Water Based Hydraulic Fluids,” Journal of Synthetic Lubrication, 4(2), 1987, pp. 115-35.

[7] Anuradha, P. and Kumar, P., “EHL Line Contacts Central and Minimum Film Thickness Equations for Lubricants with Linear piezo-viscous Behavior” Tribology International, 44, 2011, pp. 1257-1260. [8] Barus, C., “Notes on the Dependence of Viscosity on Pressure and Temprature,” Proc American Arts Science, 27, 1891, pp.13-8.

[9] Barus, C., “Isothermal, Isopiestics and Isometrics Relative to Viscosity” Am Journal of Science, 266, 1893, pp. 87-96.

[10] Chevalier, F., Lubrechet, A.A., Cann, P.M.E., Colin, F., and Dalmaz, G., “Starvation Phenomena in EHL Point Contacts: Influence of Inlet Flow Distribution”, Proceedings of 22nd Leeds-Lyon Symposium on Tribology, Elseviers Tribology Series, 31, 1995, pp. 213-223.

[11] Chevalier, F., Lubrechet, A.A., Cann, P.M.E., Colin, F., and Dalmaz, G., “Film Thickness in Starved in Starved EHL Point Contacts”, ASME Journal of Tribolgy, 120, 1998, pp. 126-133.

[12] Dumont, M.L., Lugt, P.M., and Tripp, J.H., “Surface Feature Effects in Starved Circular EHL Contacts”, ASME Journal of Tribolgy, 124(2), 2001, pp. 358-366.

[13] Venner, C.H., Berger, G., and Lugt, P.M., “Waviness Deformation in Starved EHL Circular Contacts”, ASME Journal of Tribolgy, 126(2), 2004, pp. 248-257.

[14] Kumar, P. and Khonsari, M.M., “Effect of Starvation on Traction and Film Thickness in Thermo-EHL Line Contacts with Shear-Thinning Lubricants”, Tribology Letters, 32, 2008, pp. 171-177.

[15] Kumar P. and Khonsari M.M., “Combined Effects of Shear Thinning and Viscous Heating on EHL Characteristics of Rolling/Sliding Line Contact”, ASME Journal of Tribolgy, 130, 2008, pp. 041505-1-13.

Select your language of interest to view the total content in your interested language

- ACS Nano
- AI Algorithm
- Advanced Materials
- Animation
- Antenna Design
- Architect
- Architectural Drawing
- Architectural Engineering
- Architecture
- Artificial Intelligence
- Automated Mining
- Automation
- Behavior-based systems
- Bioinformatics Techniques
- Biomechanics
- Biomechanics and Robotics
- Bionics
- Building design
- Communication Systems
- Computational Neuroscience
- Computer Networking
- Computer Networks
- Concrete
- Construction
- Construction Engineering
- Construction Estimating Software
- Cordless
- Cyber Physical System
- Cyber Security
- Cyber Technologies
- Data Engineering
- Datamining
- Digital Image Processing
- Digital Signal Processing
- Engineering Drawing
- Fabric Formwork
- Fluid Engineering
- Fuzzy Logic
- Gaming Technology
- Genomics
- Graphic design
- Graphic designing
- Image Processing
- Industrial Robotics
- Information Science
- Information Theory
- Intelligent Robotics
- Interior Design
- Interior Designing
- Landscape Architecture
- Lovotics
- Machine Vision
- Manufacturing
- Medical Device
- Medical Robotics
- Mobile Device
- Mobile Repots
- Modeling
- Motion picture
- Multimedia Network
- Neural Networks
- Neurorobotics
- Product Modeling
- Radar Imaging
- Radio Waves
- Robotic Rehabilitation
- Satellite Communication
- Search Engine Design
- Sensors
- Simulation
- Simulation and Visualization
- Simulation and Visualization Effects
- Smart Grid
- Social Robots
- Sociology of Architecture
- Soft Computing
- Spatial Database
- Telecommunications Journal
- Telerobotics
- Tribology
- Urban Design
- Urban Planner
- Virtual Reality
- Virtual Reality Technology
- Wave Length
- Web Developer
- Wireless Journals
- swarm intelligence and robotics

- International Conference on Aerospace & Aerodynamics, August 02-03, 2018 Barcelona, Spain

August 02-03, 2018 Barcelona, Spain - Big Data Innovation, Data Mining and Analytics Summit, August 20-21, 2018 Holiday Inn Singapore Atrium, Singapore

August 20-21, 2018 Singapore City, Singapore - International Conference on Artificial Intelligence, Robotics & IoT, August 21-22, 2018 Paris, France

August 21-22, 2018 Paris, France

- Total views:
**11779** - [From(publication date):

July-2014 - Jul 20, 2018] - Breakdown by view type
- HTML page views :
**8014** - PDF downloads :
**3765**

Peer Reviewed Journals

International Conferences 2018-19