Improvement of Thrust Bearing Calculation Considering the Convectional Heating within the Space between the Pads

A modern thrust bearing tool is used to estimate the behavior of tilting pad thrust bearings not only in the oil film between pad and rotating collar, but also in the space between the pads. The oil flow in the space significantly influences the oil film inlet temperature and the heating of pad and collar. For that reason, it is necessary to define an oil mixing model for the space between the pads. In the bearing tool, the solutions of the Reynolds equation including a cavitation model, the energy equation and the heat transfer equation are done iteratively with the finite volume method by considering a constant flow rate. Both effects—laminar/turbulent flow and centrifugal force—are considered. The calculation results are compared with measurements done for a flooded thrust bearing with nominal eight tilting pads with an outer diameter of 180 mm. The heat convection coefficients for the pad surfaces mainly influence the pad temperature field and are adjusted to the measurement results. In the following paper, the calculation results for variable space distances, influence of different parameters on the bearing behavior and operating condition at high load are presented.


Introduction
The thrust bearing tool used in GE Renewable Energy Hydro for vertically-aligned hydro machines considered the distance between the pads in a range near 20% of the pad length.GE Renewable Energy Hydro started an investigation to study the distance influence for a higher range.Measurements on a test rig at Gdansk University of Technology were done for a flooded tilting pad thrust bearing.For the recalculation of the measurement, the tool Combros-A, developed at Clausthal University of Technology [1], was used.Finally, the influences of selected parameters, such as the distance between the pads, support system, speed, load, oil flow and the operating condition at high load were investigated.Loads higher than load during stationary operating conditions could occur for transient conditions such as a shut-down of a machine.The quasi-stationary solution is used.
The oil flow in the space has a major influence on the oil film inlet temperature T A and the heating of pad and collar.To estimate the oil film inlet temperature T A [2,3], the hot oil carry-over coefficient is used, which is the ratio inlet temperature rise (T A -T 0 ) to collar temperature rise (T C -T 0 ).The temperature T 0 is the oil supply temperature: the bearing housing temperature in a flooded bearing.The coefficient depends on the ratio oil film inlet gap to oil film outlet gap, varies between 0.6 and 0.9 and is based on experiments for several thrust bearing types.The temperatures T A and T C are constant average values and only a small variation of the distance between the pads is investigated.Further investigations of [4,5] are based on experiments wherein the distance between the pads was changed by removing pads or decreasing the pad number, for example from eight to four.For the calculation, the oil film temperature is a constant average value.The difference between the oil film outlet and inlet temperature, T E and T A , is estimated by power loss and flow.For the estimation of the average collar temperature T C , the Peclet number, which is the ratio of the convective to conductive heat transport, is used.The bearing tool [6] uses the warm oil mixing factor in the range between 0 and 1 to describe the heat transport of the oil film outlet flow to oil film inlet of the following pad.The oil film temperature is a 3D field whereas the warm oil mixing factor depends on thrust bearing type and is an input value in the bearing tool.The influence of the oil supply on the oil mixing in the space between the pads is studied for a radial bearing in [7].The three basic equations, Reynolds, energy and heating equations, are simultaneously solved.Numerical analysis of thermo-elastic hydrodynamic thrust bearings by using computational fluid dynamics (CFD) was done by [8,9].In this paper, a modern thrust bearing tool is described which estimates the heating of the oil in the space between the pads using the energy and heating equation.

Thrust Bearing Calculation Tool with the Oil Mixing Model
The tool Combros-A is a TEHD (thermo-elastic hydrodynamic) bearing calculation software [1,10] in which the Reynolds, energy and heating equations are solved iteratively.The oil film gap and pressure are 2D fields (ϕ,r).The load capacity is primary defined by the minimum oil film gap.The heating of oil, pad and runner is estimated.Among wide range of physical effects considered in the calculation, a few key features should be mentioned: 3D fields (ϕ,r,z) for velocity, temperature, viscosity and density including the effect of the centrifugal force, approximation of thermal and mechanical pad and collar deformation, local laminar and turbulent flow transition, cavitation model in divergent gap zones based on [11].
The model of the oil mixing in the space between the tilting pads is shown in Figure 1.The basic theory of the mixing model was published by [12].For the replacement of the oil film outlet flow at trailing edge Q E by fresh oil with the temperature T 0 , the oil replacing factor f r is used.Q T is the oil flow with the fresh oil temperature T 0 and Q A is the oil film inlet flow at leading edge with the oil film inlet temperature T A .The coefficient depends on the ratio oil film inlet gap to oil film outlet gap, varies between 0.6 and 0.9 and is based on experiments for several thrust bearing types.The temperatures TA and TC are constant average values and only a small variation of the distance between the pads is investigated.Further investigations of [4,5] are based on experiments wherein the distance between the pads was changed by removing pads or decreasing the pad number, for example from eight to four.For the calculation, the oil film temperature is a constant average value.The difference between the oil film outlet and inlet temperature, TE and TA, is estimated by power loss and flow.For the estimation of the average collar temperature TC, the Peclet number, which is the ratio of the convective to conductive heat transport, is used.The bearing tool [6] uses the warm oil mixing factor in the range between 0 and 1 to describe the heat transport of the oil film outlet flow to oil film inlet of the following pad.The oil film temperature is a 3D field whereas the warm oil mixing factor depends on thrust bearing type and is an input value in the bearing tool.The influence of the oil supply on the oil mixing in the space between the pads is studied for a radial bearing in [7].The three basic equations, Reynolds, energy and heating equations, are simultaneously solved.Numerical analysis of thermo-elastic hydrodynamic thrust bearings by using computational fluid dynamics (CFD) was done by [8,9].In this paper, a modern thrust bearing tool is described which estimates the heating of the oil in the space between the pads using the energy and heating equation.

Thrust Bearing Calculation Tool with the Oil Mixing Model
The tool Combros-A is a TEHD (thermo-elastic hydrodynamic) bearing calculation software [1,10] in which the Reynolds, energy and heating equations are solved iteratively.The oil film gap and pressure are 2D fields (φ,r).The load capacity is primary defined by the minimum oil film gap.The heating of oil, pad and runner is estimated.Among wide range of physical effects considered in the calculation, a few key features should be mentioned: 3D fields (φ,r,z) for velocity, temperature, viscosity and density including the effect of the centrifugal force, approximation of thermal and mechanical pad and collar deformation, local laminar and turbulent flow transition, cavitation model in divergent gap zones based on [11].
The model of the oil mixing in the space between the tilting pads is shown in Figure 1.The basic theory of the mixing model was published by [12].For the replacement of the oil film outlet flow at trailing edge QE by fresh oil with the temperature T0, the oil replacing factor fr is used.QT is the oil flow with the fresh oil temperature T0 and QA is the oil film inlet flow at leading edge with the oil film inlet temperature TA. ( The oil temperatures at oil inlet and outlet TA and TE are 2D fields (r,z) like the flow fields.The fresh oil replacing factor fr depends on the thrust bearing type.Heating of the oil in the space and the resulting 3D oil temperature field (φ,r,z) are estimated.The flows at oil inlet and outlet Q A and Q E are 2D fields (r,z).The hydraulic integrated oil flow balance is calculated as The oil temperatures at oil inlet and outlet T A and T E are 2D fields (r,z) like the flow fields.The fresh oil replacing factor f r depends on the thrust bearing type.Heating of the oil in the space and the resulting 3D oil temperature field (ϕ,r,z) are estimated.

Bearing Test Rig for Measurement
Experimental investigations with different spaces between the pads (number of pads) and two support positions (removable ribs) at different speeds and loads were done on the bearing test rig [13].A cross section of the bearing test rig is shown in Figure 2. The inlet of the rated oil flow Q in is arranged on both bearing sides and there is one pipe for oil outlet flow Q out above the bearing housing.The configuration of pads and support systems are presented in Figure 3.The most important geometrical and operating data are summarized in Table 1.

Bearing Test Rig for Measurement
Experimental investigations with different spaces between the pads (number of pads) and two support positions (removable ribs) at different speeds and loads were done on the bearing test rig [13].A cross section of the bearing test rig is shown in Figure 2. The inlet of the rated oil flow Qin is arranged on both bearing sides and there is one pipe for oil outlet flow Qout above the bearing housing.The configuration of pads and support systems are presented in Figure 3.The most important geometrical and operating data are summarized in Table 1.Lubricants 2018, 6, x FOR PEER REVIEW 3 of 12

Bearing Test Rig for Measurement
Experimental investigations with different spaces between the pads (number of pads) and two support positions (removable ribs) at different speeds and loads were done on the bearing test rig [13].A cross section of the bearing test rig is shown in Figure 2. The inlet of the rated oil flow Qin is arranged on both bearing sides and there is one pipe for oil outlet flow Qout above the bearing housing.The configuration of pads and support systems are presented in Figure 3.The most important geometrical and operating data are summarized in Table 1.The pads were equipped with twelve temperature sensors.Nine of them were installed in the pad, seven (red dots) in the holes drilled through the whole pad thickness to monitor temperature next to hydrodynamic oil film gap on the pad surface and two others (blue dots) about 3 mm from the pad bottom.Another three sensors (green dots) were located in the collar near the inner and outer radius and in the middle of thrust bearing about 2 mm below the collar surface.The arrangement of temperature sensors for pad and collar is shown in Figure 4.The pads were equipped with twelve temperature sensors.Nine of them were installed in the pad, seven (red dots) in the holes drilled through the whole pad thickness to monitor temperature next to hydrodynamic oil film gap on the pad surface and two others (blue dots) about 3 mm from the pad bottom.Another three sensors (green dots) were located in the collar near the inner and outer radius and in the middle of thrust bearing about 2 mm below the collar surface.The arrangement of temperature sensors for pad and collar is shown in Figure 4. Based on increase of the oil temperature in the bearing housing ΔT = Tout − Tin and oil flow rate Qin, the total power loss Pbear is determined.Constant values for oil density ρ = 867 kg/m 3 and oil specific heat cp = 2112.5J/kg•K were assumed.The conductive power loss over the bearing housing to the environment is neglected.

Combros-A Calculation
The thrust bearing is fully flooded with oil, the angular pad pivot positions are 0.6 (offset) and 0.5 (central) and the pad number varies between four and nine.The variables in Table 2 are estimated with measurement results at speed 4000 rpm and a specific pressure of 3 MPa.The pad surface overlapping factor fsurf is the ratio pad area to ring area.As shown in the following Section 5, the fresh oil replacing factor corresponds to the bearing type and not to operating data such as speed and specific pressure.Based on increase of the oil temperature in the bearing housing ∆T = T out − T in and oil flow rate Q in , the total power loss P bear is determined.Constant values for oil density ρ = 867 kg/m 3 and oil specific heat c p = 2112.5J/kg•K were assumed.The conductive power loss over the bearing housing to the environment is neglected.

Combros-A Calculation
The thrust bearing is fully flooded with oil, the angular pad pivot positions are 0.6 (offset) and 0.5 (central) and the pad number varies between four and nine.The variables in Table 2 are estimated with measurement results at speed 4000 rpm and a specific pressure of 3 MPa.The pad surface overlapping factor f surf is the ratio pad area to ring area.As shown in the following Section 5, the fresh oil replacing factor corresponds to the bearing type and not to operating data such as speed and specific pressure.In [8], the effect of the heat convection coefficient at the pad surfaces of a tilting pad thrust bearing is studied by using CFD investigations.The tool Combros-A defines average heat transfer coefficients at the five surfaces at the leading edge, trailing edge, outer radius, inner radius and bottom.The pad heat convection coefficients from Table 3 are estimated by comparison of measurement at speed of 4000 rpm and specific pressure of 3 MPa.The value of the heat convection coefficient on the bottom side is high and the physical background was not finally clarified.Only with this high value do the calculated pad temperatures near pad bottom agree with the measurements.The surrounding surface temperatures depend on the temperatures in the bearing housing and therefore from the heating of the oil in the housing caused by bearing power loss.The change of the heat convection coefficients at different speeds seems to be low as shown in the following Section 6.The values from Table 3 are used in this paper for all operating conditions investigated.The pads are arranged on the both sides of the collar, Figure 2. The symmetrical collar temperature curve from the middle of the collar to both pads allows the collar to divide in two parts and no heat transfer to be estimated at the middle surface.The heat transfers near the outer and inner collar diameter are low because of their small surfaces in comparison to the middle collar surface.The values used in the calculation tool are shown in Table 4; the surrounding surface temperatures are adapted to the temperatures in the bearing housing.The heating of the oil in the film, the oil in the space between the pads, the pad and collar are shown in Figure 5 for 8 and 4 pads at rotational speed of 4000 rpm and specific pressure of 3 MPa.The graphical outputs show the temperature fields at the mean diameter D mean = 135 mm.More pads are creating higher dissipation in the oil film and therefore higher heating of space, pad and collar.Calculated power losses in oil film, the space and the heat flow to pad and collar are summarized in Table 5.The percentage value in the parentheses is the ratio heat to sum of both calculated power losses, P f,pad + P f,space .
Table 5.The percentage value in the parentheses is the ratio heat to sum of both calculated power losses, Pf,pad + Pf,space.4 (right).For the pad, three temperatures at r = 78.75mm in the circumferential direction were chosen.The results show good agreement between the measured and predicted values, especially in case of the eight pads with higher temperatures.The higher collar and pad heating for more pads could be observed.The pad temperature TA1 near the oil film inlet (small angle α), Figure 7, is lower for a reduced number of pads.Lower pad temperature TA1 relates to lower oil film inlet temperature Tin which is estimated by using the fresh oil replacing factor fr.      4 (right).For the pad, three temperatures at r = 78.75mm in the circumferential direction were chosen.The results show good agreement between the measured and predicted values, especially in case of the eight pads with higher temperatures.The higher collar and pad heating for more pads could be observed.The pad temperature TA1 near the oil film inlet (small angle α), Figure 7, is lower for a reduced number of pads.Lower pad temperature TA1 relates to lower oil film inlet temperature T in which is estimated by using the fresh oil replacing factor f r .
Table 5.The percentage value in the parentheses is the ratio heat to sum of both calculated power losses, Pf,pad + Pf,space.4 (right).For the pad, three temperatures at r = 78.75mm in the circumferential direction were chosen.The results show good agreement between the measured and predicted values, especially in case of the eight pads with higher temperatures.The higher collar and pad heating for more pads could be observed.The pad temperature TA1 near the oil film inlet (small angle α), Figure 7, is lower for a reduced number of pads.Lower pad temperature TA1 relates to lower oil film inlet temperature Tin which is estimated by using the fresh oil replacing factor fr.   Table 5.The percentage value in the parentheses is the ratio heat to sum of both calculated power losses, Pf,pad + Pf,space.4 (right).For the pad, three temperatures at r = 78.75mm in the circumferential direction were chosen.The results show good agreement between the measured and predicted values, especially in case of the eight pads with higher temperatures.The higher collar and pad heating for more pads could be observed.The pad temperature TA1 near the oil film inlet (small angle α), Figure 7, is lower for a reduced number of pads.Lower pad temperature TA1 relates to lower oil film inlet temperature Tin which is estimated by using the fresh oil replacing factor fr.

Bearing Power Loss
Measured bearing power loss P bear estimated with the oil temperature difference in the bearing housing ∆T = T out − T in and oil flow rate Q in is higher than the calculated sum of the power loss in the oil film gap and the space between the pads.The free collar rotational surfaces outside the range of outer and inner pad diameter, Figure 2, create additional rotating power loss P f,rot .The bearing tool does not consider free rotating surfaces, therefore in Table 6 the rotational power loss, P f,rot , is calculated by subtraction of the measured and calculated power losses: P f,rot = P bear − P f,pad − P f,space . ( GE Renewable Energy Hydro calculates similar values for rotational power loss, P f,rot by using an experienced formula.

Collar and Pad Heating Depending on Distance between the Pads
The influence of the collar and pad temperatures on the distance between the pads l distance is shown in Figure 8 for offset pivot and in Figure 9 for central pivot.The x-axis represents the surface factor f surf , the dimensionless ratio l pad to l arc with the pad length l pad = D mean π α pad /360 • and the arc length l arc = D mean π/N pad .The distance between the pads l distance = l arc − l pad .
does not consider free rotating surfaces, therefore in Table 6 the rotational power loss, Pf,rot, is calculated by subtraction of the measured and calculated power losses: Pf,rot = Pbear − Pf,pad − Pf,space. ( GE Renewable Energy Hydro calculates similar values for rotational power loss, Pf,rot by using an experienced formula.

Collar and Pad Heating Depending on Distance between the Pads
The influence of the collar and pad temperatures on the distance between the pads ldistance is shown in Figure 8 for offset pivot and in Figure 9 for central pivot.The x-axis represents the surface factor fsurf, the dimensionless ratio lpad to larc with the pad length lpad = Dmean π αpad/360° and the arc length larc = Dmean π/Npad.The distance between the pads ldistance = larc -lpad.
The heating of the collar is shown with the temperature at radius r = 78.75mm (sensor TSRz).For the pad, two locations at the same radius were chosen near the leading and trailing edge at the angles α = 3.2° and 34.8° (sensor TA1 and TA4).The single quadrats in the graphics represent the measured sensor temperatures.
There is an increase mainly of collar temperature TSRz for an increasing surface factor fsurf.For small number of pads (low surface factor) the collar temperature TSRz is closer to the pad temperature at oil film inlet side TA1 and increases to pad temperature at oil film outlet side TA4 for higher number of pads, as shown in [4].There are the same trends of the curves for offset and central pivot.The three power loss components in oil film, space and on free collar rotational surfaces, Pf,pad,, Pf,space, Pf,rot, are shown in Figure 10 for the offset pivot and in Figure 11 for the central pivot depending on the surface factor fsurf.The low values of the space power loss Pf,space are significant.The measured bearing power loss Pbear is not significantly different between offset and central pivot.The lower losses due to higher temperatures are compensated by the higher losses due to smaller oil film gaps.The heating of the collar is shown with the temperature at radius r = 78.75mm (sensor TSRz).For the pad, two locations at the same radius were chosen near the leading and trailing edge at the angles α = 3.2 • and 34.8 • (sensor TA1 and TA4).The single quadrats in the graphics represent the measured sensor temperatures.
There is an increase mainly of collar temperature TSRz for an increasing surface factor f surf .For small number of pads (low surface factor) the collar temperature TSRz is closer to the pad temperature at oil film inlet side TA1 and increases to pad temperature at oil film outlet side TA4 for higher number of pads, as shown in [4].There are the same trends of the curves for offset and central pivot.
The three power loss components in oil film, space and on free collar rotational surfaces, P f,pad, , P f,space , P f,rot , are shown in Figure 10 for the offset pivot and in Figure 11 for the central pivot depending on the surface factor f surf .The low values of the space power loss P f,space are significant.The measured bearing power loss P bear is not significantly different between offset and central pivot.The lower losses due to higher temperatures are compensated by the higher losses due to smaller oil film gaps.The three power loss components in oil film, space and on free collar rotational surfaces, Pf,pad,, Pf,space, Pf,rot, are shown in Figure 10 for the offset pivot and in Figure 11 for the central pivot depending on the surface factor fsurf.The low values of the space power loss Pf,space are significant.The measured bearing power loss Pbear is not significantly different between offset and central pivot.The lower losses due to higher temperatures are compensated by the higher losses due to smaller oil film gaps.

Sensitivity Analysis
In the following chapter, the influence of different parameters, such as pad pivot position, rotational speed, thrust load, oil flow rate, on the bearing behavior and the operating condition of a high load is presented.

Sensitivity Analysis
In the following chapter, the influence of different parameters, such as pad pivot position, rotational speed, thrust load, oil flow rate, on the bearing behavior and the operating condition of a high load is presented.

Pad Pivot Position
For an eight-pad design at 4000 rpm and 3 MPa, the influence on both support systems, offset and central pivot, was investigated.Oil film pressure curves are different, as shown in Figure 12, with the consequence that the pad tilting for the central pivot is lower than for offset pivot.Lower pad tilting creates higher temperatures in space, pad and collar, as shown in Figure 13.The oil film inlet and outlet gap as well as temperature fields are shown in Figure 14.The replacing factor is higher for the central pivot, as shown in Table 2, and seems to depend on the ratio inlet to outlet oil film gap.The dependence on this ratio was also found in [2] for the hot oil carry-over coefficient.

Pad Pivot Position
For an eight-pad design at 4000 rpm and 3 MPa, the influence on both support systems, offset and central pivot, was investigated.Oil film pressure curves are different, as shown in Figure 12, with the consequence that the pad tilting for the central pivot is lower than for offset pivot.Lower pad tilting creates higher temperatures in space, pad and collar, as shown in Figure 13.The oil film inlet and outlet gap as well as temperature fields are shown in Figure 14.The replacing factor is higher for the central pivot, as shown in Table 2, and seems to depend on the ratio inlet to outlet oil film gap.The dependence on this ratio was also found in [2] for the hot oil carry-over coefficient.

Pad Pivot Position
For an eight-pad design at 4000 rpm and 3 MPa, the influence on both support systems, offset and central pivot, was investigated.Oil film pressure curves are different, as shown in Figure 12, with the consequence that the pad tilting for the central pivot is lower than for offset pivot.Lower pad tilting creates higher temperatures in space, pad and collar, as shown in Figure 13.The oil film inlet and outlet gap as well as temperature fields are shown in Figure 14.The replacing factor is higher for the central pivot, as shown in Table 2, and seems to depend on the ratio inlet to outlet oil film gap.The dependence on this ratio was also found in [2] for the hot oil carry-over coefficient.

Rotational Speed
The speed was reduced from 4000 to 2000 rpm at a specific pressure of 3 MPa.Calculated and measured pad and collar temperatures are shown for offset pivot and all number of pads: see Figure 15.

Rotational Speed
The speed was reduced from 4000 to 2000 rpm at a specific pressure of 3 MPa.Calculated and measured pad and collar temperatures are shown for offset pivot and all number of pads: see Figure 15.
All temperatures are lower in case of the reduced speed in good agreement with the measured data.The fresh oil replacing factor does not depend on operating condition.All temperatures are lower in case of the reduced speed in good agreement with the measured data.The fresh oil replacing factor does not depend on operating condition.

Thrust Load
To analyze thrust load sensitivity, the specific pressure was reduced from 3 to 2 MPa at a rotational speed of 4000 rpm.The calculated and measured pad and collar temperatures are shown for offset pivot and all investigated number of pads in Figure 16.All temperatures decrease with lower pressure in good agreement with measured data.

Oil Flow Rate
To analyze the sensitivity of the oil flow rate at 4000 rpm and 3 MPa on the test rig, Qin was increased from 15.1 to 22.0 L/min.This causes a lower increase of the oil temperature in the bearing housing (ΔT = Tout − Tin).The calculated pad temperatures are approximately 3 K lower as shown in

Thrust Load
To analyze thrust load sensitivity, the specific pressure was reduced from 3 to 2 MPa at a rotational speed of 4000 rpm.The calculated and measured pad and collar temperatures are shown for offset pivot and all investigated number of pads in Figure 16.All temperatures decrease with lower pressure in good agreement with measured data.All temperatures are lower in case of the reduced speed in good agreement with the measured data.The fresh oil replacing factor does not depend on operating condition.

Thrust Load
To analyze thrust load sensitivity, the specific pressure was reduced from 3 to 2 MPa at a rotational speed of 4000 rpm.The calculated and measured pad and collar temperatures are shown for offset pivot and all investigated number of pads in Figure 16.All temperatures decrease with lower pressure in good agreement with measured data.

Oil Flow Rate
To analyze the sensitivity of the oil flow rate at 4000 rpm and 3 MPa on the test rig, Qin was increased from 15.1 to 22.0 L/min.This causes a lower increase of the oil temperature in the bearing housing (ΔT = Tout − Tin).The calculated pad temperatures are approximately 3 K lower as shown in

Oil Flow Rate
To analyze the sensitivity of the oil flow rate at 4000 rpm and 3 MPa on the test rig, Q in was increased from 15.1 to 22.0 L/min.This causes a lower increase of the oil temperature in the bearing housing (∆T = T out − T in ).The calculated pad temperatures are approximately 3 K lower as shown in Figure 17.

Oil Flow Rate
To analyze the sensitivity of the oil flow rate at 4000 rpm and 3 MPa on the test rig, Qin was increased from 15.1 to 22.0 L/min.This causes a lower increase of the oil temperature in the bearing housing (ΔT = Tout − Tin).The calculated pad temperatures are approximately 3 K lower as shown in Figure 17.

High Thrust Load
High thrust load can occur for example during an emergency stop of a vertically-aligned hydro generator.For the eight-pad design with offset pivot, the thrust load was increased from 3 to 6 MPa at rotational speed of 4000 rpm, as shown in Figure 18.As this load value is outside the measured range of the test rig, only calculation results are shown by using the fresh oil replacing factor which is independent from operating condition.There is a decrease of oil film thickness and an increase of the temperature for higher load.

High Thrust Load
High thrust load can occur for example during an emergency stop of a vertically-aligned hydro generator.For the eight-pad design with offset pivot, the thrust load was increased from 3 to 6 MPa at rotational speed of 4000 rpm, as shown in Figure 18.As this load value is outside the measured range of the test rig, only calculation results are shown by using the fresh oil replacing factor which is independent from operating condition.There is a decrease of oil film thickness and an increase of the temperature for higher load.

Conclusions
Recalculation of measurements was carried out for a fully flooded tilting pad thrust bearing by using Combros-A calculation tool, which solves the Reynolds, energy and heating equations iteratively.The variation of pad number in a thrust bearing leads to different distances between the pads.Its influence on the thermal behavior is the focus of this publication.The overlapping factor, a dimensionless value, is used to describe the pad distance.An oil mixing model for the space between the pads is defined, and the fresh oil replacing factor depends on the bearing types.The fresh oil replacing factor changes for different pad distances but is independent from operating conditions such as speed and load.A good agreement between calculation and measurements was reached by using the adjusted values shown in Table 2.
Author Contributions: Michal Wodtke and Michal Wasilczuk designed the experiment and analyzed the experimental data; Hubert Schwarze and Christian Kraft developed the tool, Monika Chmielowiec-Jablczyk and Andreas Schubert analyzed the calculated data and wrote the paper.

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

Conclusions
Recalculation of measurements was carried out for a fully flooded tilting pad thrust bearing by using Combros-A calculation tool, which solves the Reynolds, energy and heating equations iteratively.The variation of pad number in a thrust bearing leads to different distances between the pads.Its influence on the thermal behavior is the focus of this publication.The overlapping factor, a dimensionless value, is used to describe the pad distance.An oil mixing model for the space between the pads is defined, and the fresh oil replacing factor depends on the bearing types.The fresh oil replacing factor changes for different pad distances but is independent from operating conditions such as speed and load.A good agreement between calculation and measurements was reached by using the adjusted values shown in Table 2.

Figure 1 .
Figure 1.Model of the oil mixing in the space between the pads.

Figure 1 .
Figure 1.Model of the oil mixing in the space between the pads.

Figures 6
Figures 6 and 7 show a comparison of measured and calculated temperatures for collar and pad.The radii of the collar, where the temperature was measured, are shown in Figure4(right).For the pad, three temperatures at r = 78.75mm in the circumferential direction were chosen.The results show good agreement between the measured and predicted values, especially in case of the eight pads with higher temperatures.The higher collar and pad heating for more pads could be observed.The pad temperature TA1 near the oil film inlet (small angle α), Figure7, is lower for a reduced number of pads.Lower pad temperature TA1 relates to lower oil film inlet temperature Tin which is estimated by using the fresh oil replacing factor fr.

Figures 6 and 7
Figures 6 and 7 show a comparison of measured and calculated temperatures for collar and pad.The radii of the collar, where the temperature was measured, are shown in Figure4(right).For the pad, three temperatures at r = 78.75mm in the circumferential direction were chosen.The results show good agreement between the measured and predicted values, especially in case of the eight pads with higher temperatures.The higher collar and pad heating for more pads could be observed.The pad temperature TA1 near the oil film inlet (small angle α), Figure7, is lower for a reduced number of pads.Lower pad temperature TA1 relates to lower oil film inlet temperature T in which is estimated by using the fresh oil replacing factor f r .

Figures 6 and 7
Figures 6 and 7 show a comparison of measured and calculated temperatures for collar and pad.The radii of the collar, where the temperature was measured, are shown in Figure4(right).For the pad, three temperatures at r = 78.75mm in the circumferential direction were chosen.The results show good agreement between the measured and predicted values, especially in case of the eight pads with higher temperatures.The higher collar and pad heating for more pads could be observed.The pad temperature TA1 near the oil film inlet (small angle α), Figure7, is lower for a reduced number of pads.Lower pad temperature TA1 relates to lower oil film inlet temperature Tin which is estimated by using the fresh oil replacing factor fr.

Figures 6 and 7
Figures 6 and 7 show a comparison of measured and calculated temperatures for collar and pad.The radii of the collar, where the temperature was measured, are shown in Figure4(right).For the pad, three temperatures at r = 78.75mm in the circumferential direction were chosen.The results show good agreement between the measured and predicted values, especially in case of the eight pads with higher temperatures.The higher collar and pad heating for more pads could be observed.The pad temperature TA1 near the oil film inlet (small angle α), Figure7, is lower for a reduced number of pads.Lower pad temperature TA1 relates to lower oil film inlet temperature Tin which is estimated by using the fresh oil replacing factor fr.

Figure
Figure Calculated oil film inlet/outlet gap and temperature fields, pivot: offset (a/b), central (c/d).

Table 2 .
Fresh oil replacing factor.

Table 3 .
Heat convection coefficient at the pad surfaces for calculation.

Table 4 .
Heat convection coefficient at the collar surfaces for calculation.

Table 5 .
Calculated power losses and heat flows: rotational speed 4000 rpm, specific pressure 3 MPa.

Table 5 .
Calculated power losses and heat flows: rotational speed 4000 rpm, specific pressure 3 MPa.
4.2.Comparison with Collar and Pad Temperature Measurement

Table 5 .
Calculated power losses and heat flows: rotational speed 4000 rpm, specific pressure 3 MPa.

Table 5 .
Calculated power losses and heat flows: rotational speed 4000 rpm, specific pressure 3 MPa.

Table 6 .
Rotational power loss P f,rot rotational speed 4000 rpm, specific pressure 3 MPa.