In this paper, the hydrodynamic characteristics and flow field around rectangular and delta hydrofoils, moving with a constant speed beneath the free surface are numerically studied by means of isoparametric boundary element method (IBEM). The quantities (source and dipole strengths) and the geometry of the elements are represented by a linear distribution. Two types of three-dimensional hydrofoils (rectangular and delta) are selected with NACA4412 and symmetric Joukowski sections. Some numerical results of pressure distribution, lift, wave-making drag coefficients and velocity field around the hydrofoils are presented. Also, the wave pattern due to moving hydrofoil is predicted at different operational conditions. Comparisons are made between computational results obtained through this method and those from the experimental measurements and other numerical results which reveal good agreement.