Les transferts radiatifs jouent un rôle important dans les chambres de combustion des installations industrielles. En effet, il existe un couplage fort entre la combustion turbulente et le rayonnement. Dans le but d'étudier ce couplage, le code Rainier est développé pour les calculs de pertes par rayonnement dans un écoulement réactif dans des géométries complexes. Ce code repose sur des simulations aux grandes échelles (LES) de la combustion turbulente. Il est basé sur les maillages tétraédriques non structurés. Le modèle de rayonnement appliqué à la modélisation des propriétés radiatives des gaz est le modèle CK (Correlated-k). La méthode statistique de Monte-Carlo (ERM) est utilisée pour résoudre l'équation de Transfert du Rayonnement (ETR). Le code de rayonnement est parallélisé et il montre une réponse linéaire en fonction du nombre de processeurs très proche de la réponse idéale. Une méthode de couplage de code de rayonnement avec le code de combustion LES est développée. Chacun des codes a sa propre logique d'architecture et de développement. En conséquence, le couplage entre les deux domaines d'étude est réalisé de telle façon que les échanges des données et les synchronisations entre eux soient assurés. Les résultats obtenus à partir du couplage des sur une chambre de combustion d'hélicoptère sont présentés. Nous avons montré que le rayonnement modifie les champs instantanés de température et d'espèces à l'intérieur de la chambre de combustion.

Les transferts radiatifs jouent un rôle important dans les chambres de combustion des installations industrielles. En effet, il existe un couplage fort entre la combustion turbulente et le rayonnement. Dans le but d’étudier ce couplage, le code Rainier est développé pour les calculs de pertes par rayonnement dans un écoulement réactif dans des géométries complexes. Ce code repose sur des simulations aux grandes échelles (LES) de la combustion turbulente. Il est basé sur les maillages tétraédriques non structurés. Le modèle de rayonnement appliqué à la modélisation des propriétés radiatives des gaz est le modèle CK (Correlated-k). La méthode statistique de Monte-Carlo (ERM) est utilisée pour résoudre l’équation de Transfert du Rayonnement (ETR). Le code de rayonnement est parallélisé et il montre une réponse linéaire en fonction du nombre de processeurs très proche de la réponse idéale. Une méthode de couplage de code de rayonnement avec le code de combustion LES est développée. Chacun des codes a sa propre logique d’architecture et de développement. En conséquence, le couplage entre les deux domaines d’étude est réalisé de telle façon que les échanges des données et les synchronisations entre eux soient assurés. Les résultats obtenus à partir du couplage des sur une chambre de combustion d’hélicoptère sont présentés. Nous avons montré que le rayonnement modifie les champs instantanés de température et d’espèces à l’intérieur de la chambre de combustion. / Radiation plays an important role in industrial combustion chambers. In fact, there is a strong coupling between combustion and heat transfer in these turbulence chambers. The Rainier code is developed for the calculation of radiation in a reactive flow within a complex geometry. This code is dedicated to Large Eddy Simulation (LES) of turbulent reactive flows. The code is based on unstructured tetrahedral mesh. The correlated k-distribution method (CK) is applied to the modelling of radiative properties of gases. A statistical method of Monte Carlo (ERM) is used to resolve the Radiation Transfer Equation (RTE). The code is parallelized and it shows a linear response to the number of processors, which is very close to the ideal response. A coupling method of the Rainier code and the turbulent combustion code LES is implemented. Each code has its own logic architecture and development. For this, the coupling between these two different fields of study is achieved in such a way that the data exchange and synchronization between them is assured. The results obtained by applying the coupled codes to the combustion chamber of helicopter are presented. We have shown that the radiative transfer modifies the temperature and species fields inside the combustion chamber.

Lors des études de désenfumage, le recours à la mécanique des fluides numérique, par l’emploi de codes tridimensionnels (CFD), est de plus en plus employé. Cependant, suite aux études de Dalmarnock (2007) et d’Arcueil (2009), il a été montré que les résultats numériques peuvent s’éloigner de façon significative des résultats expérimentaux. Ces différences peuvent être notamment induites par des faibles erreurs de certains paramètres clés, ce qui limite l’utilisation prédictive de ces codes dans les études d’ingénierie de la sécurité incendie. Dans ce cadre, cette thèse est consacrée à la simulation numérique d’écoulement induit par une source de chaleur. Les simulations sont réalisées avec le code FDS, largement utilisé par la communauté incendie, afin d’étudier sa capacité à reproduire la dynamique de ces écoulements. L’étude se porte principalement sur la simulation d’un panache, créé par une plaque chauffante, avec une comparaison détaillée à des résultats expérimentaux. Les résultats des simulations montrent que, les profils radiaux des statistiques en un point de vitesse et température d’ordre 1 et 2, sont en très bon accord avec les résultats de la littérature. Les évolutions des flux de masse, de flottabilité et de quantité de mouvement, ainsi que celle du nombre de Richardson du panache, sont également bien reproduites. Après avoir validé la simulation, l’évolution du coefficient d’entraînement avec la hauteur a pu être examinée en détail, montrant que ce coefficient diminue en s’éloignant de la source de chaleur, passant de 0.4 en partie proche de la source à 0.1 en partie éloignée. Cette évolution a pu être expliquée par une dépendance du coefficient d’entraînement avec le nombre de Richardson local du panache et le rapport des masses volumiques entre le panache et l’air ambiant. La thèse a également permis de montrer que les variables étudiées sont sensibles à différents modèles de sous maille. Pour finir, une méthodologie a été proposée pour étudier la sensibilité de différents paramètres d’entrée de simulation et appliquée sur un cas réaliste d’incendie dans un bâtiment. Cette analyse a permis de mettre en évidence plusieurs paramètres explicites (paramètres définis en concertation avec le propriétaire du bâtiment, les services d’incendie et les autorités compétentes) qui modifient les résultats de la simulation : la puissance du foyer, le débit d’extraction des gaz chauds, mais également certains paramètres implicites (paramètres définis par le modélisateur) comme le débit d’insufflation au niveau des entrées d’air frais. / Fire Safety Engineering studies often rely on the use of Computational Fluid Dynamics, but several recent studies (e.g. Dalmarnock 2007, Arcueil Numérique 2009) have highlighted significant differences between experimental data and the results of numerical simulations. The numerical models are often rather sensitive to the choice of boundary conditions, and may depend on key parameters which are poorly understood. The research presented in this thesis addresses some of these questions, using the CFD code FDS to simulate a set of experimental test cases, and to test the sensitivity of the simulations to different input parameters. The general configuration studied is that of a thermal plume generated by a heated plate, either rising freely or else impacting on a ceiling. Radial profiles of 1st and 2nd order statistics of temperature and velocity agree well with published data, as does the axial variation of Mass, Momentum and Buoyancy fluxes, and the local Richardson number. A detailed study of the entrainment coefficient shows that this varies with distance from the heat source, between about 0.4 and 0.1; this is explained by the dependence of entrainment on the local Richardson number and the density difference between the plume and the ambient air. The comparison between simulations has also been used to investigate the applicability of different sub-grid scale models for this type of flow. Finally, a methodology has been developed for evaluating the sensitivity of the calculations to different input parameters, and the method has been applied to a test configuration representing a fire in a building. The results show the relative importance of both explicitly-defined parameters that characterize the situation (thermal power of the fire, ventilation regime inside the building) and implicit model parameters.

La prédiction des flux aux parois joue un rôle déterminant dans le cycle de vie des chambres de combustion. Le transfert de chaleur de la flamme aux parois est entraîné, outre la convection, également par le rayonnement des gaz chauds au sein de la chambre. Afin d’intégrer les contributions convectives et radiatives au flux pariétal il est nécessaire de résoudre simultanément l’équation de transfert radiatif et les équations régissant l’écoulement réactif. Quand les méthodes de Monte Carlo sont couplées aux simulations aux grandes échelles (LES), de telles simulations deviennent très coûteuses. L’objectif de cette thèse est donc d’investiguer une technique pour améliorer l’efficacité de la méthode MC, basée sur un mécanisme alternatif d’échantillonnage appelée intégration Quasi-Monte Carlo (QMC). Au cours de cette thèse, la méthode QMC a été couplée à une simulation LES dans une configuration où le rayonnement joue un rôle très important : la flamme méthane-air de la chambre Oxytec. La comparaison entre les simulations couplées et non couplées avec les données expérimentales montre que le rayonnement thermique a un impact sur la topologie de l’écoulement et de la flamme. Enfin, un bon accord est trouvé entre le flux de chaleur pariétal prédit par la simulation et les données expérimentales. / The prediction of wall fluxes is a significant aspect in the life cycle of combustors, since it allows to prevent eventual wall damages. Heat transfer from flame to the walls is driven, apart from convection, also by radiation of burnt gases inside the chamber. In order to correctly account for both convective and radiative contributions to wall fluxes, the simultaneous solution of the radiative transfer equation (RTE) and the governing equations for reactive flows is required. However, multi-physics simulations where MC methods are coupled to Large Eddy Simulation (LES), remain very costly. The purpose of this study is then to investigate improvements of MC methods, by using an alternative sampling mechanism for numerical integration usually referred to as Quasi-Monte Carlo (QMC) integration. In this study, QMC method is coupled to Large Eddy Simulation (LES) of a configuration where the radiation plays an important role: the methane-air flame investigated during the experimental campaign Oxytec. Coupled and non-coupled simulations are compared and their comparison with experimental data shows that thermal radiation has an impact on both flow and flame topology. Finally a good agreement is found between numerical wall fluxes and experimental conductive fluxes.

Nous cherchons à modéliser numériquement un jet flottant air-hélium dans une cavité avec deux ouvertures à partir de simulations aux grandes échelles (LES) et de simulations numériques directes (DNS). La configuration considérée est basée sur une étude expérimentale menée au CEA de Saclay reproduisant une fuite d'hydrogène en environnement confiné. La dimension de la cavité a été choisie pour permettre une transition laminaire-turbulent intervenant environ à la mi-hauteur de la cavité. Cette étude porte principalement sur trois points majeurs : l'influence des conditions aux limites sur le développement du jet et son interaction avec l'environnement extérieur, la validité du modèle numérique qui est analysée en comparant la distribution de vitesse obtenue numériquement aux mesures expérimentales (PIV) et, enfin, la compréhension de la distribution air-hélium et le phénomène de stratification qui s'établit à l'intérieur de la cavité. Nous observons dans un premier temps que des conditions limites de pression constante appliquées directement au ras des évents conduisent à une sous-estimation du débit volumique d'air entrant dans la cavité et donc à une surestimation de la masse de l'hélium à l'intérieur de la cavité, ce qui n'est pas acceptable dans un contexte d'évaluation du risque hydrogène. En revanche, la prise en compte, dans le domaine de calcul, d'une région extérieure à la cavité prédit correctement le flux d'air entrant. Les résultats numériques sont alors en bon accord avec les données PIV. Il a été montré que les prédictions de la DNS, par rapport à la LES, concordent mieux avec les mesures de vitesse par PIV. Le champ de concentration prédit numériquement présente une couche homogène en haut de la cavité, dont la concentration est en accord avec le modèle théorique de Linden et al. 1990. Cependant, sa position et son épaisseur ne correspondent pas au modèle. Ceci est principalement dû aux interactions directes entre le jet flottant et, d'une part, avec les limites solides de la cavité et d'autre part, avec l'environnement extérieur. L'analyse statistique concernant la production de la flottabilité de l'énergie cinétique turbulente (TKE) a permis d'identifier les limites du jet flottant. / We present numerical results from large eddy simulations (LES) and coarse direct numerical simulations (DNS) of an air-helium buoyant jet rising in a two vented cavity. The geometrical configuration mimics the helium release experimental set-up studied at CEA Saclay in the framework of security assessment of hydrogen-based systems with an indoor usage. The dimension of the enclosure was chosen to ensure a laminar-turbulent transition occurring at about the middle height of the cavity. This study focuses mainly on three key points : the influence of the boundary conditions on the jet development and its interaction with the exterior environment, the validity of the numerical model which is analyzed by comparing the numerical velocity distribution versus the measured particle image velocimetry (PIV) ones, and finally understanding the distribution of air-helium and the stratification phenomenon that takes place inside the cavity. We observe at first that applying constant pressure outlet boundary conditions directly at the vent surfaces underestimates the volumetric flow rate of air entering the enclosure and thus overestimate the helium mass inside the cavity. On the contrary, modelling an exterior region in the computational domain better predicts the air flow-rate entrance and numerical results matches better with the experimental PIV data. It has been figured out that the coarse DNS predictions match better with the velocities PIV measurements, compared to the LES. Numerical prediction of the helium field depicts a homogeneous layer formed at the top of the cavity, with a concentration in good agreement with the theoretical model of Linden et al. 1990. However, the position and the thickness of the layer do not correspond to the theory. This is mainly due to the direct interactions between the buoyant jet and both the solid boundaries of the cavity and the exterior environment. Statistical analysis regarding the buoyancy production of the turbulent kinetic energy (TKE) served to identify the limits of the buoyant jet.

L'impact environnemental du trafic aérien fait maintenant l'objet d'une réglementation qui tend à se sévériser. Dans ce contexte, les industriels misent sur l'amélioration des technologies afin de réduire la consommation de carburant et l'émission de polluants. Ces phénomènes dépendent en grande partie des chemins réactionnels sous-jacents, qui peuvent s'avérer très complexes. La Simulation aux Grandes Échelles (SGE) est un outil intéressant afin d'étudier ces phénomènes pour un coût de calcul qui reste raisonnable. Cependant, les processus chimiques, s'ils sont considérés sans simplification, font intervenir des centaines d'espèces aux temps caractéristiques très différents au sein de processus non-linéaires qui induisent une forte raideur dans le système d'équations, et un coût de calcul prohibitif. Permettant de s'absoudre de ces problèmes tout en conservant une bonne capacité de prédiction des polluants, les Chimies Analytiquement Réduites (CAR) font l'objet d'une attention grandissante au sein de la communauté. Les CAR permettent de conserver la physique du problème considéré, en conservant les espèces et voies réactionnelles les plus importantes. Grâce à l'évolution toujours croissante des moyens de calculs, les CAR sont appliqués dans des configurations de plus en plus complexes. Les travaux de thèse ont principalement portés sur deux sujets. Premièrement, une étude poussée des techniques et outils permettant une réduction efficace et systématique de chimies détaillées. L'outil de réduction multiétapes YARC est retenu et exhaustivement employé dans la dérivation et la validation d'une série de CAR préservant la description de la structure de flamme. Ensuite, une investigation de la faisabilité et des bénéfices qu'apportent l'utilisation de CAR en LES, comparé a des approches plus classiques, sur des cas tests de complexité croissante. La première configuration étudiée est une chambre de combustion partiellement pré-mélangée brûlant de l'éthylène, étudiée expérimentalement au DLR. Différentes modélisations de la chimie sont considérées, dont un CAR développé spécifiquement pour ce cas test, et les résultats démontrent qu'une prise en compte des interactions flamme-écoulement est cruciale pour une prédiction juste de la structure de la flamme et des niveaux de suies. La seconde configuration est un brûleur diphasique, avec une injection directe pauvre, brûlant du Jet-A2. Dans cette étude, une approche novatrice pour la prise en compte de la complexité du fuel réel (HyChem) est considérée, permettant la dérivation d’un CAR. Les résultats sont excellents et valident la méthodologie tout en fournissant une analyse précieuse des interactions flamme-spray et de la formation de polluants (NOx) dans des flammes à la structure complexe. / Recent implementation of emission control regulations has resulted in a considerable demand from industry to improve the efficiency while minimizing the consumption and pollutant emissions of the next generation of aero-engine combustors. Those phenomena are shown to strongly depend upon the underlying complex chemical pathways and their interaction with turbulence. Large Eddy Simulation (LES) is an attractive tool to address those issues with high accuracy at a reasonable computing cost. However, the computation of accurate combustion chemistry remains a challenge. Indeed, combustion proceeds through complex and highly non-linear processes that involve up to hundreds of different chemical compounds, which significantly increases the computational time and often induces stiffness in the resolved equations. As a mean to circumvent these drawbacks while retaining the necessary kinetics for the prediction of pollutants, Analytically Reduced Chemistry (ARC) has recently received high interest in the Computational Fluid Dynamics (CFD) community. ARC is a strategy for the description of combustion chemistry where only the most important species and reactions are retained, in a "physically-oriented way". ARC is on the verge of becoming affordable at a design stage, thanks to the continuously increasing available computational resources. The goal of the present work is twofold. A first objective is to test and validate efficient techniques and tools by which detailed chemistries are reduced to an LES-compliant format. To do so, the multi-step reduction tool YARC is selected and employed to derive and validate a series of ARC specifically designed to retrieve correct flame structures. A second objective is to investigate the overall feasibility and benefits of using ARC, combined to the Thickened Flame model (DTFLES), in performing LES of configurations of increasing complexity. The first configuration is a sooting swirl-stabilized non-premixed aero-engine combustor experimentally studied at DLR, burning ethylene. LES of this configuration is performed with the AVBP solver, in which ARC has been implemented. By comparison with global chemistry and tabulated chemistry, results highlight the importance of accurately capturing the flow-flame interactions for a good prediction of pollutants and soot. The second configuration is a swirled twophase flow burner featuring a lean direct injection system and burning Jet-A2. A novel methodology to real fuel modeling (HyChem approach) is employed, which allows subsequent ARC derivation. The excellent results in comparison with measurements constitute an additional validation of the methodology, and provide valuable qualitative and quantitative insights on the flame-spray interactions and on the pollutant formation (NOx) mechanisms in complex flame configurations.

La régulation des polluants a mené à la création de nouveaux systèmes de combustion. Le carburant étant stocké sous forme liquide, sa transformation jusqu’à sa combustion est complexe. La capacité de la Simulation aux grandes échelles à simuler des écoulements turbulents réactifs a été montrée sur des cas académiques comme sur des configurations industrielles, tout en prenant en compte les phénomènes multiphysiques intervenant dans ces configurations, mais les études sur la structure de flamme diphasique sont encore trop peu nombreuses. La présence de deux solveurs pour la simulation d’une phase liquide étant disponible dans le code AVBP, leur utilisation permet une comparaison et une compréhension des phénomènes en jeu combinant dispersion, évaporation, et combustion. La première partie de l’étude relate la validation du modèle d’injection FIM-UR. Ce modèle est capable de reconstruire les profils de vitesses et de granulométrie à l’injecteur sans avoir à simuler les phénomènes d’atomisation primaire et secondaire. Une validation en régime turbulent avait déjà été réalisée, et on propose ici de valider le modèle dans un cas laminaire. Des comparaisons entre simulations monodisperses et polydisperse et des expériences sont effectuées. La simulation monodisperse Lagrangienne donne une bonne structure globale mais la simulation polydisperse Lagrangienne permet de retrouver le comportement au centre du cône avec la présence des petites gouttes et à la périphérie du cône par la présence des grosses gouttes. De plus, des améliorations sont apportées au modèle pour le formalisme Eulérien et montrent de bons résultats. La partie suivante s’intéresse à caractériser un spray polydisperse par une distribution monodisperse. En effet, au cas où une approche polydisperse n’est pas possible, le choix du diamètre moyen à prendre pour une simulation monodisperse est délicat. On propose donc d’analyser le comportement d’un spray polydisperse en le comparant à ceux de sprays monodisperses. Deux configurations académiques sont choisies : des cas de Turbulence Homogène Isotrope chargée en particules pour étudier la dynamique, et des calculs d’évaporation 0D. Trois paramètres sont étudiés pour la dynamique : la concentration préférentielle (ou ségrégation), la traînée moyenne et la traînée réduite moyenne. Cette dernière et la ségrégation de la distribution polydisperse semblent affectées par les tailles de goutte les plus faibles, et la concentration préférentielle apparait alors comme la moyenne des ségrégations des classes qui la composent pondérées par l’inverse du nombre de Stokes associé à chacune de ces classes. La traînée moyenne de la simulation polydisperse possède un comportement proche des diamètres moyens D10 et D20. Ces analyses nous poussent donc à choisir le D10 pour caractériser la dynamique d’un spray polydisperse. Les calculs d’évaporation 0D ne permettent pas dans un premier temps de caractériser efficacement la masse évaporée d’un spray polydisperse par celle d’un spray monodisperse équivalent, mais la définition de nouveaux diamètres issus de la littérature des lits fluidisés comme le D50% le permet, ce qui le place autour du D32. On propose donc de caractériser l’évaporation d’un spray polydisperse par ce diamètre. Enfin, la dernière partie étudie la structure de flamme diphasique dans la chambre MERCATO, à l’aide du formalisme Lagrangien, monodisperse et polydisperse, mais aussi en utilisant le formalisme Eulérien. La validation du modèle FIM-UR du premier chapitre et ses améliorations sont utilisées pour représenter les conditions d’injection liquide. En plus d’un calcul polydisperse, deux simulations monodisperses Lagrangiennes sont réalisées en prenant les diamètres moyens D10 et D32, suite à la partie précédente. Des comparaisons qualitatives et des validations sont réalisées, en comparant des profils de vitesses gazeuses axiale et fluctuante et vitesse axiale liquide issus de l’expérience. / Regulations on pollutants have led to the creation of new combustion systems. Giving that fuel is stored in a liquid form, its evolution until combustion is complex. The ability of Large Eddy Simulation has been demonstrated on academic cases, as well as on industrial configurations, by taking into account the multi-physics phenomena, but there is a lack of studies about two-phase flow flame structures. Two solvers for the simulation of two-phase flows are available in the AVBP code, hence both simulations are performed to compare and increase understanding of the phenomena involved such as dispersion, evaporation and combustion. The first part of the study focuses on the validation of the FIM-UR injection model. This model is able to build velocity and droplet profiles at the injector, without simulating primary and secondary break up. A validation in a turbulent case has already been done, and this study validates the model in a laminar case. Comparisons between monodisperse and polydisperse simulations, and experiments are performed. The monodisperse Lagrangian simulation shows good results but the polydisperse simulation is able to represent profiles in the center of the cone by small droplets and at the peripheral part of the cone, by big ones. Moreover, improvements in the Eulerian model exhibit good results. The next section tries to evaluate the impact of polydispersion. Indeed, when a polydisperse approach is not available, choosing the mean diameter can be tricky. A comparison between the behavior of polydisperse spray and monodisperse sprays ones is realised. Two academic cases are studied: Homogeneous Isotropic Turbulence with particles to analyze the dynamics, and 0D evaporation cases. For the dynamics, preferential concentration, mean drag and reduced mean drag are studied. The latter and preferential concentration are affected by small droplets, and the preferential concentration of a polydisperse spray is equivalent to the average of preferential concentration of classes, extracted from the polydisperse distribution, weighted by the inverse of the Stokes number of each class. The mean drag behaves like the D10 and D20 mean drags. This analysis allows us to choose the D10 to characterize a polydisperse distribution for the dynamics. Zero-D evaporation simulations cannot characterize the polydisperse spray evaporated mass by the evaporated mass of monodisperses sprays. New definitions of diameters from fluidized bed literature enable the use of D50%, which is close to D32. We propose to use this diameter to characterize the evaporation of a polydisperse spray. Finally, the last section studies the structure of two-phase flames in the MERCATO bench, using the Lagrangian formalism, monodisperse and polydisperse but also using the Eulerian formalism. The validation of FIM-UR model and improvements from the first section are used to represent liquid injection conditions. A polydisperse simulation is realized and two monodisperse simulations are computed using mean diameters D10 and D32, thanks to the previous section. Qualitative comparisons and validations are realized, comparing gaseous velocity profiles and liquid velocity profiles. Good agreements are found and the mean diameter D32 seems to be close to the polydisperse spray. A comparison between mean flames is done with an Abel transform of the flame from the experiments. The flame has an "M shape", anchored by small recirculation zones out of the swirler, and by a point at the tip of the central recirculation zone. Then, the impact of droplet distributions is analyzed. Even if few bigger droplets from the polydisperse distribution are convected in the hot gases due to bigger particular time and evaporation time, two-phase flow flame structures are equivalent. Different combustion regimes appeared with premixed flames and pockets of fuel burning in the hot gases.

La simulation de la combustion turbulente connait un nouvel essor avec l'introduction de la Simulation aux Grandes Échelles (SGE) qui permet de prédire l'évolution in stationnaire de l'écoulement réactif turbulent. Dans ce contexte la prise en compte du rayonnement soulève des questions d'ordre a la fois fondamental et pratique. En effet les processus physiques du rayonnement et de la combustion sont de nature radicalement différente : la combustion est contrôlée par des échanges locaux sur une durée finie, alors que le rayonnement est instantané et fait intervenir des échanges a distance. En premier lieu il convient de s'interroger sur l'impact de la modélisation SGE de la combustion turbulente sur le rayonnement. Cette question est traitée dans le cadre plus général de l'interaction rayonnement-turbulence. A partir d'études théoriques et numériques, il est montre que cette interaction est faible et qu'une solution SGE peut être directement utilisée pour un calcul radiatif, sans modélisation supplémentaire. Il s'agit ensuite de mettre en place de façon pratique le couplage in stationnaire rayonnement-combustion turbulente. Un point clé est la réduction du temps de calcul pour le rayonnement, et diverses stratégies sont proposées. En particulier un nouveau modèle spectral est introduit, utilisant une technique de tabulation et garantissant un niveau de précision suffisant. Le temps de calcul radiatif a ainsi été réduit de deux ordres de grandeur, permettant la réalisation d'un calcul couple sur une configuration de flamme pré-melangée turbulente. / Simulation of turbulent combustion has gained high potential with the Large Eddy Simulation (LES) approach, allowing to predict unsteady turbulent reactive flows. In this context, taking into account radiation rises new fundamental and practical questions. Indeed the physics involved in radiation and in combustion are completely different : combustion is controlled by local exchanges and finite times whereas radiation is instantaneous and is based on non-local exchanges. In a first step, the impact of LES modelling of turbulent combustion on radiation is regarded. This question is treated in the more general frame of the turbulence-radiation interaction. From theoretical and numerical studies, it is shown that this interaction is weak in the LES context so that LES solutions can be directly coupled to radiative calculations, without further modelling. Then the unsteady coupling of radiation and turbulent combustion is realised. A key point is the reduction of calculation time of radiation, and several strategies are proposed. In particular a new global spectral model is introduced, using a tabulation technique and ensuring a sufficient level of accuracy. The radiative time calculation is finally decreased by two orders of magnitude, enabling the realization of a coupled calculation of a turbulent premixed flame

L’étude présentée dans cette thèse relève de la mécanique des fluides expérimentale et numérique appliquée aux écoulements pariétaux de refroidissement de chambres de combustion aéronautiques. En présence de phénomènes thermo-acoustiques, comme les instabilités de combustion, il est important d’évaluer si les capacités de l’écoulement pariétal à protéger les parois de chambre restent suffisantes. C’est ainsi que nous nous sommes intéressés aux écoulements de paroi multiperforée soumis à une excitation acoustique. Dans ce but, le banc d’essais MAVERIC a été amélioré grâce à l’installation d’un système qui permet de forcer acoustiquement l’écoulement transverse dans lequel les jets pariétaux débouchent. Nous avons pu alors mettre en évidence la forte sensibilité de ce type d’écoulements à l’excitation acoustique. Le bon accord entre les résultats expérimentaux et les simulations numériques aux grandes échelles (LES) effectuées est très encourageant dans le cas d’un forçage par onde stationnaire. Le forçage par onde progressive, étudié uniquement par simulations numériques, s’est révélé être capable de modifier significativement la topologie de l’écoulement. Enfin, à partir de l’outil numérique AVBP-AVTP qui permet le couplage de calculs fluide-solide, nous avons réalisé une étude de l’influence de la présence d’une excitation acoustique sur le comportement thermique de l’écoulement autour d’une paroi multiperforée de chambre de combustion. / This experimental and numerical study in the field of fluid mechanics deals with jets-in cross flow configurations that are relevant for the cooling of aero engine combustion chambers. Indeed, in presence of instabilities it is important to determine to which extent the film cooling is able to do its job of preserving the combustion chamber walls from the thermal load. The test facility MAVERIC has been upgraded in order to acoustically force the crossflow in which the jets are discharging. The strong sensitivity of the overall flow unsteady properties to the presence of the acoustic forcing has been clearly evidenced. The agreement between the experimental results and large-eddy simulations proved to be quite encouraging for a stationary acoustic wave whereas the case of a propagating acoustic wave investigated only numerically reveals also quite a significant change of the flow topology. In this context, the effect of the acoustic forcing on the wall thermal behavior has been analyzed thanks to the use of the fluid-solid coupled AVBP-AVTP solver.