Lean premixed combustion is at present one of the most promising methods to reduce emissions and to maintain high efficiency in combustion systems. As the emission legislation becomes more stringent, modelling of turbulent premixed combustion has become an important tool for designing efficient and environmentally friendlier combustion systems. However, in order to predict these emissions reliable predictive models are required. One of the methods used for predicting pollutants is the conditional moment closure (CMC), which is suitable to predict pollutants with slow time scales. Despite the fact that CMC has been successfully applied to various non-premixed combustion systems, its application to premixed flames is not fully tested and validated. The main difficulty is associated with the modelling of the conditional scalar dissipation rate (CSDR) of the conditioning scalar, the progress variable. In premixed CMC, this term is an important quantity and represents the rate of mixing at small scales of relevance for combustion. The numerical accuracy of the CMC method depends on the accuracy of the CSDR model. In this study, two different models for CSDR, an algebraic model and an inverse problem model, are validated using two different DNS data sets. The algebraic model along with standard k-ε turbulence modelling is used in the computations of stoichiometric and very lean pilot stabilized Bunsen flames using the RANS-CMC method. A first order closure is used for the conditional mean reaction rate. The computed nonreacting and reacting scalars are in reasonable agreement with the experiments and are consistent with earlier computations using flamlets and transported PDF methods for the stoichiometric flames, and transported PDF methods for the very lean flames. Sensitivity to chemical kinetics mechanism is also assessed.

<p> In premixed turbulent combustion, the normal speed of propagation of the flame front is enhanced by the turbulent velocity field. This project will focus on the method of computing the normal speed of propagation of the flame front in the Majda-Souganidis model of turbulent combustion. Solving this problem involves computing the eigenvalue of a nonlinear cell problem. Discussed in this thesis is a new, simple and direct numerical method for approximating the eigenvalue, also called the effective Hamiltonian.</p> / Thesis / Master of Science (MSc)

Producer gas derived from biomass is one of the most environment friendly substitutes to the fossil fuels. Usage of producer gas for power generation has effect of zero net addition of CO2 in atmosphere. The engines working on producer gas have potential to decrease the dependence on conventional fuels for power generation. However, the combustion process is governed by complex interactions between chemistry and fluid dynamics, some of which are not completely understood. Improved knowledge of combustion is, therefore, of vital importance for both direct use in the design of engines, and for the evolution of reliable simulation tools for engine development. The present work is related to the turbulent combustion of producer gas in closed vessels and engine cylinders. The main objective of the work was multi-dimensional simulation of turbulent combustion in the bowl-in-piston engine operating on producer gas fuel and to observe the flame and flow field interaction. First, the combustion model was validated in constant volume combustion chamber with experimental results. Experimental turbulent combustion data of producer gas (composition matching with engine operating conditions) was presented. The required data of laminar burning velocity of producer gas was computed and used in the simulation of turbulent combustion in closed vessel. The effect of squish and reverse squish flow on flame propagation in the bowl-in-piston engine cylinder was described. Laminar burning velocity of unstretched flame was computed using flame code which was developed earlier in this laboratory. One dimensional computations of unstretched planar flame were made to calculate laminar burning velocity of the producer gas-air mixture at pressures (1-10 bar) and temperatures (300-600 K). A correlation of laminar burning velocity of producer gas as a function of pressure and temperature was fitted and compared with experiments. A fixed composition and equivalence ratio of producer gas-air mixture, typical of the engine operating conditions, was considered. The correlation was used in simulation of turbulent combustion in closed vessel. The turbulent combustion experiments with producer gas-air mixture were conducted in a closed vessel. The aim of experiments was to generate pressure-time data, in closed vessel during turbulent flame propagation, which was required to validate turbulent combustion models. Determination of (ST /SL) was made from pressure-time data which requires corresponding laminar combustion data with same initial conditions. For this purpose a set of laminar combustion experiments was conducted. Experimental setup consists of a constant volume combustion chamber of cubical shape and size 80 x 80 x 80 mm3 . The initial mixtures pressure and temperature were 1 bar and 300 K respectively. A fixed composition and equivalence ratio of producer gas-air mixture, typical of the engine operating conditions, was used. The composition of producer gas was H2 -19.61%, CO2 -19.68%, CH4 -2.52%, CO2 -12.55% and N2 -45.64% on volume basis. Fuel-air mixture was ignited with electric spark at the center of the cube. Initial turbulence in the chamber was created by moving a perforated plate with specified velocity. Perforated plate was placed in chamber so that the central hole in the plate passes over the spark electrodes as it sweeps across the chamber. Two geometrically similar plates with hole diameter of 5 and 10 mm were used. The new experimental setup constructed as a part of this work was first tested with one set of experiments each with methane and propane data of SL and ST /SL from the literature. Maximum turbulent intensity (u’) achieved was 1.092 ms−1 . The ratios of turbulent to laminar burning velocity (ST /SL) values were determined at six different turbulence intensity levels. Laminar combustion experiments were extended to elevated initial pressures 2-5 bar and temperature 300 K. The value of SL was calculated from the pressure-time history recorded during laminar stretching flame propagation inside closed vessel. These SL values were compared with computed SL,∞ after accounting for stretch. Turbulent combustion simulations were carried out to validate combustion models suitable for multi-dimensional CFD simulation of combustion in constant volume closed chamber. Two models proposed by Choi and Huh, based on Flame Surface Density (FSD) were tested with the present experimental results. User FORTRAN code for the source terms in transport equation of FSD was implemented in ANSYS-CFX 10.0 software. First model called CFM1, grossly under-predicted the rate of combustion. The second model called CFM2, predicted the results satisfactorily after replacing the arbitrary length scale with turbulent integral length scale (lt) having a limiting value near the wall. The modified CFM2 model was able to predict the propagation phase of the developed flame satisfactorily, though the duration for initial flame development was over-predicted by the model. CFD simulation of producer gas engine combustion process was carried out using ANSYSCFX software. Mesh deformation option was used to take care of moving boundaries such as piston and valve surfaces. The fluid domain expands during suction process and contracts during compression process. In order to avoid excessive distortion of the mesh elements, a series of meshes at different crank angle positions were generated and checked for their quality during mesh motion in the solver. For suction process simulation, unstructured meshes having 0.1 to 0.3 million cells were used. During the compression and combustion process simulations, structured meshes having 40,000 to 0.1 million cells were used. k-ε model was used for turbulence simulation. The suction, compression and combustion processes of an SI engine were simulated. Initial flame kernel was given by providing high flame surface density in a small volume comparable to the spark size at the time of ignition. The flame surface density model, CFM-2, was adapted with the modification of length scale tested against constant volume experiments. A suitable limiting value was used to avoid abnormal flame propagation near the wall. The limiting value of integral length scale (lt) near the wall was determined by linear extrapolation of the integral length scale in the domain to the wall. Engine p - θ curves of three different ignition timings 26°, 12° and 6°before top dead center (TDC) were simulated and compared with earlier experimental results. The effects of flow field on flame propagation have been observed. A comparison of the simulated and experimental p - θ diagram of the engine for all above cases gave mixed results. For the ignition timing at 26° before TDC case, predicted peak pressure value was 17% higher and at 3° earlier than those of the experimental peak. For the other two cases, the predicted peak pressure value was 28% lower and 5° later than those of the experimental peak. The reason for under-prediction of the pressure values could be due to the delay in development of initial flame kernel. Simulated pressure curves have offset about 3-4° compared to the experimental pressure curves. It was observed that in all predicted p - θ cases, there was a delay in the initial flame development. It is evident from the under-prediction of pressure values, especially in the initial flame kernel development phase and it also affects the p - θ curve at later stage. The delay was about 3-4° of crank angle rotation in various cases. The delay in predicting the initial flame development needs to be corrected in order to predict the combustion process properly. The proposed FSD model seems to have capability to predict p - θ values fairly in the propagation phase of developed flame. Reasonably good match was obtained by advancing the ignition timing in the computation by about 3-4° compared to the experimental setting. In the bowl-in-piston engine cylinders, the flow in the cylinder is characterised by squish and reverse squish when the piston is moving towards and away from the top dead center (TDC) respectively. The effect of squish and reverse squish flow on flame propagation has been assessed. For the more advanced ignition case, i.e., 26° before TDC, The flame propagation did not have favorable effect by the flow field. The direction of flame propagation was against the squish and reverse squish flow. This resulted in suppressed peak velocities in the cylinder compared the motoring process. Hence the burning rate was not augmented by the turbulence inside the cylinder. For the ignition 12° before TDC case, the flame propagation did have favorable effect by the flow field. During the reverse squish period, the flame had reached the bowl wall. At this stage, the flame was pushing the reactants out and this augments the reverse-squish flow, and hence the maximum reverse-squish velocity was increased to 2.03 times the peak reverse-squish velocity of motoring case. The reverse-squish flow was distorting the flame from spherical shape and the flame gets stretched. Flame surface enters the cylindrical region faster compared to the previous case. The stretched flame in the reverse-squish flow may be considered as reverse squish flame, as was proposed earlier by Sridhar G. The burn rate during the reverse squish period may be 2 to 2.5 times the normal burn rate. For the ignition 6° before TDC case, the flame was very small in size and it did not affect the flow in squish period. During the reverse squish period, the flame radius was moderate compared to the bowl radius. The flame was pushing the reactants out and it increased the maximum reverse-squish velocity to 1.3 times by the flame. In this case, the reverse-squish flow moderately affecting the flame shapes. The results of this study could give an idea of what ignition timing must be kept for favorable use of flow field inside the engine cylinder. Main contributions from the present work are: Multi-dimensional simulation of combustion process inside the engine cylinder operating on producer gas was carried out to examine flame/flow field interactions. Two models based on FSD were first tested against present experimental results in constant volume combustion chamber. In CFM2 model; a modification of replacing the arbitrary length scale by integral length scale with a limiting value near the wall was suggested to avoid prediction of abnormally large turbulent burning velocity near the wall. This combustion model has been implemented in ANSYS-CFX10. The required data of laminar and turbulent burning velocities of producer gas-air mixture has been determined by experiments and computations at varied initial pressures and turbulent intensities. Finally, the simulated engine pressure data has been compared with earlier experimental data of the engine operating on producer gas. The proposed FSD model has the capability to match well with the experimental results except for the initial flame kernel development phase. Even though this issue needs to be resolved, the work has brought out the important interaction between the flame propagation and flow field within the bowl-in-piston engine cylinder.

In the last 10-15 years, large-eddy simulation (LES) has become well established for non-reacting flows, and several successful models have been developed for the transfer of momentum and kinetic energy to the subfilter-scales (SFS). However, for reacting flows, LES is still undergoing significant development. In particular, for many premixed combustion applications, the chemical reactions are confined to propagating surfaces that are significantly thinner than the computational grids used in practical LES. In these situations, the chemical kinetics and its interaction with the turbulence are not resolved and must be entirely modelled. There is, therefore, a need for accurate and robust physical modelling of combustion at the subfilter-scales. In this thesis, modelled transport equations for progress variable and flame surface density (FSD) were implemented and coupled to the Favre-filtered Navier-Stokes equations for a compressible reactive thermally perfect mixture. In order to reduce the computational costs and increase the resolution of simulating combusting flows, a parallel adaptive mesh (AMR) refinement finite-volume algorithm was extended and used for the prediction of turbulent premixed flames. The proposed LES methodology was applied to the numerical solution of freely propagating flames in decaying isotropic turbulent flow and Bunsen-type flames. Results for both stoichiometric and lean flames are presented. Comparisons are made between turbulent flame structure predictions for methane, propane, hydrogen fuels, and other available numerical results and experimental data. Details of subfilter-scale modelling, numerical solution scheme, computational results, and capabilities of the methodology for predicting premixed combustion processes are included in the discussions. The current study represents the first application of a full transport equation model for the FSD to LES of a laboratory-scale turbulent premixed flame. The comparisons of the LES results of this thesis to the experimental data provide strong support for the validity of the modelled transport equation for the FSD. While the LES predictions of turbulent burning rate are seemingly correct for flames lying within the wrinkled and corrugated flamelet regimes and for lower turbulence intensities, the findings cast doubt on the validity of the flamelet approximation for flames within the thin reaction zones regime.

Avec l’accroissement considérable de la puissance de calcul, les simulations aux grandes échelles (SGE) sont maintenant utilisées de façon routinière dans de nombreuses applications d’ingénierie. Les modèles de combustion usuels utilisés dans les SGE sont le plus souvent basés sur une hypothèse d’équilibre entre le mouvement des structures turbulentes et le plissement de la surface de la flamme. Ils s’écrivent alors sous forme d’expressions algébriques fonctions de grandeurs connues aux échelles résolues ainsi que de paramètres dont l’ajustement est à la charge de l’utilisateur selon la configuration étudiée et les conditions opératoires. Le modèle dynamique récemment développé ajuste automatiquement au cours du calcul les paramètres de modélisation qui peuvent alors dépendre du temps et de l’espace. Cette thèse présente une étude détaillée d’un modèle dynamique pour la simulation aux grandes échelles de la combustion turbulente prémélangée. L’objectif est de caractériser, explorer les avantages et les inconvénients, appliquer et valider le modèle dynamique dans plusieurs configurations. / Large eddy simulation (LES) is currently applied in a wide range of engineering applications. Classical LES combustion models are based on algebraic expressions and assume equilibrium between turbulence and flame wrinkling which is generally not verified in many circumstances as the flame is laminar at early stages and progressively wrinkled by turbulent motions. In practice, this conceptual drawback has a strong consequence: every computation needs its own set of constants, i.e. any small change in the operating conditions or in the geometry requires an adjustment of model parameters. The dynamic model recently developed adjust automatically the flame wrinkling factor from the knowledge of resolved scales. Widely used to describe the unresolved turbulent transport, the dynamic approach remains underexplored in combustion despite its interesting potential. This thesis presents a detailed study of a dynamic wrinkling factor model for large eddy simulation of turbulent premixed combustion. The goal of this thesis is to characterize, unveil pros and cons, apply and validate the dynamic modeling in different flow configurations.

L’industrie automobile se trouve aucentre des préoccupations environnementalesactuelles. Les moteurs essence à injection directeopérés en condition pauvres offrent un fortpotentiel en terme de réduction des émissions depolluants. En contrepartie, ils sont sujets à uneforte variabilité cyclique de combustion (CCV)qui ne peut être que partiellement étudiéeexpérimentalement. La simulation aux grandeséchelles (SGE) apparait comme une approchenumérique adaptée pour étudier de telsphénomènes, du fait de sa capacité naturelle àcapter les phénomènes instationnaires. Laprésente thèse se propose d’une part d’estimer lacapacité de la SGE à reproduire les CCVobservées expérimentalement, et d’autre part decontribuer à une meilleure compréhension deleur apparition. Dans ce contexte, un effortparticulier est mis sur la modélisation desphénomènes proche paroi. En particulier, unmodèle de paroi adapté à l’étude del’aérodynamique interne de configurationindustrielles est proposé. Il est validé sur desconfigurations de complexités variées. Enfin,l’étude porte sur le cas du moteur M256 qui estétudié en s’appuyant sur une solide base dedonnées expérimentales. Les causes de CCV sontexplorées, et notamment, l’impact de lavariabilité de l’écoulement généré pendant laphase d’admission sur la propagation du front deflamme est clarifié. / The automotive industry finds itselfat the center of current environmental concerns.Modern direct injection engines, operated underlean condition have the potential to reducepollutant emissions. As a drawback, they aresubject to large cyclic combustion variability(CCV), that can be explained only partially byexperimental measurements. Large-EddySimulation (LES) appears as an adapted tool tocomplement experiments, due to its naturalability to capture unsteady phenomena. Thepresent PhD thesis first aims at reproducing theCCV, and at contributing achieving a betterunderstanding of their occurrence.In this context, a special effort is put on thereproduction of near-wall phenomena, throughthe proposal of a new wall boundary conditionthat is validated on cases of differentcomplexity. Then the focus is put on the M256case, for which an extensive experimentaldatabase is available. The causes of CCV areexplored, and in particular, the impact of thevariability of the intake flow on the flame frontpropagation is clarified.

In the present work, the processes of steady combustion and autoignition of hydrogen are investigated using the Conditional Moment Closure (CMC) model with a Reynolds Averaged Navier-Stokes (RANS) Computational Fluid Dynamics (CFD) code. A study of the effects on the flowfield of changing turbulence model constants, specifically the turbulent Schmidt number, Sct, and C epsilon 1 of the k − epsilon model, are investigated. The effects of two different mixing models are explored: the AMC model, which is commonly used in CMC implementations, and a model based on the assumption of inhomogeneous turbulence. The background equations required for implementation of the CMC model are presented, and all relevant closures are discussed. The numerical implementation of the CMC model, in addition to other techniques aimed at reducing computational expense of the CMC calculations, are provided. The CMC equation is discretised using finite volume (FV) method. The CFD and CMC calculations are fully coupled, allowing for simulations of steady flames or flame development after the occurrence of autoignition. Through testing of a steady jet flame, it is observed that the flowfield calculations follow typical k − epsilon model trends, with an overprediction of spreading and an underprediction of penetration. The CMC calculations are observed to perform well, providing good agreement with experimental measurements. Autoignition simulations are conducted for 3 different cases of turbulence constants and 7 different coflow temperatures to determine the final effect on the steady flowfield. In comparison to the standard constants, reduction of Sct results in a reduction of the centreline mixing intensity within the flowfield and a corresponding reduction of ignition length, while reducing C 1 results in an increase of centreline mixing intensity and an increase in the ignition length. All scenarios tested result in an underprediction of ignition length in comparison to experimental results; however, good agreement with the experimental trends is achieved. At low coflow temperatures, the effects of mixing intensity within the flowfield are seen to have the largest influence on ignition length, while at high coflow temperatures, the chemical source term in the CMC equation increases in magnitude, resulting in very little difference between predictions for different sets of turbulence constants. The inhomogeneous mixing model is compared using the standard turbulence constants. A reduction of ignition lengths in comparison to the AMC model is observed. In steady state simulation of the autoigniting flow, the inhomogeneous model is observed to predict both lifted flames and fully anchored flames, depending on coflow temperature.

This work is intended to compare Large Eddy Simulation and Implicit Large Eddy Simulation (LES and ILES) for a turbulent, non-reacting or reacting high speed H2 jet in co-flowing air, typical of scramjet engines. Numerical simulations are performed at resolutions ranging from 32&times / 32&times / 128 to 256&times / 256&times / 1024, using a 5th order WENO scheme. Physical LES are carried out with the Smagorinsky and the Selective Structure Function models associated to molecular diffusion. Implicit LES are performed with and without molecular diffusion, by solving either the Navier-Stokes or the Euler equations. In the nonreacting case, the Smagorinsky model is too dissipative. The Selective Structure Function leads to better results, but does not show any superiority compared to ILES, whatever the grid resolution. In the reacting case, a molecular viscous cut-off in the simulation is mandatory to set a physical width for the reaction zone in the ILES approach, hence to achieve grid-convergence. It is also found that ILES/LES are less sensitive to the inlet conditions than the RANS approach. The first chapter is an introduction to the context of this study. In the second chapter, the governing equations for multispecies reacting flows are presented, with emphasis on the thermodynamic and transport models. In the third chapter, physical LES equations and explicit sub-grid modeling strategies iv are detailed. Some properties of the numerical scheme are also investigated. In chapter four, the numerical scheme and some aspects of the solver are explained. Finally, non-reacting and reacting numerical experiments are presented and the results are discussed.

The staged combustion aft-injected hybrid (SCAIH) rocket motor is a promising design for the future of hybrid rocket propulsion. Advances in computational fluid dynamics and scientific computing have made computational modelling an effective tool in design and development. The focus of this thesis is the numerical modelling of the SCAIH rocket motor in a turbulent combustion, high-speed, reactive flow accounting for solid soot transport and radiative heat transfer. The SCAIH motor has a shear coaxial injector with liquid oxygen injected centrally at sub-critical conditions: 150K, 150m/s (Mach≈0.9), and a gas-generator gas-solid mixture of one-third carbon soot by mass injected in the annual opening at 1175K, and 460m/s (Mach≈0.6). Flow conditions in the near injector region and the flame anchoring mechanism are of particular interest. Overall, the flow is shown to exhibit instabilities and the flame is shown to anchor directly on the injector faceplate with temperatures in excess of 2700K.

In the last 10-15 years, large-eddy simulation (LES) has become well established for non-reacting flows, and several successful models have been developed for the transfer of momentum and kinetic energy to the subfilter-scales (SFS). However, for reacting flows, LES is still undergoing significant development. In particular, for many premixed combustion applications, the chemical reactions are confined to propagating surfaces that are significantly thinner than the computational grids used in practical LES. In these situations, the chemical kinetics and its interaction with the turbulence are not resolved and must be entirely modelled. There is, therefore, a need for accurate and robust physical modelling of combustion at the subfilter-scales. In this thesis, modelled transport equations for progress variable and flame surface density (FSD) were implemented and coupled to the Favre-filtered Navier-Stokes equations for a compressible reactive thermally perfect mixture. In order to reduce the computational costs and increase the resolution of simulating combusting flows, a parallel adaptive mesh (AMR) refinement finite-volume algorithm was extended and used for the prediction of turbulent premixed flames. The proposed LES methodology was applied to the numerical solution of freely propagating flames in decaying isotropic turbulent flow and Bunsen-type flames. Results for both stoichiometric and lean flames are presented. Comparisons are made between turbulent flame structure predictions for methane, propane, hydrogen fuels, and other available numerical results and experimental data. Details of subfilter-scale modelling, numerical solution scheme, computational results, and capabilities of the methodology for predicting premixed combustion processes are included in the discussions. The current study represents the first application of a full transport equation model for the FSD to LES of a laboratory-scale turbulent premixed flame. The comparisons of the LES results of this thesis to the experimental data provide strong support for the validity of the modelled transport equation for the FSD. While the LES predictions of turbulent burning rate are seemingly correct for flames lying within the wrinkled and corrugated flamelet regimes and for lower turbulence intensities, the findings cast doubt on the validity of the flamelet approximation for flames within the thin reaction zones regime.