<div>Studies of turbulent combustion in the past mainly focus on problems with single-regime combustion. In practical combustion systems, however, combustion rarely occurs in a single regime, and different regimes of combustion can be observed in the same system. This creates a significant gap between our existing knowledge of combustion in single regime and the practical need in multi-regime combustion. In this work, we aim to extend the traditional single-regime combustion models to problems involving different regimes of combustion. Among the existing modeling methods, Transported Probability Density Function (PDF) method is attractive for its intrinsic closure of treating detailed chemical kinetics and has been demonstrated to be promising in predicting low-probability but practically important combustion events like local extinction and re-ignition. In this work, we focus on the model assessment and advancement of the Large Eddy Simulation (LES)/ PDF method in predicting turbulent multi-regime combustion.</div><div><br></div><div><div>Two combustion benchmark problems are considered for the model assessment. One is a recently designed turbulent piloted jet flame that features statistically transient processes, the Sydney turbulent pulsed piloted jet flame. A direct comparison of the predicted and measured time series of the axial velocity demonstrates a satisfactory prediction of the flow and turbulence fields of the pulsed jet flame by the employed LES/PDF modeling method. A comparison of the PLIF-OH images and the predicted OH mass fraction contours at a few selected times shows that the method captures the different combustion stages including healthy burning, significant extinction, and the re-establishment of healthy burning, in the statistically transient process. The temporal history of the conditional PDF of OH mass fraction/temperature at around stoichiometric conditions at different axial locations suggests that the method predicts the extinction and re-establishment timings accurately at upstream locations but less accurately at downstream locations with a delay of burning reestablishment. The other test case is a unified series of existing turbulent piloted flames. To facilitate model assessment across different combustion regimes, we develop a model validation framework by unifying several existing pilot stabilized turbulent jet flames in different combustion regimes. The characteristic similarity and difference of the employed piloted flames are examined, including the Sydney piloted flames L, B, and M, the Sandia piloted flames D, E, and F, a series of piloted premixed Bunsen flames, and the Sydney/Sandia inhomogeneous inlet piloted jet flames. Proper parameterization and a regime diagram are introduced to characterize the pilot stabilized flames covering non-premixed, partially premixed, and premixed flames. A preliminary model assessment is carried out to examine the simultaneous model performance of the LES/PDF method for the piloted jet flames across different combustion regimes.</div><div><br></div><div>With the assessment work in the above two test cases, it is found that the LES/PDF method can predict the statistically transient combustion and multi-regime combustion reasonably well but some modeling limitations are also identified. Thus, further model advancement is needed for the LES/PDF method. In this work, we focus on two model advancement studies related to the molecular diffusion and sub-filter scale mixing processes in turbulent combustion. The first study is to deal with differential molecular diffusion (DMD) among different species. The importance of theDMD effects on combustion has been found in many applications. However, in most previous combustion models equal molecular diffusivity is assumed. To incorporate the DMD effects accurately, we develop a model called Variance Consistent Mean Shift (VCMS) model. The second model advancement focuses on the sub-filter scale mixing in high-Karlovitz (Ka) number turbulent combustion. We analyze the DNS data of a Sandia high-Ka premixed jet flame to gain insights into the modeling of sub-filter scale mixing. A sub-filter scale mixing time scale is analyzed with respect to the filter size to examine the validity of a power-law scaling model for the mixing time scale.</div></div>

[pt] O presente trabalho apresenta um estudo numérico da transferência de energia por radiação e sua interação com as propriedades radiantes cinzas e espectrais dos gases produtos da combustão e da fuligem que são formados em um processo de combustão turbulenta. Assim, utilizam-se sistemas de forno/- queimador que operam em regime de chama não pré-misturada de maneira a avaliar, através da dinâmica dos fluidos computacional (CFD), a influência que exercem os diversos modelos de propriedades radiantes sobre a representação da termoquímica do escoamento reativo. Com o objetivo de identificar as principais características e deficiências que apresentam cada um destes modelos, foram considerados dois cenários. O primeiro, correspondente a um problema de radiação unidimensional de um sistema homogêneo e não isotérmico onde são estudados, modelos de propriedades radiantes (i) disponíveis em um software comercial de CFD e (ii) aqueles que foram implementados neste trabalho. Além disso, foi empregado um código numérico que determina as propriedades radiantes espectrais de gases produtos da combustão e da fuligem através de uma abordagem de banda estreita. Para este fim, este código foi acoplado com o software de CFD. Em seguida, dois queimadores de porte laboratorial são empregados de forma a avaliar a capacidade preditiva dos modelos de propriedades radiantes: o primeiro queima propano gasoso e ar enriquecido com oxigênio e o segundo utiliza querosene líquido e oxigênio como reagentes. Dados experimentais de fluxo de calor radiante e de fração volumétrica da fuligem são utilizados para comparação com os resultados obtidos da simulação. Para ambas as configurações de queimador foi também estudado o modelo de Moss-Brookes para previsão da formação/consumo da fuligem. Os resultados obtidos demonstraram o bom desempenho da maioria dos modelos de propriedades radiantes estudados. Em particular, a abordagem de banda estreita foi o que melhor previu a radiação térmica. Além disso, a sua utilização com o modelo de Moss-Brookes levou à melhor previsão da fração volumétrica da fuligem. / [en] This work presents a numerical study of radiation heat transfer and its interaction with gray and spectral radiation of combustion products and soot that are formed in a turbulent combustion process. Different burner/furnace systems operating in a non-premixed combustion regime were used in order to evaluate, through computational fluid dynamics (CFD), the influence of several radiant properties models. Aiming to identify the key features and shortcomings that exhibit each of these models, two scenarios were considered. The first corresponds to a 1-D radiation problem where radiative properties models of a homogeneous non isothermal system are studied as (i) available CFD commercial software and (ii) those implemented in this work. Moreover, a numerical code was used in order to determine, through a narrow band approach, the spectral radiative properties of soot and combustion products. For this purpose, this code was coupled with the CFD software. Then, two laboratory-scale burners are used to assess the predictive capacity of radiative properties models: the first, burning propane and enriched air oxygen, and the second uses kerosene and oxygen as reactants. Measurements of radiant heat flux and soot volumetric fraction are used for comparison with simulation results. For both configurations, the performance of the Moss-Brookes model for predicting the soot production was also studied. The results of this study demonstrated the good performance of the majority of the radiant properties models studied. Particularly, the narrow band approach was the model that provided the best thermal radiation prediction. Moreover, the combination of the narrow band approach with the Moss-Brookes model lead to the best prediction of soot volume fraction.

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.

On s’intéresse ici aux erreurs de modélisation liées à l’usage de modèles de flammelette sous-maille en combustion turbulente non prémélangée. Le but de cette thèse est de développer une stratégie d’estimation d’erreur a posteriori pour déterminer le meilleur modèle parmi une hiérarchie, à un coût numérique similaire à l’utilisation de ces mêmes modèles. Dans un premier temps, une stratégie faisant appel à un estimateur basé sur les résidus pondérés est développée et testée sur un système d’équations d’advection-diffusion-réaction. Dans un deuxième temps, on teste la méthodologie d’estimation d’erreur sur un autre système d’équations, où des effets d’extinction et de réallumage sont ajoutés. Lorsqu’il n’y a pas d’advection, une analyse asymptotique rigoureuse montre l’existence de plusieurs régimes de combustion déjà observés dans les simulations numériques. Nous obtenons une approximation des paramètres de réallumage et d’extinction avec la courbe en «S», un graphe de la température maximale de la flamme en fonction du nombre de Damköhler, composée de trois branches et d’une double courbure. En ajoutant des effets advectifs, on obtient également une courbe en «S» correspondant aux régimes de combustion déjà identifiés. Nous comparons les erreurs de modélisation liées aux approximations asymptotiques dans les deux régimes stables et établissons une nouvelle hiérarchie des modèles en fonction du régime de combustion. Ces erreurs sont comparées aux estimations données par la stratégie d’estimation d’erreur. Si un seul régime stable de combustion existe, l’estimateur d’erreur l’identifie correctement ; si plus d’un régime est possible, on obtient une fac˛on systématique de choisir un régime. Pour les régimes où plus d’un modèle est approprié, la hiérarchie prédite par l’estimateur est correcte. / We are interested here in the modeling errors of subgrid flamelet models in nonpremixed turbulent combustion. The goal of this thesis is to develop an a posteriori error estimation strategy to determine the best model within a hierarchy, with a numerical cost at most that of using the models in the first place. Firstly, we develop and test a dual-weighted residual estimator strategy on a system of advection-diffusion-reaction equations. Secondly, we test that methodology on another system of equations, where quenching and ignition effects are added. In the absence of advection, a rigorous asymptotic analysis shows the existence of many combustion regimes already observed in numerical simulations. We obtain approximations of the quenching and ignition parameters, alongside the S-shaped curve, a plot of the maximal flame temperature as a function of the Damköhler number, consisting of three branches and two bends. When advection effects are added, we still obtain a S-shaped curve corresponding to the known combustion regimes. We compare the modeling errors of the asymptotic approximations in the two stable regimes and establish new model hierarchies for each combustion regime. These errors are compared with the estimations obtained by using the error estimation strategy. When only one stable combustion regime exists, the error estimator correctly identifies that regime; when two or more regimes are possible, it gives a systematic way of choosing one regime. For regimes where more than one model is appropriate, the error estimator’s predicted hierarchy is correct.

This thesis describes the application of the first-order Conditional Moment Closure (CMC) to the autoignition of high-pressure fuel jets, and to piloted and lifted turbulent jet flames using classical and advanced CMC submodels. A Doubly-Conditional Moment Closure (DCMC) formulation is further proposed. In the first study, CMC is applied to investigate the impact of C₂H₆, H₂ and N₂ additives on the autoignition of high-pressure CH₄ jets injected into lower pressure heated air. A wide range of pre-combustion air temperatures is considered and detailed chemical kinetics are employed. It is demonstrated that the addition of C₂H₆ and H₂ does not change the main CH₄ oxidisation pathways. The decomposition of these additives provides additional ignition-promoting radicals, and therefore leads to shorter ignition delays. N₂ additives do not alter the CH₄ oxidisation pathways, however, they reduce the amount of CH₄ available for reaction, causing delayed ignition. It is further shown that ignition always occurs in lean mixtures and at low scalar dissipation rates. The second study is concerned with the modelling of a piloted CH₄/air turbulent jet flame. A detailed assessment of several Probability Density Function (PDF), Conditional Scalar Dissipation Rate (CSDR) and Conditional Velocity (CV) submodels is first performed. The results of two β-PDF-based implementations are then presented. The two realisations differ by the modelling of the CSDR. Homogeneous (inconsistent) and inhomogeneous (consistent) closures are considered. It is shown that the levels of all reactive scalars, including minor intermediates and radicals, are better predicted when the effects of inhomogeneity are included in the modelling of the CSDR. The two following studies are focused on the consistent modelling of a lifted H₂/N₂ turbulent jet flame issuing into a vitiated coflow. Two approaches are followed to model the PDF. In the first, a presumed β-distribution is assumed, whereas in the second, the Presumed Mapping Function (PMF) approach is employed. Fully consistent CV and CSDR closures based on the β-PDF and the PMF-PDF are employed. The homogeneous versions of the CSDR closures are also considered in order to assess the effect of the spurious sources which stem from the inconsistent modelling of mixing. The flame response is analysed over a narrow range of coflow temperatures (Tc). The stabilisation mechanism is determined from the analysis of the transport budgets in mixture fraction and physical spaces, and the history of radical build-up ahead of the stabilisation height. The β-PDF realisations indicate that the flame is stabilised by autoignition irrespective of the value of Tc. On the other hand, the PMF realisations reveal that the stabilisation mechanism is susceptible to Tc. Autoignition remains the controlling stabilisation mechanism for sufficiently high Tc. However, as Tc is decreased, stabilisation is achieved by means of premixed flame propagation. The analysis of the spurious sources reveals that their effect is small but non-negligible, most notably within the flame zone. Further, the assessment of several H₂ oxidation mechanisms show that the flame is very sensitive to chemical kinetics. In the last study, a DCMC method is proposed for the treatment of fluctuations in non-premixed and partially premixed turbulent combustion. The classical CMC theory is extended by introducing a normalised Progress Variable (PV) as a second conditioning variable beside the mixture fraction. The unburnt and burnt states involved in the normalisation of the PV are specified such that they are mixture fraction-dependent. A transport equation for the normalised PV is first obtained. The doubly-conditional species, enthalpy and temperature transport equations are then derived using the decomposition approach and the primary closure hypothesis is applied. Submodels for the doubly-conditioned unclosed terms which arise from the derivation of DCMC are proposed. As a preliminary analysis, the governing equations are simplified for homogeneous turbulence and a parametric assessment is performed by varying the strain rate levels in mixture fraction and PV spaces.

This thesis describes the application of the first-order Conditional Moment Closure (CMC) to the autoignition of high-pressure fuel jets, and to piloted and lifted turbulent jet flames using classical and advanced CMC submodels. A Doubly-Conditional Moment Closure (DCMC) formulation is further proposed. In the first study, CMC is applied to investigate the impact of C₂H₆, H₂ and N₂ additives on the autoignition of high-pressure CH₄ jets injected into lower pressure heated air. A wide range of pre-combustion air temperatures is considered and detailed chemical kinetics are employed. It is demonstrated that the addition of C₂H₆ and H₂ does not change the main CH₄ oxidisation pathways. The decomposition of these additives provides additional ignition-promoting radicals, and therefore leads to shorter ignition delays. N₂ additives do not alter the CH₄ oxidisation pathways, however, they reduce the amount of CH₄ available for reaction, causing delayed ignition. It is further shown that ignition always occurs in lean mixtures and at low scalar dissipation rates. The second study is concerned with the modelling of a piloted CH₄/air turbulent jet flame. A detailed assessment of several Probability Density Function (PDF), Conditional Scalar Dissipation Rate (CSDR) and Conditional Velocity (CV) submodels is first performed. The results of two β-PDF-based implementations are then presented. The two realisations differ by the modelling of the CSDR. Homogeneous (inconsistent) and inhomogeneous (consistent) closures are considered. It is shown that the levels of all reactive scalars, including minor intermediates and radicals, are better predicted when the effects of inhomogeneity are included in the modelling of the CSDR. The two following studies are focused on the consistent modelling of a lifted H₂/N₂ turbulent jet flame issuing into a vitiated coflow. Two approaches are followed to model the PDF. In the first, a presumed β-distribution is assumed, whereas in the second, the Presumed Mapping Function (PMF) approach is employed. Fully consistent CV and CSDR closures based on the β-PDF and the PMF-PDF are employed. The homogeneous versions of the CSDR closures are also considered in order to assess the effect of the spurious sources which stem from the inconsistent modelling of mixing. The flame response is analysed over a narrow range of coflow temperatures (Tc). The stabilisation mechanism is determined from the analysis of the transport budgets in mixture fraction and physical spaces, and the history of radical build-up ahead of the stabilisation height. The β-PDF realisations indicate that the flame is stabilised by autoignition irrespective of the value of Tc. On the other hand, the PMF realisations reveal that the stabilisation mechanism is susceptible to Tc. Autoignition remains the controlling stabilisation mechanism for sufficiently high Tc. However, as Tc is decreased, stabilisation is achieved by means of premixed flame propagation. The analysis of the spurious sources reveals that their effect is small but non-negligible, most notably within the flame zone. Further, the assessment of several H₂ oxidation mechanisms show that the flame is very sensitive to chemical kinetics. In the last study, a DCMC method is proposed for the treatment of fluctuations in non-premixed and partially premixed turbulent combustion. The classical CMC theory is extended by introducing a normalised Progress Variable (PV) as a second conditioning variable beside the mixture fraction. The unburnt and burnt states involved in the normalisation of the PV are specified such that they are mixture fraction-dependent. A transport equation for the normalised PV is first obtained. The doubly-conditional species, enthalpy and temperature transport equations are then derived using the decomposition approach and the primary closure hypothesis is applied. Submodels for the doubly-conditioned unclosed terms which arise from the derivation of DCMC are proposed. As a preliminary analysis, the governing equations are simplified for homogeneous turbulence and a parametric assessment is performed by varying the strain rate levels in mixture fraction and PV spaces.

La combustion cryogénique dans les moteurs de fusée dits à propulsion liquide utilise généralement un couple d'ergols, le plus couramment composé d'hydrogène/oxygène (H2/O2). Privilégiée pour le fort pouvoir calorifique du dihydrogène, cette combustion à haute pression, induit des températures de fonctionnement très élevées et nécessite l'intégration d'un système de refroidissement. La prédiction des flux thermiques aux parois est donc un élément essentiel de la conception d'une chambre de combustion de moteur fusée. Ces flux sont le résultat d'écoulements fortement turbulents, compressibles, avec une cinétique chimique violente induisant de forts gradients d'espèces et de température. La simulation de ces phénomènes nécessite des approches spécifiques telles que la Simulation aux Grandes Echelles (SGE) qui réalise un très bon compromis entre précision et coût de calcul. Cette thèse a ainsi pour objectif la simulation par SGE des transferts de chaleur aux parois dans les chambres de combustion de moteurs fusée opérant en régime sous-critique. Le régime sous-critique implique un état liquide pour un des ergols, dont il faut traiter l'injection et l'atomisation. Dans un premier temps ce travail s'intéresse à plusieurs éléments de modélisation nécessaire pour réaliser les simulations visées. Le comportement des flammes H2/O2 est décrit par un schéma cinétique réduit et validé sur des configurations académiques. La prédictivité de ce schéma est évaluée sur une large gamme de fonctionnement dans des conditions représentatives des moteurs fusée. La simulation de l'injection de l'oxygène liquide (LOx) est un autre point critique qui nécessite de décrire l'atomisation et la phase dispersée ainsi que son couplage avec la phase gazeuse. La déstabilisation et l'atomisation primaire du jet liquide, trop complexe à simuler en SGE 3D, sont omises ici pour injecter directement un spray paramétré grâce à des corrélations empiriques. Enfin, la prédiction des flux thermiques utilise un modèle de loi de paroi spécifiquement dédiée aux écoulements à fort gradient de température. Cette loi de paroi est validée sur des configurations de canaux turbulents par comparaison avec des simulations avec résolution directe de la couche limite. La méthodologie basée sur les modèles développés est ensuite employée pour la simulation d'une chambre de combustion représentative du fonctionnement des moteurs cryogéniques. Il s'agit de la configuration CONFORTH testée sur le banc MASCOTTE (ONERA) et pour laquelle des mesures de température de paroi et de flux thermiques sont disponibles. Les résultats des SGE montrent un bon accord avec l'expérience et démontrent la capacité de la SGE à prédire les flux thermiques dans une chambre de combustion de moteur fusée. Enfin, dans un dernier chapitre ce travail s'intéresse à une méthode d'augmentation des transferts thermiques via une expérience de JAXA utilisant des parois rainurées dans la direction axiale. Par comparaison avec une chambre à parois lisses, les résultats démontrent la bonne prédiction par la SGE de l'augmentation du flux de chaleur grâce aux rainures et confirment la validité de la méthode développée pour des géométries de paroi complexes. / Combustion in cryogenic engines is a complex phenomenon, involving either liquid or supercritical fluids at high pressure, strong and fast oxidation chemistry, and high turbulence intensity. Due to extreme operating conditions, a particularly critical issue in rocket engine is wall heat transfer which requires efficient cooling of the combustor walls. The concern goes beyond material resistance: heat fluxes extracted through the chamber walls may be reused to reduce ergol mass or increase the power of the engine. In expander-type engine cycle, this is even more important since the heat extracted by the cooling system is used to drive the turbo-pumps that feed the chamber in fuel and oxidizer. The design of rocket combustors requires therefore an accurate prediction of wall heat flux. To understand and control the physics at play in such combustor, the Large Eddy Simulation (LES) approach is an efficient and reliable numerical tool. In this thesis work, the objective is to predict wall fluxes in a subcritical rocket engine configuration by means of LES. In such condition, ergols may be in their liquid state and it is necessary to model liquid jet atomization, dispersion and evaporation.The physics that have to be treated in such engine are: highly turbulent reactive flow, liquid jet atomization, fast and strong kinetic chemistry and finally important wall heat fluxes. This work first focuses on several modeling aspects that are needed to perform the target simulations. H2/O2 flames are driven by a very fast chemistry, modeled with a reduced mechanism validated on academic configurations for a large range of operating conditions in laminar pre- mixed and non-premixed flames. To form the spray issued from the atomization of liquid oxygen (LOx) an injection model is proposed based on empirical correlations. Finally, a wall law is employed to recover the wall fluxes without resolving directly the boundary layer. It has been specifically developed for important temperature gradients at the wall and validated on turbulent channel configurations by comparison with wall resolved LES. The above models are then applied first to the simulation of the CONFORTH sub-scale thrust chamber. This configuration studied on the MASCOTTE test facility (ONERA) has been measured in terms of wall temperature and heat flux. The LES shows a good agreement compared to experiment, which demonstrates the capability of LES to predict heat fluxes in rocket combustion chambers. Finally, the JAXA experiment conducted at JAXA/Kakuda space center to observe heat transfer enhancement brought by longitudinal ribs along the chamber inner walls is also simulated with the same methodology. Temperature and wall fluxes measured with smooth walls and ribbed walls are well recovered by LES. This confirms that the LES methodology proposed in this work is able to handle wall fluxes in complex geometries for rocket operating conditions.

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

La Simulation aux Grandes échelles (LES) est de plus en plus présentée comme un outil à part entière dans le développement des chambres de combustion des turbomachines. Dans ce contexte, les écoulements réactifs considérés sont complexes et, dans un souci de validation, la LES doit montrer ses capacités sur des configurations modèles. Le but de cette thèse est de démontrer le potentiel de la LES pour la prédiction des écoulements vrillés réactifs non prémélangés de chambres de combustion modèles. - La LES est tout d'abord appliquée sur une configuration turbulente avec une hypothèse de prémélange parfait, afin d'étudier l'influence de la modélisation de la cinétique chimique, des modèles de combustion turbulente et de leur paramètres internes. Dans ces conditions, chacun de ces modèles montre ses avantages et désavantages. - L'hypothèse de prémélange parfait est ensuite retirée et l'étude réalisée permet d'évaluer l'influence de la prise en compte du mélange air/carburant dans un injecteur vrillé, des pertes thermiques et des conditions limites acoustiques. - Enfin, une chambre de combustion non prémélangée est simulée afin de démontrer les capacités du modèle de flamme épaissie sur ce type de flamme, pour lequel il n'a pas été initialement développé. Les résultats obtenus sont encourageants et démontrent, entre autres, la bonne représentation du positionnement de la flamme. / Large Eddy Simulation (LES) is considered as the next generation tool for the development of turbomachinery combustion chambers. In this specific context, reactive flows are of very complex nature and, as a validation goal, LES needs to prove its capabilities on academic configurations. This dissertation aims at demonstrating LES capabilities for the simulation of non-premixed reactive flows that can be found in swirled academic combustion chambers. - LES is first applied to a turbulent reacting configuration with a perfect premixing assumption. Chemical kinetics, turbulent combustion models and their internal parameters are studied. For this flow condition, each model shows his advantages and disadvantages. - Then, the perfect premixed hypothesis is removed, allowing the evaluation of mixing, thermal losses and acoustic boundary conditions for this swirled injector. - Finally, a non premixed combustion chamber is simulated with the dynamically thickened flame model, which was not developped for this kind of reactive flow. However, results are encouraging and demonstrate that the flame localisation is well represented by LES.

Une méthode de frontières immergées est développée pour la simulation d’écoulements compressibles et validée au travers de cas-tests spécifiques (réflexion d’ondes acoustiques et quantification de la conservation de la masse dans des canaux inclinés). La simulation aux grandes échelles (LES) d’une cavité transsonique est ensuite présentée. Le bouclage aéro-acoustique, très sensible aux conditions aux limites, est reproduit avec précision par la LES dans le cas où les parois sont immergées dans un maillage structurée. La comparaison des stratégies de modélisation de sous-maille pour cet écoulement transsonique et l’adaptation des filtres en présence de frontières immergées sont également discutées. Le rôle, souvent sous-estimé, du schéma de viscosité artificiel, est quantifié.Dans la dernière partie du manuscrit, des études sont réalisées pour aider au dimensionnement d’un nouveau concept de chambre de combustion où la flamme est stabilisée par la recirculation de gaz brûlés dans une cavité (chambre TVC pour Trapped Vortex Combustor). La modélisation de la combustion turbulente est basée sur une chimie tabulée, couplée à une fonction densité de probabilité présumée (PCM-FPI). L’étude de la dynamique de la flamme est réalisée pour diverses conditions de fonctionnement (débit de l’écoulement principal et présence ou non d’un swirl). Les spécificités de mise en œuvre de la simulation d’un écoulement de ce type sont discutées et un soin particulier est apporté au traitement de la condition de sortie, qui constitue un point sensible de la chaîne de modélisation. Les phénomènes d’instabilités et de retour de la flamme sont mis en évidence ainsi que les modifications à apporter au dispositif afin de minimiser ces effets. L’existence d’un cycle limite acoustique est souligné et une formule permettant d’anticiper le niveau des fluctuations de pression est proposée et validée. Une correction au modèle PCM-FPI est présentée afin de préserver la vitesse de flamme et d’assurer une reproduction plus précise de la dynamique de flamme. / An immersed boundary method has been developed for the simulation of compressible flow and validated with reference test cases (pressure wave reflection and quantification of mass conservation for various inclined channels). Large Eddy Simulation (LES) of a transonic cavity is then presented. The aeroacoustic feedback loop, which is highly sensitive to the boundary conditions, was accurately reproduced where the walls are immersed inside a structured grid. The comparison between the modeling approaches for this transonic flow and the correction of the filtering operation near immersed boundaries are also discussed. The often underestimated role of the numerical artificial dissipation is also quantified.In the last part of this manuscript, many studies are realized to help in the design of a new combustion chamber for Trapped Vortex Combustor (TVC). The turbulent combustion model is based on tabulated chemistry and a presumed probability density function (PCM-FPI) method.The flame dynamics is studied for various operating conditions (flowrate of the main flow and presence of swirl motion). Details concerning the realization of such a flow are discussed and special care is taken for the treatment of the most sensitive outlet boundary condition. The phenomena of combustion instabilities and of flame backflow are highlighted along with the modifications to be made for the device to minimize these effects. The existence of a acoustic limit cycle is emphasized and a formula is proposed and validated to anticipate the level of pressure fluctuations. Finally a correction to the PCM-FPI model is suggested to preserve the flame front speed and to ensure a more accurate description of the flame dynamics.