Analyse et contrôle des écoulements en compresseur centrifuge avec diffuseur aspiré

L'étude effectuée au cours de cette thèse s'inscrit dans le contexte du contrôle des écoulements en turbomachines. Suite à l'intérêt avéré de l'aspiration de la couche limite en compresseurs axiaux pour en améliorer les performances, elle vise à évaluer si cette méthodologie de contrôle est pertinente en géométrie radiale. Elle s'appuie pour cela sur un compresseur centrifuge transsonique Turbomeca. Le fonctionnement du compresseur est d'abord analysé sur la base de simulations numériques stationnaires et instationnaires chorochroniques, réalisées grâce au code de calcul elsA développé par l'Onera et en utilisant le modèle de turbulence k-l de Smith. La validité des résultats numériques est vérifiée par comparaison avec les données expérimentales disponibles : performances mono-dimensionnelles et signaux de pressions instationnaires. Ces données sont issues à la fois de la caractérisation expérimentale du compresseur effectuée par Turbomeca et de campagnes de mesures menées a l'ISAE, ou le compresseur est monté sur un banc de recherche dédié à l'étude des interactions entre le rouet et le diffuseur radial. L'analyse détaillée des résultats numériques révèle le développement d'un décollement de coin dans le diffuseur, entre la paroi moyeu et la face en dépression des aubages, lorsque le point de fonctionnement se rapproche du pompage. L'étude topologique du spectre du frottement pariétal permet d'en préciser le lieu et la structure, et de localiser les points singuliers ainsi que la ligne de séparation principale. En stationnaire, la croissance du décollement aboutit au décrochage du diffuseur, et compromet les performances de l'étage de compression et finalement la stabilité numérique. L'exploitation des résultats instationnaires permet de préciser le comportement temporel de l'écoulement : l'étendue du décollement fluctue sous l'effet du défilement des ondes de pression générées par l'interaction rouet-diffuseur, mais les trajectoires des particules fluides impliquées dans le décollement sont en accord avec les lignes de courant du champ de l'écoulement moyenné temporellement. Le décollement instationnaire est fixe, et sa topologie correspond à celle prévue par les simulations stationnaires. En particulier, le lieu du décollement est prévu de façon similaire par les modèles stationnaires et instationnaires. Une stratégie d'aspiration est ensuite développée sur la base de ces observations : une fente de prélèvement est positionnée au voisinage du col de séparation principal, dont le lieu correspond à celui du maximum du gradient de pression adverse stationnaire. Cette stratégie est implémentée au sein des modèles numériques stationnaires et instationnaires. En stationnaire, un prélèvement de 1 % de la valeur totale du débit traversant le rouet permet un contrôle total du décollement dans le diffuseur, et conduit à une augmentation significative de la plage de stabilité numérique. La réduction du débit de prélèvement à 0,3 % du débit total permet un contrôle partiel du décollement, et conduit également à une augmentation de la plage de stabilité numérique. En instationnaire, le décollement de coin initial est contrôlé. Mais les résultats mettent en évidence le rôle majeur joue par le défilement des ondes de pression le long des aubages du diffuseur. Celles-ci se renforcent au passage du col du diffuseur, et engendrent l'existence d'un important gradient de pression adverse instantané. Ce maximum du gradient de pression instantané conduit au développement d'un nouveau décollement de coin, en aval de la fente de prélèvement, et la plage de stabilité des calculs instationnaires n'est pas augmentée. Ces résultats mettent en évidence la possibilité d'agir sur les décollements se produisant dans les diffuseurs radiaux à l'aide de la technique d'aspiration de la couche limite. […]The study presented in this manuscript takes place in the context of flow control within tur - bomachinery. Following the demonstration of the effectiveness of the boundary layer suction technique in axial compressors for improving their performance, the present work aims at determining if this same technique may be relevant in radial compressors. Boundary layer suction is then applied on a centrifugal compressor stage, designed and built by Turbomeca, Safran group. The working of the compressor is first analyzed thanks to steady-state and unsteady numerical simulations, performed using the elsA solver developed by Onera, the French Aerospace Laboratory. The turbulence is modeled with the two equations k-l model of Smith. The validity of the numerical results is ensured by comparison with available experimental measurements results: one-dimensional performance coefficients and time-dependent pressure signals. This data were obtained by both Turbomeca, during the characterization of the compressor, and the fluid mechanics laboratory of ISAE, Université de Toulouse, were the compressor is mounted in an experimental test rig dedicated for studying the impeller-diffuser interactions. The detailed analysis of the numerical results reveals the growth of a corner separation within the diffuser between the hub endwall and the vane suction side when the operating point moves toward surge. The precise location and the internal structure of that boundary layer separation are then elucidated thanks to a topological study, which allows to identify the singular points and the separating lines of the skin-friction pattern. In steady-state numerical simulations, the development of that corner separation leads to the stall of the diffuser, which compromises the compressor stage performance and finally the stability of the numerical model. Unsteady numerical simulations results allow to specify the temporal behavior of the corner separation: the extent of the separated zone is modulated by the scrolling of pressure waves created by the impeller-diffuser interaction, but the trajectories of fluid particles within the separation match with the streamlines of the time-averaged flow field. The unsteady separation is then fixed, and its topology is in agreement with the prediction of the steady-state numerical simulations. In particular, the location of the separation is predicted similarly by the steady-state and the unsteady numerical models. Afterward, a control strategy using suction technique is developed thanks to the previous conclusions. The suction slot is set in the neighborhood of the main saddle of the separation, which corresponds to the location of the maximum of the time-averaged adverse pressure gradient. This strategy is implemented within both the steady-state and unsteady numerical models. The steady-state numerical model predicts the complete control of the separated zone with a removal of 1% of the total massflow through the compressor. This leads to a significant increase of the numerical stable range. With a removal of 0.3% of the compressor total massflow, the numerical model predicts only a partial control of the corner separation, but it also leads to a significant increase of the numerical stable range. With the unsteady numerical model, the initial hub corner separation is again controlled thanks to aspiration. But the results also highlight the major role played by the scrolling of pressure waves. They reinforce when crossing the diffuser throat, and generates a strong instantaneous adverse pressure gradient. This maximum provokes a new boundary layer separation, further downstream of the suction slot. The stable range of the unsteady numerical model is not increased. These results put into evidence the possibility to act on boundary layer separation that occur in radial diffusers thanks to the boundary layer suction technique. [...