Numerical simulations have been carried out for a 138K LNG Carrier (KLNG) model ship with free surface, using WAVIS 2.0 and Fluent 6.3.26 with various $y^+$ values and different grid densities. Level-set method for free surface capturing was adopted in WAVIS, while VOF has been used in Fluent. The calculated results were compared with the experiment data. Resistance coefficient, wave pattern, wave profile along the hull surface, axial velocity contours and transverse vectors had been analyzed. When the first $y^+$ value was fixed at 60, the simulation results from both WAVIS and Fluent were improved as the number of grids increased. The convergence time of WAVIS was much shorter than that of Fluent. Furthermore, WAVIS predicted the velocity field and the wave profile along the hull surface better than Fluent. However, Fluent gave better wave patterns.