In hydraulic turbines, the tip-leakage vortex is responsible for flow instabilities and for promoting erosion due to cavitation. To better understand the tip vortex flow, Reynolds-averaged Navier–Stokes (RANS) and large eddy simulation (LES) computations are carried out to simulate the flow around a NACA0009 blade including the gap between the tip and the wall. The main focus of the study is to understand the influence of the gap width on the development of the tip vortex, as for instance its trajectory. The RANS computations are performed using the open source solver OpenFOAM 2.1.0, two incidences and five gaps are considered. The LESs are achieved using the YALES2 solver for one incidence and two gaps. The validation of the results is performed by comparisons with experimental data available downstream the trailing edge. The position of the vortex core, the mean velocity and the mean axial vorticity fields are compared at three different downstream locations. The results show that the mean behaviour of the tip vortex is well captured by the RANS and LES computations compared to the experiment. The LES results are also analysed to bring out the influence of the gap width on the development of the tip-leakage vortex. Finally, a law that matches the vortex trajectory from the leading edge to the mid-chord is proposed. Such a law can be helpful to determine, in case of cavitation, if the tip vortex will interact with the walls and cause erosion.

Due to the penetration of alternative renewable energies, the stabilization of the electrical power network relies on the off-design operation of turbines and pump-turbines in hydro-power plants. The occurrence of cavitation is however a common phenomenon at such operating conditions, often leading to critical flow instabilities which undercut the grid stabilizing capacity of the power plant. In order to predict and extend the stable operating range of hydraulic machines, a better understanding of the cavitating flows and mainly of the transition between stable and unstable flow regimes is required. In the case of Francis turbines operating at full load, an axisymmetric cavitation vortex rope develops at the runner outlet. The cavity may enter self-oscillation, with violent periodic pressure pulsations. The flow fluctuations lead to dangerous electrical power swings and mechanical vibrations, dictating an inconvenient and costly restriction of the operating range. The present paper reports an extensive numerical and experimental investigation on a reduced scale model of a Francis turbine at full load. For a given operating point, three pressure levels in the draft tube are considered, two of them featuring a stable flow configuration and one of them displaying a self-excited oscillation of the cavitation vortex rope. The velocity field is measured by two-dimensional (2D) particle image velocimetry (PIV) and systematically compared to the results of a simulation based on a homogeneous unsteady Reynolds-averaged Navier–Stokes (URANS) model. The validation of the numerical approach enables a first comprehensive analysis of the flow transition as well as an attempt to explain the onset mechanism.

Due to the integration of new renewable energies, the electrical grid undergoes instabilities. Hydroelectric power plants are key players for grid control thanks to pumped storage power plants. However, this objective requires extending the operating range of the machines and increasing the number of start-up, stand-by, and shut-down procedures, which reduces the lifespan of the machines. CFD based on standard URANS turbulence modeling is currently able to predict accurately the performances of the hydraulic turbines for operating points close to the Best Efficiency Point (BEP). However, far from the BEP, the standard URANS approach is less efficient to capture the dynamics of 3D flows. The current study focuses on a hydraulic turbine, which has been investigated at the BEP and at the Speed-No-Load (SNL) operating conditions. Several “advanced” URANS models such as the Scale-Adaptive Simulation (SAS) SST $k-\omega$ and the BSL- EARSM have been considered and compared with the SST $k-\omega$ model. The main conclusion of this study is that, at the SNL operating condition, the prediction of the topology and the dynamics of the flow on the suction side of the runner blade channels close to the trailing edge are influenced by the turbulence model.

Hydropower is a key source of electricity production for allowing the integration of intermittent renewable energy resources. Among the various hydraulic power plants around the world, the ones equipped with Pelton turbines already provide large flexibility that is still enhanced with the development, for instance, of the hydraulic short circuit operating mode. However, the knowledge of the flow inside Pelton turbines is still a challenging task, both numerically and experimentally, despite progress in the last two decades. One key feature of the Pelton efficiency is the jet quality, i.e., the jet velocity needs to be uniform, not perturbed by secondary flows and compact. The compactness of the jet is mainly dependent o nthe location of the jet detachment at the nozzle outlet, which is challenging for computational fluid dynamics simulations mainly due to numerical diffusion. Even if this point has already been mentioned in previous papers, the present paper focuses on all the parameters that can influence the jet detachment: the nozzle geometry, the mesh and the numerical scheme used to discretize the convective fluxes. The simulations of an existing Pelton injector are performed using the OpenFOAM toolbox. It is noticed that, in addition to the nozzle geometry and the mesh resolution at the nozzle outlet, the choice of the numerical schemes influences the jet detachment and, consequently, the jet diameter and discharge. The use of an anti-diffusive scheme such as the “SUPERBEE” limiter improves the prediction of the jet in accordance with the on-site measurements.

The flexibility of the FMHL+ pumped storage power plants can be improved by extending the hydraulic short-circuit operating mode. CFD simulations of the flow in three bifurcations are performed to calculate the head losses and to investigate the flow topology in the pipes. A specific attention is paid to the influence of the curvature correction that has been developed for two-equation RANS turbulence models. For the T-junction considered, the activation of the curvature correction influences the head losses whereas for the two Y-junctions computed, no effect is observed. By comparing with the Y-junctions, the T-junction leads to higher head losses and helicity in the pipes downstream of the bifurcation. Compared to the current the intragroup hydraulic short circuit operation permitted, the intergroup and interplant hydraulic short circuit mode should provide better performances with possible gains until of −55% in head losses and −94% in helicity upstream of the turbines.