Fluids slip on superhydrophobic surfaces. The slip velocity is modeled by Navier's slip-length. A user defined function (UDF) in ANSYS Fluent 15.0 was developed to implement the slip boundary condition. The UDF was validated for different values of slip length by two benchmark solutions in laminar and turbulent flows. We utilized a periodic approach to model the VLH turbine in steady-state condition. For modeling, single-phase is assumed and the shear stress transport model was used for turbulence modeling. The results of employing partial slip at different positions of the runner blade and other turbine components with a constant slip-length of 50 mu m showed an approximately 4% improvement in the turbine hydraulic efficiency, at the design point and when only the runner blade is superhydrophobic. In addition to the design point analysis, the turbine efficiency values were studied at part and high load ranges. The entropy generation method was applied to the output results of the simulated cases to determine the energy dissipating zones and to detect the effect of superhydrophobic walls on the dissipation mechanisms. Since this work is the first numerical evaluation of superhydrophobic surfaces for VLH, the conclusions can be meaningful for future numerical and experimental studies. (C) 2019 Elsevier Ltd. All rights reserved.