Modelos algébricos para prever o comprimento de uma chama turbulenta têm sido foco de estudo de diversos grupos de pesquisa por suas aplicações na área de engenharia. O método experimental para obter o modelo varia desde visualizações simples, até técnicas fotográficas, este último com parâmetros fotográficos variando entre os autores. Técnicas fotográficas são usadas para estimar a altura de levantamento da base da chama, (Lift-Off) e o comprimento médio visível de chama (Visible Flame Length, VFL). Duas técnicas comuns que podem ser encontradas na literatura: por imagens de chama com baixo tempo de exposição e longo tempo de exposição, são comparados com um terceiro que se baseia na intensidade luminosa e na frequência de imagens de chama que ocupam um pixel. O melhor método foi utilizado para caracterizar o comportamento das chamas turbulentas de gás natural para diferentes regimes de velocidade do escoamento. Modelos algébricos que preveem o comprimento de chama, altura de levantamento e a velocidade crítica de extinção de chama são avaliados com os novos resultados experimentais. Logo após, os coeficientes numéricos dos melhores modelos algébricos são reajustados Finalmente, foram obtidos mapas de estabilidade relacionados à altura de levantamento e à velocidade crítica de extinção de chama para cada diâmetro em função da diluição com CO2 e do número adimensional de Reynolds. A terceira parte deste trabalho está focada no estudo da distribuição de radiação térmica. Em particular, foram consideradas três distâncias radiais medidas em comprimentos de chama (0,5 Lf, 1 Lf, 2 Lf) visando obter a distribuição do fluxo radiante experimental ao longo de um eixo vertical adjacente às chamas. Finalmente, os dados experimentais foram utilizados como dados de entrada em uma análise inversa com o objetivo de calcular os fatores de ponderação do modelo das múltiplas fontes ponderadas (por suas siglas em inglês WMPS). Nesta última parte, são apresentados frações radiantes e distribuições de fluxo de calor radiante de chamas de gás natural diluídas para diversas diluições com dióxido de carbono e diâmetros do queimador. / Predicting models for turbulent diffusion flame lengths have several applications driven the attention of many research groups. Since several studies use photographs to measure the flame length, with photographic parameters varying among authors, in other cases simple visualizations were used. It is important to explore possible discrepancies among measurement technics that could affect the results. Optical visualizations of turbulent diffusion flames are used to estimate the visible average flame length (VFL) and the lift-off. The study presents a study of three different methods to measure the VFL using optical techniques. The effect on the image of the main optic parameters such as focus, exposure time and ISO sensibility are analyzed. The VFL obtained with images in low exposure time and long exposure time are compared with a third optical method that is based on the luminous intensity and the frequency of flame images occupying a pixel. One method was used to characterize the behavior of turbulent diffusion flames of natural gas for a range of flames in function of the flow velocity. Universal non-dimensional models that describe the VFL, lift-off and the blow-out stability limit of gaseous jet diffusion flames in the still air have been compared with new experimental data. The numerical coefficients of the best models are adjusted. Finally, maps of stability related to lift-off and blow-out were obtained for each diameter in function of the dilution with CO2 and flow exit velocity expressed in non-dimensional Reynolds number The third part of this work focuses on the estimation of the thermal distribution of radiative flux from turbulent diffusion flames in laboratory-scale. The experimental measurements were gotten from the previous stability study. In particular, was considered three radial distances measured in flame lengths (0,5 Lf, 1 Lf, 2 Lf) aiming at obtaining the experimental radiant flux along a vertical axis adjacent to the flames. Finally, the experimental data was used as input data in an inverse analysis with the purpose of computing weight coefficients of the weighted multi-point source (WMPS) model. Then, experimental data that include: radiant fractions and radiative heat flux are presents for several flames with different dilutions with carbon dioxide and burner´s diameters.

Modelos algébricos para prever o comprimento de uma chama turbulenta têm sido foco de estudo de diversos grupos de pesquisa por suas aplicações na área de engenharia. O método experimental para obter o modelo varia desde visualizações simples, até técnicas fotográficas, este último com parâmetros fotográficos variando entre os autores. Técnicas fotográficas são usadas para estimar a altura de levantamento da base da chama, (Lift-Off) e o comprimento médio visível de chama (Visible Flame Length, VFL). Duas técnicas comuns que podem ser encontradas na literatura: por imagens de chama com baixo tempo de exposição e longo tempo de exposição, são comparados com um terceiro que se baseia na intensidade luminosa e na frequência de imagens de chama que ocupam um pixel. O melhor método foi utilizado para caracterizar o comportamento das chamas turbulentas de gás natural para diferentes regimes de velocidade do escoamento. Modelos algébricos que preveem o comprimento de chama, altura de levantamento e a velocidade crítica de extinção de chama são avaliados com os novos resultados experimentais. Logo após, os coeficientes numéricos dos melhores modelos algébricos são reajustados Finalmente, foram obtidos mapas de estabilidade relacionados à altura de levantamento e à velocidade crítica de extinção de chama para cada diâmetro em função da diluição com CO2 e do número adimensional de Reynolds. A terceira parte deste trabalho está focada no estudo da distribuição de radiação térmica. Em particular, foram consideradas três distâncias radiais medidas em comprimentos de chama (0,5 Lf, 1 Lf, 2 Lf) visando obter a distribuição do fluxo radiante experimental ao longo de um eixo vertical adjacente às chamas. Finalmente, os dados experimentais foram utilizados como dados de entrada em uma análise inversa com o objetivo de calcular os fatores de ponderação do modelo das múltiplas fontes ponderadas (por suas siglas em inglês WMPS). Nesta última parte, são apresentados frações radiantes e distribuições de fluxo de calor radiante de chamas de gás natural diluídas para diversas diluições com dióxido de carbono e diâmetros do queimador. / Predicting models for turbulent diffusion flame lengths have several applications driven the attention of many research groups. Since several studies use photographs to measure the flame length, with photographic parameters varying among authors, in other cases simple visualizations were used. It is important to explore possible discrepancies among measurement technics that could affect the results. Optical visualizations of turbulent diffusion flames are used to estimate the visible average flame length (VFL) and the lift-off. The study presents a study of three different methods to measure the VFL using optical techniques. The effect on the image of the main optic parameters such as focus, exposure time and ISO sensibility are analyzed. The VFL obtained with images in low exposure time and long exposure time are compared with a third optical method that is based on the luminous intensity and the frequency of flame images occupying a pixel. One method was used to characterize the behavior of turbulent diffusion flames of natural gas for a range of flames in function of the flow velocity. Universal non-dimensional models that describe the VFL, lift-off and the blow-out stability limit of gaseous jet diffusion flames in the still air have been compared with new experimental data. The numerical coefficients of the best models are adjusted. Finally, maps of stability related to lift-off and blow-out were obtained for each diameter in function of the dilution with CO2 and flow exit velocity expressed in non-dimensional Reynolds number The third part of this work focuses on the estimation of the thermal distribution of radiative flux from turbulent diffusion flames in laboratory-scale. The experimental measurements were gotten from the previous stability study. In particular, was considered three radial distances measured in flame lengths (0,5 Lf, 1 Lf, 2 Lf) aiming at obtaining the experimental radiant flux along a vertical axis adjacent to the flames. Finally, the experimental data was used as input data in an inverse analysis with the purpose of computing weight coefficients of the weighted multi-point source (WMPS) model. Then, experimental data that include: radiant fractions and radiative heat flux are presents for several flames with different dilutions with carbon dioxide and burner´s diameters.

This dissertation addresses the analysis of partially premixed flame configurations and the detection and characterization of their local flame regimes. First, the identification of flame regimes in experimental data is intensively discussed. Current methods for combustion regime characterization, such as the flame index, rely on 3D gradient information that is not accessible with available experimental techniques. Here, a method is proposed for reaction zone detection and characterization, which can be applied to instantaneous 1D Raman/Rayleigh line measurements of major species and temperature as well as to the results of laminar and turbulent flame simulations, without the need for 3D gradient information. Several derived flame markers, namely the mixture fraction, the heat release rate and the chemical explosive mode, are combined to detect and characterize premixed versus non-premixed reaction zones. The methodology is developed and evaluated using fully resolved simulation data from laminar flames. The fully resolved 1D simulation data are spatially filtered to account for the difference in spatial resolution between the experiment and the simulation, and experimental uncertainty is superimposed onto the filtered numerical results to produce Raman/Rayleigh equivalent data. Then, starting from just the temperature and major species, a constrained homogeneous batch reactor calculation gives an approximation of the full thermochemical state at each sample location. Finally, the chemical explosive mode and the heat release rate are calculated from this approximated state and compared to those calculated directly from the simulation data. After successful validation, the approach is applied to Raman/Rayleigh line measurements from laminar counterflow flames, a mildly turbulent lifted flame and turbulent benchmark cases. The results confirm that the reaction zones can be reliably detected and characterized using experimental data. In contrast to other approaches, the presented methodology circumvents uncertainties arising from the use of limited gradient information and offers an alternative to known reaction zone identification methods. Second, this work focuses on the flame structure of partially premixed dimethyl ether (DME) flames. DME flames form significant intermediate hydrocarbons in the reaction zone and are classified as the next more complex fuel candidate in research after methane. To simulate DME combustion processes, accurate predictions by computational combustion models are required. To evaluate such models and to identify appropriate flame regimes, numerical simulations are necessary. Therefore, fully resolved simulations of laminar dimethyl ether flames, defined by different levels of premixing, are performed. Further, the qualitative two-dimensional structures of the partially premixed DME flames are discussed and analyses are carried out at selected slices and compared to each other as well as to experimental data. Further, the flamelet/progress variable (FPV) approach is investigated to predict the partially premixed flame structures of the DME flames. In the context of the FPV approach, a rigorous analysis of the underlying manifold is carried out based on the newly developed regime identification approach and an a priori analysis. The most promising flamelet look-up table is chosen for the fully coupled tabulated chemistry simulations and the results are further compared to the fully resolved simulation data.

Fossil fuel reserves are projected to be decreasing, and emission regulations are becoming more stringent due to increasing atmospheric pollution. Alternative fuels for power generation in industrial gas turbines are thus required able to meet the above demands. Examples of such fuels are synthetic gas, blast furnace gas and coke oven gas. A common characteristic of these fuels is that they are multi-component fuels, whose composition varies greatly depending on their production process. This implies that their combustion characteristics will also vary significantly. Thus, accurate and yet flexible enough combustion sub-models are required for such fuels, which are used during the design stage, to ensure optimum performance during practical operating conditions. Most combustion sub-model development and validation is based on Direct Numerical Simulation (DNS) studies. DNS however is computationally expensive. This, has so far limited DNS to single-component fuels such as methane and hydrogen. Furthermore, the majority of DNS conducted to date used one-step chemistry in 3D, and skeletal chemistry in 2D only. The need for 3D DNS using skeletal chemistry is thus apparent. In this study, an accurate reduced chemical mechanism suitable for multi-component fuel-air combustion is developed from a skeletal mechanism. Three-dimensional DNS of a freely propagating turbulent premixed flame is then conducted using both mechanisms to shed some light into the flame structure and turbulence-scalar interaction of such multi-component fuel flames. It is found that for the multi-component fuel flame heat is released over a wider temperature range contrary to a methane flame. This, results from the presence of individual species reactions zones which do not all overlap. The performance of the reduced mechanism is also validated using the DNS data. Results suggest it to be a good substitute of the skeletal mechanism, resulting in significant time and memory savings. The flame markers commonly used to visualize heat release rate in laser diagnostics are found to be inadequate for the multi-component fuel flame, and alternative markers are proposed. Finally, some popular mean reaction rate closures are tested for the multi-component fuel flame. Significant differences are observed between the models’ performance at the highest turbulence level considered in this study. These arise from the chemical complexity of the fuel, and further parametric studies using skeletal chemistry DNS would be useful for the refinement of the models.

Cette thèse est consacrée à l’amélioration des capacités de prédiction des émissions d’oxydes d’azote (NO et NO2) des foyers de combustion aéronautiques. Les travaux, exclusivement numériques, consistent d’abord dans une étude de la cinétique chimique responsable des émissions polluantes. Cetteétude conduit à l’écriture d’un modèle, nommé NOMANI (pour Nitrogen Oxide emission model with one-dimensional MANIfold), basé sur l’approche PCM-FPI (pour Presumed Conditional Moments - Flame Prolongation of ILDM) avec une variable de progrès additionnelle afin calculer l’avancement de la chimie azotée une fois la chimie carbonée à l’équilibre. Différentes validations sur des configurations laminaires simples puis des flammes de laboratoire de Sandia sont présentées. Les résultats en terme de structure de flamme et d'émission de monoxyde d’azote sont confrontés aux mesures expérimentales. Le dernier volet de ces travaux, disponible uniquement dans la version confidentielle du manuscrit, consiste dans le développement d’un modèle de prédiction de polluants associé au modèle TF-LES (pour Thickening Flame for Large Eddy Simulation). Le modèle développé est ensuite appliqué à des calculs d’une chambre de combustion aéronautique. / This thesis is focused on the prediction capabilities of nitrogen oxides (NO and NO2) for numerical tools applied to aeronautical combustion chambers. The modeling work is based on a study of the chemical kinetic that produced the pollutant emissions. This study leads to a model, called NOMANI (Nitrogen Oxide emission model with one-dimensional MANIfold), based on PCM-FPI (Presumed Conditional Moments - Flame Prolongation of ILDM) with an additional progress variable to compute the NO evolution once the carbon chemistry is at the equilibrium. Several benchmarks and test-cases (laminar and turbulent flames) are gathered in this study : Sandia flame have been computed and satisfactory comparisons with measurements are obtained. The last part of this work, only available in the confidential version of the manuscript, is the development of a model to predict pollutant associated with the model TF-LES (for Thickening Flame for Large Eddy Simulation). This model is then applied to computations of a aeronautical combustion chambers.

La combustion stratifiée, qui consiste à brûler un mélange carburant/oxydant inhomogène, et la combustion diluée, consistant à ajouter une quantité limitée de gaz brûlés, sont deux technologies utilisées dans les moteurs à piston pour réduire leur consommation. Cette thèse est dédiée à l’étude de l’allumage dans ces deux types de milieux en régimes laminaire et turbulent. Un nouveau schéma cinétique pour la combustion propane/air a été dérivé et combiné à deux approches de modélisation différentes : la chimie complexe et une approche de chimie tabulée de type FPI. Dans le cas laminaire, les deux approches de modélisation donnent des résultats similaires et un modèle simple a mis en évidence l’importance de la dynamique des gaz frais et des gaz brûlés sur le développement du noyau. Dans le cas turbulent, plusieurs techniques d’analyse ont montré la dépendance de la vitesse absolue de la flamme au champ de vitesse moyen et la décorrélation des fluctuations locales de richesse. / Stratified combustion, which consists in burning an inhomogeneous fuel/air mixture, and diluted combustion, which consists in adding a limited quantity of burnt gases, are two technologies used in internal combustion engines to reduce fuel consumption. This Ph.D is devoted to the study of ignition in these two types of combustion in laminar and turbulent regimes. A new kinetic scheme for propane/air combustion has been derived and combined to two modeling approaches: finite-rate chemistry and an FPI tabulated chemistry approach. In the laminar case, both approaches give similar results and a simplified model has highlighted the importance of fresh and burnt gases dynamics on the kernel development. In the turbulent case, several techniques of analysis have shown the dependency of absolute flame speed on the mean fluid velocity and the lack of correlation to the local equivalence ratio.

This dissertation addresses the analysis of partially premixed flame configurations and the detection and characterization of their local flame regimes. First, the identification of flame regimes in experimental data is intensively discussed. Current methods for combustion regime characterization, such as the flame index, rely on 3D gradient information that is not accessible with available experimental techniques. Here, a method is proposed for reaction zone detection and characterization, which can be applied to instantaneous 1D Raman/Rayleigh line measurements of major species and temperature as well as to the results of laminar and turbulent flame simulations, without the need for 3D gradient information. Several derived flame markers, namely the mixture fraction, the heat release rate and the chemical explosive mode, are combined to detect and characterize premixed versus non-premixed reaction zones. The methodology is developed and evaluated using fully resolved simulation data from laminar flames. The fully resolved 1D simulation data are spatially filtered to account for the difference in spatial resolution between the experiment and the simulation, and experimental uncertainty is superimposed onto the filtered numerical results to produce Raman/Rayleigh equivalent data. Then, starting from just the temperature and major species, a constrained homogeneous batch reactor calculation gives an approximation of the full thermochemical state at each sample location. Finally, the chemical explosive mode and the heat release rate are calculated from this approximated state and compared to those calculated directly from the simulation data. After successful validation, the approach is applied to Raman/Rayleigh line measurements from laminar counterflow flames, a mildly turbulent lifted flame and turbulent benchmark cases. The results confirm that the reaction zones can be reliably detected and characterized using experimental data. In contrast to other approaches, the presented methodology circumvents uncertainties arising from the use of limited gradient information and offers an alternative to known reaction zone identification methods. Second, this work focuses on the flame structure of partially premixed dimethyl ether (DME) flames. DME flames form significant intermediate hydrocarbons in the reaction zone and are classified as the next more complex fuel candidate in research after methane. To simulate DME combustion processes, accurate predictions by computational combustion models are required. To evaluate such models and to identify appropriate flame regimes, numerical simulations are necessary. Therefore, fully resolved simulations of laminar dimethyl ether flames, defined by different levels of premixing, are performed. Further, the qualitative two-dimensional structures of the partially premixed DME flames are discussed and analyses are carried out at selected slices and compared to each other as well as to experimental data. Further, the flamelet/progress variable (FPV) approach is investigated to predict the partially premixed flame structures of the DME flames. In the context of the FPV approach, a rigorous analysis of the underlying manifold is carried out based on the newly developed regime identification approach and an a priori analysis. The most promising flamelet look-up table is chosen for the fully coupled tabulated chemistry simulations and the results are further compared to the fully resolved simulation data.

<div>This research has been studied a laser-induced spark in flowing gases. The relationship between the minimum ignition energy (MIE), the turbulence intensity, and the flame kernel propagation speed is considered. Plasma emission, produced by the laser-induced spark, and flame kernel generation by the plasma are investigated. The energy balance equation between an ignition energy and energy losses by heat transfer is studied at laminar flows and turbulent flows. Hydrogen and air mixtures were used in a premixed jet burner for ignition experiments. Particle image velocimetry (PIV) examined the velocity and the turbulence intensity under the turbulent flows. The flame kernel development was visualized using Schlieren imaging and infrared images (IR camera). Flame kernel temperatures were measured through Rayleigh scattering and infrared images (IR camera). Plasma evaluations were captured through an intensified CCD camera (ICCD camera). Minimum ignition energies were measured at the laminar flows and the turbulent flows. The MIE decreases with an increase in the turbulence intensity which changed by ignition locations and perforated plates at the constant bulk velocity. Improved mixing rates due to the ignition locations or the geometry of the perforated plates decrease the MIE at the constant bulk velocity. The turbulence intensity increases wrinkles in the flame kernel surface, thus the contact between the flame kernel and reactants increases due to the wrinkles. Therefore, the flame kernel propagation speed increases as the turbulence intensity is higher since the increased reaction by the wrinkles and the contact. Thus, the MIE decreases as the turbulence intensity increases at the constant ignition condition, including bulk velocities and ignition heights, since the high turbulence intensity increases the flame kernel propagation speed. Laser energy differences affect the plasma expansions by the laser absorption. Laser-supported radiation (LSR) wave speeds were measured and calculated using energy balance equations. Velocity does not affect the flame kernel temperature distribution during the early reaction steps because the plasma generates a flame kernel and determines the flame kernel temperature distribution. The MIE increases with increasing the bulk velocity. The energy losses considering convection, conduction, and radiation were calculated using the flame kernel radius, the flame kernel temperature, mixture properties, and the flame speed. The energy balance equation in the ignition of flowing gases is newly written at the laminar flows and the turbulent flows.</div>

The present thesis deals with the fundamental aspects of turbulent mixingand non-premixed combustion in wall-jet flows. Direct numerical simulations(DNS) of compressible turbulent flows are performed in a wall-jet configura-tion, which has a close resemblance to many industrial combustion applica-tions. The triple ”turbulence-chemistry-wall” interactions are also present inthis flow set-up. These interactions have been addressed by first focusing onturbulent flow effects on the isothermal reaction, including the near-wall issues.Then, by adding heat-release to the simulations, it has been concentrated onheat-release effects on various phenomena that occur in the reacting turbulentwall-jet flow. In the computational domain, fuel and oxidizer enter separatelyin a non-premixed manner and the flow is fully turbulent and subsonic in allsimulations. In the first phase of this study, the case of a turbulent wall-jetincluding an isothermal reaction without heat release is addressed in order toisolate the near-wall effects and the mixing characteristics of the flow and thekey statistics for combustion are studied in the absence of thermal effects. Adeeper insight into three-dimensional mixing and reaction characteristics in aturbulent wall-jet has been gained through investigation of the probability den-sity functions, higher order moments of velocities and reacting scalars and thescalar dissipation rates of different species. In the second phase, DNS of turbu-lent reacting wall-jets including heat release is performed, where a single-stepglobal exothermic reaction with an Arrhenius-type reaction rate is considered.The main target was to identify the heat-release effects on different mixingscales of turbulent wall-jet flow. The scalar dissipation rates, time scale ratios,two-point correlations, one and two-dimensional premultiplied spectra are usedto illustrate the heat release induced modifications. It is observed that heatrelease effects delay the transition process in the chemically reacting cases andenlarge the fluctuation intensities of density and pressure, but have a dampingeffect on all velocity fluctuation intensities. Finer small mixing scales were ob-served in the isothermal simulations and larger vortical structures formed afteradding significant amounts of heat-release. Simulations with different Damk ̈h-  oler numbers, but comparable temperature-rise are performed and the expectedbehavior, a thinner flame with increasing Damk ̈hler number, is observed. Finally, some heat transfer related quantities are examined. The wall heat fluxand the corresponding Nusselt numbers are addressed. The near-wall reactioneffects on the skin friction coefficient are studied and further the reaction char-acteristics are investigated throughout the domain. / QC 20110908

[pt] Chamas pré-misturadas em escoamentos turbulentos com rotação são encontradas em diversos sistemas de engenharia, como turbinas a gás e motores a jato. Neste trabalho, regimes de chamas característicos de tais sistemas são estudados numericamente num queimador de escala laboratorial. O estado da arte dos estudos numéricos de tais chamas é revisado, com respeito a simulações de grandes escalas, bem como o de modelos computacionais baseados em médias de Reynolds. Um estudo isotérmico é feito no escoamento turbulento, num domínio computacional que consiste de um swirler radial e uma câmara de combustão. O impacto de diferentes modelos de turbulência, níveis de refinamentos de malha e condições de contorno no número de swirl e na estrutura do escoamento é investigado. Os resultados revelam que os três modelos de turbulência propostos resultam em campos de escoamento e número de swirl similares, enquanto o nível de refinamento de malha e a condição de contorno de parede deslizante alteram o número de swirl significativamente. Utilizando as equações de média de Reynolds, com o fechamento do modelo k − E realizável, acoplado a um modelo de duas equações para chamas pré-misturadas de metano e ar, dois regimes de chamas são analisados. Estes regimes correspondem à chama de recirculação externa (chama tipo M) e um regime de instabilidade, que ocorre na transição entre a chama tipo V e chama tornado. A estrutura do escoamento é caracterizada em termos de velocidade e propriedades de turbulência e combustão. Uma comparação entre variáveis de progresso também é feita, utilizando resultados experimentais prévios, levando a boa concordância qualitativa para os dois regimes estudados. / [en] Lean premixed turbulent swirling flames are found in many engineering systems, such as gas turbines and jet engines. This work aims to numerically study flame regimes, representative of such systems, stabilized in a laboratory scale burner. The state of the art of the numerical studies concerning these types of flames is reviewed, with respect to Reynolds-Averaged NavierStokes and Large Eddy Simulations. A turbulent, isothermal flow study is performed within the radial swirler and the combustion chamber. The impact of different turbulence models (realizable k − E, RNG k − E and SST k−W), mesh refinement levels and boundary conditions on the swirl number and overall flow structure is investigated. The results show that the three tested turbulence models yield similar results, with respect to the obtained flow field, whereas the mesh refinement level and slip wall boundary condition alter the swirl number significantly. Using Reynolds-Averaged NavierStokes transport equations, closed by the realizable k − E model, coupled with a two-equation premixed combustion model for methane/air mixtures, two combustion regimes are analyzed. These regimes correspond to the outer recirculation zone flame (M-shaped flame) and an unstable regime, which occurs at the transition between the V-shaped flame and tornado-flame. The flow structure is characterized in terms of velocity fields, turbulence and combustion properties. A reaction progress variable comparison is also performed, using existing experimental results, yielding qualitatively similar results for both studied regimes.