A Numerical Investigation of the Effects of Groove Texture on the Dynamics of a Water-Lubricated Bearing–Rotor System

: This paper aims to investigate the combined effects of working condition and structural parameters of groove texture on the dynamic characteristics, stability and unbalance response of a water-lubricated hydrodynamic bearing–rotor system to avoid instability and excessive vibration of the rotor. The Navier–Stokes equation, standard K-ε model with enhanced wall treatment and Zwart– Gerber–Belamri cavitation model are considered using the commercial software Fluent to calculate the stiffness and damping coefﬁcients of a groove-textured, water-lubricated bearing based on the dynamic mesh method; the critical mass to express the stability and the unbalance response solved by the fourth order Runge–Kutta method of the rotor are calculated based on dynamic equations. The results indicate that shallower and longer groove textures can improve the direct stiffness along the load direction k yy , weaken the stiffness in the orthogonal direction k xx , improve stability and decrease the unbalance response amplitude of the water-lubricated bearing–rotor system at a greater rotational speed and smaller eccentricity ratio; however, the impact of grooves on damping parameters is not as great as it is on stiffness—there exists an optimum groove width to achieve a best dynamic performance.


Introduction
The water-lubricated hydrodynamic journal bearing (WLHB) has gained attention in the area of high-speed spindle systems and marine shafts due to the combined advantages of low power loss, being free of pollution at large rotary speed, good damping and long lifetime. Numerous researches about the static performance [1][2][3][4] and dynamic performance [5] of WLHBs have been carried out based on the Reynolds equation [2,5], Navier-Stokes using the CFD analysis method [1,3,4] or experimental methods [6]. Despite the numerous advantages of water as a lubricant, further researches to improve the load capacity and reduce friction still need to be carried out.
Recently, surface texture has been successfully used in the field of water-lubricated bearings to further improve load capacity and tribological performance [7,8]. However, a higher rotary speed and smaller radial clearance usually adopted in WLHBs easily leads to a reduction in stability and the tendency to unbalance the bearing-rotor system [9]. To avoid instability and excessive vibration of the rotor, the effect of surface texture on water-lubricated bearing-rotor systems should be fully studied.
In the past few years, the effects of shape, size and distribution of textures on the static performance such as friction, wear and carrying capacity of water-lubricated bearings have been thoroughly investigated [7,8,[10][11][12][13]. It can be concluded that the wear could be apparently reduced and the bearing capacity can be significantly enhanced with specific structural parameters for different working conditions. The mathematical models adopted usually are the Reynolds equation and Navier-Stokes equations. Cavitation, turbulence or inertia are taken into consideration partially or entirely to account for the mechanisms of textured bearings [14][15][16]. Results show that the Navier-Stokes equations are needed when film heights are in the same order or greater than pocket depth, and convective inertia effects are important [17].
To numerically evaluate the influence of surface texture on the dynamics of a rotor, the dynamic characteristics of the bearings including stiffness and damping coefficients should be determined first, which will then be substituted into dynamic equations [18][19][20][21]. Xie et al. [12] pointed out that speed has a significant influence on the equivalent stiffness. Sharma et al. [19] numerically investigated the influences of a triangular shape texture on the stability and found that lower texture depth is beneficial for direct stiffness and threshold speed. Singh and Awasthi [20] found that using the optimal distribution and aspect ratio could increase the stability apparently. However, it should be noted that for the available researches, the stiffness and damping coefficients being substituted into the dynamic rotor equations are solved by the laminar Reynolds equations.
Recently, some researchers carried out experimental investigations on the effect of texture on the stability of bearing-rotor systems [14,[21][22][23] and proved that the vibration of the water-lubricated bearing-rotor system could be apparently reduced by reasonable design of the texture.
However, a comprehensive numerical study to explore how textures affect the stability and unbalance response of the WLHB rotor system based on Navier-Stokes equation considering cavitation, turbulence and inertia has not been undertaken. This research is a follow-up study based on the previous study about static performance of grooved WLHB we made before [7] that aims to fill the gap and investigate the combined influences of working conditions (rotary speed and eccentricity ratio) and structural parameters (width, depth and length of groove texture) on the dynamic characteristics, stability and unbalance response of the water-lubricated hydrodynamic bearing-rotor system used in high-speed spindles.

Governing Equations
The groove-textured WLHB is depicted in Figure 1. Half of the bearing is modeled due to its axial symmetry perpendicular to the y-O-x plane along the z axis.
anisms of textured bearings [14][15][16]. Results show that the Navier-Stokes needed when film heights are in the same order or greater than pocket dept tive inertia effects are important [17].
To numerically evaluate the influence of surface texture on the dynam the dynamic characteristics of the bearings including stiffness and dampi should be determined first, which will then be substituted into dynamic e 21]. Xie et al. [12] pointed out that speed has a significant influence on the e ness. Sharma et al. [19] numerically investigated the influences of a triangu ture on the stability and found that lower texture depth is beneficial for d and threshold speed. Singh and Awasthi [20] found that using the optim and aspect ratio could increase the stability apparently. However, it should for the available researches, the stiffness and damping coefficients being su the dynamic rotor equations are solved by the laminar Reynolds equations Recently, some researchers carried out experimental investigations o texture on the stability of bearing-rotor systems [14,[21][22][23] and proved tha of the water-lubricated bearing-rotor system could be apparently reduced design of the texture.
However, a comprehensive numerical study to explore how textures aff ity and unbalance response of the WLHB rotor system based on Navier-St considering cavitation, turbulence and inertia has not been undertaken. Th follow-up study based on the previous study about static performance of g we made before [7] that aims to fill the gap and investigate the combined working conditions (rotary speed and eccentricity ratio) and structural para depth and length of groove texture) on the dynamic characteristics, stabil ance response of the water-lubricated hydrodynamic bearing-rotor system speed spindles.

Governing Equations
The groove-textured WLHB is depicted in Figure 1. Half of the beari due to its axial symmetry perpendicular to the y-O-x plane along the z axis In the present work, the flow behavior is solved by CFD and calculate of the Navier-Stokes equation for incompressible flow and the isothermal c mass conservation equation coupled with the momentum conservation eq follows [7,24]: The mass conservation equation is In the present work, the flow behavior is solved by CFD and calculated on the bases of the Navier-Stokes equation for incompressible flow and the isothermal conditions. The mass conservation equation coupled with the momentum conservation equations are as follows [7,24]: The momentum conservation equation is where ρ is the density of the fluid, u, v, w is the velocity vector, P is the static pressure, = τ is the stress tensor, ρ → g is the gravitational body force and → F is the external body force. For groove-textured water-lubricated bearings, turbulence and cavitation might occur owning to the low viscosity of water, high rotational speed, large film thickness and divergent region. The standard K-ε model with enhanced wall treatment will be adopted. The equation for the kinetic energy of turbulence and the equation for the dissipation rate are as follows: where k is the turbulent kinetic energy and ε is the turbulent dissipation rate. C 1ε , C 2ε , σ k and σ ε are empirical constants with values of 1.44, 1.92, 1.0 and 1.3, respectively. The Zwart-Gerber-Belamri cavitation model, which is shown as follows, is adopted [25]: where R e and R c are the mass transfer source terms connected to the vaporization and condensation processes, respectively; the vaporization pressure P v of water is 2340 Pa; and the bubble radius R B, the nucleation site volume fraction a nue , the evaporation coefficient F vap and F cond are set to default values.

Dynamic Characteristics
According to the above models, the distribution of pressure across a bearing and the equilibrium position can be calculated through iteration. Superimposed on the equilibrium position, the journal center disturbs with small displacement and velocity in each direction (x,y) when small disturbances exist. The relationship between the disturbed displacement, velocity and pressure in the film can be expressed as follows: where F x0 , F y0 are the forces at the equilibrium position; F xi , F yi are disturbed forces in each direction; ∆x, ∆y are disturbed displacements and ∆ .
x, ∆ . y are disturbed velocity. The linearized bearing water film forces using stiffness and damping coefficients are illustrated in Figure 2.

Bearing-Rotor Dynamic Model
The groove-textured WLHBs studied here are used to support the motorized s Since a motorized spindle usually works under its first critical speed, the bearin system can be regarded as rigid. To eliminate the gyroscopic effects, two identical W symmetrically support a rotor, which is shown in Figure 3.

Stability Analysis
The bearing-rotor system will start oscillating from the equilibrium positio there is a small disturbance. The stiffness and damping coefficients of the beari then substituted into the transient dynamic equations of the bearing-rotor system mx + c xx x + c xy y + k xx x + k xy y = 0 my + c yx x + c yy y + k yx x + k yy y = 0 When the rotor is disturbed by a small disturbance, the rotor will whirl aro steady position. When it is at the threshold of instability, the solution of Equation pure imaginary number, which can be written as where w is the whirling frequency of the rotor center.
Substitute solution (9) into Equation (8) The determinate of the matrix should be zero. Since both the real part and im part should be equal to zero, there will be two equations. The critical whirling fre and critical mass can be solved. By substituting the non-dimensional dynamic coe

Bearing-Rotor Dynamic Model
The groove-textured WLHBs studied here are used to support the motorized spindle. Since a motorized spindle usually works under its first critical speed, the bearing-rotor system can be regarded as rigid. To eliminate the gyroscopic effects, two identical WLHBs symmetrically support a rotor, which is shown in Figure 3.

Bearing-Rotor Dynamic Model
The groove-textured WLHBs studied here are used to support the moto Since a motorized spindle usually works under its first critical speed, the system can be regarded as rigid. To eliminate the gyroscopic effects, two ide symmetrically support a rotor, which is shown in Figure 3.

Stability Analysis
The bearing-rotor system will start oscillating from the equilibrium p there is a small disturbance. The stiffness and damping coefficients of the then substituted into the transient dynamic equations of the bearing-rotor s mx + c xx x + c xy y + k xx x + k xy y = 0 my + c yx x + c yy y + k yx x + k yy y = 0 When the rotor is disturbed by a small disturbance, the rotor will wh steady position. When it is at the threshold of instability, the solution of Eq pure imaginary number, which can be written as where w is the whirling frequency of the rotor center.
Substitute solution (9) into Equation (8) The determinate of the matrix should be zero. Since both the real part a part should be equal to zero, there will be two equations. The critical whirl

Stability Analysis
The bearing-rotor system will start oscillating from the equilibrium position when there is a small disturbance. The stiffness and damping coefficients of the bearings are then substituted into the transient dynamic equations of the bearing-rotor system: When the rotor is disturbed by a small disturbance, the rotor will whirl around its steady position. When it is at the threshold of instability, the solution of Equation (7) is a pure imaginary number, which can be written as where w is the whirling frequency of the rotor center. Substitute solution (9) into Equation (8), the expression can be obtained as follows: −mw 2 st +iw st c xx +k xx iw st c xy +k xy iw st c yx +k yx −mw 2 st +iw st c yy +k yy The determinate of the matrix should be zero. Since both the real part and imaginary part should be equal to zero, there will be two equations. The critical whirling frequency and critical mass can be solved. By substituting the non-dimensional dynamic coefficients (K = kψ 3 µΩL , C = cψ 3 µL , where ψ = h 0 /L, h 0 is radial clearance and L is the length of the bearing), the non-dimensional critical mass M st = m st ψ 3 Ω µL can be determined to assess the stability [26,27]. When the rotor mass is no greater than the critical mass, the rotor is stable; the greater the critical mass is, the more stable the rotor is.

Unbalance Response Analysis
Assuming that the residual unbalance is located in the plane of the mass center of the bearing-rotor system, the dynamic equations of motion can be expressed as follows: [M] ..
where X represents the unbalance response of the rotor center in the orthogonal direction of x and y, [C] is the damping matrix and [K] is stiffness matrix. The unbalance force ∆F is shown as follows: Equation (10) can be solved by using the fourth-order Runge-Kutta method to obtain the unbalance response of the rotor mass center (x center , y cener ).

Calculation Process
The static performances of groove-textured water-lubricated bearings can be investigated thoroughly using Fluent, as in our previous study [7]. Based on the static results, further calculation including the dynamic coefficients of the bearings, stability and unbalance response of the WLHB rotor using Fluent and Matlab are discussed in the following part. The computational procedure is depicted in Figure 4.
When the equilibrium position is obtained, the dynamic mesh method is adopted. The User Defined Function (UDF) in Fluent should be used to deal with the linear whirl motion of the rotor. In this process, since the position of the journal center changes due to the combined motions of translational and rotational motions, the traditional dynamic mesh method is no longer suitable due to the mesh distortion leading by rotational motion. To this point, two UDF macros are used [28,29]: the DEFINE_CG_MOTION macro is used to define the displacement and velocity whirling motion; the DEFINE_PROFILE macro is used to update the coordinates of the journal meshes that change according to the position of the journal.
To solve the stiffness and damping coefficients, the calculation is carried out in two steps: firstly, the disturbed velocities are assumed to be unchanged to calculate the stiffness, thus, the differences of velocities are zero (∆ . x= 0; ∆ . y = 0); secondly, the disturbed displacements are assumed to be unchanged to calculate the damping coefficients, thus, the differences of displacements are zero (∆x = 0; ∆y = 0).
In the first step, the dynamic equation can be expressed as follows: where F xi is the force in the x direction with a disturbed displacement ∆x in the +x direction; F xj is the force in the x direction with a disturbed displacement ∆x in the −x direction. Similarly, F yi , F yj are the forces in the y direction with a disturbed displacement ∆y in the +y and −y directions, respectively. When the disturbed displacement in the x direction ∆x is applied and the disturbed displacement in the y direction ∆y is zero, the stiffness k xx can be solved. Similarly, when the disturbed displacement in the y direction ∆y is applied and disturbed displacement in the x direction ∆x is zero, the stiffness k yy can be solved. Then, both of the disturbed displacements in the x and y directions are applied by substituting the stiffness k xx and k yy into Equation (13); thus, the coupled stiffness k xy and k yx can be solved.  Similarly, in the second step, the dynamic equation can be expressed as follows: The disturbed velocities in the x, y direction ∆ , ∆ are changed successively to solve the damping coefficients and . Then, both of the disturbed velocities in the x and y directions are applied by substituting the damping coefficients and into Equation (14); thus, the coupled damping coefficients and can be solved. Similarly, in the second step, the dynamic equation can be expressed as follows: Lubricants 2023, 11, 242 7 of 24 The disturbed velocities in the x, y direction ∆ .
x, ∆ . y are changed successively to solve the damping coefficients c xx and c yy . Then, both of the disturbed velocities in the x and y directions are applied by substituting the damping coefficients c xx and c yy into Equation (14); thus, the coupled damping coefficients c xy and c yx can be solved.
Then, by substituting the dynamic characteristics of the bearings into dynamic equations of rotor, the critical mass M st to express the stability and the unbalance response using the fourth-order Runge-Kutta method will be solved by Matlab codes.
The water film is divided into meshes of 220 × 150 × 5, and the sensitive analyses for meshes have been fully studied in our previous study [7].

Validation of the Calculation Process
A test rig was manufactured to verify the calculation of static performance of the WLHB in our previous study [7], which can partially validate the calculation of the CFD analysis since the static performance is the foundation of the dynamic characteristics. To verify the correctness and effectiveness of the dynamic mesh method to calculate the dynamic coefficients, two comparisons are carried out. Firstly, the stiffness coefficients of a smooth WHJB available from literature are calculated by the dynamic mesh method. Comparisons with the available results from the literature are used to partially show the correctness of calculation of the dynamic characteristics for a smooth bearing. Then, the stiffness coefficients of the grooved WHJB are calculated by both the UDF method and the difference of the load capacity method (DLC) according to the definition of stiffness to partially show the validation of the dynamic characteristics for a grooved bearing.
Since there are no available measurement data for validation of a smooth WHJB, the dynamic simulation procedure is verified with the simulation results from Ref. [30]. Comparisons of the stiffness coefficients k xx and k yy of a smooth WHJB calculated by the UDF method with the results obtained by Zhang, Yin et al. [30] are accomplished. As shown in Figure 5, according to the comparisons, the maximum error of the results by the present method and the results of Zhang, Yin et al. [30] are not greater than 16%. Then, by substituting the dynamic characteristics of the bearings into dynamic equations of rotor, the critical mass Mst to express the stability and the unbalance response using the fourth-order Runge-Kutta method will be solved by Matlab codes.
The water film is divided into meshes of 220 × 150 × 5, and the sensitive analyses for meshes have been fully studied in our previous study [7].

Validation of the Calculation Process
A test rig was manufactured to verify the calculation of static performance of the WLHB in our previous study [7], which can partially validate the calculation of the CFD analysis since the static performance is the foundation of the dynamic characteristics. To verify the correctness and effectiveness of the dynamic mesh method to calculate the dynamic coefficients, two comparisons are carried out. Firstly, the stiffness coefficients of a smooth WHJB available from literature are calculated by the dynamic mesh method. Comparisons with the available results from the literature are used to partially show the correctness of calculation of the dynamic characteristics for a smooth bearing. Then, the stiffness coefficients of the grooved WHJB are calculated by both the UDF method and the difference of the load capacity method (DLC) according to the definition of stiffness to partially show the validation of the dynamic characteristics for a grooved bearing.
Since there are no available measurement data for validation of a smooth WHJB, the dynamic simulation procedure is verified with the simulation results from Ref. [30]. Comparisons of the stiffness coefficients and of a smooth WHJB calculated by the UDF method with the results obtained by Zhang, Yin et al. [30] are accomplished. As shown in Figure 5, according to the comparisons, the maximum error of the results by the present method and the results of Zhang, Yin et al. [30] are not greater than 16%. Since the definition of stiffness is the ratio of the difference of load capacities to the difference of displacements, to further validate the correctness of UDF codes for grooved WHJB, results calculated by UDF using dynamic mesh method are compared with the ratio of the difference of load capacities to the difference of displacements using the static analysis method. Table 1 provides a list of the structural parameters of the WLHB. The eccentricity ratio is 0.3. There are ten grooves located in the region of 210°~240° and these are oriented axially through the bearing.

Item
Value Radial clearance h0 (um) 50 Since the definition of stiffness is the ratio of the difference of load capacities to the difference of displacements, to further validate the correctness of UDF codes for grooved WHJB, results calculated by UDF using dynamic mesh method are compared with the ratio of the difference of load capacities to the difference of displacements using the static analysis method. Table 1 provides a list of the structural parameters of the WLHB. The eccentricity ratio is 0.3. There are ten grooves located in the region of 210 •~2 40 • and these are oriented axially through the bearing. Table 1. Initial structural parameters of grooved WLJB.

Item Value
Radial clearance h 0 (um) 50 10 Groove length (mm) 50 As illustrated in Figure 6, the stiffness coefficients obtained by the two methods match quite well, with a maximum difference of 6.6% at a rotary speed of 10,000 rpm for the direct stiffness k yy . Since the damping coefficients are defined as the response of load capacities to disturbed velocities, which cannot be solved by the difference of load capacity using the static method, the damping coefficient calculation cannot directly be verified. However, the calculation processes using the UDFs of stiffness and damping coefficients are quite similar, only replacing the disturbed displacement by a disturbed velocity. So, the applicability and reliability of the UDF codes is verified to obtain the dynamic coefficients for the groove-textured water-lubricated bearings.
ubricants 2023, 11, x FOR PEER REVIEW Length L (mm) 50 Groove width (°) 0.1 Groove depth (µm) 10 Groove length (mm) 50 As illustrated in Figure 6, the stiffness coefficients obtained by the two match quite well, with a maximum difference of 6.6% at a rotary speed of 10,00 the direct stiffness . Since the damping coefficients are defined as the respon capacities to disturbed velocities, which cannot be solved by the difference of lo ity using the static method, the damping coefficient calculation cannot directly b However, the calculation processes using the UDFs of stiffness and damping co are quite similar, only replacing the disturbed displacement by a disturbed ve the applicability and reliability of the UDF codes is verified to obtain the dynam cients for the groove-textured water-lubricated bearings.

Results and Discussion
The combined influences of working condition and structural parameters texture on the dynamic coefficients, stability and unbalance response of a rotor s by both the smooth and grooved WLHBs are thoroughly studied. As shown in the groove length ratio λl is calculated by dividing the axial length of the groo half of the axial length of the bearing (lhb). When λl is 1, the grooves are axial th bearing. Initial parameters of the grooved WLJB are listed in Table 1. In the follo cussion, when the increase or decrease in the dynamic coefficients is beneficial namics of the bearing-rotor system, the groove texture is considered to be able to the dynamic coefficients.

Results and Discussion
The combined influences of working condition and structural parameters of groove texture on the dynamic coefficients, stability and unbalance response of a rotor supported by both the smooth and grooved WLHBs are thoroughly studied. As shown in Figure 7, the groove length ratio λ l is calculated by dividing the axial length of the groove (l g ) by half of the axial length of the bearing (l hb ). When λ l is 1, the grooves are axial through the bearing. Initial parameters of the grooved WLJB are listed in Table 1. In the following discussion, when the increase or decrease in the dynamic coefficients is beneficial to the dynamics of the bearing-rotor system, the groove texture is considered to be able to improve the dynamic coefficients. Length L (mm) 50 Groove width (°) 0.1 Groove depth (µm) 10 Groove length (mm) 50 As illustrated in Figure 6, the stiffness coefficients obtained by the two methods match quite well, with a maximum difference of 6.6% at a rotary speed of 10,000 rpm for the direct stiffness . Since the damping coefficients are defined as the response of load capacities to disturbed velocities, which cannot be solved by the difference of load capacity using the static method, the damping coefficient calculation cannot directly be verified. However, the calculation processes using the UDFs of stiffness and damping coefficients are quite similar, only replacing the disturbed displacement by a disturbed velocity. So, the applicability and reliability of the UDF codes is verified to obtain the dynamic coefficients for the groove-textured water-lubricated bearings.

Results and Discussion
The combined influences of working condition and structural parameters of groove texture on the dynamic coefficients, stability and unbalance response of a rotor supported by both the smooth and grooved WLHBs are thoroughly studied. As shown in Figure 7, the groove length ratio λl is calculated by dividing the axial length of the groove (lg) by half of the axial length of the bearing (lhb). When λl is 1, the grooves are axial through the bearing. Initial parameters of the grooved WLJB are listed in Table 1. In the following discussion, when the increase or decrease in the dynamic coefficients is beneficial to the dynamics of the bearing-rotor system, the groove texture is considered to be able to improve the dynamic coefficients.

Dynamic Characteristics of Groove-Textured WLHB
The combined influences of working condition and structural parameters of groove texture on the dynamic coefficients, stability and unbalance response of a rotor supported

Dynamic Characteristics of Groove-Textured WLHB
The combined influences of working condition and structural parameters of groove texture on the dynamic coefficients, stability and unbalance response of a rotor supported by both the smooth and grooved WLHBs were thoroughly studied. Figure 8a-h show the effects of rotary speed on stiffness and damping coefficients of both the smooth and groove-textured, water-lubricated, hydrodynamic bearings. The width and depth of a groove are 0.1 • and 10 µm, respectively.

Effects of Rotary Speed
As shown in the figure, compared with the smooth bearing, the direct stiffness coefficients k xx are smaller for the groove-textured bearing, and the differences become increasingly prominent with higher rotary speeds. However, the direct stiffness coefficients k yy along the main load direction of the groove-textured bearing are greater than those of the smooth bearing; the higher the rotational speed becomes, the greater the differences become. The groove texture can improve the stiffness along the load direction and weaken the stiffness in the orthogonal direction, especially at greater rotary speed. The reason for this is probably that larger rotary speed could promote an extra hydrodynamic effect along the load direction; thus, the stiffnesses k yy of the groove-textured bearing are larger than those of the smooth bearing. Furthermore, the greater the rotary speed is, the more significant the hydrodynamic effect it generates. However, in the orthogonal direction, the reduction exceeds the improvement introduced by the groove texture.
For the coupled stiffness k xy and k yx illustrated in Figure 8b,c, the groove texture decreases the absolute values compared with those of the smooth bearing. The reason for this is that the presence of the groove texture can produce several small convergence regions, leading to changes of the cavitation area and distribution of pressure. The reduction in cavitation area can improve the hydrodynamic effect along the load direction and weaken the hydrodynamic effect along the orthogonal direction; furthermore, instead of one convergence region, several small convergence regions may break the generation of whirling motion. As shown in the figure, compared with the smooth bearing, the direct stiffness coefficients are smaller for the groove-textured bearing, and the differences become increasingly prominent with higher rotary speeds. However, the direct stiffness coefficients along the main load direction of the groove-textured bearing are greater than those of the smooth bearing; the higher the rotational speed becomes, the greater the differences become. The groove texture can improve the stiffness along the load direction and weaken the stiffness in the orthogonal direction, especially at greater rotary speed. The reason for this is probably that larger rotary speed could promote an extra hydrodynamic effect along the load direction; thus, the stiffnesses of the groove-textured bearing are larger than those of the smooth bearing. Furthermore, the greater the rotary speed is, the more significant the hydrodynamic effect it generates. However, in the orthogonal direction, the reduction exceeds the improvement introduced by the groove texture.
For the coupled stiffness and illustrated in Figure 8b,c, the groove texture decreases the absolute values compared with those of the smooth bearing. The reason for this is that the presence of the groove texture can produce several small convergence regions, leading to changes of the cavitation area and distribution of pressure. The reduction in cavitation area can improve the hydrodynamic effect along the load direction and weaken the hydrodynamic effect along the orthogonal direction; furthermore, instead of one convergence region, several small convergence regions may break the generation of whirling motion.  Figure 8e-h illustrate the influences of the groove texture on damping coefficients of both the smooth and grooved bearings. The study demonstrates that the damping coefficients increase with respect to rotary speed. The variation amplitudes of direct damping coefficients and are smaller than those of coupled ones and the of the groovetextured bearing gradually exceeds that of the smooth one. However, the effect of the groove texture on with rotary speed is quite opposite to that in the orthogonal direction . As illustrated in Figure 8f,g, the groove texture increases the coupled damping coefficients and . It should be noted that no matter the improvement or weakness, with the increased rotational speed, the groove texture has little effect on the damping coefficients compared to stiffness.

Effects of Eccentricity Ratio
The influences of eccentricity ratio on the stiffness coefficients are shown in Figure  9a-d. The width, depth and groove length ratio λl of a groove are 0.1°, 10 µm and 1, respectively. As shown in Figure 9, by comparing the stiffness coefficients of groove-textured WLHB and smooth WLHB, it is found that the stiffness along the direction of load is significantly improved by the groove texture, especially at higher rotary speed and smaller eccentricity ratio. Operating at 20,000 rpm with an eccentricity ratio of 0.1, the value of stiffness of groove-textured WLHB is almost double that of smooth WLHB. Besides, the magnitudes of coupled stiffness coefficients , are decreased by the groove texture. The reason for this is probably that along the axis of the applied force, the existence of the groove texture generates additional hydrodynamic pressure and weakens the coupling motion of water. Figure 9e-h show the influence of coupled eccentricity ratio and groove texture on the damping coefficients. According to the results, at a larger eccentricity ratio and greater rotary speed, the direct damping coefficient in the direction of the applied load is Figure 8e-h illustrate the influences of the groove texture on damping coefficients of both the smooth and grooved bearings. The study demonstrates that the damping coefficients increase with respect to rotary speed. The variation amplitudes of direct damping coefficients c xx and c yy are smaller than those of coupled ones and the c yy of the groove-textured bearing gradually exceeds that of the smooth one. However, the effect of the groove texture on c yy with rotary speed is quite opposite to that in the orthogonal direction c xx . As illustrated in Figure 8f,g, the groove texture increases the coupled damping coefficients c xy and c yx . It should be noted that no matter the improvement or weakness, with the increased rotational speed, the groove texture has little effect on the damping coefficients compared to stiffness.

Effects of Eccentricity Ratio
The influences of eccentricity ratio on the stiffness coefficients are shown in Figure 9a-d. The width, depth and groove length ratio λ l of a groove are 0.1 • , 10 µm and 1, respectively. As shown in Figure 9, by comparing the stiffness coefficients of groove-textured WLHB and smooth WLHB, it is found that the stiffness along the direction of load k yy is significantly improved by the groove texture, especially at higher rotary speed and smaller eccentricity ratio. Operating at 20,000 rpm with an eccentricity ratio of 0.1, the value of stiffness k yy of groove-textured WLHB is almost double that of smooth WLHB. Besides, the magnitudes of coupled stiffness coefficients k xy , k xy are decreased by the groove texture. The reason for this is probably that along the axis of the applied force, the existence of the groove texture generates additional hydrodynamic pressure and weakens the coupling motion of water. Figure 9e-h show the influence of coupled eccentricity ratio and groove texture on the damping coefficients. According to the results, at a larger eccentricity ratio and greater rotary speed, the direct damping coefficient in the direction of the applied load c yy is greater than that in the perpendicular direction c xx . It is noted that with an increased eccentricity ratio, the direct damping coefficients c xx and c yy are little affected by the groove texture while the coupled damping coefficients are improved by the groove texture by no more than 10%.
bricants 2023, 11, x FOR PEER REVIEW 11 greater than that in the perpendicular direction . It is noted that with an increased centricity ratio, the direct damping coefficients and are little affected by the gr texture while the coupled damping coefficients are improved by the groove texture b more than 10%.   Figure 10a-h present the dynamic coefficients of WLHBs versus groove dept 15,000 rpm, 17,500 rpm and 20,000 rpm and an eccentricity ratio of 0.1. The groove w is 0.1° and the groove length ratio λl is 1.

Effects of Groove Depth
As shown in Figure 10a,b, with an increase in groove depth from 10 µm to 30 µm stiffness is significantly improved; however, when the groove depth becomes dee the improvement by the groove texture on remains almost unchanged. For exam for the value at a rotary speed of 20,000 rpm and eccentricity ratio of 0.1, with a gro depth of 10 µm is approximately double that of of the smooth bearing. Besides absolute values of coupled stiffness and of the groove-textured WLHB are sm than those of the smooth bearing. The results indicate that a shallower groove texture improve the direct stiffness in the direction of the applied load more significantly a deeper one. The reason for this is probably that with the groove depth increase from µm to 30 µm, the flow vortex generated at the bottom of groove is enhanced significa leading to a sharp variation in the stiffness; however, with further increased groove de the effect of flow vortex at the bottom of the groove weakens gradually, leading to a ne unchanged variation in stiffness.
As shown in Figure 10e-h, compared with smooth WLHB, for most of the cases groove texture can improve the damping coefficients by no more than 10% with incre groove depth. When the groove depth increases from 10 µm to 20 µm, the damping c ficients of groove-textured WLHB change slightly; however, with increased groove de the damping coefficients remain almost unchanged. The reason for this is that dampin a function of the pressure variation for a given squeeze velocity, which is not depen on the groove depth. As shown in Figure 10a,b, with an increase in groove depth from 10 µm to 30 µm, the stiffness k yy is significantly improved; however, when the groove depth becomes deeper, the improvement by the groove texture on k yy remains almost unchanged. For example, for the value at a rotary speed of 20,000 rpm and eccentricity ratio of 0.1, k yy with a groove depth of 10 µm is approximately double that of k yy of the smooth bearing. Besides, the absolute values of coupled stiffness k yx and k yx of the groove-textured WLHB are smaller than those of the smooth bearing. The results indicate that a shallower groove texture can improve the direct stiffness in the direction of the applied load k yy more significantly than a deeper one. The reason for this is probably that with the groove depth increase from 10 µm to 30 µm, the flow vortex generated at the bottom of groove is enhanced significantly, leading to a sharp variation in the stiffness; however, with further increased groove depth, the effect of flow vortex at the bottom of the groove weakens gradually, leading to a nearly unchanged variation in stiffness.     As shown in Figure 10e-h, compared with smooth WLHB, for most of the cases, the groove texture can improve the damping coefficients by no more than 10% with increased groove depth. When the groove depth increases from 10 µm to 20 µm, the damping coefficients of groove-textured WLHB change slightly; however, with increased groove depth, the damping coefficients remain almost unchanged. The reason for this is that damping is a function of the pressure variation for a given squeeze velocity, which is not dependent on the groove depth.

Effects of Groove Width
The effects of groove width on stiffness are illustrated in Figure 11a-d. The eccentricity ratio and groove length ratio λ l are 0.1 and 1, respectively. As shown in Figure 11, with the increased groove width, the stiffness either shows an increase-decrease-increase trend or decrease-increase-decrease trend. These tendencies are likely to be caused by a combination of turbulence and flow vortex. As the width of the groove is increased, the cavitation area will change and the flow vortex generated at the bottom of the texture is reduced. Shrinking of the cavitation surface and reduced flow vortex can improve both the load capacity and the stiffness k yy . Therefore, the overall load capacity and stiffness in the direction of load application can only be improved by the texture if the positive impact is greater than the negative impact caused by the confluence of flow vortex and cavitation.
As for the damping coefficients presented in Figure 11e-h, c yy , c xy and c yx of the surface-textured WLHB are generally greater than those of the smooth WLHB. This is an indication that the groove texture weakens the whirling motion of the bearing. It is noted that the increase-decrease-increase trend or the opposite trend is also found. Two peaks for c yy are obtained at the groove widths of 0.1 and 1.1 degrees for bearings with a groove depth of 10 µm as well as at the groove width of 0.3 and 0.9 degrees for bearings with a groove depth of 20 µm. It should be noted that for bearings with a greater rotary speed, the effect of groove width on damping coefficients is more significant than at a lower rotary speed.

Effects of Groove Width
The effects of groove width on stiffness are illustrated in Figure 11a-d. The eccentricity ratio and groove length ratio λl are 0.1 and 1, respectively. As shown in Figure 11, with the increased groove width, the stiffness either shows an increase-decrease-increase trend or decrease-increase-decrease trend. These tendencies are likely to be caused by a combination of turbulence and flow vortex. As the width of the groove is increased, the cavitation area will change and the flow vortex generated at the bottom of the texture is reduced. Shrinking of the cavitation surface and reduced flow vortex can improve both the load capacity and the stiffness . Therefore, the overall load capacity and stiffness in the direction of load application can only be improved by the texture if the positive impact is greater than the negative impact caused by the confluence of flow vortex and cavitation.
As for the damping coefficients presented in Figure 11e-h, , and of the surface-textured WLHB are generally greater than those of the smooth WLHB. This is an indication that the groove texture weakens the whirling motion of the bearing. It is noted that the increase-decrease-increase trend or the opposite trend is also found. Two peaks for are obtained at the groove widths of 0.1 and 1.1 degrees for bearings with a groove depth of 10 µm as well as at the groove width of 0.3 and 0.9 degrees for bearings with a groove depth of 20 µm. It should be noted that for bearings with a greater rotary speed, the effect of groove width on damping coefficients is more significant than at a lower rotary speed.

Effect of Groove Length Ratio
Figure 12a-h show the effect of groove length ratio on the dynamic coefficients. Two symmetrical grooves are distributed along the axis when the ratio of the groove length is less than 1. According to the results, compared with the stiffness coefficients of the smooth bearing, the of the groove-textured WLHB is greater, while the absolute values of and of the groove-textured WLHB are smaller. For the three cases with a groove length ratio of 1, the improvement of the direct stiffness coefficient in the load application direction is 91.7%, 96.0% and 98.4%, respectively, indicating the significant improvement of grooves on the dynamic coefficients. The additional hydrodynamic effect produced by the groove texture is probably responsible for this. The longer the groove texture is, the more significant the extra hydrodynamic effect that is generated.
As for the damping coefficients illustrated in Figure 12e-h, and increase while the orthogonal damping coefficients decrease with increasing groove length ratio. Except for , all the other damping coefficients of the textured bearing are greater than those of the smooth bearing, indicating that the groove texture can suppress the whirling motion. Two symmetrical grooves are distributed along the axis when the ratio of the groove length is less than 1. According to the results, compared with the stiffness coefficients of the smooth bearing, the k yy of the groove-textured WLHB is greater, while the absolute values of k xy and k yx of the groove-textured WLHB are smaller. For the three cases with a groove length ratio of 1, the improvement of the direct stiffness coefficient in the load application direction is 91.7%, 96.0% and 98.4%, respectively, indicating the significant improvement of grooves on the dynamic coefficients. The additional hydrodynamic effect produced by the groove texture is probably responsible for this. The longer the groove texture is, the more significant the extra hydrodynamic effect that is generated.
As for the damping coefficients illustrated in Figure 12e-h, c yx and c yy increase while the orthogonal damping coefficients decrease with increasing groove length ratio. Except for c xx , all the other damping coefficients of the textured bearing are greater than those of the smooth bearing, indicating that the groove texture can suppress the whirling motion.    Figure 13 presents the dimensionless critical mass of a rotor supported by two identical WLHBs as a function of rotary speed. As can be seen in Figure 13, as the rotational speed increases, the critical masses of both the grooved and smooth WLHBs decrease, indicating that the stability of the rotor gradually deteriorates. When the eccentricity ratio is 0.1, compared with the smooth bearing, the stability of the groove-textured WLHB is improved when the rotary speed varies from 12,500 rpm to 22,500 rpm. The stability of the rotor supported by the groove-textured bearings improves significantly, from approximately 2.4% to 20.9% with an eccentricity ratio of 0.1, when the rotary speed increases to 22,500 rpm. The groove texture has a more significant stabilizing effect at higher rotary speeds. However, when the rotary speed is 12,500 rpm, the critical mass of the grooved bearing is approximately 2.0% and 5.9% smaller than that of the smooth one at eccentricity ratios of 0.3 and 0.5, respectively, indicating that the groove texture plays a negative role under these conditions. With the increased rotary speed, the difference in critical mass between the grooved and smooth bearings grows steadily, illustrating that grooves can play a more obvious positive effect on the stability of the grooved WLHB rotor system at greater rotary speeds. However, for a low rotary speed and large eccentricity ratio, the WLHB might be weakened by the groove texture.

Effects of Rotary Speed
indicating that the stability of the rotor gradually deteriorates. When the e is 0.1, compared with the smooth bearing, the stability of the groove-tex improved when the rotary speed varies from 12,500 rpm to 22,500 rpm. the rotor supported by the groove-textured bearings improves significant imately 2.4% to 20.9% with an eccentricity ratio of 0.1, when the rotary sp 22,500 rpm. The groove texture has a more significant stabilizing effect speeds. However, when the rotary speed is 12,500 rpm, the critical mass bearing is approximately 2.0% and 5.9% smaller than that of the smooth on ratios of 0.3 and 0.5, respectively, indicating that the groove texture play under these conditions. With the increased rotary speed, the difference between the grooved and smooth bearings grows steadily, illustrating t play a more obvious positive effect on the stability of the grooved WLHB greater rotary speeds. However, for a low rotary speed and large eccen WLHB might be weakened by the groove texture.
The cause of this can be attributed to the improvement of direct stiff damping coefficient in the direction of the applied load, as well as the wea ing motion due to the existence of the groove texture with proper parame decrease in coupled stiffness and increase in coupled damping coefficien tured WLHB are also important factors influencing the stability of the rot  Figure 14 presents the influences of eccentricity ratio on the stabili rotor system. It can be deduced that the stability becomes more pronounc tricity ratio increases. Compared with the smooth bearing, the stability 2.7~20.5% when the bearing is rotating at 20,000 rpm with an eccentricity 0.6; however, when the rotary speed decreases to 15,000 rpm, only when ratio is not greater than 0.4 can the stability be improved by the groove tex for this is probably that the promotion of whirling motion induced by the pled stiffness coefficients exceeds the suppression of whirling motion by efficients at a low rotary speed and large eccentricity ratio. Furthermore, speed, the extra hydrodynamic effect generated by the groove texture c the stability of the rotor. The cause of this can be attributed to the improvement of direct stiffness and direct damping coefficient in the direction of the applied load, as well as the weakening of whirling motion due to the existence of the groove texture with proper parameter settings. The decrease in coupled stiffness and increase in coupled damping coefficients of groovetextured WLHB are also important factors influencing the stability of the rotor. Figure 14 presents the influences of eccentricity ratio on the stability of the WLHB rotor system. It can be deduced that the stability becomes more pronounced as the eccentricity ratio increases. Compared with the smooth bearing, the stability is improved by 2.7~20.5% when the bearing is rotating at 20,000 rpm with an eccentricity ratio from 0.1 to 0.6; however, when the rotary speed decreases to 15,000 rpm, only when the eccentricity ratio is not greater than 0.4 can the stability be improved by the groove texture. The reason for this is probably that the promotion of whirling motion induced by the increase in coupled stiffness coefficients exceeds the suppression of whirling motion by the damping coefficients at a low rotary speed and large eccentricity ratio. Furthermore, at a high rotary speed, the extra hydrodynamic effect generated by the groove texture can also improve the stability of the rotor.  Figure 15 shows the stability with respect to groove depth. Th and the groove length ratio is 1. As shown in Figure 15, the critica textured bearing are greater than those of the smooth bearing whe creases from 10 µm to 50 µm. As the depth of the groove grows, the Specifically, when the groove depth is 50 µm and rotary speed is mass of the grooved bearing is only 0.53% greater than that of the critical mass of the grooved bearing is 20.53% larger when the groo the rotary speed is 20,000 rpm. The results indicate that the stabi cated bearing-rotor can be improved at a high rotary speed with a

Effect of Groove Width
The effect of groove width on rotor stability is shown in Figure ratio is 1. According to Figure 16, the critical masses in all cases of the smooth bearing. When the rotary speed is 20,000 rpm and g the stability is improved by a maximum of 21%, approximately. W  Figure 15 shows the stability with respect to groove depth. The groove width is 0.1 • and the groove length ratio is 1. As shown in Figure 15, the critical masses of the groove-textured bearing are greater than those of the smooth bearing when the groove depth increases from 10 µm to 50 µm. As the depth of the groove grows, the critical mass decreases. Specifically, when the groove depth is 50 µm and rotary speed is 15,000 rpm, the critical mass of the grooved bearing is only 0.53% greater than that of the smooth one, while the critical mass of the grooved bearing is 20.53% larger when the groove depth is 10 µm and the rotary speed is 20,000 rpm. The results indicate that the stability of the water-lubricated bearing-rotor can be improved at a high rotary speed with a shallow groove.  Figure 15 shows the stability with respect to groove depth. The g and the groove length ratio is 1. As shown in Figure 15, the critical ma textured bearing are greater than those of the smooth bearing when th creases from 10 µm to 50 µm. As the depth of the groove grows, the criti Specifically, when the groove depth is 50 µm and rotary speed is 15,0 mass of the grooved bearing is only 0.53% greater than that of the sm critical mass of the grooved bearing is 20.53% larger when the groove the rotary speed is 20,000 rpm. The results indicate that the stability cated bearing-rotor can be improved at a high rotary speed with a sha

Effect of Groove Width
The effect of groove width on rotor stability is shown in Figure 16 ratio is 1. According to Figure 16, the critical masses in all cases are of the smooth bearing. When the rotary speed is 20,000 rpm and groo the stability is improved by a maximum of 21%, approximately. Whe is 10 µm, the critical mass reaches its peak value at groove widths of degrees for both 15,000 rpm and 20,000 rpm; however, when groove d µm, the best groove width to improve the stability is 1.1 degrees. Gen an optimum groove width to improve stability.

Effect of Groove Width
The effect of groove width on rotor stability is shown in Figure 16. The groove length ratio λ l is 1. According to Figure 16, the critical masses in all cases are greater than those of the smooth bearing. When the rotary speed is 20,000 rpm and groove depth is 10 µm, the stability is improved by a maximum of 21%, approximately. When the groove depth is 10 µm, the critical mass reaches its peak value at groove widths of 1.3 degrees and 1.1 degrees for both 15,000 rpm and 20,000 rpm; however, when groove depth increases to 20 µm, the best groove width to improve the stability is 1.1 degrees. Generally, there exists an optimum groove width to improve stability. the stability is improved by a maximum of 21%, approximately. When the g is 10 µm, the critical mass reaches its peak value at groove widths of 1.3 deg degrees for both 15,000 rpm and 20,000 rpm; however, when groove depth in µm, the best groove width to improve the stability is 1.1 degrees. Generally an optimum groove width to improve stability.

Effect of Groove Length
The effect of groove length on the stability of the water-lubricated bearing-rotor is shown in Figure 17. With the increased length ratio, the stability of the grooved bearingrotor system increases gradually. For the rotor at a rotary speed of 15,000 rpm, the groove texture can weaken the stability when the length ratio is less than 0.5. However, for the rotor at a rotary speed of 20,000 rpm, the groove texture improves the stability significantly. Further, the greater the ratio of the length is, the more significant the contribution of the groove texture to stability is. When the length ratio is 1, the stability is improved by about 20.5% with groove widths of 0.1 degrees and 1.1 degrees at a rotary speed of 20,000 rpm. However, when the length ratio decreases to 0.2 for a bearing with a groove width of 0.9 degrees and depth of 20 µm at 15,000 rpm, the critical mass of the groove-textured bearing is only about 0.78% larger than that of the smooth bearing. Results show that the improvement by groove texture is more significant at a larger groove length ratio and greater rotary speed.

Effect of Groove Length
The effect of groove length on the stability of the water-lubricated bea shown in Figure 17. With the increased length ratio, the stability of the groo rotor system increases gradually. For the rotor at a rotary speed of 15,000 rpm texture can weaken the stability when the length ratio is less than 0.5. How rotor at a rotary speed of 20,000 rpm, the groove texture improves the sta cantly. Further, the greater the ratio of the length is, the more significant the of the groove texture to stability is. When the length ratio is 1, the stability is about 20.5% with groove widths of 0.1 degrees and 1.1 degrees at a rotary sp rpm. However, when the length ratio decreases to 0.2 for a bearing with a g of 0.9 degrees and depth of 20 µm at 15,000 rpm, the critical mass of the gro bearing is only about 0.78% larger than that of the smooth bearing. Results s improvement by groove texture is more significant at a larger groove leng greater rotary speed.

Unbalance Response of Groove-Textured, Water-Lubricated Bearing-Rotor Sys
To evaluate the effects of groove texture on the unbalance response of t amplitude of unbalance response is studied. The residual unbalance is 1 × The time interval is 1 × 10 s and 10,000 time steps are run to ensure a stead

Effect of Rotary Speed
The combined effects of rotary speed and groove texture on the unbala of the water-lubricated bearing-rotor system is illustrated in Figure 18. The g is 0.1 degrees, the depth is 10 µm and the ratio of length is 1. As shown in the the unbalance responses of the groove-textured and smooth water-lubricated

Unbalance Response of Groove-Textured, Water-Lubricated Bearing-Rotor System
To evaluate the effects of groove texture on the unbalance response of the rotor, the amplitude of unbalance response is studied. The residual unbalance is 1 × 10 6 (kg∆m). The time interval is 1 × 10 −6 s and 10,000 time steps are run to ensure a steady state.

Effect of Rotary Speed
The combined effects of rotary speed and groove texture on the unbalance response of the water-lubricated bearing-rotor system is illustrated in Figure 18. The groove width is 0.1 degrees, the depth is 10 µm and the ratio of length is 1. As shown in the figure, both the unbalance responses of the groove-textured and smooth water-lubricated bearingrotor system increase with the rotary speed. It should be noted that at speeds above 17,500 rpm and eccentricity ratio of 0.1, grooves can decrease the amplitude in both the x and y directions. When the rotary speed is 20,000 rpm and eccentricity ratio is 0.1, the amplitudes of unbalance in the x and y directions are reduced by approximately 10.51 and 3.79%; when the rotary speed increases to 22,500 rpm, the amplitudes of unbalance in the x and y directions are reduced by approximately 13.55 and 4.65%. However, at a lower rotary speed or larger eccentricity ratio, the groove will increase the amplitude of the unbalance response.

Effect of Eccentricity Ratio
The unbalance responses of the smooth bearing-rotor and grooved bearing-rotor in terms of eccentricity ratio are presented in Figure 19. As shown in Figure 19a, the amplitudes of unbalance in both the x and y directions decrease with increased eccentricity ratio due to the increased direct stiffness coefficients and damping coefficients. It should be noted that only with an eccentricity ratio of 0.1 and a rotational speed of 20,000 rpm, the groove texture can decrease the unbalance amplitude by 10.51% and 3.79%, respectively. When the eccentricity ratio is 0.1, the groove texture can reduce the amplitude at rotational speeds of 15,000 rpm, 17,500 rpm and 20,000 rpm; when it increases to 0.3, the groove texture will be beneficial to the dynamics of the rotor with a rotational speed of 17,500 rpm or 20,000 rpm. However, when the eccentricity ratio increases to 0.5, the amplitude can only be reduced by the groove texture at a high rotational speed of 20,000 rpm.
In general, at a small eccentricity ratio and great rotational speed, the groove texture is beneficial to the dynamics of rotor; however, at a small eccentricity ratio and lower rotary speed, the groove texture significantly increases the unbalance response; with the further increased eccentricity ratio, the rotor dynamics are not significantly affected by the groove texture.

Effect of Eccentricity Ratio
The unbalance responses of the smooth bearing-rotor and grooved bearing-rotor in terms of eccentricity ratio are presented in Figure 19. As shown in Figure 19a, the amplitudes of unbalance in both the x and y directions decrease with increased eccentricity ratio due to the increased direct stiffness coefficients and damping coefficients. It should be noted that only with an eccentricity ratio of 0.1 and a rotational speed of 20,000 rpm, the groove texture can decrease the unbalance amplitude by 10.51% and 3.79%, respectively. When the eccentricity ratio is 0.1, the groove texture can reduce the amplitude at rotational speeds of 15,000 rpm, 17,500 rpm and 20,000 rpm; when it increases to 0.3, the groove texture will be beneficial to the dynamics of the rotor with a rotational speed of 17,500 rpm or 20,000 rpm. However, when the eccentricity ratio increases to 0.5, the amplitude can only be reduced by the groove texture at a high rotational speed of 20,000 rpm.

Effect of Eccentricity Ratio
The unbalance responses of the smooth bearing-rotor and grooved bearing-rotor in terms of eccentricity ratio are presented in Figure 19. As shown in Figure 19a, the amplitudes of unbalance in both the x and y directions decrease with increased eccentricity ratio due to the increased direct stiffness coefficients and damping coefficients. It should be noted that only with an eccentricity ratio of 0.1 and a rotational speed of 20,000 rpm, the groove texture can decrease the unbalance amplitude by 10.51% and 3.79%, respectively. When the eccentricity ratio is 0.1, the groove texture can reduce the amplitude at rotational speeds of 15,000 rpm, 17,500 rpm and 20,000 rpm; when it increases to 0.3, the groove texture will be beneficial to the dynamics of the rotor with a rotational speed of 17,500 rpm or 20,000 rpm. However, when the eccentricity ratio increases to 0.5, the amplitude can only be reduced by the groove texture at a high rotational speed of 20,000 rpm.
In general, at a small eccentricity ratio and great rotational speed, the groove texture is beneficial to the dynamics of rotor; however, at a small eccentricity ratio and lower rotary speed, the groove texture significantly increases the unbalance response; with the further increased eccentricity ratio, the rotor dynamics are not significantly affected by the groove texture. Eccentricity ratio  The effects of groove depth on unbalance response are shown in Figure 20. The groove texture can decrease the unbalance response in both the x and y directions with respect to groove depth when the rotary speed is 20,000 rpm. The most significant de- In general, at a small eccentricity ratio and great rotational speed, the groove texture is beneficial to the dynamics of rotor; however, at a small eccentricity ratio and lower rotary speed, the groove texture significantly increases the unbalance response; with the further increased eccentricity ratio, the rotor dynamics are not significantly affected by the groove texture.

Effect of Groove Depth
The effects of groove depth on unbalance response are shown in Figure 20. The groove texture can decrease the unbalance response in both the x and y directions with respect to groove depth when the rotary speed is 20,000 rpm. The most significant decrease in groove texture on unbalance response exists at a rotary speed of 20,000 rpm and groove depth of 10 µm, which is 10.51% and 3.79% in the directions of x and y, respectively. It can be seen that the reduction in unbalance response is greater at shallower depth. However, at speeds below 20,000 rpm, the groove texture gradually increases the amplitude with groove depth. respectively. It can be seen that the reduction in unbalance response is greater at shallower depth. However, at speeds below 20,000 rpm, the groove texture gradually increases the amplitude with groove depth. Groove depth (μm) Smooth Figure 21 shows the influence of groove width on unbalance response. According to the results, the amplitudes in both the x and y directions increase first, then decrease and then increase with the increased eccentricity ratio. Two valleys exist when the groove width increases from 0.1 degrees to 1.3 degrees. The minimum unbalance amplitudes of the grooved WLHB rotor for the three cases are 5.1%, 6.5% and 12.6% smaller than those of the smooth WLHB rotor in the x direction and 5.5%, 6.3% and 14.8% smaller than those of the smooth WLHB rotor in the y direction.   Figure 21 shows the influence of groove width on unbalance response. According to the results, the amplitudes in both the x and y directions increase first, then decrease and then increase with the increased eccentricity ratio. Two valleys exist when the groove width increases from 0.1 degrees to 1.3 degrees. The minimum unbalance amplitudes of the grooved WLHB rotor for the three cases are 5.1%, 6.5% and 12.6% smaller than those of the smooth WLHB rotor in the x direction and 5.5%, 6.3% and 14.8% smaller than those of the smooth WLHB rotor in the y direction. respectively. It can be seen that the reduction in unbalance response is greater at shallower depth. However, at speeds below 20,000 rpm, the groove texture gradually increases the amplitude with groove depth.  Figure 21 shows the influence of groove width on unbalance response. According to the results, the amplitudes in both the x and y directions increase first, then decrease and then increase with the increased eccentricity ratio. Two valleys exist when the groove width increases from 0.1 degrees to 1.3 degrees. The minimum unbalance amplitudes of the grooved WLHB rotor for the three cases are 5.1%, 6.5% and 12.6% smaller than those of the smooth WLHB rotor in the x direction and 5.5%, 6.3% and 14.8% smaller than those of the smooth WLHB rotor in the y direction. The effect of groove length ratio on the unbalance response is shown in Figure 22. According to the results, generally, the unbalance response in both directions of x and y decrease with respect to groove length ratio. The most significant improvement exists The effect of groove length ratio on the unbalance response is shown in Figure 22. According to the results, generally, the unbalance response in both directions of x and y decrease with respect to groove length ratio. The most significant improvement exists when the groove length ratio is 1.0, groove depth is 10 µm and rotary speed is 20,000 rpm. The effect of the groove texture on unbalance response changes from negative to positive, indicating that a larger groove length ratio can be beneficial for the dynamics of the waterlubricated bearing-rotor system.

Conclusions
(1) The groove texture can enhance the direct stiffness along the load direction and weaken the stiffness in the orthogonal direction , especially at greater rotational speed, smaller eccentricity ratio, and a shallower and longer groove. When a bearing with a groove depth of 10 µm, a groove length of 1 and width of 1.1 degrees is working at 20,000 rpm and a 0.1 eccentricity ratio, the stiffness can be improved by 98.4%. With the increased width of the groove, the dynamic coefficients show an increase-decrease-increase trend or decrease-increase-decrease trend; with increased eccentricity ratio and rotary speed, the groove texture has little effect on the direct damping coefficients but improves the coupled damping coefficients; for most of the cases, the damping coefficients can be improved with increased groove depth; with the increased width of the groove, the damping coefficients increase slightly. These complex variations may be attributed to the combination of cavitation and flow vortex.
(2) The groove texture can play a more obvious positive effect on the stability of a rotor supported by WLHBs with shallow and longer grooves at a greater rotary speed and smaller eccentricity ratio; there exists an optimum groove width to achieve the largest critical mass. When the length ratio is 1, the stability is improved by about 20.5% with groove widths of 0.1 degrees and 1.1 degrees at a rotary speed of 20,000 rpm. The main reason for this phenomenon is due to the improvement of direct stiffness and direct damping coefficient in the direction of the applied load, as well as the weakening of whirling motion due to the existence of a groove texture with proper parameter settings. (3) As for the unbalance response, the groove texture is beneficial for the dynamics of the rotor at a small eccentricity ratio and great rotary speed. A shallower and longer groove can decrease the unbalance response more apparently; there exists an optimal groove width to decrease the amplitude of the unbalance response. The most significant decrease in groove texture on unbalance response exists at a rotary speed of 20,000 rpm, groove depth of 10 µm, groove length ratio of 1, groove width of 0.1 degree and 1.1 degrees. However, it should be noted that the improvement of the groove texture on the unbalance response of the bearing-rotor system is not as sig-

Conclusions
(1) The groove texture can enhance the direct stiffness along the load direction k yy and weaken the stiffness in the orthogonal direction k xx , especially at greater rotational speed, smaller eccentricity ratio, and a shallower and longer groove. When a bearing with a groove depth of 10 µm, a groove length of 1 and width of 1.1 degrees is working at 20,000 rpm and a 0.1 eccentricity ratio, the stiffness k yy can be improved by 98.4%. With the increased width of the groove, the dynamic coefficients show an increase-decrease-increase trend or decrease-increase-decrease trend; with increased eccentricity ratio and rotary speed, the groove texture has little effect on the direct damping coefficients but improves the coupled damping coefficients; for most of the cases, the damping coefficients can be improved with increased groove depth; with the increased width of the groove, the damping coefficients increase slightly. These complex variations may be attributed to the combination of cavitation and flow vortex. (2) The groove texture can play a more obvious positive effect on the stability of a rotor supported by WLHBs with shallow and longer grooves at a greater rotary speed and smaller eccentricity ratio; there exists an optimum groove width to achieve the largest critical mass. When the length ratio is 1, the stability is improved by about 20.5% with groove widths of 0.1 degrees and 1.1 degrees at a rotary speed of 20,000 rpm. The main reason for this phenomenon is due to the improvement of direct stiffness and direct damping coefficient in the direction of the applied load, as well as the weakening of whirling motion due to the existence of a groove texture with proper parameter settings. (3) As for the unbalance response, the groove texture is beneficial for the dynamics of the rotor at a small eccentricity ratio and great rotary speed. A shallower and longer groove can decrease the unbalance response more apparently; there exists an optimal groove width to decrease the amplitude of the unbalance response. The most significant decrease in groove texture on unbalance response exists at a rotary speed of 20,000 rpm, groove depth of 10 µm, groove length ratio of 1, groove width of 0.1 degree and 1.1 degrees. However, it should be noted that the improvement of the groove texture on the unbalance response of the bearing-rotor system is not as significant as that on the dynamic characteristics of the bearing. As a result, in actual application, considering the manufacturing cost of texture, the effects of texture should be comprehensively considered. Data Availability Statement: The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.