Ce travail de recherche s’intéresse au problème de la stabilisation de flammes, et du transfert de chaleur résultant, dans les moteurs fusées cryogéniques. La stabilisation de flamme dans un moteur fusée est un phénomène critique, et toute instabilité peut conduire à des dégâts importants, voire à la destruction du lanceur et des satellites embarqués. Les moteurs (Vulcain 2 et Vinci) qui équipent Ariane 5, et la future Ariane 6, utilisent le couple hydrogène / oxygène, dont la grande impulsion spécifique permet, en minimisant la masse des ergols par rapport à la masse de la charge, une meilleure performance du lanceur. Pour réduire le volume de stockage, les ergols sont refroidis à des températures très basses de l’ordre de quelques dizaines de Kelvin. Ils alimentent une flamme dont la température maximale peut atteindre 3500K, générant de très forts gradients de température dans le fluide, et des flux de chaleur extrêmes dans les parties solides de l’injecteur. Pour prédire les flux de chaleur entre la flamme, l’injecteur et les ergols froids, l’approche de Simulation aux Grandes Echelles (SGE), pour reproduire l’écoulement réactif turbulent instationnaire, est couplée au calcul de thermique du solide dans l’injecteur. Cette approche est d’abord validée par comparaison à une expérience en conditions ambiantes, menée au Laboratoire EM2C (Paris). L’interaction flamme-paroi en présence de transfert de chaleur, qui est un mécanisme de base de la stabilisation de flamme, est ensuite étudiée pour différents niveaux de pression. Finalement une configuration représentative d’un injecteur coaxial de moteur fusée est simulée pour étudier la structure et les mécanismes de stabilisation de la flamme, ainsi que les flux de chaleur reçus par l’injecteur, en vue d’évaluer la fatigue thermique du système. / This research work deals with the problem of the flame stabilization in the context of high pressure liquid rocket engines. Flame stabilization in a rocket engine is a critical feature. An instability can lead to important damages of the engine or the destruction of the launcher and the satellite. The engines (Vulcain 2 and Vinci) of the Ariane 5, and the future Ariane 6, use the hydrogen/oxygen propellants. One characteristic of this couple is its high specific impulse. The launcher performance is linked to the ratio of the payload to the total mass of propellants. For volume reasons the propellants are stored at low temperature of the order of a few tens of Kelvin. When injected in the combustion chamber, their combustion releases a huge amount of heat leading to temperature of 3500K. In order to predict the heat transfer between the flame, the solid injector and the cold propellants the Large Eddy Simulation, which allows to capture the unsteady features of the flow, is used in association with a thermal solver for the injector. This approach is validated with a low pressure experiment conducted at Centrale Paris, then a basic 1D configuration allows to understand the phenomena of high pressure flame-wall interaction. Finally a configuration representative of a coaxial rocket engine injector allows to evaluate the structure and the stabilization mechanisms of a cryogenic flame, the heat flux and the temperature of the injector.

Afin de répondre aux enjeux environnementaux, les fabricants de turbine à gaz ont mis au point de nouveaux concepts de chambre de combustion plus propres et moins consommateurs. Ces technologies sont cependant plus sensibles aux instabilités de combustion, un couplage entre acoustique et flamme pouvant conduire à des niveaux dangereux de fluctuations de pression et de dégagement de chaleur. Les solveurs de Helmholtz sont une méthode numérique efficace pour prédire ces instabilités de combustion. Ils reposent sur la description d'un fluide non visqueux au repos, dont le comportement acoustique est régi par une équation d'Helmholtz thermoacoustique, résolue dans le domaine fréquentiel comme un problème aux valeurs propres. Le couplage flamme/acoustique est modélisé par une fonction de transfert du premier ordre entre les perturbations de dégagement de chaleur et la vitesse acoustique en un point de référence. Bien que performants, les solveurs de Helmholtz négligent l'interaction entre acoustique et vorticité aux coins, car celle-ci dépend d'effets visqueux. Cette interaction pourrait fortement amortir l'acoustique d'une chambre de combustion et la négliger revient à faire des prédictions trop pessimistes voire erronées. Par conséquent, une méthodologie a été mise au point afin d'inclure dans un solveur de Helmholtz l'effet d'interactions complexes entre acoustique et écoulement. Ces interactions étant compactes, elles sont modélisées par des matrices 2x2 et ajoutées au solveur comme des paires de conditions limites : les conditions limites de matrice. Grâce à cette méthodologie, les fréquences et modes d'une configuration académique non-réactive sont correctement calculées en présence de deux éléments où une telle interaction est forte: un orifice et un tourbilloneur. Afin d'être applicable aux chambres industrielles, deux extensions sont nécessaires. Premièrement, les surfaces de matrices doivent pouvoir être non-planes, afin de s'adapter aux géométries industrielles complexes. Pour cela, une procédure d'ajustement a été mise en place. La matrice est mesurée sur des surfaces planes et des transformations nondissipatives lui sont appliquées afin de la déplacer sur les surfaces non planes. Ces transformations peuvent être déterminées analytiquement ou calculées avec un solveur de propagation acoustique. Le deuxième point concerne le point de référence du modèle de flamme. En effet, celui-ci est souvent choisi à l'intérieur de l'injecteur ce qui pose problème si celui-ci est retiré du domaine de calcul et remplacé par sa matrice. Dans cette thèse, le point de référence est remplacé par une surface de référence. La méthodologie étendue est validée sur des configurations académiques puis appliquée à une chambre annulaire de Safran. Cette nouvelle méthodologie permet de constater que l'interaction écoulement/acoustique au niveau des trous de dilution et des injecteurs joue un effet important sur la stabilité de la chambre mais aussi sur la structure des modes. Les premiers résultats avec une surface de référence pour la flamme sont encourageants. / Environmental concerns have motivated turbine engine manufacturers to create new combustor designs with reduced fuel consumption and pollutant emissions. These designs are however more sensitive to a mechanism known as combustion instabilities, a coupling between flame and acoustics that can generate dangerous levels of heat release and pressure fluctuations. Combustion instabilities can be predicted at an attractive cost by Helmholtz solvers. These solvers describe the acoustic behavior of an inviscid fluid at rest with a thermoacoustic Helmholtz equation, that can be solved in the frequency domain as an eigenvalue problem. The flame/acoustics coupling is modeled, often with a first order transfer function relating heat release fluctuations to the acoustic velocity at a reference point. One limitation of Helmholtz solvers is that they cannot account for the interaction between acoustics and vorticity at sharp edges. Indeed, this interaction relies on viscous processes at the tip of the edge and is suspected to play a strong damping role in a combustor. Neglecting it results in overly pessimistic stability predictions but can also affect the spatial structure of the unstable modes. In this thesis, a methodology was developed to include the effect of complex flow-acoustic interactions into a Helmholtz solver. It takes advantage of the compactness of these interactions and models them as 2-port matrices, introduced in the Helmholtz solver as a pair of coupled boundary conditions: the Matrix Boundary Conditions. This methodology correctly predicts the frequencies and mode shapes of a nonreactive academic configuration with either an orifice or a swirler, two elements where flowacoustic interactions are important. For industrial combustors, the matrix methodology must be extended for two reasons. First, industrial geometries are complex, and the Matrix Boundary Conditions must be applied to non-plane surfaces. This limitation is overcome thanks to an adjustment procedure. The matrix data on non-plane surfaces is obtained from the well-defined data on plane surfaces, by applying non-dissipative transformations determined either analytically or from an acoustics propagation solver. Second, the reference point of the flame/acoustics model is often chosen inside the injector and a new reference location must be defined if the injector is removed and replaced by its equivalent matrix. In this work, the reference point is replaced by a reference surface, chosen as the upstream matrix surface of the injector. The extended matrix methodology is successfully validated on academic configurations. It is then applied to study the stability of an annular combustor from Safran. Compared to standard Helmholtz computations, it is found that complex flow-acoustic features at dilution holes and injectors play an important role on the combustor stability and mode shapes. First encouraging results are obtained with surfacebased flame models.

Les progrès du calcul scientifique ont permis des avancées importantes dans la simulation et la compréhension de problèmes complexes tels que les différents phénomènes physiques qui ont lieu dans des turbines à gaz industrielles. Cependant' l'essentiel de ces avancées portent sur la résolution d'un seul problème à la fois. En effet on résout soit les équations de la phase fluide d'un côté' de la thermique d'un autre' du rayonnement' etc... Pourtant' dans la réalité tous ces différents problèmes physiques interagissent entre eux: on parle de problèmes couplés. Ainsi en réalisant des calculs couplés on peut continuer à améliorer la qualité des simulations et donc donner aux concepteurs de turbines à gaz des outils supplémentaires. Aujourd'hui' des logiciels récents permettent de résoudre plusieurs physiques simultanément grâce à des solveurs génériques. En revanche' la contrepartie de cette généricité est qu'ils se révèlent peu efficaces sur des problèmes coûteux tels que la Simulation aux Grandes Echelles (SGE). Une autre solution consiste à connecter des codes spécialisés en leur faisant échanger des informations' cela s'appelle le couplage de codes. Dans cette thèse on s'intéresse au couplage d'un domaine fluide dans lequel on simule une SGE réactive (combustion) avec un domaine solide dans lequel on résout la conduction thermique. Pour réaliser ce couplage une méthodologie est mise en place en abordant différentes problématiques. Tout d'abord' la problématique spécifique au couplage de la SGE et de la thermique : l'impact de la fréquence d'échange sur la convergence du système ainsi que sur les problèmes de repliement de spectre et la stabilité du système couplé. Ensuite les problèmes d'interpolation et de géométrie sont traités avec notamment le développement d'une méthode d'interpolation conservative et la mise en évidence des difficultés spécifiques au couplage de géométries industrielles. Finalement la problématique du calcul haute performance (HPC) est traitée avec le développement d'une méthode permettant de réaliser efficacement l'échange des données et l'interpolation entre différents codes parallèles. Ces travaux ont été appliqués sur une configuration de chambre de combustion aéronautique industrielle. / Progress in scientific computing has led to major advances in simulation and understanding of the different physical phenomena that exist in industrial gas turbines. However' most of these advances have focused on solving one problem at a time. Indeed' the combustion problem is solved independently from the thermal or radiation problems' etc... In reality all these problems interact: one speaks of coupled problems. Thus performing coupled computations can improve the quality of simulations and provide gas turbines engineers with new design tools. Recently' solutions have been developed to handle multiple physics simultaneously using generic solvers. However' due to their genericity these solutions reveal to be ineffective on expensive problems such as Large Eddy Simulation (LES). Another solution is to perform code coupling: specialized codes are connected together' one for each problem and they exchange data periodically. In this thesis a conjugate heat transfer problem is considered. A fluid domain solved by a combustion LES solver is coupled with a solid domain in which the conduction problem is solved. Implementing this coupled problem raises multiple issues which are addressed in this thesis. Firstly' the specific problem of coupling an LES solver to a conduction solver is considered: the impact of the inter-solver exchange frequency on convergence' possible temporal aliasing' and stability of the coupled system is studied. Then interpolation and geometrical issues are addressed: a conservative interpolation method is developed and compared to other methods. These methods are then applied to an industrial configuration' highlighting the problems and solutions specific to complex geometry. Finally' high performance computing (HPC) is considered: an efficient method to perform data exchange and interpolation between parallel codes is developed. This work has been applied to an aeronautical combustion chamber configuration.

La nature conflictuelle des contraintes de performances, d'opérabilité et de respect des normes environnementales conduit les motoristes à optimiser finement la géométrie du brûleur afin d'identifier le meilleur design.La Simulation aux Grande Echelles (SGE) est aujourd'hui un outil performant et est déployé de manière courante dans les Bureaux d'Etudes pour la prédiction des propriétés macroscopiques de l'écoulement.Toutefois, de nombreux phénomènes influencés par les effets de chimie complexe, tels que la stabilisation, l'extinction de flamme et la formation des polluants, reste un problème crucial.En effet, la description des effets de chimie complexe nécessite l'utilisation de modèles cinétiques détaillés imposant des coûts de calculs prohibitifs, des problèmes de raideurs numérique et des difficultés de couplage avec les échelles non résolues turbulentes.Afin d'inclure une description des processus chimiques, dans les simulations numériques de chambres de combustion réelles, des modèles réduits doivent être proposés.Dans cette thèse, une méthode originale, appelée chimie virtuelle optimisée, est développée.Cette stratégie a pour objectif la description de la structure chimique de la flamme et la formation des polluants dans des configurations de flamme représentatives.Les schémas cinétiques virtuels optimisés, composés de réactions virtuelles et d'espèces virtuelles, sont construits par optimisation des paramètres réactionnels et des propriétés thermochimiques des espèces virtuelles afin de capturer les propriétés de flamme d'intérêt. / The conflicting nature of performance, operability and environmental constraints leads engine manufacturers to perform a fine optimization of the burner geometry to find the best design compromise.Large Eddy Simulation (LES) is an attractive tool to achieve this challenging task, and is routinely used in design office to capture macroscopic flow features.However, the prediction of phenomena influenced by complex kinetic effects, such as flame stabilization, extinction and pollutant formation, is still a crucial issue.Indeed, the comprehensive description of combustion chemistry effects requires the use of detailed models imposing prohibitive computational costs, numerical stiffness and difficulties related to model the coupling with unresolved turbulent scales.Reduced-cost chemistry description strategies must then be proposed to account for kinetic effects in LES of real combustion chambers.In this thesis an original modeling approach, called virtual optimized chemistry, is developed.This strategy aims at describing the chemical flame structure and pollutant formation in relevant flame configurations, at a low computational cost.Virtual optimized kinetic schemes, composed by virtual reactions and virtual species, are built through optimization of both kinetic rate parameters and virtual species thermo-chemical properties so as to capture reference target flame quantity.

L’application de cette thèse est le contrôle actif de la combustion dans les brûleurs industriels à combustible liquide. Il s’agit d’explorer les possibilités de contrôle d’un spray par des jets gazeux auxiliaires. Deux familles d’actionneurs utilisant ce procédé ont été testées sur un atomiseur coaxial assisté par air. Le premier dispositif est appelé (Dev). Composé d’un unique jet actionneur, il vise à dévier le spray. La seconde configuration, appelée (Sw), est équipée de 4 jets auxiliaires tangents au spray afin de lui conférer un effet de swirl et d’en augmenter le taux d’expansion. Les mesures de granulométrie par PDA et les visualisations du spray par strioscopie démontrent un effet important du contrôle sur l’atomisation et la forme du spray. On observe en outre une déviation pouvant atteindre 30°avec l’actionneur (Dev) et une augmentation du taux d’expansion de 80% dans le cas (Sw). Des simulations du banc expérimental ont de plus été menées avec le code AVBP. L’écoulement de gaz est calculé par simulation aux grandes échelles (SGE ou LES en Anglais). L’approche lagrangienne est utilisée pour simuler la phase dispersée. Une attention particulière a été portée aux conditions d’injection du gaz et des gouttes dans le calcul. Ceci a abouti au développement d’une nouvelle condition limite caractéristique non réfléchissante (VFCBC) destinée à l’injection d’écoulements turbulents en LES compressible. Les résultats de LES présentent un bon accord avec les mesures expérimentales. Les effets du contrôle sur la dynamique des gouttes et sur la topologie du spray (forme, déviation, expansion) sont correctement décrits. / The present work focuses on active control of two-phase combustion in industrial burners. The generic method explored in this thesis consists in controlling the injected fuel spray with transverse air jets. Two families of these jet actuators are tested on a coaxial airblast atomizer. The first system (Dev) is used to modify the trajectory of the spray, while the second one (Sw) introduces swirl into the spray to modify its spreading rate and mixing with the surrounding air. Experimental characterisations of the controlled flow with Schlieren visualisations and Phase Doppler Anemometry (PDA) show that actuators induce important effect on the spray. The deviation angle reaches 30° for the actuator (Dev) and the expansion rate increases of 80 % in the swirl case (Sw). Simulations of the experiment are then performed with the CFD code AVBP. The gas flow is computed with Large Eddy Simulation (LES). A Lagrangian formulation is used to simulate droplets trajectories. A particular attention is given to the injection of the gas flow and the droplets in the calculations. Therefore, a new non-reflecting characteristic boundary condition (VFCBC) has been derived to inject turbulent flows in compressible LES. A good agreement is observed between simulation and experiment. Control effects on the spray topology ( features, deviation, spread rate) and on the droplets velocities and diameters are correctly described by the Lagrangian LES.

Dans la première partie de cette thèse de doctorat une méthodologie est présentée qui permet de prédire les niveaux de suies produits dans des flammes laminaires monodimensionnelles, ou un modèle semi-empirique de suies est utilisé en combinaison avec une chimie complexe et un solveur radiatif détaillé. La méthodologie est appliquée au calcul de suies dans une série de flammes de diffusion à contre-courant d'éthylène/air. Plusieurs modèles d'oxydation de suies sont testés et les constantes du modèle sont ajustées afin de retrouver un meilleur accord avec les expériences. L'effet des pertes thermiques radiatives sur la formation de suies et la structure des flammes est évalué. Finalement, la performance du modèle de suies est évalué sur des flammes prémélangées monodimensionnelles, ou une expression alternative du terme de croissance de surface est proposée pour reproduire les résultats expérimentaux. Dans la deuxième partie de cette thèse, des outils de Simulation aux Grandes Échelles (SGE) et d'analyse acoustique sont appliqués à la prédiction des oscillations de cycle limite (OCL) d'une instabilité thermo-acoustique qui apparaît dans un brûleur académique partiellement prémélangé de méthane/air à pression atmosphérique. La SGE prédit bien l'apparition et le développement des OCL est un bon accord est trouvé entre simulations et expériences en termes d'amplitude et fréquence des OCL. La simulation permet de révéler certains aspects clés responsables du comportement instable de la flamme. Ensuite, une analyse préliminaire de la quantification des incertitudes est fait, ou l'effet des paramètres tels que l'impédance des entrées, le degré de raffinement du maillage ou les pertes thermiques sur les caractéristiques des OCL est évalué. Aussi, la SGE prédit bien la dépendance de la stabilité de la flamme du point d'opération et de la géométrie du brûleur / In the first part of the present PhD. thesis a methodology is presented that allows to predict the soot produced in one-dimensional academic flames, where a semi-empirical soot model is used in combination with a complex chemistry and a detailed radiation solver. The methodology is applied to the computation of soot in a set of ethylene/air counterflow diffusion flames. Several oxidation models are tested and the constants of the model were adjusted to retrieve the experimental results. Also, the effect of radiative losses on soot formation and the flame structure is evaluated. Finally, the performance of the soot model is evaluated on 1D premixed flames, where an alternative expression for the surface growth term is proposed to better reproduce the experimental findings. In the second part of the thesis, Large-Eddy Simulation (LES) and acoustic analysis tools are applied to the prediction of limit cycle oscillations (LCO) of a thermo-acoustic instability appearing in a partially premixed methane/air academic burner operating at atmospheric pressure. The LES captures well the appearance and development of the LCO and a good agreement is found between simulations and experiments in terms of amplitude and frequency of the LCO. Some light is shed on the mechanisms leading to the existence of such instability. Then, a preliminar uncertainty quantification (UQ) analysis is performed, where the effect on the features of the LCO of several computational parameters such as the inlets impedances, mesh refinement or heat losses is assessed. Also, the LES captures well the flame stability behaviour dependence on the operating point and the burner geometry

La Simulation aux Grandes Echelles (SGE) est appliquée à des brûleurs industriels pour prédire de nombreux phénomènes physiques complexes, tel que l’allumage ou la formation de polluants. La prise en compte de réactions chimiques détaillées est alors indispensable pour obtenir des résultats précis. L’amélioration des moyens de calculs permet de réaliser des simulations de brûleurs avec une chimie de plus en plus détaillée. La principale problématique est le couplage entre les réactions chimiques et l’écoulement turbulent. Bien que la dynamique de flamme soit souvent bien reproduite avec les modèles actuels, la prédiction de phénomènes complexes comme la formation de polluants reste une tâche difficile. En particulier, des études ont montré que l’influence du plissement de sous-maille sur la structure chimique des flammes n’était pas prise en compte de manière précise. Deux modèles basés sur le filtrage explicite des fronts de flammes sont étudiés dans cette thèse afin d’améliorer la prédiction de polluants en combustion turbulente prémélangée : (i) le premier modèle met en jeu une méthode de déconvolution des variables filtrées ; (ii) le second modèle implique l’optimisation de la chimie pour obtenir des flammes turbulentes filtrées. L’objectif de la thèse est d’obtenir une prédiction précise des polluants à coût de calcul réduit. / Large Eddy Simulation (LES) is applied to industrial burners to predict a wide range of complex physical phenomena, such as flame ignition and pollutants formation. The prediction accuracy is tightly linked to the ability to describe in detail the chemical reactions and thus the flame chemical structure. With the improvement of computational clusters, the simulation of industrial burners with detailed chemistry becomes possible. A major issue is then to couple detailed chemical mechanisms to turbulent flows. While the flame dynamics is often correctly simulated with stateof- the-art models, the prediction of complex phenomena such as pollutants formation remains a difficult task. Several investigations show that, in many models, the impact of flame subgrid scale wrinkling on the chemical flame structure is not accurately taken into account. Two models based on explicit flame front filtering are explored in this thesis to improve pollutants formation in turbulent premixed combustion: (i) a model based on deconvolution of filtered scalars; (ii) a model involving the optimization of chemistry to reproduce filtered turbulent flames. The objective of the work is to achieve high accuracy in pollutants formation prediction at low computational costs.

Une étude tridimensionnelle du déferlement plongeant des vagues est proposée sur la base de résultats obtenus par Simulation des Grandes Echelles des équations de Navier-Stokes en formulation diphasique. Après une étude bibliographique situant le problème, la première partie du travail est consacrée au développement et à la validation des méthodes numériques employées pour la résolution des équations. Dans une seconde partie de travail, l'outil numérique est exploité afin d'étudier différents problèmes bidimensionnels de propagation d'ondes solitaires. La dernière partie du travail consiste à appliquer les méthodes numériques au déferlement tridimensionnel d'ondes sinusoïdales périodiques de larges amplitudes. Une description complète du phénomène est proposée, incluant la génération du splash-up ainsi que la dissipation de l'énergie par les différentes structures turbulentes présentes dans l'écoulement.

Ces travaux, initiés par le CNES et la SNECMA, sont principalement consacrés à la simulation numérique d'écoulements turbulents supersoniques en présence d'interactions ondes de choc/couche limite. L'étude de ces écoulements, sièges de phénomènes complexes, se révèle être d'une importance particulière en vue du dimensionnement d'organes de propulsion de lanceurs spatiaux. Il a été choisi, dans ce travail, d'investir les aspects instationnaires au moyen de la simulation des grandes échelles turbulentes. Les contraintes liées au traitement numérique ont aboutit au développement d'un schéma numérique hybride, permettant de minimiser les effets de dissipation inhérents aux schémas dédiés à la capture de chocs. Le problème de la génération de conditions aux limites turbulentes réalistes est également abordé. Une technique due à Lund, fondée sur un principe de renormalisation et permettant à l'écoulement de générer lui-même ses conditions d'entrée, a été retenue. Les résultats des simulations se sont révélées être en très bon accord avec les mesures expérimentales ainsi qu'avec les DNS de référence. En particulier, l'existence de basses fréquences associées au mouvement du choc réfléchi a pu être constatée, en accord avec les observations expérimentales. De plus, la simulation a révélé la présence de fréquences similaires à l'intérieur du bulbe de recirculation, venant ainsi conforter l'hypothèse selon laquelle les instabilités du choc de décollement sont conditionnées par celles de la zone décollée.

Contexte du sujet : La climatisation est un poste de dépense énergétique important dans le secteur du bâtiment, qui pourrait être réduit de manière drastique par l'utilisation de systèmes passifs de rafraîchissement. Dans les climats chauds et humides, le rafraîchissement passif des bâtiments est une solution éprouvée, qui s'ordonne autour de quatre principes : minimiser les apports de chaleur interne et externe, apporter de l'inertie au bâtiment, humidifier l'air et assurer une bonne ventilation pour favoriser les échanges convectifs. Objectifs : La description des transferts thermo-convectifs (évaluation des débits massiques, des transferts de chaleur) mis en jeu dans des cavités ouvertes (pièces avec ventilation traversante, cheminées solaires, doubles peaux, espaces sous-toiture) est encore mal connue et constitue un enjeu pour l'amélioration des systèmes passifs. L'étude de ces phénomènes peut être évaluée par la mécanique des fluides numérique. Les objectifs de ce travail de thèse sont : de réaliser des simulations numériques fines du comportement thermo-aéraulique dans des configurations typiques de systèmes passifs en climat tropical humide, afin d'améliorer nos connaissances sur la ventilation naturelle, d'approfondir et d'apporter des éléments de réponses en ce qui concerne le choix des conditions limites numériques à appliquer pour les systèmes ouverts. Modélisation numérique : L'approche numérique adoptée dans ce travail, pour étudier la convection naturelle turbulente, est la simulation des grandes échelles (SGE ou LES en anglais). Cette approche se situe à mi-parcours entre la méthode de calcul direct et la réso- lution des équations moyennées en temps, de type RANS. L'avantage d'une telle technique est la réduction appréciable du nombre de points de discrétisation nécessaire par rapport à ce- lui exigé par la méthode de calcul direct, tout en conservant l'aspect dynamique des écoulements. Résultats : Les résultats obtenus lors de ce travail concernent l'étude des conditions limites dynamiques à imposer pour des géométries ouvertes avec une SND et l'application de la SGE à différentes configurations de cavités ouvertes en régime turbulent, afin de caractériser les champs de température et de vitesse et d'en déduire les grandeurs intégrales d'intérêt (débit massique, débit enthalpique, renouvellement d'air, ...). Les résultats de ces calculs ont été comparés soit à d'autres résultats numériques dans le cadre de benchmarks nationaux (benchmark numérique AmeTh et ADNBâti) ou à des résultats expérimentaux.