Reynolds-averaged Navier-Stokes equations are numerically solved using a secondorder finite difference method for the analysis of turbulent flows around single and tandem hydrofoils advancing under the free surface. The location of the free surface, not known a priori, is computed from the kinematic free surface condition and the computational grid is conformed at each iteration to the free surface deformation. The eddy viscosity model of Baldwin-Lomax is employed for the turbulence closure. The method is validated through the comparision of the numerical results with the experimental data for a single hydrofoil of a Joukowski foil section. A computational study is also carried out to investigate the effect of the submergence depth and the Froude number on the lift and the drag of the hydrofoil. For tandem hydrofoils, computations are performed for several separation distances between the forward and aft foils to see the interference effect. The result shows clearly how the lift and drag change with the separation distance.