Knowing the characteristics of the flow around the ship hull is a very important element in the process of ship design. Nowadays, the most widely used method for the determination of the hydrodynamic characteristics is to conduct experiments with ship models in the towing tank, which is a very expensive and time consuming method. In order to determine the hydrodynamic characteristics in a much shorter period of time and to reduce the costs of ship design it is possible to use computational fluid dynamics (CFD). With this approach it is possible to gain insight into the details of the flow around the ship hull at an early design stage and to get a recommendation on how to improve a specific project or to choose the most suitable project for model testing. In this thesis, numerical simulation of viscous flow around the ship hull was carried out by use of software package STAR-CCM+. A mathematical model based on Reynolds averaged Navier – Stokes equations, k - ε turbulence model and Volume of Fluid method (VOF) for describing the motion of two-phase media are given. Also, necessary boundary conditions for the mathematical model and the method of implementation of discretization are described. The influence of the grid density on the results for the total ship resistance was investigated using three different grid densities. Results of total ship resistance, obtained by numerical simulations, were compared to available experimental results. It has been shown that for all three grid densities satisfactory agreement with experimental results can be achieved.