A numerical simulation code was developed for estimating the surface temperature of POM wear test with different shapes. By using the developed code, the effect of sliding geometry with various contact areas on the surface temperature was numerically investigated. The simulation results show that the limiting sliding velocity depends on contact area and reaches a peak value at a certain geometry. It was also found that the simulation code leads an exact prediction of limiting velocity in the range of small contact area. This is important since the exact temperature measurement is very difficult in this range.