Next Article in Journal
Tribological Properties of Cu-Based Self-Lubricating Materials Composed with Cu- or Ni-Plated Graphite
Next Article in Special Issue
Research on Internal Flow Field Characteristics of Straight-Groove Texture Using Three-Dimensional Modeling
Previous Article in Journal
Tire Wear Sensitivity Analysis and Modeling Based on a Statistical Multidisciplinary Approach for High-Performance Vehicles
Previous Article in Special Issue
Investigation of the Static Performance of Hydrostatic Thrust Bearings Considering Non-Gaussian Surface Topography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Effect of Fluid Cavitation and Inertia on the Pressure Buildup of Parallel Textured Surfaces

College of Petrochemical Engineering, Lanzhou University of Technology, Lanzhou 730050, China
Lubricants 2023, 11(7), 270; https://doi.org/10.3390/lubricants11070270
Submission received: 6 May 2023 / Revised: 7 June 2023 / Accepted: 17 June 2023 / Published: 21 June 2023
(This article belongs to the Special Issue Hydrodynamic Lubrication of Textured Surfaces)

Abstract

:
A mathematical model is developed to investigate the combined effect of fluid cavitation and inertia on the fluid pressure buildup of parallel textured surfaces. The fluid cavitation is analyzed using the Rayleigh–Plesset model, and the fluid inertia is analyzed with an averaged method. The finite element method and Newton-downhill method are employed to solve the governing equations. The numerical model is validated by comparing the experimental and numerical results, and the combined effect of fluid cavitation and inertia on the fluid pressure buildup is analyzed and discussed. The research indicates that the cavitation weakens the fluid inertia effect on the pressure distribution at the inlet area of textures. The fluid inertia greatly enhances the hydrodynamic effect and effectively limits the excessive extension of the low-pressure zone caused by cavitation. The fluid cavitation and inertia, especially their interaction, significantly affect the fluid pressure buildup and generate a net load-carrying capacity (LCC). The numerical model with the fluid inertia and cavitation is more time saving than the commercial CFD tools in solutions, which gives a novel and optional HD foundation for developing a more efficient and accurate THD or TEHD model by numerical programming.

1. Introduction

In hydrodynamic lubrication, the fluid cavitation and inertia are two common special effects, especially at the location where the fluid film thickness changes suddenly. The fluid film cavitation is complicated and attracts researchers’ attention more often. In the review of fluid film cavitation [1], Braun et al. outlined three recognized forms of fluid cavitation: gaseous cavitation, pseudo-cavitation and vaporous cavitation. Several studies [2,3,4,5,6] have shown that it affects the fluid pressure buildup, so that the fluid film does not appear negative pressure or too low pressure. Studies [7,8,9] on the fluid inertia of textured surfaces showed that it can significantly affect the fluid pressure buildup due to the side wall. Their effects on the fluid pressure buildup of textured surface have been often investigated separately. In some cases, fluid cavitation and inertia, however, are simultaneously encountered, and they affect the pressure distribution of fluid film together as well as influence each other. Therefore, it is necessary to take into account the above two effects simultaneously for these cases.
Nowadays, the Reynolds equation (RE) combined with JFO boundary conditions (RE-JFO) [10,11,12] is a widely accepted mass-conserving model in analyzing textured surfaces. Singhal et al. [13] developed a full cavitation model (FCM) and incorporated it into the CFD software. Then, Bakir et al. [14] outlined the volume fraction (VOF) model in the CFD software and analyzed the cavitation behavior in an inducer. The simulation results are in good agreement with the experimental results. Both of the cavitation models above are based on the Rayleigh–Plesset equation (RPE) in the dynamics theory of cavitation bubbles [15]. Remarkably, they satisfy the mass conservation. Compared with the JFO theory, it can better approximate the physical process of fluid cavitation. Several researchers [16,17,18,19,20,21] have successfully analyzed the hydrodynamic lubrication in journal bearings by the similar cavitation models based on RPE. Recently, Wang et al. [22] presented a homogeneous phase change model by combining the RPE and RE to investigate the boiling flow in wavy face seals. It can be observed that the RPE model is a feasible means to analyze the cavitation flow of fluid film.
To study the fluid inertia effect on the pressure distribution of textured surfaces, several researchers [7,8] have solved the Navier–Stokes equations (NSEs) by the CFD technology. Despite its advantage in the detailed investigation of fluid inertia, the numerical calculation based on NSEs is very time consuming. Another way to deal with fluid inertia is based on the hypothesis proposed by Constantinescu and Galetuse [23]. They assume that fluid inertia only affects the velocity magnitude and does not affect the velocity profile. This method has been successfully used to study the fluid inertia effect on the hydrostatic seals [24], in which the seal face is a two-dimensional axisymmetric structure. In general, it has a better approximation to the fluid inertia for the shallow textures and grooves due to the limitation of the assumption.
The above short review emphasizes the problem in the studies on the fluid cavitation and inertia of textured surfaces as well as their treatments and solutions. Moreover, the research objects are often tribological components with surface texture, such as seals [25,26], bearings [27,28], piston rings [29] and so on. Recesses of various shapes, as ones of the typical surface textures, have been widely used in the above fields. In order to study the mechanism conveniently, the sliding cell shown in Figure 1 is taken as the research object instead of a specific tribological component.
In this work, a mathematical model with fluid cavitation and inertia based on the RPE [15] and hypothesis of Constantinescu and Galetuse [23] is proposed for the parallel surfaces with a shallow recess. A numerical procedure is developed to solve the governing equations by the finite element method and Newton-downhill technology. The numerical model is then validated by confronting it with the previously published experimental results and the numerical ones by RE-JFO and CFD. The combined effect of fluid cavitation and inertia on the fluid pressure buildup and their interaction mechanism are analyzed and discussed.

2. Theoretical Models

2.1. Geometrical Model

The geometrical model of the parallel sliding surfaces with a shallow recess is shown in Figure 1a,b. In the model, the lower surface moves along the x direction with a velocity u0, while the upper surface cell with a side length W0 is stationary. A shallow recess is located in the center of the upper surface cell, where its shapes include a square, circle, and triangle (shown in Figure 1c), it has the same length W1 along the x and y directions, and its depth is denoted by h1. The gap between two surfaces is filled with the lubricant #15 aviation hydraulic oil with the film thickness h0 in the non-textured area. The film thickness between the parallel surfaces is given as
h = h 0 non - textured   region h 0 + h 1 textured   region
The geometric and operating parameters used in the paper are listed in Table 1.

2.2. Mathematical Model

The following assumptions are made in the present model.
  • The roughness of parallel surfaces is neglected, and the slip flow between the surfaces and the fluid film are not considered.
  • The flow regime of the fluid film is laminar, as calculated from Re.
  • The fluid inertia only affects the velocity magnitude and does not affect the velocity profile [23,24].
  • The fluid film in the cavitation zone is a homogeneous mixture of liquid and vapor, and the mixture is assumed to be isothermal and incompressible [14,22,30].
  • The density of liquid is considered constant, the vapor is assumed to be the ideal gas, and the non-condensable gas in the fluid film is not considered. The relative motion between them is neglected and they share the same flow rates and pressure [14,22,30].

2.2.1. Mixture Mass and Momentum Equations

The equation of mass conservation for the steady flow can be expressed as follows.
u x + v y + w z = 0
where u, v and w are the velocity components in the x, y and z coordinate systems. Since the film thickness in the z direction is in the micrometer level, much smaller than the scales in the x and y directions, the pressure variation along the film thickness is neglected. The NSEs can be written as
ρ u u x + v u y + w u z = p x + μ 2 u z 2
ρ u v x + v v y + w v z = p y + μ 2 v z 2
0 = p z
where p, ρ and μ are the pressure, density and dynamic viscosity of the mixture, respectively. The physical property parameters ρ and μ of the mixture are to be defined in Section 2.2.2. Equations (3) and (4) are rewritten in the following forms based on the mass conservation shown in Equation (2).
ρ u 2 x + v u y + w u z = p x + μ 2 u z 2
ρ u v x + v 2 y + w v z = p y + μ 2 v z 2
Integrating Equations (2), (6) and (7) in the z direction, where the following forms can be obtained (note that the velocity component w along the z axis is 0),
x 0 h u d z + y 0 h v d z = 0
ρ x 0 h u 2 d z + y 0 h u v d z = p x h + μ u z 0 h
ρ x 0 h u v d z + y 0 h v 2 d z = p y h + μ v z 0 h
According to Refs. [23,24], the velocity profiles of inertial flow can be deduced as
u = 6 q x h 2 z + u 0 1 z h v = 6 q y h 2 z 1 z h
The variables qx and qy are the flow rates of unit length along the x and y directions. Introducing Equation (11) into Equations (8)–(10) gives
q x x + q y y + u 0 2 h x = 0
ρ I x x x + I x y y + h p x + 12 μ h 2 q x = 0
ρ I x y x + I y y y + h p y + 12 μ h 2 q y = 0
where Ixx, Ixy, and Iyy are the integrals of the products of velocity components along the z direction, which are defined as follows.
I x x = 0 h u 2 d z = 6 5 h q x 2 + q x u 0 + 1 3 h u 0 2
I x y = 0 h u v d z = 6 5 h q x q y + 1 2 q y u 0
I y y = 0 h v 2 d z = 6 5 h q y 2

2.2.2. Mass Transfer Equations

According to the theory of homogeneous two-phase flow [22,30], the physical property parameters ρ and μ of the mixture are defined as follows.
ρ = F L ρ L + F G ρ G
μ = F L μ L + F G μ G
where ρL and ρG are the densities of the liquid and vapor, ρL is constant and ρG is given by the ideal gas state equation. FL and FG are the volume fractions of the liquid and vapor per unit volume in the mixture, which can be defined as
F L = V L V L + V G F G = V G V L + V G
where VL and VG are liquid and vapor volumes, respectively.
The liquid flow and vapor flow are governed by the following split phase mass transfer equations:
x F L ρ L u + y F L ρ L v + z F L ρ L w = S x F G ρ G u + y F G ρ G v + z F G ρ G w = S
where S is the cavitation source term and is defined in Section 2.2.3. According to the assumption (5) and Section 2.2.1, the w value is 0. Integrating Equation (21) along the film thickness direction (z axis), the integral forms can be expressed as
x F L ρ L 0 h u d z + y F L ρ L 0 h v d z = S h x F G ρ G 0 h u d z + y F G ρ G 0 h v d z = S h
The integral forms of the mixture velocity components in Equation (11) across the film thickness are as follows
0 h u d z = 0 h 6 q x h 2 z + u 0 1 z h d z = q x + u 0 h 2 0 h v d z = 0 h 6 q y h 2 z 1 z h d z = q y
Introducing Equation (23) into Equation (22), we can obtain
x F L ρ L q x + u 0 h 2 + y F L ρ L q y = S h x F G ρ G q x + u 0 h 2 + y F G ρ G q y = S h
Equation (24) can be rewritten as
F L x q x + u 0 h 2 + F L y q y 1 F L F L p p x q x + u 0 h 2 + p y q y 1 F L ρ L + F L ρ G S h = 0
where there is only one unknown variable FL, because FG is replaced by 1 − FL.

2.2.3. Cavitation Source Term: RPE Model

The RPE in the dynamics theory of cavitation bubbles [15] can well describe the growth and rupture of cavitation bubbles. Without considering the effects of viscous damping and surface tension, the first-order approximation of RPE [13,14,15] can be written as
Ŕ b = 2 3 p p c ρ L
where pc is the cavitation pressure of the fluid film, and Ŕb is the dynamic radius of cavitation bubbles. It is assumed that bubbles grow from an initial average radius of Rb when the liquid film ruptures and the bubbles change to the initial ones when the liquid film regenerates. The initial average radius Rb can be assumed to be 1 μm [14]. When the liquid film ruptures and regenerates, cavitation source terms can be written as
S = C 1 N ρ G 4 π R b 2 Ŕ b p < p c C 2 N ρ G 4 π R b 2 Ŕ b p > p c
where C1 and C2 are the empirical coefficients of vaporization and liquefaction, respectively. The values provided in Ref. [22] are used.
C 1 = C 2 = 0.05 p p c p c
N is the number of initial bubbles per unit volume in the mixture. It has been given as follows in Ref. [14].
N = 3 F L 4 π R b 3 p < p c 3 F G 4 π R b 3 p > p c
Then, the source term in Equation (27) can be written uniformly as
S = T ρ G 4 π R b 2 2 3 p p c ρ L sgn p p c
where T is a comprehensive coefficient, and it makes vaporization and liquefaction source terms to be a unified formula.
T = 0.05 p p c p c p p c + p p c 2 p c 3 F L 4 π R b 3 + p p c + p p c 2 p c 3 1 F L 4 π R b 3
Equation (31) has already been presented in Ref. [22] and applied directly in the present model.

2.2.4. Boundary Conditions

The governing equations of the present model include Equations (12)–(14) and Equation (25). The unknown variables include p, qx, qy and FL. To solve the governing equations, the periodic boundary conditions shown in Equation (32) and the imposed boundary conditions shown in Equation (33) are employed.
q x x = 3 = q x x = 3
q y x = 3 = q y x = 3
p x = 3 = p x = 3
F L x = 3 = F L x = 3
p y = 3 = p y = 3 = p a
F L y = 3 = F L y = 3 = 1

3. Numerical Procedure

The numerical model is solved by the numerical programming with the finite element method. The Newton-downhill method is applied to solve the nonlinear Equation (12)–(14) and (25). The weighted residual forms of governing equations are as follows.
R c i = Ω w i q x x + q y y + u 0 2 h x d x d y
R m x i = Ω w i ρ I x x x + I x y y + h p x + 12 μ h 2 q x d x d y
R m y i = Ω w i ρ I x y x + I y y y + h p y + 12 μ h 2 q y d x d y
R F L i = Ω w i F L x q x + u 0 h 2 + F L y q y 1 F L F L p p x q x + u 0 h 2 + p y q y 1 F L ρ L + F L ρ G S h d x d y
where wi is a weight function on the two-dimensional computational domain Ω. When the numerical solution converges, the above residuals tend to zero. Numerical iterative systems with the Newton-downhill method are as follows.
R m x i q x j R m x i q y j R m x i p j R m y i q x j R m y i q y j R m y i p j R c i q x j R c i q y j 0 δ q x j δ q y j δ p j = λ q p R m x i R m y i R c i
R F L i F L j δ F L j = λ F R F L i
Then, the variables are expressed as
q x ( k ) = q x ( k 1 ) + δ q x ( k ) q y ( k ) = q y ( k 1 ) + δ q y ( k ) p ( k ) = p ( k 1 ) + δ p ( k )
F L ( k ) = F L ( k 1 ) + δ F L ( k )
where k is the number of iterations. The procedure convergence criterion is as follows, where X includes p, qx, qy and FL.
X ( k ) X ( k - 1 ) X ( k ) < 10 4
The LCC and flow rate can be expressed as
F o = Ω p d x d y
Q = l ( q x l x + q y l y ) d l
where l is the integral curve, and lx and ly are the unit vectors along the x and y directions.
The numerical procedure is described as follows.
(1).
Input geometric and operating parameters.
(2).
Give initial values X(0) (p(0), qx(0), qy(0), FL(0)) to unknown variables X (p, qx, qy, FL), and impose pre-boundary conditions.
(3).
Calculate the physical property parameters ρ, μ and flow characteristic parameters Ixx, Ixy, and Iyy.
(4).
Solve Equation (38) and obtain the pressure and flow rate distributions.
(5).
Solve Equation (39) based on step (4) and obtain the liquid volume fraction distribution.
(6).
Check the convergence criterion. If the unknowns meet the tolerance, end the iteration; otherwise, go to step (3) and repeat the procedure.
(7).
Output results and post-processing.

4. Results and Discussion

4.1. Validation

To validate the mathematical model and computational procedure, the results given by the present model are compared with the published experimental results of the journal bearing in Ref. [31]. They are further compared with the numerical results from the RE-JFO model, which is solved with the method in Ref. [32]. Figure 2 shows the pressure profiles at the centerline of the journal bearing. It shows that the pressure peak from the present model is higher than that from the RE-JFO model. This may be the effect of the fluid inertia. The cavitation zone obtained from the present model is narrower than that from the RE-JFO model. It seems that the pressure distribution from the present model is more consistent with the experimental results than the RE-JFO model. In order to further verify the cavitation prediction of textured surfaces, Figure 3 shows the pressure distribution of the circle recess with a side length W1 = 3 mm at a very small Re = 0.265 and a small gap thickness h0 = 10 μm, where the Reynolds number Re is defined as ρLu0h0/μL. It is worth noting that the pressure in the low-pressure zone given by the present model changes continuously near the cavitation pressure pc (0.08 MPa) rather than the fixed value. The main reason is that in the present model, the pc is only used to calculate the source term S, and the pressure in the cavitation region is solved by the mass conservation equation and momentum equations. It is in accordance to many experimental research papers on cavitation [1,33]. In addition, according to the Bernoulli equation, when the fluid flows through the recess, the pressure should first decrease and then increase. Although the fluid cavitation makes the pressure in the cavitation region greatly increase, its trend should be consistent. In general, the pressure distribution given by the present model is in good agreement with the one given by the RE-JFO model. The above verification analysis shows that the present model can well solve the cavitation problem of textured surfaces.
Figure 4 shows the fluid pressure distribution of the square recess with a side length W1 = 3 mm at y = 0 when the fluid cavitation is neglected for the case of Re = 148.82. The pressure distribution from the present model is compared with the results by solving N-S equations from ANSYS FLUENT (NSEs) and by solving the Reynolds equation (RE). Due to the change of fluid film thickness, a low-pressure zone is generated around the zone of x = −1.5, and a high one is generated around the zone of x = 1.5. When the fluid inertia is not considered, that is the result from the RE, and the fluid pressure distribution shows triangular profiles. For the case of NSEs and the present model, there are pressure jumps at x = −1.5 and x = 1.5 due to the fluid inertia. The pressure peak values are larger and valley values are smaller than those from the RE. That means the fluid inertia plays a remarkable role in the fluid pressure distribution. It is worth noting that the results from the present model and the NSEs agree well with each other. It proves that the present model is validated.
In this paper, the triangular elements are used to discretize the computational domain. The meshes in the recess region are denser than those in the other region because the cavitation may occur there (shown in Figure 5a). The mesh independence of the present model is checked, and the results are shown in Figure 5b. It indicates that with the increase in the node number, the LCC of the fluid film and the flow rates on the periodic boundaries gradually tend to their stable values.

4.2. Effect of Re

Figure 6 shows the pressure distribution of fluid film considering only fluid inertia or both fluid inertia and cavitation at Re = 148.82. It can be seen from Figure 6a that due to the fluid inertia and side wall [7,8,9], there are significant pressure jumps at the inlet and outlet of texture, that is, the absolute values of the pressure gradients of these two regions are much larger than those of the other regions. Figure 6b shows that due to the fluid cavitation, the lower pressure at the inlet is greatly enhanced, and there is no change in the peak pressure. Figure 7 further illustrates the effect of fluid cavitation and Re on the pressure distribution of fluid film (where “I + C” is inertia and cavitation; “I” is only inertia, and “C” is only cavitation). It is worth noting that the fluid cavitation deceases the pressure jump at the inlet of the texture. The main reason is the decrease in fluid density in the cavitation region, where the fluid is a mixture of the liquid and vapor. It is well known that the mass is the only measure of inertia. As a result, a smaller fluid inertia causes a smaller pressure jump. In addition, with the increase in Re, the pressure in the low-pressure zone is gradually reduced. It is in accordance to experimental research papers on cavitation [1,33]; according to Braun et al. and Etsion et al., the greater the speed, the lower the pressure in the cavitation region. Moreover, it also can be found that the fluid cavitation improves the valley pressure value; however, the peak value is hardly affected.
Figure 8 presents the effects of fluid inertia and Re on the pressure distribution. When the fluid inertia is not considered, with the increase in Re, the hydrodynamic effect is enhanced slowly, while the low-pressure zone caused by the fluid cavitation increases gradually. When the fluid inertia is considered, with the increase in Re, the hydrodynamic effect is rapidly enhanced. It shows that the fluid inertia significantly affects the local pressure distribution of fluid film. Compared with the results without the fluid inertia, when the Re is 115.75, 148.82, and 181.89 in turn, the incremental percentage ∆p of the pressure peak is 18%, 28%, and 40%, respectively. When the fluid inertia is not considered, the peak pressure is much smaller. The greater the Re, the more serious the deviation is. At the same position of the longitudinal axis shown in Figure 8, the incremental percentage ∆x of the high-pressure zone width is 133%, 161% and 196%, respectively. When the fluid inertia is not considered, the width of the high-pressure zone is greatly reduced; that is, the width of the low one is greatly increased. The greater the Re, the more obvious the change. The main reason is that the excessive expansion of the low-pressure zone caused by the fluid cavitation greatly weakens the hydrodynamic effect of texture.
Figure 9 shows the effect of Re on the LCC Fo; there are five cases: flat surface (Reynolds equation, RE-Flat); only texture (Reynolds equation, RE); texture and fluid cavitation (RE-JFO, C); texture and fluid inertia (present model only considering inertia, I); texture, fluid inertia and cavitation (present model considering both inertia and cavitation, I + C). When the fluid inertia is not considered, the LCCs of RE-Flat, RE and C are almost identical and unaffected by the Re. The main reason is that the pressure distribution is anti-symmetric for the case RE, and for the case C, although fluid cavitation avoids extremely small fluid pressure, it makes the high-pressure zone be reduced greatly. When the fluid inertia is only considered (I), with the increasing Re, the fluid inertia makes the pressure distribution of fluid film asymmetrical and then produces a small net LCC. This is consistent with many research papers on fluid inertia [7,8,9]. When both the fluid inertia and cavitation are considered, the fluid cavitation avoids extremely small fluid pressure, and the fluid inertia effectively enhances the hydrodynamic effect of texture and limits the excessive expansion of the low-pressure zone. As a result, a net LCC is generated, and it increases with the Increase in Re. When the Re value is 181.89, the net increment of the LCC is about 7.1% compared with the cases without fluid inertia. It is worth noting that the combined effect of fluid cavitation and inertia on the LCC is more significant than the case with only one of them. The main reason is the interaction between fluid inertia and cavitation, that is, their combined effect on the fluid pressure buildup.

4.3. Effect of pc

Figure 10 presents the effect of cavitation pressure pc on the fluid film pressure distribution for the case of Re = 165.35. For the current Re and gap thickness, the fluid inertia has a significant influence on the pressure distribution. It also be seen that the pressure jump at the inlet of the texture is far less than that of the texture outlet due to the smaller density of mixture at the inlet zone. For these two cases, with the increase in pc, the pressure in the low-pressure zone increases rapidly. When the fluid inertia is not considered, with the increase in pc, the low-pressure zone increases and the peak pressure decreases rapidly. The main reason is that the cavitation becomes more severe with the increasing pc, and the excessive expansion of the low-pressure zone caused by the fluid cavitation greatly weakens the hydrodynamic effect of the texture. When the fluid inertia is considered, as the pc increases, the low-pressure zone increases slowly. Another important finding is that the pc has little effect on the pressure peak or the pressure jump at the texture outlet. As a result, when pc is 0.065, 0.08, and 0.095, in turn, the incremental percentage ∆p of the pressure peak with the fluid inertia is 10%, 15%, and 33%, respectively. At the same position of the longitudinal axis shown in Figure 10, the incremental percentage ∆x of the high-pressure zone width is 13%, 29% and 149%, respectively. The main reason is that the fluid inertia significantly enhances the hydrodynamic effect of texture; thus, it greatly limits the excessive expansion of the low-pressure zone. The above analysis shows that the pressure jump at the texture outlet significantly affects the whole fluid film pressure buildup, and the one at the texture inlet mainly affects the local pressure distribution. According to the experimental results [34,35] of Sun et al., the pressure in the cavitation region is close to 0, and according to the experimental results [1,33] of Braun et al. and Etsion et al., it is near the atmospheric pressure. It indicates that the cavitation pressure is uncertain for actual working conditions. Therefore, for a working tribological component, the cavitation pressure can affect the fluid pressure buildup and LCC, especially when it has a larger value.
Figure 11 shows the effect of pc on the LCC Fo. For these cases without the fluid cavitation (I, RE, RE-Flat), there is only one data point. When the fluid inertia is not considered (RE-Flat, RE, C), the LCC is almost fixed. Considering only the fluid inertia, the LCC has a small net increment. When both the fluid cavitation and inertia are considered, a considerable net increment of LCC is generated: the greater the pc, the more significant the enhancement of LCC. When the pc value is 0.095 (near the atmospheric pressure), the net increment of the LCC is about 6%. It is the contribution of the interaction of fluid cavitation and inertia; that is, the pressure in the low-pressure zone is significantly enhanced due to a larger pc and the excessive expansion of the low-pressure zone is effectively limited due to the fluid inertia.

4.4. Effect of Sr, h1 and Shapes

Figure 12 shows the effects of texture area ratio Sr on the pressure distribution and LCC of the fluid film. The Sr value is a ratio of the texture area to the area of the upper surface cell, in which the texture area is adjusted by changing the texture size W1. As shown in Figure 12a, with the increase in Sr, the hydrodynamic effect of the texture first increases and then decreases. The main reason is the location changes of the texture inlet and outlet. It can be seen from Figure 12b that when the fluid inertia is not considered (RE, C), the Sr has little effect on the LCC. When the fluid inertia is only considered (I), the LCC increases slowly with the increase in Sr. When both the fluid cavitation and inertia are considered (I + C), the LCC increases firstly and then decreases with the increasing Sr, and when the Sr is about 0.44, the net increment of LCC is about 7.1%. Figure 13 shows the effect of texture depth h1 on the pressure distribution and LCC of fluid film. As shown in Figure 13a, for the two cases, the hydrodynamic effect is gradually enhanced with the increase in h1. It is worth noting that when the fluid inertia is not considered, the enhancement of hydrodynamic effect caused by the increase in h1 is very weak. In addition, as the h1 increases, the pressure jump generated by the fluid inertia becomes more significant. It is consistent with many studies [7,8,9] on the fluid inertia. The main reason is that the increase in h1 enhances the effect of side walls. As shown in Figure 13b, it is also found that the h1 value has the most significant effect on the LCC when both the fluid cavitation and inertia are considered. In the range of calculation, the maximum net increment of LCC is about 8%. Figure 14a–c show the effect of texture shape on the fluid film pressure distribution, where the texture shape includes a square, circle, and triangle, respectively. The pressure distributions of the square and the circle are similar, because they are symmetrical in the sliding direction. For the triangular shape, the pressure jump at the texture inlet is less than that at the outlet. The main reason is that the two side walls of the texture inlet are not perpendicular to the sliding direction. For these three shapes, the combined effect of fluid cavitation and inertia greatly enhances the hydrodynamic effect. As shown in Table 2, the combined effect has the most significant influence on the pressure buildup and LCC. Moreover, the square provides the maximum net LCC. The above analysis indicates that when only the fluid inertia or cavitation, or neither of them is considered, the enhancement of hydrodynamic effect is smaller and the contribution to LCC is very small or almost 0 for Sr, h1 and recess shapes. Only when they are considered simultaneously can the adjustment of Sr, h1 and recess shapes maximize the combined effect of fluid cavitation and inertia on the fluid pressure buildup and obtain the maximum LCC.

5. Conclusions

A mathematical model with fluid cavitation and inertia was developed to investigate the fluid pressure buildup mechanism of parallel textured surfaces. The fluid cavitation was dealt with by the Rayleigh–Plesset model, and the fluid inertia was treated by an averaged method. The mass conservation equation, the momentum equations and the liquid volume fraction equation were solved with the finite element method. The Newton-downhill method was employed to solve the nonlinear equations. The numerical model was validated by comparing with the previously published experimental results and the numerical results by the RE-JFO and ANSYS FLUENT models. The combined effect of fluid cavitation and inertia on the fluid pressure buildup and LCC of hydrodynamic lubrication in parallel textured surfaces was analyzed and discussed.
Due to the fluid inertia and side wall, the pressure jumps occur at the texture inlet and outlet. The fluid cavitation significantly reduces the pressure jump at the texture inlet; thus, the fluid pressure distribution at the inlet is less affected by the fluid inertia.
The pressure jump at the outlet greatly increases the peak pressure and effectively limits the excessive extension of the low-pressure zone caused by the fluid cavitation, which greatly increases the width of the high-pressure zone and significantly enhances the hydrodynamic effect.
The fluid cavitation and inertia, especially their interaction, significantly affect the fluid film pressure buildup and generate a net LCC. The cavitation pressure can affect the fluid pressure buildup and LCC, especially when it has a larger value. The adjustment of recess area ratio Sr, recess depth h1 and recess shapes can maximize the combined effect of fluid cavitation and inertia on the fluid pressure buildup and obtain the maximum LCC.

Funding

This research was funded by the National Natural Science foundation of China (No. 52005236), the Natural Science foundation of Gansu Province (No. 22JR5RA289), and the Hongliu Distinguished Young Talent Support Program of Lanzhou University of Technology (No. 02/062205).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The author declares no conflict of interest.

Nomenclature

ŔbDynamic radius of cavitation bubbles, mReReynolds number, Re = ρLu0h0/μL
C1, C2Empirical coefficient of vapor liquid transition SCavitation source term, kg/(m3∙s)
FG, FLVapor and liquid volume fraction per unit volumeTComprehensive coefficient
FoLoad-carrying capacity, NuVelocity in x direction, m/s
hFluid film thickness, mu0Sliding velocity in x direction, m/s
h0, h1Gap thickness and square recess depth, mvVelocity in y direction, m/s
Ixx, Ixy, IyyFlow characteristic parameters, m3/s2 v0Sliding velocity in y direction, m/s
lx, lyUnit vectors in x and y directionsVG, VLVapor and liquid volume
NNumber of initial bubbles per unit volumewVelocity in z direction, m/s
pFluid film pressure, PaW0, W1Slider and square recess length, m
pa, pcAmbient pressure and cavitation pressure, Pax, y, zCartesian coordinate, m
QFlow rate, m3/sλF, λqpDownhill factors
qx, qyFlow rates of unit length in x and y directions, m2/sμ, μG, μ LMixture, vapor and liquid dynamic viscosity, Pa∙s
RbBubbles initial average radius, mρ, ρG, ρ LMixture, vapor and liquid density, kg/m3

References

  1. Braun, M.J.; Hannon, W.M. Cavitation formation and modelling for fluid film bearings: A review. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2010, 224, 839–863. [Google Scholar] [CrossRef]
  2. Salant, R.F.; Homiller, S.J. Stiffness and leakage in spiral groove upstream pumping mechanical seals. Tribol. Trans. 1993, 36, 55–60. [Google Scholar] [CrossRef]
  3. Etsion, I.; Burstein, L. A model for mechanical seals with regular microsurface structure. Tribol. Trans. 1996, 39, 677–683. [Google Scholar] [CrossRef]
  4. Qiu, Y.Y.; Khonsari, M.M. Performance analysis of full-film textured surfaces with consideration of roughness effects. J. Tribol. 2011, 133, 021704. [Google Scholar] [CrossRef]
  5. Ma, X.; Meng, X.; Wang, Y.; Peng, X. Suction effect of cavitation in the reverse-spiral-grooved mechanical face seals. Tribol. Int. 2019, 132, 142–153. [Google Scholar] [CrossRef]
  6. Lin, Q.; Bao, Q.; Li, K.; Khonsari, M.M.; Zhao, H. An investigation into the transient behavior of journal bearing with surface texture based on fluid-structure interaction approach. Tribol. Int. 2018, 118, 246–255. [Google Scholar] [CrossRef]
  7. Han, J.; Fang, L.; Sun, J.P.; Ge, S.R. Hydrodynamic lubrication of microdimple textured surface using three-dimensional CFD. Tribol. Trans. 2010, 53, 860–870. [Google Scholar] [CrossRef]
  8. Han, J.; Fang, L.; Sun, J.P.; Wang, Y.Q.; Ge, S.R.; Zhu, H. Hydrodynamic lubrication of surfaces with asymmetric microdimple. Tribol. Trans. 2011, 54, 607–615. [Google Scholar] [CrossRef]
  9. Ma, X.; Meng, X.; Wang, Y.; Liang, Y.; Peng, X. Fluid inertia effect on spiral-grooved mechanical face seals considering cavitation effects. Tribol. Trans. 2021, 64, 367–380. [Google Scholar] [CrossRef]
  10. Ausas, R.; Ragot, P.; Leiva, J.; Jai, M.; Bayada, G.; Buscaglia, G.C. The impact of the cavitation model in the analysis of microtextured lubricated journal bearings. J. Tribol. 2007, 129, 868–875. [Google Scholar] [CrossRef]
  11. Qiu, Y.Y.; Khonsari, M.M. On the prediction of cavitation in dimples using a mass-conservative algorithm. J. Tribol. 2009, 131, 041702. [Google Scholar] [CrossRef]
  12. Zhang, J.Y.; Meng, Y.G. Direct observation of cavitation phenomenon and hydrodynamic lubrication analysis of textured surfaces. Tribol. Lett. 2012, 46, 147–158. [Google Scholar] [CrossRef]
  13. Singhal, A.K.; Athavale, M.M.; Li, H.; Jiang, Y. Mathematical basis and validation of the full cavitation model. J. Fluids Eng. 2002, 124, 617–624. [Google Scholar] [CrossRef]
  14. Bakir, F.; Rey, R.; Gerber, A.G.; Belamri, T.; Hutchinson, B. Numerical and experimental investigations of the cavitating behavior of an inducer. Int. J. Rotating Mach. 2004, 10, 15–25. [Google Scholar] [CrossRef] [Green Version]
  15. Plesset, M.S. The dynamics of cavitation bubbles. J. Appl. Mech. 1949, 16, 277–282. [Google Scholar] [CrossRef]
  16. Morris, N.J.; Shahmohamadi, H.; Rahmani, R.; Rahnejat, H.; Garner, C.P. Combined experimental and multiphase computational fluid dynamics analysis of surface textured journal bearings in mixed regime of lubrication. Lubr. Sci. 2018, 30, 161–173. [Google Scholar] [CrossRef] [Green Version]
  17. Shahmohamadi, H.; Rahmani, R.; Rahnejat, H.; Garner, C.P.; Dowson, D. Big end bearing losses with thermal cavitation flow under cylinder deactivation. Tribol. Lett. 2015, 57, 2. [Google Scholar] [CrossRef] [Green Version]
  18. Dhande, D.Y.; Pande, D.W. Multiphase flow analysis of hydrodynamic journal bearing using CFD coupled fluid structure interaction considering cavitation. J. King Saud Univ.-Eng. Sci. 2018, 30, 345–354. [Google Scholar] [CrossRef] [Green Version]
  19. Cupillard, S.; Glavatskih, S.; Cervantes, M.J. Computational fluid dynamics analysis of a journal bearing with surface texturing. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2008, 222, 97–107. [Google Scholar] [CrossRef]
  20. Mishra, S.; Choudhury, A.; Sahu, S. CFD investigation of influences of reverse textures on bearing surface of a journal bearing. J. Appl. Fluid Mech. 2014, 7, 395–399. [Google Scholar]
  21. Song, Y.; Gu, C.W. Development and validation of a three-dimensional computational fluid dynamics analysis for journal bearings considering cavitation and conjugate heat transfer. J. Eng. Gas Turbines Power 2015, 137, 122502. [Google Scholar] [CrossRef]
  22. Wang, T.; Huang, W.F.; Liu, Y.; Liu, X.F.; Wang, Y.M. A homogeneous phase change model for two-phase mechanical seals with three-dimensional face structures. J. Tribol. 2014, 136, 041708. [Google Scholar] [CrossRef]
  23. Constantinescu, V.N.; Galetuse, S. Operating characteristics of journal bearings in turbulent inertial flow. J. Lubr. Technol. 1982, 104, 173–179. [Google Scholar] [CrossRef]
  24. Brunetière, N.; Tournerie, B. Finite element solution of inertia influenced flow in thin fluid films. J. Tribol. 2007, 129, 876–886. [Google Scholar] [CrossRef]
  25. Wang, T.; Huang, W.F.; Liu, X.F.; Li, Y.J.; Wang, Y.M. Experimental study of two-phase mechanical face seals with laser surface texturing. Tribol. Int. 2014, 72, 90–97. [Google Scholar] [CrossRef]
  26. Li, W.; Stephens, L.S.; Wenk, J.F. Experimental benchmarking of the numerical model of a radial lip seal with a surface textured shaft. Tribol. Trans. 2013, 56, 75–87. [Google Scholar] [CrossRef]
  27. Rao, T.V.V.L.N.; Rani, A.M.A.; Nagarajan, T.; Hashim, F.M. Analysis of slider and journal bearing using partially textured slip surface. Tribol. Int. 2012, 56, 121–128. [Google Scholar] [CrossRef]
  28. Cross, A.T.; Sadeghi, F.; Rateick, G.R.; Rowan, S.; Laboda, D. Temperature distribution in pocketed thrust washers. Tribol. Trans. 2015, 58, 31–43. [Google Scholar] [CrossRef]
  29. Tomanik, E. Modelling the hydrodynamic support of cylinder bore and piston rings with laser textured surfaces. Tribol. Int. 2013, 59, 90–96. [Google Scholar] [CrossRef]
  30. Morgut, M.; Nobile, E.; Bilus, I. Comparison of mass transfer models for the numerical prediction of sheet cavitation around a hydrofoil. Int. J. Multiph. Flow 2011, 37, 620–626. [Google Scholar] [CrossRef]
  31. Jakobsson, B.; Floberg, L. The Finite Journal Bearing, Considering Vaporization; Chalmers Tekniska Högskolas Handlingar: Gothenburg, Sweden, 1957. [Google Scholar]
  32. Meng, X.K.; Bai, S.X.; Peng, X.D. Lubrication film flow control by oriented dimples for liquid lubricated mechanical seals. Tribol. Int. 2014, 77, 132–141. [Google Scholar] [CrossRef]
  33. Etsion, I.; Ludwig, L.P. Observation of pressure variation in the cavitation region of submerged journal bearings. ASME J. Lubr. Tech. 1982, 104, 157–163. [Google Scholar] [CrossRef]
  34. Sun, D.C.; Brewe, D.E. A high-speed photography study of cavitation in a dynamically loaded journal bearing. ASME J. Tribol. 1991, 113, 287–292. [Google Scholar] [CrossRef]
  35. Sun, D.C.; Brewe, D.E.; Abel, P.B. Simultaneous pressure measurement and high-speed photography study of cavitation in a dynamically loaded journal bearing. ASME J. Tribol. 1993, 115, 88–95. [Google Scholar] [CrossRef]
Figure 1. Geometrical model: (a) Parallel sliding surfaces, (b) Textured surface cell and (c) Recess shapes.
Figure 1. Geometrical model: (a) Parallel sliding surfaces, (b) Textured surface cell and (c) Recess shapes.
Lubricants 11 00270 g001
Figure 2. Pressure distribution of the journal bearing [31].
Figure 2. Pressure distribution of the journal bearing [31].
Lubricants 11 00270 g002
Figure 3. Pressure distribution of the circle recess.
Figure 3. Pressure distribution of the circle recess.
Lubricants 11 00270 g003
Figure 4. Pressure distribution of the square recess.
Figure 4. Pressure distribution of the square recess.
Lubricants 11 00270 g004
Figure 5. (a) Mesh and (b) Validation of mesh independence.
Figure 5. (a) Mesh and (b) Validation of mesh independence.
Lubricants 11 00270 g005
Figure 6. Pressure distribution comparison with and without cavitation.
Figure 6. Pressure distribution comparison with and without cavitation.
Lubricants 11 00270 g006
Figure 7. Effect of cavitation and Re on pressure distribution.
Figure 7. Effect of cavitation and Re on pressure distribution.
Lubricants 11 00270 g007
Figure 8. Effect of inertia and Re on pressure distribution.
Figure 8. Effect of inertia and Re on pressure distribution.
Lubricants 11 00270 g008
Figure 9. Effect of Re on Fo.
Figure 9. Effect of Re on Fo.
Lubricants 11 00270 g009
Figure 10. Effect of pc on pressure distribution.
Figure 10. Effect of pc on pressure distribution.
Lubricants 11 00270 g010
Figure 11. Effect of pc on Fo.
Figure 11. Effect of pc on Fo.
Lubricants 11 00270 g011
Figure 12. Effect of Sr.
Figure 12. Effect of Sr.
Lubricants 11 00270 g012
Figure 13. Effect of h1.
Figure 13. Effect of h1.
Lubricants 11 00270 g013
Figure 14. Effect of recess shapes on pressure distribution.
Figure 14. Effect of recess shapes on pressure distribution.
Lubricants 11 00270 g014
Table 1. Geometric parameters and operating parameters.
Table 1. Geometric parameters and operating parameters.
ParametersValues
Textured surface cell side length, W0/(mm)6
Recess size length, W1/(mm)4
Gap thickness, h0/(μm)50
Recess depth, h1/(μm)4
Ambient pressure, pa/(MPa)0.1
Cavitation pressure, pc/(MPa)0.095
Dynamic viscosity of liquid, μL/(Pa·s)0.0127
Liquid density, ρL/(kg/m3)840
Reynolds number, Re165.35
Table 2. Effect of recess shapes on peak pressure and LCC.
Table 2. Effect of recess shapes on peak pressure and LCC.
SquareCircleTriangle
p(pmax ICpmax C)/ pmax C 39.42%38.88%39.10%
(pmax Ipmax RE)/ pmax RE8.38%8.94%11.11%
LCCIC7.3%7%4.7%
I0.6%0.43%1.4%
C≈00.4%≈0
RE≈0≈0≈0
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

Ma, X. Combined Effect of Fluid Cavitation and Inertia on the Pressure Buildup of Parallel Textured Surfaces. Lubricants 2023, 11, 270. https://doi.org/10.3390/lubricants11070270

AMA Style

Ma X. Combined Effect of Fluid Cavitation and Inertia on the Pressure Buildup of Parallel Textured Surfaces. Lubricants. 2023; 11(7):270. https://doi.org/10.3390/lubricants11070270

Chicago/Turabian Style

Ma, Xuezhong. 2023. "Combined Effect of Fluid Cavitation and Inertia on the Pressure Buildup of Parallel Textured Surfaces" Lubricants 11, no. 7: 270. https://doi.org/10.3390/lubricants11070270

APA Style

Ma, X. (2023). Combined Effect of Fluid Cavitation and Inertia on the Pressure Buildup of Parallel Textured Surfaces. Lubricants, 11(7), 270. https://doi.org/10.3390/lubricants11070270

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