A three dimensional thermohydrodynamic lubrication model which couples the Reynolds and energy equations is developed. The model uses the streamline upwind Petrov-Galerkin (SUPG) method to solve the nonsymmetric stiffness ...