Adaptive high-order Discontinuous Galerkin solution of elastohydrodynamic lubrication point contact problems

This paper describes an adaptive implementation of a high order Discontinuous Galerkin (DG) method for the solution of Elastohydrodynamic Lubrication (EHL) point contact problems. These problems arise when modelling the thin lubricating film between contacts which are under sufficiently high pressure that the elastic deformation of the contacting elements cannot be neglected. The governing equations are highly non-linear and include a second order partial differential equation that is derived via the thin-film approximation. Furthermore, the problem features a free boundary, which models where cavitation occurs, and this is automatically captured as part of the solution process. The need for spatial adaptivity stems from the highly variable length scales that are present in typical solutions. Results are presented which demonstrate both the effectiveness and the limitations of the proposed adaptive algorithm.

[1]  R. Bosma,et al.  Multigrid, An Alternative Method for Calculating Film Thickness and Pressure Profiles in Elastohydrodynamically Lubricated Line Contacts , 1986 .

[2]  S. R. Wu.  A penalty formulation and numerical approximation of the Reynolds-Hertz problem of elastohydrodynamic lubrication , 1986 .

[3]  D. Dowson,et al.  A Numerical Solution to the Elasto-Hydrodynamic Problem , 1959 .

[4]  C. Venner.  Multilevel solution of the EHL line and point contact problems , 1991 .

[5]  Henry Peredur Evans,et al.  A novel method for integrating first- and second-order differential equations in elastohydrodynamic lubrication for the solution of smooth isothermal, line contact problems , 1999 .

[6]  Martin Berzins,et al.  Parallelization and scalability issues of a multilevel elastohydrodynamic lubrication solver , 2007, Concurr. Comput. Pract. Exp..

[7]  Seppo Santavirta.  Biotribology , 2005, Acta orthopaedica.

[8]  O. Reynolds.  I. On the theory of lubrication and its application to Mr. Beauchamp tower’s experiments, including an experimental determination of the viscosity of olive oil , 1886, Proceedings of the Royal Society of London.

[9]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[10]  G. Evans.  Practical Numerical Integration , 1993 .

[11]  Peter K. Jimack,et al.  On the use of adjoint-based sensitivity estimates to control local mesh refinement , 2009 .

[12]  D. Eyheramendy,et al.  A Full-System Approach of the Elastohydrodynamic Line/Point Contact Problem , 2008 .

[13]  Jing Wang,et al.  Simplified multigrid technique for the numerical solution to the steady-state and transient EHL line contacts and the arbitrary entrainment EHL point contacts , 2001 .

[14]  Martin Berzins,et al.  Using adjoint error estimation techniques for elastohydrodynamic lubrication line contact problems , 2011 .

[15]  Martin Berzins,et al.  Calculation of friction in steady-state and transient EHL simulations , 2003 .

[16]  Cornelis H. Venner,et al.  Higher-Order Multilevel Solvers for the EHL Line and Point Contact Problem , 1994 .

[17]  H P Evans,et al.  Evaluation of deflection in semi-infinite bodies by a differential method , 2000 .

[18]  Wassim Habchi,et al.  A full-system finite element approach to elastohydrodynamic lubrication problems , 2008 .

[19]  R W Snidle,et al.  Elastohydrodynamic Lubrication of Heavily Loaded Circular Contacts , 1989 .

[20]  A. Brandt,et al.  Multilevel matrix multiplication and fast solution of integral equations , 1990 .

[21]  A. Lubrecht,et al.  MultiLevel Methods in Lubrication , 2013 .

[22]  C. E. Goodyer,et al.  Adaptive Numerical Methods for Elastohydrodynamic Lubrication , 2001 .

[23]  Wassim Habchi,et al.  Stabilized fully-coupled finite elements for elastohydrodynamic lubrication problems , 2012, Adv. Eng. Softw..

[24]  D. Griffin,et al.  Finite-Element Analysis , 1975 .

[25]  Philippe Vergne,et al.  A finite element approach of thin film lubrication in circular EHD contacts , 2007 .

[26]  Hongqiang Lu,et al.  High order finite element solution of elastohydrodynamic lubrication problems , 2006 .

[27]  Henry Peredur Evans,et al.  On the coupling of the elastohydrodynamic problem , 1998 .

[28]  Leszek Demkowicz,et al.  A fully automatic hp-adaptivity for elliptic PDEs in three dimensions , 2007 .

[29]  Martin Berzins,et al.  High-order discontinuous Galerkin method for elastohydrodynamic lubrication line contact problems , 2005 .

[30]  T. Stolarski.  Tribology in Machine Design , 1990 .

[31]  Ethan J. Kubatko,et al.  hp Discontinuous Galerkin methods for advection dominated problems in shallow water flow , 2006 .

[32]  I. Babuska,et al.  A DiscontinuoushpFinite Element Method for Diffusion Problems , 1998 .

[33]  Martin Berzins,et al.  Adaptive high-order finite element solution of transient elastohydrodynamic lubrication problems , 2006 .

[34]  J. Oden,et al.  A discontinuous hp finite element method for convection—diffusion problems , 1999 .

[35]  Chi-Wang Shu,et al.  Discontinuous Galerkin Methods: Theory, Computation and Applications , 2011 .

[36]  David L. Darmofal,et al.  DEVELOPMENT OF A HIGHER-ORDER SOLVER FOR AERODYNAMIC APPLICATIONS , 2004 .