Next Article in Journal
Hybrid Impulsive Feedback Control for Drive–Response Synchronization of Fractional-Order Multi-Link Memristive Neural Networks with Multi-Delays
Previous Article in Journal
Fractional Factor Model for Data Transmission in Real-Time Monitoring Network
Previous Article in Special Issue
Thermal Conductivity Model of Porous Media Embedded with a Damaged Tree-like Branching Network Considering the Influence of Roughness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Total Leakage Characteristics of Finger Seal Considering Fractal Wear and Fractal Porous Media Seepage Effects

Faculty of Mechanical and Electrical Engineering, Kunming University of Science and Technology, Kunming 650500, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(7), 494; https://doi.org/10.3390/fractalfract7070494
Submission received: 31 May 2023 / Revised: 15 June 2023 / Accepted: 21 June 2023 / Published: 22 June 2023
(This article belongs to the Special Issue Transport Phenomena in Porous Media and Fractal Geometry)

Abstract

:
As an advanced flexible dynamic sealing technology, the leakage characteristics of a finger seal (FS) is one of the key research areas in this technology field. Based on the fractal theory, this paper establishes a mathematical model of the FS main leakage rate considering the fractal wear effect by taking into account the influence of the wear height on the basis of the eccentric annular gap flow equation. Based on the Hagen-Poiseuille law and the fractal geometry theory of porous media, a mathematical model of the FS side leakage rate considering the fractal porous media seepage effect is developed. Then, a mathematical model of the FS total leakage rate is established. The results show that the mathematical model of the FS total leakage rate is verified with the test results, the maximum error rate is less than 5%, and the mathematical model of the FS total leakage rate is feasible. With the gradual increase in working conditions and eccentricity, the FS main leakage rate gradually increases. In addition, the effects of the fractal dimension, fractal roughness parameters and porosity after loading on the FS main leakage rate are negligible. As the fractal dimension of tortuosity after loading gradually decreases, the fractal dimension of porosity after loading gradually increases, and the FS side leakage rate gradually increases. As the porosity after loading gradually increases, the FS side leakage rate gradually increases. Under different working conditions, different fractal characteristic parameters and different porosities after loading, the weight of the FS main leakage rate is much greater than that of the FS side leakage rate by more than 95%.

1. Introduction

As a green technology in the field of high-end equipment manufacturing, sealing technology has a decisive impact on the stable operation, energy consumption reduction, fatigue life and green cycle of high-end precision rotating machinery systems [1,2,3]. A finger seal (FS) is a kind of contact flexible dynamic sealing technology, which has excellent sealing performance and is widely used in the exhaust cavity of aircraft engine and the inlet and outlet of compressor [4]. The dynamic leakage rate of an FS is about 1/2 of that of standard carbon seal and 4/5 of that of standard brush seal. Its manufacturing cost is significantly lower than that of carbon seal and brush seals [5,6,7,8]. Since the FS can meet the strict requirements of high cycle parameters of aircraft engines, many researchers have carried out a lot of research on the sealing performance and mechanism of the FS.
As one of the critical indicators to measure the sealing performance, the flow leakage characteristic has been the research hotspot in the FS field. Wang et al. [9] found that the pressure fluid has a great influence on the sealing performance of an FS. Du et al. [10] found that when the pressure difference is more than 0.3 MPa, the leakage coefficient tends to be flat, and the pressure difference has no significant effect on the leakage coefficient in the interference fit. Bai et al. [11] first proposed a two-dimensional porous media computational fluid dynamics (CFD) model of an FS. Hu et al. [12] and Wang et al. [13] used a two-dimensional porous media CFD model to study the leakage characteristics of the FS. However, the two-dimensional porous media CFD model cannot accurately reflect the influence of the hysteresis gap on the leakage characteristics of an FS, so the prediction of the leakage performance of an FS may not be accurate. Chen et al. [14] first proposed a finite element analysis model that can properly reflect the hysteresis gap between the FS and the rotor. Su [15], Zhang et al. [16], Yin et al. [17], and Zhao et al. [18] used the finite element analysis model of the FS to study the influence of the hysteresis gap on the leakage performance of the FS. However, most of the researchers did not consider the effect of wear characteristics on the leakage performance of the FS.
In fact, the wear characteristics will change the assembly state of the FS/rotor, causing the seal structure to fail and the leakage rate to shift. Zhang et al. [19] found that the leakage rate is closely related to the amount of wear. Persson [20] found that the plastic deformation of metallic materials has a significant effect on the leakage rate of metallic seals. Du et al. [21] found that wear has a significant effect on the leakage characteristics in the early stage of the experiment; the effect of wear on the leakage characteristics is gradually reduced in the middle and late stages of the experiment. Therefore, it is crucial to construct the mathematical model of FS leakage characteristics considering the wear characteristics. Wang et al. [22] established a finite element analysis model of dynamic leakage of an FS considering the progressive wear effect based on the Archard model, and they found that the leakage rate considering the progressive wear effect was closer to the experimental results. Wang et al. [23] established a numerical analysis model for the transient leakage characteristics of an FS considering the wear effect based on the Archard model, and they numerically studied the leakage performance of the FS by using the two-dimensional porous media CFD model. The research showed that the leakage in the finger foot area was mainly in the early stage of wear, and the leakage in the finger beam area was mainly in the late stage of wear. However, the Archard model cannot accurately reflect the changing characteristics of the friction surface morphology in the wear process, and there is a certain error in the prediction of the wear rate, which may lead to the deviation of the prediction of the leakage rate.
Fractal theory can exclusively and quantitatively describe the contact behaviour of the contact surface, and it can effectively characterise the morphological changes of the contact surface and the leakage channel changes of the sealing interface during the operation of the sealing structure [24,25,26]. Maleki et al. [27] found that surface topography has a significant effect on the leakage rate of sealed parts. A study by Lv et al. [28] found that fractal feature parameters can effectively characterise the morphological features of the sealing surface and have a significant effect on leakage performance. Zhao et al. [29] established a fractal model for the mechanical seal leakage rate based on fractal theory, taking into account the characteristics of contact surface morphology. Li et al. [30] established a fractal model for mechanical seal end face leakage based on fractal theory and analysed the effect of changes in seal end face morphology on leakage performance from a microscopic point of view. Sun [31] used fractal feature parameters to describe in detail the morphological characteristics of the dynamic and static ring end faces of mechanical seals, established a fractal model for mechanical seal leakage, and investigated the relationship between leakage characteristics, fractal feature parameters, and working conditions. Feng et al. [32] established a fractal model for metal gasket leakage based on fractal theory, which accurately revealed the relationship between leakage characteristics and sealing surface morphology. Ni et al. [33] found that the fractal dimension has a significant effect on the contact surface throughout the working phase, and the fractal characteristic parameters have an effect on the contact surface only in the initial phase. The above research shows that fractal theory can effectively describe the morphological changes of contact surfaces and accurately reveal the relationship between the morphological changes of friction surfaces and leakage characteristics.
In recent years, most researchers have paid great attention to the fluid leakage caused by hysteresis and wear of FSs and rotors. In fact, the leakage channels of FSs are mainly divided into main leakage channels and side leakage channels. The main leakage channel is the hysteresis gap formed between the finger beam and the rotor after deformation, as well as the radial wear height caused by the wear of the FS. The radial gap is formed by the combination of the above two reasons; the side leakage channel is a microchannel formed by the interconnection of micropores on the rough contact surface of the adjacent finger seal annulus (FSA). Bai et al. [34,35] through experimental research and numerical analysis, found that the side leakage channel has a significant effect on the total leakage rate of the FS. Zhao et al. [36] considered adjacent FSAs as porous media structures and established a mathematical model for the side leakage channel rate of FSs based on Hagen-Poisson’s law and Darcy’s law. Through numerical research, they found that the side leakage channel dominates when the contact surface roughness is steep. The above research indicates that the side leakage channel can improve the total leakage rate of the FS, and it has a significant effect on the total leakage rate of the FS, which is of great significance for the composition of the total leakage channel of the FS. In a complete finger sealing system, several layers of FSAs are overlapped, staggered and stacked, tightly joined by riveting, and the adjacent FSAs are subjected to high axial pressure for a long time during the working process. As a result, the gap between adjacent FSAs is particularly narrow. Lv et al. [37] found through numerical research that rough surface characteristics have a significant effect on the air-tightness of sealing structures. Bao et al. [38] established a calculation method for the leakage rate of mechanical seals using the porous medium theory, and their research found that the surface roughness characteristics have a significant effect on the leakage rate. Huang et al. [39] established a prediction model for the leakage rate of metal gaskets based on the theory of fractal porous media, and they found that the geometric morphology characteristics of a rough surfaces are crucial factors affecting the leakage characteristics. The above research indicates that the morphological characteristics of rough surfaces have a significant impact on the side leakage channels and can be analysed from both microscopic and fractal perspectives.
Based on the fractal theory, this paper uses fractal characteristic parameters to characterise the wear height and establishes a mathematical model of the main leakage rate of the FS considering the fractal wear effect. Then, based on the fractal porous media transport theory, a mathematical model of the side leakage rate of the FS considering the fractal surface morphological characteristics is established. In addition, the mathematical model of the total leakage rate of the FS is established. The mathematical model of the total leakage rate not only reveals the influence of the fractal wear effect on the leakage rate from the macroscopic point of view, but also reveals the influence of the microscopic morphological features of the seal surface on the leakage rate from the microscopic point of view. The parameters in each of the mathematical models constructed above have specific physical meanings.

2. Basic Situation of Finger Seal

A standard type FS system, consisting of front and back baffles, a front and back fixed distance annulus, and a multiple contact FSA, is shown in Figure 1a. The contact FSA consists of a series of curved beams machined in the circumferential direction according to a certain rule, as shown in Figure 1b. We applied the following FSA principles: adjacent FSAs are staggered stacked, and have close contact; we always ensure that the gap between the former finger beam can be covered by the latter finger beam, effectively preventing the fluid through the gap between the finger beam axial leakage; the plates are always located on the high- and low-pressure air flow sides; the fixed distance annulus prevents the front and back plate and FSA from having direct contact to reduce the hysteresis effect; the FSA and rotor contact form an interference fit to achieve a critical sealing interface, effectively preventing fluid flow from high to low pressure. The FSA compensates for the shortcomings of the grate seal with high leakage, and the finger beam is similar to the brush wire bundle in the brush seal in the structure, which can adapt to the vortex effect formed by the rotor due to uneven material and the thermal effect. In addition, the rigidity and strength of the finger beam is greater than that of the brush wire bundle, overcoming the defect of the brush wire breakage.
There are many random and irregular asperities and gaps on the contact surfaces of the adjacent FSA. In this paper, the asperities are considered as skeletons and the gaps are considered as pores, so that the contact surfaces of the adjacent FSA can be reduced to microscopic pore structures. In this paper, the microscopic leakage channels on the contact surfaces of the adjacent FSA are assumed to be parallel capillaries of the same diameter and bent. The total flow space in the thin tube is equal to the volume of all the leakage channels in the sealing contact surface. The internal surface area of the thin tube is equal to the area of all the leakage channels in the rough contact surface. The asperity of the contact surface is circular in cross section in any direction perpendicular to the depth of roughness.
In this paper, the arc FS structure is used for the analysis. The main structural parameters are shown in Figure 2 and Table 1. Here, Dw is the outer circle radius; Df is the root circle radius; Dr is the inner circle radius; Dcc is the base circle radius; Rc is the arc radius; r is the rotor radius; Ia is the gap width; α is the clearance angle; gd is the downstream protection height; xg is the heel height. In this paper, the material used for the FSA is GH605, with a density of 9.13 × 10−9 t/mm3; a Young’s modulus of 2.31 × 105 MPa; and a Poisson’s ratio of 0.286. The rotor is GH4169 with a density of 8.24 × 10−9 t/mm3; a Young’s modulus of 2.04 × 105 MPa; and a Poisson’s ratio of 0.3. The back plate is 1Cr13 with a density of 7.75 × 10−9 t/mm3; a Young’s modulus of 2.17 × 105 MPa and a Poisson’s ratio of 0.3.

3. Analysis Method of Total Leakage Characteristics of Finger Seal

3.1. Mathematical Model for the Main Leakage Rate of Finger Seal

The FS main leakage is mainly divided into leakage caused by the radial hysteresis gap and leakage caused by the radial wear height, as shown in Figure 3. The FS and the rotor are usually assembled by interference assembly, and the rotor generates a vortex effect due to uneven material. Thus, it is difficult to maintain strict concentricity between the rotor and the FS, which often has a certain amount of eccentricity, and the FS and the rotor may form an eccentric annulus gap. Therefore, the qm mathematical model for predicting eccentric annulus gap leakage is [40].
q m = 2 π r h 3 12 η t s r Δ p 1 + 1.5 ε 2 + u 0 2 π r h
where r is the rotor radius, mm; h is the leakage gap, mm; Δp is the axial pressure difference, MPa; ε is the eccentricity; u0 is the flow rate of the seal medium; η is the dynamic viscosity of the seal gas; and tsr is the total thickness of the FSA, mm.
According to the research results of Zhang et al. [41], it is known that the circumferential flow velocity u0 at the rotor surface due to friction is
u 0 = π n r 60
where n is the rotor speed, r/min.
t s r = t m t b
where tm is the number of layers of the FSA; tb is the thickness of a single FSA, mm.
h = h g + h w
where hg is the radial hysteresis clearance, mm, obtained by finite element calculation; and hw is the mean radial wear height, mm.
The average radial wear height of the FS is mainly represented by the radial length variation of the finger foot, so the mathematical model for the wear height is [42]
h w = Δ V A a
where ΔV is the wear volume of the FS, mm3; and Aa is the nominal effective contact area between the FS and the rotor, mm2.
Fractal theory can effectively describe the complex behaviour of frictional wear and contribute to the quantitative analysis of tribological behaviour [25]. Therefore, in this paper, a fractal wear model is used to describe the wear volume of the FS. In the FS system, micromotion friction exists between adjacent FSAs; sliding friction exists between the rotor and the finger boots of each FSA; and the FSA operates in a very harsh environment with high-pressure intensity for a long time. As a result, the contact surfaces of the FSAs eventually tend to become completely elastically deformed when repeatedly subjected to the above loads [43]. In this paper, the research results of Lei et al. [44] are improved, and the improved mathematical model for the volume wear of the FS is shown in Equation (6).
Δ V = 1 + ψ μ 2 0 . 5 K e A r e Δ S
where ψ is the shear stress influence coefficient on the actual contact area [45]; μ is the friction coefficient; Ke is the elastic contact wear coefficient; Are is the elastic contact area; and ΔS is the sliding distance.
Δ S = v 0 Δ t  
where v0 is the relative linear sliding linear velocity between the FS and the rotor; and Δt is the relative sliding time.
v 0 = π n r 30
where n is the rotor speed; and r is the rotor radius.
A r e = λ D 2 D a l a l 0 . 5 D a c 1 0 . 5 D
where λ is the contact coefficient between the FS and the rotor [44]; D is the fractal dimension; al is the maximum contact area of the asperity; and ac is the critical contact area for plastic deformation of the asperity [46].
λ = 720 N t s r F r * 0 . 5 E π 3 0 . 5 r + D r 360 N f b α N n t s r r * 1
where Nfb is the total number of finger beams in a single FSA; F is the total pressure on the contact surface; r* is the comprehensive radius of curvature; E is the composite Young’s modulus; Dr fis the inner radius of the FSA; α is the clearance angle between the finger feet; and Nn is the number of finger beams in contact with a single FSA and rotor.
r * 1 = r i r d r
where ri is the interference amount.
E = 1 v 1 2 E 1 + 1 v 2 2 E 2 1
where E1 and E2 are the Young’s modulus of the two materials, respectively; v1 and v2 are Poisson’s ratio of the two materials, respectively.
a c = G 2 33 π 0 . 5 φ k μ 40 2 1 D
where G is the fractal roughness parameter that reflects the height of the asperity contour; φ is the characteristic parameter of the material, φ = ơy/E, where ơy is the yield strength of the softer material, and E is the composite Young’s modulus; and kμ is the correction factor for friction force.
k μ = 1 0 . 228 μ 0 μ 0 . 3 0 . 932 e 1 . 58 ( μ 0 . 3 ) 0 . 3 < μ 0 . 9
where μ is the friction coefficient.
A a = 360 α N f b π r 180 N f b N n t sr
Therefore, if Equations (6) and (15) are substituted into Equation (5), the wear height hw is
h w = 180 N f b K e A r e Δ S ( 1 + γ μ 2 ) 1 2 360 N f b α π r N n t s r
Therefore, by substituting Equation (16) into Equation (4), it can be seen that the leakage gap h is
h = h g + 180 N   f b K e A r e Δ S ( 1 + γ μ 2 ) 1 2 360 N   f b α π r N n t s r
In summary, the mathematical model Qm for the FS main leakage rate, taking into account the fractal wear effect, is as follows
Q m = π r h 3 Δ p 6 η t m t b 1 + 1 . 5 ε 2 + n h r 2 π 2 120

3.2. Mathematical Model for the Side Leakage Rate of Finger Seal

The volumetric flow of fluid through the microscopic channels of the adjacent FSA rough surface is referred to as the FS side leakage. The G-W model [47] first considered the contour height distribution of the contact surface as a random variable and assumed that the mean radius of curvature is closely related to the experimental instrument resolution. Sayles et al. [48] found that the contour height distribution of the contact surface has a non-stationary random characteristic and the standard deviation of its height distribution is the sampling length. Therefore, the statistical parameters of the surface topography are closely related to the instrument resolution and sampling length, and the characterisation of the contact surface is not unique and cannot fully reflect all the information of the surface roughness. Mandelbrot et al. [49] showed that the rough surface has self-similarity and self-simulation, the local and overall characteristics of the profile are self-similar, and this feature can be characterised by the fractal characteristic parameters, independent of the experimental instrument resolution and scale. This fractal characteristic parameter can provide full roughness information on the fractal surface for all the scale ranges.
The materials used for the FS are mostly cobalt-based superalloys. According to reference [42], cobalt-based superalloys have self-similarity and fractal characteristics, so the FS contact surface is a contact surface with fractal characteristics. The W-M function is an ideal fractal curve which is non-differentiable and has self-radiating fractal characteristics. It has uniqueness and certainty in characterising the height of the contact surface profile and its mathematical model is [50]
Z x 0 = G D 1 n = n m i n cos 2 π γ n x 0 γ 2 D n
Here, Z(x) is the contour curve of the rough surface; x0 is the coordinate of the contour; G is the fractal roughness parameter which reflects the magnitude of the rough contour curve; D is the fractal dimension which can describe the irregularity of the rough contour curve; γ is a constant, and for the rough surface obeying the normal distribution, γ = 1.5; γn is the spatial frequency of the contour curve; and nmin is the lowest cut-off frequency of the contour curve corresponding to the order number.
In this paper, based on Equation (19), a 2D rough surface scheme with profile height varying with different fractal dimensions (D = 1.4, 1.5, 1.6) and the same fractal roughness parameter (G = 1 × 10−12) is generated as shown in Figure 4, and a 2D rough surface scheme with profile height varying with different fractal roughness parameters (G = 1 × 10−11, 1 × 10−10, 1 × 10−9) and the same fractal dimension (D = 1.3) is generated as shown in Figure 5.
From Figure 4, it can be seen that as the fractal dimension gradually increases, the contour height of the rough surface gradually decreases, because the larger the fractal dimension, the smaller the height of the asperity on the rough surface, and therefore, the lower the contour height of the rough surface. This conclusion can be further extended to indicate that the larger the fractal dimension, the lower the profile height of the rough surface, the smoother the contact surface, the greater the number of asperities in contact in the contact surface, and the larger the contact area, which is consistent with the conclusion in the reference [30,51]. From Figure 5, it can be seen that the height of the profile of the rough surface gradually increases with the increase in the fractal roughness parameter, because the larger the fractal roughness parameter is, the higher the height of the asperity on the rough surface is, and therefore the height of the profile of the rough surface is larger. The conclusion can be further extended that the greater the fractal roughness parameter, the greater the contour height of the rough surface, the rougher the contact surface, the smaller the number of asperities in the contact surface that are in contact with each other, and the smaller the contact area. Therefore, the fractal dimension D and the characteristic roughness parameter G can affect the contour size of the rough surface and are some of the important parameters to characterise the contour features of the fractal surface.
If each asperity in the contact surface satisfies Equation (19), the profile curve before deformation can be described as [50]
z x 0 = G D 1 l 2 D cos π x 0 l l 2 < x 0 < l 2
where, l is the diameter of the bottom surface of the asperity profile.
For an asperity satisfying Equation (20), the height of the contact gap hb before loading can be expressed as [52]
h b = G D 1 l 2 D
The Roth model [53] provides the relationship between the contact gap height ha and the axial pressure difference Δp after loading as follows
h a = h a e Δ p R c
where, Rc is the sealing performance coefficient [37], which depends on the yield strength of the material and the characteristic parameters of the rough surface.
R c = 64 . 14 Δ p 0 . 3406 Δ p 90 MPa
Therefore, the deformation of the asperity after being loaded ω by
ω = h a h b
According to the research results of Yu et al. [54], the initial porosity of porous media before loading ϕ0 is
ϕ 0 = λ m i n λ m a x D E D f b
where, λmin is the minimum pore diameter [54]; According to reference [55], (λmin/λmax) < 10−2; DE is the Euclidean dimension, which takes the value of 2 in this paper; and Dfb is the pore fractal dimension before loading.
According to the research results of Ji et al. [52], the porosity of porous media after loading ϕ by
ϕ = h b ϕ 0 ω h a
Therefore, the pore fractal dimension Df0 after loading is
D f 0 = D E l n ϕ ln λ m i n / λ m a x
The FS side leakage channel is shown in Figure 6. From Figure 6, it can be seen that the FS side gap leakage channel is mainly composed of interconnected microscopic leakage channels on the contact surfaces of adjacent FSAs, the width of the leakage channel is related to the thickness of the FSA, and its length is the outer radius of the FS minus the inner radius.
In this paper, the microscopic pore structure on the contact surface of adjacent FSAs is considered as a fractal porous medium, and the capillary bundle model is used to describe the fractal porous structure. Extensive references [55,56,57,58] show that the pore area and pore size distribution of porous media satisfy the fractal scalar law, and the relationship between the number of pores N and the pore diameter λ0 is [54].
N L λ 0 = λ m a x λ 0 D f b
where L is the measured length; and λmax is the maximum pore diameter.
Therefore, from Equation (28), the total number of pores Nt is
N t L λ m i n = λ m a x λ m i n D f b
Differentiating Equation (29) gives the number of pores −dN in the interval λ to λ + .
d N = D f b λ m a x D f λ D f + 1 d λ
where −dN > 0.
Considering the fractal geometry of the microscopic pore structure of the contact surface, the leakage channel can be viewed as a bundle of curved capillaries with different cross-sectional areas. According to the Hagen-Poiseuille equation, the fluid flow rate q(λ) in a single curved capillary is [59]
q λ = π Δ p λ 4 128 L t λ η
where Lt(λ) is the actual length of the bent capillary before loading [59]; and η is the viscosity of the fluid.
L t λ = λ 1 D T b L 0 D T b
where L0 is the length of the straight line along the macroscopic pressure gradient direction, mm; DTb is the fractal dimension of tortuosity before loading; and the greater the DTb value, the more tortuous the capillary [55].
L 0 = D w D r
where Dw is the outer circle radius of the FS; and Dr is the inner circle radius of the FS.
D T b = 1 + ln τ a v e 0 ln L 0 / 2 λ a v e 0
where τave0 is the average tortuosity before loading [60]; and λave0 is the average pore radius before loading [55,61].
τ a v e 0 = 1 2 1 + 1 2 1 ϕ 0 + 1 ϕ 0 1 1 ϕ 0 1 2 + 1 4 1 1 ϕ 0
where ϕ0 is the initial porosity of the porous medium before loading.
λ a v e 0 = D f b λ m i n D f b 1
where Dfb is the fractal dimension of the pore before loading.
Therefore, the fractal dimension of the tortuosity DT after loading is
D T = 1 + ln τ a v e ln L 0 / 2 λ a v e
where τave is the average tortuosity after loading; and λave is the average pore radius after loading.
τ a v e = 1 2 1 + 1 2 1 ϕ + 1 ϕ 1 1 ϕ 1 2 + 1 4 1 1 ϕ
where ϕ is the initial porosity of the porous medium after loading.
λ a v e = D f 0 λ m i n D f 0 1
where Df0 is the fractal dimension of the pore after loading.
Integrating Equation (31), the FS side leakage rate Qb, taking into account the fractal porous media seepage effect before loading, is given by
Q b = λ m i n λ m a x q λ d N λ = π Δ p D f b λ m a x 3 + D T b 128 η L 0 D T b 3 + D T b D f b
According to the actual working conditions of the FS system, it is known that the contact surface of the adjacent FSA is subjected to high-strength pressure for a long time, so it is necessary to improve Equation (40) by replacing Dfb with Df0 and DT with DTb; then, the mathematical model of the FS side leakage rate Qs, taking into account the fractal porous media seepage effect after loading, is
Q s = λ m i n λ m a x q λ d N λ = π Δ p D f 0 λ m a x 3 + D T 128 η L 0 D T 3 + D T D f 0

3.3. Mathematical Model for Total Leakage Rate of Finger Seal

According to the above analysis, the FS total leakage rate Q
Q = Q m + Q s = π r h 3 Δ p 6 η t m t b 1 + 1 . 5 ε 2 + n h r 2 π 2 120 + π Δ p D f 0 λ m a x 3 + D T 128 η L 0 D T 3 + D T D f 0

4. Numerical Calculation Method for Dynamic Performance of Finger Seal

4.1. Finite Element Calculation Model for Dynamic Performance of Finger Seal

Most of the performance analyses of the FS belong to static performance analyses, and fewer FS dynamic performance analyses are considered, so most of the theoretical results have large errors with the actual test results [62,63,64]. In this paper, the numerical calculation method of the dynamic performance of the FS is used to obtain the value of the hysteresis gap in Equation (18) by the finite element method. The FS structure has the characteristics of cyclic symmetry and staggered superposition, so the finite element calculation model of the FS dynamic performance can be simplified into four parts. In the first part, there is one complete finger beam in the circumferential direction of the FSA and two half-finger beams, which is referred to as the high-pressure FS. In the second part, there are two complete finger beams in the circumferential direction of the FS, which is referred to as the low-pressure FS. In the third part includes part of the back plate. The fourth part includes part of the rotor. The finite element calculation model of the dynamic performance of the FS is shown in Figure 7.

4.2. Boundary Condition

According to the actual situation of the FS structure in operation, fixed constraint boundary conditions must be applied to the circumferential surfaces of the high-pressure FSA and low-pressure FSA, the circumferential surface and back of the back plate, and the axial direction of the rotor when performing finite element calculations. The fixed constraint boundary condition is usually understood as a displacement of 0, which means that the three components u, v and w along the x, y and z axes are 0.
q = u x , y , z v x , y , z w x , y , z = u , v , w x , y , z T = 0
where {q} is the displacement.
Then it is also necessary to apply the cyclic symmetry constraint to the high-pressure FSA, the low-pressure FSA cross-section and the rotor cross-section, and the mathematical model of the cyclic symmetry constraint boundary condition is
u A u B = cos k α n sin k α n sin k α n cos k α n u A u B
where ua and ub are the basic sector displacement and the replica sector displacement of the low-angle edge, respectively; u’a and u’b are the basic sector displacement and the replica sector displacement of the high-angle edge, respectively; k denotes the harmonic index and αn is the sector angle.
α n = 360 N s
where Ns is the number of sectors.

4.3. Setting of Load Step and Extraction of Calculation Results

The calculation process of the FS needs to be divided into three load steps: the first load step, which applies an axial pressure to the high-pressure FSA; the second load step, which applies a periodic displacement function in the radial direction of the rotor [65], which is used to simulate the vortex effect of the rotor due to unbalance effects; and the third load step, which resets the rotor.
The mathematical model of periodic displacement function y(t) is [65]
y t = Δ r sin n π t 30 + Δ r x sin n π t 30 ± r i
where Δr is the radial displacement excitation of the rotor; Δrx is the inclination of the rotor; “+” is the is the interference amount; and “−” is the amount of clearance.
At the end of the third load step, the gap between the deformed FSA and the rotor, the hysteresis gap, is extracted from the result file and substituted into Equation (18) to obtain the main leakage rate of the FS.

5. Results and Discussion

5.1. Verification of the Accuracy of Theoretical Models

It should be noted that the FS is a shaft seal component, as shown in Figure 1. Zhao et al. [36] first divided the FS total leakage (Q) into the FS main leakage (Qm) and the FS side leakage (Qs). In this paper, the FS main leakage (Qm) is divided into the leakage caused by the hysteresis gap and the leakage caused by the radial wear height, and in this paper, the FS side leakage (Qs) is considered as the flow rate of the adjacent FSA microscopic leakage channel after loading. Due to the limited testing equipment, only the FS total leakage (Q) can be measured at present [65,66,67]. Therefore, the test value of the main leakage volume of the FSA cannot be measured by the test alone and cannot be compared individually; similarly, the test value of the side leakage volume of the FSA cannot be measured by the test alone and cannot be compared individually. Furthermore, the mathematical model of the FS main leakage rate (Qm) and the FS total leakage rate (Q) established in this paper is only applicable to the fully elastic deformation stage of the asperity.
The test results in reference [67] were used to compare the analysis with the mathematical models Q (Equation (42)) for the FS total leakage, Qm (Equation (18)) for the FS main leakage, Qs (Equation (41)) for the FS side leakage and the mathematical model of the FS total leakage developed in reference [68] and established in this paper, as shown in Figure 8. The finite element calculation model used in the theoretical model validation was consistent with the structural parameters of the FS specimen for the FS bench test in the reference [67].
The calculation results show that the leakage gradually decreases with the increase in rotor speed, and the numerical calculation results keep the same trend and order of magnitude of change compared with the experimental results of the reference [67]. The maximum error between the numerical calculation results of the reference [68] and the experimental results of the reference [67] is 4.56%. The maximum error between the numerical calculation results of the FS total leakage rate mathematical model (Equation (42)) established in this paper and the experimental results is 2.195%. The maximum error between the numerical calculation results of the FS main leakage rate mathematical model (Equation (18)) and the experimental results is 2.196%. Meanwhile, for the FS side leakage rate mathematical model, the error rate between the numerical calculation results and the experimental results is very large. The numerical calculation results of the mathematical model of the FS total leakage rate established in this paper are closer to the experimental results than the numerical calculation results of the reference [68], because the eccentricity of the rotor, the wear height of the FS/rotor, and the FS side leakage are considered in this paper. The above three factors have an influence on the calculation of the FS total leakage rate, so considering the influence of these three factors in the calculation of the leakage rate will be more consistent with the actual operating condition of the FS. The comparison of the above numerical calculation results proves the reasonableness and correctness of the numerical calculation method of the FS considering the change of fractal surface shape.

5.2. Analysis of Factors Influencing the Main Leakage Rate of Finger Seals

5.2.1. Working Conditions

The influence law of the FS main leakage rate was analysed under different axial pressure differences, different rotor speeds and different radial displacement excitations, as shown in Figure 9. When n = 13,000 r/min, D = 1.3, G = 1 × 10−9, from Figure 9a, it can be seen that the FS main leakage rate gradually increases with the continuous increase in the axial pressure difference; when the axial pressure difference is constant, the FS main leakage rate also gradually increases with the increase of the radial displacement excitation. When Δp = 0.3 MPa, D = 1.3, and G = 1 × 10−9, it can be seen from Figure 9b that the FS main leakage rate gradually Increases with the gradual increase in the radial displacement excitation; when the radial displacement excitation is constant, the FS main leakage rate gradually decreases with the increase in the rotor speed. When Δr = 0.03 mm, D = 1.3, G = 1 × 10−9, it can be seen from Figure 9c that the FS main leakage rate gradually decreases with the increase in the rotor speed; when the rotor speed is constant, the FS main leakage rate gradually increases with the increase in the axial pressure difference. The variation laws in Figure 9a–c are consistent with the trend of the leakage rate under the actual working conditions of the FS.

5.2.2. Fractal Dimension

Under different fractal dimension D, different axial pressure differences Δp and different rotor speeds n, the influence law of the main leakage rate of FS is analysed, as shown in Figure 10. When Δr = 0.03 mm, n = 13,000 r/min, and G = 1 × 10−9, it can be seen from Figure 10a that the FS main leakage rate gradually increases with the increasing axial pressure difference; when the axial pressure difference is constant, the FS main leakage rate gradually increases with the increasing fractal dimension, because the larger the fractal dimension, the larger the contact area is, and according to Equation (6), the wear rate of the FS is larger, and then according to Equation (5), the wear height of the FS is higher, so the FS main leakage rate is greater than the leakage rate of the FS. When Δr = 0.03 mm, Δp = 0.3 MPa, and G = 1 × 10−9, it can be seen from Figure 10b that the FS main leakage rate decreases as the rotor speed continues to increase; at constant rotor speed, the FS main leakage rate increases as the fractal dimension increases, and the reason for this phenomenon is the same as in Figure 10a.

5.2.3. Fractal Roughness Parameter

The influence law of the FS main leakage rate was analysed for different fractal roughness parameters at different axial differential pressures and different rotor speeds, as shown in Figure 11. When Δr = 0.03 mm, n = 13,000 r/min, and D = 1.3, it can be seen from Figure 11a that the FS main leakage rate gradually increases with the increasing axial pressure difference; when the axial pressure difference is constant, the FS main leakage rate gradually decreases with the increasing fractal roughness parameter, because the larger the fractal roughness parameter, the smaller the contact area, and according to Equation (6), the smaller the wear rate of the FS, and then according to Equation (5), the smaller the wear height of the FS, so the FS main leakage rate is smaller. When Δr = 0.03 mm, Δp = 0.3 MPa, and D = 1.3, from Figure 11b, it can be seen that the FS main leakage rate gradually decreases as the rotor speed continues to increase; at constant rotor speed, the FS main leakage rate gradually becomes smaller as the fractal roughness parameter increases, and the reason for this phenomenon is consistent with Figure 11a.

5.2.4. Eccentricity

Under different values of eccentricity ε, axial pressure difference Δp and rotor speed n, the influence law of the FS main leakage rate is analysed, as shown in Figure 12. When Δr = 0.03 mm, n = 13,000 r/min, D = 1.3, and G = 1 × 10−9, the effect of axial pressure difference Δp and eccentricity ε on the main leakage rate is shown in Figure 12a. From Figure 12a, it can be observed that the main leakage rate gradually increases as the axial pressure difference continues to increase. When the axial pressure difference is constant, the main leakage rate gradually increases with the increase in eccentricity, because the larger the eccentricity, the larger the radial displacement excitation of the rotor, and the larger the leakage channel between the finger foot and the rotor, so the main leakage rate gradually increases. When Δr = 0.03 mm, Δp = 0.3 MPa, D = 1.3, and G = 1 × 10−9, the effect of rotor speed n and eccentricity ε on the main leakage rate is shown in Figure 12b. It can be seen from Figure 12b that as the rotor speed continues to increase, the main leakage rate gradually decreases; when the rotor speed is constant, the main leakage rate gradually increases with the increase in eccentricity. The reason for this phenomenon is consistent with Figure 12a.

5.3. Analysis of Factors Influencing the Side Leakage Rate of Finger Seals

5.3.1. Axial Pressure Difference and Fractal Characteristic Parameters

Under different values of axial pressure difference Δp, tortuosity fractal dimension DT and pore fractal dimension Df0, the influence law of the FS side leakage rate is analysed, as shown in Figure 13. When Df0 = 1.6, D = 1.3, and G = 1 × 10−9, the effects of axial pressure difference Δp and tortuosity fractal dimension DT on the side leakage rate are shown in Figure 13a. From Figure 13a, it can be observed that as the axial pressure difference continues to increase, the side leakage rate gradually increases; when the axial pressure difference is constant, the side leakage rate gradually decreases as the tortuosity fractal dimension continues to increase, because the larger the tortuosity fractal dimension, the higher the actual leakage channel bends, the longer the actual leakage channel, the more obvious the change trend of the side leakage rate, and the smaller the flow of the side leakage rate. When DT = 1.3, D = 1.3, G = 1 × 10−9, the axial pressure difference Δp and pore fractal dimension Df0 on the side leakage rate are shown in Figure 13b. From Figure 13b, it can be seen that the side leakage rate gradually increases as the axial pressure difference gradually increases; when the axial pressure difference is constant, the side leakage rate gradually increases with the increase in the pore fractal dimension, because the larger the pore fractal dimension, the larger the porosity, so the larger the side leakage channel, the larger the side leakage rate.

5.3.2. Porosity

When DT = 1.3, D = 1.3, and G = 1 × 10−9, the influence law of the FS side leakage rate was analysed under different axial pressure differences Δp and porosity ϕ, as shown in Figure 14. From Figure 14, it can be observed that as the axial pressure difference continues to increase, the side leakage rate gradually increases; when the axial pressure difference remains constant, as the porosity continues to increase, the side leakage rate gradually increases and decreases. This is because the larger the porosity, the larger the micro-leakage channel between the contact surfaces, resulting in a larger flow rate and a higher side leakage rate.

5.4. Analysis of Influencing Factors on Leakage Performance of Finger Seal

5.4.1. Working Conditions

Under different values of axial pressure difference, radial displacement excitation and rotor speed, the influence laws of the total leakage, main leakage and side leakage performance of the FS are analysed, as shown in Figure 15. When D = 1.3, G = 1 × 10−9, ϕ = 0.5, Δr = 0.03 mm, and n = 13,000 r/min, the effects of total leakage, main leakage, and side leakage under different axial pressure differences Δp are shown in Figure 15a. It can be seen from Figure 15a that with the continuous increase in axial pressure difference, the increase in the total leakage rate and main leakage rate is large, and the increase in the side leakage rate is small. The above phenomenon is consistent with the actual working conditions of the FS. When D = 1.3, G = 1 × 10−9, ϕ = 0.5, Δp = 0.3 MPa, and n = 13,000 r/min, the effects of total leakage, main leakage, and side leakage under different radial displacement excitation Δr are shown in Figure 15b. It can be seen from Figure 15b that with the continuous increase in the radial displacement excitation, the total leakage rate and the main leakage rate increase significantly, and the side leakage rate remains unchanged because the side leakage rate is not affected by the radial displacement excitation. When D = 1.3, G = 1 × 10−9, ϕ = 0.5, Δp = 0.3 MPa, and Δr = 0.03 mm, the effects of total leakage, main leakage, and side leakage under different values of rotor speed n are shown in Figure 15c. It can be seen from Figure 15c that with the continuous increase in the rotor speed, the total leakage rate and the main leakage rate gradually decrease, and the side leakage rate remains unchanged, because the side leakage rate is not affected by the rotor speed.
Corresponding to Figure 15, Table 2 shows the ratio of the main leakage rate and the side leakage rate of the FS under different values of axial pressure difference Δp, and it can be seen that the axial pressure difference Δp has a significant effect on the ratio of the main leakage rate. Table 3 shows the ratio of the main leakage rate and side leakage rate of the FS under different radial displacement excitation Δr. It can be seen that the radial displacement excitation Δp has a significant effect on the ratio of the main leakage rate ratio. Table 4 shows the ratio of the main leakage rate and the side leakage rate of the FS at different rotor speeds n. It can be seen that the rotor speed n has a significant effect on the ratio of the main leakage rate.

5.4.2. Fractal Feature Parameters

Under different fractal dimensions and fractal roughness parameters, the influence law of the total leakage, main leakage and side leakage performance of the FS is analysed, as shown in Figure 16. When G = 1 × 10−9, ϕ = 0.5, Δp = 0.3 MPa, n = 13,000 r/min, and Δr = 0.03 mm, the effects of total leakage, main leakage, and side leakage under different values of fractal dimension D are shown in Figure 16a. It can be seen from Figure 16a that the total leakage rate and the main leakage rate are increasing with the increasing of fractal dimension, and the reasons for this conclusion are consistent with those in Figure 10a. The side leakage rate has practically no shift, because the influence of the fractal dimension on the side leakage rate is extremely small and can be ignored. When D = 1.3, ϕ = 0.5, Δp = 0.3 MPa, n = 13,000 r/min, and Δr = 0.03 mm, the effects of total leakage, main leakage, and side leakage under different values of fractal roughness parameter G are shown in Figure 16b. From Figure 16b, it can be observed that as the feature scale coefficient continues to increase, the total leakage rate and the main leakage rate continue to decrease, and the reasons for this conclusion are consistent with those in Figure 11a. The side leakage rate remains practically unchanged because the influence of the fractal roughness parameter on the side leakage rate is extremely small and can be ignored.
Corresponding to Figure 16, Table 5 shows the ratio of the main leakage rate to the side leakage rate of the FS under different values of fractal dimension D. It can be seen that the fractal dimension has a significant effect on the ratio of the main leakage rate. Table 6 shows the ratio of the main leakage rate to the side leakage rate of the FS under different values of fractal roughness parameter G, and it can be found that the fractal roughness parameter has a significant effect on the ratio of the main leakage rate.

5.4.3. Porosity

When D = 1.3, G = 1 × 10−9, Δp = 0.3 MPa, n = 13,000 r/min, and Δr = 0.03 mm, the influence laws of the total leakage, main leakage, and side leakage performance of the FS under different values of porosity ϕ were analysed as shown in Figure 17. From Figure 17, it can be seen that as the porosity increases, the total leakage rate and the side leakage rate continue to increase. This is because the larger the porosity, the larger the leakage channel, and therefore the larger the leakage rate. Meanwhile, the main leakage rate remains virtually unchanged as it is not affected by porosity.
Corresponding to Figure 17, Table 7 shows the ratio of the main leakage rate to the side leakage rate of the FS under different porosities, and it can be seen that porosity has a significant effect on the ratio of the side leakage rate.

6. Conclusions

In this study, based on the fractal theory and the fractal geometry theory of porous media, a numerical analysis method is proposed for the finger seal (FS) total leakage characteristics considering the fractal wear and fractal porous media seepage effects, and a mathematical model for the FS main leakage rate considering the fractal wear effect and a mathematical model for the FS side leakage rate considering the fractal porous media seepage effect are proposed as well. The influence of working conditions, fractal characteristic parameters and other parameters on the total, main and side leakage rates of the FS is numerically investigated. The main conclusions are as follows:
  • The maximum error rate between the numerical calculation and the experimental value of the mathematical model of the FS total leakage rate is 2.2%, which is less than 5%, and the mathematical model of the FS total leakage rate is feasible. However, the model is only applicable to the study of the leakage rate during the fully elastic deformation phase of the asperity, and the use of this model assumes that the contact surface has fractal characteristics. The FS main leakage rate mathematical model is used in the same way as the FS total leakage rate mathematical model.
  • A mathematical model for the FS side leakage rate considering the seepage effect of fractal porous media takes into account the physical characteristics of isotropic porous media before and after loading, and the model can better reflect the FS side leakage characteristics, but the model does not take into account the physical characteristics of each anisotropic porous media and the non-Darcy seepage characteristics of porous media.
  • In the analysis of the FS total leakage performance under various parameters, the FS main leakage rate is above 95% and the FS side leakage rate is below 5%. Thus, the FS’s main leakage rate is much larger than the FS’s side leakage rate, so the density of the material used for the FS can be improved to reduce the porosity, or non-metallic materials (such as C/C composite materials, etc.) can be used to reduce the FS wear rate and thus improve the sealing performance.
  • In this study, the fractal characteristic parameters are used to characterise the contact surface morphology, which in turn reflects the leakage characteristics of the FS. It is recommended that in subsequent studies, the arithmetic mean roughness Ra or the effective root mean square surface roughness σ be used to characterise the contact surface topography and then Ra or σ be used to characterise the leakage performance.

Author Contributions

Conceptualisation, J.L. and M.L.; methodology, J.L. and M.L.; software, J.L. and Y.W.; validation, J.L., W.C. and Y.W.; formal analysis, J.L.; investigation, J.L.; resources, J.L.; data curation, W.C.; writing—original draft preparation, J.L.; writing—review and editing, W.C.; visualisation, J.L.; supervision, M.L.; project administration, M.L. 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 (grant number 51765024), and Analysis and Measurement Fund of Kunming University of Science and Technology (grant number 2021P20193103004).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Wang, Z.D.; Wang, Z.G. Development and Prospect of Brush Seal in Aero-engine. Lubr. Eng. 2005, 5, 203–205. [Google Scholar]
  2. Chupp, R.; Farshad, G.; Moore, G. Applying Abradable Seals to Industrial Gas Turbines. In Proceedings of the 38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Indianapolis, IN, USA, 7–10 July 2002. [Google Scholar]
  3. Chupp, R.E.; Hendricks, R.C.; Lattime, S.B.; Steinetz, B.M. Sealing in turbomachinery. J. Propuls. Power 2006, 22, 313–349. [Google Scholar] [CrossRef]
  4. Chen, G.D.; Su, H.; Zhang, Y.C. Analysis and Design of Finger Seal; Northwestern Polytechnical University Press: Xi’an, China, 2012; pp. 40–55. [Google Scholar]
  5. Zhang, Y.C.; Zhang, Y.T.; Wang, T.; Cui, Y.H.; Liu, K. Experimental Study on Performances of Carbon Seal and Finger Seal under High-speed and High-pressure Condition. In Proceedings of the IOP Conference Series: Materials Science and Engineering, Bristol, UK, 18–20 May 2018. [Google Scholar]
  6. Proctor, M.P.; Kumar, A.; Delgado, I.R. High-Speed, High-Temperature Finger Seal Test Results. J. Propuls. Power 2004, 20, 312–318. [Google Scholar] [CrossRef] [Green Version]
  7. Zhou, K.; Li, N.; Guo, W.; Pan, J.; Tan, J. Experimental Investigation on Static and Dynamic Characteristics of Finger Seals. Lubr. Eng. 2018, 43, 132–136. [Google Scholar]
  8. Proctor, M.P.; Delgado, I.R. Preliminary Test Results of Non-Contacting Finger Seal on Herringbone-Grooved Rotor. In Proceedings of the 44th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Hartford, CT, USA, 21–23 July 2008. [Google Scholar]
  9. Wang, L.N.; Chen, G.D.; Su, H.; Lu, F. Transient Performance Analysis of Finger Seal Considering Compressed Fluid in the Leakage Gap Effect. J. Aerosp. Power 2015, 30, 2004–2010. [Google Scholar]
  10. Du, C.H.; Ji, H.H.; Hu, Y.P.; Luo, J.; Ma, D.; Tang, L.P.; Liao, K. Experiment on Axis Locus of Rotor and Leakage Characteristics of Finger Seal. J. Aerosp. Power 2016, 31, 2575–2584. [Google Scholar]
  11. Bai, H.L. Experimental and Numerical Study on Leakage Characteristics of Finger Seal. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2007; pp. 32–47. (In Chinese). [Google Scholar]
  12. Hu, T.X.; Zhou, K.; Li, N.; Pan, J. Numerical Calculation and Experimental Study on Leakage Characteristic of Finger Seal. J. Propuls. Technol. 2020, 41, 1089–1096. [Google Scholar]
  13. Wang, Q.; Hu, Y.P.; Ji, H.H. An Anisotropic Porous Media Model for Leakage Analysis of Finger Seal. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2020, 234, 280–292. [Google Scholar] [CrossRef]
  14. Chen, G.D.; Xu, H.; Yu, L.; Xiao, Y.X.; Zhou, L.J. On Selecting Proper Shape of Finger Seal. J. Northwestern Polytech. Univ. 2002, 20, 218–221. [Google Scholar]
  15. Su, H. Structural Design Performance Analysis and Tests of Finger Seal. Ph.D. Thesis, Northwestern Polytechnical University, Xi’an, China, 2007; pp. 35–43. (In Chinese). [Google Scholar]
  16. Zhang, Y.C.; Yin, M.H.; Zeng, Q.R.; Wang, T.; Wang, R. Theoretical and Experimental Investigation of Variable Stiffness Finger Seal. Tribol. Trans. 2020, 63, 634–646. [Google Scholar] [CrossRef]
  17. Yin, M.H.; Zhang, Y.C.; Zhou, R.M.; Zhai, Z.Y.; Wang, J.L.; Cui, Y.H.; Li, D.S. Friction Mechanism and Application of PTFE Coating in Finger Seal. Tribol. Trans. 2021, 65, 260–269. [Google Scholar] [CrossRef]
  18. Zhao, H.L.; Su, H.; Chen, G.D. Analysis of Total Leakage of Finger Seal with Side Leakage Flow. Tribol. Int. 2020, 150, 106371. [Google Scholar] [CrossRef]
  19. Zhang, Y.C.; Liu, K.; Zhou, L.J.; Hu, H.T. Analysis of Leakage Characteristics of Finger Seal Based on System Responses. J. Aerosp. Power 2013, 28, 205–210. [Google Scholar]
  20. Persson, B.N.J. Leakage of Metallic Seals: Role of Plastic Deformations. Tribol. Lett. 2016, 63, 42. [Google Scholar] [CrossRef]
  21. Du, C.H.; Ji, H.H.; Hu, Y.P.; Luo, J.; Ma, D.; Tang, L.P.; Liao, K. Experiment of Wearing Characteristics and its Effect on Leakage of Finger Seal. J. Aerosp. Power 2017, 32, 53–59. [Google Scholar]
  22. Wang, L.N.; Sun, H.C.; Sun, W.; Zhou, X.Q.; Cui, Y.H.; Sun, L.C.; Meng, D.H. Fluid-solid Coupling Performance Analysis of Finger Seal under Considering Progressive Wear Effect. Lubr. Eng. 2021, 46, 47–53. [Google Scholar]
  23. Wang, Q.; Ji, H.H.; Hu, Y.P.; Du, C.H. Transient Leakage Characteristics Analysis of Finger Seal Considering Wear Effects. J. Propuls. Technol. 2020, 41, 2815–2826. [Google Scholar]
  24. Sun, J.J.; Gu, B.Q.; Wie, L. Leakage Model of Contacting Mechanical Seal Based on Fractal Geometry Theory. J. Chem. Ind. Eng. 2006, 57, 1626–1631. [Google Scholar]
  25. Ge, S.R.; Zhu, H. Fractal in Tribology; China Machine Press: Beijing, China, 2005; pp. 126–154. [Google Scholar]
  26. Persson, B.N.J.; Yang, C. Theory of the leak-rate of seals. J. Phys. Condens. Matter 2008, 20, 315011. [Google Scholar] [CrossRef]
  27. Maleki, I.; Wolski, M.; Woloszynski, T.; Podsiadlo, P.; Stachowiak, G. A Comparison of Multiscale Surface Curvature Characterization Methods for Tribological Surfaces. Tribol. Online 2019, 14, 8–17. [Google Scholar] [CrossRef] [Green Version]
  28. Lv, B.; Han, K.; Wang, Y.Z.; Li, X.L. Analysis and Experimental Verification of the Sealing Performance of PEM Fuel Cell Based on Fractal Theory. Fractal Fract. 2023, 7, 401. [Google Scholar] [CrossRef]
  29. Zhao, Y.X.; Ding, X.X.; Wang, S.P. Prediction of Leakage Rate and Film Thickness of Mechanical Seal Based on Fractal Contact Theory. Lubr. Eng. 2022, 47, 156–163. [Google Scholar]
  30. Li, X.P.; Yang, Z.M.; Wang, L.L.; Yang, Y.X. Leakage Model of Contacting Mechanical Seal Based on Fractal Theory. J. Northeast. Univ. (Nat. Sci.) 2019, 40, 526–530. [Google Scholar]
  31. Sun, J.J. Research on Mechanical Seal Leakage Prediction Theory and Its Application. Ph.D. Thesis, Nanjing Tech University, Nanjing, China, 2006. [Google Scholar]
  32. Feng, X.; Gu, B.Q. Research on Leakage Model of Metallic Gasket Seal. Lubr. Eng. 2006, 8, 78–80. [Google Scholar]
  33. Ni, X.Y.; Sun, J.J.; Ma, C.B.; Zhang, Y.Y. Wear Model of a Mechanical Seal Based on Piecewise Fractal Theory. Fractal Fract. 2023, 7, 251. [Google Scholar] [CrossRef]
  34. Bai, H.L.; Wang, W.; Zhang, Z.S.; Hu, G.Y. Analysis of Leakage Flow through Finger Seal Based on Porous Medium. J. Aerosp. Power 2016, 31, 1303–1308. [Google Scholar]
  35. Bai, H.L.; Ji, H.H.; Ji, G.J.; Cao, G.Z. Experimental Investigation on Leakage Characteristics of Finger Seal. J. Aerosp. Power 2009, 24, 532–536. [Google Scholar]
  36. Zhao, H.L.; Chen, G.D.; Su, H. Analysis of Total Leakage Performance of Finger Seal Considering the Rough Seepage Effect. J. Mech. Eng. 2020, 56, 152–161. [Google Scholar]
  37. Lv, X.K.; Yang, W.J.; Xu, G.L.; Huang, X.M. The Influence of Characteristic of Rough Surface on Gas Sealing Performance in Seal Structure. J. Mech. Eng. 2015, 51, 110–115. [Google Scholar] [CrossRef]
  38. Bao, C.Y.; Meng, X.K.; Li, J.Y.; Peng, X.D. The Leakage Performance Analysis of Mechanical Seals under Hydrostatic Pressures Based on Porous Media Model. Lubr. Eng. 2015, 40, 57–63. [Google Scholar]
  39. Huang, X.M.; Yao, B.; Xu, G.L.; Lu, X.K. Research on Leakage of Metallic Gasket Based on Fractal Porous Seepage. J. Huazhong Univ. Sci. Technol. (Nat. Sci. Ed.) 2016, 44, 1–5. [Google Scholar]
  40. Ning, C.X. Hydraulic and Pneumatic Technology; Chemical Industry Press: Beijing, China, 2017; pp. 24–27. [Google Scholar]
  41. Zhang, Y.C.; Liu, K.; Hu, H.T.; Song, F. Quasi-Dynamic Performances Analysis of Finger Seal Based on Finite Element Simulation. J. Propuls. Technol. 2016, 37, 2352–2358. [Google Scholar]
  42. Zhang, X.L. Study on Wear Performance of Finger Seal Contact Pair. Master’s Thesis, Xi’an University of Technology, Xi’an, China, 2020; pp. 46–52. (In Chinese). [Google Scholar]
  43. Archard, J.F. Elastic Deformation and the Laws of Friction. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1957, 243, 190–205. [Google Scholar]
  44. Lei, J.J.; Liu, M.H.; Hu, X.P.; Wang, J.; Wang, X.L. Analysis of Wear Characteristics of Elliptic Arc Finger Seal Based on Fractal Theory. J. Mech. Eng. 2023, 59, 144–157. [Google Scholar]
  45. Dong, L. A Research for Fractal Geometry Models of Elastic-Plastic Contact and Wear Prediction. Master’s Thesis, Sichuan University, Chengdu, China, 2000; pp. 20–21. (In Chinese). [Google Scholar]
  46. Wei, L.; Gu, B.Q.; Liu, Q.H.; Zhang, P.G.; Fang, G.M. Correction of contact fractal model for friction faces of mechanical seals. CIESC J. 2013, 64, 1723–1729. [Google Scholar]
  47. Greenwood, J.A.; Williamson, J.B.P. Contact of nominally flat surfaces. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1966, 295, 300–319. [Google Scholar]
  48. Sayles, R.S.; Thomas, T.R. Surface topography as a nonst ationary random process. Nature 1978, 271, 431–434. [Google Scholar] [CrossRef]
  49. Majumdar, A.; Bhushan, B. Role of Fractal Geometry in Roughness Characterization and Contact Mechanics of Surfaces. Trans. Asme J. Tribol. 1990, 112, 205–216. [Google Scholar] [CrossRef]
  50. Majumdar, A.; Bhushan, B. Fractal Model of Elastic-Plastic Contact Between Rough Surfaces. J. Tribol. 1991, 113, 1–11. [Google Scholar] [CrossRef]
  51. Xu, G.L.; Zhu, Y.P.; Fang, L.; Du, Y.; Huang, X.M. Three-dimensional Fractal Reconstruction Technique and Leakage Characteristics of Micro-pore Sealing Interfaces. Chin. J. Comput. Phys. 2019, 36, 440–448. [Google Scholar]
  52. Ji, Z.B.; Sun, J.J.; Lu, J.H.; Ma, C.B.; Yu, Q.P. Predicting Method for Static Leakage of Contacting Mechanical Seals Interface Based on Percolation Theory. Tribology 2017, 37, 734–742. [Google Scholar]
  53. Zhou, X.; Pang, H.W.; Liu, H.Y. Leak Rate Prediction of Orbicular Seal Jonint. J. Astronaut. 2007, 28, 762–766. [Google Scholar]
  54. Yu, B.M.; Li, J.H. Some Fractal Characters of Porous Media. Fractals-Complex Geom. Patterns Scaling Nat. Soc. 2001, 9, 365–372. [Google Scholar] [CrossRef]
  55. Yu, B.M.; Xu, P.; Zou, M.Q.; Cai, J.C.; Zheng, Q. Transport Physics in Fractal Porous Media; Science Press: Beijing, China, 2014; pp. 7–8. [Google Scholar]
  56. Li, C.; Shen, Y.Q.; Qiu, S.X.; Yu, B.M.; Xu, P. Mapping Correlation Between Heat Conduction and Fluid Flow Through Unsaturated Porous Media. J. Eng. Thermophys. 2022, 43, 2742–2750. [Google Scholar]
  57. Qiu, S.X.; Xu, P.; Yang, M. The Effective Gas Permeability of Porous Media with Multi-scale Pore Structure. J. Eng. Thermophys. 2019, 40, 1375–1379. [Google Scholar]
  58. Xu, P.; Li, C.H.; Liu, H.C.; Qiu, S.X.; Yu, B.M. Fractal Features of the Effective Gas Transport Coefficient for Multiscale Porous Media. Earth Sci. 2017, 42, 1373–1378. [Google Scholar]
  59. Yu, B.M.; Cheng, P. A Fractal Permeability Model for Bi-dispersed Porous Media. Int. J. Heat Mass Transf. 2002, 45, 2983–2993. [Google Scholar] [CrossRef]
  60. Yu, B.M.; Li, J.H. A Geometry Model for Tortuosity of Flow Path in Porous Media. Chin. Phys. Lett. 2004, 21, 1569–1571. [Google Scholar]
  61. Yu, B.M. Fractal Character for Tortuous Streamtubes in Porous Media. Chin. Phys. Lett. 2005, 22, 158–160. [Google Scholar]
  62. Zhang, Y.C.; Chen, G.D.; Shen, X.L. Analysis of Dynamic Performance and Leakage for Finger Seal. Acta Aeronaut. Astronaut. Sin. 2009, 30, 2193–2199. [Google Scholar]
  63. Wang, H.Y.; Chen, G.D.; Zhang, Y.C. Sensitivity Analysis of Finger Seal Performances and Construction Parameters. Lubr. Eng. 2008, 33, 11–13. [Google Scholar]
  64. Zhang, M.C. Performance Analysis and Optimization of the Extended Rotary Finger Seal. Master’s Thesis, Xi’an University of Technology, Xi’an, China, 2020; pp. 11–36. (In Chinese). [Google Scholar]
  65. Wang, L.N. Research on Dynamic Performance of C/C Composite Finger Seal under Considering Work Status. Ph.D. Thesis, Northwestern Polytechnical University, Xi’an, China, 2016; pp. 99–108. (In Chinese). [Google Scholar]
  66. Lu, F. Study on Wear and Dynamic Performance of C/C Composite Finger Seal. Ph.D. Thesis, Northwestern Polytechnical University, Xi’an, China, 2017; pp. 111–126. (In Chinese). [Google Scholar]
  67. Zhi, B.W. Study on the Effect of Surface Texture on the Tribological Characteristics of Contact Finger Seal. Master’s Thesis, Xi’an University of Technology, Xi’an, China, 2022; pp. 57–63. (In Chinese). [Google Scholar]
  68. Zong, Z.K.; Su, H. Geometric Feature and Structure Optimization of Arc Finger Seal. Acta Aeronaut. Astronaut. Sin. 2010, 31, 393–399. [Google Scholar]
Figure 1. Schematic diagram of the finger seal (FS): (a) general structure, (b) single finger seal annulus (FSA).
Figure 1. Schematic diagram of the finger seal (FS): (a) general structure, (b) single finger seal annulus (FSA).
Fractalfract 07 00494 g001
Figure 2. Design method of arc FS. The red dashed rectangular box and the blue dashed rectangular box are the partially enlarged area; the position where the red dashed arrow and the blue dashed arrow point is the final partially enlarged area.
Figure 2. Design method of arc FS. The red dashed rectangular box and the blue dashed rectangular box are the partially enlarged area; the position where the red dashed arrow and the blue dashed arrow point is the final partially enlarged area.
Fractalfract 07 00494 g002
Figure 3. Main leakage channel of FS.
Figure 3. Main leakage channel of FS.
Fractalfract 07 00494 g003
Figure 4. Roughness surface contour height with different fractal dimension (D): (a) D = 1.4, G = 1 × 10−12; (b) D = 1.5, G = 1 × 10−12; (c) D = 1.6, G = 1 × 10−12.
Figure 4. Roughness surface contour height with different fractal dimension (D): (a) D = 1.4, G = 1 × 10−12; (b) D = 1.5, G = 1 × 10−12; (c) D = 1.6, G = 1 × 10−12.
Fractalfract 07 00494 g004
Figure 5. Roughness surface contour height with different fractal roughness parameter (G): (a) D = 1.3, G = 1 × 10−11; (b) D = 1.3, G = 1 × 10−10; (c) D = 1.3, G = 1 × 10−9.
Figure 5. Roughness surface contour height with different fractal roughness parameter (G): (a) D = 1.3, G = 1 × 10−11; (b) D = 1.3, G = 1 × 10−10; (c) D = 1.3, G = 1 × 10−9.
Fractalfract 07 00494 g005
Figure 6. Side leakage channel of FS.
Figure 6. Side leakage channel of FS.
Fractalfract 07 00494 g006
Figure 7. The finite element calculation model of dynamic performance of FS: (a) The finite element calculation model; (b) Low-pressure FS; (c) High-pressure FS.
Figure 7. The finite element calculation model of dynamic performance of FS: (a) The finite element calculation model; (b) Low-pressure FS; (c) High-pressure FS.
Fractalfract 07 00494 g007
Figure 8. Comparison of calculated and experimental values of leakage. The black line with a black solid square is the experimental result of reference [67], the red line with a red circle is the numerical calculation of reference [68], the blue line with a blue square triangle is the numerical calculation of Equation (42), the green line with a green inverted triangle is the numerical calculation of Equation (18), and the purple line with a purple diamond is the numerical calculation of Equation (41).
Figure 8. Comparison of calculated and experimental values of leakage. The black line with a black solid square is the experimental result of reference [67], the red line with a red circle is the numerical calculation of reference [68], the blue line with a blue square triangle is the numerical calculation of Equation (42), the green line with a green inverted triangle is the numerical calculation of Equation (18), and the purple line with a purple diamond is the numerical calculation of Equation (41).
Fractalfract 07 00494 g008
Figure 9. The influence of working conditions on the FS main leakage rate (Qm): (a) The influence of axial pressure difference Δp on the FS main leakage rate (Qm); (b) The influence of radial displacement excitation Δr on the FS main leakage rate (Qm); (c) The influence of rotor speed n on the FS main leakage rate (Qm).
Figure 9. The influence of working conditions on the FS main leakage rate (Qm): (a) The influence of axial pressure difference Δp on the FS main leakage rate (Qm); (b) The influence of radial displacement excitation Δr on the FS main leakage rate (Qm); (c) The influence of rotor speed n on the FS main leakage rate (Qm).
Fractalfract 07 00494 g009
Figure 10. The influence of fractal dimension (D) on the FS main leakage rate (Qm): (a) The influence of axial pressure difference (Δp) and fractal dimension (D) on the FS main leakage rate (Qm); (b) The influence of rotor speed (n) and fractal dimension (D) on the FS main leakage rate (Qm).
Figure 10. The influence of fractal dimension (D) on the FS main leakage rate (Qm): (a) The influence of axial pressure difference (Δp) and fractal dimension (D) on the FS main leakage rate (Qm); (b) The influence of rotor speed (n) and fractal dimension (D) on the FS main leakage rate (Qm).
Fractalfract 07 00494 g010
Figure 11. The influence of fractal roughness parameter (G) on the FS main leakage rate (Qm): (a) The influence of axial pressure difference (Δp) and fractal roughness parameter (G) on the FS main leakage rate (Qm); (b) The influence of rotor speed (n) and fractal roughness parameter (G) on the FS main leakage rate (Qm).
Figure 11. The influence of fractal roughness parameter (G) on the FS main leakage rate (Qm): (a) The influence of axial pressure difference (Δp) and fractal roughness parameter (G) on the FS main leakage rate (Qm); (b) The influence of rotor speed (n) and fractal roughness parameter (G) on the FS main leakage rate (Qm).
Fractalfract 07 00494 g011
Figure 12. The influence of eccentricity ε on the FS main leakage rate (Qm): (a) The influence of axial pressure difference (Δp) and eccentricity ε on the FS main leakage rate (Qm); (b) The influence of rotor speed (n) and eccentricity (ε) on the FS main leakage rate (Qm).
Figure 12. The influence of eccentricity ε on the FS main leakage rate (Qm): (a) The influence of axial pressure difference (Δp) and eccentricity ε on the FS main leakage rate (Qm); (b) The influence of rotor speed (n) and eccentricity (ε) on the FS main leakage rate (Qm).
Fractalfract 07 00494 g012
Figure 13. The influence of axial pressure difference (Δp), tortuosity fractal dimension (DT) after loading and the pore fractal dimension (Df0) after loading on the FS side leakage rate (Qs): (a) The influence of axial pressure difference (Δp) and tortuosity fractal dimension (DT) on the FS side leakage rate (Qs); (b) The influence of axial pressure difference (Δp) and the pore fractal dimension (Df0) after loading on the FS side leakage rate (Qs).
Figure 13. The influence of axial pressure difference (Δp), tortuosity fractal dimension (DT) after loading and the pore fractal dimension (Df0) after loading on the FS side leakage rate (Qs): (a) The influence of axial pressure difference (Δp) and tortuosity fractal dimension (DT) on the FS side leakage rate (Qs); (b) The influence of axial pressure difference (Δp) and the pore fractal dimension (Df0) after loading on the FS side leakage rate (Qs).
Fractalfract 07 00494 g013
Figure 14. The influence of axial pressure difference (Δp) and porosity (ϕ) on the FS side leakage rate (Qs).
Figure 14. The influence of axial pressure difference (Δp) and porosity (ϕ) on the FS side leakage rate (Qs).
Fractalfract 07 00494 g014
Figure 15. Influence of working conditions on leakage performance of FS: (a) Influence of axial pressure difference (Δp) on leakage performance of FS; (b) Influence of radial displacement excitation (Δr) on leakage performance of FS; (c) Influence of rotor speed (n) on leakage performance of FS.
Figure 15. Influence of working conditions on leakage performance of FS: (a) Influence of axial pressure difference (Δp) on leakage performance of FS; (b) Influence of radial displacement excitation (Δr) on leakage performance of FS; (c) Influence of rotor speed (n) on leakage performance of FS.
Fractalfract 07 00494 g015
Figure 16. Influence of fractal feature parameters on leakage performance of FS: (a) Influence of fractal dimension (D) on leakage performance of FS; (b) Influence of fractal roughness parameter (G) on leakage performance of FS.
Figure 16. Influence of fractal feature parameters on leakage performance of FS: (a) Influence of fractal dimension (D) on leakage performance of FS; (b) Influence of fractal roughness parameter (G) on leakage performance of FS.
Fractalfract 07 00494 g016
Figure 17. Influence of porosity (ϕ) on leakage performance of FS.
Figure 17. Influence of porosity (ϕ) on leakage performance of FS.
Fractalfract 07 00494 g017
Table 1. Structural parameters of arc finger seal (FS).
Table 1. Structural parameters of arc finger seal (FS).
ParametersUnitsValues
Radius of the outer circle/Dwmm77
Radius of the root circle/Dfmm68.75
Radius of the inner circle/Drmm60.5
Radius of the base circle/Dccmm11
Arc radius/Rcmm63
Rotor radius/rmm59.5
Gap width/Iamm2
Clearance angle/α°0.2
Downstream protection height/gdmm1.2
Height of the heel/xgmm1.2
Number of finger beams/Nfb 42
Table 2. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different axial pressure differences.
Table 2. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different axial pressure differences.
Axial Pressure Differences
(MPa)
Main Leakage Rate Ratio
(%)
Side Leakage Rate Ratio
(%)
0.197.902.10
0.1598.631.37
0.298.711.29
0.2598.771.23
0.399.090.91
0.3599.250.75
0.499.390.61
0.4599.510.49
0.599.590.41
Table 3. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different radial displacement excitation values.
Table 3. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different radial displacement excitation values.
Radial Displacement Excitation
(mm)
Main Leakage Rate Ratio
(%)
Side Leakage Rate Ratio
(%)
0.0196.213.79
0.0296.963.04
0.0396.973.03
0.0497.542.46
0.0598.031.97
0.0698.311.69
0.0798.571.43
0.0898.781.22
0.0998.781.22
Table 4. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different rotor speed values.
Table 4. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different rotor speed values.
Rotor Speed
(r/min)
Main Leakage Rate Ratio
(%)
Side Leakage Rate Ratio
(%)
10,00098.721.28
10,50098.711.29
11,00098.601.40
11,50098.491.51
12,00098.351.65
12,50098.181.82
13,00097.972.03
13,50097.702.30
14,00097.352.65
14,50096.893.11
15,00096.303.70
Table 5. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different fractal dimensions.
Table 5. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different fractal dimensions.
Fractal DimensionMain Leakage Rate Ratio
(%)
Side Leakage Rate Ratio
(%)
1.397.972.03
1.497.982.02
1.597.992.01
1.698.011.99
1.798.051.95
1.898.111.89
1.998.291.71
Table 6. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different fractal roughness parameters.
Table 6. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different fractal roughness parameters.
Fractal Roughness ParameterMain Leakage Rate Ratio
(%)
Side Leakage Rate Ratio
(%)
1 × 10−797.882.12
1 × 10−897.972.03
1 × 10−997.972.03
1 × 10−1097.972.03
1 × 10−1197.972.03
1 × 10−1297.972.03
Table 7. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different porosity values.
Table 7. The proportion of the FS main leakage rate and the FS side leakage rate to the FS total leakage rate under different porosity values.
PorosityMain Leakage Rate Ratio
(%)
Side Leakage Rate Ratio
(%)
0.197.252.75
0.297.092.91
0.396.993.01
0.496.923.08
0.596.863.14
0.696.813.19
0.796.773.23
0.896.743.26
0.996.703.30
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

Lei, J.; Liu, M.; Chang, W.; Wan, Y. Analysis of the Total Leakage Characteristics of Finger Seal Considering Fractal Wear and Fractal Porous Media Seepage Effects. Fractal Fract. 2023, 7, 494. https://doi.org/10.3390/fractalfract7070494

AMA Style

Lei J, Liu M, Chang W, Wan Y. Analysis of the Total Leakage Characteristics of Finger Seal Considering Fractal Wear and Fractal Porous Media Seepage Effects. Fractal and Fractional. 2023; 7(7):494. https://doi.org/10.3390/fractalfract7070494

Chicago/Turabian Style

Lei, Junjie, Meihong Liu, Wei Chang, and Yongneng Wan. 2023. "Analysis of the Total Leakage Characteristics of Finger Seal Considering Fractal Wear and Fractal Porous Media Seepage Effects" Fractal and Fractional 7, no. 7: 494. https://doi.org/10.3390/fractalfract7070494

APA Style

Lei, J., Liu, M., Chang, W., & Wan, Y. (2023). Analysis of the Total Leakage Characteristics of Finger Seal Considering Fractal Wear and Fractal Porous Media Seepage Effects. Fractal and Fractional, 7(7), 494. https://doi.org/10.3390/fractalfract7070494

Article Metrics

Back to TopTop