In this paper, the wave run-up height and depression depth around air-water interface-piercing circular cylinder have been numerically studied. The Reynolds Averaged Navier-Stokes equations (RANS) and continuity equations are solved with Reynolds Stress model (RSM) and volume of fluid (VOF) method as turbulence model and free surface modeling, respectively. A commercial Computational Fluid Dynamics (CFD) software “Star-CCM+” has been used for the current simulations. Various Froude numbers ranged from 0.2 to 1.6 are used to investigate the change of air-water interface structures around the cylinder and experimental data and theoretical values by Bernoulli are compared. The present results showed a good agreement with other studies. Kelvin waves behind the cylinder were generated and its wave lengths are longer as Froude numbers increase and they have good agreement with theoretical values. And its angles are smaller with the increase of Froude numbers.