Numerical Study of Turbulent Jet Ignition in a Lean Premixed Configuration

Direct numerical simulations (DNS) of a hot combustion product jet interacting with a lean premixed hydrogen-air coflow are conducted to fundamentally investigate turbulent jet ignition (TJI) in a three-dimensional configuration. TJI is an efficient method for initiating and controlling combustion in ultra-lean combustion systems. Fully compressible gas dynamics and species equations are solved with high order finite difference methods. The hydrogen-air reaction is simulated with a reliable detailed chemical kinetics mechanism. The physical processes involved in the TJI-assisted combustion are investigated by considering the flame heat release, temperature, species concentrations, vorticity, and Baroclinc torque. The complex turbulent flame and flow structures are delineated in three main: i) hot product jet, ii) burned-mixed, and iii) flame zones. In the TJI-assisted combustion, the flow structures and the flame features such as flame speed, temperature, and species distribution are found to be quite different than those in “standard” turbulent premixed combustion due to the existence of a high energy turbulent hot product jet. The flow structures and statistics are also found to be different than those normally seen in non-isothermal non-reacting jets.     

A pore scale study on turbulent combustion in porous media

This paper presents pore scale simulation of turbulent combustion of air/methane mixture in porous media to investigate the effects of multidimensionality and turbulence on the flame within the pores of porous media. In order to investigate combustion in the pores of porous medium, a simple but often used porous medium consisting of a staggered arrangement of square cylinders is considered in the present study. Results of turbulent kinetic energy, turbulent viscosity ratio, temperature, flame speed, convective heat transfer and thermal conductivity are presented and compared for laminar and turbulent simulations. It is shown that the turbulent kinetic energy increases from the inlet of burner, because of turbulence created by the solid matrix with a sudden jump or reduction at the flame front due to increase in temperature and velocity. Also, the pore scale simulation revealed that the laminarization of flow occurs after flame front in the combustion zone and turbulence effects are important mainly in the preheat zone. It is shown that turbulence enhances the diffusion processes in the preheat zone, but it is not enough to affect the maximum flame speed, temperature distribution and convective heat transfer in the porous burner. The dimensionless parameters associated with the Borghi–Peters diagram of turbulent combustion have been analyzed for the case of combustion in porous media and it is found that the combustion in the porous burner considered in the present study concerns the range of well stirred reactor very close to the laminar flame region.     

Leading-Edge Reaction Zones in Lifted-Jet Gas and Spray Flames

An investigation of the leading edge characteristics in lifted turbulent methane-air (gaseous) and ethanol-air (spray) diffusion flames is presented. Both combustion systems consist of a central nonpremixed fuel jet surrounded by low-speed air co-flow. Non-intrusive laser-based diagnostic techniques have been applied to each system to provide information regarding the behavior of the combustion structures and turbulent flow field in the regions of flame stabilization. Simultaneous sequential CH-PLIF/particle image velocimetry and CH-PLIF/Rayleigh scattering measurements are presented for the lifted gaseous flame. The CH-PLIF data for the lifted gas flame reveals the role that ``leading-edge' combustion plays as the stabilization mechanism in gaseous diffusion flames. This phenomenon, characterized by a fuel-lean premixed flame branch protruding radially outward at the flame base, permits partially premixed flame propagation against the incoming flow field. In contrast, the leading edge of the ethanol spray flame, examined using single-shot OH-PLIF imaging and smoke-based flow visualization, does not exhibit the same variety of leading-edge combustion structure, but instead develops a dual reaction zone structure as the liftoff height increases. This dual structure is a result of the partial evaporation (hence partial premixing) of the polydisperse spray and the enhanced rate of air entrainment with increased liftoff height (due to co-flow). The flame stabilizes in a region of the spray, near the edge, occupied by small fuel droplets and characterized by intense mixing due to the presence of turbulent structures. This revised version was published online in July 2006 with corrections to the Cover Date.     

Turbulence,vortex and external explosion induced by venting

The process of explosion venting to air in a cylindrical vent vessel connected to a duct, filling with a stoichiometric methane-oxygen gas mixture, was simulated numerically by using a colocated grid SIMPLE scheme based on k-epsilon turbulent model and Eddydissipation combustion model. The characteristics of the combustible cloud, flame and pressure distribution in the external flow field during venting were analyzed in terms of the predicted results. The results show that the external explosion is generated due to violent turbulent combustion in the high pressure region within the external combustible cloud ignited by a jet flame. And the turbulence and vortex in the external flow field were also discussed in detail. After the jet flame penetrating into the external combustible cloud, the turbulent intensity is greater in the regions with greater average kinetic energy gradient, rather than in the flame front ; and the vortex in the external flow field is generated primarily due to the baroclinic effect, which is greater in the regions where the pressure and density gradients are nearly perpendicular.     

Estimation of three-dimensional flame surface densities from planar images in turbulent premixed combustion

Turbulence motions are, by nature, three-dimensional while planar imaging techniques, widely used in turbulent combustion, give only access to two-dimensional information. For example, to extract flame surface densities, a key ingredient of some turbulent combustion models, from planar images implicitly assumes an instantaneously two-dimensional flow, neglecting the unresolved flame front wrinkling. The objective here is to estimate flame surface densities from two-dimensional measurements assuming that (1) the flow is statistically two dimensional; (2) the measuring plane is a plane of symmetry of the mean flow, either by translation (homogeneous third direction as in slot burners for example) or by rotation (axi-symmetrical flows such as jets) and (3) flame movements in transverse directions are similar. The unknown flame front wrinkling is then modelled from known quantities. An excellent agreement is achieved against direct numerical simulation (DNS) data where all three-dimensional quantities are known, but validations in other conditions (larger DNS, experiments) are required.     

Influence of flow field scaling on flashback of swirl flames

In this paper, the effect of geometrical scaling on the onset of flashback into a cylindrical premixing zone of a swirl flame is investigated. We discriminate two types of flashback. In the first type of flashback the flame propagates upstream inside an already present axial recirculation zone. This flashback is caused by turbulent burning along the vortex axis (TBVA¹) and is controlled by flame extinction inside the recirculation zone. The second type of flashback is caused by combustion induced vortex breakdown (CIVB²). This type of flashback is characterised by the aerodynamic influence of the combustion heat release that leads to propagation of the axial recirculation zone and the flame in upstream direction.To study the effects of geometrical scaling on the flow fields and the two types of flashback, the operation of two geometrically scaled burners are compared at equal Reynolds number. By this method it is possible to observe the flashback phenomena in similar swirl flow fields but with different turbulent scales affecting the combustion process. To check flow field similarity and to indentify the flashback type, the non-reacting and reacting flow fields have been examined by planar particle imaging velocimetry and simultaneous recording of the flame luminescence.It is shown that geometrical scaling of the burner shifts the equivalence ratio at which flashback occurs and that this shift is different for the two types of flashback. Consistency and inconsistency with known scaling and stability criterions is discussed. Analysing the fluid dynamics and turbulent combustion gives a first explanation of why CIVB and TBVA are affected differently by geometrical scaling at constant Reynolds number which is in good agreement with the experimental observations.     

Turbulent Combustion of Hydrogen–CO Mixtures

Laminar and turbulent burning velocities were measured in a closed-volume fan-stirred vessel for H₂–CO mixtures using two independent methods of flame definition. It has been shown that the unsteady flame development is an important factor and it needs to be taken into account for comparison of the burning rates obtained in different experiments. For the atmospheric pressure flames, the mixtures with faster laminar flame velocities burnt faster in turbulent flow despite the fact that the lean flames exhibit cellular structures. However, even a modest increase of the initial pressure promotes strongly cellularity and causes a significant acceleration of a lean laminar flame. The same lean flame burns faster in turbulent flow as well and this increase in the rate of combustion is greater that can be deduced from variation of the molecular heat diffusivity and laminar flame speed.     

Effects of Small- and Large-Scale Structures in a Piloted Jet Diffusion Flame

Laser Doppler Anemometry (LDA) and Planar Laser-Induced Fluorescence (PLIF) measurements have been performed in a turbulent nonpremixed jet flame. One of the features of this configuration is a central co-axial fuel jet surrounded by a turbulent annular air flow. The whole is placed within a low-speed coflowing air stream. This three-flow system with turbulent primary air differs from flow systems used for nonpremixed jet flames reported in the literature and is very useful for obtaining information on the mixing process between fuel and primary air. Next to the characterization of the velocity field, special attention has been paid to the conditional seeding of the central fuel jet and of the annular air flow. Together with visualizations of the OH radical, an important combustion intermediate which is formed during combustion, and the NO radical, which is seeded to the central jet flow, the resulting statistics reveal the properties of small- and large-scale structures in the flame.     

Flame Structure and Propagation in Turbulent Flame-Droplet Interaction: A Direct Numerical Simulation Analysis

Three-dimensional Direct Numerical Simulations (DNS) in canonical configuration have been employed to study the combustion of mono-disperse droplet-mist under turbulent flow conditions. A parametric study has been performed for a range of values of droplet equivalence ratio ?_d, droplet diameter a_d and root-mean-square value of turbulent velocity u^′. The fuel is supplied entirely in liquid phase such that the evaporation of the droplets gives rise to gaseous fuel which then facilitates flame propagation into the droplet-mist. The combustion process in gaseous phase takes place predominantly in fuel-lean mode even for ?_d>1. The probability of finding fuel-lean mixture increases with increasing initial droplet diameter because of slower evaporation of larger droplets. The chemical reaction is found to take place under both premixed and non-premixed modes of combustion: the premixed mode ocurring mainly under fuel-lean conditions and the non-premixed mode under stoichiometric or fuel-rich conditions. The prevalence of premixed combustion was seen to decrease with increasing droplet size. Furthermore, droplet-fuelled turbulent flames have been found to be thicker than the corresponding turbulent stoichiometric premixed flames and this thickening increases with increasing droplet diameter. The flame thickening in droplet cases has been explained in terms of normal strain rate induced by fluid motion and due to flame normal propagation arising from different components of displacement speed. The statistical behaviours of the effective normal strain rate and flame stretching have been analysed in detail and detailed physical explanations have been provided for the observed behaviour. It has been found that the droplet cases show higher probability of finding positive effective normal strain rate (i.e. combined contribution of fluid motion and flame propagation), and negative values of stretch rate than in the stoichiometric premixed flame under similar flow conditions, which are responsible for higher flame thickness and smaller flame area generation in droplet cases.     

Large-Scale Turbulent Vortical Structures inside a Sudden Expansion Cylinder Chamber

A large eddy simulation (LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model and the Lagrangian dynamic subgrid-scale model are employed and tested. The calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structures behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent behavior behind the flame holder is analyzed by visualizing the sectional views of vortical structure. This revised version was published online in July 2006 with corrections to the Cover Date.     

Joint RANS/LES Approach to Premixed Flame Modelling in the Context of the TFC Combustion Model

We present an original timesaving joint RANS/LES approach to simulate turbulent premixed combustion. It is intended mainly for industrial applications where LES may not be practical. It is based on successive RANS/LES numerical modelling, where turbulent characteristics determined from RANS simulations are used in LES equations for estimation of the subgrid chemical source and viscosity. This approach has been developed using our TFC premixed combustion model, which is based on a generalization of the Kolmogorov’s ideas. We assume existence of small-scale statistically equilibrium structures not only of turbulence but also of the reaction zones. At the same time, non-equilibrium large-scale structures of reaction sheets and turbulent eddies are described statistically by model combustion and turbulence equations in RANS simulations or follow directly without modelling in LES. Assumption of small-scale equilibrium gives an opportunity to express the mean combustion rate (controlled by small-scale coupling of turbulence and chemistry) in the RANS and LES sub-problems in terms of integral or subgrid parameters of turbulence and the chemical time, i.e. the definition of the reaction rate is similar to that of the mean dissipation rate in turbulence models where it is expressed in terms of integral or subgrid turbulent parameters. Our approach therefore renders compatible the combustion and turbulent parts of the RANS and LES sub-problems and yields reasonable agreement between the RANS and averaged LES results. Combining RANS simulations of averaged fields with LES method (and especially coupled and acoustic codes) for simulation of corresponding nonstationary process (and unsteady combustion regimes) is a promising strategy for industrial applications. In this work we present results of simulations carried out employing the joint RANS/LES approach for three examples: High velocity premixed combustion in a channel, combustion in the shear flow behind an obstacle and the impinging flame (a premixed flame attached to an obstacle).     

Efficient Thermo-Chemistry Tabulation for Non-Premixed Combustion at High-Pressure Conditions

Propellant injection and turbulent combustion in high-pressure engines is often dominated by real-gas effects. However, previous studies suggested that the departure of the fluid properties from an ideal gas behavior has only a limited effect on the laminar flame structure. This is due to the fact that chemical reactions take place in the flame zone where the temperature is sufficiently high and molecular interactions are negligible, i.e., the ideal gas assumption is valid. On the other hand, various experimental and numerical studies of injection processes at high-pressure conditions demonstrated that real-gas effects can have a strong impact on the turbulent flow. Mixing is influenced by the rapid change of fluid properties. In this work, we exploit the gap in the fidelity of the thermodynamics model needed to describe the laminar flame structure and that needed to describe the turbulent flow field. We then propose a new real-gas flamelet model with increased numerical performance. The computational cost of the new formulation is not significantly higher than that of an ideal gas simulation. The performance of the method is analyzed and the error that is introduced by our assumptions is assessed by comparison to more complete modeling. Finally, the method is used to simulate a turbulent jet flame emanating from a coaxial injector at supercritical pressure and cryogenic oxidizer temperature. The results are compared with experimental OH^? images giving evidence of the suitability of the present method.     

Numerical study of turbulent structures for lean premixed prevaporized combustion

In this paper we demonstrate the advantages of the LES-WALE model coupled with the PDF approach to resolve a set of aerothermochemistry equations for turbulent lean premixed prevaporized combustion. The main issue is the modeling of the closure in the turbulent combustion equations. So, combustion problems involve a strong coupling between dynamic and scalar parameters. The validation is based on comparisons of three parameters: mean longitudinal velocity, fluctuation of longitudinal velocity, and length of recirculation zones. In line with what was observed by an experimental reference study, the simulation succeeds to detect the flame zone and to model the flow morphology for different equivalence ratios and inlet mass flow rates.     

Effects of Turbulence on Self-sustained Combustion in Premixed Flame Kernels: A Direct Numerical Simulation (DNS) Study

The effects of mean flame radius and turbulence on self-sustained combustion of turbulent premixed spherical flames in decaying turbulence have been investigated using three-dimensional direct numerical simulations (DNS) with single step Arrhenius chemistry. Several flame kernels with different initial radius or initial turbulent field have been studied for identical conditions of thermo-chemistry. It has been found that for very small kernel radius the mean displacement speed may become negative leading ultimately to extinction of the flame kernel. A mean negative displacement speed is shown to signify a physical situation where heat transfer from the kernel overcomes the heat release due to combustion. This mechanism is further enhanced by turbulent transport and, based on simulations with different initial turbulent velocity fields, it has been found that self-sustained combustion is adversely affected by higher turbulent velocity fluctuation magnitude and integral length scale. A scaling analysis is performed to estimate the critical radius for self-sustained combustion in premixed flame kernels in a turbulent environment. The scaling analysis is found to be in good agreement with the results of the simulations.     

Schlieren measurement of turbulent structure in a diffusion flame

The regular and random mixing structures in a turbulent diffusion flame were investigated using the quantitative, dynamic crossed-beam schlieren method. Evidence was found close to the nozzle relating to the vortexlike structure of eddies surrounding the central fuel jet flow. The observations also make possible resolution of turbulent intensity, scales, convection, and spectra within the diffusion flame without the use of seeding or intrusion of measuring probes. It is found that length scales and other turbulence parameters in the diffusion flame progressively revert to values similar to those expected and observed in scalar passive mixing as the combustion reaction intensity reduces with axial distance from the nozzle system.     

Development of 3D Pocket Tracking Algorithm from Volumetric Measured Turbulent Flames

The flame pocket formation, including reactant pocket, product pocket, soot pocket, and fluid parcel, is a common phenomenon in turbulent combustion occurred as a response of the flame to flow straining and shearing. Understanding pocket behavior is vital to study the flames in such a regime. This work addresses the research need to experimentally measure and track multiple flame pockets in 3D. For this purpose, volumetric measurements were performed to measure the high-speed turbulent flame structure at 15 kHz based on emission tomography. With the 3D flame structures, a new tracking algorithm was developed to identify and track the multiple flame pockets simultaneously in 3D. The instantaneously tracked 3D flame pockets enabled the extraction of key properties of pocket dynamics, including the favorable formation location, 3D3C movement speed, and pocket expanding/shrinking speed. The developed methods were evidently able to resolve the detailed behavior of flame pockets in highly turbulent flames.

     

LES of Low to High Turbulent Combustion in an Elevated Pressure Environment

A subgrid scale flame surface density combustion model for the Large Eddy Simulation (LES) of premixed combustion is derived and validated. The model is based on fractal characteristics of the flame surface, assuming a self similar wrinkling of the flame between smallest and largest wrinkling length scales. Experimental and direct numerical simulation databases as well as theoretical models are used to derive a model for the fractal parameters, namely the cut-off lengths and the fractal dimension suitable in the LES context. The combustion model is designed with the intent to simulate low as well as high Reynolds number premixed turbulent flame propagation and with a focus on correct scaling with pressure. The combustion model is validated by simulations of turbulent Bunsen flames with methane and propane fuel at pressure levels between 0.1 MPa and 2 MPa and at turbulence levels of $0 < u^{\prime }/s_{L}^{0} < 11$ , conditions typical for spark ignition engines. The predicted turbulent flame speed is in a very good agreement with the experimental data and a smooth transition from resolved flame wrinkling to fully modelled, nearly subgrid-only wrinkling is realized. Evaluating the influence of mesh resolution shows a predicted mean flame surface and turbulent flame speed independent of mesh resolution for cases with 9–86 % resolved flame surface. Additional simulations of a highly turbulent jet flame at 0.1 MPa and 0.5 MPa and the comparison with experimental data in terms of flame shape, velocity field and turbulent fluctuations validates the model also at conditions typical for gas turbines.     

Subgrid Linear Eddy Mixing and Combustion Modelling of a Turbulent Nonpremixed Piloted Jet Flame

A linear eddy model for subgrid mixing and combustion has been coupled to a large eddy simulation of the turbulent nonpremixed piloted jet flame (Sandia Flame D). For the combustion reaction, simplified, single-step, irreversible, Arrhenius kinetics are used. The large scale and the subgrid structure of the flow are compared with experimental observations and, where appropriate, with a flamelet model of the flame. The main objective of this work is to demonstrate the feasibility of the LES-LEM approach for determining the structure of the subgrid scalar dissipation rate and the turbulence-chemistry interactions. The results for the large- and subgrid-scale structure of the flow show a reasonable agreement with the experimental observations.     

New Perspectives on Turbulent Combustion: Multi-Parameter High-Speed Planar Laser Diagnostics

Over the past three decades laser combustion diagnostics have guided an improved understanding of turbulent combustion processes. Until recently, this was based on statistically independent sampling using sampling rates much slower than typical integral time-scales of turbulent flames. Recent developments in laser and camera technology enabled an increase in sampling rates by more than three orders of magnitudes. Using these new instruments for particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) at high sampling rates (high-speed diagnostics) allowed the resolution of integral time-scales of turbulent flames. This statistically dependent sampling is increasingly used to temporally track transients in turbulent combustion, such as flame extinction, ignition, flashback and cycle-to-cycle variations in IC engines. The simultaneous application of flow and scalar field measurements makes insights into these transients possible that were not when using statistically independent sampling with low data acquisition rates. Conditioning on distinct flame features with high-speed diagnostics enables the inclusion of time as an additional dimension. This paper reviews the emerging field of multi-parameter, high-speed, planar laser diagnostics in combustion applications. The benefit of high data acquisition rates in turbulent combustion applications is discussed in detail as well as requirements and constraints imposed by the time-scales of the investigated phenomenon are addressed. Recent developments in laser and detector hardware are highlighted, as these are the limiting factors of the sampling rate. Finally, multi-parameter high-speed measurements in combustion are summarized, with a few examples discussed in more detail.     

Study of a Filtered Flamelet Formulation for Large Eddy Simulation of Premixed Turbulent Flames

Despite significant advances in the understanding and modelling of turbulent combustion, no general model has been proposed for simulating flames in industrial combustion devices. Recently, the increase in computational possibilities has raised the hope of directly solving the large turbulent scales using large eddy simulation (LES) and capturing the important time-dependant phenomena. However, the chemical reactions involved in combustion occur at very small scales and the modelling of turbulent combustion processes is still required within the LES framework. In the present paper, a recently presented model for the LES of turbulent premixed flames is presented, analysed and discussed. The flamelet hypothesis is used to derive a filtered source term for the filtered progress variable equation. The model ensures proper flame propagation. The effect of subgrid scale (SGS) turbulence on the flame is modelled through the flame-wrinkling factor. The present modelling of the source term is successfully tested against filtered direct numerical simulation (DNS) data of a V-shape flame. Further, a premixed turbulent flame, stabilised behind an expansion, is simulated. The predictions agree well with the available experimental data, showing the capabilities of the model for performing accurate simulations of unsteady premixed flames.