Next Article in Journal
Effect of Feedstock Powder Intrinsic Characteristics on the Tribological Behavior of Inductively Remelted NiCrBSi Flame-Sprayed Coatings
Next Article in Special Issue
Partial Surface Texturing in Hydrodynamic Lubrication: A CFD-Based Investigation
Previous Article in Journal
Skidding Analysis of Exhaust Cam-Roller Unit in the Steady/Startup Operation of Internal Combustion Engine
Previous Article in Special Issue
Research on Internal Flow Field Characteristics of Straight-Groove Texture Using Three-Dimensional Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of the Turbulent Lubrication of a Textured Hydrodynamic Journal Bearing

1
School of Mechanical Engineering, Henan University of Engineering, Zhengzhou 451191, China
2
College of Vehicle and Traffic Engineering, Henan University of Science and Technology, Luoyang 471003, China
*
Author to whom correspondence should be addressed.
Lubricants 2023, 11(9), 362; https://doi.org/10.3390/lubricants11090362
Submission received: 14 July 2023 / Revised: 9 August 2023 / Accepted: 21 August 2023 / Published: 25 August 2023
(This article belongs to the Special Issue Hydrodynamic Lubrication of Textured Surfaces)

Abstract

:
In order to investigate the turbulent lubrication performance of a textured hydrodynamic journal bearing (THJB), a model of turbulent lubrication was established in this paper. The variations in the Reynolds number, oil film thickness, oil film pressure, bearing capacity, attitude angle, and side leakage flow with structural and working parameters were studied, and the axis whirl orbit was further analyzed. The results show that turbulent lubrication is suitable for the actual operating conditions of THJBs. The Reynolds number decreases with the eccentricity ratio in the pressure-bearing zone but increases with rotational speeds, whereas the variation in the maximum oil film pressure increases and the minimum oil film thickness decreases with the eccentricity ratio under various Reynolds numbers. The bearing capacity decreases with the dimple diameter, depth, oil film thickness, and clearance ratio but increases with the length/diameter ratio and dimple spacing. As the eccentricity ratio increases, the attitude angle decreases, but the side leakage flow increases. In addition, the system tends to be unstable as the rotational speed and length/diameter ratio increase, and the friction and wear on the surface are three-body friction. This work not only helps in analyzing the characteristics of a THJB under actual operating conditions but also provides support for research on the simulation of THJB’s lubrication mechanism of THJB via computational fluid dynamics.

1. Introduction

With the development of highly precise rotational machinery that can achieve high speed and tolerate heavy loads [1,2,3], as a key part that supports the rotating spindle, the lubrication performance of a hydrodynamic journal bearing (HJB) directly affects the system of rotating machinery. The surface texture can effectively make the lubrication performance of HJBs better [2,4,5], and the lubrication performance of textured bearings will be affected by the flow regime; laminar and turbulence are the common flow regimes that affect the properties of lubrication.
Until now, the development of laminar lubrication theory has been quite comprehensive. Turbulence is often encountered in natural science and engineering technology, but laminar lubrication, which has been studied well, is seldom encountered. In 1923, the viscous flow between concentric cylinders was studied, an expression for calculating the critical Reynolds number between concentric cylinders was proposed [6], and the application of turbulent lubrication theory in the design of friction pairs was deeply studied [7,8]. Since 1923, the lubrication theory has made a breakthrough. In particular, a turbulent Reynolds equation was established based on the Ng-Pan theory [9], and the lubrication performance parameters under a greater range of Reynolds numbers were obtained. The turbulent model was established to study the performance of turbulent lubrication, and the results are in good agreement with the corresponding experimental data [10,11,12]. While the finite difference method was used to solve a modified Reynolds equation of couple stress fluids considering turbulence [13,14,15], the dynamic properties of water-lubricated bearings were analyzed under contact conditions, and it was found that the dynamic characteristics were significantly affected by speed and pressure [16]. Jiang’s team conducted a series of studies, and the results showed that the influence of factors such as film cavitation, inertial force, and turbulence on the bearing lubrication performance could not be ignored under high-speed working conditions [17,18,19] and that surface topography has an impact on bearing performance [20,21]. Moreover, the lubrication performances of bearings with groove angles of 36° and 18° in a turbulent regime were analyzed, and the load capacity improved as the misalignment angle increased [22]. Further studies have shown that local turbulence has a significant impact on the friction coefficient [23,24], and the influence of micro-groove on lubrication characteristics and bearing capacity varies with variations of the eccentricity ratio [25,26]. The effect of the micro-groove on the lubrication performance of water-lubricated bearings during start-up was also studied [27,28,29].
A series of reports on lubrication properties have been analyzed by scholars, but there are few reports on the performance of micro pit THJBs under turbulent regimes. In this work, a systematic study of the influences of geometry and working condition parameters on the characteristics of micro pit THJBs was carried out by the authors. This research has a certain significance for our understanding of the performance laws and engineering applications of THJBs serving in turbulent regimes.

2. Geometry Model and Governing Equations

2.1. Geometry Model

As schematic and section of a THJB are shown in Figure 1. In Figure 1, O1 and O2 represent the center of THJB and the rotating shaft, respectively, e is the eccentric distance, R and r represent the bearing radius and spindle radius, respectively, c is the clearance c = Rr, L is the bearing width, and hmin and hmax are the minimum and maximum oil film thickness, respectively.
The geometric parameters selected for the investigation of the lubrication characteristics of THJBs within a turbulent model are shown in Table 1.

2.2. Flow Regime Criterion

The flow regime of a lubricating oil regime in the clearance of a bearing is related to its Reynolds number. The Reynolds numbers ReA and ReB of the THJB in region A and region B can be calculated using the following expression [6]:
Re A = ρ × u × h / μ
Re B = ρ × u × h × ( 1 - a × ε ) / μ
where ρ is the density, u is the shaft velocity, h is the oil film thickness, μ is the viscosity, and a is the coefficient relating to the bearing’s structural and working conditions, which can be obtained according to experimental data.
Expressions (1) and (2) show that the relationship between the Reynolds number, velocity, and oil film thickness, and the expression of the critical Reynolds number under turbulence can be obtained via estimation as follows [30]:
Re cc = 41.1 1 + 1.5 ε 2 / ψ
where ε is the eccentricity ratio, ε = e/c, ψ is the clearance ratio, and ψ = c/R. The critical Reynolds number Recc is obtained according to the expression in (3).

2.3. Reynolds Governing Equation

According to the flow regime analysis of the lubricant mentioned above, the Reynolds governing equation is as follows:
x h 3 G θ μ p x + z h 3 G λ μ p z = 6 R ω h x + 12 h t
where x and z represent two coordinates of a cylindrical coordinate system, ω represents the shaft angular velocity, p represents the oil film pressure, and Gθ and Gλ represent turbulence factors.
Considering the convenience of the analysis process and the simplicity of the governing equation, the dimensionless, unified Reynolds governing equation, which is suitable for analyzing the turbulent lubrication of THJBs, is provided below:
θ H 3 G θ P θ + D L 2 λ H 3 G λ P λ = H θ + 2 e ˙ cos θ
where H and P represent the normalized oil film thickness and pressure, respectively, D and L represent the diameter and length of the bearing, respectively, and λ and θ represent the normalized circumferential coordinates and axial coordinates, respectively.
In the dimensionless process of the Reynolds equation [4], the following dimensionless parameters are introduced:
x = R θ h = Hc z = λ L / 2 p = Pp 0 p 0 = 6 μ ω / ψ 2
Based on the Ng-Pan mathematical model of lubrication [10], the turbulence factors Gθ and Gλ in Equation (4) are Gθ = 1/(12 + 0.0136Recc0.9) and Gλ = 1/(12 + 0.0043Recc0.98), respectively, while the oil film thickness of the THJB mainly consists of untextured and textured oil film thicknesses, and its expression for the THJB is shown as follows:
H = 1 + ε cos θ φ + i = 1 n H i
where Hi represents the oil film thickness of the micro-pits and i represents the number of micro-dimples.
The Reynolds boundary condition of the Reynolds governing equation was adopted. In this work, the oil film pressure of the first column and column (m + 1) were set to zero as a starting boundary. The process of realizing the oil film’s pressure rupture boundary was as follows: The pressure along the circumferential direction is calculated point by point from j = 2~n rows and i = 2~m columns (see Figure 2). If the pressure value is negative, it is determined that oil film rupture occurs at this point, and the pressure values at this point and the point after its row are set to zero. After an iterative calculation, the rupture boundary gradually approaches the natural rupture boundary, and the pressure distribution meets the Reynolds boundary condition.

2.4. Bearing Capacity and Attitude Angle Equation

As the oil film force comes into balance with the normal force, the oil film component forces of the THJB are Wx and Wy, respectively. Therefore, the expression of the bearing’s capacity is as follows:
W x = 0 2 π 0 L P cos θ λ d θ d λ W y = 0 2 π 0 L P sin θ λ d θ d λ W = W x 2 + W y 2
Thus, the attitude angle of the THJB is shown as follows:
φ = arctg W y W x

2.5. Lubricant Side Leakage Flow Equation

For THJBs, the velocity components in the circumferential direction and the axial direction are as follows:
V x = 1 2 μ p G θ x y 2 hy + U h y V z = 1 2 μ p G λ z y 2 hy
According to the flow equation, Equation (10) is integrated, and the flow on any given section is obtained as follows:
q x = h 3 p 12 μ G θ x + Uh 2 q z = - h 3 p 12 μ G λ z
Thus, the dimensionless expression of the obtained lubricant side leakage flow is shown as follows:
Q = RL 0 2 π H 3 c 3 μ p 0 P G λ ( λ L ) λ = 1 d φ 6 r 2 cU = - 0 2 π H 3 0 2 π P G λ λ λ = 1 d φ

2.6. Axis Whirl Orbit Governing Equation

Axis whirl orbit analysis is an important way to determine the stability of the rotor. The governing equation of the axis whirl orbit for a rotor–THJB system is shown as follows.
m x ¨ = W x + m ω n 2 sin ω t m y ¨ = W y + me ω n 2 cos ω t + mg
where e is the centrifugal distance, ω is the speed, Wx and Wy are the oil film forces in x-axis and y-axis directions, respectively, g is the acceleration of gravity, and x and y are the displacements of journal in x axis and y axis directions, respectively.
The expressions Wx and Wy of the oil film force components on the journal could be obtained by solving the Reynolds governing equation, and the instantaneous acceleration of the journal was obtained via the following expression.
x ¨ = W x / m y ¨ = W y W / m
Then, the expressions for the transient velocity and displacement of the axis were
x ˙ i = x ¨ i t + x ˙ i 1 y ˙ i = y ¨ i t + y ˙ i 1 x i = x ˙ i t + x i 1 y i = y ˙ i t + y i 1
where Δt is a tiny time interval, Δt = 120/50n, and n is the rotational speed. The axis whirl orbit of the journal at any time was obtained via Expressions (13)–(15).

3. Numerical Computation Method and Effectiveness Analysis

In this work, the effectiveness was verified prior to conducting a numerical analysis. In order to verify the correctness of the finite difference solution, the turbulent lubrication performance of a THJB with the geometric parameters ψ = 0.003 and L/D = 1 was analyzed, and the Sommerfeld number, a characteristic parameter, was obtained via calculation and compared with results in the literature [31]. The Sommerfeld number is:
S = μ Lu π W R c 2
The results of comparing the solution obtained via FDM with test values are shown in Figure 3. In Figure 3, the calculations agree well with test values, which proves the correctness of the FDM applied in this work.
Equation (5) was solved by FDM, and the grid number was important when ensuring the correctness of the calculation result. In order to identify a suitable differential grid, the properties of a THJB were analyzed using various differential grids based on the geometric parameters in Table 1. The results are listed in Table 2. Relative errors of oil film pressure are 0.54% and 1.48%, respectively, for the 390 × 125 differential grid compared with the 628 × 200 differential grid, while the relative errors are 0.03% and 0.15%, respectively, in [32,33] compared with the FDM result achieved using the same differential grid, 314 × 100. Therefore, the 314 × 100 differential grid is used in this paper.
In order to further verify the accuracy of calculation results based on the turbulent lubrication model, a comparison of the theoretical result, test value result, and the results from other studies in the literature [34,35] for the oil film pressure distribution are shown in Figure 4. The results show that the difference between the theoretical result presented in this paper, the test value, and the result from other studies are relatively small, which proves the accuracy of the calculation result based on the turbulent model.

4. Results and Discussion

4.1. Lubrication Characteristics of a Textured Bearing under the Turbulent Model

4.1.1. Reynolds Number Varies with Eccentricity Ratios and Rotational Speeds

The variation in the Reynolds number with eccentricity ratios (ε = 0.25, 0.50, 0.75, and 0.95, n = 1000 rpm) and rotational speeds (ε = 0.25, n = 200, 400, 600, 800, 1000, 1200, and 1400 rpm) is shown in Figure 5. In Figure 5a, the Reynolds number decreases with an increase in the eccentricity ratios in the pressure-bearing region. In contrast, the Reynolds number increases with eccentricity ratios, and the Reynolds number for the pressure-bearing region is higher than the critical Reynolds number; however, the opposite is true for non-pressure-bearing regions. This is because the oil film thickness in the pressure-bearing zone is relatively thinner, and the mainstream region where the Reynolds number is lower than the critical Reynolds number is laminar. Conversely, the mainstream region of the lubricating oil is turbulent where the flow states are more complex. In Figure 5b, the Reynolds number increases with the increase in the rotational speed n, and the actual Reynolds number is less than the critical Reynolds number when n = 200–800 rpm. When n = 1000–1200 rpm, the actual Reynolds number is partially less than the critical Reynolds number, and the others are higher than the critical Reynolds number. However, when n is 1400 rpm, the Reynolds number is higher than the critical Reynolds number. The reason for this is that the flow regimes of lubricating oil are laminar and the mixture flows when n = 200–800 rpm and 1000–1200 rpm, respectively. In addition, the flow regime of lubricating oil is turbulent when n is 1400 rpm.

4.1.2. Variations in Oil Film Thickness and Maximum Oil Film Pressure

Variations in oil film thickness and pressure with eccentricity ratio and angular coordinate are shown in Figure 6. In Figure 6a, the trend in the variation in the oil film thickness is basically consistent and shows a sine wave shape under six eccentricity ratios. As the eccentricity ratio increases, the minimum oil film thickness decreases gradually. The reason for this is that as the eccentricity ratio increases, the load on the rotor increases and the offset of the axial center position increases, which causes the oil film thickness to decrease. In Figure 6b, the minimum oil film thickness decreases with the increase in the eccentricity ratio and follows the rule of H = −ε, while the maximum oil film pressure increases gradually with eccentricity ratios. In addition, the minimum oil film thickness is not affected by the variation in the Reynolds numbers under the same eccentricity, and the maximum oil film pressure is almost unaffected by the eccentricity ratios (<0.25) and Reynolds numbers. The larger the Reynolds number at the same eccentricity ratio, the greater the corresponding maximum oil film pressure.

4.1.3. Variation in Oil Film Pressure Distribution

The oil film distribution of a THJB with a length-to-diameter ratio L/D = 1 and an eccentricity ratio ε = 0.3 is shown in Figure 7. In Figure 7a, the circumferential oil film pressure distribution is basically the same, but the oil film pressure under a turbulent regime is much higher than under the laminar regime, which indicates that a turbulent regime can effectively make the oil film pressure better. In order to observe the pressure distribution of the THJB under a turbulent regime, the variation in the oil film pressure is shown in Figure 7b. In Figure 7b, the smallest and largest THJB oil film pressures are located in the leaking position and the axial center position X = 0, respectively. Moreover, the oil film pressure is symmetrically distributed along the axial direction, with X = 0 as the axis of symmetry.
The circumferential and axial oil film pressure distributions of the THJB are shown in Figure 8. In Figure 8a, the variation in the oil film pressure distribution becomes serrated due to the existence of surface texture, and arbitrary serrations correspond to the texture point. The dynamic lubrication of the THJB system is superimposed, which further enhances the oil film pressure as the surface texture is in the convergence region. Therefore, the serrated change is not obvious. In contrast, in the divergent region, the negative pressure effect of the surface texture will truncate the dynamic pressure action in the convergence region, which will cause the cavitation phenomenon to be more obvious and the dynamic pressure lubrication effect caused by the gap contraction to be more prominent so that the serration presents a significant change. In Figure 8b, the axial oil film pressure distribution shows a continuous phenomenon of alternating changes between high-pressure and low-pressure regions. The maximum and minimum oil film pressure values occur at the axial center coordinate (X = 0) and the axial leaking position (X = ±0.5), respectively.
In order to further explain the variation in the pressure distribution of THJB, the local oil film pressure distribution is shown in Figure 9. In Figure 9a, there is a zone in which the relative oil film pressure increases around the drop zone of the THJB under a turbulent regime. This is because the distribution interval of surface texture will lead to the local film thickness increasing as shown in Figure 9b. According to the hydrodynamic lubrication mechanism, an increase in the oil film thickness in the local area will inevitably lead to a decrease in the oil film pressure; thus, the situation in which the oil film pressure is reduced is shown in Figure 9a. However, as the lubricating oil in the micro-pit texture flows out, the clearance becomes smaller and overlaps with the hydrodynamic effects of the THJB. Compared with the area in which the inner oil film pressure of the micro pit decreases, the outside has a higher oil film pressure increase area, that is, the low-pressure area appears in front of the surface texture the and high-pressure area appears in the rear.

4.2. The Influence of Structural and Working Condition Parameters

The bearing capacity distribution of a THJB is influenced by the micro dimple diameter, depth, spacing, and oil film thickness, as shown in Figure 10. In Figure 10a,b, the bearing capacity increases gradually as the Reynolds number increases with the same dimple diameter and depth. Moreover, a higher Reynolds number is beneficial to improve the bearing capacity. Therefore, the Reynolds number is a very important factor affecting the bearing capacity.
In Figure 10c, the bearing capacity of the THJB increases gradually as the dimple spacing increases. The reason for this is that the adjacent micro pits do not interfere with each other when the spacing is small. However, the adjacent micro pits experience a mutual disturbance with an increase in the spacing, which causes the change rule to increase slowly. Additionally, the bearing capacity increases gradually as the Reynolds number increases with the same dimple spacing. Thus, a larger load capacity is obtained by changing the dimple spacing and Reynolds number.
Moreover, in Figure 10d, there is an approximately negative linear relationship between the oil film thickness and bearing capacity. The bearing capacity is approximately twice that of 1.44 as the oil film thickness is 1.26. Additionally, the smaller the eccentricity ratio, Reynolds number, and length/diameter ratio, the smaller the bearing capacity. However, the maximum bearing capacity of the THJB occurs at a length-to-diameter ratio of L/D = 1, an eccentricity ratio ε = 0.5, and a Reynolds number of Recc = 800. Therefore, the load capacity will be affected by the geometric parameter and Reynolds number.
The variation in the bearing capacity of the THJB with the length/diameter ratio and clearance ratio is shown in Figure 11. In Figure 11a, the bearing capacity increases as the length/diameter ratio increases. If the length/diameter ratio is too large, its influence on the bearing capacity is almost unchanged. The reason for this is that with the increase in the length/diameter ratio, although the axial width increases and the bearing range increases, the pressure per unit area decreases, so the bearing capacity decreases as the length-to-diameter ratio increases. Therefore, the length/diameter ratios of THJBs should be selected in the range of 1–2.4. The bearing capacity increases with an increase in the eccentricity ratio under the same Reynolds number and the Reynolds number under the same eccentricity ratio. In addition, the bearing capacity of a THJB with a Reynolds number of Recc0.3 = 1500 is better than that of a textured bearing with Reynolds numbers of Recc0.25 = 1000 and Recc0.25 = 1200. However, in Figure 11b, the bearing capacity decreases with an increase in the clearance ratio but increases gradually as the Reynolds number increases under the same clearance ratio, and a higher Reynolds number is beneficial to improve the bearing capacity.
The distribution of the bearing capacity of a THJB with rotational speed is shown in Figure 12. In Figure 12a, the linear increase in the bearing capacity with an increasing rotational speed occurs because the increase in rotational speed changes the THJB’s oil film pressure distribution. However, at lower rotational speeds and at the same rotational speed, the bearing capacity of the THJB under a turbulent regime is approximately the same as that of a THJB under a laminar regime because, at this speed, the flow regime gradually changes from laminar to turbulent. In Figure 12b, the bearing capacity becomes greater as the length/diameter ratio increases with the increase in the rotational speed under the same Reynolds number and eccentricity ratio. However, under the same length/diameter ratio and eccentricity ratio, the bearing capacity increases as the Reynolds number increases with the increase in the rotational speed, for example, Recc = 600–1200. In addition, the greater eccentricity, the greater the bearing capacity as the rotational speed increases under the same Reynolds number and length/diameter ratio.

4.3. Variations in the Attitude Angle and Lubricant Side Leakage Flow with the Eccentricity Ratio

The variations in the attitude angle and lubricant side leakage flow with the eccentricity ratio are shown in Figure 13. In Figure 13a, the attitude angle decreases as the eccentricity ratio increases, which means that bearing capacity increases as the eccentricity ratio increases, and the oil film’s bearing capacity increases as the attitude angle decreases. In addition, when extreme conditions occur, since the eccentricity ratio ε tends toward 1, the attitude angle φ tends toward 0°.This is because as eccentricity ε = 1, the shaft and bearing bush are in direct contact, the direction of the bearing capacity is along the connecting line between the shaft diameter and bearing bush, and the corresponding attitude angle is 0°.
In Figure 13b, the dimensionless lubricant’s side leakage flow increases with the increase in the eccentricity ratio. This occurs because the lubricating oil’s side leakage flow is determined by the oil film pressure gradient in the axial direction of the THJB. As the eccentricity ratio increases, the oil film’s pressure gradient also increases at both ends of the THJB. Therefore, the dimensionless lubricating oil’s side leakage flow increases as the eccentricity ratio increases, whereas when the eccentricity ratio tends toward 1 due to dry friction occurring between the shaft diameter and bearing bush, and when good hydrodynamic lubrication cannot be established, no oil film pressure gradient will be generated. Thus, the dimensionless lubricating oil’s side leakage flow tends toward zero at this time. Additionally, as the length/diameter ratio increases, the THJB’s lubricating oil’s side leakage flow decreases. This is because the axial flow resistance increases as the length/diameter ratio increases, which leads to decreases in the pressure gradient.

4.4. The Axis Whirl Orbit Variation of the Rotor

The axis whirl orbit variation of the rotor calculated based on Section 2.6 is shown in Figure 14. In Figure 14a,b, the rotor is initially located at the center of the THJB, and offset has not occurred due to the rotational spindle. The rotor begins to move downward due to the action of gravity, and the pressure generated by the oil film gradually increases as the eccentricity ratio increases. When the pressure generated by the oil film is higher than the rotor’s gravity, the rotor starts to move upward, and the running area gradually converges and tends toward a stable state. In Figure 14a,d, with the increase in the length/diameter ratio, the axis whirl orbit of the rotor is a limit cycle, and the rotor is in a critical state at this time. Additionally, in Figure 14a,c, as the speed increases, the divergent state is occupied by the axis whirl orbit. Compared with the nonlinear axis whirl orbit when the rotational speed is 1000 rpm, the rotor axis whirl orbit diverges significantly outward. The reason for this is that the oil film flow regime tends to become more complex when the rotational speed increases, the oil film force on the rotor increases, and the attitude angle of the bearing capacity is larger, so the axis whirl orbit of the rotor is in a divergent state.
Figure 15 shows the friction and wear surface morphology of the THJB when the load and rotational speed were 4 kN, 1100 rpm, and 2000 rpm, respectively. In Figure 15, obvious scratches can be seen on the surface of the THJB. Apart from the micro pits being smoothed out, most of the micro pit textures were filled with wear particles. Because the rotor-bearing system was working at high speeds and under a heavy load, the lubricating oil was not supplied in time, and the particles produced by friction were not removed effectively, which caused obvious grinding and filling of furrows.

5. Conclusions

In this present work, a model of the turbulent lubrication properties of a THJB was investigated. A few conclusions are summarized as follows:
(1)
A turbulent lubrication model was established for a performance analysis based on the Ng-Pan model, and the correctness of the model and the numerical method was verified. For the THJB in this work, the flow regimes were laminar and mixed flow at rotational speeds of n = 200~800 rpm and 1000~1200 rpm, whereas a completely turbulent regime occurred at n ≥ 1400 rpm. The difference between the values calculated for the laminar and turbulent model shows that the turbulent effect enhanced the dynamic pressure effect of lubricating oil, which cannot be ignored when designing textured bearings.
(2)
The Reynolds number decreases with the eccentricity in the pressure-bearing zone but increases with the rotational speed. The variation in the maximum oil film pressure increases and the minimum oil film thickness decreases with the eccentricity ratio under various Reynolds numbers. Furthermore, the bearing capacity of the THJB under a turbulent regime decreases with the dimple diameter, depth, oil film thickness, and clearance ratio but increases with the dimple spacing and length/diameter ratio. In addition, the attitude angle decreases but the side leakage flow increases as the eccentricity ratio increases.
(3)
The axis trajectory converged to an equilibrium point as the eccentricity values were ε = 0.25 and 0.95 under the same rotational speed n = 1000 rpm. Additionally, the axis trajectory spread outward and gradually became unstable as the rotational speed increased under the same eccentricity ε = 0.25. Further study showed that the mechanism of friction and wear during the axis trajectory operation was mainly three-body friction. Moreover, the axis trajectory approached the cycle limit gradually as the length/diameter ratio changed from 1 to 2, and the length/diameter ratio should be properly selected to ensure the stable operation of the rotor when designing the THJB.

Author Contributions

Conceptualization, Y.M. and L.L.; methodology, Y.M., L.L. and D.L.; software, J.Z. and D.L.; formal analysis, L.L., D.L. and J.Z.; investigation, L.L. and J.Z.; resources, D.L. and J.Z.; data curation, Y.M., L.L. and D.L.; writing—original draft preparation, Y.M. and L.L.; writing—review and editing, Y.M. and D.L.; supervision, L.L.; project administration, Y.M., D.L. and J.Z.; funding acquisition, D.L. and J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Henan Provincial Department of Science and Technology Research Project (No. 232102220041) and the Young Elite Scientists Sponsorship Program (No. 2023HYTP012).

Data Availability Statement

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

Acknowledgments

The authors would also like to express their sincere for their constructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mao, Y.; Yang, J.; Ji, J.; Xu, W.; Guo, Q. An analytical solution of Reynolds equation for evaluation the characteristics of surface textured bearing. Ind. Lubr. Tribol. 2020, 72, 1075–1085. [Google Scholar] [CrossRef]
  2. Mao, Y.; Yang, J.; Xu, W.; Liu, Y. Study on the influence of round pits arrangement patterns on tribological properties of journal bearings. Ind. Lubr. Tribol. 2019, 71, 931–941. [Google Scholar] [CrossRef]
  3. Lin, C.-G.; Yang, Y.-N.; Chu, J.-L.; Sima, C.; Liu, P.; Qi, L.-B.; Zou, M.-S. Study on nonlinear dynamic characteristics of propulsion shafting under friction contact of stern bearings. Tribol. Int. 2023, 183, 108391. [Google Scholar] [CrossRef]
  4. Grützmacher, P.G.; Profito, F.J.; Rosenkranz, A. Multi-scale surface texturing in tribology-Current knowledge and future perspectives. Lubricants 2019, 7, 95. [Google Scholar] [CrossRef]
  5. Rosenkranz, A.; Costa, H.L.; Profito, F.; Gachot, C.; Medina, S.; Dini, D. Influence of surface texturing on hydrodynamic friction in plane converging bearings-An experimental and numerical approach. Tribol. Int. 2019, 134, 190–204. [Google Scholar] [CrossRef]
  6. Taylor, G.I. Stability of a viscous liquid contained between two rotating cylinders. Proc. R. Soc. A 1923, 102, 541–542. [Google Scholar]
  7. Mariot, A.; Arghir, M.; Hélies, P. Experimental analysis of floating ring annular seals and comparisons with theoretical predictions. J. Eng. Gas Turbines Power 2015, 138, 42503. [Google Scholar] [CrossRef]
  8. Li, G.Q.; Zhang, Q.; Huang, E.L. Leakage performance of floating ring seal in cold/hot state for aero-engine. Chin. J. Aeronaut. 2019, 32, 2085–2094. [Google Scholar] [CrossRef]
  9. Das, S.; Guha, S.K. Non-linear stability analysis of micropolar fluid lubricated journal bearings with turbulent effect. Ind. Lubr. Tribol. 2019, 71, 31–39. [Google Scholar] [CrossRef]
  10. Zhu, S.; Sun, J.; Li, B.; Zhao, X.; Wang, H.; Teng, Q.; Ren, Y.; Zhu, G. Stochastic models for turbulent lubrication of bearing with rough surfaces. Tribol. Int. 2019, 136, 224–233. [Google Scholar] [CrossRef]
  11. Zhu, S.; Sun, J.; Li, B.; Zhu, G. Thermal turbulent lubrication analysis of rough surface journal bearing with journal misalignment. Tribol. Int. 2019, 144, 106109. [Google Scholar] [CrossRef]
  12. Du, Y.; Lan, J.; Quan, H.; Sun, C.; Liu, X.; Yang, X. Effect of different turbulent lubrication models on the lubrication characteristics of water-lubricated rubber bearings at a high Reynolds number. Phys. Fluids 2021, 33, 065118. [Google Scholar] [CrossRef]
  13. Chetti, B. The effect of turbulence on the performance of a two-lobe journal bearing lubricated with a couple stress fluid. Ind. Lubr. Tribol. 2016, 68, 336–340. [Google Scholar] [CrossRef]
  14. Cai, J.; Xiang, G.; Li, S.; Guo, J.; Wang, J.; Chen, S.; Yang, T. Mathematical modeling for nonlinear dynamic mixed friction behaviors of novel coupled bearing lubricated with low-viscosity fluid. Phys. Fluids 2022, 34, 093612. [Google Scholar] [CrossRef]
  15. Singh, A.; Sharma, S.C. Influence of geometric imperfection of journal on the performance of porous hybrid journal bearing under turbulent condition. Proc. Inst. Mech. Eng. Part J 2021, 235, 1879–1896. [Google Scholar] [CrossRef]
  16. Ouyang, W.; Zhang, X.; Jin, Y.; Yuan, X. Experimental study on the dynamic performance of water lubricated rubber bearings with local contact. Shock Vib. 2018, 2018, 6309727. [Google Scholar] [CrossRef]
  17. Feng, H.H.; Jiang, S.Y.; Ji, A.M. Investigations of the static and dynamic characteristics of water-lubricated hydrodynamic journal bearing considering turbulent, thermohydrodynamic and misaligned effects. Tribol. Int. 2019, 130, 245–260. [Google Scholar] [CrossRef]
  18. Lin, X.H.; Jiang, S.Y.; Zhang, C.B. Thermohydrodynamic analysis of high speed water-lubricated spiral groove thrust bearing considering effects of cavitation, inertia and turbulence. Tribol. Int. 2018, 119, 645–665. [Google Scholar] [CrossRef]
  19. Du, Y.Y.; Li, M. Effects on lubrication characteristics of water-lubricated rubber bearings with journal tilting and surface roughness. Proc. Inst. Mech. Eng. Part J 2019, 234, 161–171. [Google Scholar] [CrossRef]
  20. Zhang, X.H.; Xu, Y.; Jackson, R. A mixed lubrication analysis of a thrust bearing with fractal rough surfaces. Proc. Inst. Mech. Eng. Part J 2020, 234, 608–621. [Google Scholar] [CrossRef]
  21. Zhou, G.W.; Wu, K.P.; Pu, W.; Li, P.; Han, Y. Tribological modification of hydrogenated nitrile rubber nanocomposites for water-lubricated bearing of ship stern shaft. Wear 2022, 504, 204432. [Google Scholar] [CrossRef]
  22. Mallya, R.; Shenoy, S.B.; Pai, R. Static characteristics of misaligned multiple axial groove water-lubricated bearing in the turbulent regime. J. Eng. Tribol. 2017, 231, 385–398. [Google Scholar] [CrossRef]
  23. Lv, F.R.; Zou, D.L.; Ta, N.; Rao, Z.S. Influence of local turbulent flow on the performance of a mixed-lubrication bearing. Proc. Inst. Mech. Eng. Part J 2019, 233, 1029–1035. [Google Scholar] [CrossRef]
  24. Zhang, L.; Xu, H. Effect of radial clearance on the dynamic behavior of adjustable journal bearings in turbulent and laminar flow regimes. AIP Adv. 2023, 13, 025362. [Google Scholar] [CrossRef]
  25. Gong, J.Y.; Jin, Y.; Liu, Z.L. Study on influencing factors of lubrication performance of water-lubricated micro-groove bearing. Tribol. Int. 2019, 129, 390–397. [Google Scholar] [CrossRef]
  26. Yang, T.Y.; Cai, J.L.; Wang, L.W.; Tang, D.; Chen, S.; Wang, J. Numerical analysis of turbulence effect for coupled journal-thrust water-lubricated bearing with micro grooves. J. Tribol. 2023, 145, 084101. [Google Scholar] [CrossRef]
  27. Guo, X.; Han, Y.F. Study on the tribo-dynamic performances of water-lubricated microgroove bearings during start-up. Tribol. Int. 2020, 151, 106395. [Google Scholar]
  28. Liang, X.X.; Yan, X.P.; Ouyang, W. Experimental research on tribological and vibration performance of water-lubricated hydrodynamic thrust bearings used in marine shaft-less rim driven thrusters. Wear 2019, 426–427, 778–791. [Google Scholar] [CrossRef]
  29. Wang, L.; Pei, S.Y.; Xiong, X.Z. Study on the static performance and stability of a water-lubricated hybrid bearing with circumferential grooves and stepped recesses considering the influence of recess sizes. Tribol. Trans. 2014, 57, 36–45. [Google Scholar] [CrossRef]
  30. Constantinescu, V.N. On gas lubrication in turbulent regime. J. Basic Eng. 1964, 86, 475–482. [Google Scholar] [CrossRef]
  31. Frene, J.; Arghir, M.; Constantinescu, V. Combined thin-film and Navier–Stokes analysis in high Reynolds number lubrication. Tribol. Int. 2006, 39, 734–747. [Google Scholar] [CrossRef]
  32. Sun, J.; Gui, C.L. Hydrodynamic lubrication analysis of journal bearing considering misalignment caused by shaft deformation. Tribol. Int. 2004, 37, 841–848. [Google Scholar] [CrossRef]
  33. Lv, F.R.; Ta, N.; Rao, Z.S. Analysis of equivalent supporting point location and carrying capacity of misaligned journal bearing. Tribol. Int. 2017, 116, 26–38. [Google Scholar] [CrossRef]
  34. Gertzos, K.; Nikolakopoulos, P.; Papadopoulos, C. CFD analysis of journal bearing hydrodynamic lubrication by Bingham lubricant. Tribol. Int. 2008, 41, 1190–1204. [Google Scholar] [CrossRef]
  35. Chasalevris, A.; Sfyris, D. Evaluation of the finite journal bearing characteristics, using the exact analytical solution of the Reynolds equation. Tribol. Int. 2013, 57, 216–234. [Google Scholar] [CrossRef]
Figure 1. Schematic of textured hydrodynamic journal bearing.
Figure 1. Schematic of textured hydrodynamic journal bearing.
Lubricants 11 00362 g001
Figure 2. Grid division of a textured hydrodynamic journal bearing.
Figure 2. Grid division of a textured hydrodynamic journal bearing.
Lubricants 11 00362 g002
Figure 3. The variation of the eccentricity ratio with Sommerfeld number [31].
Figure 3. The variation of the eccentricity ratio with Sommerfeld number [31].
Lubricants 11 00362 g003
Figure 4. Comparison of oil film pressure distribution theoretical result, test value, and results from other studies in the literature. (a) Oil film pressure test bench; (b) Comparison of oil film pressure distribution [34,35].
Figure 4. Comparison of oil film pressure distribution theoretical result, test value, and results from other studies in the literature. (a) Oil film pressure test bench; (b) Comparison of oil film pressure distribution [34,35].
Lubricants 11 00362 g004
Figure 5. Reynolds number varies with eccentricity ratios and rotational speeds. (a) Reynolds number varies with eccentricities; (b) Reynolds number varies with rotational speeds.
Figure 5. Reynolds number varies with eccentricity ratios and rotational speeds. (a) Reynolds number varies with eccentricities; (b) Reynolds number varies with rotational speeds.
Lubricants 11 00362 g005
Figure 6. Variations in the oil film thickness and pressure with eccentricity and angular coordinate. (a) The variation of oil film thickness; (b) The variation of minimum oil film thickness and maximum oil film pressure with eccentricity ratio.
Figure 6. Variations in the oil film thickness and pressure with eccentricity and angular coordinate. (a) The variation of oil film thickness; (b) The variation of minimum oil film thickness and maximum oil film pressure with eccentricity ratio.
Lubricants 11 00362 g006
Figure 7. The variation in oil film pressure distribution with L/D = 1 and ε = 0.3. (a) The variation of oil film pressure under various flow regimes; (b) The distribution of oil film pressure.
Figure 7. The variation in oil film pressure distribution with L/D = 1 and ε = 0.3. (a) The variation of oil film pressure under various flow regimes; (b) The distribution of oil film pressure.
Lubricants 11 00362 g007
Figure 8. Sectional view oil film pressure of a textured bearing; (a) Circumferential oil film pressure; (b) axial oil film pressure.
Figure 8. Sectional view oil film pressure of a textured bearing; (a) Circumferential oil film pressure; (b) axial oil film pressure.
Lubricants 11 00362 g008
Figure 9. Local oil film pressure distribution of textured bearing: (a) Local oil film pressure distribution; (b) local oil film thickness distribution.
Figure 9. Local oil film pressure distribution of textured bearing: (a) Local oil film pressure distribution; (b) local oil film thickness distribution.
Lubricants 11 00362 g009
Figure 10. Variation in bearing capacity with dimple diameter, depth, spacing, and film thickness. (a) The variation of load capacity with dimple diameter; (b) The variation of load capacity with dimple depth; (c) The variation of load capacity with dimple spacing; (d) The variation of load capacity with oil film thickness.
Figure 10. Variation in bearing capacity with dimple diameter, depth, spacing, and film thickness. (a) The variation of load capacity with dimple diameter; (b) The variation of load capacity with dimple depth; (c) The variation of load capacity with dimple spacing; (d) The variation of load capacity with oil film thickness.
Lubricants 11 00362 g010
Figure 11. Variation in bearing capacity with length/diameter ratio and clearance ratio. (a) Bearing capacity with length to diameter ratio; (b) Bearing capacity with clearance ratio.
Figure 11. Variation in bearing capacity with length/diameter ratio and clearance ratio. (a) Bearing capacity with length to diameter ratio; (b) Bearing capacity with clearance ratio.
Lubricants 11 00362 g011
Figure 12. Variation in bearing capacity with rotational speed. (a) Bearing capacity with rotational peeds under various flow regimes; (b) Bearing capacity with rotational speeds under various parameters.
Figure 12. Variation in bearing capacity with rotational speed. (a) Bearing capacity with rotational peeds under various flow regimes; (b) Bearing capacity with rotational speeds under various parameters.
Lubricants 11 00362 g012
Figure 13. Variation in the attitude angle and leakage with eccentricity ratio. (a) Attitude angle with eccentricity ratio; (b) Leakage with eccentricity ratio.
Figure 13. Variation in the attitude angle and leakage with eccentricity ratio. (a) Attitude angle with eccentricity ratio; (b) Leakage with eccentricity ratio.
Lubricants 11 00362 g013
Figure 14. The variation in the axis whirl orbit of the rotor. (a) The axis whirl orbit variation of the rotor under L/D = 1, ε = 0.25 and N = 1000 rpm; (b) The axis whirl orbit variation of the rotor under L/D = 1, ε = 0.95 and N = 1000 rpm; (c) The axis whirl orbit variation of the rotor under L/D = 1, ε = 0.25 and N = 2000 rpm; (d) The axis whirl orbit variation of the rotor under L/D = 2, ε = 0.25 and N = 1000 rpm.
Figure 14. The variation in the axis whirl orbit of the rotor. (a) The axis whirl orbit variation of the rotor under L/D = 1, ε = 0.25 and N = 1000 rpm; (b) The axis whirl orbit variation of the rotor under L/D = 1, ε = 0.95 and N = 1000 rpm; (c) The axis whirl orbit variation of the rotor under L/D = 1, ε = 0.25 and N = 2000 rpm; (d) The axis whirl orbit variation of the rotor under L/D = 2, ε = 0.25 and N = 1000 rpm.
Lubricants 11 00362 g014
Figure 15. Friction and wear surface morphology of the THJB following operation as various rotational speeds; (a) SEM morphology of surface after operation at 1000 rpm; (b) SEM morphology of surface after operation at 2000 rpm.
Figure 15. Friction and wear surface morphology of the THJB following operation as various rotational speeds; (a) SEM morphology of surface after operation at 1000 rpm; (b) SEM morphology of surface after operation at 2000 rpm.
Lubricants 11 00362 g015
Table 1. Geometric parameters of a textured hydrodynamic journal bearing.
Table 1. Geometric parameters of a textured hydrodynamic journal bearing.
ParametersValues
Bearing radius R/mm100
Radial clearance, c/mm0.05–0.20
Length/diameter ratio, L/D1
Eccentricity ratio, ε0.25–0.95
Lubricant viscosity, μ/Pa·s−10.0048
Density, ρ/kg·m−3890
Mass, m/kg13.43
Rotational speed, n/rpm1000–20,000
Texture radius, rp/μm16
Texture depth, hp/μm25
Texture density, sp/%25
Table 2. Comparison results for various differential grids.
Table 2. Comparison results for various differential grids.
Differential GridMaximum Oil Film Pressure Pmax/MPa
Sun et al. [32]Lv et al. [33]FDMErrors 1Errors 2
314 × 10033.0633.033.050.03%0.15%
390 × 12533.0633.032.870.58%0.39%
628 × 20033.0633.032.561.54%1.35%
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

Mao, Y.; Li, L.; Li, D.; Zheng, J. Analysis of the Turbulent Lubrication of a Textured Hydrodynamic Journal Bearing. Lubricants 2023, 11, 362. https://doi.org/10.3390/lubricants11090362

AMA Style

Mao Y, Li L, Li D, Zheng J. Analysis of the Turbulent Lubrication of a Textured Hydrodynamic Journal Bearing. Lubricants. 2023; 11(9):362. https://doi.org/10.3390/lubricants11090362

Chicago/Turabian Style

Mao, Yazhou, Lilin Li, Daqing Li, and Jingyang Zheng. 2023. "Analysis of the Turbulent Lubrication of a Textured Hydrodynamic Journal Bearing" Lubricants 11, no. 9: 362. https://doi.org/10.3390/lubricants11090362

APA Style

Mao, Y., Li, L., Li, D., & Zheng, J. (2023). Analysis of the Turbulent Lubrication of a Textured Hydrodynamic Journal Bearing. Lubricants, 11(9), 362. https://doi.org/10.3390/lubricants11090362

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