Numerical Modeling of A Newtonian and A Bingham Fluid in A Rayleigh Step Bearing
YP Hong,DR Chen,JD Wang,XD Li
DOI: https://doi.org/10.1115/1.1506321
2002-01-01
Journal of Tribology
Abstract:Contributed by the Tribology Division for publication in the ASME JOURNAL OF TRIBOLOGY. Manuscript received by the Tribology Division January 3, 2002 revised manuscript received June 14, 2002. Associate Editor: M. M. Khonsari. The Reynolds lubrication equation is derived from fundamental equations governing fluid flows 1. There are several assumptions in Reynolds lubrication theory, one of which is on characteristic length ratio. Let the length scale in the “plane” of the film be denoted by Lxz and let Ly be the length scale across the film (Fig. 1(a)). We have Ly/Lxz⩽O10−2 for common bearing geometries. The characteristic length relationship in the inlet and outlet regions of a Rayleigh step bearing (Fig. 1(b)) is Ly≪Lxz. However, this relationship changes in the vicinity of the step as Ly=H1−H2 and Lxz→0. As a result, when ∂ρh/2∂x→∞, the behavior of Eq. (1) becomes paradoxical. (1)∂∂x ρh312η∂p∂x+∂∂z ρh312η∂p∂z=u ∂ρh2∂xThis paper presents an investigation of the flow of a Newtonian fluid in a Rayleigh step bearing using the Navier-Stokes equations. Calculation results show some important phenomena neglected by the Reynolds equation. These phenomena are important in designing a bearing, especially in the study of wear characteristics. The transport properties of a Bingham fluid in a Rayleigh step bearing are also investigated in this paper. As many papers suggest that electro-rheological fluids, which are helpful for bearings controls in mechanical engineering 23, can be considered as a continuum with Bingham plastic constitutive properties 45. Over the last 15–20 years, a great multitude of studies have been dedicated to Bingham fluids 456. However, for this purpose the thin layer lubrication theory with two-viscosity model is used in the general case 7. This paper presents results on a Bingham fluid with more fundamental equations, in which the effects of boundary curvature (a 90 deg step in the context) on the pressure change are described. The performance of viscoelastic model with regard to different characteristic length (aspect) ratios and different yield stresses are investigated. Detailed information near the step can be obtained for better control of these bearings. The steady flow of an incompressible fluid is governed by the conservation equations of mass and momentum (or Navier-Stokes equations) 8. (2)∇s˙ρv=0(3)∇s˙ρvv=−∇p+∇s˙τHere ρ is the fluid density, v is the velocity vector, p is pressure and τ is the deviatoric stress tensor. The constitutive relation of two-viscosity model 7 is (4a)(4b)τ=η1dτ:τ<2τ12η2d+τ01/2d:d1/2 dτ:τ⩾2τ12Here d is the rate of stain tensor. The unyielded material, described by Eq. (4a), is treated as a Newtonian liquid of viscosity η1. Eq. (4b) is the expression for a Bingham-type fluid above the yield point. There are two distinguished stresses, τ1 and τ0, and a dimensionless parameter ε=η2/η1. We have τ0=τ11−ε; when ε=1, we obtain a Newtonian model and as ε→0 we recover the Bingham model. Some papers adopt the following expression instead of Eq. (4): η=η1τ:τ<2τ12η2+τ01/2d:d1/2τ:τ⩾2τ12A finite volume method developed for environmental hydraulics calculations 9 was modified to solve lubrication problems. Variables were set on staggered grids. SIMPLEC 10 procedure (an improved semi-implicit method for pressure linked equations) was adopted. The governing Eq. (2) and Eq. (3) were written in terms of the following general transport equation of primitive variable shown in Eq. (5a). Convective terms were discretized by the power-law scheme 11. Algebraic expression after discretization on each control volume is given by Eq. (5b). (5a)∇s˙ρvϕ−Γ∇ϕ=S(5b)APϕP=∑nbAnbϕnb+SΔxΔyThe definitions of ϕ, Γ, and S for all the governing equations are given in Table 1. AP=∑nbAnb,Anb=Dnb max0,1−0.1|Fnb/Dnb|0.5+max−Fnb,0. Here Dnb and Fnb are diffusion and flux parameter respectively on grid faces; Δx and Δy are distance between central nodes and neighbor nodes. Algebraic equations are solved by ADI algorithm (alternating-direction implicit iteration) 12. The program checks for degree of freedom (velocities and pressure) convergence after every global iteration. Convergence criterion for cases in this context is: ∑i=1N|ϕik−ϕik−1|/∑i=1N|ϕik|⩽10−7, where ϕik is the variable value at node i of the current kth iteration. The summation is performed over all nodes (N), using the absolute values of the differences. A Rayleigh bearing of infinite width in the direction Oz is shown in Fig. 1(b). The geometry and operation conditions are characterized by H1,B0,s=B1/B0,a=H1/H2 and the translation velocity u0 of the lower plate. From the analytic solution of Reynolds lubrication theory, the optimum values of s and a corresponding to maximum load for a Newtonian fluid are: s=0.718 and a=1.866. The load on unit width is: 0.5B0pm, where pm is the maximum pressure. The pressure difference along direction Oy is omitted. The relationship between pm and bearing parameters is 13: (6)pm=6ηu0B0s1−sa−1a2H12[s+1−sa3]The values of the bearing parameters for the case of the Newtonian lubricant are: H1=2.5×10−4m,B0=50×H1,u0=1m/s,s=0.718,η=η1=η2=0.188Pas,τ0=0, and a=1.866. The maximum pressure pm=54.0KPa can be obtained with the analytic expression of Eq. (6). With Eqs. (2) and (3) used to model this Rayleigh bearing, pressure curves on up and lower plates are plotted in Fig. 2(a) and Fig. 2(b) respectively. Fig. 2b shows that pressure changes about 5.4 percent from point R to S along the y direction. This difference disappears in the inlet and outlet regions. It should be mentioned that, when the flow is still in a laminar condition, the pressure change from point R to S is approximately 5.7 percent for u0=10m/s.The results of shear stresses from Eq. (2) and (3) are shown in Fig. 3(a) and Fig. 3(b). In Fig. 3(a), τyx=η∂u/∂y on RR′ is almost zero while that on SS′ is very high. The shear stress τyx from Navier-Stokes equation is approximately 29 percent greater than that from Reynolds theory at step tip (S). In Fig. 3(b), τxy=η∂ν/∂x on RS is too strong to be omitted. There is a stress peak on RS, followed by a rapid and sharp drop, which is helpful to supply more accurate boundary condition for a better solid analysis of the step tip. These jumping phenomena are not predicted if the Reynolds lubrication equation is used. Results from the above calculation case show that with more fundamental equations, more exact phenomena can be obtained. But it is just a case of a Newtonian fluid. Calculations on a Bingham lubricant are as follows. The Rayleigh step parameters are: H1=2.5×10−4m,B0=50×H1,u0=1m/s,s=0.718. As pressure control of Electro-rheological lubricants is quite important, emphases are put on pressure results. The parameters for a Bingham lubricant are: η2=0.188Pas,η1=105η2≫η2, and τ0=148.7Pa. Pressure contours around the steps with different aspect ratios (a=1.866 and a=1.866×5) are plotted in Fig. 4. The maximum pressures for a Newtonian fluid, pm=54.0KPa and pm=144.2KPa, can be obtained for these aspect ratios with Eq. (6). As the yield stress is not great here, maximum pressures for this Bingham lubricant are MX=57.06KPa and MX=165.08KPa in Fig. 4(a) and Fig. 4(b). The film pressure near the step increases evidently as the value of a=H1/H2 increases. To further test the validity of the computation program, the case of a=1.866×5 was also calculated with commercial ANSYS software installed at Virtual Manufacturing Workstation of Tsinghua University. The maximum pressure from ANSYS software is 165824 Pa. The difference between our model and the finite element software is less than 0.46 percent for this case. ANSYS FEM software is powerful in solving different engineering problems, while our computer program can conveniently run on a personal computer as less memory is required. Electro-rheological fluids are often used under strong electric field and thus take on high yield stress. Pressure properties on yield stresses are to be discussed below. When H2 in Fig. 1(b) tends to zero, the hydrodynamic pressures generated using an electro-rheological fluid in a Rayleigh step bearing at the limiting condition of zero flow rate were studied by Leek 5. To analyze pressure characteristics, the dimensionless pressure gradient Pg is defined by Pg=dp/dxH12/η2u0. For a Bingham plastic fluid it is convenient to employ the Bingham number, Bn, defined by Bn=τ0H1/η2u0. At the limiting condition of zero net flow rate, the approximate expression between Pg and Bn from Leek is (7)Pg=6+2.69Bn−0.0253Bn2This relation for an incompressible Bingham fluid is used under the following assumptions: (i) the flow is fully developed, one-dimensional, laminar and isothermal; (ii) curvature effects are negligible. Equation (7) was obtained from the developed thin layer lubrication theory of two-viscosity model. With Eq. (7), it is easy to get the expression for a Newtonian fluid in the zero flow rate case, Pg=dp/dxH12/η2u0=6.Equations (2) and (3) are used to take curvature effects into account. For Bingham lubricants in a Rayleigh bearing of zero net flow rate (η1≫η2,η2=0.188Pas and a→∞), three yield stresses τ0=148.7Pa,τ0=446.1Pa, and τ0=1487Pa are considered. Pressure distributions of numerical results near the step are shown in Fig. 5. For these three yield stresses, pressure values at corner R (Fig. 1(b)) calculated by approximate expression (7) are 176.3, 204.8, and 302.9 KPa respectively. While pressure values of contour B at the same corner are 166.0, 226.1, and 482.7 KPa respectively in Fig. 5. The larger Bingham number (the higher yield stress here) leads to greater difference between numerical results and Eq. (7). For the case of τ0=1487Pa, the difference is 179.8 KPa, which is about 50 percent of the result obtained from the developed thin layer lubrication theory. In the inlet region, at least 1.5H1 distant from the step, the velocity profile of a Newtonian fluid is in a parabolic form along the lubrication film. For a Bingham plastic flow, a plug of unyielded fluid occurs where τ:τ⩽2τ02. Velocity profiles for zero flow rate of different τ0 are shown in Fig. 6(a). The larger Bingham number (the higher yield stress here) leads to greater velocity gradient or shear stress along the translating plate. For both Newtonian and Bingham lubricants, velocity vector near the step is apparently two-dimensional, shown in Fig. 6(b). Velocity along film thickness direction (Oy) is quite large in this recirculation zone. Jumping phenomena in Fig. 3 is apprehensible by analyzing velocity vectors in Fig. 6(b). The numerical solutions of the conservation equations of mass and momentum show some special jumping phenomena around a Rayleigh step. Compared with the results of commonly used Reynolds lubrication equation, the increment of shear stress is up to 29 percent for a Newtonian fluid, which may aggravate the wear process of the step tip. Much attention should be paid to this increment when designing a Rayleigh bearing. Pressure differences between the results of Navier-stokes equations of two-viscosity model and those of developed thin layer lubrication theory are investigated. These differences are great for Bingham flows with large aspect ratios and high yield stresses. For the case of τ0=1487Pa in a Rayleigh bearing of zero flow rate, 50 percent of the pressure value might be omitted if boundary curvature is not taken into account. Compared with the thin layer lubrication theory, finite volume model in this paper can present more accurate data for electro-rheological fluids in a bearing controlling system. The authors wish to acknowledge the help of Prof. Yang Peiran, Prof. Hu Yuanzhong and Prof. Kong Xianmei. This research is supported by the National Natural Science Foundation of China under the grant number 59875046. a= aspect ratio of step heights, a=H1/H2B0= total bearing width B1= bearing width of inlet region Bn= Bingham number Bn=τ0H1/η2u0d= rate of stain tensor H1= bearing height of inlet region H2= bearing height of outlet region MX= maximum pressure in numerical results Pg= dimensionless pressure gradient dp/dxH12/η2u0pm= maximum pressure in Reynolds (thin layer) lubrication theory s= ratio of step width s=B1/B0τ = deviatoric stress tensor τ0= yield stress of a Bingham fluid τ1= shear stress in 2-viscosity model ε = dimensionless parameter ε=η2/η1η = apparent viscosity η1= newtonian fluid viscosity in 2-viscosity model η2= plastic viscosity of a Bingham fluid u0= translation velocity of the lower plate