Le présent travail s'inscrit dans le cadre d'une collaboration entre le Laboratoire de Mécanique des Fluides et d'Acoustique à l'École Centrale de Lyon (LMFA-ECL), Turbomeca et l'ONERA. Le sujet de recherche porte sur l'interaction rouet-diffuseur et sur l'entrée en régime de pompage d'un étage de compresseur centrifuge transsonique à fort taux de compression étudié à vitesse de rotation de croisière (0,927Nn). L'alimentation des analyses est réalisée par trois méthodes de mesures et deux types de simulations numériques. Le module d'essai est installé sur le banc d'essai 1 MW du LMFA. La caractérisation expérimentale du compresseur TM est réalisée par le biais de mesures de pression et température auxquelles sont adjointes des mesures du débit et de la vitesse de rotation de la roue mobile. La description de l'écoulement interne au compresseur s'appuie sur les résultats de sondages par Anémométrie Laser à effet Doppler (LDA) et de mesures de pression à haute fréquence. Les simulations numériques sont réalisées par l'intermédiaire du code de calcul elsA développé par l'ONERA, qui permet de résoudre le système d'équations de Navier-Stokes couplé à un modèle de turbulence k-l de Smith. Deux modélisations de l'interaction rouet-diffuseur permettent de générer des champs aérodynamiques stationnaires (modèle plan de mélange) et instationnaires (modèle chorochronique). La comparaison entre les données expérimentales et numériques est très satisfaisante et permet alors de profiter pleinement de la richesse des informations numériques. L'examen détaillé de l'écoulement interne au rouet pour trois points de fonctionnement (à débit bloqué, à rendement maximum et à proximité du pompage) révèle que, du blocage vers le pompage, l'évolution de l'intensité et de la taille du tourbillon de jeu est le point de départ d'un enchaînement de mécanismes conduisant à la dilatation du sillage de la structure jet-sillage. Pour l'écoulement en amont du diffuseur, ceci ce traduit en moyenne temporelle par une augmentation de l'incidence principalement au voisinage du moyeu. L'onde de choc en amont des aubes du diffuseur remonte à mesure que le débit du compresseur diminue. La trajectoire de l'écoulement principal bascule du côté de la face en dépression vers le côté de la face en pression du canal inter-aubes du diffuseur. La cartographie des nombreux décollements de couches limites est également modifiée à l'approche du pompage. Les structures instationnaires majeures sont produites par l'interaction de l'onde de choc en amont des aubes du diffuseur avec les pales du rouet. Des ondes de pression progressives et des poches à faible nombre de Mach sont ainsi générées. Les ondes pression impriment d'intenses fluctuations au champ de vitesse qui favorisent le processus de mélange. En conséquence, les couches limites sont plus robustes vis-à-vis des décollements (en moyenne temporelle). Le défilement instationnaire des poches à faible nombre de Mach engendre une dissymétrie marquée des conditions d'alimentation du diffuseur dans la direction azimutale. Au cours du changement de point de fonctionnement en allant du débit bloqué vers le pompage, les ondes de pression se renforcent et la taille des poches à faible nombre de Mach diminue. De ce fait, les conditions sont plutôt favorables à retarder l'entrée en pompage du compresseur qui est localement initié en amont du diffuseur aubé. Le pompage du compresseur est provoqué par un changement rapide de la structure supersonique de l'écoulement en entrée de diffuseur qui est alors déstabilisée par les fluctuations de pression des ondes progressives. / The present work is in line with a collaboration between the Laboratoire de Mécanique des Fluides et d'Acoustique at École Centrale de Lyon (LMFA-ECL), Turbomeca and ONERA. The subject is focused on impeller-diffuser interaction and phenomena occurring during surge ignition of a transonic centrifugal stage with high-compression ratio at cruise rotation speed (0,927Nn). The analysed data come from three measurement devices and two kinds of numerical simulations. The 1MW LMFA-ECL test rig was used for carrying out the experiments on the centrifugal compressor stage. The global performances were obtained by pressure, temperature, mass flow rate and rotation speed measurements. The internal flow field properties were probed by Laser Doppler Anemometry (LDA) and high frequency pressure measurements. The computations were performed with the elsA software developed at ONERA. The code solves the compressible Reynolds Averaged Navier-Stokes equations associated with the two-equations (k-l) turbulence model of Smith. Two models of impeller-diffuser interaction were used to simulate the flow within the compressor. The first model is based on the Averaged Passage equations, gives a stationary description of the internal flow fields. The second model is based on the phase lagged approach and allows access to the unsteady phenomena. A good agreement between experiments and simulations was obtained, which justifies the use of the CFD results for the flow field analysis. Analysis of the flow development within the impeller were carried out for three operating points (choke, peak efficiency and close to surge). This study shows that, from choke to near surge, the development of the tip clearance vortex is the starting point of a sequence of physical mechanisms which lead to an extend of the wake of the jet-wake structure. The change in flow field at impeller exit tends to increase the upstream incidence of the vaned diffuser, especially close to the hub. The vane bow shock wave progresses in the impeller passages when the operating point moves from choke to near surge. The location of the main flow moves from suction side to pressure side of the vaned diffuser passage. The topology of the boundary layer separations within the diffuser passage is deeply affected when approaching surge. The main unsteady structures are generated by the interaction between the bow shock wave and the impeller blades. Progressive pressure waves and low Mach number flow bubbles are generated by this locally time-spaced interaction phenomenon. The pressure waves lead to strong fluctuations of the velocity field. As a consequence, the boundary layer becomes more resistant in relation to the separations (in term of time-averaged point of view). In time-averaged point of view, the low Mach number flow bubbles lead to inhomogeneous conditions at the vaned diffuser inlet along the azimuthal direction . From choke to near surge operating point, the strength of the pressure wave increases whereas the size of the low Mach number flow bubbles decreases. These conditions are quite favourable to push back the surge limit which is locally initiated in the inlet part of the vaned diffuser. The surge of the compressor is produced by a rapid change in supersonic flow structure at diffuser inlet. The pressure fluctuations due to the progressive waves lead to destabilise the new supersonic flow configuration and degenerate in the global instability of the compressor flow.

En configuration croisière d'un avion de transport civil, où les conditions de l'écoulement sont transsoniques, une augmentation de la vitesse ou de l'incidence, nécessaire par exemple lors d'un virage ou d'un changement d'altitude, peut générer des décollements sur l'extrados de l'aile par la recompression de l'écoulement ou par l'interaction onde de choc/couche limite plus forte. Ces décollements modifient la répartition de charge sur la voilure, donc la portance, et deviennent instables par les échappées tourbillonnaires dans le sillage lorsqu'ils sont étendus jusqu'au bord de fuite. Ces instabilités aérodynamiques, correspondent à la répartition spatiale et temporelle des décollements et de la position du choc. Elles modifient la charge de la voilure et par couplage excitent la structure sur ces modes propres, tremblement. Par conséquent, il apparaît intéressant d'approfondir l'étude des instabilités qu'on vient d'identifier à travers d'une nouvelle approche. De ce fait, ce mémoire s'orientera vers une analyse de stabilité linéaire globale d'un écoulement compressible turbulent soumis au phénomène de tremblement. En particulier, la démarche de cette thèse s'inspirera des travaux effectués par l'équipe de Crouch. Dans ce contexte, l'objectif de cette thèse est le développement d'un outil numérique permettant d'étudier la stabilité globale d'un écoulement compressible turbulent pour des géométries quelconques. Ce travail est essentiellement numérique, il mettra en oeuvre l'implémentation du code dans le logiciel elsA de l'ONERA, pour calculer le champ de base et réaliser l'extraction explicite de la jacobienne des équations RANS linéarisées. Cet objectif, a pour but de prouver qu'on est capable de retrouver les acteurs responsables du phénomène de tremblement qui apparaît à l'extrados d'une voilure bidimensionnelle. En vue de l'étude de la stabilité de l'écoulement, l'implémentation des méthodes numériques ainsi que leurs validations au sein du logiciel elsA, ont été réalisé avec succès. Une fois ce résultat est acquis, on s'intéressera, en particulier, à déterminer le spectre des modes propres, grâce aux techniques spécifiques fournies par la théorie de stabilité.

Les solutions exactes autosemblables de l'hydrodynamique, avec conduction thermique non-linéaire, pour des gaz parfaits en milieu semi-infini, constituent une approche possible pour l'étude de la stabilité d'écoulements en fusion par confinement inertiel. Ces solutions et leurs perturbations linéaires, calculées 'a l'aide d'une méthode spectrale de Chebyshev multidomaine, permettent de rendre compte, sans autre approximation, du caractère compressible et instationnaire de ces écoulements. Suite aux résultats (Clarisse et al., 2008; Lombard, 2008) représentatifs de l'ablation précoce d'une cible par un flux laser non-uniforme (conduction électronique, front thermique subsonique précédé d'un choc quasi-parfait), nous explorons ici d'autres configurations. Pour cette même phase d'ablation, mais avec un rayonnement X incident non-uniforme (conduction radiative), nous traitons un écoulement compressible et un autre faiblement compressible. Dans les deux cas, nous retrouvons les comportements des écoulements compressibles obtenus en conduction électronique avec une instabilité maximale pour un nombre d'onde nul. D'autre part, la méthode spectrale est étendue au calcul de solutions autosemblables en tenant compte de l'onde thermique supersonique en amont du choc. Basée sur une analyse des singularités des équations réduites (front infiniment raide), cette méthode permet d'accéder au régime d'onde thermique supersonique propre à l'irradiation initiale d'une cible et de retrouver les solutions ablatives calculées antérieurement dans l'approximation de précurseur thermique négligeable.