A Triangle Based Finite Volume Method for the Integration of Lubrication’s Incompressible Bulk Flow Equations

It is well known that for a reduced Reynolds number (Re * =ρVH/μ.H/L) greater than unity, inertia forces have a dominant effect in the transport equations, thus rendering the classical lubrication equation inapplicable. The so called bulk flow system of equations is then the appropriate mathematical model for describing the flow in bearing and seals operating at Re * ≥1. The difficulty in integrating this system of equations is that one has to deal with coupled pressure and velocity fields. Analytic methods have a very narrow application range so a numerical method has been proposed by Launder and Leschziner in 1978. It represents a natural extrapolation of the successful SIMPLE algorithm applied to the bulk flow system of equations. The algorithm used rectangular, staggered control volumes and represented the state of the art at that moment. In the present work we introduced a method using triangular control volumes. The basic advantage of triangles versus rectangles is that non rectangular domains can be dealt without any a priori limitation. The present paper is focused on the description of the discretized equations and of the solution algorithm. Validations for bearings and seals operating in incompressible, laminar and turbulent flow regime are finally proving the accuracy of the method.

[1]  V. Venkatakrishnan Convergence to steady state solutions of the Euler equations on unstructured grids with limiters , 1995 .

[2]  W. Kahlert Der Einfluß der Trägheitskräfte bei der hydrodynamischen Schmiermitteltheorie , 1948 .

[3]  A Quasi-Two-Dimensional Method for the Rotordynamic Analysis of Centered Labyrinth Liquid Seals , 1999 .

[4]  Luis San Andrés,et al.  Analysis of Variable Fluid Properties, Turbulent Annular Seals , 1991 .

[5]  C. Rhie,et al.  A numerical study of the turbulent flow past an isolated airfoil with trailing edge separation , 1982 .

[6]  R. Glowinski,et al.  FINITE ELEMENT METHODS FOR NAVIER-STOKES EQUATIONS , 1992 .

[7]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[8]  J. Murthy,et al.  A PRESSURE-BASED METHOD FOR UNSTRUCTURED MESHES , 1997 .

[9]  S. Majumdar Role of underrelaxation in momentum interpolation for calculation of flow with nonstaggered grids , 1988 .

[10]  P. L. George,et al.  Automatic Mesh Generation: Application to Finite Element Methods , 1992 .

[11]  V. Constantinescu Laminar viscous flow , 1995 .

[12]  J. Donea Generalized Galerkin Methods for Convection Dominated Transport Phenomena , 1991 .

[13]  Michael A. Leschziner,et al.  Flow in Finite-Width, Thrust Bearings Including Inertial Effects: I—Laminar Flow , 1978 .

[14]  S. Patankar Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[15]  Vivette Girault,et al.  Finite Element Methods for Navier-Stokes Equations - Theory and Algorithms , 1986, Springer Series in Computational Mathematics.

[16]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[17]  G. G. Hirs A Bulk-Flow Theory for Turbulence in Lubricant Films , 1973 .

[18]  Y. Lai AN UNSTRUCTURED GRID METHOD FOR A PRESSURE-BASED FLOW AND HEAT TRANSFER SOLVER , 1997 .

[19]  M. Peric A finite volume method for the prediction of three-dimensional fluid flow in complex ducts , 1985 .