Accuracy and Grid Convergence of the Numerical Solution of the Energy Equation in Fluid Film Lubrication: Application to the 1D Slider

Abstract

The present work is focused on the numerical solution of the complete energy equation used in fluid film lubrication. The work was motivated by the fact that the complete energy equation has no analytical solution that can be used for validations. Its accuracy and computation time are related to the employed numerical method and to the grid resolution. The natural discretization method (NDM) applied on different grids is systematically compared with the spectral method (the Lobatto Point Colocation Method or LPCM) with different polynomial degrees. A one dimensional inclined slider is used for the numerical tests, and the energy equation is artificially decoupled from the Reynolds equation. This approach enables us to focus all the attention on the numerical solution of the energy equation. The results show that the LPCM is one or two orders of magnitude more efficient than the NDM in terms of computation time. The energy equation is then coupled with the Reynolds equation in a thermo-hydrodynamic analysis of the same 1D slider; the numerical results confirm again the efficiency of the LPCM. A thermo-hydrodynamic analysis of a two-lobe journal bearing is then presented as a practical application.

Keywords

• energy equation
• numerical solution
• lobatto point quadrature
• legendre polynomials