Simulation aux Grandes Échelles pour la modélisation aérothermique des aubages de turbines refroidies

This PhD dissertation, conducted as part of a CIFRE research project between TURBOMECA and CERFACS in partnership with the VKI, deals with improving performance of axial turbines from helicopter engines. One of the most critical design points of such engines is the control of the high pressure turbine blade lifetime which face the high temperatures from the combustor. Today, industrial numerical aerothermal predictions of the flows around the blade (in the vein and in its cooling system) are performed with the Reynolds Averaged Navier-Stokes (RANS). Thanks to the increasing computational power, Large Eddy Simulation (LES) becomes affordable to offer further flow predictions. Therefore, this thesis focuses on the capabilities of the LES to estimate the flow in turbine blade internal cooling channels. To simplify this analysis where several physical phenomenon are present, the problem is described in three parts with increasing complexity. The first part addresses simplified typical geometries of cooling channel (U-bend and ribbed channel) in a static configuration. Considering the flow regime, a wall-resolved approach using a hybrid unstructured mesh is proposed in view of the application on an industrial case. The second part extends the study of the ribbed channel in rotation using an inertial reference frame. LES provides mean and unsteady results in good agreement with the available experimental data and previous works, for the flow dynamic and the heat transfer. Finally, the third part presents the application of the method to an industrial case with conjugate heat transfer between a complex cooling channel and the blade. This last section is not present in the public manuscrit for confidential reasons. Results of the use of the wall-resolved approach in rotation in an inertial frame of reference are compared to RANS predictions and show the potential of the method with high local differences.Ce travail de thèse, réalisé dans le cadre d'une convention CIFRE entre TURBOMECA et le CERFACS et en partenariat avec l'IVK, se place dans un contexte d'amélioration des performances des turbines axiales équipant les turboréacteurs d'hélicoptère. Un des points critiques du dimensionnement de tels moteurs est la maitrise de la durée de vie des pales de la turbine haute pression qui font face à de très hautes températures provenant de la chambre de combustion. Les prédictions numériques de l'environnement aérothermique des pales (écoulements dans la veine et système de refroidissement) sont réalisées aujourd'hui dans le milieu industriel à l'aide de la modélisation Reynolds Averaged Navier-Stokes (RANS). Grâce à des capacités de calculs grandissantes, l'approche Simulation aux Grandes Echelles (SGE) offre désormais un nouveau potentiel de prédictions d'écoulements. Les travaux de cette thèse s'intéressent ainsi à la capacité de la SGE à prédire l'écoulement du circuit de refroidissement interne d'une pale de turbine. Pour simplifier l'analyse de ce problème ou plusieurs phénomènes physiques sont en jeu, une progression en trois parties est proposée. La première s'intéresse à l'étude aérothermique de géométries simplifiées de canaux de refroidissement (coude à 180° et canal avec promoteurs de turbulence) en configuration statique. Aux régimes d'écoulement considérés, une approche résolue en paroi avec maillage non-structuré hybride est proposée et validée en vue d'une application industrielle facilitée. La seconde partie étend l'analyse de l'écoulement à un cas de canal avec promoteurs de turbulence en rotation utilisant une méthode de résolution numérique dans un repère absolu. Les investigations des résultats de la SGE fournissent des prédictions moyennes et instationnaires en bon accord avec les expériences disponibles et les travaux précédents aussi bien pour la dynamique de l'écoulement que les transferts de chaleur. Enfin, une troisième partie présente une application de la méthode sur un cas de pale réelle avec couplage thermique entre le circuit de refroidissement et le solide de la pale. Cette dernière partie classée confidentielle n'est pas présente dans le manuscrit disponible publiquement. Les résultats de l'approche résolue en paroi et de la rotation dans le repère absolu comparés aux résultats RANS disponibles pour le cas applicatif montrent d'importante différences locales et ainsi le potentiel de la méthode proposée