Adaptive high-order finite element solution of transient elastohydrodynamic lubrication problems

This article presents a new numerical method to solve transient line contact elastohydrodynamic lubrication (EHL) problems. A high-order discontinuous Galerkin (DG) finite element method is used for the spatial discretization, and the standard Crank-Nicolson method is employed to approximate the time derivative. An h-adaptivity method is used for grid adaptation with the time-stepping, and the penalty method is employed to handle the cavitation condition. The roughness model employed here is a simple indentation, which is located on the upper surface. Numerical results are presented comparing the DG method to standard finite difference (FD) techniques. It is shown that micro-EHL features are captured with far fewer degrees of freedom than when using low-order FD methods.