Numerical evidence for Kosterlitz—Thouless transition in the 2D XY Villain model