dynamic characteristics of high-speed water- lubricated
TRANSCRIPT
![Page 1: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/1.jpg)
Dynamic Characteristics of High-speed Water-lubricated Spiral Groove Thrust Bearing based onTurbulent Cavitating Flow Lubrication ModelXiaohui Lin ( [email protected] )
Southeast University https://orcid.org/0000-0001-7849-5734shun Wang
Southeast Universityshuyun Jiang
Southeast Universityshaowen Zhang
Southeast University
Original Article
Keywords: spiral groove thrust bearings, water lubrication, turbulent lubrication, cavitating οΏ½ow, dynamiccharacteristics
Posted Date: May 29th, 2020
DOI: https://doi.org/10.21203/rs.3.rs-31639/v1
License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License
Version of Record: A version of this preprint was published at Chinese Journal of Mechanical Engineeringon February 22nd, 2022. See the published version at https://doi.org/10.1186/s10033-021-00671-3.
![Page 2: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/2.jpg)
Title page
Dynamic Characteristics of High-speed Water-lubricated Spiral Groove Thrust Bearing based on Turbulent Cavitating Flow Lubrication Model
Xiaohui Lin, born in 1960, is currently an associate professor at School of Mechanical Engineering, Southeast
University, China.He received his bachelor degree from Wuhan University of Technology, China, in 1982. His
research interests include two-phase flow transport theory. born in 1960,Associate Professor of School of Mechanical
Engineering, Southeast University.
TelοΌ13851760077οΌE-mailοΌ[email protected]
Shun WangοΌborn in 1995, is currently a master at Southeast University, China.
Tel: 18251985252; E-mail: [email protected]
Shuyun Jiang, born in 1966, is currently a PhD supervisor at Southeast University, China.
Tel: 13605161056; E-mail: [email protected]
Shaowen Zhang, is currently PhD at Southeast University, China.
TelοΌ15051867239οΌE-mailοΌ[email protected]
Corresponding authorοΌXiaohui Lin E-mailοΌ[email protected]
![Page 3: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/3.jpg)
ORIGINAL ARTICLE
Dynamic Characteristics of High-speed Water-lubricated Spiral Groove Thrust
Bearing based on Turbulent Cavitating Flow Lubrication Model
Xiaohui Lin1, Shun Wang, Shuyun Jiang, Shaowen Zhang
School of Mechanical Engineering Southeast University, Nanjing 211189, China
Abstract: In this paper, a turbulent cavitating flow lubrication model based on two-
phase fluid and population balance equation of bubbles was established to analyze
dynamic characteristics of high speed water-lubricated spiral groove thrust
bearings(SGTB). Stiffness and the damping coefficients of the SGTB were calculated
using the perturbation pressure equations. An experimental apparatus was developed to
verify the theoretical model. Simulating and experimental results show that the small-
sized bubbles tend to generate in the turbulent cavitating flow when at a high rotary
speed, and the bubbles mainly locate at the edges of the spiral groove. The simulating
results also show that the direct stiffness coefficients are increased due to cavitation
effect, and cross stiffness coefficients and damping coefficients are hardly affected by
the cavitation effect. Turbulent effect on the dynamic characteristics of SGTB is much
stronger than the cavitating effect.
Keywords: spiral groove thrust bearings; water lubrication; turbulent lubrication;
cavitating flow; dynamic characteristics
1. Introduction
Water-lubricated bearing has been widely applied due to its low viscosity, low
temperature rise, no pollution and etc. Spiral groove thrust bearing (SGTB) has
advantages in hydrodynamic effect and stability. Nowadays, the water-lubricated spiral
groove thrust bearing has been used in high-speed spindle. However, since inception
cavitation number of water is larger than that of oil, cavitation phenomenon in water-
lubricated bearings is more serious than the oil-lubricated bearing at high speed, and
the water-lubricated bearing tends to fall into turbulent fluid condition when rotary
speed exceeds a certain value. Hence, both cavitating effect and turbulent effect should
be considered to model the high-speed water-lubricated bearing.
The cavitating and turbulent effects of water or oil bearings with different
![Page 4: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/4.jpg)
configurations have been investigated in previous studies. In the early studies on
bearing cavitation, the different boundary conditions for lubrication equation including
the Swift-Stiebe condition [1], the JFO condition [2], the Elrod condition [3] have been
adopted to describe the cavitation effect of the lubricating film. With those conditions,
solution of a two-phase flow in cavitated region can be avoided. However, the cavitating
flow is not involved in those models. To overcome this problem, the cavitation lubricant
is treated as two-phase mixed fluid. The two-phase mixed fluid models have been
developed, those models can be divided into three types based on calculation
approaches of gas phase volume fractionοΌ(1) the model based on the R-P equation [4,5]
(2) the model based on gas solubility and surface tension of bubble [4,6-9] (3) the model
based on transport equation of the gas phase volume fraction [10-15]. In addition, with
the rapid development of computational fluid dynamics (CFD) technology, some
researchers have proposed performance analysis of bearings with the cavitation using
the CFD software [7,11,16-19]. However, the interfacial effect of the two-phase flow is
not included in the above lubrication models of cavitating flow. Actually, the
momentum and energy exchange between the two phases produces at the interface, so
the interface effect is a very important feature for two-phase flow, which can not be
ignored when modeling the high-speed water-lubricated bearing. In addition, the size
distribution of the bubbles cannot be obtained using the present mixed fluid model or
the CFD model. In our previous study, lubrication models considering the interface
effect of two-phase flow have been proposed[20-23]. However, those models are
established under heat isolation, and the bubble population equilibrium equation is
solved directly. Since the bubble size probability density function is an unknown
quantity, so it is difficult to determine the bubble size probability density function at the
initial moment and at the interface. Furthermore, this kind of algorithm tends to unstable.
This study aims to establish a lubrication model for the water-lubricated spiral groove
thrust bearings including both two-phase interfacial and fluid turbulent and bearing heat
conduction effects. A generalized Reynolds equation considering the cavitation
interface and turbulence, and an energy equation considering bearing heat conduction
are derived. A population balance equation considering the breakage and coalescence
of bubbles is introduced to describe the evolution of the size distribution of bubbles. A
discrete interval method is employed to solve the bubble population balance equation.
The influence of cavitation and turbulence on the dynamic characteristics of the SGTB
![Page 5: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/5.jpg)
is analyzed using the established model. The size distribution of cavitation bubbles and
stiffness coefficient of water film are measured using a self-developed experimental
device. A comparison of the simulating results with experimental ones is given to verify
the proposed model.
2.Theoretical model2.1 The generalized Reynolds equation based on the two-phase flow theoryThis paper studies the water-lubricated SGTB with the turbulent cavitating flow film.
The classical Reynolds equation is not suitable for describing the bearing lubrication
state, a Reynolds equation which contains turbulent and two-phase interfacial effects is
needed, so the following basic assumptions are made for the modeling:
The bubbles in the cavitation flow be regarded as spheres.
The flow field is a small perturbation of the turbulent Couette flow [24].
The turbulent shear stress is defined by the law of wall, Reichardt [25] empirical
formula is used to define the eddy viscosity coefficient.
The average viscosity along the film thickness is adopted in generalized turbulent
Reynolds equation.
The spiral area is an irregular area in the polar coordinate system, a spiral
coordinates system (π, π) are adopted to improve numerical calculation accuracy. The
polar coordinates (π, π) are converted into the spiral coordinates system (π, π) as
follows: π = ππ = π β π(π)π(π) = 1cotπ½ ππ ( πππ) (1)
The spiral patterns before and after the coordinates conversion are shown in Fig. 2.
![Page 6: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/6.jpg)
Fig. 1 Structure schematic of inward-pumping SGTB
Fig. 2 Coordinates transformation of spiral patterns
A generalized turbulent Reynolds equation with cavitation interface effects and
inertial effects in the (π, π) coordinates is derived based on two-phase flow theory[26]
πππ ( β3ποΏ½Μ οΏ½π€ππ ππΆπ€πππ ) + πππ ((π2πβ²2(π)ππ + 1ππ) β3οΏ½Μ οΏ½π€π ππΆπ€πππ )β πππ ( β3ποΏ½Μ οΏ½π€ππ πβ²(π) ππΆπ€πππ ) β πππ ( β3ποΏ½Μ οΏ½π€ππ πβ²(π) ππΆπ€πππ )= πππ (βππΆπ€2 ) + πππ ( β3οΏ½Μ οΏ½π€πππ(π)) +πππ (β3ππ€π’2πΆπ€οΏ½Μ οΏ½π€ππ ) β πππ (πβ²(π) β3ππ€οΏ½Μ οΏ½2πΆπ€οΏ½Μ οΏ½π€ππ ) (2)
where οΏ½Μ οΏ½π€ = 1β β« ππ€β0 ππ§οΏ½Μ οΏ½ = 1β β« π’β0 π’ππ§ (3)
![Page 7: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/7.jpg)
1ππ = β« β« 1ππ (1 β ππππ) (12 β π§ββ²) ππ§ββ²ππ§βπ§β0
101ππ = β« β« 1ππ (12β π§ββ²) ππ§ββ²π§β
01
0 ππ§βππ = 0.4 [π§ββπβ β 10.7π‘β (π§ββπβ10.7)] + 1ππ = 0.2π§ββπβπ‘β2 (π§ββπβ10.7)π§β = π§β
(4)
βπβ = βππ€β ππππ€ (5)
The following expression can be given for the logarithmic spirals πβ²(π) = 1πcotπ½ (6)
The circumferential velocity distribution of the water film is written asπ’ = π2 + (βπβ)2π π πββ πΊ1(π§β) β 1πΆπ€ β2οΏ½Μ οΏ½π€ (1π π(πΆπ€π)ππ βπ(π))πΊ2(π§β) (7)
where πΊ1(π§β) = β« 1ππ ππ§ββ²π§β12πΊ2(π§β) = β« 1ππ (1 β ππππ) ππ§ββ²π§β
0(8)
πΆπ€ = 1 β πΆπ (9)πΆπ = π6 β« π3πππ(π,π)β0 ππ1+π6 β« π3πππ(π,π)β0 ππ (10)
where, πππ(π, π) is the equilibrium distribution function of the cavitation bubble
diameter; πππ(π, π)ππ represents the number of bubbles with diameter between (π, π + ππ) per unit volume of water at spatial coordinate Ξ· under the equilibrium state.
The first term on the right-hand side of Eq. (2) is the hydrodynamic effect, the second
term on the right side of Eq. (2) is the interface hydrodynamic effect, the third and
fourth terms on the right side of Eq. (2) is the inertia hydrodynamic effect. The interface
momentum transfer function includes the following three parts:π(π) = π1(π) + π2(π) + π3(π) (11)
where, π1(π) is the interface momentum transfer term due to mass transfer, π2(π) is the interface momentum transfer term due to viscous drag, π3(π) is the interface
![Page 8: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/8.jpg)
momentum transfer term due to surface tension. Because the film thickness is thin, the
liquid velocity in π1(π),π2(π) are approximated by the average
velocity. π1(π),π2(π) and π3(π) are written as:π1 = β[ππ€ π3 β« π2β0 πππ(π, π)ππ] οΏ½Μ οΏ½(οΏ½Μ οΏ½ β π’π)π2 = [π3 πΆπ·π ππ€ β« π2β0 πππ(π, π)ππ] (οΏ½Μ οΏ½ β π’π)2π3 = 2ππ β« ππππ(π, π)β0 ππ (12)
where πΆπ·π = 24π πAssume that the range of bubble diameter is divided into M smaller intervals, the
diameter distribution of the bubbles is uniform in a smaller interval. So the equilibrium
distribution function of bubble diameter in the kth smaller diameter interval (ππ, ππ+1) can be expressed as. πππ(π, π) β ππππ (π)πΏ(π β π₯π) (13)
where ππππ (π) = β« πππππ+1ππ (π, π)ππ (14)π₯π = 12 (ππ + ππ+1) (15)
where ππππ (π) is the number of bubbles in the kth diameter interval (ππ, ππ+1). And
the integral in Eq. (12) can be approximated as
β« π2β0 πππ(π, π)ππ β β π₯π2ππππ (π)ππ=1β« πβ0 πππ(π, π)ππ β β π₯πππππ (π)ππ=1 (16)
The boundary conditions of Eq. (2) are
![Page 9: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/9.jpg)
π|π=πππ’π‘ = 0π|π=πππ = ππππ(π) = π(π + 2ππ ) (17)
The pressure at the groove-ridge boundary can be solved by using continuous condition
of flow at the groove-ridge boundary. The finite volume at the groove-ridge boundary
is shown in Fig. 3.
Fig.3 Finite volume separated by grooveβridge boundary
In a finite volume at the groove-ridge boundary, the continuous condition of the flow
can be written as β ππππ + β ππππ = 0 (18)
where
πππ =β« (ππβ²(π) ( β3πΆπ€ππ (ππΆπ€πππ β 1πcotπ½ ππΆπ€πππ ) β β3ππ€ ππ’2π ) 1ππ+ πβ2 β β3ππ ( 1ππΆπ€ ππΆπ€πππ β π(π)πΆπ€ ) 1ππ )π2π1
ππ (19)
(π = π΅πΉ, πΉπΆ, π·πΈ, πΈπ΄)πππ = β« π (β β3πΆπ€ππ (ππΆπ€πππ β πβ²(π) ππΆπ€πππ ) 1ππ + β3ππ ππ’2π 1ππ)π2π1 ππ (20)(π=π΄π΅, πΆπ·)2.2 Turbulent cavitating flow energy equation
The temperature field of the cavitating flow water film can be obtained by solving
the two-phase flow energy equation under high-speed conditions, and the following
assumptions are adopted to simplify the energy equation.
![Page 10: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/10.jpg)
(1) The circumferential convection term is only considered in the energy equation,
because the circumferential velocity is much larger than the radial velocity.
(2) Turbulent pulsation heat conduction term is ignored.
(3) The viscous dissipation term along the film thickness is only considered in the
energy equation.
According to the above assumption, the energy equation in the (π, π) coordinates
in the spiral coordinate system can be simplified toππ€ππ£ (π’π π(πΆπ€π)ππ ) = πΆπ€ππ€ π2πππ§2 + πΆπ€(ππ‘(π§)) (ππ’ππ§) + πΈπππ‘ (21)
where
(ππ‘(π§)) = (βπβ)2οΏ½Μ οΏ½π€ππ πβ + 1πΆπ€ (1π π(πΆπ€π)ππ βπ(π)) (π§ β β2) (22)
ππ’ππ§ = (βπβ)2π πββ ππΊ1(π§β)ππ§β πβ β βπΆπ€οΏ½Μ οΏ½π€ (1π π(πΆπ€π)ππ βπ(π)) ππΊ2(π§β)ππ§β (23)
The interface term πΈintin the energy equation is expressed asπΈπππ‘ β πΆπ·π π(οΏ½Μ οΏ½ β π’π)2οΏ½Μ οΏ½β π₯π2ππππ (π)ππ=1 β ππ€ππ£π(οΏ½Μ οΏ½ β π’π)πβ π₯π2ππππ (π)ππ=1+2ποΏ½Μ οΏ½πβ π₯πππππ (π)ππ=1(24)
The following boundary conditions are adopted for the energy Eq. (21)
(1) At the inlet of the spiral grooveπ(0, π§) = π0 (25)
(2) On the interface of the water film and the thrust diskπ(π, β) = ππ(π, 0) (26)
(3) On the interface of the water film and the stationary ringπ(π, β) = ππ (π, πΏ) (27)
The viscosity temperature relationship is written as ππ€ = ππ€0π(βπ(πβπ0)) (28)
The water film temperature field can be obtained solved by simultaneously solving
the energy equation (21), the heat conduction equation of the stationary ring and the
heat conduction equation of the thrust disk. The surface temperature ππ(π, 0) of the
stationary ring and surface temperatureππ (π, πΏ) of the thrust plate are set as the
boundary condition of the energy equation (21).
2.3 The heat conduction equation in stationary ring
![Page 11: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/11.jpg)
The circumferential and radial heat conduction of the stationary ring are ignored, the
heat conduction equation is simplified as follows
π2πsππ§π 2 = 0 (0 β€ π§π β€ πΏπ ) (29)
The boundary conditions of Eq. (29) are given as follows
1) At the surface between the stationary ring and the fluid (π§π = πΏπ )βππ πππ ππ§π |π§π =πΏπ = πΌπ€(ππ |π§π =πΏπ β ππ)ππ = 1β β« πβ0 ππ§ (30)
2) At the surface between the stationary ring and the ambient (π§π = 0)ππ πππ ππ§π |π§π =0 = πΌπ(ππ |π§π =0 β π0) (31)
The expression of the temperature at the surface of the stationary ring can be
obtained from Eq. (29) and boundary conditions (30), (31). The temperature at the
surface of the stationary ring is written as
ππ |π§π =πΏπ = ππβπ0πΌπ€πΌπβπΌπ€ππ πΏπ β1+ ππ (32)
2.4 The heat conduction equation in thrust diskFor the rotary thrust disk, the heat transfer effects in circumferential and radial
directions are also ignored, the heat conduction equation is simplified as followsππππππ π2ππππ§π2 = π πππππ (0 β€ π§π β€ πΏπ) (33)
The boundary conditions of Eq. (33) are given as follows
1) At the surface between the thrust disk and the fluid (π§π = 0)ππ πππππ§π|π§π=0 = πΌπ€(ππ|π§π=0 β ππ) (34)
2) At the surface between the thrust disk and the ambient (π§π = πΏπ)ππ πππππ§π|π§π=πΏπ = βπΌπ(ππ|π§π=πΏπ β π0) (35)
2.5 The Force Balance Equation of Bubble
The force balance equation of bubbles is expressed as πΉπ΅ + πΉπ· + πΉπΊ = 0 (36)
where
![Page 12: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/12.jpg)
πΉπ΅ = β¨πβ©(ππ€ β ππ)π (37)πΉπΊ = (ππ€ β ππ)β¨πβ© ππ’ππ‘ (38)ππ’ππ‘ = β 1ππ€ πππππ (39)β¨πβ© β π6β π₯π3ππππ (π)ππ=1 (40)
πΉπ· = 12ποΏ½Μ οΏ½π€β¨π πβ©(π’π β οΏ½Μ οΏ½) 1βπΆπ53(1βπΆπ)2 (41)
β¨π πβ© β 12β π₯πππππ (π)ππ=1 (42)
The bubble velocity π’π in the cavitating flow is calculated by the force balance Eq.
(36). Substitute πΉπ΅, πΉπΊ, πΉπ· into Eq. (36), the bubble velocity bu can be expressed as
[27] π’π = οΏ½Μ οΏ½ + (1βπΆπ)2(ππ€βππ)12ποΏ½Μ οΏ½π€β¨π πβ©(1βπΆπ53) β¨πβ© (π + 1ππ€ πππππ) (43)
2.6 The population balance equation of bubbles
2.6.1 The population balance equation
Interface effect is an important phenomenon for cavitating flow, the momentum,
mass and energy transfer occur through the interface between the gas-liquid. The
interfacial area per unit volume liquid is positively correlated with the size distribution
of the bubble, the bubble size distribution can be described by defining a probability
density functionπ(π«, π, π‘), and the internal coordinates π are taken as bubble diameter
in the model, and the π(π«, π, π‘)ππ represents the number of bubbles with between (π, π + ππ) at per volume liquid. Thus, integration of π over bubble diameter results
in the total number of bubbles per volume liquid. The breakage, coalescence are two
main factors affecting bubble size distribution in the cavitating flow, and bubble size
distribution is predicted by the model of the breakage and coalescence processes of
bubbles, this leads to the so-called population balance equation (PBE).
The PBE is a transport equation, the evolution of function π can be described by the
PBE in the spiral coordinates (π, π), the PBE is written as follows
![Page 13: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/13.jpg)
ππ(π, π)ππ‘ + π’π ππππ (π(π, π)) = βπ(π)π(π, π)+β« βπ(π, π)π(π)π(π, π)πmaxπ ππ β π(π, π)β« π(π, π)π(π, π)πmax0
ππ+ π22 β« π((π3 β π3)1 3β , π)(π3 β π3)2 3β π(π, (π3 β π3)1 3β )π(π, π)π
0 ππ + ππ(π) (44)
where, the first term on the right-hand side of Eq. (44) represents breakup sink term of
bubbles with diameter π per unit time, the second term on the right-hand side of Eq.
(44) represents breakup source term of bubbles with diameter π per unit time, the third
term on the right-hand side of Eq. (44) represents coalescence sink term of bubbles with
diameter π per unit time, the fourth term on the right-hand side of Eq. (44) represents
coalescence source term of bubbles with diameter π per unit time, the fifth term on the
right-hand side of Eq. (44) represents bubbles with diameter π source term.
The initial condition for Eq. (44) is givenπ(π, π, 0) = 0 (45)
The breakage frequency π(π) is given as [28]π(π) = π 1π1 3βπ2 3β exp(β ππ 2πππ2 3β π5 3β ) (46)π = πΆππ’ππ (47)
The daughter bubble size redistribution function βπ(π, π)is given as [28]
βπ(π, π) = 4.8π ππ₯π (β4.5 (2πβπ)π2 2) (48)
The coalescence closure is given as [28]π(π, π) = 0.05 π4 (π + π)2[π½π(ππ)2 3β + π½π(ππ)2 3β ]1 2βππ₯π (β (ππ3ππ€ 16πβ )1 2β π1 3β ln(βπ βπβ )ππ2 3β ) (49)
ππ = 14 (1π + 1π)β1 (50)
The value of parameter π½π is taken as 2.48. The expression of source term ππ(π) is
![Page 14: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/14.jpg)
ππ(π) = πΆπβ(ππ£ β π)π(π) (51)
where π(π) = β« (π(π, ππ)ππ(ππ))ππππ0 (52)πΆπ = ( 23ππ€)1 2β(53)
The redistribution function of bubble diameter π is given as π(ππ, π) = π1ππ ππ₯π (β π2(πβππ)2πmax2 ) (54)
The diameter distribution density function of the gas nucleus is given by
experimental data. [29] ππ(ππ) β (πΌ+1)πππππ₯ππ ( πππππππ₯)πΌ (55)
where, the value of parameter Ξ± is taken as -10/3.
3. Calculation of dynamic characteristics of the SGTB
3.1 Perturbation generalized Reynolds equations of cavitation water film
As shown in Fig.4, the origin of Cartesian coordinate is located at the center of the
stationary ring. The possible independent motions of the thrust disk are axial movement
along z axis and angular movements around x and y axes. The perturbed displacements
and velocities of the thrust disk with respect to the steady equilibrium position are
defined as (π₯π§, π₯ππ₯, π₯ππ¦) and (π₯οΏ½ΜοΏ½, π₯οΏ½ΜοΏ½π₯, π₯οΏ½ΜοΏ½π¦) . The assume initial position of the
thrust disk is (β0, ππ₯0, ππ¦0), the film thickness in quasi-equilibrium can be expressed
as β = β0 + π₯π§ + πcosππ₯ππ¦ β πsinππ₯ππ₯ (56)
The perturbations will influence the pressure distribution in the water film through
Reynolds equation. The high order terms are ignored, the transient pressure of water
film can be expressed as π = π0 + ππ§π₯π§ + πππ₯π₯ππ₯ + πππ¦π₯ππ¦ + ποΏ½ΜοΏ½π₯οΏ½ΜοΏ½ + ποΏ½ΜοΏ½π₯π₯οΏ½ΜοΏ½π₯ + ποΏ½ΜοΏ½π¦π₯οΏ½ΜοΏ½π¦ (57)
![Page 15: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/15.jpg)
Substituting Eqs. (56) and (57) into Eq. (2), the zeroth and first-order perturbation
generalized Reynolds equations in the (π β π) coordinates system are obtained for the
SGTB respectively,
π ππ¦(){ πΆπ€ππ§πΆπ€πππ₯πΆπ€πππ¦πΆπ€ποΏ½ΜοΏ½πΆπ€ποΏ½ΜοΏ½π₯πΆπ€ποΏ½ΜοΏ½π¦}
={ ππ§πππ₯πππ¦ποΏ½ΜοΏ½ποΏ½ΜοΏ½π₯ποΏ½ΜοΏ½π¦}
(58)
where, the operator π ππ¦() is defined asπππ ( β03οΏ½Μ οΏ½π€π (π2πβ²2(π)ππ +1ππ) πππ)+
πππ ( πβ03οΏ½Μ οΏ½π€ππ πππ)β πππ ( β03οΏ½Μ οΏ½π€ππ ππβ²(π) πππ) β πππ ( β03οΏ½Μ οΏ½π€ππ ππβ²(π) πππ)ππ(π = π§, ππ₯, ππ¦, οΏ½ΜοΏ½, ποΏ½ΜοΏ½, ποΏ½ΜοΏ½)can be found in Appendix A.
The stiffness and damping coefficients in the (π β π) coordinates system can be
found by the following integrals
[ πΎπ§π§ πΎπ§ππ₯ βπΎπ§ππ¦πΎππ₯π§ πΎππ₯ππ₯ πΎππ₯ππ¦πΎππ¦π§ πΎππ¦ππ₯ πΎππ¦ππ¦] =β« β« [ β1πΊ1(π, π)βπΊ2(π, π)] {ππ§ ββπππ₯ ββπππ¦}ππππππππ’π‘πππ
2π0
(59)
[ πΆπ§π§ πΆπ§ππ₯ βπΆπ§ππ¦πΆππ₯π§ πΆππ₯ππ₯ πΆππ₯ππ¦πΆππ¦π§ πΆππ¦ππ₯ πΆππ¦ππ¦] =β« β« [ β1πΊ1(π, π)βπΊ2(π, π)] {ποΏ½ΜοΏ½ ββποΏ½ΜοΏ½π₯ ββποΏ½ΜοΏ½π¦}ππππππππ’π‘πππ
2π0
(60)
where πΊ1(π, π) = πsin(π + π(π))πΊ2(π, π) = πcos(π + π(π)) (61)
![Page 16: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/16.jpg)
Fig. 4 Motion of a SGTB
4. Numerical calculation method of theoretical model
The generalized Reynolds equation (2), energy equation (21), and PBE (44) are the
main three governing equations for the bearing. The pressure field, temperature field,
and bubble diameter distribution are obtained by solving these three equations
respectively, and pressure field; temperature field and bubble diameter distribution are
coupled to each other, and an iterative algorithm needs to be employed to solve these
distributions. Eq. (2) and Eq. (21) are discretized using finite difference, the pressure
field is obtained by solving discretization equations of Eq. (2) with the SOR iterative
method, the temperature field is obtained by solving discretization equations of Eq. (21)
with the stepping method.
4.1 Calculation of parameterππβThe shear stress parameter βπβ of the Couette flow is included in turbulent velocity
distribution Eq. (2), The governing equation for parameter βπβ is written as [24](βπβ)2β« ππ§βππ(π§β,βπβ)10 β π πβββ = 0 (62)
Eq. (62) is a nonlinear equation with respect to βπβ, the parameter βπβ can be solved
using the Newton iteration method. The Newton iteration formula is expressed asβπβ
![Page 17: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/17.jpg)
(βπβ)(π+1) = (βπβ)(π) β πΉ(βπβ(π))πΉβ²(βπβ(π)) (63)
where πΉ(βπβ) = (βπβ)2β« ππ§βππ(π§β,βπβ)10 β π πβββ (64)
4.2 Numerical solution of bubble PBEA discrete interval method [30] is employed to solve the bubble population balance
Eq. (44). The procedure of the discrete interval method is as follows:
(1) Assume that the entire bubble diameter range is divided into M smaller intervals.
The bubble diameter is approximated uniform distribution at kth representative
diameter range(ππ, ππ+1), then the π(π, π) can be expressed by the bubble number
density ππ(π) at kth representative diameter range (ππ, ππ+1). Thus,π(π, π) β ππ(π)πΏ(π β π₯π) (65)
where ππ(π) = β« π(π, π)ππ+1ππ ππ (66)
(2) The integrating the continuous Eq. (44) over a smaller range (ππ, ππ+1). Thus,πππ‘ β« π(π, π)ππ+1ππ ππ + π’π ππππ (β« π(π, π)ππ+1ππ ππ) = ββ« π(π, π)π(π, π)ππ+1ππ ππ+β« ππππ+1ππ β« βπ(π, π)π(π, π)π(π, π)ππππ₯π ππββ« π(π, π)ππ+1ππ ππ β« π(π, π)π(π, π)ππππ₯0
ππ+ 12β« π2ππππ+1ππ β« π((π3βπ3)1 3β ,π)(π3βπ3)2 3β π(π, (π3 β π3)1 3β )π(π, π)π
0 ππ+πΆπβ(ππ£ β π)β« π(π)ππ+1ππ ππ
(67)
(3) Since coalescence or breakage of bubbles of these diameters results in the formation
of new bubbles whose diameters do not match with any of the node diameters, the
new bubble diameter will be assigned to the adjoining node by introducing weight
factor πΎ, the weight factor πΎ can be determined by the total mass conservation of
the bubble and the number of bubbles conservation. The expression is written asπΎππ₯π3 + πΎπ+1π₯π+13 = π3πΎπ + πΎπ+1 = 1 (68)
Substituting Eq. (65) into Eq. (67), and the integral term at the right end of Eq.
(67) is reconstructed using the mean value theorem on breakage frequency and
![Page 18: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/18.jpg)
coalescence frequency, so Eq. (67) can be closed, a closed set of equations with
regard to ππ(π) can be written as πππ(π)ππ‘ + π’π ππππ (ππ(π)) = βπ(π, π₯π)ππ(π)+β πΌπ,πππ=π π(π, π₯π)ππ(π) β ππ(π)β π(π, π₯π , π₯π)ππ(π)ππ=1+β (1 β 12 πΏπ,π) πΎπ,ππ π(π, π₯π , π₯π)ππ(π)πβ₯π π,ππ₯πβ1β€(π₯π3+π₯π3)1 3β β€π₯π+1 ππ(π)+πΆπβ(ππ£ β π)β« π(π)ππ+1ππ ππ(69)
where
πΎπ,ππ = { π₯π+13 β(π₯π3+π₯π3)(π₯π+13 βπ₯π3) βββββββββπ₯π3 β€ (π₯π3 + π₯π3) β€ π₯π+13(π₯π3+π₯π3)βπ₯π3(π₯π+13 βπ₯π3) ββββββββββπ₯πβ13 β€ (π₯π3 + π₯π3) β€ π₯π30ββββββββββββββββββββββπππ π
(70)
π = 1,2, . . . . . . π β 2 , π3 = ππ3 + ππ3 πΎπ,ππβ1 = {
(π₯π3+π₯π3)π₯πβ13 ββββββββββββββββββββββββπ₯π-13 β€ (π₯π3 + π₯π3) β€ ππ3(π₯π3+π₯π3)βπ₯πβ23(π₯πβ13 βπ₯πβ23 ) ββββββββββπ₯πβ23 β€ (π₯π3 + π₯π3) β€ π₯πβ130ββββββββββββββββββββββββββββπππ πββ
(71)
πΏπ,π = {1βββββββββπ = π0ββββββββπ β π (72)
Coefficient πΌπ,π is given asπΌπ,π = β« (π3βπ₯πβ13 )(π₯π3βπ₯πβ13 )βπ(π₯π, π)πππ₯ππ₯πβ1 +β« (π₯π+13 βπ3)(π₯π+13 βπ₯π3)βπ(π₯π, π)πππ₯π+1π₯π (73)
Eq. (69) is a set of first-order differential equations. Suppose that at a fixed time t, a
set of nonlinear equations with regard to ππ,ππ‘ can be obtained by discretizing the
second term on the left side of Eq. (69). The discrete equation at node π is written as
follows οΏ½ΜοΏ½(ππ,ππ‘ ) = 0 (74)
where, the operator οΏ½ΜοΏ½(ππ,π)are defined as
![Page 19: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/19.jpg)
οΏ½ΜοΏ½(ππ,ππ‘ ) = π’πππ ππ,ππ‘ βππ,πβ1π‘π₯π + π(ππ, π₯π)ππ,ππ‘ββ πΌπ,ππ(ππ, π₯π)ππ,ππ‘ππ=π + ππ,ππ‘ β π(ππ, π₯π , π₯π)ππ,ππ‘ππ=1ββ (1 β 12 πΏπ,π) πΎπ,ππ π(ππ , π₯π , π₯π)ππ,ππ‘πβ₯π π,ππ₯πβ1β€π₯π+π₯πβ€π₯π+1 ππ,ππ‘ β πΆπβππ£ β πβ« π(π)ππ+1ππ ππ(75)
Eq. (74) can be solved using the Newton-SOR iterative method. In addition, since βπ(π₯π, π) is a given function, the coefficient πΌπ,π can be calculated using the
Gaussian integration method, and the value of the coefficient πΌπ,π is stored in a matrix,
it will be called directly in the iterative calculation process, and the calculation time
can be greatly saved by using the iterative process.
The term of the time derivative in the transient equation (69) is specially processed
with finite difference, and a weighted averaging method is adopted for the time-
differencing. It can be written as [31]ππ,ππ‘+1βππ,ππ‘π₯π‘ + ((1 β ππ‘)οΏ½ΜοΏ½(ππ,ππ‘ ) + ποΏ½ΜοΏ½(ππ,ππ‘+1)) = 0 (76)
where, ππ,ππ‘+1 = ππ,π(π‘ + π₯π‘), weight factor 0 β€ ππ‘ β€ 1.
A set of nonlinear equations with regard to ππ,ππ‘+1 can be rewritten asπΉπ(ππ,ππ‘+1) = ππ,ππ‘+1 + π₯π‘ππ‘οΏ½ΜοΏ½(ππ,ππ‘+1) + π₯π‘(1 β ππ‘)οΏ½ΜοΏ½(ππ,ππ‘ ) β ππ,ππ‘ = 0β(π = 1, . . . ππ) (77)
The nonlinear equation (77) is also solved using the Newton-SOR iterative method.
And downhill condition is adopted in numerical iterative procedure to improve
convergence. The downhill condition is expressed asβπΉπ(ππ,ππ‘+1(π+1))β β€ βπΉπ(ππ,ππ‘+1(π))β (78)
where π is iteration number. The overall algorithm flow diagram is given by Fig.
5.
![Page 20: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/20.jpg)
![Page 21: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/21.jpg)
Fig. 5 Flow diagram of overall algorithm
5. Results and Discussion
5.1 Experimental verification of theoretical modelsTable 1 list the related parameters of the water-lubricated SGTB for the numerical
calculation. In order to verify this theoretical model, the simulating results of bubble
distribution and axial stiffness coefficient are compared with experimental ones from
our previous test study in Refs.[22,23].
Table 1 The parameters of the water-lubricated SGTB
Item Value
Inner radius πππ(ππ) 7.5
Outer radius πππ’π‘(ππ) 20
Base circle radius ππ(ππ) 12
Water viscosity π(ππ. π ) 0.001(20Β°)
Groove number π 12
Groove depth βπ(ππ) 40
Spiral angles π½(β) 20
Pressure of supply water πππ(πππ) 0.1
Specific heat at constant volume of water ππ£β(π½ ππβ . πΆβ) 4200
Density of water ππ€ β(ππ π3β ) 1000
Bearing capacity π(π) 20
Surface tension of bubble π(π πβ ) 0.3
Groove-to-land ratio ππ 0.5
Groove-to-dam ratio ππ 0.6
1οΌBubble distribution
Fig.6 shows the photographs of bubble distributions taken during the experiments
under the rotational speeds of 15000 rpm and 18000 rpm.
![Page 22: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/22.jpg)
(a) 15000rpm
(b) 18000rpm
Fig.6 Photographs of bubble distributions under both rotational speeds: (a)
15000 rpm, (b) 18000 rpmThe bubble size distribution experimental data are counted from the photos. A
comparison of the theoretical results of the bubble size distribution with the
experimental statistics results is shown in Fig.7. It can be seen that, for a specified
rotational speed of 15000 rpm, the results predicted by the model are in good agreement
with the experimental statistics results over the entire bubble size range; while, for a
specified rotational speed of 18000 rpm, the results predicted by the model are in good
agreement with the experimental statistics results in the large size bubble range, but the
deviation between the theoretical calculation and the experimental statistics results is
large in the small size bubble range, and the experimental statistical values are all larger
than the theoretical ones in the small size bubble range. This may be explained by the
fact that the number of bubbles is excessively counted in the small-sized bubble range
from the experimental photos, since there exists an error unavoidably when the edge
contours of small size bubbles were identified from the photos.
![Page 23: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/23.jpg)
0 2 4 6 8 10
0.0
0.1
0.2
0.3
0.4
Experiment
Theory
Dimensionsless bubble diameter *
Pro
ba
bil
ity
den
sity
of
bu
bb
le v
olu
me,
PV
(a) 15000rpm
0 2 4 6 8 10
0.0
0.1
0.2
0.3
0.4
0.5
Experiment
Theory
Dimensionsless bubble diameter *
Pro
ba
bil
ity
den
sity
of
bu
bb
le v
olu
me,
PV
(b) 18000rpm
Fig. 7 Comparison of the theoretical probability density distributions with the
experimental ones: (a)15000rpmοΌ(b)18000rpm
(π½=20β, βπ = 15ππ, βπ = 40ππ)
2οΌAxial stiffness coefficient
A comparison of the predicting axial stiffness coefficients with the experimental
ones is shown in Fig. 8. It can be seen that the predicting stiffness coefficients with
axial load are generally in agreement with the experimental ones, but the experimental
values are larger than the predicting ones in the entire range of applied load. This may
be explained by the fact that the actual water flowrate of the bearing is less than the
![Page 24: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/24.jpg)
theoretical one, because the water is tended to be prevented from entering the clearance
of inward-pumping SGTB due to the centrifugal effect.
30 60 90 120 150 180 210
5
10
15
20
25
W (N)
Kzz
(N
/m
)
Experiment
Theory
Fig. 8 Comparison of the theoretical stiffness coefficients with the
experimental ones at 15000rpm
5.2 Influence of cavitation effect on dynamic characteristics
1οΌStiffness coefficient
Fig. 9 shows the influence of cavitation effect on the stiffness coefficients of the
SGTB with different spiral angles. The results show the direct stiffness
coefficients(πΎπ§π§β , πΎππ₯ππ₯β , πΎππ¦ππ¦β )of the cavitating flow are larger than those with non-
cavitating flow when the spiral angle is less than 50 degrees. This may be explained by
the fact that the bubbles cause the interfacial hydrodynamic effect for the the cavitating
flow. In addition, the relative deviation between the cavitating flow model and the non-
cavitating flow model decreases with the increasing of spiral angle, so the interface
hydrodynamic effect of the bubbles decreases with the spiral angle. It means that the
influence of the cavitation effect on the water film stiffness is weakened with the
increasing of spiral angle. The cross-coupled stiffness coefficients (πΎππ₯ππ¦β , πΎππ¦ππ₯β )are
![Page 25: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/25.jpg)
of opposite sign, and πΎππ₯ππ¦β and πΎππ¦ππ₯β of cavitating flow with the angle generally
coincide with those of non-cavitating flow. It shows that the influence of cavitation
effect on cross-coupled stiffness coefficients (πΎππ₯ππ¦β , πΎππ¦ππ₯β )are negligible. In addition,
the cross-coupled stiffness (πΎππ₯π§β , πΎπ§ππ₯β , πΎππ¦π§β , πΎπ§ππ¦β ) with different spiral angles are
approximately equal to zero for the cavitating flow or the non-cavitating flow.
5 10 15 20 25 30 35 40 45 50 55 60 65-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8 cavitating turbulent flow
non-cavitating turbulent flow
K*
zz
K*yz
Spiral angle (Β°)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
5 10 15 20 25 30 35 40 45 50 55 60 65
-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
K*zx
K*xx
cavitating turbulent flow
non-cavitaing turbulent flow
Spiral angle (Β°)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(a) (b)
5 10 15 20 25 30 35 40 45 50 55 60 65-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
K*zy
K*yy
Spiral angle (Β°)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
cavitaing turbulent flow
non-cavitating turbulent flow
5 10 15 20 25 30 35 40 45 50 55 60 65-1.5
-1.0
-0.5
0.0
0.5
1.0
1.5
K*xz
K*yx
K*xy
Spiral angle (Β°)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
cavitaing turbulent flow
non-cavitating turbulent flow
(c) (d)
Fig. 9 Influence of cavitation effect on the stiffness coefficient for different
spiral angles
(π = 18000πππ, βπ = 15ππ, βπ = 40ππ)
Fig. 10 shows the influence of cavitation effect on the stiffness coefficient of the
SGTB with different groove depth. The result also show the direct stiffness coefficients (πΎπ§π§β , πΎππ₯ππ₯β , πΎππ¦ππ¦β )of cavitating flow are larger than those of non-cavitating flow at a
![Page 26: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/26.jpg)
specific groove depth. This is also due to the fact that the bubbles cause the interfacial
hydrodynamic effect for the cavitating flow. However, the differences of direct stiffness
coefficients (πΎπ§π§β , πΎππ₯ππ₯β , πΎππ¦ππ¦β ) between the cavitating flow and the non-cavitating
flow increase first and then decrease with the groove depth. The step hydrodynamic
effect becomes weak when the groove depth exceeds a given value, so the velocity
difference between fluid and bubble reduces and the bubble density increases due to the
weakening of hydrodynamic effect. The interface hydrodynamic effect is positively
correlated with the velocity difference and the bubble density. Therefore, the interface
hydrodynamic effect changes with the groove depth, meaning that the influence of
cavitation effect on stiffness depends on the groove depth.
In addition, the variation of the cross-coupled stiffness coefficients(πΎππ₯ππ¦β , πΎππ¦ππ₯β )
with the groove depth is identical with that with the spiral angle. Similarly, it also shows
that the influence of cavitation effect on cross-coupled angular stiffness coefficients (πΎππ₯ππ¦β , πΎππ¦ππ₯β ) are negligible for all the groove depths. The cross-coupled stiffness
coefficients (πΎππ₯π§β , πΎπ§ππ₯β , πΎππ¦π§β , πΎπ§ππ¦β ) with different groove depths are approximately
equal to zero for the cavitating flow or the non-cavitating flow.
15 20 25 30 35 40 45 50
-0.2
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
1.8
2.0
2.2
2.4
2.6
K*yz
K*
zz
cavitating turbulent flow
non-cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
15 20 25 30 35 40 45 50-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
K*zx
K*xx
cavitating turbulent flow
non-cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(a) (b)
![Page 27: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/27.jpg)
15 20 25 30 35 40 45 50-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
K*zy
K*yy
cavitating turbulent flow
non-cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
15 20 25 30 35 40 45 50-4
-3
-2
-1
0
1
2
3
4
K*xz
K*yx
K*xy
Groove Depth hg (m)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
cavitating turbulent flow
non-cavitating turbulent flow
(c) (d)
Fig. 10 Influence of cavitation effect on stiffness coefficients for different
groove depth
(π = 18000πππ, π½ = 20β, βπ = 15ππ)
Fig. 11 shows the influence of cavitation effect on the stiffness coefficient of the
SGTB with different rotating speed. It can be seen that the direct stiffness coefficients (πΎπ§π§β , πΎππ₯ππ₯β , πΎππ¦ππ¦β ) with cavitating flow are larger than those with non-cavitating
flow at a specific rotating speed. This is also due to the interfacial hydrodynamic effect
caused by the bubbles. The cross-coupled stiffness coefficients (πΎππ₯ππ¦β , πΎππ¦ππ₯β ) of
cavitating flow are approximately equal to those of non-cavitating flow for different
rotating speeds. The difference of direct stiffness coefficients (πΎπ§π§β , πΎππ₯ππ₯β , πΎππ¦ππ¦β )between the cavitating flow and the non-cavitating flow gradually
increase as the rotating speed increases. This also can be explained by the fact that the
hydrodynamic effect enhances as the rotational speed increases, and the velocity
difference between the fluid and the bubble increases with this hydrodynamic effect. It
means that the influence of the cavitation effect on stiffness is increased with the
increase of rotating speed. The cross-coupled stiffness (πΎπ§ππ₯β , πΎππ₯π§β , πΎπ§ππ¦β , πΎππ¦π§β ) for
different rotating speed are also approximately equal to zero for the cavitating flow or
![Page 28: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/28.jpg)
the non-cavitating flow.
12000 15000 18000 21000 24000 27000 30000-0.2
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
K*yz
K*
zz
cavitating turbulent flow
non-cavitating turbulent flow
rpm
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
12000 15000 18000 21000 24000 27000 30000
-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
K*zy
K*yy
cavitating turbulent flow
non-cavitating turbulent flow
rpm
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(a) (b)
12000 15000 18000 21000 24000 27000 30000-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
K*zx
K*xx
cavitating turbulent flow
non-cavitating turbulent flow
rpm
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
12000 15000 18000 21000 24000 27000 30000-1.5
-1.2
-0.9
-0.6
-0.3
0.0
0.3
0.6
0.9
1.2
1.5
K*xz
K*yx
K*xy
cavitating turbulent flow
non-cavitating turbulent flow
rpm
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(c) (d)
Fig. 11 Influence of cavitation effect on stiffness coefficients for different rotating
speed
(π½ = 20β, βπ = 15ππ, βπ = 40ππ)
2οΌDamping coefficient
Fig. 12 shows the influence of cavitation effect on the water film damping
coefficients for different spiral angles. It can be seen that the cavitation effect on
damping coefficients is very weak. Because the damping coefficient is determined by
the perturbed pressure, and the perturbed pressure with cavitating flow is only affected
by the fraction πΆπ of bubble volume, and the variation of bubble volume fraction is
![Page 29: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/29.jpg)
very small for the cavitating flow. In addition, the direct damping coefficients (πΆπ§π§β , πΆππ₯ππ₯β , πΆππ¦ππ¦β ) decreases with the increasing of the spiral angle, and the cross-
coupled damping coefficients (πΆπ§ππ₯β , πΆππ₯π§β , πΆπ§ππ¦β , πΆππ¦π§β , πΆππ₯ππ¦β , πΆππ¦ππ₯β ) with the
cavitating flow and the non-cavitating flow are also approximately equal to zero for
different spiral angles.
5 10 15 20 25 30 35 40 45 50 55 60 65-0.2
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
C*yzC
*xz
C*
zz
cavitating turbulent flow
non-cavitating turnulent flow
Spiral angle (Β°)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
5 10 15 20 25 30 35 40 45 50 55 60 65-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
0.45
C*yx C
*zx
C*xx
cavitating turbulent flow
non-cavitating turnulent flow
Spiral angle (Β°)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
(a) (b)
5 10 15 20 25 30 35 40 45 50 55 60 65-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
C*xy
C*zy
C*yy
Spiral angle (Β°)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents cavitating turbulent flow
non-cavitating turnulent flow
(c)Fig. 12 Influence of cavitation effect on the damping coefficients for different
spiral angles
(π = 18000πππ, βπ = 15ππ, βπ = 40ππ)
Fig. 13 shows the influence of cavitation effect on the water film damping
coefficients for different groove depth. It can be seen that the cavitation effect on
damping coefficients is very weak for different groove depth. And the cross-coupled
damping coefficients (πΆπ§ππ₯β , πΆππ₯π§β , πΆπ§ππ¦β , πΆππ¦π§β , πΆππ₯ππ¦β , πΆππ¦ππ₯β ) with the cavitating flow
![Page 30: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/30.jpg)
and the non-cavitating flow are also approximately equal to zero for different groove
depth.
20 25 30 35 40 45
-0.4
0.0
0.4
0.8
1.2
1.6
2.0
2.4
2.8
C*
zz
Groove Depth hg (m)
Dim
en
sion
less
da
mp
ing
coeff
icie
nts
cavitating turbulent flow
non-cavitating turnulent flow
C*yzC
*xz
20 25 30 35 40 45-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
C*xx
cavitating turbulent flow
non-cavitating turnulent flow
Groove Depth hg (m)
Dim
ensi
on
less
da
mp
ing
co
effi
cien
ts
C*zxC
*yx
(a) (b)
20 25 30 35 40 45-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
C*yy
Groove Depth hg (m)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
cavitating turbulent flow
non-cavitating turnulent flow
C*zyC
*xy
(c)Fig. 13 Influence of cavitation effect on damping coefficients for different
groove depth.
(π = 18000πππ, π½=20β, βπ = 15ππ)
5.3 Influence of turbulence effect on dynamic characteristics
1οΌStiffness coefficient
Fig.14 shows the influence of turbulence effect on the water film stiffness
coefficients for different spiral angles. It can be seen that the direct stiffness
coefficients (πΎπ§π§β , πΎππ₯ππ₯β , πΎππ¦ππ¦β ) and the cross-coupled stiffness coefficients (πΎππ₯ππ¦β , πΎππ¦ππ₯β ) of turbulent cavitating flow are larger than those of the laminar
cavitating flow at a specified spiral angle. This may be explained by the fact that
additional Reynolds shear stress generates due to pulsating velocity in the turbulent
![Page 31: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/31.jpg)
flow, and equivalent viscosity of the turbulent water film increases, and the force of the
turbulent water film needed to produce unit displacement increases, so the stiffness of
the water film increases. In addition, it can be seen that the deviation of stiffness
coefficient between turbulence model and laminar model gradually decreases with the
increase of the spiral angle. This is due to the fact that the outward pumping
hydrodynamic effect increases first and then decreases with the increase of the spiral
angle. Therefore, when the spiral angle exceeds a certain value, the outward pumping
hydrodynamic effect becomes weak and the additional Reynolds shear stress of the
turbulent flow reduces. This indicates that the influence of turbulence effect on the
water film stiffness coefficients is weakened when the spiral angle exceeds a certain
value.
5 10 15 20 25 30 35 40 45 50 55 60 65
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
K*yz
K*
zz
Spiral angle (Β°)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents cavitating turbulent flow
cavitating laminar flow
5 10 15 20 25 30 35 40 45 50 55 60 65
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
K*zx
K*xx
Spiral angle (Β°)
cavitating laminar flow
cavitaing turbulent flow
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(a) (b)
5 10 15 20 25 30 35 40 45 50 55 60 65
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
K*yy
K*zy
Spiral angle (Β°)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
cavitating laminar flow
cavitating turbulent flow
5 10 15 20 25 30 35 40 45 50 55 60 65-1.5
-1.0
-0.5
0.0
0.5
1.0
1.5
K*xz
K*yx
K*xy
cavitating laminar flow
cavitating turbulent flow
Spiral angle (Β°)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(c) (d)
Fig. 14 Influence of turbulence effect on the water film stiffness coefficient for
different spiral angles
![Page 32: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/32.jpg)
(π = 18000πππ, βπ = 15ππ, βπ = 40ππ)
Fig. 15 shows the influence of turbulence effect on the water film stiffness coefficients
for different groove depth. It can also be seen that the direct stiffness coefficients (πΎπ§π§β , πΎππ₯ππ₯β , πΎππ¦ππ¦β ) and the cross-coupled stiffness coefficients of the turbulent
cavitating flow are larger than those of the laminar cavitating flow at a specified groove
depth. And the deviation of stiffness coefficients between turbulence model and laminar
model also reduce with the increase of groove depth. Similarly, the step hydrodynamic
effect becomes weak when the groove depth exceeds a given groove depth. This also
show that the influence of turbulence effect on water film stiffness coefficients are also
weakened when the groove depth exceeds a certain value.
15 20 25 30 35 40 45 50
-0.4
0.0
0.4
0.8
1.2
1.6
2.0
2.4
2.8
K*yz
K*
zz
cavitating laminar flow
cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
15 20 25 30 35 40 45 50-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
K*zx
K*xx
cavitating laminar flow
cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(a) (b)
15 20 25 30 35 40 45 50-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
K*zy
K*yy
cavitating laminar flow
cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
15 20 25 30 35 40 45 50-4
-3
-2
-1
0
1
2
3
4
K*xz
K*xy
K*yx
cavitating laminar flow
cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(c) (d)
Fig. 15 Influence of turbulence effect on the stiffness coefficient for different
groove depths
(π = 18000πππ, π½=20β, βπ = 15ππ)
Fig. 16 shows the influence of turbulence effect on the water film stiffness
![Page 33: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/33.jpg)
coefficients for different rotating speed. It can also be seen that the direct stiffness
coefficients and the cross-coupled stiffness coefficients of turbulent cavitating flow are
larger than those of the laminar cavitating flow at a specified rotating speed. However,
the deviation between turbulence model and laminar model is gradually increased with
the increasing of rotating speed, the reason for this result is that the additional Reynolds
shear stress of the turbulent flow increases with the rotational speed, leading to an
enhancement of turbulence intensity and a larger deviation between turbulence model
and laminar model.
12000 15000 18000 21000 24000 27000 30000
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
K*yz
K*
zz
cavitating laminar flow
cavitating turbulent flow
rpm
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
12000 15000 18000 21000 24000 27000 30000-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
K*zx
K*xx
cavitating laminar flow
cavitating turbulent flow
rpm
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(a) (b)
12000 15000 18000 21000 24000 27000 30000-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
K*zy
K*yy
cavitating laminar flow
cavitating turbulent flow
rpm
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
12000 15000 18000 21000 24000 27000 30000-1.5
-1.2
-0.9
-0.6
-0.3
0.0
0.3
0.6
0.9
1.2
1.5
K*xz
K*yx
K*xy
cavitating laminar flow
cavitating turbulent flow
rpm
Dim
ensi
on
less
sti
ffn
ess
coef
fici
ents
(c) (d)
Fig. 16 Influence of turbulence effect on the stiffness coefficient for different
rotating speed
(π½ = 20β, βπ = 15ππ, βπ = 40ππ)
![Page 34: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/34.jpg)
2οΌDamping coefficient
Fig.17 shows the influence of turbulence effect on the water film damping
coefficients for different spiral angles. It can also be seen that the direct damping
coefficients with turbulent cavitating flow are larger than those with the laminar
cavitating flow at a specified spiral angle. The equivalent viscosity of water film is
increased, and the force of the turbulent water film needed to produce unit velocity is
increased, so the damping of the turbulent water film is increased. Similar to Fig. 14,
the deviation of damping coefficients between turbulence model and laminar model is
also gradually decreased with the increase of the spiral angle, correspondingly, the
influence of the turbulence effect on the disturbance pressures s ποΏ½ΜοΏ½ , ποΏ½ΜοΏ½π₯ , ποΏ½ΜοΏ½π¦ is also
weakened. This indicates that the influence of turbulence effect on damping coefficient
of water film is weakened when the spiral angle exceeds a certain value.
5 10 15 20 25 30 35 40 45 50 55 60 65-0.2
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
C*yzC
*xz
C*
zz
cavitating laminar flow
cavitating turbulent flow
Spiral angle (Β°)
Dim
en
sio
nle
ss d
am
pin
g c
oeff
icie
nts
5 10 15 20 25 30 35 40 45 50 55 60 65-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
C*zxC
*yx
C*xx
cavitating laminar flow
cavitating turbulent flow
Spiral angle (Β°)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
(a) (b)
5 10 15 20 25 30 35 40 45 50 55 60 65-0.05
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
C*xyC
*zy
C*yy
cavitating laminar flow
cavitating turbulent flow
Spiral angle (Β°)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
(c)
Fig. 17 Influence of turbulence effect on the damping coefficients for different
![Page 35: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/35.jpg)
spiral angles.
(π = 18000πππ, βπ = 15ππ, βπ = 40ππ)
Fig.18 shows the influence of turbulence effect on the water film damping coefficients
for different groove depth. It can also be seen that the direct damping coefficients with
turbulent cavitating flow are larger than those with the laminar cavitating flow at a
specified groove depth. However, the deviation of damping coefficient between
turbulence model and laminar model decreases slowly with the increase of the groove
depth. This indicates that the turbulence effect on damping coefficient is not sensitive
to the groove depth.
15 20 25 30 35 40 45 50-0.4
0.0
0.4
0.8
1.2
1.6
2.0
2.4
C*yzC
*xz
C*
zz
cavitating laminar flow
cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
15 20 25 30 35 40 45 50
-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
C*zxC
*yx
C*xx
cavitating laminar flow
cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
(a) (b)
15 20 25 30 35 40 45 50-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
C*zyC
*xy
C*yy
cavitating laminar flow
cavitating turbulent flow
Groove Depth hg (m)
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
(c)Fig. 18 Influence of turbulence effect on the damping coefficients for different
groove depths.
(π = 18000πππ, π½=20β, βπ = 15ππ)
Fig.19 shows the influence of turbulence effect on the water film damping
coefficients for different rotating speed. It can be seen that the direct damping
coefficients with turbulent cavitating flow are larger than those with the laminar
![Page 36: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/36.jpg)
cavitating flow at a specified rotating speed. And the deviation of water film damping
coefficients between turbulence model and laminar model is also gradually increased
with the increasing of rotating speed.
12000 15000 18000 21000 24000 27000 30000-0.2
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
1.8
2.0
C*
zz
cavitating laminar flow
cavitating turbulent flow
rpm
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
C*xz C
*yz
12000 15000 18000 21000 24000 27000 30000-0.06
0.00
0.06
0.12
0.18
0.24
0.30
0.36
0.42
0.48
C*xx
cavitating laminar flow
cavitating turbulent flow
rpm
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
C*yx C
*zx
(a) (b)
12000 15000 18000 21000 24000 27000 30000-0.06
0.00
0.06
0.12
0.18
0.24
0.30
0.36
0.42
0.48
C*yy
cavitating laminar flow
cavitating turbulent flow
rpm
Dim
ensi
on
less
da
mp
ing
coef
fici
ents
C*xyC
*zy
(c)Fig. 19 Influence of turbulence effect on the damping coefficients for different
rotating speed.
(π½ = 20β, βπ = 15ππ, βπ = 40ππ)
3οΌBubble density
Fig. 20 shows the influence of turbulence effect on the bubble number for different
radius. It can be seen that the bubbles mainly locate at the groove-ridge edges and the
outer rim of bearing. This may be explained by the fact that the bubbles are formed due
to the pressure drop near the groove-ridge edges. The result also show that the bubbles
are approximately evenly distributed in dam region of bearing. In addition, the number
![Page 37: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/37.jpg)
of bubbles with laminar flow is more than the number of bubbles with turbulent flow
at a specified circumferential position. Because the Reynolds shear stress due to the
turbulent flow causes the equivalent viscosity of the turbulent liquid to increase,
consequently, the turbulent flow effect on dynamic pressure is enhanced, which
contributes to the generation of the bubbles.
-5 0 5 10 15 20 25 30 355
10
15
20
25
30
35
Nu
mb
er o
f b
ub
ble
s in
th
e co
ntr
ol
bo
dy
, n
Cavitation turbulent flow model
Cavitation laminar flow model
Circumferential direction (Β°)
(a)
-5 0 5 10 15 20 25 30 35
20
21
22
23
24
25
26
27
28
29
30
Cavitation turbulent flow model
Cavitation laminar flow model
Nu
mb
er o
f b
ub
ble
s in
th
e co
ntr
ol
bo
dy
n
Circumferential direction (Β°)
(b)
Fig. 20 Influence of turbulence effect on bubbles number for different
radius.(a) π=0.2 (b) π=0.5(π = 18000πππ, π½=20β,βπ = 15ππ, βπ = 40ππ)
![Page 38: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/38.jpg)
6. Conclusions
In this paper, a theoretical model of water lubrication based on two-phase flow
considering both turbulence and cavitation bubble effect is established. In the turbulent
state, the evolution of the size distribution of the bubbles is described by the bubble
transport equation considering the mechanism of breakage and coalescence. The
interval discretization method is used to solve the bubble transport equation. The size
and spatial equilibrium distribution of the bubbles under turbulent flow are predicted
by solving the generalized Reynolds equation and the bubble transport equation
simultaneously. The dynamic characteristics of water-lubricated SGTB under turbulent
cavitating flow is analyzed. Based on discussion, the following conclusion can be
drawn:
(1) The bubbles are mainly concentrated near the spiral-groove edges of the SGTB.
(2) The small-sized bubbles is much more than the large-sized bubbles for the water-
lubricated SGTB under the cavitation flow condition.
(3) The cavitation has a significant effect on the direct stiffness coefficients, and the
cavitation hardly affect the cross-coupled stiffness coefficients and the damping
coefficients of SGTB.
(4) Compared to the cavitation, the turbulence has a much greater effect on the
dynamic characteristics of the high-speed water-lubricated SGTB.
7. DeclarationFunding
The authors gratefully acknowledge the support provided by the National Science
Foundation of China through grant Nos. 51635004, 11472078.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article.
Authorsβ contributionsThe authorβ contributions are as follows: Xiaohui Lin was in charge of the whole trial;
Xiaohui Lin wrote the manuscript; Shun Wang is responsible for model calculation and
drawing;Shuyun Jiang revises and reviews the manuscript; Shaowen Zhang assisted
with laboratory analyses.
![Page 39: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/39.jpg)
Competing interests
The authors declare no competing financial interests.
Consent for publication
Not applicable
Ethics approval and consent to participate
Not applicable
AcknowledgementsNot applicable
Nomenclature π=temperature viscosity coefficient (1πΎ)πΆππ(π = π§, ππ₯, ππ¦; π = π§, ππ₯, ππ¦)nine damping coefficients ((
πβ π π , πβ π πππ , π β π , πβ πβ π πππ ))πΆππβ (π = π§, ππ₯, ππ¦; π = π§, ππ₯, ππ¦) nine dimensionless damping coefficients
(dimensionless)πΆπ€=the water volume fraction (dimensionless)πΆπ·π =drag coefficient (dimensionless)
= the bubble volume fraction (dimensionless)ππ£=the specific heat capacity of water ( π½ππβ πΎ) ππ=the specific heat capacity of thrust disk ( π½ππβ πΎ) π1, π2=experimental constant (dimensionless)ππ=the function of π§β containing the parameter βπβ (dimensionless)πππ(π, π) =the equilibrium distribution function of bubble diameter at (π, π) coordinates system ( 1π6)π(π, π) =the distribution function of bubble diameter at (π, π) coordinates system
( 1π6)
gC
![Page 40: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/40.jpg)
π= gravity acceleration (ππ 2)ππ= the function of π§β containing the parameter βπβ (dimensionless)β=the film thickness (π)βπ=the groove depth (π)β0=Static film thickness (π)βπ =the film thickness of the ridge (π)βπ= initial film thickness of bubble (π)βπ= critical rupture film thickness of bubble (π)βπβ= the shear stress parameter βπβ of the Couette flow (dimensionless)ππ =heat transfer coefficient of the stationary ring ( ππβ πΎ)ππ=heat transfer coefficient of the thrust disk ( ππβ πΎ)ππ€= heat transfer coefficient of water ( ππβ πΎ)ππ=the turbulence correction coefficient of radial direction (dimensionless)ππ= the turbulence correction coefficient of circumferential direction (dimensionless)πΎππ(π = π§, ππ₯, ππ¦; π = π§, ππ₯, ππ¦) nine stiffness coefficients (ππ , ππππ , π, πβ ππππ)πΎππβ (π = π§, ππ₯, ππ¦; π = π§, ππ₯, ππ¦)nine dimensionless stiffness coefficients
(dimensionless)π= the circumferential discrete node(dimensionless)π(π)=the interface item at (π, π) coordinates system ( ππ3)π =the number of spiral groove (dimensionless)ππ = the number density of nuclei ( 1π3)ππππ (π)= the number of bubbles in the kth volume range(ππ, ππ+1) at equilibrium
![Page 41: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/41.jpg)
(dimensionless)ππ(π)= the number of bubbles in the kth volume range(ππ, ππ+1) (dimensionless)π=the water film pressure ( ππ2)π0=the static pressure ( ππ2)ππ§ , πππ₯ , πππ¦ , ποΏ½ΜοΏ½ , ποΏ½ΜοΏ½π₯ , ποΏ½ΜοΏ½π¦=perturbed pressure components (ππ3 , ππ2 , πβ π π3 , πβ π π2)ππ£=the vaporization pressure of water ( ππ2)πππ=the pressure at inter radius ( ππ2)πππ, πππ=the volume flow (π3π )π = the radial coordinate in polar coordinates (π)ππ= the base circle radius of spiral line (π)
πππ= thrust disk inner radius (π)πππ’π‘= thrust disk outer radius (π)π π=the statistical average radius of bubbles (π)π π=Reynolds number (dimensionless)ππ=the source terms ( 1π6π )π‘ =the time (π )π = Water film temperature (πΎ)ππ=the average temperature of water film (πΎ)ππ =the temperature of stationary ring (πΎ)ππ=the temperature of thrust plate (πΎ)π0= the ambient temperature (πΎ)π= the max circumferential velocity of thrust disk (ππ )οΏ½Μ οΏ½= the circumferential average velocity of water (ππ )
![Page 42: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/42.jpg)
π’ = the circumferential velocity of water at (π, π) coordinates system (ππ )π’π= the bubble velocity (ππ )
π₯π=the diameter of bubbles (π)π§ = the coordinate of the film thickness (π)π§β= the coordinate of the dimensionless film thickness (π§β = π§β)
π§π = the coordinate of thickness of the stationary ring (π)
π§π= the coordinate of thickness of the thrust disk (π)
Greek LettersπΌ =experimental constant (dimensionless)πΌπ= air convection heat transfer coefficient ( ππ2β πΎ) πΌπ€=forced convective heat transfer coefficient of water ( ππ2β πΎ)π½= the helix angle of bearing (deg)π½π= the parameter in coalescence kernel (dimensionless)πΏ=thickness coordinate of stationary ring (π)πΏπ =the thickness of the stationary ring (π)πΏπ=the thickness of thrust plate (π)π=the spiral coordinates (π)ππ= the base circle radius of spiral line at (π, π) coordinates system (π) πππ= radius of the oil inlet at (π, π) coordinates system (π)πππ’π‘= radius of the oil outlet at (π, π) coordinates system (π)π = turbulent kinetic energy dissipation rate per unit mass (π2π 3 )π =the spiral coordinates (deg)
![Page 43: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/43.jpg)
ππππ₯= the maximum value of spiral coordinates (deg)π 1, π 2= the adjustable parameters (dimensionless)ππ=the groove-to-land ratio (dimensionless)ππ=the groove-to-dam ratio (dimensionless)ππ€β = dimensional dynamic viscosity of water (dimensionless)
οΏ½Μ οΏ½π€= the average dynamic viscosity of water (πβ π π2)ππ€= the dynamic viscosity of water (πβ π π2)ππ€0= dynamic viscosity of water at π0 (πβ π π2)ππ= the density of bubble (πππ3)ππ€=the density of water (πππ3)ππ=the density of thrust disk (πππ3)π=the surface tension of bubble (ππ)π = the diameter of bubbles (π)ππ = the diameter of gas nucleus (π) ππmax =the maximum diameter of gas nucleus (π)π = the angular velocity of bearing rotation (ππππ )
π = the diameter of bubbles (π)ππππ₯= the maximum diameter of bubbles (π)π₯π‘ = The time step (π )π₯π§, π₯ππ₯, π₯ππ¦= perturbations quantity (π)
References
[1] B C Majumdar, R Pai, D J Hargreaves. Analysis of water-lubricated journal
![Page 44: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/44.jpg)
bearings with multiple axial grooves. Proc Inst Mech Eng Part J: J Eng Tribol,
2004, 218: 135β146.
[2] Y M Zhao, C Wei, S Yuan, et al. Theoretical and Experimental Study of Cavitation
Effects on the Dynamic Characteristic of Spiral-Groove Rotary Seals(SGRSs).
Tribology letters, 2016, 64(3): 1-18.
[3] H G Elrod. A Cavitation Algorithm. ASME J. Lubr. Technol,1981, 103(3): 350β
354.
[4] Y Song, C Gu. Development and validation of a three-dimensional computational
fluid dynamics analysis for journal bearings considering cavitation and conjugate
heat transfer. Journal of Engineering for Gas Turbines and Power, 2015,
137(12):122502.
[5] Y Wang, Z Yin, D Jiang, G Gao, et al. Study of the lubrication performance of water-
lubricated journal bearings with CFD and FSI method. Industrial Lubrication and
Tribology, 2016, 68(3): 341-348.
[6] E P Grando, M Priest, A T Prata. A two-phase flow approach to cavitation modelling
in journal bearings. Tribology Letters, 2006, 21(3): 233-244.
[7] H Liu, H Xu, P J Ellison, et al. Application of computational fluid dynamics and
fluidβ structure interaction method to the lubrication study of a rotorβbearing
system. Tribology Letters, 2010, 38(3): 325-336.
[8] J L Nikolajsen. The Effect of Aerated Oil on the Load Capacity of a Plain Journal
Bearing. Tribology Transactions, 1999, 42(1): 58β62.[9] S Choi, K W Kim. Analysis of Bubbly Lubrication in Journal Bearings. JSME
International Journal Series C, 2002, 45(3): 802β808.
[10] Q Li, S L Liu, X H Pan, et al. A new method for studying the 3D transient flow of
misaligned journal bearings in flexible rotor-bearing systems. Journal of Zhejiang
University Science A, 2012, 13(4): 293-310.
[11] Q Li, G Yu, S Liu, S Zheng. Application of computational fluid dynamics and fluid
structure interaction techniques for calculating the 3D transient flow of journal
bearings coupled with rotor systems. Chinese Journal of Mechanical Engineering,
2012, 25(5): 926-932.
![Page 45: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/45.jpg)
[12] J Hu, W Wei, M Wu, S Yuan. Numerical investigation of the airβoil two-phase
flow inside an oil-jet lubricated ball bearing. International Journal of Heat and
Mass Transfer, 2014, 68: 85-93.
[13] X Zhang, Z Yin, G Gao, Z Li. Determination of stiffness coefficients of
hydrodynamic water-lubricated plain journal bearings. Tribology International.,
2015, 85: 37-47.
[14] Q Li, S Zhang, L Ma, W Xu, S Zheng. Stiffness and damping coefficients for
journal bearing using the 3D transient flow calculation. Journal of Mechanical
Science and Technology, 2017, 31(5): 2083-2091.
[15] Y Chen, Y Sun, Q He, J Feng. Elastohydrodynamic Behavior Analysis of Journal
Bearing Using FluidβStructure Interaction Considering Cavitation,β Arabian
Journal for Science and Engineering, 2019, 44: 1305β1320.
[16] D Sun, S Y Li, C W Fei, et al. Investigation of the effect of cavitation and journal
whirl on static and dynamic characteristics of journal bearing. Journal of
Mechanical Science and Technology, 2019, 33(1): 77-86.
[17] H Liu, H Xu, Y Zhang. The influence of sea water in oil emulsion on bearing
performance. Part J: Journal of Engineering Tribology, 2009, 223(3): 457-468.
[18] Q Y Lin, Z Y Wei, N Wang, W Chen. Analysis on the lubrication performances of
journal bearing system using computational fluid dynamics and fluidβstructure
interaction considering thermal influence and cavitation. Tribology International,
2013, 64: 8β15.
[19] S B Shenoy, R Pai, D Rao, R B Pai. Elasto-hydrodynamic lubrication analysis of
full 360 journal bearing using CFD and FSI techniques. World Journal of
Modelling and Simulation, 2009, 5(4).
[20] X H Lin, S Y Jiang, C B Zhang. Thermohydrodynamic analysis of high speed
water-lubricated spiral groove thrust bearing considering effects of cavitation,
inertia and turbulence. Tribology International, 2018, 119: 645-658.
[21] X H Lin, S Y Jiang, C B Zhang, X Liu. Thermohydrodynamic Analysis of High-
Speed Water-Lubricated Spiral Groove Thrust Bearing Using Cavitating Flow
Model. Journal of Tribology, 2018, 140: 0517031-05170312.
[22] X H Lin, R Q Wang, S W Zhang, C B Zhang, S Y Jiang. Study of Cavitation
Bubbles Evolution for High-Speed Water-Lubricated Spiral Groove Thrust
Bearings. ASME Journal of Tribology, 2019, 141: 051703-1.
[23] X H Lin, R Q Wang, S W Zhang, S Y Jiang. Study on dynamic characteristics for
![Page 46: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/46.jpg)
high speed water-lubricated spiral groove thrust bearing considering cavitating
effect. Tribology International, 2019, 143:106022.
[24] C W Ng, C H T Pan. A Linearized Turbulent Lubrication Theory. ASME J. Fluids
Eng, 1965, 87(3): 675β688.
[25] J O Hinze. Turbulence. McGraw-Hili Book Company., Inc., New York, 1959.
[26] J R Ni, G Q Wang, H W Zhang. The Basic Theory of Solid Liquid Two-Phase
Flow and Its Latest Applications. Science Press, 1991.
[27] L J Guo. Two-Phase and Multiphase Flow Dynamics. Xiβan Jiaotong University
Press, 2002.
[28] J Solsvik, H A Jakobsen. Evaluation of Weighted Residual Methods for the
Solution of a Population Balance Model Describing Bubbly Flows: The Least
Squares, Galerkin, Tau, and Orthogonal Collocation Methods. Industrial &
Engineering Chemistry Research, 2013, 52(45): 15988-16013.
[29] M L Billet. Cavitation Nuclei Measurements. In: New Orleans, Louisiana, USA.
Proc. of 2nd Intcrn Symp on Cavitation Inception, 1984, 33-42.
[30] S Kumar, D Ramkrishna. On The Solution of Population Balance Equations By
Discretization--I. A Fixed Pivot Technique. Chemical Engineering Science, 1996,
51(8): 1311-1332.
[31] C A Dorao, H A Jakobsen. A Least Squares Method for The Solution of Population
Balance Problems. Computers and Chemical Engineering, 2006, 30: 535β547.
![Page 47: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/47.jpg)
Appendix A. The derivation of perturbation generalized Reynolds equations
The generalized Reynolds equation in (π, π) coordinate system,
πππ ( β3πππππ (ππΆπ€πππ )) + πππ ( β3πππππ (ππΆπ€πππ ))= πππ (βππΆπ€2 + β3πππππ(π)) + πππ (β3ππ€οΏ½Μ οΏ½2πΆπ€ππππ ) + πΆπ€π πβππ‘ (A1)
Ignoring the infinitesimal value of higher order, the expression of β3 is as follows,β3 β β03 + 3β02(π₯π§ + πcosππ₯ππ¦ β πsinππ₯ππ₯) (A2)
Bring Eq. (A2) and Eq. (57) into the right side of Eq. (A1),
πππ (βππΆπ€2 + β3οΏ½Μ οΏ½π€πππ(π)) + πππ (β3ππ€οΏ½Μ οΏ½2πΆπ€οΏ½Μ οΏ½π€ππ ) + πΆπ€π πβππ‘= πππ (πΆπ€πβ0) + πππ ( β03οΏ½Μ οΏ½π€ππ) + πππ (β03ππ€π’2πΆπ€οΏ½Μ οΏ½π€ππ β03)+ [ πππ (πΆπ€π)+ πππ (π(π)οΏ½Μ οΏ½π€ππ 3β02) + πππ (β3ππ€οΏ½Μ οΏ½2πΆπ€οΏ½Μ οΏ½π€ππ 3β02)] π₯π§+ [ πππ (πΆπ€ππcosπ) + πππ (π(π)οΏ½Μ οΏ½π€ππ 3β02πcosπ) + πππ (β3ππ€οΏ½Μ οΏ½2πΆπ€οΏ½Μ οΏ½π€ππ 3β02πcosπ)] π₯ππ¦β[ πππ (πΆπ€ππsinπ) + πππ (π(π)οΏ½Μ οΏ½π€ππ 3β02πsinπ) + πππ (β3ππ€π’2πΆπ€οΏ½Μ οΏ½π€ππ 3β02πsinπ)] π₯ππ₯+πΆπ€ππ₯οΏ½ΜοΏ½ + πΆπ€π2cosππ₯οΏ½ΜοΏ½π¦ β πΆπ€π2sinππ₯οΏ½ΜοΏ½π₯
(A3)
Convert the coordinate system of Eq. (A3) to the (π, π) coordinate system, and
compare similar terms,ππ§ = πππ (πΆπ€π) + πππ (π(π)ππ€ππ 3β02) + πππ (ππ€π’2πΆπ€ππ€ππ 3β02)β πππ ( ππ€οΏ½Μ οΏ½2πΆπ€πππ€ππcotπ½ 3β02) β πππ ( 3β02ππ€πππ ππΆπ€π0ππ )β [ πππ (3β02πππ€ππ ππΆπ€π0ππ ) β πππ ( 3β02cotπ½ππ€ππ ππΆπ€π0ππ )β πππ ( 3β02cotπ½ππ€ππ ππΆπ€π0ππ ) + πππ ( 3β02πcot2π½ππ€ππ ππΆπ€π0ππ )]
(A4)
πππ₯ = πππ (πΆπ€ππΊ1(π, π)) + πππ (π(π)ππ€ππ 3β02πΊ1(π, π)) + πππ (ππ€π’2πΆπ€ππ€ππ 3β02πΊ1(π, π))β πππ (ππ€π’2πΆπ€πππ€ππ 3β02cotπ½πΊ1(π, π)) + πππ ( 3β02ππ€πππ ππΆπ€π0ππ πΊ1(π, π))+ [ πππ (3β02πππ€ππ ππΆπ€π0ππ πΊ1(π, π)) β πππ ( 3β02cotπ½ππ€ππ ππΆπ€π0ππ πΊ1(π, π))β πππ ( 3β02cotπ½ππ€ππ ππΆπ€π0ππ πΊ1(π, π)) + πππ ( 3β02πcot2π½ππ€ππ ππΆπ€π0ππ πΊ1(π, π))]
(A5)
![Page 48: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/48.jpg)
πππ¦ = πππ (πΆπ€ππΊ2(π, π)) + πππ (π(π)ππ€ππ 3β02πΊ2(π, π)) + πππ (ππ€π’2πΆπ€ππ€ππ 3β02πΊ2(π, π))β πππ (ππ€π’2πΆπ€πππππ 3β02cotπ½πΊ2(π, π)) β πππ ( 3β02ππ€πππ ππΆπ€π0ππ πΊ2(π, π))β [ πππ (3β02πππ€ππ ππΆπ€π0ππ πΊ2(π, π)) β πππ ( 3β02cotπ½ππ€ππ ππΆπ€π0ππ πΊ2(π, π))β πππ ( 3β02cotπ½ππ€ππ ππΆπ€π0ππ πΊ2(π, π)) + πππ ( 3β02πcot2π½ππ€ππ ππΆπ€π0ππ πΊ2(π, π))]
(A6)ποΏ½ΜοΏ½ = πΆπ€π (A7)
ποΏ½ΜοΏ½π₯ = βπΆπ€ππΊ1(π, π) (A8)
ποΏ½ΜοΏ½π¦ = πΆπ€ππΊ2(π, π) (A9)
where πΊ1(π, π) = πsin(π + π(π))πΊ2(π, π) = πcos(π + π(π)) (A10)
![Page 49: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/49.jpg)
Figures
Figure 1
Structure schematic of inward-pumping SGTB
Figure 2
Coordinates transformation of spiral patterns
![Page 50: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/50.jpg)
Figure 3
Finite volume separated by grooveβridge boundary
![Page 51: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/51.jpg)
Figure 4
Motion of a SGTB
![Page 52: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/52.jpg)
Figure 5
Flow diagram of overall algorithm
![Page 53: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/53.jpg)
Figure 6
Photographs of bubble distributions under both rotational speeds: (a) 15000 rpm, (b) 18000 rpm
![Page 54: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/54.jpg)
Figure 7
Comparison of the theoretical probability density distributions with the experimental ones: (a)15000rpm(b)18000rpm
![Page 55: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/55.jpg)
Figure 8
Comparison of the theoretical stiffness coeοΏ½cients with the experimental ones at 15000rpm
![Page 56: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/56.jpg)
Figure 9
InοΏ½uence of cavitation effect on the stiffness coeοΏ½cient for different spiral angles
![Page 57: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/57.jpg)
Figure 10
InοΏ½uence of cavitation effect on stiffness coeοΏ½cients for different groove depth
![Page 58: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/58.jpg)
Figure 11
InοΏ½uence of cavitation effect on stiffness coeοΏ½cients for different rotating speed
![Page 59: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/59.jpg)
Figure 12
InοΏ½uence of cavitation effect on the damping coeοΏ½cients for different spiral angles
![Page 60: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/60.jpg)
Figure 13
InοΏ½uence of cavitation effect on damping coeοΏ½cients for different groove depth.
![Page 61: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/61.jpg)
Figure 14
InοΏ½uence of turbulence effect on the water οΏ½lm stiffness coeοΏ½cient for different spiral angle
![Page 62: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/62.jpg)
Figure 15
InοΏ½uence of turbulence effect on the stiffness coeοΏ½cient for different groove depths
![Page 63: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/63.jpg)
Figure 16
nοΏ½uence of turbulence effect on the stiffness coeοΏ½cient for different otating speed
![Page 64: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/64.jpg)
Figure 17
InοΏ½uence of turbulence effect on the damping coeοΏ½cients for different spiral angles.
![Page 65: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/65.jpg)
Figure 18
InοΏ½uence of turbulence effect on the damping coeοΏ½cients for different groove depths.
![Page 66: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/66.jpg)
Figure 19
nοΏ½uence of turbulence effect on the damping coeοΏ½cients for different rotating speed.
![Page 67: Dynamic Characteristics of High-speed Water- lubricated](https://reader034.vdocument.in/reader034/viewer/2022042421/6261144db221fb42ed5ac31e/html5/thumbnails/67.jpg)
Figure 20
InοΏ½uence of turbulence effect on bubbles number for different radius.(a) =0.2 (b) =0.5 ( = 18000,=20, = 15, = 40)