Contact mechanics of rough surfaces is becoming increasingly important in understanding the real behaviours of machine elements in contact. Due to the complicated physical, chemical and mechanical phenomena occurring at surfaces, especially in boundary lubrication, a multiphysics numerical model is essential to capture the behaviour. Boundary Element Method is a well-known numerical approach to model such a problem because of several advantages. Firstly, it is far faster than Finite Element Method since only the boundaries of the solids are discretised. In addition, there is no problem of remeshing the contacting bodies due to plastic deformation and wear. Conventional Boundary Element models simulate movement of contacting surfaces by shif...