The turbulent free surface flow around a self-propelled KRISO 138K LNG Carrier is numerically simulated using the finite volume based multi-block RANS code, WAVIS developed at HRISO. The realizable k-$varepsilon$ turbulence model with a wail function is employed for the turbulence closure. The free surface is captured with the Level-Set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. In order to obtain an accurate free surface solution and stable convergence, the computations are executed with a proper fine grid refinement around the free surface and with an adoption of implicit discretization scheme for the Level-Set formulation. The computed velocity vectors at the several stations and wave patterns show a good agreement with the experimental results measured at the KRISO towing tank.