Two-dimensional (2D) XY model subject to three different types of dynamics, namely Monte Carlo, resistivity shunted junction (RSJ), and relaxational dynamics, is numerically simulated. From the comparisons of the current-voltage (I-V ) characteristics, it is found that up to some constants I-V curves at a given temperature are identical to each other in a broad range of external currents. Simulations of the Villain model and the modified 2D XY model allowing stronger thermal vortex fluctuations are also performed with RSJ type of dynamics. The finite-size scaling suggested in Medvedyeva et al. [Phys. Rev. B (in press)] is confirmed for all dynamic models used, implying that this finite-size scaling behaviors in the vicinity of the Kosterlitz-Thouless transition are quite robust.