1. Introduction
Contact between rough surfaces is the direct cause of friction and the main cause of wear failure of face seals [
1,
2]. With the continuous application of gas face seals in aero-power machinery, the problem of thermoelastic contact has become an increasingly restrictive factor for the long life and reliable design of seals in the face of complex and variable conditions such as frequent start–stop, high temperature and high speed.
At present, the mixed-lubrication model based on Greenwood–Tripp (GT) model [
3] and Patir–Cheng (PC) model [
4,
5] is widely used in the thermal effect analysis of the mixed lubrication of the seal face, and the two rough surfaces of the friction pair are equivalent to a smooth rigid surface and a rough elastic surface. The GT model is close to the actual engineering situation and lays a foundation for research on contact problems in contemporary Tribology. In 1978, Pair and Cheng [
4,
5] proposed an average flow model (PC model for short) considering the contact conditions of micro-surface roughness peaks, and a calculation method for the influence of surface roughness on the flow of lubricating fluid was given by introducing the flow factor into the Reynolds equation. In 1997, Ruan [
6] analyzed the mixed lubrication characteristics of mechanical end seals based on the mixed lubrication model established by combining the PC model and the GT model.
With the application of groove technology in face sealing, the superposition of surface roughness and geometric groove on the face of the seal has further affected the bearing capacity of the sealing fluid. In a previous study, Lebeck carried out systematic theoretical and experimental research on the circumferential surface waviness in liquid lubricated contact mechanical seals, focusing on their thermoelastic instability [
7], hydrodynamic lubrication [
8], two-dimensional model [
9] and waviness effect [
10], showing that the circumferential waviness of the carbon ring surface after operation can be generated with an amplitude between 0.5 and 0.8 μm, resulting in an obvious hydrodynamic pressure effect. For the aqueous medium under pressure of 2.87 MPa and 93 °C, when the wave number is 3 and the dimensionless wave amplitude increases from 1 to 5, the proportion of liquid film bearing capacity can increase to 95%. When the wave number is nine, the hydrodynamic pressure effect is enhanced, and the proportion of liquid film bearing capacity can be increased to 100%, which can significantly reduce the contact wear.
Surface contact friction heating forms a temperature gradient in the radial direction of the seal ring, and the uneven distribution of thermal deformation on the seal face directly affects the contact state of the seal face, which changes from the theoretical surface contact to linear contact, resulting in poor fitting of the end face, an increased leakage rate and radial contact bias wear. In 1976, Li’s research pointed out that the thermal deformation of the liquid contact mechanical seal would cause the deformation of the seal end face in the amplitude of 2 μm, resulting in the convergence of the clearance and resulting in the contact wear of the seal ring at the inner diameter [
11]. In 1986, Doust and Parmar were the first researchers to experimentally measure the thermal deformation of the non-contact wedge seal, which was in good agreement with the theoretical results calculated by the boundary integral equation method [
12]. In 2015, Blasiak’s research on heat transfer and thermal deformation in a non-contact liquid end seal showed that when the rotating speed was 800 rad/s, the temperature of the sealing fluid lubrication film could rise by 16 K from the outer diameter to the inner diameter, resulting in the convergence and deformation of 0.7 μm seal clearance [
13].
For gas seals, the thermal deformation is also prominent (as observed in radial tapered face seals studied by Thomas in 2007 [
14] and grooved faced seals by Xie in 2020 [
15] and Zhu in 2021 [
16,
17]). Xie’s work [
15] shows that with the increase in ambient temperature, the thermal deformation reduces the opening force slightly, while the thermal viscosity effect significantly reduces the leakage rate. When the ambient temperature is increased from 300 K to 700 K, the opening force and leakage rate are reduced by 3.7% and 76%, respectively. The temperature difference between inside and outside diameters has a greater influence on the thermal deformation of the sealing surface. With the decrease in opening force and the increase in the leakage rate, the sealing surface will have a sharp divergence deformation. From 0 to 75 K, with the increase in temperature difference, the opening force decreases by 12% and the leakage rate increases by 224%. Theoretically, the thermal deformation problem not only leads to the complexity of the contact distribution, but also the temperature change leads to a change in the contact force. However, the effect of rough contact on seal properties such as temperature rise is not clearly quantified, especially considering coupling factors such as deformation of seal face.
In this paper, based on the elastic contact hypothesis of rough and wavy surfaces and the influence of temperature on the elastic modulus of materials, a thermoelastic contact lubrication model of a gas-lubricated end seal is established. The gas film pressure, temperature distribution, contact force and deformation of the spiral groove in the end face seal are numerically solved. The influence of surface roughness on the contact distribution, deformation and temperature of the seal end face at different speeds and pressures is analyzed.
2. Model Development
Figure 1 illustrates a typical gas spiral face seal, which consists of a smooth ring and a spiral-grooved one. In the analysis, a pressure–periodic boundary condition is applied, so a periodic computing region is defined by the dashed lines, as illustrated in the figure.
2.1. Thermoelastic Contact of Rough Surfaces
Here, we describe the method used to create the one-dimensional topography parameters, contour arithmetic mean deviation or center line mean
Ra, shown in
Figure 2. It is the arithmetic average of the height of each point on the contour within the measured length range [
18]:
where
z(
x) is the contour height of each point and
L is the measured length.
The contour root-mean-square deviation or root-mean-square value
σ is defined as [
18]
For a rough surface with Gaussian distribution, the relationship between the
Ra value and
σ is [
18]:
σ′ is the standard deviation of the combined roughness of two surfaces of
σ1 and
σ2 [
18].
E′ is the plane-stress modulus for the material, calculated from the elastic modulus
E1,
E2 and Poisson ratio
υ1,
υ2 of two rings [
18].
The influence of temperature,
T, on the elastic modulus,
E, can be expressed as follows [
18]:
Here, for the austenitic stainless steel, E0 = 31.673 and α = 3.9929 × 10−4.
When the mechanical seal is running, the heat has two main sources: the shear action of the fluid and the contact action of the solid. The frictional heat generated by surface contact is introduced into the energy equation of the gas film as the gas film heat source term.
According to the GT contact model, the nominal pressure caused by rough peak contact is obtained [
3]:
where
µβσ is the roughness characterization parameter, between 0.04~0.08 [
19,
20];
λ is the ratio of film thickness to roughness,
λ =
h/
σ.
Assuming that the height of the surface roughness peak conforms to the Gaussian distribution, the calculation formula of
F2.5(
λ) is as follows [
3]:
The actual contact area is as follows [
1]:
where
is the nominal contact area and the calculation formula of
F2(
λ) is as follows [
3]:
According to the furrow effect [
21], the relationship between the cutting force
Fe and the actual contact area is calculated as follows.
where
σs is the yield limit of the material,
σs = 206.8 MPa in the flowing analysis. Here, the cutting force is often considered as the friction force.
So, the friction coefficient
f can be calculated as follows:
Further, according to the frictional heat formula, the heat generated by solid contact is obtained:
2.2. Governing Equations
Since the magnitude of the film thickness is close to the magnitude of the roughness peak, the influence of the surface roughness on the flow of the medium must be considered. The average Reynolds equation with a contact coefficient and flow direction coefficient is introduced [
22]:
where
Qθ is the circumferential pressure flow factor,
Qr is radial pressure flow factor, and
Qs is shear flow factor.
where
γ represents the directionality of asperity.
A1,
A2,
C,
α1,
α2 and
α3 are the calculation parameters [
4,
5].
The clearance between two rings can be given as follows:
where
h0 is the minimum film clearance and
hg is the groove depth;
hwave is the wave amplitude and
Nw is the number of waves.
By adding the friction heat to the energy equation of the gas, we can obtain the following equation [
22]:
where
where
cv is the constant specific heat capacity of gas;
kgs1 and
kgs2 is the convection heat transfer coefficient between gas and ring;
Ts1 and
Ts2 is the ring’s surface temperature.
For the stator ring, the heat conduction equation is given as follows [
22]:
Since there is motion in the direction of
θ, the rotor ring’s heat conduction equation is obtained as shown in [
22]
where
kc2 is thermal conductivity of the rotor ring;
ρs2 is the density of the rotor ring;
cs2 is the specific heat capacity of rotor ring.
Here, assuming the gas molecules are rigid spheres, according to the principle of equal energy portions, the energy per degree of motion freedom is equal to
Em. Further, we assume that the gas temperature only represents the macroscopic inner energy of the gas molecular, so the following equation can be obtained [
22].
where
id is the freedom number of gas motions and
cv is the specific heat at a constant volume. Here,
id = 5 for the oxygen gas.
Further, it is assumed that gas pressure is determined only by both translation energy of gas molecules and gas density. Hence, the pressure component induced by the gas can be expressed as follows [
22]:
where
Ru is the universal ideal gas constant and the value is 8.31434 J/(mol·K).
2.3. Boundary Conditions of the Fluid Film
Equation (1) can be solved for the density and pressure distribution of the gas film. In the analysis, pressure–periodic boundary condition is applied as follows [
22]:
Equation (2) can be solved for the temperature distribution of the gas film. The dynamic temperature boundary condition is applied as follows [
22]:
where
and mean that the gas flows from ambient into the lubricating region at the boundaries.
In addition, temperature–periodic boundary condition is expressed as follows [
22]:
For a compressible fluid, the exit pressure is not equal to the surrounding inner pressure when the flow is choked [
23,
24]. The exit condition is expressed as follows [
22]:
where
Mexit is the Mach number at the seal exit with gas pressure
pexit and density
ρexit, and
M1 is the Mach number at any point near the exit with gas pressure
p1 and density
ρ1. Then, the expression of parameter
γ1 can also be obtained [
22].
At the inlet of the seal, pressure and temperature losses can be expressed as follows [
22]:
where the loss coefficient
CL is equal to 0.65 [
16].
Minlet is the Mach number at the seal inlet with gas pressure
pinlet and density
ρinlet. Then, the expression of parameter
γ2 can be obtained [
22].
where
p2 and
ρ2 are gas pressure and density, respectively, at any point near the inlet.
2.4. Boundary Conditions Relating to the Solids
Here, convective, adiabatic and imposed heat boundary conditions are used on surfaces as shown in
Figure 1, and the heat exchange condition between solid and ambient is as follows [
22]:
Figure 3 presents the geometry of the solids and the mechanical boundary conditions. The fluid film between the faces is delimited by inner and outer radii (
ri and
ro). On the rear face of the stator (the rotor), the area bounded by the inner radius (the balance radius) and the outer radius,
ri and
ro, respectively (
rb and
ro, respectively), can be in contact with the ring support. Meanwhile, the O-ring location on the stator,
horing, is also of concern, as it determines the influence of sealed pressure acting on the outside cylindrical surface of the stator. Non-penetration conditions with a rigid surface are imposed on rear surfaces of the stator and the rotor, with the ring supports being assumed to be rigid.
2.5. Numerical Method
The finite difference method is utilized to obtain the film pressure, film temperature and ring temperature. The finite element method is used for the coupling calculation of the face elastic and thermal distortions.
In numerical analysis, the opening force is expressed as follows:
The closing force can be calculated as follows:
The solid contact force can be calculated as follows:
Average gas temperature,
Tav, is defined as follows [
22]:
The convergence criterion is defined as follows [
22]:
where
and
j is the iterative number.
The mesh density of gas film is 80 × 80, and that of rings is 80 × 80 × 30. The value of error limit,
ε, for the convergence criterion is 10
−5. The film pressure, film temperature, ring temperature, face distortions and seal clearance are successively calculated into four overlapping loops as shown in
Figure 4. The new equilibrium position of the seal clearance is then established, and the entire iterative process is repeated until the convergence criterion on the equilibrium position is satisfied.
In the following numerical analysis, the dimensions of the seal, the parameters of asperity and the sealing materials are shown in
Table 1,
Table 2 and
Table 3 respectively. The proposed model can not only apply for graphite and steel, but also for other sealing pairs if the material parameters are changed.
3. Result and Discussion
For the gas face seal, the rough contact directly dependents on clearance between two sealing faces, which mainly affected by operating parameters and surface roughness. In the following analysis, we put points on the influence of rotational speed, seal pressure and surface roughness on the frictional performance.
3.1. Lubrication Status Transition with Rotational Speed
Generally, lubrication status transits from mixed lubrication to full-film lubrication with the increase in rotational speed. Meanwhile, the film thickness increases, which leads to a sharp decrease in contact pressure. In theory, the face deformation affects the surface contact distribution and changes with the sealing speed and film thickness. In this section, we calculate the influence of rough surface contact on lubrication state after coupling face deformation.
Figure 5 presents the film thickness and contact pressure distributions under different rotational speeds at
Ra = 0.20 μm and
po = 1.1 MPa. As can be seen, the film thickness increases as the rotational speed increases from 1 rpm to 2000 rpm, while the contact pressure sharply decreases from 0.25 kPa to 0. The more important thing is that the roughness contact mainly happens on the inner side of the rings, which means there will be a convergent distortion of seal faces.
Further, the pressure and temperature distributions of gas film between the seal surfaces can be obtained at different rotational speeds, as shown in
Figure 6. Clearly, under low-speed conditions smaller than 500 rpm, there is no significant temperature rise, although there is a contact pressure of 0.25 MPa. It can also be found that the spiral grooves on a sealing surface can produce an obvious hydrodynamic pressure effect due to the function of shear speed when the speed increases to 2000 rpm, while the film temperature increases from 293.3 K to about 306 K.
Figure 7 shows the influence of rotational speed on frictional performance. According to
Figure 7a, it can be found that the minimum film thickness demonstrates a tendency to increase with the growth of rotational speed from 0 to 10,000 rpm. Meanwhile, the contact pressure presents a law of first increasing and then decreasing. When the rotational speed increases to a speed exceeding 100 rpm, the value of contact pressure sharply decreases to zero.
As expected in
Figure 7b, the frictional torque presents obvious Stribeck curve with increase of rotational speed, especially in the case of
Ra = 0.4 μm. In another words, the lubrication status transfers from mix-lubrication to full film lubrication when the rotational speed reaches 100 rpm.
Another important conclusion can be found in
Figure 7b, which indicates that the greater surface roughness results in a larger temperature rise under low-rotational-speed conditions. For the case of
Ra = 0.4 μm, the maximum film temperature reaches about 650 K at a rotational speed of 50 rpm. So, the larger surface roughness may result in serious wear at the start-up stage.
3.2. Effect of Surface Roughness
It is well known that the surface roughness may lead to a significant temperature rise and risk of surface wear. In particular, the amount of surface roughness is the key factor that directly affects the temperature rise. But it is still not very clear how the magnitude of roughness contact affects the temperature rise. In this part of the paper, the effect of surface roughness on film pressure and temperature is further analyzed.
Figure 8 shows the film thickness and contact pressure distributions at different surface roughnesses at
ω = 5000 rpm and
po = 1.1 MPa. Clearly, the contact pressure sharply increases from 0 to 0.4 MPa with the increase in surface roughness from 0.05 μm to 0.60 μm. As a result, there is a significant convergent distortion of about 1 μm in the radius direction for the film thickness, as shown in the figure. The reason is that the increase in surface roughness leads to greater frictional heat, which makes the temperature increase significantly. As shown in
Figure 9, when the surface roughness increases from 0.05 μm to 0.60 μm, the maximum film temperature increases from 293.3 K to about 368 K.
Figure 10 show the influence of surface roughness on frictional performance. It can be found that, for the spiral groove face seal, the lubricating condition will deteriorate rapidly with increases in the surface roughness after 0.4 μm. Both the contact pressure and frictional torque increase significantly, which leads to a rapid increase in film temperature, exceeding 500 K. This means that excessive roughness increases the risk of wear failure of the seal face.
3.3. Effect of Seal Pressure
Seal pressure is another important factor affecting surface contact because of the face distortion. In theory, the greater the seal pressure, the greater the deformation. Due to the uneven deformation of the seal end face, the contact distribution will be very complicated. When designing high-pressure gas face seals, the coupling effect of end deformation and surface contact becomes more important.
Figure 11 shows the film thickness and contact pressure distributions at different seal pressures at
Ra = 0.20 μm,
ω = 5000 rpm. Clearly, the contact pressure sharply decreases from 0.14 kPa to 0 kPa as the seal pressure increases from 0.10 MPa to 2.3 MPa due to the increase in film thickness. The more important result is that, due to smaller face distortion at low seal pressure, the contact pressure presents radial uniform distribution. With the increase in seal pressure, the contact occurs mainly in the inner diameter due to the convergent distortion.
Additionally, since there is no obvious contact at higher seal pressure, the film temperature presents no obvious increase, as shown in
Figure 12.
Figure 13 show the influence of seal pressure on frictional performance. It can be found that, for the spiral groove face seal at lower seal pressure, there is an obvious temperature rise due to smaller film thickness and higher contact pressure. With the increase in seal pressure and film thickness, the contact pressure and frictional torque drastically decrease.