Analytical modelling of mechanical face seals: post #4

Analytical model for full film lubrication

In this post, a first model for full film lubrication is proposed. It will allow to determine the seal faces temperature, the thickness of the film between the faces, the friction torque, etc.

If the distance between the faces is assumed to be equal to the mean film thickness, the dissipated power by viscous friction is:

P=\mu \omega^2 S_f \frac{\left( r_o+r_i\right)^2}{4 h_m}

This power is transferred to the seal rings by conduction. According to post# 1, we have:

E_t \Delta T=\mu \omega^2 S_f \frac{\left( r_o+r_i\right)^2}{4 h_m}

According to post# 3, the film thickness is given by:

h_m=\frac{\beta \Delta R}{4 \left(B_t-0.5 \right)}

The coning angle \beta is due to the faces deformations as explained in post #2:

\beta=N \Delta T+\beta_e

Finally, the fluid viscosity  \mu is a function of the local temperature. By using an exponential law, we have:

\mu \left(\Delta T\right)=\mu_0 \exp \left(-\alpha \Delta T\right)

where \mu_o is the viscosity at the reference temperature and \alpha the thermoviscosity coefficient.

A non linear equation of the seal faces temperature is obtained:

E_t \Delta T=\mu_0 \exp{\left( -\alpha \Delta T\right)} \omega^2 S_f \frac{ \left( r_o+r_i \right)^2 \left( B_t-0.5 \right) }{ \left( N \Delta T+\beta_e \right) \Delta R }

It is convenient to use a dimensionless form of this equation. For that a temperature scale based on the thermoviscosity coefficient is used:

\bar{T}=\Delta T \alpha

Then the previous non linear equation temperature becomes:

\bar{T}=\exp{\left(-\bar{T}\right)}\frac{Se}{\bar{T} + Co}

In this equation Se is the sealing number and is a kind of thermal loading number of the seal face:

Se=\mu_o \pi \omega^2\frac{\left( r_o+r_i \right)^3\left( B_t-0.5 \right) }{E_t N }

The second dimensionless number, Co, is the coning number. It represents an initial coning of the faces due to mechanical loading (fluid pressure) or wear of the faces:

Co=\frac{\beta_e \alpha}{N}

According to governing equation, the seal faces temperature rises with the sealing number, Se. However, the magnitude of this temperature is reduced by the decrease in viscosity with temperature (exponential term) and the increase in film thickness with temperature (denominator). A positive initial coning number, Co , will increase the film thickness and thus reduce temperature whereas a negative value will produce the invert.

The governing equation can be easily solved by using Excel for example. When the temperature of the faces is known, it is then possible to compute the film thickness, the leakage flow, the friction torque, etc. The next figure presents the dimensionless temperature rise as a function of the sealing number for different values of the coning number.

In the case of negative coning number, the temperature can not be lower than -Co to avoid negative coning and film thickness. This situation is not realistic and will lead to faces contact. To analyze

this effect, the minimum film thickness can be expressed in a dimensionless form:

\bar{h}_{min}=Co+\bar{T}

 

The film film thickness is presented as a function of the sealing number; Se, for different values of the coning number Co in the next figure. As expected the film thickness is an increasing function of Se. For positive coning number, the minimum film thickness can not decrease below a threshold avoiding faces contact. For negative or null values of Co, the film thickness can reach zero at low sealing number values. In this case face contact will be experienced. According to my experience, the dimensionless film thickness corresponding to the appearance of contact (three times the r.m.s roughness height), is approximately equal to 0.5. Above this value, the present model is correct. Below this threshold mixed lubrication is expected to occur. This boundary between full film and mixed lubrication regimes has been added on the temperature map (dashed curve):

0.5=Co+\bar{T}

It is now necessary to consider contact between the surfaces of the seal.

References

Brunetière, N. & Apostolescu, A. A Simple Approach to the ThermoElastoHydroDynamic Behavior of Mechanical Face Seals Tribology Transactions, 2009, 52, 243-255

 

 

Leave a Reply