Les méthodes hybrides Moments/PDF sont bien adaptées pour la description des écoulements diphasiques turbulents, polydispersés, hors équilibre thermodynamique. Ces méthodes permettent d'avoir une description assez fine de la polydispersion, de la convection et des termes sources non-linéaires. Cependant, les approximations issues de telles simulations sont bruitées ce qui, dans certaines situations, occasionne un biais. L'approche alternative étudiée dans ce travail consiste à coupler une description Eulerienne des moments avec une description stochastique Lagrangienne à l'intérieur de la phase dispersée, permettant de réduire l'erreur statistique et d'éliminer le biais. La mise en oeuvre de cette méthode nécessite le développement de schémas numériques robustes. Les approches proposées sont basées sur l'utilisation simultanée des techniques de relaxation et de décentrement, et permettent d'obtenir des approximations stables des solutions instationnaires du système d'équations aux dérivées partielles, avec des données peu régulières provenant du couplage avec le modèle stochastique. Une comparaison des résultats de la méthode hybride Moments-Moments/PDF avec ceux issus de la méthode hybride "classique'' est présentée en termes d'analyse des erreurs numériques sur un cas de jet co-courant gaz-particules. / Hybrid Moments/PDF methods have shown to be well suitable for the description of polydispersed turbulent two-phase flows in non-equilibrium which are encountered in some industrial situations involving chemical reactions, combustion or sprays. hey allow to obtain a fine enough physical description of the polydispersity, non-linear source terms and convection phenomena. However, their approximations are noised with the statistical error, which in several situations may be a source of a bias. An alternative hybrid Moments-Moments/PDF approach examined in this work consists in coupling the Moments and the PDF descriptions, within the description of the dispersed phase itself. This hybrid method could reduce the statistical error and remove the bias. However, such a coupling is not straightforward in practice and requires the development of accurate and stable numerical schemes. The approaches introduced in this work rely on the combined use of the upwinding and relaxation-type techniques. They allow to obtain stable unsteady approximations for a system of partial differential equations containing non-smooth external data which are provided by the PDF part of the model. A comparison of the results obtained using the present method with those of the ``classical'' hybrid approach is presented in terms of the numerical errors for a case of a co-current gas-particle wall jet.

Le polymer flooding est une des techniques de récupération assistée des hydrocarbures (RAH) ; elle consiste à injecter une solution de polymères de forte masse moléculaire afin de déplacer plus efficacement le pétrole emprisonné dans la roche. Cependant, une limite importante de cette technique est la possible dégradation mécanique des polymères au cours de l'injection et dans le réservoir, due à une scission des chaînes macromoléculaires induite par l’écoulement. Ce travail de thèse a pour objectif de mieux comprendre les mécanismes et scénarios de scission, mais aussi leur impact sur le procédé de polymer flooding. Nous nous sommes intéressés au seuil de dégradation mécanique de solutions de poly(oxyde d’éthylène) et de de polyacrylamide partiellement hydrolysé, pour différents régimes de concentration (solutions diluées et semi-diluées) en régime laminaire et inertiel, et pour des solvants de différentes qualités. L’étude de la dégradation mécanique des solutions et de leur impact sur les propriétés rhéologiques a été menée à l’aide de différents dispositifs de dégradation et de différents rhéomètres, dont un dispositif microfluidique en élongation ; ces techniques de mesure ont été combinées à des mesures de distribution de masses moléculaires par chromatographie d’exclusion stérique couplée à la diffusion de lumière. L’étude montre en premier lieu qu’une composante élongationnelle est indispensable pour dégrader les chaînes macromoléculaires en solution. Les résultats mettent aussi clairement en évidence que les mécanismes de dégradation sont très différents en régime dilué et semi-dilué. En régime dilué, la dégradation mécanique des solutions de polymères est indépendante du régime d’écoulement et affecte préférentiellement les macromolécules de fortes masses, avec une scission en milieu de chaîne. En revanche, en régime semi-dilué, la dégradation mécanique dépend du régime de l’écoulement : en écoulement laminaire, la dégradation est gouvernée par le réseau d’enchevêtrements et la scission des chaînes est aléatoire, tandis qu’en régime inertiel, les chaînes se dégradent comme en régime dilué, avec le même scénario de scission en milieu de chaîne. Par ailleurs, les résultats montrent que les propriétés rhéologiques en élongation peuvent être très fortement impactées par la dégradation mécanique. Enfin, les résultats de l’étude préliminaire des propriétés d’injectivité dans un milieu poreux d’une solution de polymère semi-diluée faiblement dégradée montrent que la dégradation mécanique améliore l’injectivité du polymère aux abords du puits. / Polymer flooding is a technique used in enhanced oil recovery; it consists in injecting high molecular weight polymer solutions in order improve oil sweep efficiency in the reservoir. However, polymer flooding is challenged by possible mechanical degradation of polymer solutions during injection and in the reservoir, due to the flow induced scission of macromolecules. This work aims at better understanding the scission mechanisms and scenarios, but also their impact on polymer flooding. We investigated the onset of mechanical degradation of poly(ethylene oxide) and partially hydrolysed polyacrylamide solutions, for different concentration regimes (dilute and semi-dilute regimes), under laminar or inertial conditions, but also under good or bad solvent conditions. The study of mechanical degradation of polymer solutions and their impact on the rheological properties was performed using different degradation devices and different rheometers, including a microfluidic extensional device; these investigation techniques were combined with measurements of the molecular weight distributions using size exclusion chromatography coupled with light scattering experiments. The study first shows that an extensional component is needed to get a mechanical degradation of polymer chains. The results also clearly show that the degradation mechanisms are very different in dilute and semi-dilute regime. In dilute regime, the mechanical degradation of polymer solutions does not dependent on flow regime and mainly affects the macromolecules with high molecular weights, with a mid-chain scission scenario. On the other hand, in semi-dilute regime, mechanical degradation depends on flow regime: in laminar flows, degradation is governed by the entanglement network and chain scission is random, whereas in inertial flows, chain degradation is similar to that observed in dilute regime, with the same mid-chain scission scenario. Besides, the results show that the extensional rheological properties can be very strongly affected by mechanical degradation. At last, the results of a preliminary study of the injectivity of a slightly degraded semi-dilute polymer solution in porous media show that mechanical degradation improves polymer injectivity near the wellbore.

Cette thèse est consacrée à la simulation d'écoulements turbulents, incompressibles ou à faible nombre de Mach pour des applications touchant à la sûreté nucléaire. En particulier, nous nous concentrons sur le développement et l'analyse mathématique de schémas numériques performants pour la méthode dite de Simulation des Grandes Echelles. Ces schémas sont basés sur des méthodes à pas fractionnaires de type correction de pression et des éléments finis non conformes de bas degré. Deux arguments semblent essentiels à la construction de tels schémas: le contrôle de l'énergie cinétique et la précision pour des écoulements à convection dominante. Concernant la discrétisation en temps, nous proposons un schéma de type Crank-Nicolson et nous montrons qu'il satisfait un contrôle de l'énergie cinétique. Ce schéma présente de plus l'avantage d'être peu dissipatif numériquement (résidu d'ordre deux en temps). Concernant le défaut de précision de la discrétisation par l'élément fini de Rannacher-Turek, nous envisageons deux approches. La première consiste à construire un schéma pénalisé contraignant les degrés de liberté tangents aux faces des cellules à s'écrire comme combinaison linéaire des degrés de liberté normaux alentour. La deuxième approche repose sur l'enrichissement de l'espace discret d'approximation pour la pression. Enfin, différents tests numériques sont présentés en dimensions deux et trois et pour des maillages généraux, afin d'illustrer les capacités des schémas étudiés et de confronter les résultats théoriques et expérimentaux. / This thesis is devoted to the simulation of incompressible or low Mach turbulent flows, for nuclear safety applications. In particular, we focus on the development and analysis of performing numerical schemes for the Large Eddy Simulation technique. These schemes are based on fractional step methods of pressure correction type and on nonconforming low degree finite elements. Two requirements seems essential to build such schemes, namely a control of kinetic energy and the accuracy for convection dominated flows. Concerning the time marching algorithm, we propose a Crank-Nicolson like scheme for which we prove a kinetic energy control. This scheme has the advantage to be numerically low dissipative (numerical dissipation residual is second order in time). Concerning the low accracy of the Rannacher-Turek discretization, two approaches are investigated in this work. The first one consists in building a penalized scheme constraining the velocity degrees of freedom tangent to the faces to be written as a linear combination of the normal ones. The second approach relies on the enrichment of the pressure approximation discrete space. Finally, various numerical tests are presented in both two and three dimensions and for general meshes, to illustrate the capacity of the schemes and compare theoretical and experimental results.

Cette thèse est consacrée à l'étude de méthodes numériques pour la simulation des écoulements polyphasiques en milieu poreux, en vue de leur application à des problèmes d'ingénierie pétrolière ou environnementale. Nous présentons une formulation générique du modèle d'écoulements à nombre quelconque de composants présents dans un nombre quelconque de phases. Dans notre approche l'approximation des flux diffusifs (issus, par exemple, de la loi de Darcy) s'appuie sur de nouveaux schémas, appelés schémas gradient, qui ont plusieurs avantages sur les schémas industriels standard : ces derniers, qui sont des schémas volumes finis multi-points centrés aux mailles, ne sont généralement pas symétriques et convergent difficilement sur des cas à forts rapports d'anisotropie. Nous montrons en revanche que les schémas gradient conduisent naturellement à des approximations symétriques et convergentes. Parmi cette classe de schémas, nous étudions plus particulièrement le schéma "VAG" qui fait intervenir des inconnues au centre des mailles et aux sommets du maillage. Ce schéma conduit à la définition de flux entre le centre d'une maille et ses sommets, qui sont utilisés pour généraliser la méthode "VAG" au contexte polyphasique. Des tests numériques montrent alors que ce schéma est robuste, et conduit à un très bon compromis précision/coût, ce qui en fait un candidat idoine pour les applications industrielles. Nous présentons notamment un cas test, basé sur des observations de terrains, d'injection et de dissolution de CO2 dans la région proche d'un puits foré dans un aquifère salin. Nous montrons alors que le schéma numérique permet de simuler l'assèchement et la précipitation de minéral observée en pratique. Un chapitre de la thèse est enfin consacré à l'étude pratique et théorique d'une méthode numérique générique pour contrôler l'effet d'axe lors de l'utilisation de schémas industriels / This thesis is focused on numerical methods dedicated to the simulation of multiphase flow in porous media, involved in petroleum or environmental engineering. We present a generic formulation of the flow model which is able to take into account any number of components within any number of phases. In our approach the approximation of the diffusive fluxes (mainly resulting from Darcy's law) is based on new schemes, called gradient schemes, which show several advantages over the standard industrial numerical schemes : these schemes, which belong to the class of the cell-centred MultiPoint Flux Approximation finite volume schemes, are not symmetric and may lead to difficulties of convergence in the case of high anisotropy ratios. We indeed show that gradient schemes are naturally providing symmetric and convergent approximations. We particularly study the Vertex Approximate Gradient scheme (called the VAG scheme), whose discrete unknowns are the values at the cell centres and at the vertices of the mesh. This scheme implies the definition of fluxes between the centre of a given cell and its vertices, used for the extension of the scheme to multiphase flow. Numerical tests show the robustness and the accuracy of the method for a low computational cost, which enables the use of the VAG scheme in an industrial framework. A test case, based on experimental data of injection and dissolution of CO2 in the near-well region within a saline aquifer, shows the aptitude of the scheme for reproducing drying and salt precipitation, which are practically observed. Finally, a chapter of the thesis is devoted to the theoretical and practical study of a general numerical method for controlling Grid Orientation Effect in industrial simulators

Le modèle numérique que nous avons développé au cours de cette thèse présente deux caractéristiques principales : un modèle dilatable pour l'eau et la prise en compte de domaines ouverts. Les difficultés associées au premier aspect concernent l'adaptation de la loi d'état de l’eau au modèle dilatable sous l’approximation à faibles nombres de Mach, tandis que celles associées au second sont relatives à la mise en œuvre de conditions aux limites numériques de sortie compatibles avec l'algorithme de projection utilisé. Les résultats de simulations d'écoulement de convection mixte en canal horizontal chauffé par le bas ont été confrontés à celles utilisant l'approximation de Boussinesq et aux expériences. / The 3D numerical model which we developed in this thesis presents two main features: a Low-Mach-Number approximation for water along with an open boundary condition formulation. Indeed, the difficulties related to the former point stand in a computationally efficient adaptation of the water equation of state in the framework of Low Mach number approximation, whereas the difficulties related to the latter concern the introduction of Open Boundary Conditions in the projection algorithm used. We have computed a mixed convection flow in a horizontal channel uniformly heated from below and compared the results obtained with both the Boussinesq approximation and experimental results.

Une grande quantité d'énergie est présente dans les mouvements de rotations propre et orbitale des planètes. Des forçages harmoniques tels que les déformations de marées, la précession ou la libration peuvent en convertir une partie pour générer des écoulements dans les couches fluides d'une planète. Ces écoulements restent largement méconnus même s'ils sont importants pour contraindre des modèles d'intérieur planétaire ou expliquer la présence de champs magnétiques dans certains astres.Dans cette thèse, nous étudions les mécanismes engendrés par ces forçages en combinant une approche théorique, expérimentale et numérique et soulignons la généricité des phénomènes observés. L'étude d'un forçage de libration longitudinale, i.e. des oscillations de la vitesse de rotation d'un astre, montre la présence d'un écoulement zonal généré par des interactions non-linéaires dans les couches visqueuses. Nous étudions ensuite l'instabilité qui apparaît à la paroi pour des amplitudes de libration suffisantes et qui peut transférer de l'énergie vers le volume du fluide. Finalement, une étude expérimentale de forçage de marées dans une sphère met en évidence que l'excitation directe d'ondes inertielles induit un écoulement zonal intense et localisé. Cet écoulement peut se déstabiliser par une instabilité de cisaillement et générer un écoulement turbulent dans tout le volume.Pour finir, nous considérons la pertinence de ces résultats pour des applications géo-/astrophysiques, telles que l'étude des océans internes sous la surface de glace des satellites joviens Ganymède, Encelade et Europe. / A huge amount of energy is stored in the spin and orbital motions of any planet. Harmonic forcings such as libration, precession and tides are capable of conveying a portion of this energy to drive intense three-dimensional flows in liquid layers of planetary bodies. The generated flows remain largely unknown even if they are important to constraint model of planetary interior or to explain the presence of magnetic fields in some astrophysical bodies.In this thesis, we study the mechanisms induced by these forcings by combining theoretical, experimental and numerical approaches and we highlight the genericity of the observed phenomena. The study of a longitudinal libration forcing, corresponding to oscillations of the rotation rate of a planet, shows the presence of a mean zonal flow generated by non-linear interactions in the viscous layers. We then study the instability which appears at the outer boundary at sufficiently large libration amplitude or small Ekman number and which can transfer energy to the bulk of the fluid. Finally, an experimental study of tidal forcing in a sphere shows that the nonlinear self-interaction of excited inertial waves may drive an intense and localised axisymmetric jet, which becomes unstable at low Ekman number following a shear instability, generating space-filling turbulence.To conclude, we consider the relevance of these results to geo-/astrophysical applications, such as the subsurface oceans of the icy satellites Ganymede, Enceladus or Europa.

Cette thèse s'intéresse à la simulation numérique de l'explosion vapeur. Ce phénomène correspond à une vaporisation instantanée d'un volume d'eau liquide entraînant un choc de pression. Nous nous y intéressons dans le cadre de la sûreté nucléaire. En effet, lors d'un accident entraînant la fusion du cœur du réacteur, du métal fondu pourrait interagir avec de l'eau liquide et entraîner un tel choc. On voudrait alors connaître l'ampleur de ce phénomène et les risques d'endommagements de la centrale qu'il implique. Pour y parvenir, nous utilisons pour modèle les équations d'Euler dans un cadre Lagrangien. Cette description a l'avantage de suivre les fluides au cours du temps et donc de parfaitement conserver les interfaces entre l'eau liquide et sa vapeur. Pour résoudre numériquement les équations obtenues, nous développons un nouveau schéma de type Godunov utilisant des flux nodaux. Le solveur nodal développé durant cette thèse ne dépend que de la répartition angulaire des variables physiques autour du nœud. De plus, nous nous intéressons aux changements de phase liquide-vapeur. Nous proposons une méthode pour les prendre en compte et mettons en avant les avantages qu'il y a à l'implémentation de ce phénomène dans un algorithme Lagrangien. / This thesis studies numerical simulation of steam explosion. This phenomenon correspond to a fast vaporization of a liquid leading to a pressure shock. It is of interest in the nuclear safety field. During a core-meltdown crisis, molten fuel rods interacting with water could lead to steam explosion. Consequently we want to evaluate the risks created by this phenomenon.In order to do it, we use Euler equations written in a Lagrangian form. This description has the advantage of following the fluid motion and consequently preserves interfaces between the liquid and its vapor. To solve these equations, we develop a new Godunov type scheme using nodal fluxes. The nodal solver developed here only depends on the angular repartition of the physical variables around the node.Moreover, we study liquid-vapor phase changes. We describe a method to take it into account and highlight the advantages of using this method into a Lagrangian framework.

Nous nous sommes intéressés dans la présente thèse à l’étude de l’écoulement gravitaire de nappesliquides non guidées qui s’écoulent verticalement dans l’air ambiant. Après une synthèse bibliographique,nous avons réalisé une double étude, expérimentale et numérique, en considérant des liquidesde viscosités différentes (allant de 1 à 50 fois celle de l’eau) et de tension superficielle proche de cellede l’eau. Le nombre de Reynolds Rel du liquide a varié de quelques unités à quelques milliers alorsque le nombre de Weber du liquide allait approximativement de 0.1 à 10. Le dispositif expérimentalque nous avons mis en place nous a permis de créer et d’étudier les formes géométriques des nappesliquides. Nous avons employé une méthode expérimentale originale pour mesurer le champ d’épaisseurdes nappes. Nous avons mené, parallèlement aux expériences, des simulations numériques 2D et 3Dinstationnaires et diphasiques (VOF), utilisant le calcul parallèle. Nous avons trouvé que les nappesexpérimentales et numériques deviennent plus courtes (verticalement) et moins épaisses lorsque le débitdiminue. Expérimentalement, lorsque le débit du liquide devient suffisamment faible, des filamentsliquides commencent à apparaître à coté d’une nappe moins large qu’auparavant. Dans le cas desnappes d’eau, cette transition de régime d’écoulement a été précédée de l’apparition systématique detrous dans la partie inférieure des nappes. Pour les autres liquides newtoniens plus visqueux (solutionsaqueuses de glycérine), l’apparition des filaments liquides a été précédée d’une déstabilisation des bourreletsqui délimitent la partie plane des nappes. Nous avons étendu par la suite l’étude expérimentaleà celle de fluides au comportement rhéologique plus complexe en utilisant un liquide non newtonienrhéofluidifiant. Nous avons montré pour ce fluide rhéofluidifiant que le débit associé à la transition versle régime des filaments diminue en comparaison avec celui associé à un liquide newtonien de viscositésimilaire. Cela nous a conduit à suggérer que la présence des propriétés rhéofluidifiantes des nappesliquides peut représenter une solution pour les applications de "coating" pour lesquelles on cherche àproduire des nappes stables et sans percement dans des configurations d’écoulement de faibles débits. / In this thesis, unguided plane liquid sheets flowing vertically by gravity in an ambient air atmosphereare studied experimentally and numerically. First of all a litterature survey clearly identified themain issues regarding the dynamics and modelling of such flows. Subsequently, different liquids exhibitinga wide range of viscosity (1 to 50 times that of water) and a surface tension close to that of waterwere selected. The liquid flow regimes were characterized by a Reynolds number Rel ranging from afew units to a few thousand while the Weber number Wel was varied between 0.1 to 10. A dedicatedexperimental system was designed and operated to study the relevant sheet features (geometry, thickness)by means of an original optical method. In parallel, finite volume based 2D and 3D simulationsof the flows were undertaken. All rely on the volume of fluid method (VOF) combined with adaptivemeshing. The experimental and numerical sheets became shorter (vertically) and thinner as the massflow rate decreased. Experimentally, when the mass flow rate of the liquid becomes sufficiently low,liquid threads begin to appear next to a narrower sheet than before. In the case of water, this flowregime transition was preceded by the systematic appearance of holes in the lower part of the sheets.For the other more viscous Newtonian liquids (mixtures of water and glycerin), the appearance of theliquid threads was preceded by a destabilization of the rims which delimited the flat part of the sheets.The experimental study was then extended to fluids featuring more complex rheological behavior e.g.by the use of a non-Newtonian shear-thinning fluid. For such a fluid, it was shown that the criticalmass flow rate associated with the transition towards the threads regime was lower than its Newtoniancounterpart of similar viscosity. It is suggested that the presence of shear-thinning properties in liquidsheets may represent a solution for "coating" applications for which stable and non-pierced curtainsin flow configurations of low mass flow rates are targeted.

La thèse a pour but de développer une méthodologie de calcul instationnaire permettant une étude qualitative et quantitative de l’aérodynamique externe d’une automobile. La première partie de l’étude est consacrée au développement de la méthodologie numérique sur les différents corps d’Ahmed à 25°, géométries simplifiées d’une automobile réelle, afin de valider les choix stratégiques de maillages et de méthodes numériques ainsi que s’assurer de la bonne résolution de l’écoulement. Les résultats numériques sont comparés aux données expérimentales obtenues au cours d’essais réalisés à la soufflerie La Ferté Vidame lors de travaux de thèse précédents. A l’issue de ce travail, la méthode hybride Delayed Detached Eddy Simulation Shear-Stress Transport (DDES SST) est sélectionnée pour la suite de l’étude au vu des meilleures performances obtenues (torseur aérodynamique, coefficient de pression, topologie d’écoulement) par rapport aux méthodes Reynolds Averaged Navier-Stokes (RANS), Scale Adaptive Simulation (SAS) et Stress-Blended Eddy Simulation (SBES). La seconde phase de l’étude consiste à adapter la méthodologie précédemment développée sur un véhicule réel, la Peugeot 308 SW. Au préalable, une base de données expérimentales conséquente a été réalisée au sein du Groupement d’Intérêt Economique Souffleries Aéroacoustiques Automobiles (GIE S2A) au cours de ces travaux. La géométrie est tout d’abord simplifiée afin de faciliter la mise en place de la méthodologie numérique : entrées d’air fermées, soubassement lissé, roues remplacées par des carénages. Les résultats obtenus sont encourageants et démontrent globalement la supériorité de la DDES par rapport aux méthodes RANS classiques. La topologie d’écoulement est mieux prédite (soubassement et sillage), même si la prédiction du coefficient de portance reste une difficulté majeure pour ce type de méthode hybride. / The main goal of this PhD is to develop an unsteady numerical method to study the external aerodynamic flow around real vehicles. The first part of the study focuses on the flow around simplified geometries, such as 25° Ahmed bodies (with sharp and rounded edges on the back of the body), in order to determine the optimal turbulence model, mesh setup and numerical parameters. Computational Fluid Dynamics (CFD) results are compared to experimental data reported in literature conducted in the La Ferté Vidame wind tunnel. Based on this study, Shear-Stress Transport Delayed Detached Eddy Simulation (SST DDES) demonstrates superiority over Reynolds Averaged Navier-Stokes (RANS), Scale Adaptive Simulation (SAS) and Stress Blended Eddy Simulation (SBES) turbulence models, regarding both drag and lift coefficients predictions, and flow topology.Secondly, the numerical procedure is adapted for a real vehicle, the Peugeot 308 SW estate car. A substantial experimental campaign was carried out in the Groupement d’Intérêt Economique Souffleries Aéroacoustiques Automobiles (GIE S2A) wind tunnel to provide data against which the numerical results are compared. Given the geometric complexity of a real vehicle, the car is simplified for this study as follows: the front air inlets are closed, the underbody is smoothed with additional panels and the wheels are replaced by fairings. DDES computations show encouraging results. A significant improvement of the flow topology is obtained with DDES compared to RANS models. However, the prediction of the lift coefficient remains a major difficulty with these hybrid methods.

Afin de faire face aux fortes températures rencontrées par les composantsen aval de la chambre de combustion, des prélèvements d’air plus frais sont réalisésau niveau du compresseur. Cet air alimente les cavités en pied de turbine et refroidiles disques rotor permettant d’assurer le bon fonctionnement de la turbine.Ce manuscrit présente une étude numérique de l’effet de ces écoulements de cavitéau pied de la turbine sur ses performances aérodynamiques. Les phénomènesd’interaction entre l’air de cavité en pied de turbine et l’air de veine principal est unphénomène encore difficilement compris. L’étude de ces phénomènes est réalisée autravers de différentes approches numériques (RANS, LES et LES-LBM) appliquéesà deux configurations pour lesquelles des résultats expérimentaux s ont disponibles.Une première configuration en grille d’aube linéaire en amont de laquelle différentesgéométries d’entrefer (interface entre plateforme rotor et stator) et débits de cavitépouvaient être variés. Une seconde configuration annulaire composée de deux étagesde turbine comprenant les cavités en pied et plus proche d’une configuration industrielle.Les pertes additionnelles associées à l’écoulement de cavité sont mesurées etétudiées à l’aide d’une méthode basée sur l’exergie (bilans d’énergie dans l’objectifde générer du travail). / In order to deal with high temperatures faced by the components downstreamof the combustion chamber, some relatively cold air is bled at the compressor.This air feeds the cavities under the turbine main annulus and cool down the rotordisks ensuring a proper and safe operation of the turbine. This thesis manuscriptintroduces a numerical study of the effect of the cavity flow close to the turbine hubon its aerodynamic performance. The interaction phenomena between the cav-ity andmain annulus flow are not currently fully understood. The study of these phenomenais performed based on different numerical approaches (RANS, LES and LES-LBM)applied to two configurations for which experimental results are avail-able. A linearcascade configuration with an upstream cavity and various rim seal geometries(interface between rotor and stator platform) and cavity flow rate avail-able. Arotating configuration that is a two stage turbine including cavities close to realisticindustrial configurations. Additional losses incurred by the cavity flow are measuredand studied using a method based on exergy (energy balance in the purpose togenerate work).