The propagation of premixed flames in turbulent flows is a problem of wide physical and technological interest, with a significant literature on their propagation speed and front topology. While certain scalings and parametric dependencies are well understood, a variety of problems remain. One major challenge, and focus of this thesis, is to model the influence of fuel/oxidizer composition on turbulent burning rates. Classical explanations for augmentation of turbulent burning rates by turbulent velocity fluctuations rely on global arguments - i.e., the turbulent burning velocity increase is directly proportional to the increase in flame surface area and mean local burning rate along the flame. However, the development of such global approaches is complicated by the abundance of phenomena influencing the propagation of turbulent premixed flames. Emphasizing key governing processes and cutting-off interesting but marginal phenomena appears to be necessary to make further progress in understanding the subject. An alternative approach to understand turbulent augmentation of burning rates is based upon so-called "leading points", which are intrinsically local properties of the turbulent flame. Leading points concepts suggest that the key physical mechanism controlling turbulent burning velocities of premixed flames is the velocity of the points on the flame that propagate farthest out into the reactants. It is postulated that modifications in the overall turbulent combustion speed depend solely on modifications of the burning rate at the leading points since an increase (decrease) in the average propagation speed of these points causes more (less) flame area to be produced behind them. In this framework, modeling of turbulent burning rates can be thought as consisting of two sub-problems: the modeling of (1) burning rates at the leading points and of (2) the dynamics/statistics of the leading points in the turbulent flame. The main objective of this thesis is to critically address both aspects, providing validation and development of the physical description put forward by leading point concepts. To address the first sub-problem, a comparison between numerical simulations of one-dimensional laminar flames in different geometrical configurations and statistics from a database of direct numerical simulations (DNS) is detailed. In this thesis, it is shown that the leading portions of the turbulent flame front display a structure that on average can be reproduced reasonably well by results obtained from model geometries with the same curvature. However, the comparison between model laminar flame computations and highly curved flamelets is complicated by the presence of negative (i.e., compressive) strain rates, due to gas expansion. For the highest turbulent intensity investigated, local consumption speeds, curvatures, strain rates and flame thicknesses approach the maximum values obtained by the laminar model geometries, while other cases display substantially lower values. To address the second sub-problem, the dynamics of flame propagation in simplified flow geometries is studied theoretically. Utilizing results for Hamilton-Jacobi equations from the Aubry-Mather theory, it is shown how the overall flame front progation under certain conditions is controlled only by discrete points on the flame. Based on these results, definitions of leading points are proposed and their dynamics is studied. These results validate some basic ideas from leading points arguments, but also modify them appreciably. For the simple case of a front propagating in a one-dimensional shear flow, these results clearly show that the front displacement speed is controlled by velocity field characteristics at discrete points on the flame only when the amplitude of the shear flow is sufficiently large and does not vary too rapidly in time. However, these points do not generally lie on the farthest forward point of the front. On the contrary, for sufficiently weak or unsteady flow perturbations, the front displacement speed is not controlled by discrete points, but rather by the entire spatial distribution of the velocity field. For these conditions, the leading points do not have any dynamical significance in controlling the front displacement speed. Finally, these results clearly show that the effects of flame curvature sensitivity in modifying the front displacement speed can be successfully interpreted in term of leading point concepts.

Turbulent premixed propane/air and methane/air flames were studied using planar Rayleigh scattering and particle image velocimetry on a stabilized Bunsen type burner. The fuel-air equivalence ratio was varied from Φ=0.7 to 1.0 for propane flames, and from Φ=0.6 to 1.0 for methane flames. The non-dimensional turbulence intensity, u'/SL (ratio of fluctuation velocity to laminar burning velocity), covered the range from 3 to 24, equivalent to conditions of corrugated flamelets and thin reaction zones regimes. Temperature gradients decreased with the increasing u'/SL and levelled off beyond u'/SL > 10 for both propane and methane flames. Flame front thickness increased slightly as u'/SL increased for both mixtures, although the thickness increase was more noticeable for propane flames, which meant the thermal flame front structure was being thickened. A zone of higher temperature was observed on the average temperature profile in the preheat zone of the flame front as well as some instantaneous temperature profiles at the highest u'/SL. Curvature probability density functions were similar to the Gaussian distribution at all u'/SL for both mixtures and for all the flame sections. The mean curvature values decreased as a function of u'/SL and approached zero. Flame front thickness was smaller when evaluated at flame front locations with zero curvature than that with curvature. Temperature gradients and FSD were larger when the flame curvature was zero. The combined thickness and FSD data suggest that the curvature effect is more dominant than that of the stretch by turbulent eddies during flame propagation. Integrated flame surface density for both propane and methane flames exhibited no dependance on u'/SL regardless of the FSD method used for evaluation. This observation implies that flame surface area may not be the dominant factor in increasing the turbulent burning velocity and the flamelet assumption may not be valid under the conditions studied. Dκ term, the product of diffusivity evaluated at conditions studied and the flame front curvature, was a magnitude smaller than or the same magnitude as the laminar burning velocity.

Turbulent premixed propane/air and methane/air flames were studied using planar Rayleigh scattering and particle image velocimetry on a stabilized Bunsen type burner. The fuel-air equivalence ratio was varied from Φ=0.7 to 1.0 for propane flames, and from Φ=0.6 to 1.0 for methane flames. The non-dimensional turbulence intensity, u'/SL (ratio of fluctuation velocity to laminar burning velocity), covered the range from 3 to 24, equivalent to conditions of corrugated flamelets and thin reaction zones regimes. Temperature gradients decreased with the increasing u'/SL and levelled off beyond u'/SL > 10 for both propane and methane flames. Flame front thickness increased slightly as u'/SL increased for both mixtures, although the thickness increase was more noticeable for propane flames, which meant the thermal flame front structure was being thickened. A zone of higher temperature was observed on the average temperature profile in the preheat zone of the flame front as well as some instantaneous temperature profiles at the highest u'/SL. Curvature probability density functions were similar to the Gaussian distribution at all u'/SL for both mixtures and for all the flame sections. The mean curvature values decreased as a function of u'/SL and approached zero. Flame front thickness was smaller when evaluated at flame front locations with zero curvature than that with curvature. Temperature gradients and FSD were larger when the flame curvature was zero. The combined thickness and FSD data suggest that the curvature effect is more dominant than that of the stretch by turbulent eddies during flame propagation. Integrated flame surface density for both propane and methane flames exhibited no dependance on u'/SL regardless of the FSD method used for evaluation. This observation implies that flame surface area may not be the dominant factor in increasing the turbulent burning velocity and the flamelet assumption may not be valid under the conditions studied. Dκ term, the product of diffusivity evaluated at conditions studied and the flame front curvature, was a magnitude smaller than or the same magnitude as the laminar burning velocity.

The persistent thrust for a cleaner, greener environment has prompted air pollution regulations to be enforced with increased stringency by environmental protection bodies all over the world. This has prompted gas turbine manufacturers to move from non-premixed combustion to lean, premixed combustion. These lean premixed combustors operate quite fuel-lean compared to the stochiometric, in order to minimize CO and NOx productions, and are very susceptible to oscillations in any of the upstream flow variables. These oscillations cause the heat release rate of the flame to oscillate, which can engage one or more acoustic modes of the combustor or gas turbine components, and under certain conditions, lead to limit cycle oscillations. This phenomenon, called thermoacoustic instabilities, is characterized by very high pressure oscillations and increased heat fluxes at system walls, and can cause significant problems in the routine operability of these combustors, not to mention the occasional hardware damages that could occur, all of which cumulatively cost several millions of dollars. In a bid towards understanding this flow-flame interaction, this research works studies the heat release response of premixed flames to oscillations in reactant equivalence ratio, reactant velocity and pressure, under conditions where the flame preheat zone is convectively compact to these disturbances, using the G-equation. The heat release response is quantified by means of the flame transfer function and together with combustor acoustics, forms a critical component of the analytical models that can predict combustor dynamics. To this end, low excitation amplitude (linear) and high excitation amplitude (nonlinear) responses of the flame are studied in this work. The linear heat release response of lean, premixed flames are seen to be dominated by responses to velocity and equivalence ratio fluctuations at low frequencies, and to pressure fluctuations at high frequencies which are in the vicinity of typical screech frequencies in gas turbine combustors. The nonlinear response problem is exclusively studied in the case of equivalence ratio coupling. Various nonlinearity mechanisms are identified, amongst which the crossover mechanisms, viz., stoichiometric and flammability crossovers, are seen to be responsible in causing saturation in the overall heat release magnitude of the flame. The response physics remain the same across various preheat temperatures and reactant pressures. Finally, comparisons between the chemiluminescence transfer function obtained experimentally and the heat release transfer functions obtained from the reduced order model (ROM) are performed for lean, CH4/Air swirl-stabilized, axisymmetric V-flames. While the comparison between the phases of the experimental and theoretical transfer functions are encouraging, their magnitudes show disagreement at lower Strouhal number gains show disagreement.

This dissertation introduces a general, predictive and cost-efficient reduced-order modeling (ROM) technique for characterization of flame response under acoustic modulation. The model is built upon the kinematic flame model–G-equation to describe the flame topology and dynamics, and the novelties of the ROM lie in i) a procedure to create the compatible base flow that can reproduce the correct flame geometry and ii) the use of a physically-consistent acoustic modulation field for the characterization of flame response. This ROM addresses the significant limitations of the classical kinematic model, which is only applicable to simple flame configurations and relies on ad-hoc models for the modulation field. The ROM is validated by considering the acoustically-excited premixed methane/air flames in conical and M-shape configurations. To test the model availability to practical burners, a confined flame configuration is also employed for model evaluation. Furthermore, to investigate the generality of the ROM to the burner flame, the performance of the ROM with respect to the V-shape and the swirled V-shape is investigated. The model accuracy is evaluated concerning flame geometrical features and flame describing function, and assessed by comparing the ROM results with both experimental measurements and direct- numerical-simulation results. It is found that the flame describing/transfer functions predicted by the ROM compare well with reference data, and are more accurate than those obtained from the conventional kinematic model built upon heuristically-presumed modulation fields.

Utsläpp från bilar har spelat stor roll de senaste decennierna. Detta har lett till ökad användning av Otto gasmotorer som använder naturgas som bränsle. Nya motordesigner behöver optimeras för att förbättra motorens effektivitet. Ett effektivt sätt att göra detta på är genom användningen av simuleringar för att minska ledtiden i motorutvecklingen. Verifiering och validering av simuleringarna spelar stor roll för att bygga förtroende för och förutsägbarhet hos simuleringsresultaten. Syftet med detta examensarbete är att föreslå förbränningsmodellparametrarna efter utvärdering av olika kombinationer av förbrännings- och tändmodeller för Otto förbränning, vad gäller beräkningstid och noggrannhet. In-cylindertrycksspår från simulering och mätning jämförs för att hitta den bästa kombinationen av förbrännings- och tändmodell. Inverkan av tändtid, antal motorcykler och randvillkor för simuleringsresultatet studeras också. Resultaten visar att ECFM-förbränningsmodellen förutsäger simuleringsresultaten mer exakt när man jämför med mätningarna. Effekten av tändningstiden på olika kombinationer av förbrännings- och tändningsmodell utvärderas också. Stabiliteten hos olika förbränningssimuleringsmodeller diskuteras också under körning för fler motorcykler. Jämförelse av beräkningstid görs även för olika kombinationer av förbrännings- och tändmodeller. Resultaten visar också att flamspårningsmetoden med Euler är mer känslig för cellstorlek och kvalitet hos simuleringsnätet, jämfört med övriga studerade modeller. Rekommendationer och förslag ges om nät- och simulerings-inställningar för att prediktera förbränningen på ett så bra sätt som möjligt. Några möjliga förbättringsområden ges som framtida arbete för att förbättra noggrannheten i simuleringsresultaten. / Emission from automobiles has been gaining importance for past few decades. This has gained a lot of impetus in search for alternate fuels among the automotive manufacturers. This led to the increase usage of Otto gas engine which uses natural gas as fuel. New engine designs have to be optimized for improving the engine efficiency. This led to usage of virtual simulations for reducing the lead time in the engine development. The verification and validation of actual phenomenon in the virtual simulations with respect to the physical measurements was quite important.  The aim of this master thesis is to suggest the combustion model parameters after evaluating various combination of combustion and ignition models in terms of computational time and accuracy. In-cylinder pressure trace from the simulation is compared with the measurement in order to find the nest suited combination of combustion and ignition models. The influence of ignition timing, number of engine cycles and boundary conditions on the simulation results are also studied. Results showed that ECFM combustion model predicts the simulation results more accurately when compare to the measurements. Impact of ignition timing on various combination of combustion and ignition model is also assessed. Stability of various combustion simulation models is also discussed while running for more engine cycles. Comparison of computational time is also made for various combination of combustion and ignition models. Results also showed that the flame tracking method using Euler is dependent on the mesh resolution and the mesh quality.  Recommendations and suggestions are given about the mesh and simulation settings for predicting the combustion simulation accurately. Some possible areas of improvement are given as future work for improving the accuracy of the simulation results.