The unresolved definition of the pressure-viscosity coefficient

In the classical approach to elastohydrodynamic lubrication (EHL) a single parameter, the pressure-viscosity coefficient, quantifies the isothermal pressure dependence of the viscosity for use in prediction of film thickness. Many definitions are in current use. Progress toward a successful definition of this property has been hampered by the refusal of those working in classical EHL to acknowledge the existence of accurate measurements of the piezoviscous effect that have existed for nearly a century. The Hamrock and Dowson pressure-viscosity coefficient at high temperature requires knowledge of the piezoviscous response at pressures which exceed the inlet pressure and may exceed the Hertz pressure. The definition of pressure-viscosity coefficient and the assumed equation of state must limit the use of the classical formulas, including Hamrock and Dowson, to liquids with high Newtonian limit and to low temperature. Given that this problem has existed for at least fifty years without resolution, it is reasonable to conclude that there is no definition of pressure-viscosity coefficient that will quantify the piezoviscous response for an analytical calculation of EHL film thickness at temperatures above ambient.

The unresolved definition of the pressure-viscosity coefficient Scott Bair In the classical approach to elastohydrodynamic lubrication (EHL) a single parameter, the pressureviscosity coefficient, quantifies the isothermal pressure dependence of the viscosity for use in prediction of film thickness. Many definitions are in current use. Progress toward a successful definition of this property has been hampered by the refusal of those working in classical EHL to acknowledge the existence of accurate measurements of the piezoviscous effect that have existed for nearly a century. The Hamrock and Dowson pressure-viscosity coefficient at high temperature requires knowledge of the piezoviscous response at pressures which exceed the inlet pressure and may exceed the Hertz pressure. The definition of pressure-viscosity coefficient and the assumed equation of state must limit the use of the classical formulas, including Hamrock and Dowson, to liquids with high Newtonian limit and to low temperature. Given that this problem has existed for at least fifty years without resolution, it is reasonable to conclude that there is no definition of pressure-viscosity coefficient that will quantify the piezoviscous response for an analytical calculation of EHL film thickness at temperatures above ambient. www.nature.com/scientificreports/ This is not to say that exponential response does not occur. Real piezoviscous response is slower than exponential (at times even linear) at low pressure and high temperature, becoming faster than exponential at high pressure as the glass transition pressure is approached. Therefore, Eq. (1) can be applied near the inflection in the log viscosity versus pressure curve which may include ambient pressure if the temperature is low. The inflection in the log viscosity versus pressure curve occurs at a characteristic viscosity that is roughly independent of temperature 4 . Here, μ is the limiting low shear (Newtonian) viscosity in contrast to η, the generalized Newtonian viscosity which varies with shear rate.

Abbreviations
One impediment to EHL becoming a quantitative field is the lack of a consensus for the definition of a pressure-viscosity coefficient which sufficiently characterizes the strength of the piezoviscous response in the inlet zone to allow for an accurate prediction of film thickness. A significant hindrance to this goal is the use of unrealistic assumptions in the classical film thickness formulas. Two of the unrealistic assumptions are the Newtonian inlet zone and the Dowson and Higginson equation of state 5 . Many definitions of pressure-viscosity coefficient are in current use. A few are listed below.

Some definitions of alpha
The qualitative prediction of the EHL thickness profile is rightly regarded as one of the greatest achievements in tribology 6 and the Hamrock and Dowson formula 7 for central film thickness, h C , of a point contact is widely used to estimate the thickness of the oil film in machinery in spite of a lack of quantitative validation.
Notice that parameters for shear-thinning and compressibility, which have a measurable effect on film thickness, are absent from this formula.
In classical EHL a single property, the pressure-viscosity coefficient or "alpha value", quantifies the isothermal pressure dependence of the viscosity of the liquid lubricant for prediction of film thickness. Progress toward a successful definition of this property has been hampered by the refusal of those working in classical EHL to acknowledge the existence of the accurate measurements of the piezoviscous effect available for nearly a century 2 . See, for example, a tribology textbook 8 . One definition of the pressure-viscosity coefficient, α B (T), invokes the "Barus equation" (1), requiring exponential response for some interval in pressure near ambient.
This definition and another, the tangent pressure-viscosity coefficient, α T require that a pressure be specified.
The tangent pressure-viscosity coefficient, evaluated at the Hertz pressure of a concentrated contact, influences the slope of the logarithmic portion of a traction curve, a slope which has been incorrectly associated with Eyring theory. The "conventional" pressure-viscosity coefficient is 10 α 0 .
The Hamrock and Dowson formula (2) employs yet another definition of pressure-viscosity coefficient, the reciprocal asymptotic isoviscous pressure coefficient, α*, a sort of weighted average of local coefficients. This definition was suggested by Blok 11 considering the behavior of a mineral oil at room temperature. The isoviscous pressure is found from integration 11,12 .
The isoviscous pressure provides a (Weibull) transformation whereby a piezoviscous solution for the pressures in the Reynolds equation can be found from an isoviscous solution 12 . The reciprocal asymptotic value is found by setting the argument of the isoviscous pressure function to infinity.
It would seem that the integration to infinite pressure would make this definition of the pressure-viscosity coefficient untenable. However, the integration may quickly converge to a good approximation of the asymptotic isoviscous pressure for many oils within the inlet pressure range when the temperature is low.

A problem for the Hamrock and Dowson coefficient, a*
One pressing problem for EHL film thickness prediction is the film forming capability at high temperature in a roller bearing of a jet engine oil, meeting specification PRF 7808 13 . The limiting low shear viscosity, μ(T, p) has been measured in a falling cylinder viscometer of the type introduced by Bridgman 2 nearly a century ago. These  Fig. 1 along with two correlations for the slower than exponential response at low pressures. The McEwen equation 14 is used by many laboratories that do accurate measurements of viscosity at elevated pressures in viscometers. It is significantly better suited for the slower than exponential response than the third Roelands 15 model often used in classical EHL and can be written with either α 0 or α* as a parameter.
The parameters of a fit to the data for pressures up to 400 MPa are listed in Table 2 along with those of the isothermal form of the Modulus equation 9 .
This equation is identical to the second Roelands model 15 and has a limiting viscosity of μ 0 exp(1/B). Although there is no real asymptotic limit to viscosity, the agreement with the data in Fig. 1 is quite satisfactory. The average absolute relative deviation, AARD, is listed in Table 2 for the two correlations and the fit to the data is excellent for both. www.nature.com/scientificreports/ Surprisingly, for the Modulus equation, α* is identically zero This is due to the limiting value for viscosity which causes p isovis (∞) to be unbounded. The McEwen equation produces values of α* that are in reasonable agreement with those obtained from the algorithm recommended by Bair 16 .
Equation (10) results in α* = 14.6 and 7.6 GPa −1 for 23 and 165 °C, respectively. The great difference in values of α* for the McEwen and the Modulus correlations, which give equivalent fits to the data, clearly show that there is insufficient information contained in the viscosities up to pressure of 400 MPa to evaluate α*. Therefore, the evaluation of the reciprocal asymptotic isoviscous pressure coefficient, α*, requires an understanding of how the viscosity depends upon pressure for pressures very much greater than the inlet zone pressure of a contact. For example, the algorithm (10) assumes exponential response for pressures greater than the highest measured pressure, p N , and the McEwen Eq. (8) assumes that the slower than exponential response simply continues, whereas the Modulus Eq. (9) reaches a limiting value of viscosity.

A problem for glass/steel EHD rigs
The contact between a glass disc and a steel ball is the working component of an EHL rig, an instrument that has become the standard for studies of EHL film thickness. In some ways, this is unfortunate because the Hertz pressure cannot exceed about 0.8 GPa and is typically no more than 0.5 GPa in most commercially available test rigs because of the strength of glass. The pressure in the inlet zone must be much less. The inlet zone where the film thickness is established ends at a pressure of 50-100 MPa according to Dench et al. 17 . A bit higher inlet pressure limit of 150 MPa is given by Spikes 18 . This means that, in many cases, the viscosity in the inlet zone may not reach sufficiently high pressure to contain the information needed to evaluate the Hamrock and Dowson pressure-viscosity coefficient, α*.
Performing the integration in Eq. (7) for the McEwen Eq. (8) yields an analytical expression for the isoviscous pressure.
The reciprocal of the isoviscous pressure is plotted in Fig. 2 for the jet engine oil meeting specification PRF 7808 at two temperatures, 23 and 165 °C. The maximum pressure of the inlet, assumed to be 100 MPa, is shown as the gray vertical line. It can be seen that, for low temperature, the reciprocal of the isoviscous pressure, 19.8 GPa −1 , is 35% greater than the asymptotic value, α* = 14.6 GPa −1 , at the pressure of 100 MPa. For the higher temperature, at 100 MPa the reciprocal of the isoviscous pressure is 16.2 GPa −1 whereas the asymptotic value is 6.8 GPa −1 , less than half. Clearly this definition is not appropriate for a jet engine oil at temperatures above ambient. It is useful, however, for a naphthenic mineral oil at room temperature, as shown in Fig. 3 where the reciprocal of the isoviscous pressure reaches to within 5% of α* at 100 MPa. However, for a helicopter oil, meeting specification PRF-85734, at 100 °C the reciprocal of the isoviscous pressure at 100 MPa is 16.6 GPa −1 while the asymptotic value is 9.4 GPa −1 .
The Hamrock and Dowson coefficient, α*, cannot be the operative pressure-viscosity coefficient in the standard glass/steel EHL rig. There is not enough information available regarding the piezoviscous effect within the inlet pressure range. The liquid cannot be responding to pressures greater than those existing in the film.

Discussion
It has repeatedly been shown that the pressure-viscosity coefficient cannot be accurately derived from EHL film thickness measurements [19][20][21][22] ; however, the method continues in use today, promoted primarily by Imperial College. The errors may be attributed to non-Newtonian inlet behavior, ignoring the temperature dependence of compressibility, operation of the EHL rig outside of the fully piezoviscous-elastic regime, and inaccurate film measurements. Another problem for the method has been shown above to be the lack of a suitable definition of www.nature.com/scientificreports/ the pressure-viscosity coefficient that is reported from the method. In particular, the film pressures in a glass/ steel contact are insufficient to reach the asymptotic value of the isoviscous pressure. The classical Hamrock and Dowson film thickness formula, which is used in this method, requires that p isovis (p) ≈ p isovis (∞) within the inlet zone and, therefore, the inlet pressure would have to reach to well over 500 MPa at high temperature. The definition of pressure-viscosity coefficient and the assumed equation of state must limit the use of the classical formulas, including Hamrock and Dowson, to liquids with high Newtonian limit and to low temperature. The severity of the situation can be illustrated in Fig. 4 by comparing the pressure-viscosity response of the jet oil at 165 °C to that of hypothetical oils with the same value of the Hamrock and Dowson pressure-viscosity coefficient, α*, but with different values of the McEwen exponent, q. For the jet oil, q = 2.3. Values of q for typical lubricating oils are in the range of 1.2-20. For traction fluids q may exceed 100. In Fig. 3, the hypothetical oils have q = 1.5 and 10 with the same α* = 6.8 GPa −1 as the jet oil. If the inlet pressure is less than 300 MPa, the jet oil will clearly have a smaller film thickness than the hypothetical q = 1.5 oil given the same compressibility and Newtonian response although they have equal α*. The conventional pressure viscosity coefficients are α 0 = 20, 12.1, and 7.5 GPa −1 for the hypothetical q = 1.5 oil, the jet oil, and the hypothetical q = 10 oil, respectively. However, it has already been shown that the conventional pressure viscosity coefficient does not correctly rank the  www.nature.com/scientificreports/ film forming capability 23 . Two oils with the same α 0 and μ 0 can generate very different film thicknesses in both experiments and simulations 23 . Fifteen years ago, this author suggested an alternative definition of pressure-viscosity coefficient, which in simulations appeared to offer a suitable replacement for α* 23 and provides a set of curves more compact at low pressures than does α* in Fig. 4 . However, this definition also invoked the integral of Eq. (7) requiring viscosity information at pressures which may exceed the contact pressure.

Conclusion
The Hamrock and Dowson pressure-viscosity coefficient at high temperature requires knowledge of the piezoviscous response at pressures which exceed the inlet pressure and may exceed the Hertz pressure. The definition of pressure-viscosity coefficient and the assumed equation of state must limit the use of the classical formulas, including Hamrock and Dowson, to liquids with high Newtonian limit and to low temperature. Given that this problem has existed for at least fifty years without resolution, it is reasonable to conclude that there is no suitable definition of pressure-viscosity coefficient that will quantify the piezoviscous response for an analytical calculation of EHL film thickness for temperatures above ambient. Consent to participate. The author consents to participate.

Data availability
Will make all data available on request.