Abstract: An unequal distance multigrid method is developed to solve the thermal elastohydrodynamic lubrication of line contact problem. The results show that the pressure, film thickness and temperature distributions can be more accurately and efficiently calculated by using the unequal distance multigrid method without very dense grid, though the stabilization of the solving process is not improved. Moreover, the pressure spike calculated with the new numerical method is smooth under moderate load.