Numerical simulations have proven itself as a significant and powerful tool for accurate prediction of turbulent premixed flames in practical engineering devices. The work presented in this thesis concerns the development of simulation techniques for premixed turbulent combustion of three different fuels, namely, CNG, LPG and Hydrogen air mixtures. The numerical results are validated against published experimental data from the newly built Sydney combustion chamber. In this work a newly developed Large Eddy Simulation (LES) CFD model is applied to the new Sydney combustion chamber of size 50 x 50 x 250 mm (0.625 litre volume). Turbulence is generated in the chamber by introducing series of baffle plates and a solid square obstacle at various axial locations. These baffles can be added or removed from the chamber to adapt various experimental configurations for studies. This is essential to understand the flame behaviour and the structure. The LES numerical simulations are conducted using the Smagorinsky eddy viscosity model with standard dynamic procedures for sub-grid scale turbulence. Combustion is modelled by using a newly developed dynamic flame surface density (DFSD) model based on the flamelet assumption. Various numerical tests are carried out to establish the confidence in the LES based combustion modelling technique. A detailed analysis has been carried out to determine the regimes of combustion at different stages of flame propagation inside the chamber. The predictions using the DFSD combustion model are evaluated and validated against experimental measurements for various flow configurations. In addition, the in-house code capability is extended by implementing the Lewis number effects. The LES predictions are identified to be in a very good agreement with the experimental measurements for cases with high turbulence levels. However, some disagreement were observed with the quasi-laminar case. In addition a data analysis for experimental data, regarding the overpressure, flame position and the flame speed is carried out for the high and low turbulence cases. Moreover, an image processing procedure is used to extract the flame rate of stretch from both the experimental and numerical flame images that are used as a further method to validate the numerical results. For the grids under investigation, it is concluded that the employed grid is independent of the filter width and grid resolution. The applicability of the DFSD model using grid-independent results for turbulent premixed propagating flames was examined by validating the generated pressure and other flame characteristics, such as flame position and speed against experimental data. This study concludes that the predictions using DFSD model provide reasonably good results. It is found that LES predictions were slightly improved in predicting overpressure, flame position and speed by incorporating the Lewis number effect in the model. Also, the investigation demonstrates the effects of placing multiple obstacles at various locations in the path of the turbulent propagating premixed flames. It is concluded that the pressure generated in any individual configuration is directly proportional to the number of baffles plates. The flame position and speed are clearly dependent on the number of obstacles used and their blockage ratio. The flame stretch extracted from both the experimental and numerical images shows that hydrogen has the highest stretch values over CNG and LPG. Finally, the regime of combustion identified for the three fuels in the present combustion chamber is found to lie within the thin reaction zone. This finding supports the use of the laminar flamelet modelling concept that has been in use for the modelling of turbulent premixed flames in practical applications.

Le fonctionnement en régime très pauvre ou avec forts taux d'EGR des moteurs à allumage commandé (MAC) permet de réduire efficacement les émissions de CO2 et de Nox ; cependant ces stratégies se heurtent à l'augmentation des variabilités cycliques. Ces dernières sont principalement dues à la phase d'allumage qui devient critique de dilution. Le modèle ECFM-LES actuellement utilisé à IFPEn, basé sur la notion de densité de surface de flamme, est insuffisant pour décrire l'allumage dans ces conditions critiques. Dans ces travaux, l'approche TF-LES est adoptée, l'allumage étant alors décrit par un emballement cinétique des réactions chimiques lors d'une élévations locale de la température. Ces travaux définissent et évaluent une stratégie de simulation pour TF-LES en configuration moteur, qui permette une prédiction fine des allumages critiques et de la propagation turbulente de la flamme, afin de décrire le cycle moteur complet.Dans une première partie, des DNS d'allumages turbulents ont été réalisées, en modélisant la phase d'allumage par un dépôt d'énergie thermique (Lacaze et al., (2009)). Les calculs ont simulé les expériences d'allumage de Cardin et al. (2013), dans lesquelles l'énergie minimum d'allumage (MIE) d'un mélange mtéhane-air a été mesuré, pour différentes richesses pauvres et sous différentes intensités turbulentes. L'objectif principal des simulations a été de déterminer les paramètres numériques et physiques du modèle permettant de reproduire les allumages de l'expérience. Deux types de schémas cinétiques ont été évalués : un schéma simplifié et un schéma analytique (ARC), ce dernier reproduisant et les délais d'auto-allumage et la vitesse de flamme laminaire. Les résultats ont permis de définir des critères d'allumage et de mettre en évidence les différentes prédiction d'allumage avec les deux types de schémas cinétiques. Les résultats ont été également démontré que l'approche choisie permettait de prédire les bons niveaux d'énergie pour les allumages laminaires et à faible nombres de Kalovitz (Ka<10). Aux plus hauts nombres de Karlovitz, il a été montré que le modèle ED était insuffisant pour prédire les énergie d'allumage et qu'une description plus fine du dépôt d'énergie est nécessaire.Dans la seconde partie des travaux, un modèle de plissement dynamique (Wang et al., 2012) a été étudié, afin de décrire le développement hors-équilibre de la flamme dans la phase de propagation turbulente. Des études sur des flammes sphériques laminaires ont d'abord été menées. Ensuite, les premiers tests de configuration moteur ayant révélé des incompatibilités du modèle, des modifications ont été proposées. Le modèle de plissement dynamique modifié a été finalement évalué sur la configuration moteur ICAMDAC. Les résultats obtenus ont été comparés aux résultats obtenus par Robert et al. (2015) avec le modèle ECFM-LES, qui utilise une équation de transport de densité de surface de flamme décrivant le plissement hors-équilibre de la flamme. Les résultats obtenus avec le plissement dynamique sont en très bon accord avec ceux du modèles ECFM-LES, démontrant ainsi la capacité du modèle dynamique à prédire des valeurs de plissement hors-équilibre. D'autre part, le modèle dynamique s'ajustant automatiquement aux conditions de turbulence de l'écoulement, nul besoin n'est d'ajuster la constante de modélisation en fonction du régime moteur, comme c'est le cas pour l'équation de transport de la densité de surface de flamme. / The use of lean equivalence ratios or high EGR rates in spark ignition engines (SIE) enables to optimize CO2 and NOx emissions; however too important dilution rates leads to increased cycle-to-cycle variability. These latter are mostly due to the ignition phase, which becomes critical when dilution rates are important and requires high ignition energy. The ECFM-LES model currently used in IFPEN, which is based on the flame surface density concept, is not sufficient to describe ignition in these critical conditions. The TF-LES approach was chosen in this study, principally because it directly resolved chemistry and can thus model ignition via a local raise of the temperature. The present work defines and evaluates a simulation strategy for TF-LES in SIE configurations, that enables a fine prediction of critical ignitions and of the turbulent flame propagation.In the first part, DNS of turbulent ignition were performed. The ignition phase was modeled using a thermal energy deposit (ED model, Lacaze et al.). Simulations reproduced the ignition experiments of Cardin et al. who determined the minimum ignition energy (MIE) of lean premixed methane/air mixtures, for different turbulence characteristics. The main purpose of the study was to determine the numerical and physical model parameters, which enable to reproduce Cardin et al. experiments. Two types of kinetic schemes were evaluated: a simplified kinetic scheme and an analytical kinetic scheme (ARC), that can predict both the auto-ignition delays and the laminar flame speed, while keeping affordable CPU times. Results analysis enabled to define ignition criteria and to highlight the differences in terms of ignition prediction using the two kinetic schemes. Results also demonstrated that the chosen approach could recover correct levels of ignition energy for laminar and low Karlovitz number cases (Ka<10). For higher Karlovitz number cases, the ED model was found to be insufficient to predict the ignition and a finer description of the energy deposit is required.In the second part, a dynamic wrinkling model (Wang et al., 2012) was studied to describe the out-of-equilibrium behavior of the flame during the propagation phase. Studies on laminar spherical flames were first performed, to assess the laminar degeneration of the model. Then, as first tests in an engine configuration have revealed incompatibilities of the model, modifications were proposed. The modified dynamic model was finally tested in the ICAMDAC engine configuration. Results of the simulations were compared against previous results of Robert et al. obtained with the ECFM-LES model using a transport equation for the flame surface density that can describe the out-of-equilibrium wrinkling of the flame. Results obtained with the dynamic model are in very good agreement with the ones of Robert et al., thus demonstrating the ability of the dynamic model to predict out-of-equilibrium values in the engine configuration. Besides, the dynamic model self-adapts to the turbulence conditions, hence does not require any model parameter adjustment, as is it the case for models based on the flame surface density transport equation.

Combustion has been used for a long time as a means of energy extraction. However, in the recent years there has been further increase in air pollution, through pollutants such as nitrogen oxides, acid rain etc. To solve this problem, there is a need to reduce carbon and nitrogen oxides through lean burning, fuel dilution and usage of bi-product fuel gases. A numerical analysis has been carried out to investigate the effectiveness of several reduced mechanisms, in terms of computational time and accuracy. The cases were tested for the combustion of hydrocarbons diluted with hydrogen, syngas, and bi-product fuel in a cylindrical combustor. The simulations were carried out using the ANSYS Fluent 19.1. By solving the conservations equations, several global reduced mechanisms (2-5-10 steps) were obtained. The reduced mechanisms were used in the simulations for a 2D cylindrical tube with dimensions of 40 cm in length and 2.0 cm diameter. The mesh of the model included a proper fine quad mesh, within the first 7 cm of the tube and around the walls. By developing a proper boundary layer, several simulations were performed on hydrocarbon/air and syngas blends to visualize the flame characteristics. To validate the results “PREMIX and CHEMKIN” codes were used to calculate 1D premixed flame based on the temperature, composition of burned and unburned gas mixtures. Numerical calculations were carried for several hydrocarbons by changing the equivalence ratios (lean to rich) and adding small amounts of hydrogen into the fuel blends. The changes in temperature, radical formation, burning velocities and the reduction in NOx and CO2 emissions were observed. The results compared to experimental data to study the changes. Once the results were within acceptable range, different fuels compositions were used for the premixed combustion through adding H2/CO/CO2 by volume and changing the equivalence ratios and preheat temperatures, in the fuel blends. The results on flame temperature, shape, burning velocity and concentrations of radicals and emissions were observed. The flame speed was calculated by finding the surface area of the flame, through the mass fractions of fuel components and products conversions that were simulated through the tube. The area method was applied to determine the flame speed. It was determined that the reduced mechanisms provided results within an acceptable range. The variation of the inlet velocity had neglectable effects on the burning velocity. The highest temperatures were obtained in lean conditions (0.5-0.9) equivalence ratio and highest flame speed was obtained for Blast Furnace Gas (BFG) at elevated preheat temperature and methane-hydrogen fuels blends in the combustor. The results included; reduction in CO2 and NOx emissions, expansion of the flammable limit, under the condition of having the same laminar flow. The usage of diluted natural gases, syngas and bi-product gases provides a step in solving environmental problems and providing efficient energy.

No description available.

No description available.

No description available.

No description available.

No description available.

No description available.

No description available.