This paper describes a fast numerical solver for the elastohydrodynamically lubricated circular contact problem. As a result of the application of multigrid techniques the algorithm is of very low complexity, i.e. O(n ln n), where n is the number of nodes on the grid. Consequently, it enables the solution of the problem using large numbers of nodes (O(105)) on a small-capacity computer, as is demonstrated by the results presented for an example load situation.
Venner, C. H., & ten Napel, W. E. (1992). Multilevel solution of the EHL circular contact problem, part 1: theory and numerical algorithm. Wear, 1(152), 351-367. https://doi.org/10.1016/0043-1648(92)90132-R