Ce projet présente l’étude expérimentale de l’écoulement à la sortie du diffuseur d’une turbine de type bulbe. La technique de mesure d’anémométrie laser à effet Doppler (LDV) a été utilisée pour les trois montages expérimentaux, situés au-dessus et sur les deux côtés du diffuseur pour obtenir les trois composantes de vitesse. Une importante chute de rendement a été mesurée dans un modèle réduit d’une turbine de type bulbe opérant à forte charge. Des études précédentes ont relié la chute de performance avec les pertes du diffuseur, et en particulier avec la zone de séparation de l'écoulement aux parois du diffuseur. Dans la présente étude, l'écoulement a été étudié à la sortie de la trompette de la turbine, qui est une section du diffuseur qui permet de passer d'une section circulaire à une section rectangulaire. La turbine a été étudiée pour cinq conditions d'opération, qui représentent les différents phénomènes de l'écoulement à la sortie du diffuseur. En plus du champ de vitesse, l'analyse a été effectuée pour le coefficient d’intermittence de la vitesse débitante et pour la composante de la vorticité autour l’axe axial Z. Les résultats révèlent une zone contrarotative dans le diffuseur, qui s'intensifie avec l'ouverture des directrices. L'ouverture des directrices induit une modification aux phénomènes hydrauliques : à partir d'une recirculation de l'écoulement dans la zone centrale pour la condition d'opération à plus faible charge, vers un écoulement de retour induit par la séparation de l'écoulement aux parois pour les conditions d'opérations à forte charge. / This project presents the experimental study of the flow at the diffuser outlet of a bulb turbine. Measurements by laser Doppler velocimetry (LDV) were performed on three experimentals setups located above and on two sides of the diffuser to obtain the three velocity components. An important drop in turbine performances has been measured in a bulb turbine model operated at overload. Previous investigations have correlated the performance drop with diffuser losses, and particularly to the flow separation zone at the diffuser wall. In the present study, the flow has been investigated at the exit of the turbine, which is a diffuser section that transforms from a circular to a rectangular section. The turbine has been operated at five operating points, which are representative of different flow patterns at the diffuser exit. In addition to the average velocity field, the analysis has been conducted based on a backflow occurrence function and on the vorticity level. Results reveal a counter-rotating zone in the diffuser, which intensifies with the guide vanes opening. The guide vanes opening induces a modification of the flow phenomena: from a central backflow recirculation zone at the lowest flowrate to a backflow zone induced by flow separation at the wall at the highest flowrate.

Tableau d'honneur de la Faculté des études supérieures et postdorales, 2015-2016 / Les travaux réalisés dans le cadre de ce mémoire de maîtrise visent à décrire la source hydraulique de vibrations mesurées lors de la mise en service d’une roue Francis. Lors de l’essai de rejet de charge à partir de la puissance maximale, il a été relevé que la roue subissait des vibrations autoexcitées à l’une de ses fréquences naturelles. L’objectif de ces travaux est d’identifier le phénomène hydraulique causant l’entrée en vibrations auto-excitées de la roue submergée. Puisque la cavitation semble jouer un rôle dans cette problématique, les capacités du solveur fluide ANSYS CFX à simuler de la cavitation instationnaire ont été évaluées. De ce fait, deux géométries ayant été étudiées expérimentalement ont été simulées. Cette étude a permis de conclure que l’approche homogène, plus stable et couramment utilisée dans la littérature, ne permet pas de reproduire de façon fiable les fluctuations de pression causées par l’instationarité des cavités de vapeur. Il a cependant été montré que les simulations permettent de prédire la forme, la localisation ainsi que les mécanismes entraînant la présence de vapeur dans l’écoulement. De plus, il a été demontré que la cavitation est un phénomène particulièrement sensible et sujet à répondre en phase à une excitation oscillatoire, par exemple la vibration de la roue. Les simulations numériques réalisées durant différentes phases du rejet de charge transitoire ont entre autres permis d’identifier que lorsque les vibrations apparaissent, une forte région de vapeur se crée au bord de fuite de l’aubage, près du plafond. En augmentant le temps de fermeture du distributeur, le partenaire industriel a réussis à éliminer les vibrations problématiques. En réalisant des simulations avec différents temps de fermeture, il a été démontré que la solution proposée permet d’augmenter le niveau de pression dans le canal inter-aubes, réduisant ainsi la quantité de vapeur s’y trouvant. Cela laisse suggérer que la cause d’entrée en vibrations de la roue est la cavitation se formant durant le rejet de charge. Cependant, il est à noter que les différentes méthodologies proposées n’ont pas permis d’obtenir les fréquences d’excitation mesurées expérimentalement, essentiellement à cause de limitations liées à la modélisation de la cavitation. / The work presented in this master degree thesis aims to identify the hydraulic cause of mechanical vibrations measured during the commissioning of a Francis runner. During the test of load rejection from maximal power output regime, it was noticed that the runner entered a state of self-excited vibrations at one of its natural frequencies. The purposes of this work is to investigate the hydraulic phenomenon which causes the submerged runner to enter self-excited vibrations. Since cavitation is expected to play a role in this problematics, there was a need to evaluate the capabilities of the fluid solver ANSYS CFX to solve unsteady cavitating flows. Two geometries which had been investigated experimentally were thus simulated. It was concluded that the homogeneous approach, more robust and widely used in the literature, does not lead to a reliable prediction of the pressure fluctuations caused by cavitation. It was however shown that the simulations allowed to predict the shape, the location as well as the physical mechanisms responsible for the presence of vapor in the flow. It was also demonstrated that cavitation is a phenomenon particularly sensitive and subject to respond in phase to oscillatory perturbations, for instance the vibrating runner. The numerical simulations carried out at different phases of the load rejection transient have established that when the vibrations appear, a wide region of vapor forms at the trailing edge of the blade, near the crown. By increasing the distributor closing time, the industrial partner in this work had success in eliminating such problematic vibrations. In our case, when performing simulations with increased closing times, it was demonstrated that the solution proposed allows to increase the pressure level in the inter-blade channel, lowering the quantity of vapor it contains. This strongly suggests that the hydraulic cause of the vibrations is the cavitation forming during the load rejection. However, one can note that the proposed methodologies have not permitted to predict the excitation frequencies as measured experimentally, essentially due to limitations in the modeling of cavitation.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2017-2018 / Ce mémoire vise à contribuer à l'amélioration de la prédiction numérique des performances d'une turbine hydraulique par la quantification de l'impact de certains paramètres d'influence, autant physiques que méthodologiques. Pour ce faire, l'écoulement dans une turbine caractérisée par une chute de rendement près du meilleur point de fonctionnement a été étudié. Cette dernière est associée à une baisse abrupte du coefficient de récupération de l'aspirateur, une composante connue pour jouer un rôle capital dans les performances globales de la machine. L'accent a particulièrement été mis sur l'étude du débalancement de l'écoulement sous la roue, observé à la suite d'une campagne de mesure expérimentale effectuée sur le banc d'essai du CREMHyG à Grenoble. En effet, l'écoulement en sortie de roue est souvent considéré axisymétrique dans les simulations numériques, une simplification qui est en désaccord avec les données expérimentales du cas étudié dans ce projet de recherche. Afin d'identifier la nature du débalancement, une analyse critique a été menée sur les données expérimentales. Toutefois, aucune conclusion satisfaisante n'a pu être tirée puisqu'une inconsistance entre les données obtenues par LDV et celles obtenues par PIV a été décelée. Considéré statique, il a été démontré que la présence d'un tel débalancement à l'entrée de l'aspirateur altère considérablement la topologie de l'écoulement en aval ainsi que la prédiction des performances. De plus, certains aspects des simulations obtenus avec ce type de débalancement ne concordent pas avec les observations expérimentales, c'est notamment le cas pour la zone d'amorçage du décollement mais aussi pour l'écart de débit entre les deux pertuis. Il a également été montré que les simulations des composantes fixes en amont de l'aspirateur, incluant une importante partie de la conduite d'amenée, ne permettent pas de retrouver l'asymétrie de l'écoulement observée expérimentalement en sortie de roue. En effet, l'écoulement secondaire, associé à la présence d'une conduite coudée en amont de la bâche spirale, s'est avéré avoir très peu d'influence sur la distribution de l'écoulement dans l'avant-distributeur et dans le distributeur. La conduite menant jusqu'à l'entrée de la bâche spirale semble être assez longue pour que les modifications apportées au profil de vitesse deviennent insuffisantes pour perturber l'écoulement à cet endroit. / This master degree thesis aims at enhancing the numerical prediction of the performance of a hydraulic turbine by quantifying the impact of physical and methodological influencing parameters. To this end, the flow in a turbine characterised by an efficiency drop near the best efficiency point was studied. The latter is associated with a sharp drop in the recovery coefficient of the draft tube, a component known to play a key role in the global performance of the machine. Particular emphasis was placed on the study of the imbalance of the flow at the runner exit, observed following an experimental measurement campaign carried out on the CREMHyG test bench in Grenoble. Indeed, the flow at the runner exit is often considered axisymetric in numerical simulations, a simplification that is not consistent with the experimental data of the case studied in this research project. In order to identify the nature of the imbalance, a critical analysis was conducted on the experimental data. However, no satisfactory conclusion could be drawn since inconsistency between the data obtained by LDV and those obtained by PIV was identified. Considered static, it has been shown that the presence of such an imbalance at the inlet of the draft tube alters considerably the topology of the downstream flow as well as the prediction of the performances. Moreover, some aspects of the simulations obtained with this type of imbalance do not agree with the experimental observations, it is notably the case for the initiation zone of the flow detachment but also for the flow rate imbalance between the two draft tube channels. It has also been shown that the numerical simulations of fixed components upstream of the draft tube, including a large part of the penstock, do not allow to recover the asymmetry of the flow experimentally observed at the exit of the runner. Indeed, the secondary flow, associated with the presence of a bend upstream of the spiral case, proved to have a very little influence on the distribution of the flow in the distributor. The pipe leading to the entrance of the spiral case appears to be long enough for changes made in the velocity profile to become insufficient to disrupt the flow at this place.

Ce mémoire a comme objectif de prédire les performances d’une turbine modèle de type bulbe à l’aide de simulations numériques RANS k-ε. Tout d’abord, par des simulations numériques, on valide les éléments de la méthodologie les plus susceptibles d’influencer la qualité des résultats. Par la suite, on obtient la colline de rendement numérique et on analyse l’effet du jeu de bout d’aube ainsi que des jeux au moyeu sur la prédiction de performance. Les résultats indiquent que la colline de rendement numérique se compare bien à la littérature portant sur les turbines bulbes. On montre, de plus, que tenir compte du jeu de bout d’aube dans la simulation numérique diminue de manière significative le rendement prédit. L’impact des jeux au moyeu est cependant négligeable sur la prédiction de performance de la turbine. Finalement, une modification de l’aspirateur est proposée pour la suite du projet entrepris par le consortium de recherche.

L'objectif du présent mémoire est de clarifier l'influence qu'ont sur la solution divers paramètres de calcul couramment utilisés en mécanique des fluides numérique, dans le but d'améliorer la prédiction de l'écoulement dans les aspirateurs de turbines hydrauliques. Pour y arriver, la résolution eulérienne des équations de Navier-Stokes est faite par une approche RAXS à l'aide du logiciel commercial ANSYS CFX. Pour les trois cas tests sélectionnés, la turbulence est modélisée à l'aide des modèles à deux équations k ? e ou SST. Les résultats principaux font ressortir très clairement l'impact de la composante radiale de la vitesse ainsi que des quantités turbulentes imposées sur la frontière amont du domaine de calcul. On remarque aussi que l'extension ajoutée en aval du domaine d'intérêt peut quant à elle faire varier de façon significative la pression statique dans la dernière portion de la géométrie et ainsi influencer le coefficient de récupération mesuré. Finalement les deux modèles de turbulence utilisés sont eux aussi susceptibles d'influencer le développement de l'écoulement.