Les écoulements vrillés, qui sont largement utilisés dans les turbines à gaz, sont connus pour être sujet à des bifurcations entre différentes topologies (grandes reconfigurations de l'écoulement) qui peuvent affecter les performances et la sécurité du moteur. Ce travail se concentre sur l'étude de ces bifurcations en utilisant la Simulation aux Grandes Echelles (SGE). Cette étude montre qu'un petit changement dans les conditions dynamique du fluide, induite par les différents modèles de sous-maille utilisés, peut provoquer une transition entre deux régimes d'écoulement distincts lorsque l'écoulement tourbillonnaire est proche des conditions critiques de transition. La sensibilité de la SGE aux modèles de sous-maille est également identifiée comme le résultat d'un manque de résolution à certains endroits critiques, un problème qui est analysé en utilisant une méthode d'adaptation de maillage. L’adaptation de maillage est testée sur des cas académiques et industriels. Ici, par ajustement de la résolution du maillage sur la base des caractéristiques de l'écoulement étudié (raffinement et grossissement de la grille en maintenant constant le coût numérique), des améliorations substantielles peuvent être obtenues, en terme de prédictions de la SGE. Ce travail peut être considéré comme une des premières étapes vers la mise en place d'une procédure standard (reproductible et indépendante de l’utilisateur) de maillage pour la SGE. / Swirling flows, which are widely employed in gas turbines, are known to undergo bifurcation between different topologies (large reconfigurations of the flow field) affecting the engine performance and safety. This work focuses on the study of such bifurcations using Large-Eddy Simulation (LES). It shows that a small change in the fluid dynamics conditions, induced by the different Sub-Grid Scale (SGS) models used in the simulations, can cause a transition between two, distinct, flow states when the swirling flow is close to transition conditions. The sensitivity of LES to SGS modeling is also identified as the result of a lack of mesh resolution at some critical locations, a problem which is analyzed using mesh adaptation. Mesh adaptation is tested on canonical and industrial flows. Here, by adjusting the mesh resolution based on the characteristics of the flow examined (refining and coarsening the grid keeping constant the numerical cost), substantial improvements of the LES predictions can be obtained. This work can be considered as the first step toward the establishment of a standard (repeatable and user independent) meshing procedure for LES.

Les objectifs ambitieux, fixés aux acteurs du secteur automobile par les pouvoirs publics, en matière d'émission de polluants et de gaz à effet de serre rendent aujourd'hui indispensable une compréhension plus fine de la combustion dans les moteurs. La simulation 3D aux grandes échelles (LES) représente une voie prometteuse pour répondre à ces enjeux. Elle permet l'étude de phénomènes transitoires complexes inaccessibles avec des moyens expérimentaux ou des méthodes de calculs traditionnelles de type RANS. Ce travail de thèse est une première étape vers la simulation LES de l'injection de carburant liquide dans les moteurs à piston. Il a consisté à adapter le code de calcul aux particularités physiques de l'injection directe, technologie qui se généralise actuellement à tous les types de moteurs à piston. Dans un premier temps, et afin de s'affranchir du calcul 3D complexe en sortie d'injecteur, une méthodologie originale, consistant à initier le calcul en aval de l'injecteur, est proposée et validée sur différents cas. Pour la simulation 3D, l'approche Eulérienne mésoscopique, à laquelle est ajouté un modèle d'interaction particules-particules, est utilisée pour simuler le spray. Les simulations ont été premièrement validées par comparaison expérimentale dans des conditions proches de l'injection Diesel. De plus, une étude sur la dynamique du spray a permis de mieux comprendre son évolution et de dégager des points communs avec un jet de gaz turbulent. Des simulations complémentaires ont également montré la prédictivité de la LES sur des injections Diesel réalistes. Enfin, un premier calcul moteur à injection directe a été réalisé et a permis de valider les développements réalisés dans le cadre de cette thèse. / Car manufacturers are facing increasingly severe regulations on pollutant emissions and fuel consumption. To respect these regulations, a better understanding of combustion processes is needed. Large Eddy Simulation (LES) is becoming a promising tool for such issues as it allows the study of complex unsteady phenomena which can not be analysed with RANS simulations or experiments. The present work is a step towards the LES of liquid injection in piston engines. The numerical code has been adapted to the specifications of Direct Injection which is more and more used in industry. Firstly, in order to avoid the difficulties linked to the 3D simulation of cavitation, primary break-up and turbulence in the near-nozzle region, an original methodology, based on an injector model, has been proposed. The idea is to initiate the spray physics downstream to the injector exit. Then LES 3D simulations of spray have been conducted using the Eulerian Mesoscopic approach extended to dense dispersed sprays by the addition of a particle-particle interactions model. The simulation results have been validated by comparison with experimental data in Diesel conditions with a low injection pressure. Furthermore a study on the spray dynamics has permitted to better understand its development and to find similarities with a turbulent gaseous jet. Additional simulations on realistic Diesel injection conditions have shown the good predictivity of LES in such cases. Finally, a first simulation of a Direct Injection Engine has been been carried out to assess the developments achieved in this work.

L'objectif de la présente étude est d'analyser la modélisation de la turbulence de simulations en moyenne de Reynolds (RANS) dans le cadre d'écoulements de type turbomachines, en utilisant des simulations aux grandes échelles (SGE) comme référence. L'étude porte sur deux cas test: un décollement de coin dans une grille d'aubes rectiligne, et un écoulement de jeu pour un aubage isolé dans un jet. Deux lois de comportement, la loi de comportement de Boussinesq et la loi de comportement quadratique (quadratic constitutive relation ou QCR), sont analysées, avec deux versions du modèle de turbulence k-omega de Wilcox. Les lois de comportement étudiées reposent sur deux hypothèses: une hypothèse d'alignement entre le tenseur de Reynolds et un tenseur construit à partir de l'écoulement moyen, et une hypothèse sur la viscosité turbulente. L'hypothèse d'alignement est étudiée à partir de la SGE, pour laquelle les deux tenseurs sont indépendamment connus, en utilisant un indicateur construit sur le produit scalaire des tenseurs. Les résultats sont présentés sous forme d'une fonction de répartition de la valeur de l'indicateur pour le domaine complet, puis pour trois sous-domaines d'intérêt: l'entrée, une région où l'écoulement interagit fortement avec les parois, et une région où l'écoulement est fortement tourbillonnaire. L'hypothèse d'alignement n'est que rarement valide pour la loi de comportement de Boussinesq. Pour la QCR, les résultats sont meilleurs en entrée, comparé à la loi de Boussinesq. Il ne sont cependant pas meilleurs pour les régions où l'écoulement est plus tourbillonnaire. Une amélioration de la loi de comportement est nécessaire pour pouvoir faire progresser la modélisation turbulente en RANS. En revanche, l'utilisation de l'énergie cinétique turbulente et du taux de dissipation spécifique semble correcte pour estimer la valeur de la viscosité turbulente. L'analyse de la modélisation de l'équation d'énergie cinétique turbulente (ECT) est réalisée au travers d'une comparaison terme à terme avec l'équation d'ECT résolue par la SGE. Les résultats SGE présentent une turbulence qui n'est pas à l'équilibre : la production et la dissipation ne sont pas superposées, et le terme de transport est important. Pour le RANS, la turbulence est à l'équilibre : la production et la dissipation sont superposées, et le terme de transport est de faible intensité. Un modèle de turbulence qui prend en compte le déséquilibre est nécessaire pour améliorer ce point. En dernier lieu, une nouvelle formulation hybride RANS/SGE est proposée, fondée sur la distance à la paroi en unités de paroi. La formulation est validée dans un canal bi-périodique et un premier essai est réalisé sur le cas de décollement de coin, mais d'autres analyses sont nécessaires avant que cette formulation ne soit fonctionnelle. / The present study aims at analysing turbulence modelling in Reynolds-averaged Navier-Stokes (RANS) simulations, in the context of turbomachinery flows, using large-eddy simulations (LES) as references. Two test cases are considered: a corner separation (CS) flow in a linear compressor cascade, and a tip-leakage (TL) flow of a single blade in a jet. Two constitutive relations, the Boussinesq constitutive relation and the quadratic constitutive relation (QCR), are investigated, with two versions of Wilcox's $k-\omega$ turbulence model. The studied constitutive relations rely on two hypotheses: an alignment hypothesis between the Reynolds stress tensor and a mean flow tensor, and an hypothesis on the turbulent viscosity. The alignment hypothesis is investigated using LES, where both the tensors are known independently, with an indicator built on the inner product of the tensors. The results are presented as probability density functions of the indicator value for the entire domain first, and then for three specific areas of interest: the inlet area, similar to a boundary-layer flow, an area of strong interaction between the flow and the walls (CS: passage area, TL: tip clearance) and an area of highly vortical flow (CS: separation wake, TL: tip-leakage vortex). The alignment hypothesis is rarely verified in any area for the Boussinesq constitutive relation. For the QCR, the results are improved for the inlet areas compared to the Boussinesq constitutive relation, but no significant improvement is found in the highly vortical regions. An improvement of the constitutive relation is needed in order to improve the RANS turbulence modelling. In contrast, the use of the turbulent kinetic energy and the specific dissipation rate appears quite correct to estimate the turbulent viscosity. The modelling of the RANS turbulent kinetic energy (TKE) budget equation is investigated through a term to term comparison with the resolved LES TKE budget equation. The LES presents a turbulence that is not at equilibrium, with the production and the dissipation not superimposed, and an important amount of transport. This differs from the RANS models, at equilibrium: the production and the dissipation are superimposed, with a small amount of transport. The development of a non-equilibrium turbulence model for RANS simulations could improve this aspect of turbulence modelling. Finally, a new hybrid RANS-LES formulation, based on the wall distance in wall units, is also proposed. It is validated on a bi-periodical channel flow, and a first attempt is made on the corner separation case, but further investigations are still needed for the model to be fully operational.