The accurate description of the turbulence chemistry interactions that can determine chemical conversion rates and flame stability in turbulent combustion modelling is a challenging research area. This thesis presents the development and implementation of a model for the treatment of fluctuations around the conditional mean (i.e., the auto-ignition and extinction phenomenon) of realistic turbulence-chemistry interactions in computational fluid dynamics (CFD) software. The wider objective is to apply the model to advanced combustion modelling and extend the present analysis to larger hydrocarbon fuels and particularly focus on the ability of the model to capture the effects of particulate formation such as soot. A comprehensive approach for modelling of turbulent combustion is developed in this work. A direct quadrature conditional moment closure (DQCMC) method for the treatment of realistic turbulence-chemistry interactions in computational fluid dynamics (CFD) software is described. The method which is based on the direct quadrature method of moments (DQMOM) coupled with the Conditional Moment Closure (CMC) equations is in simplified form and easily implementable in existing CMC formulation for CFD code. The observed fluctuations of scalar dissipation around the conditional mean values are captured by the treatment of a set of mixing environments, each with its pre-defined weight. In the DQCMC method the resulting equations are similar to that of the first-order CMC, and the “diffusion in the mixture fraction space” term is strictly positive and no correction factors are used. Results have been presented for two mixing environments, where the resulting matrices of the DQCMC can be inverted analytically. Initially the DQCMC is tested for a simple hydrogen flame using a multi species chemical scheme containing nine species. The effects of the fluctuations around the conditional means are captured qualitatively and the predicted results are in very good agreement with observed trends from direct numerical simulations (DNS). To extend the analysis further and validate the model for larger hydrocarbon fuel, the simulations have been performed for n-heptane flame using detailed multi species chemical scheme containing 67 species. The hydrocarbon fuel showed improved results in comparison to the simple hydrogen flame. It suggests that higher hydrocarbons are more sensitive to local scalar dissipation rate and the fluctuations around the conditional means than the hydrogen. Finally, the DQCMC is coupled with a semi-empirical soot model to study the effects of particulate formation such as soot. The modelling results show to predict qualitatively the trends from DNS and are in very good agreement with available experimental data from a shock tube concerning ignition delays time. Furthermore, the findings suggest that the DQCMC approach is a promising framework for soot modelling.

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.

Cette thèse est consacrée à l’amélioration des capacités de prédiction des émissions polluantes (CO, NOx . . . ) des foyers de combustion de turboréacteurs. L’étude, exclusivement numérique, repose sur des simulations aux grandes échelles (ou LES pour Large-Eddy Simulation) basées sur des méthodes de tabulation de la chimie détaillée. L’approche PCM-FPI (pour Presumed Conditional Moments - Flame Prolongation of ILDM) a été étendue à la chimie des oxydes d’azote via la modification de la variable d’avancement. Différentes validations sur des configurations laminaires simples puis des flammes de laboratoire (Cabra, Sandia) sont présentées. Les résultats en terme de structure de flamme et de champs d’espèces chimiques sont confrontés aux mesures expérimentales. Le rôle du formaldéhyde comme marqueur de la zone réactionnelle est illustré à l’aide de calculs de flammes laminaires puis confirmé par un calcul 3D LES. Une analyse des spécificités de l’implantation de ce type de modèle sur des machines à architecture massivement parallèle est ensuite menée. Diverses modifications de la structure de la table et des méthodes d’interpolation sont réalisées, servant de base à une étude de sensibilité de maillage appliquée à la flamme Sandia D. Les difficultés relatives à la prédiction du NO dans les flammes turbulentes sont exposées : divers modèles de sous-maille sont alors employés et comparés. / This thesis is focused on the prediction capabilities of pollution (CO, NOx especially) for numerical tools applied to aeronautical combustion chambers. The modeling work is based on Large-Eddy Simulation methods coupled with a tabulated detailed chemistry approach. The PCM-FPI model, which stands for Presumed Conditional Moments - Flame Prolongation of ILDM, has been revised to take into account nitrogen chemistry through a modification of the progress variable. Several benchmarks and test-cases (laminar and turbulent flames) are gathered in this study : Cabra and Sandia flames have been computed and satisfactory comparisons with measurements are obtained. The role of CH2O as a marker of heat release is investigated, first in the frame of laminar premixed flames and then validated through LES runs. The challenges of the implementation of tabulated chemistry methods on massively parallel machines are discussed. Modifications are proposed regarding both the table structure and the interpolation methods leading to a mesh sensitivity review applied to the Sandia D flame. Difficulties arising when dealing with NOx chemistry in turbulent flows are presented : new Sub-Grid Scale models are introduced and investigated.