Surrogate-based optimization of hydrofoil shapes using RANS simulations