Location of Repository

High Order Finite Element Solution of Elastohydrodynamic Lubrication Problems

By Hongqiang Lu


In this thesis, a high-order finite element scheme, based upon the Discontinuous Galerkin (DG) method, is introduced to solve one- and two-dimensional Elastohydrodynamic Lubrication (EHL) problems (line contact and point contact). This thesis provides an introduction to elastohydrodynamic lubrication, including some history, and a description of the underlying mathematical model which is based upon a thin film approximation and a linear elastic model. Following this, typical nondimensionalizations of the equations are discussed, along with boundary conditions. Two families of problems are considered: line and point contacts. Following a review of existing numerical methods for EHL problems, a different numerical technique, known as the Discontinuous Galerkin method is described. This is motivated by the high accuracy requirement for the numerical simulation of EHL problems. This method is successfully applied to steady-state line contact problems. The free boundary is captured accurately using the moving-grid method and the penalty method respectively. Highly accurate numerical results are obtained at a low expense through the use of h-adaptivity methods based on discontinuity and high-order components respectively. Combined with the Crank-Nicolson method and other implicit schemes for the temporal discretization, highly accurate solutions are also obtained for transient line contact problems using the high order DG method for the spatial discretization. In particular, an extra pressure spike is captured, which is difficult to resolve when using low order schemes for spatial discretization. The extension of this high order DG method to the two-dimensional case (point contact) is straightforward. However, the computation in the two-dimensional case is more expensive due to the extra dimension. Hence p-multigrid is employed to improve the efficiency. Since the free boundary in the two-dimensional case is more complicated, only the penalty method is used to handle the cavitation condition. This thesis is ended with the conclusions and a discussion of future work

Publisher: School of Computing (Leeds)
Year: 2006
OAI identifier: oai:etheses.whiterose.ac.uk:1341

Suggested articles



  1. (1986). A forward iterative numerical method for steady-state elastohydrodynamically lubricated contacts at high loads.
  2. (1986). Apenaltyformulationandnumericalapproximationof theReynolds-Hertz problem of elastohydrodynamic lubrication.
  3. (1992). Elasto-Hydrodynamic Lubrication.
  4. (1994). Higher-order multilevel solvers for the EHL line and point contact problem.
  5. (1991). Multilevel Solution of the EHL Line and Point Contact Problems.
  6. (2001). Simpliļ¬ed multigrid technique for the numerical solution to the steady-state and transient EHL line contacts and the arbitrary entrainment EHL point contacts.
  7. (1994). Transient analysis of surface features in an EHL line contact in the case of sliding.

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.