AT Advances in Tribology 1687-5923 1687-5915 Hindawi 10.1155/2018/9480636 9480636 Research Article Recirculation Flow and Pressure Distributions in a Rayleigh Step Bearing http://orcid.org/0000-0003-2429-2899 Shen Feng 1 Yan Cheng-Jin 1 Dai Jian-Feng 1 http://orcid.org/0000-0002-3408-595X Liu Zhao-Miao 1 Ciulli Enrico College of Mechanical Engineering and Applied Electronics Technology & Institute for Advanced Mechanics in Engineering Beijing University of Technology Beijing 100124 China bjut.edu.cn 2018 2162018 2018 26 01 2018 16 04 2018 26 04 2018 2162018 2018 Copyright © 2018 Feng Shen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Flow characteristics in the Rayleigh step slider bearing with infinite width have been studied using both analytical and numerical methods. The conservation equations of mass and momentum were solved utilizing a finite volume approach and the whole flow field was simulated. More detailed information about the flow patterns and pressure distributions neglected by the Reynolds lubrication equation has been obtained, such as jumping phenomenon around a Rayleigh step, vortex structure, and shear stress distribution. The pressure distribution of the Rayleigh step bearing with optimum geometry has been numerically simulated and the results obtained agreed with the analytical solution of the classical Reynolds lubrication equation. The simulation results show that the maximum pressure of the flow field is at the step tip not on the lower surface and the increment of the strain rate from Navier-Stokes equation is approximately 49 percent greater than that from Reynolds theory at the step tip. It is also shown that the position of the maximum pressure of the lower surface is a little less than the length of the first region. These results neglected by the Reynolds lubrication equation are important for designing a bearing.

National Natural Science Foundation of China 11572013 Beijing Municipal Natural Science Foundation 7152012 Beijing Municipal Education Commission KM201610005002
1. Introduction

Rayleigh step bearing has been widely used in industry due to its highest load capacity among all other possible bearing geometries. Many researches on improving its load capacity were carried out using an analytical method by solving the classical Reynolds lubrication equation, assuming that the bearing length should be at least 100 times of the film thickness. In 1918, the theory of a step bearing was firstly discussed by Lord Rayleigh , determining the optimum geometry with maximum load capacity per unit width for a given film thickness and bearing length. This configuration is now referred to as Rayleigh step bearing.

Because of its high load capacity and cheap manufacturing, the Rayleigh step bearing has been widely used in industry, such as thrust and pad bearings . In addition, a series of Rayleigh steps are used in journal bearings to form a grooved bearing with higher performance. Since then, researches into optimum design and fluid dynamic characteristics of this bearing have attracted much attention [6, 7]. Rahmani et al.  comprehensively studied the Rayleigh step slider bearing including the effect of variations of pressure at the boundaries on the optimum parameters. The bearing is also optimized considering the lubricant flow rate, friction force, and friction coefficient. Auloge et al.  studied the optimum design of Rayleigh step bearing and determined the relationships between step location and height along with non-Newtonian lubricants. Artiles et al.  described the analysis and design of a 50mm diameter floating-ring seals with the Rayleigh step lift pads, considering influences of the surface speed, gas pressure, and inertia. Hong et al.  investigated the flow of a Newtonian fluid and Bingham fluid in a Rayleigh step bearing.

Through the numerical solutions of the conservation equations of mass and momentum, they found some special jumping phenomena around a Rayleigh step, which are important for designing a bearing and the study of wear characteristics. Zhu  investigated both numerically and experimentally the response of a rotor supported on Rayleigh step gas bearing using Galerkin finite element method. Constantinescu et al.  analyzed the pressure variation due to fluid inertia effects in Rayleigh step bearings. Faria and San Andrés  calculated the bearing load capacity, static stiffness coefficients, and frequency-dependent force coefficients for the gas-lubricated plane and Rayleigh step slider bearings using both finite element and finite difference methods. Lee and Kim  calculated the air film temperature of Rayleigh step air foil thrust bearing by solving the Reynolds equation and 3D energy equation with thermohydrodynamic boundary conditions at the top foil, thrust disc, and cooling air plenum. By using the non-Newtonian fluid as a lubricant, the pressure distribution and load capacity for Rayleigh step bearing have also been investigated [11, 1619].

In previous literature, most researchers used an analytical method by solving the Reynolds lubrication equation, which is derived from fundamental equations governing fluid flows and based on several simplifying assumptions. Thus, the solutions of Reynolds lubrication equation for the pressure distribution, load capacity, and frictional coefficient were obtained through a relatively simple calculation and some detailed flow information on Rayleigh step bearing has been neglected . Only a few investigations used the calculation of the Navier-Stokes equations with finite volume method  or infinite element method [12, 14]. In recent years, the full Navier-Stokes equations have been increasingly applied to solve some lubrication problems . The trend is driven by a growing interest in describing lubricant flow behavior within a whole lubricated system. The Reynolds equation is confined to just the contact region, but the Navier-Stokes equations can be applied throughout the flow field [24, 25]. In order to optimize Rayleigh step bearing design scientifically, it is beneficial to obtain detailed flow information of the whole lubricant flow field .

In this study, the flow characteristics in the infinitely wide Rayleigh step bearing have been investigated by solving Navier-Stokes equations and comparing the numerical results with those results obtained by analytical solutions deduced from the Reynolds equation. From the comparison, more detailed flow characteristics and pressure and shear stress distributions have been obtained, and there are some new differences between the results obtained by two different methods. Moreover, the effects of step geometry and velocity of the lower surface on the recirculation flow and pressure and shear stress distributions in the whole flow field are investigated for better control of these bearings.

2. Analytical Solution

The geometry of a two-dimensional infinite Rayleigh step slider bearing is shown in Figure 1. There are two opposed surfaces: the upper surface separates the bearing into two regions and the lower surface has a pure tangential sliding motion relative to the upper surface. The one-dimensional Reynolds equation for a steady flow of an incompressible fluid of constant viscosity can be expressed in the following form :(1)xh3xpxx=6μU1hxxHere p is the oil film pressure, h is the oil film thickness, μ is the fluid dynamic viscosity, U1 is the velocity of the lower surface, and x is the Cartesian position variable.

The Rayleigh step bearing.

With reference to Figure 1, the film thickness is given by(2)hx=hii=1,0xL1hx=hii=2,L1<xLHere L is the characteristic length of the bearing, given by(3)L=L1+L2

The dimensionless scales of the Rayleigh step bearing can be defined as(4)ε=L1L;0<ε<1and(5)ξ=h1h2;ξ>1

The boundary conditions for the velocity and the pressure are(6)u=U1,v=0aty=0u=0,v=0aty=hii=1,2and(7)px=c=p0px=a=pa-pcHere p0 is the ambient pressure. Double integrating of (1) with respect to x, the pressure distribution along the bearing surface can be expressed as(8)px=6μU1ac1hx2dx+c1ac1hx3dx+c2Here a and c are the coordinates of the bearing leading and trailing ends (see Figure 1); c1 and c2 are integration constants. Substituting boundary conditions (6) and (7) into (8), the integration constants were obtained as(9)c1=pa-pc-6μU1ac1/hx2dxac1/hx3dxc2=0

Employing film thickness (2), the pressure distribution for the two regions of the bearing can be expressed in the following forms, respectively:(10)pxxa,b=6μU1x-ah12+c1x-ah13,pxxb,c=6μU1L1h12+x-bh22+c1L1h13+x-bh23

In addition, substituting (2) into (9), the integration constant c1 was obtained as follows:(11)c1=pa-pc-6μU1L1/h12+L2/h22L1/h13+L2/h23

Using pressure distribution (10), the load capacity of the bearing was calculated as(12)w=acpxdx=3μU1L1L1+2L2h12+L22h22+c12L1L1+2L2h13+L22h23

When there is no pressure difference (pa=pc=p0) between two ends of the bearing, the load capacity of the bearing is entirely supported by the hydrodynamic pressure distribution inside the bearing and the hydrostatic pressure support no longer exists. Equations (10), (11), and (12) will be reduced to a simpler form.(13)c1=-6μU1h1L1+L2h12/h22L1+L2h13/h23

The derivation of (10) is constant:(14)dpxdxxa,b=6μU1h12+c1h13>0,dpxdxxb,c=6μU1h22+c1h23<0

So the maximum pressure pm of the lower surface is at x=b; substituting (4), (5), and (13) into (10), pm was obtained as(15)pm=6μU1La-1ε1-εa2h12ε+1-εa3where the total load on unit width is(16)w=12Lpm=3μU1L2a-1ε1-εa2h12ε+1-εa3

3. Governing Equation and Numerical Methods

For steady state incompressible flow, the momentum and continuity equations of film region are(17)uux+vuy=-1ρpx+μρ2ux2+2uy2(18)uvx+vvy=-1ρpy+μρ2vx2+2vy2(19)ux+vy=0Here ρ is the fluid density, u and v are the velocity, μ is the fluid viscosity, and p is pressure. The pressure at the inlet and outlet boundaries of the Rayleigh step bearing was set to be atmospheric (105Pa) and zero velocity gradient in the direction normal to sliding was assumed. This defines a fully developed flow approximation through these boundaries and it is important to set the boundaries relatively far from the region of interest so that they do not influence the numerical solution. The entire domain was fully flooded. At the solid walls, the “no-slip” boundary condition was assumed for the momentum equations. The film region was considered to be a fully developed thin film flow. The bearing working condition is specified by the velocity of the moving wall (U1), varied from 1m/s to 8m/s. The fluid is considered as being incompressible and isothermal with μ=0.188Pa·s and ρ=850kg/m3. This approach was adopted because the present work is the first part of a study aimed at enlightening the flow characteristics in the Rayleigh step bearing. Meshes with 210295, 334793, and 5604066 cells were tested. The final refined mesh of the Rayleigh step bearing was built using 334793 cells and the solutions are independent of the refined mesh.

To solve the above equations, a CFD software package CFD-ACE+ (v2010, ESI Corporation) used for multiphysics computational analysis was employed. A finite volume method is utilized to turn the governing partial differential equations (PDE) into a system of algebraic equations and numerically integrated over each of the computational cells using a collocated cell-centered variable arrangement [28, 29]. A second-order upwind scheme is used for the momentum equations. For pressure calculations, SIMPLEC scheme  was used. The equation for pressure correction is obtained from the continuity equation, and the scheme of velocity and pressure calculations is fundamentally iterative in nature. The residuals of continuity and momentum equations are required to be below 10-5.

4. Results and Discussion

The dynamics of flow through a Rayleigh step was investigated. The simulations visualized both the velocity and pressure profiles, while varying system parameters including the depth and length of the step, inlet conditions (pressure and velocity), and shear speed.

4.1. Effects of Length Ratios on Pressure

The geometry and operation conditions are h1 =250μm, L=12.5mm, h2 =133.976μm, l1=8.975mm, U1=1m/s, and μ=0.188Pas, which are the same as those reported in a previous publication . From the analytic solution of Reynolds lubrication theory, the optimum values of ε and ξ according to maximum load for a Newtonian fluid are ε=0.718 and ξ=1.866 , respectively. The maximum pressure pm of 54KPa can be obtained with the analytic expression of (14). Figure 2 shows the pressure curves on upper and lower surfaces. The maximum pressure (pm) in Figure 2(a) is 53.7KPa close to the analytical value of 54KPa, which is the same as what is reported in a previous publication . Figure 2(b) shows that pressure distribution on the upper surface changes approximately 5.5 percent at the point of x=b. The reason is that there is a vortex that appeared near the step and the pressure distribution of flow field is changed significantly, which will be discussed in next section. The results were compared with the simulations of other CFD software (ANSYS-FLUENT, V15.0) in the Supplementary Materials (available here).

Pressure distributions on lower (a) and upper surfaces (b).

The pm values with different length ratios (ε=0.56, 0.64, 0.718, 0.8, and 0.88) are plotted in Figure 3. The maximum pressure pm can be obtained for ε=0.718 by both analytical and numerical methods. Figure 3 shows that the analytical and numerical results are similar and the validity of the numerical results is further verified.

p m values with different length ratios L1/L.

4.2. Effects of Height Ratios on Pressure

The effects of height ratios on the maximum pressure were simulated. The Rayleigh step parameters are h1=250μm, L=12.5mm, ε=0.718, U1=1m/s, and μ=0.188Pas. The pm values with different height ratios (ξ=1.56, 1.786, 1.866, 2.273, 2.778, 3.57, 5, and 9.33) are plotted in Figure 4. Figure 4 shows that the maximum pressure pm increases evidently as the value of ξ increases, and both the analytic and numerical results are similar.

p m values with different height ratios h1/h2.

Pressure contours around the steps with different height ratios (ξ=1.866 and ξ=9.33) are plotted in Figure 5. Figure 5 shows that the pressure values increase evidently as the value of ξ increases. It is found that the maximum pressure of the flow field is at the step tip instead of being on the lower surface. The maximum pressures of Figures 5(a) and 5(b) are 57KPa and 154.6KPa, which are higher than the pm values of the lower surface: 53.7KPa and 152KPa. These details of pressure distribution are not predicted if the Reynolds lubrication equation is used.

Pressure contours (Pa) of ξ=1.866 (a) and ξ=9.33 (b).

4.3. Vortex and Flow Details

The flow characteristics of Rayleigh step bearing with different geometric structures were studied. Some new phenomena like the vortex and reversed flow zone which are not predicted by the Reynolds lubrication equation were investigated. Figure 6 shows the streamlines of different geometric structures in parameters of U1=1m/s and μ=0.188Pas. In Figure 6(a), a reversed flow zone appears near the step with ε=0.56. When ε=0.718 and 0.88, the reversed flow zone disappeared and vortices are formed, and the vortices become smaller with the increase of length ratios (ε). In Figure 6(b), two reversed flow zones appear with ξ=3.57 and 2.27, and the center of the reserved flow zone moves rightward with the decrease of the height ratios (ξ). When ξ=1.7857, a vortex appears near the step. In the reversed flow zone, the velocity direction near the upper surface is opposite to the velocity of the lower surface U1.

Streamlines of different geometric structures.

The simulation results show that the pressure and shear stress distributions of the Rayleigh step bearing are significantly influenced by the vortex and velocity field. In Figure 7, the strain rates of fluid in Rayleigh step bearing (ε=0.718 and ξ=1.866) are changed obviously near the surfaces. Near the upper surface, the strain rate of fluid in region I is almost zero while that in region II is very high, as shown in Figure 7(a). It is found that the highest value of strain rate appears near the step tip. The strain rate from Navier-Stokes equation is approximately 49 percent greater than that from Reynolds theory at step tip and too strong to be omitted. This jumping phenomenon is not predicted by the Reynolds lubrication equation. Near the lower surface, the strain rate of fluid in region I is higher than that in region II. These results can supply more accurate boundary condition for a better solid analysis of the step tip.

Strain rates of fluid near the upper (a) and lower (b) surfaces.

4.4. Effects of <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M89"><mml:mrow><mml:msub><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> on Pressure

The velocity of the lower surface U1 has important effects on the pressure distribution. Figure 8 shows the lower surface pressures with different U1 values in a Rayleigh step bearing with optimum geometry. The pressure values increase evidently as the value of U1 increases. The maximum pressure (pm) of the lower surface was at x=8.8269 mm a little less than L1=b, which can be explained in Figure 5. A linear increase of pm values of the lower surface was noted with increasing U1 values, shown in Figure 9.

Lower surface pressures of different U1 values.

The pm values of the lower surface of different U1 values.

5. Conclusions

The current work presents a parametric study of the Rayleigh step bearing with infinite width using two different approaches: (a) analytical solution of the Reynolds equation; (b) the numerical solution of the Navier-Stokes equations by CFD method. The numerical results agree well with the analytical results with parameters of ε=0.718 and ξ=1.866. It was shown that the optimum length rate is ε=0.718 with no pressure difference (pa=pc=0) between two ends of the bearing, and the maximum pressures (pm) of the lower surface increase evidently as the value of height ratio increases. More detailed information on flow characteristics neglected by Reynolds lubrication equation, such as vortices, reversed flow, shear stress, and pressure distributions, has been obtained. It is found that the maximum pressure of the flow field is at the step tip instead of being on the lower surface, and the highest value of strain rate appears near the step tip. The velocity of the lower surface U1 has important effects on the pressure distribution. It is also found that the position of the maximum pressure (pm) of the lower surface is a little less than x=b. These results are helpful to supply more accurate flow details in the Rayleigh step bearing when designing an optimum Rayleigh step bearing geometry.

Nomenclature a :

Coordinate of bearing leading end (mm)

a n b :

Link coefficients

b :

Coordinate of region I leading end (mm)

c :

Coordinate of bearing trailing end (mm)

c 1 :

Integration constant

c 2 :

Integration constant

h 1 :

Height of region I (μm)

h 2 :

Height of region II (μm)

L :

Length of the bearing (mm)

L 1 :

Length of region I (mm)

L 2 :

Length of region II (mm)

p :

Pressure of fluid (Pa)

p a :

Pressure at bearing leading end (Pa)

p c :

Pressure at bearing trailing end (Pa)

p m :

The maximum pressure of the lower surface at x=b (Pa)

p ( x ) :

Pressure distribution along the bearing lower surface (Pa)

U 1 :

Velocity of the lower surface (m/s)

u :

Velocity in the x direction (m/s)

v :

Velocity in the y direction (m/s)

w :

Load capacity of the bearing (N)

μ :

Viscosity (Pa.s)

ρ :

Fluid density (Kg/m3)

ε :

Dimensionless length (L1/L)

ξ :

Dimensionless height (h1/h2).

Data Availability

All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of the paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (11572013), Beijing Municipal Natural Science Foundation (7152012), and Development Project of Beijing Municipal Education Commission (KM201610005002).

Supplementary Materials

Description of each file of Supplementary Materials: Verification of the results compared with the ANSYS-FLUENT simulations. The results of pressure distributions in Figures 2 and 8 and the results of streamline patterns in Figure 6(a) are simulated and compared with the FLUENT simulations, finding the same results. Figure 1: FLUENT results of pressure distributions on lower surfaces under the same conditions as that in Figure 2(a), finding the deviation of the maximum pressure is 2.4%. Figure 2: FLUENT results of pressure distributions on upper surfaces corresponding to the result of Figure 2(b), finding similar pressure distributions. Figure 3: FLUENT results of lower surface pressures of different U1 values corresponding to the result of Figure 8, finding the same pressure distributions. Figure 4: streamlines in a Rayleigh step bearing with optimum geometry at different U1 values corresponding to Figure 6(a), finding similar flow phenomena.

Rayleigh L. Notes on the theory of lubrication Philosophical Magazine 1918 35 1 12 Sato Y. Knight J. D. Performance characteristics of shrouded rayleigh-step and spiral groove viscous pumps Journal of Tribology 1992 114 3 499 503 2-s2.0-0026899330 10.1115/1.2920911 Boldyrev Y. Y. Lupulyak S. V. Shinder Y. K. Numerical solution of the variational Rayleigh problem of gas lubrication theory Journal of Fluid Dynamics 1995 30 6 826 832 2-s2.0-0029169642 10.1007/BF02078196 Zbl0860.76074 Farmer D. G. Shepherd J. J. Slip flow in the gas-lubricated Rayleigh step-slider bearing Internal Journal of Applied Mechanics and Engineering 2006 11 593 608 Dobrica M. B. Fillon M. About the validity of Reynolds equation and inertia effects in textured sliders of infinite width Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology 2009 223 1 69 78 2-s2.0-61649101519 10.1243/13506501JET433 Oike M. Nosaka M. Watanabe Y. Kikuchi M. Kamijo K. Experimental study on high-pressure gas seals for a liquid oxygen turbopump Tribology Transactions 1988 31 1 91 97 2-s2.0-0023846740 10.1080/10402008808981803 Lipschitz A. Basu P. Johnson R. P. A bi-directional gas thrust bearing Tribology Transactions 1991 34 1 9 16 2-s2.0-0026003647 10.1080/10402009108982003 Rahmani R. Shirvani A. Shirvani H. Analytical analysis and optimisation of the Rayleigh step slider bearing Tribology International 2009 42 5 666 674 2-s2.0-60649096670 10.1016/j.triboint.2008.09.002 Auloge J. Y. Bourgin P. Gay B. The optimum design of one-dimensional bearings with non-Newtonian lubricants Journal of Lubrication Technology 1983 105 391 395 Artiles A. Shapiro W. Jones H. F. Design analysis of Rayleigh-step floating-ring seals Tribology Transactions 1984 27 321 331 Hong Y. Chen D. Wang J. Li X. Numerical modeling of a Newtonian and a Bingham fluid in a Rayleigh step bearing Journal of Tribology 2003 125 1 206 210 2-s2.0-0038585106 10.1115/1.1506321 Zhu X. Experimental response of a rotor supported on Rayleigh step gas bearings [M. S. thesis] 2004 College Station, TX, USA Texas A&M University Constantinescu V. N. Galetuse S. Pressure Drop Due to Inertia Forces in Step Bearings Journal of Lubrication Technology 1976 98 1 167 10.1115/1.3452757 Faria M. T. C. San Andrés L. On the numerical modeling of high-speed hydrodynamic gas bearings Journal of Tribology 2000 122 1 124 130 2-s2.0-0033627412 10.1115/1.555335 Lee D. Kim D. Three-dimensional thermohydrodynamic analyses of rayleigh step air foil thrust bearing with radially arranged bump foils Tribology Transactions 2011 54 3 432 448 2-s2.0-84995625035 10.1080/10402004.2011.556314 Yagi K. Sugimura J. Elastohydrodynamic simulation of Rayleigh step bearings in thin film hydrodynamic lubrication Tribology International 2013 64 204 214 2-s2.0-84877338062 10.1016/j.triboint.2013.04.005 Elkouh A. F. Yang D.-F. Flow of a power-law fluid in a rayleigh step Journal of Tribology 1991 113 3 428 433 2-s2.0-0026191282 10.1115/1.2920642 Leek T. H. Lingard S. Atkini R. J. Bullough W. A. An experimental investigation of the flow of an electro-rheological fluid in a rayleigh step bearing Journal of Physics D: Applied Physics 1993 26 10 1592 1600 2-s2.0-0027684516 10.1088/0022-3727/26/10/009 Naduvinamani N. B. Non-Newtonian effects of second-order fluids on double-layered porous Rayleigh-step bearings Japan Society of Fluid Mechanics: Fluid Dynamics Research 1997 21 6 495 507 10.1016/S0169-5983(97)00019-1 MR1600816 Zbl1051.76542 Dobrica M. Fillon M. Reynolds' model suitability in simulating rayleigh step bearing thermohydrodynamic problems Tribology Transactions 2005 48 4 522 530 2-s2.0-33646408893 10.1080/05698190500385088 Braun M. J. Dzodzo M. B. Three-Dimensional Flow and Pressure Patterns in a Hydrostatic Journal Bearing Pocket Journal of Tribology 1997 119 4 711 719 2-s2.0-0031445368 10.1115/1.2833875 Helene M. Arghir M. Frene J. Numerical three-dimensional pressure patterns in a recess of a turbulent and compressible hybrid journal bearing Journal of Tribology 2003 125 2 301 308 2-s2.0-0037564167 10.1115/1.1506327 Kim J.-S. Kim K.-W. Effects of distance between pads on the inlet pressure build-up on pad bearings Journal of Tribology 2002 124 3 506 514 2-s2.0-0036629107 10.1115/1.1456452 Brajdic-Mitidieri P. Gosman A. D. Ioannides E. Spikes H. A. CFD analysis of a low friction pocketed pad bearing Journal of Tribology 2005 127 4 803 812 2-s2.0-27544473425 10.1115/1.2032990 Radulescu A. V. Theoretical Calculus Model for the Finite Length Rayleigh Step Bearings Lubricated with Greases Proceedings of the BALTTRIB 2015 November 2015 Akademija, Lithuania Aleksandras Stulginskis University 10.15544/balttrib.2015.18 Kumar R. Azam M. S. Ghosh S. K. Khan H. Effect of surface roughness and deformation on Rayleigh step bearing under thin film lubrication Industrial Lubrication and Tribology 2017 69 6 1016 1032 2-s2.0-85038603230 10.1108/ILT-04-2017-0098 Hori Y. Hydrodynamic Lubrication 2006 Tokyo, Japan Springer-Verlag Deng D. F. A numerical and experimental investigation of Taylor flow instabilities in narrow gaps and their relationship to turbulent flow in bearings [Ph.D. thesis] 2007 Akron, Ohio, USA University of Akron Horvat F. E. Experimental and numerical analysis of flow and pressure fields inside a variable depth single pocket hydrostatic bearing [M. S. thesis] 2008 Akron, Ohio, USA University of Akron van Doormaal J. P. Raithby G. D. Enhancements of the SIMPLE Method for predicting incompressible fluid flows Numerical Heat Transfer, Part B: Fundamentals 1984 7 2 147 163 10.1080/01495728408961817 ZBL0553.76005 2-s2.0-0021412675 Zbl0553.76005