0
Research Papers: Hydrodynamic Lubrication

A New Singularity Treatment Approach for Journal-Bearing Mixed Lubrication Modeled by the Finite Difference Method With a Herringbone Mesh

[+] Author and Article Information
Yanfeng Han, Jiaxu Wang

State Key Laboratory
of Mechanical Transmission,
Chongqing University,
174 Shazhengjie, Shapingba District,
Chongqing 40044, China

Shangwu Xiong

Department of Mechanical Engineering, Northwestern University,
2145 Sheridan Road,
Evanston, IL 60208

Q. Jane Wang

Department of Mechanical Engineering,
Northwestern University,
2145 Sheridan Road,
Evanston, IL 60208;
State Key Laboratory
of Mechanical Transmission,
Chongqing University,
174 Shazhengjie, Shapingba District,
Chongqing 40044, China
e-mail: qwang@northwestern.edu

1Corresponding author.

Contributed by the Tribology Division of ASME for publication in the JOURNAL OF TRIBOLOGY. Manuscript received March 14, 2015; final manuscript received July 21, 2015; published online August 31, 2015. Assoc. Editor: Mihai Arghir.

J. Tribol 138(1), 011704 (Aug 31, 2015) (10 pages) Paper No: TRIB-15-1079; doi: 10.1115/1.4031138 History: Received March 14, 2015; Revised July 21, 2015

Steady-state mixed hydrodynamic lubrication of rigid journal bearing is investigated by using a finite difference form of the Patir–Cheng average Reynolds equation under the Reynolds boundary condition. Two sets of discretization meshes, i.e., the rectangular and nonorthogonal herringbone meshes, are considered. A virtual-mesh approach is suggested to resolve the problem due to the singularities of pressure derivatives at the turning point of the herringbone mesh. The effectiveness of the new approach is examined by comparing the predicted load with that found in the literature for a smooth-surface case solved in the conventional rectangular mesh. The effects of the skewness angles of symmetric and asymmetric herringbone meshes on the predicted parameters, such as load, friction coefficient, attitude angle, and maximum pressure, are investigated for smooth, rough, and herringbone-grooved bearing surfaces. It is found that the new approach helps to improve the computational accuracy significantly, as demonstrated by comparing the results with and without the treatment of the pressure derivative discontinuity although the latter costs slightly less computational time.

FIGURES IN THIS ARTICLE
<>
Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.

References

Pinkus, O. , and Sternlicht, B. , 1961, Theory of Hydrodynamic Lubrication, McGraw-Hill, New York.
Elrod, H. G. , and Adams, M. L. , 1974, “ A Computer Program for Cavitation and Starvation Problems,” Cavitation and Related Phenomena in Lubrication, Mechanical Engineering Publication, New York, pp. 37–41.
McCallion, H. , Yousif, F. , and Lloyd, T. , 1970, “ The Analysis of Thermal Effects in a Full Journal Bearing,” ASME J. Lub. Tech., 92(4), pp. 578–587. [CrossRef]
Dmochowski, W. , Dadouche, A. , and Fillon, M. , 2013, “ Finite Difference Method for Film Fluid Bearings,” Encyclopedia of Tribology, Vol. 2, Q. J. Wang and Y.-W. Chung , eds., Springer, New York, pp. 1137–1143.
Whipple, R. T. P. , 1949, “ Herringbone Pattern Thrust Bearing,” Paper No. AERE T/M 29.
Vohr, J. H. , and Chow, C. Y. , 1965, “ Characteristics of Herringbone Grooved Gas Lubricated Journal Bearings,” ASME J. Basic Eng., 87(3), pp. 568–578. [CrossRef]
Hirs, G. S. , 1965, “ The Load Capacity and Stability Characteristic of Hydrodynamic Grooved Journal Bearings,” ASLE Trans., 8(3), pp. 296–305. [CrossRef]
Kang, K. , Rhim, Y. , and Sung, K. , 1996, “ A Study of the Oil-Lubricated Herringbone-Grooved Journal Bearing—Part 1: Numerical Analysis,” ASME J. Tribol., 118(4), pp. 906–911. [CrossRef]
Zirkelback, N. , and San Andres, L. , 1998, “ Finite Element Analysis of Herringbone Groove Journal Bearings: A Parametric Study,” ASME J. Tribol., 120(2), pp. 234–240. [CrossRef]
Jang, G. H. , and Chang, D. I. , 1999, “ Analysis of a Hydrodynamic Herringbone Grooved Journal Bearing Considering Cavitation,” ASME J. Tribol., 122(1), pp. 103–109. [CrossRef]
Wu, J. K. , Li, A. F. , Lee, T. S. , Shu, C. , and Wan, J. M. , 2004, “ Operator-Splitting Method for the Analysis of Cavitation in Liquid-Lubricated Herringbone Grooved Journal Bearings,” Int. J. Numer. Methods Fluids, 44(7), pp. 765–775. [CrossRef]
Chao, P. C. P. , and Huang, J. S. , 2005, “ Calculating Rotor Dynamic Coefficients of a Ferrofluid-Lubricated and Herringbone-Grooved Journal Bearing Via Finite Difference Analysis,” Tribol. Lett., 19(2), pp. 99–109. [CrossRef]
Liu, J. , and Yoshihiro, M. , 2010, “ Analysis of Oil-Lubricated Herringbone Grooved Journal Bearing With Trapezoidal Cross-Section, Using a Spectral Finite Difference Method,” J. Hydrodyn., Ser. B, 22(5), pp. 408–412. [CrossRef]
Lee, T. S. , Liu, Y. G. , and Winoto, S. H. , 2004, “ Analysis of Liquid-Lubricated Herringbone Grooved Journal Bearings,” Int. J. Numer. Methods Heat Fluid Flow, 14(3), pp. 341–365. [CrossRef]
Vijayaraghavan, D. , and Keith, T. G. , 1990, “ Grid Transformation and Adaption Techniques Applied in the Analysis of Cavitated Journal Bearings,” ASME J. Tribol., 112(1), pp. 52–59. [CrossRef]
Vijayaraghavan, D. , and Keith, T. G. , 1990, “ An Efficient, Robust, and Time Accurate Numerical Scheme Applied to a Cavitation Algorithm,” ASME J. Tribol., 112(1), pp. 44–51. [CrossRef]
Patir, N. , and Chen, H. S. , 1978, “ An Average Flow Model for Determining Effects of Three-Dimensional Roughness on Partial Hydrodynamic Lubrication,” ASME J. Lubr. Technol., 100(1), pp. 12–17. [CrossRef]
Lee, S. C. , and Ren, N. , 1996, “ Behavior of Elastic-Plastic Rough Surface Contacts as Affected by Surface Topography, Load, and Material Hardness,” Tribol. Trans., 39(1), pp. 67–74. [CrossRef]
Wang, Q. , Shi, F. , and Lee, S. C. , 1997, “ A Mixed-Lubrication Study of Journal Bearing Conformal Contacts,” ASME J. Tribol., 119(3), pp. 456–461. [CrossRef]
Wu, C. , and Zheng, L. , 1989, “ An Average Reynolds Equation for Partial Film Lubrication With a Contact Factor,” ASME J. Tribol., 111(1), pp. 188–191. [CrossRef]
Christopherson, D. G. , 1941, “ A New Mathematical Method for the Solution of Film Lubrication Problem,” Proc. Inst. Mech. Eng., 146(1), pp. 126–135. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

(a) Illustration of a skewed discretization mesh and (b) the neighboring nodes of node (i, j) when β = 0 deg. x-y are the Cartesian coordinates, β indicates the skewness angle of the herringbone meshes, A is a node in λ1−φ1 skewed coordinate, B is a node in λ2−φ2 skewed coordinate, and C is the turning location.

Grahic Jump Location
Fig. 2

Node (i, j) and its neighbors. (a) Common approach without any singularity treatment and (b) new approach using virtual nodes, marked by cross, to treat the singularity. L and R mark the points neighboring virtual node (i ̃, J + 1) on the line of y=yJ+1 in the φ − λ2 grids.

Grahic Jump Location
Fig. 4

Total bearing load and attitude angle obtained in the rectangular mesh and herringbone meshes without singularity treatment at the turning points. β indicates the skewness angle of the herringbone meshes, with β = 0 for the rectangular mesh, and yc is for the turning location. (a) yC = 0.5L, symmetric, using 80 × 41 and 160 × 41 nodes and (b) yC = 0.25L and yC = 0.125L using 160 × 41 nodes.

Grahic Jump Location
Fig. 5

Predicted total bearing load and attitude angle using the rectangular mesh (β = 0) and the herringbone meshes with the first virtual-node interpolation technique. (a) yC = 0.5L, symmetric, using 160 × 41 nodes; (b) yC = 0.25L using 160 × 41 nodes; and (c) yC = 0.125L using 160 × 41 nodes.

Grahic Jump Location
Fig. 6

Predicted total bearing load and attitude angle using the rectangular mesh (β = 0) and herringbone meshes with the second and third virtual-node interpolation technique (yC = 0.5L). (a) Second virtual-node interpolation technique using 160 × 41 nodes and (b) third virtual-node interpolation technique using 160 × 41 nodes.

Grahic Jump Location
Fig. 7

Total load error using the rectangular mesh (β = 0) and herringbone meshes with 160 × 41 nodes (yC = 0.5L). (a) Without singularity treatment; (b) first virtual-node interpolation technique; and (c) comparison between without singularity treatment and virtual-node interpolation technique.

Grahic Jump Location
Fig. 8

Comparison of the total bearing load obtained by using the rectangular mesh (β = 0) and the symmetric herringbone mesh (yC = 0.5L) with the first virtual-node interpolation technique

Grahic Jump Location
Fig. 9

Predicted pressures for smooth bearing surface by using the rectangular and herringbone meshes (eccentricity ratios at 0.75 and 0.99)

Grahic Jump Location
Fig. 10

Numerical results for mixed lubrication obtained by using the rectangular and symmetric herringbone meshes, withthe new method enforced by the first singularity treatment technique for the latter. (a) Attitude angle and the total load; (b) asperity contact load and the average friction coefficient; and (c) comparison of total load error for without singularity treatment and first singularity treatment technique.

Grahic Jump Location
Fig. 11

Comparison of load capacity with Refs. [7] and [10]. The geometry and lubricant parameters are the same with Ref.[10].

Grahic Jump Location
Fig. 12

Herringbone grooves in the bearing surface

Grahic Jump Location
Fig. 13

Numerical results for the herringbone-grooved bearing surface (groove angle = mesh angle = 60 deg) obtained by using the rectangular and symmetric herringbone meshes, with the new method enforced by the first singularity treatment technique for the latter. (a) Maximum lubricant pressure and theminimum film thickness; (b) asperity contact load and the average friction coefficient; and (c) attitude angle and the total load. The β = 0 (rectangular mesh) and β = 60 deg (herringbone mesh) cases with yC = 0.5L are directly comparable, showing good agreement.

Grahic Jump Location
Fig. 14

Predicted film thickness and pressure distributions for the herringbone-grooved surface (groove angle = 60 deg) at eccentricity ratio = 0.75. (a) Rectangular mesh, β = 0 deg and yC=0.5L and (b) herringbone meshes with the first singularity treatment, β = 60 deg and yC=0.5L.

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In