This study presents a numerical investigation of NACA 0012 hydrofoil–free surface interactions using ANSYS Fluent with the k–ε turbulence model. The numerical approach is first assessed by comparing the free surface elevation with experimental and numerical data from the literature. Further, to ensure the accuracy of the numerical approach, a grid refinement study is investigated. Furthermore, the drag and lift coefficients are compared with published results. Next, the analysis of volume fraction contours for various angles of attack (α = 5°, 10°, 20°, and 25°) reveals the strong influence of incidence on free surface deformation, ranging from mild disturbances at low angles to pronounced wave breaking and splashing at higher angles. The temporal evolution of hydrodynamic forces coefficients shows that stationary values of lift and drag can be obtained at small incidence, whereas persistent oscillations dominate at larger angles due to vortex shedding and enhanced hydrofoil-free surface coupling. Overall, this work emphasizes the importance of RANS (Reynolds-Averaged Navier-Stokes) models in capturing detailed flow physics and free-surface dynamics. Nevertheless, further investigations are required to accurately estimate the associated hydrodynamic forces in order to provide insights for the design and optimization of hydrofoil systems under varying operating conditions.