Next Article in Journal
Study of Heat Flow at Substrate during Sputtering of Copper–Titanium Sandwich Target
Previous Article in Journal
Inline Hot Rolling of Al-5%Mg Strip Cast Using an Unequal Diameter Twin-Roll Caster
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Thermoelastic Contact of Gas-Lubricated Rough Sealing Faces

College of Mechanical Engineering, Zhejiang University of Technology, Hangzhou 310032, China
*
Author to whom correspondence should be addressed.
Materials 2024, 17(14), 3600; https://doi.org/10.3390/ma17143600
Submission received: 19 June 2024 / Revised: 17 July 2024 / Accepted: 19 July 2024 / Published: 21 July 2024

Abstract

:
Friction and wear are the main failure sources of face seals. When the surfaces of sealing rings exhibit greater roughness, the level of friction might increase and lead to sealing failure. Therefore, 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 novelty and advantage of this study is that it takes the effect of surface roughness into consideration during thermoelastic analysis of gas-lubricated seals. The film pressure, temperature, contact force and deformation of a gas spiral groove-faced seal are numerically determined. The influence of surface roughness on the contact distribution, deformation and temperature of the end-face seal at different speeds and pressures is analyzed. 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 analysis shows that the roughness contact mainly happens on the inner side of the rings due to convergent distortion of the seal faces, which easily causes partial wear of the seal faces. Moreover, it can also be found that the spiral grooves on the sealing surface can produce 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. The greater surface roughness results in a larger temperature rise under low-rotational-speed and lower-seal-pressure conditions, which further increases the risk of severe wear or even failure of the seal faces.

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]:
R a = 1 L 0 L z ( x ) d x
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]
σ = 1 L 0 L z x 2 d x
For a rough surface with Gaussian distribution, the relationship between the Ra value and σ is [18]:
R a = 2 π σ = 0.798 σ
σ′ is the standard deviation of the combined roughness of two surfaces of σ1 and σ2 [18].
σ = σ 1 2 + σ 2 2
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].
1 E = 1 2 1 v 1 2 E 1 + 1 v 2 2 E 2
The influence of temperature, T, on the elastic modulus, E, can be expressed as follows [18]:
E = E 0 e α T
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]:
p t = 1 + α T - T 0 2 16 2 15 π μ β σ 2 σ β E F 2.5 λ
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]:
F 2.5 ( λ ) = 1 2 π λ ( s λ ) 2.5 e s 2 2 d s
The actual contact area is as follows [1]:
A λ = π 2 μ β σ 2 A ˜ F 2 λ
where A ˜ is the nominal contact area and the calculation formula of F2(λ) is as follows [3]:
F 2 ( λ ) = 1 2 π λ ( s λ ) 2 e s 2 2 d s
According to the furrow effect [21], the relationship between the cutting force Fe and the actual contact area is calculated as follows.
F e = A σ s
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:
f = A ( d ) σ s p ˜ ( d ) = π 2 μ β σ 2 A ˜ F 2 d σ σ s 16 2 15 π μ β σ 2 E σ β A ˜ F 2.5 d σ = 15 π F 2 λ σ s 16 2 E σ β F 2.5 λ
Further, according to the frictional heat formula, the heat generated by solid contact is obtained:
w fric = f p t ω r

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]:
r θ Q θ ρ h 3 η p r θ + r Q r ρ h 3 r η p r = 6 ω Q s ρ h θ
where Qθ is the circumferential pressure flow factor, Qr is radial pressure flow factor, and Qs is shear flow factor.
Q θ = 1 C e γ ( h σ ) γ 1 Q θ = 1 + C ( h σ ) γ γ > 1
Q r h σ , γ = Q r h σ , 1 γ
Q s = A 1 h σ α 1 e α 2 h σ + α 3 h σ 2 h / σ 5 Q s = A 2 e 0.25 h σ h / σ > 5
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:
h = h 0 + h wave cos ( N w θ ) Non groove area h 0 + h g + h wave cos ( N w θ ) Groove area
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]:
q θ T r θ + q r T r = η ω 2 r 2 h ρ c v h 3 12 η ρ c v Q θ p r θ 2 + Q r p r 2 + f p t ω r ρ c v + k g , s 1 ρ c v T s 1 T + k g , s 2 ρ c v T s 2 T + R u ρ T h t
where
q θ = Q θ h 3 12 η p r θ + Q s ω r h 2 q r = Q r h 3 12 η p r
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]:
2 T s 1 r 2 θ 2 + r r r T s 1 r + 2 T s 1 z 2 = 0
Since there is motion in the direction of θ, the rotor ring’s heat conduction equation is obtained as shown in [22]
2 T s 2 r 2 θ 2 + r r r T s 2 r + 2 T s 2 z 2 = ω ρ s 2 c s 2 k c 2 T s 2 θ
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].
T = i d E m c v
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]:
p = R u ρ i d E m c v
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]:
p ( r , θ = π / N ) = p ( r , θ = π / N )
Equation (2) can be solved for the temperature distribution of the gas film. The dynamic temperature boundary condition is applied as follows [22]:
T ( r = r o , θ ) = T inlet   if   q r ( r = r o , θ ) < 0
T ( r = r i , θ ) = T i   if   q r ( r = r i , θ ) > 0
where
q r = h 3 12 η p r
q r ( r = r o , θ ) < 0 and q r ( r = r i , θ ) > 0 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]:
T ( r , θ = π / N ) = T ( r , θ = π / N )
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]:
p exit = M 1 M exit 2 + ( 1 + γ 1 ) M 1 2 2 + ( 1 + γ 1 ) M exit 2 p 1
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].
γ 1 = ln p 1 p exit ln ρ 1 ρ exit
At the inlet of the seal, pressure and temperature losses can be expressed as follows [22]:
p inlet = p o 1 + ( γ 2 1 ) M inlet 2 2 C L 2 γ 2 γ 2 1
T inlet = T o 1 + ( γ 2 1 ) M inlet 2 2 C L 2
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].
γ 2 = ln p inlet p 2 ln ρ inlet ρ 2
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]:
k c 1 ( T s n ) s = k gs 1 ( T s 1 T 0 )
k c 2 ( T s n ) s = k gs 2 ( T s 2 T 0 )
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:
F open = 0 2 π 0.5 d i 0.5 d o p r d r d θ
The closing force can be calculated as follows:
F close = π r o 2 r b 2 p o + π r b 2 r i 2 p i
The solid contact force can be calculated as follows:
F cont = 0 2 π r i r o p t r d r d θ
Average gas temperature, Tav, is defined as follows [22]:
T av = T r d θ d r r d θ d r
The convergence criterion is defined as follows [22]:
δ X = X j X j / 2 X j
where X = p , T , T s , T av 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.

4. Conclusions

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 following conclusions can be drawn.
(a)
Surface roughness often leads to significant contact pressure, which mainly happens on the inner side of the rings due to convergent distortion of seal faces, which can easily cause partial wear of the seal faces. Under low-speed conditions smaller than 500 rpm, a contact pressure of 0.25 MPa may be observed.
(b)
Greater surface roughness results in a greater temperature rise under low rotational speed and lower seal-pressure conditions, which further increases the risk of severe wear or even failure of the seal faces. In the low-speed cases of 50 rpm, the roughness contact may make the film temperature increase from 293.3 K to 650 K.
(c)
The proposed model can not only apply for graphite and steel, but also for other sealing pairs through changing the material parameters. Moreover, the simulated temperature of this model can extend to a larger scale.

Author Contributions

Methodology, S.B. and J.Y.; software, S.B.; formal analysis, Y.C.; investigation, S.B., J.Y. and Y.C.; resources, S.B.; writing—original draft preparation, S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (52375212).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data is contained with the article.

Conflicts of Interest

The authors declare no conflict on interest.

Nomenclature

AThe actual contact area [mm2].
A1, A2Calculation parameters.
CCalculation parameter.
cs2Specific heat capacity [J·kg−1·K−1].
cvThe constant specific heat capacity [J·mol−1·K−1].
CLThe loss coefficient.
EPlane-stress modulus [GPa].
E, E0 , E1 , E2Elasticity modulus [GPa].
EmMolecular energy per degree of freedom of gas [J·mol−1].
fThe friction coefficient.
FeThe cutting force [N].
hFilm thickness [m].
h0Minimum film clearance [m].
hdefromDeflection of film [m].
hgDepth of texture [m].
hwaveThe wave amplitude [m].
idFreedom of gas molecules.
kc1 , kc2Thermal conductivity of dynamic and static rings, [W·m−1·K−1].
kgs1 , kgs2Heat transfer coefficient [W·m−2·K−1].
LMeasured length [um].
M1 , MinletThe Mach number.
MfricDimensionless friction torque.
nRotation speed [rpm].
NGroove numbers.
NwWaviness numbers.
pLocal pressure [Pa].
poSeal pressure [Pa].
p1 , p2Gas pressure [Pa].
pcontContact pressure [Pa].
piInside pressure [Pa].
qLeakage rate [m3·s−1].
QDimensionless leakage rate.
QrRadial pressure flow factor.
QsShear flow factor.
QθCircumferential pressure flow factor.
rRadial coordinate [m].
rbBalanced radius [m].
rgSpiral groove radius [m].
riInside radius [m].
roOutside radius [m].
RDimensionless radius.
RaContour arithmetic mean [m].
RuThe universal ideal gas constant [J·mol−1·K−1].
TLubrication gas film temperature [K].
T0Ambient temperature [K].
T1Insider temperature [K].
T2Outside temperature [K].
TavAverage gas temperature [K].
TmaxMaximum film temperature [K].
TsSeal torus temperature [K].
wOpening force [N].
v, v1 , v2Poisson’s ratio.
α, α1 , α2 , α3Calculation parameters.
βAsperity radius, [m].
γ1 , γ2The expression of parameter.
ηVolume viscosity of gas [Pa·s].
θRotational direction [rad].
ωCircumferential velocity [rad·s−1].
λThe ratio of film thickness to roughness.
µ·β·σTopography parameters.
ρLocal density [mol·m−3].
ρs2Density of rotor ring [kg·m−3].
σThe contour root-mean-square deviation [m].
σ´ Standard deviation of combined roughness of two Surfaces of σ1 and σ2 [m].
σsThe yield limit of the material [MPa].

References

  1. Carbone, G.; Lorenz, B.; Persson, B.N.J. Contact mechanics and rubber friction for randomly rough surfaces with anisotropic statistical properties. Eur. Phys. J. E Soft Matter 2009, 29, 275. [Google Scholar] [CrossRef] [PubMed]
  2. Martin, H.; Müser Dapp, W.B.; Bugnicourt, R. Meeting the Contact-Mechanics Challenge. Tribol. Lett. 2017, 66, 118. [Google Scholar]
  3. Greenwood, J.A.; Tripp, J.H. The contact of two nominally flat rough surfaces. Proc. Inst. Mech. Eng. 1970, 185, 625–633. [Google Scholar] [CrossRef]
  4. Patir, N.; Cheng, H.S. An average flow model for determining effects of three-dimensional roughness on partial hydrodynamic lubrication. J. Lubr. Technol. 1978, 100, 12–17. [Google Scholar] [CrossRef]
  5. Patir, N.; Cheng, H.S. Application of average flow model to lubrication between rough sliding surfaces. J. Lubr. Technol. 1979, 101, 220–229. [Google Scholar] [CrossRef]
  6. Ruan, B.; Salant, R.F.; Green, I. A mixed lubrication model of liquid/gas mechanical face seals. Tribol. Trans. 1997, 40, 647–657. [Google Scholar] [CrossRef]
  7. Lebeck, A.O. Theory of Thermoelastic instability of rotating rings in sliding contact with wear. J. Lubr. Technol. Trans. ASME 1976, 98, 277. [Google Scholar] [CrossRef]
  8. Lebeck, A.O.; Tealej, L.; Piercer, E. Hydrodynamic Lubrication and Wear in Wavy Contacting Face Seals. J. Tribol. 1978, 100, 81. [Google Scholar] [CrossRef]
  9. Lebeck, A.O. Hydrodynamic lubrication in wavy contacting face seal—A two dimensional dodel. J. Tribol. 1981, 103, 578. [Google Scholar]
  10. Lebeck, A.O. A Test apparatus for measuring the effects of waviness in mechanical face seals. ASLE Trans. 1981, 24, 371–378. [Google Scholar] [CrossRef]
  11. Li, C.H. Thermal deformation in a mechanical face seal. Tribol. Trans. 1976, 19, 146–152. [Google Scholar] [CrossRef]
  12. Doust, T.G.; Parmar, A. An experimental and theoretical study of pressure and thermal distortions in a mechanical seal. Tribol. Trans. 1986, 29, 151–159. [Google Scholar] [CrossRef]
  13. Blasiak, S. An analytical approach to heat transfer and thermal distortions in non-contacting face seals. Int. J. Heat Mass Tran. 2015, 81, 90–102. [Google Scholar] [CrossRef]
  14. Thomas, S.; Brunetière, N.; Tournerie, B. Thermoelastohydrodynamic behavior of mechanical gas face seals operating at high pressure. J. Tribol. 2007, 129, 841–850. [Google Scholar] [CrossRef]
  15. Xie, J.; Ma, C.H.; Bai, S.X. Thermo-distortion characteristics of spiral groove gas face seal at high temperature. Numer. Heat Transf. 2020, 77, 242–256. [Google Scholar] [CrossRef]
  16. Zhu, D.L.; Yang, J.; Bai, S.X. Thermoelastohydrodynamic characteristics of low-temperature helium gas T-groove face seals. Materials 2021, 14, 2873. [Google Scholar] [CrossRef] [PubMed]
  17. Zhu, D.L.; Bai, S.X. Thermoelastohydrodynamic characteristics of supercritical CO2 spiral groove face seals. Ind. Lubr. Tribol. 2021, 73, 153–162. [Google Scholar] [CrossRef]
  18. Wang, Q.J.; Chung, Y.W. Encyclopedia of Tribology; Elsevier: Amsterdam, The Netherlands, 2013; ISBN 9780387928975. [Google Scholar]
  19. Etsion, I.; Front, I. A model for static sealing performance of end face seals. Tribol. Trans. 1994, 37, 111–119. [Google Scholar] [CrossRef]
  20. Hutchings, I.; Shipway, P. Tribology: Friction and Wear of Engineering Materials; Butterworth-Heinemann: Oxford, UK, 2017. [Google Scholar]
  21. Wen, S.Z.; Huang, P. Principles of Tribology, 4th ed.; Tsinghua University Press: Beijing, China, 2012. [Google Scholar]
  22. Bai, S.X.; Wen, S.Z. Gas Thermohydrodynamic Lubrication and Seals, 1st ed.; Elsevier Academic Press: Cambridge, MA, USA, 2019. [Google Scholar]
  23. Zuk, J.; Ludwig, L.P.; Johnson, R.L. Compressible flow across shaft face seals. In Proceedings of the Fifth International Conference on Fluid Sealing, BHRA, Coventry, UK, 30 March–2 April 1971; p. FICFS-H6. [Google Scholar]
  24. Zuk, J. Analysis of face deformation effects on gas film seal performance. ASLE Trans. 1973, 16, 267–275. [Google Scholar] [CrossRef]
  25. Di Schino, A. Manufacturing and Application of Stainless Steels. Metals 2020, 10, 327. [Google Scholar] [CrossRef]
Figure 1. Schematic of spiral face gas seal and thermal boundary conditions.
Figure 1. Schematic of spiral face gas seal and thermal boundary conditions.
Materials 17 03600 g001
Figure 2. GT contact model [3].
Figure 2. GT contact model [3].
Materials 17 03600 g002
Figure 3. Mechanical boundary conditions.
Figure 3. Mechanical boundary conditions.
Materials 17 03600 g003
Figure 4. Algorithm of the program.
Figure 4. Algorithm of the program.
Materials 17 03600 g004
Figure 5. Film thickness and contact pressure distributions under different rotational speed (Ra = 0.20 μm, po = 1.1 MPa).
Figure 5. Film thickness and contact pressure distributions under different rotational speed (Ra = 0.20 μm, po = 1.1 MPa).
Materials 17 03600 g005
Figure 6. Film pressure and temperature distributions under different rotational speed (Ra = 0.20 μm, po = 1.1 MPa).
Figure 6. Film pressure and temperature distributions under different rotational speed (Ra = 0.20 μm, po = 1.1 MPa).
Materials 17 03600 g006
Figure 7. Influence of rotational speed on frictional performance (po = 1.1 MPa).
Figure 7. Influence of rotational speed on frictional performance (po = 1.1 MPa).
Materials 17 03600 g007
Figure 8. Film thickness and contact pressure distributions at different surface roughness (ω = 5000 rpm, po = 1.1 MPa).
Figure 8. Film thickness and contact pressure distributions at different surface roughness (ω = 5000 rpm, po = 1.1 MPa).
Materials 17 03600 g008
Figure 9. Film pressure and temperature distributions at different surface roughness (ω = 5000 rpm, po = 1.1 MPa).
Figure 9. Film pressure and temperature distributions at different surface roughness (ω = 5000 rpm, po = 1.1 MPa).
Materials 17 03600 g009
Figure 10. Influence of surface roughness on frictional performance (po = 1.1 MPa).
Figure 10. Influence of surface roughness on frictional performance (po = 1.1 MPa).
Materials 17 03600 g010
Figure 11. Film thickness and contact pressure distributions at different seal pressure (Ra = 0.20 μm, ω = 5000 rpm).
Figure 11. Film thickness and contact pressure distributions at different seal pressure (Ra = 0.20 μm, ω = 5000 rpm).
Materials 17 03600 g011
Figure 12. Film pressure and temperature distributions at different seal pressure (Ra = 0.20 μm, ω = 5000 rpm).
Figure 12. Film pressure and temperature distributions at different seal pressure (Ra = 0.20 μm, ω = 5000 rpm).
Materials 17 03600 g012
Figure 13. Influence of seal pressure on frictional performance (Ra = 0.20 μm).
Figure 13. Influence of seal pressure on frictional performance (Ra = 0.20 μm).
Materials 17 03600 g013
Table 1. Parameter table of spiral groove seal.
Table 1. Parameter table of spiral groove seal.
Item and SymbolDimensions and Data
Inside radius, ri23 mm
Outside radius, ro31 mm
Spiral groove radius, rg26 mm
Groove numbers, N12
Groove depth, hg5 μm
Waviness height, hwave0.8 μm
Waviness numbers, Nw12
Table 2. Parameter table of asperity.
Table 2. Parameter table of asperity.
Item and SymbolDimensions and Data
Topography parameters, µβσ0.075
Standard deviation of asperity height, σ120.2 µm
Asperity radius, β{Tribology: friction and wear of engineering materials}250 µm
Young’s modulus, E125 GPa
Young’s modulus, E2204 GPa
Table 3. Characteristics of the seal ring materials.
Table 3. Characteristics of the seal ring materials.
CharacteristicsGraphiteStainless Steel [25]
Density/kg·m−319607800
Young’s modulus/G·Pa25204
Poisson’s coefficient0.20.27
Specific heat capacity/J·kg−1·K−1710460
Thermal conductivity/W·m−1·K−112916.5
Linear thermal expansion coefficient/10−6 °C44.3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bai, S.; Chen, Y.; Yang, J. Analysis of Thermoelastic Contact of Gas-Lubricated Rough Sealing Faces. Materials 2024, 17, 3600. https://doi.org/10.3390/ma17143600

AMA Style

Bai S, Chen Y, Yang J. Analysis of Thermoelastic Contact of Gas-Lubricated Rough Sealing Faces. Materials. 2024; 17(14):3600. https://doi.org/10.3390/ma17143600

Chicago/Turabian Style

Bai, Shaoxian, Yangyang Chen, and Jing Yang. 2024. "Analysis of Thermoelastic Contact of Gas-Lubricated Rough Sealing Faces" Materials 17, no. 14: 3600. https://doi.org/10.3390/ma17143600

APA Style

Bai, S., Chen, Y., & Yang, J. (2024). Analysis of Thermoelastic Contact of Gas-Lubricated Rough Sealing Faces. Materials, 17(14), 3600. https://doi.org/10.3390/ma17143600

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop