Large eddy simulation of bluff-body stabilized swirling non-premixed flames

Large eddy simulations (LES) using a subgrid mixing and combustion model are carried out to study two bluff-body stabilized swirling non-premixed flames (SM1 and SMA2). The similarities and differences between the two flames are highlighted and discussed. Flow features, such as, the recirculation zone (RZ) size and the flame structure are captured accurately in both cases. The SM1 flame shows a toroidal RZ just behind the bluff body and a vortex breakdown bubble (VBB) downstream. In addition, a highly rotational non-recirculating region in-between the RZ and VBB is observed as well. On the other hand, the SMA2 shows a single elongated recirculation zone downstream the bluff body. Flame necking is observed downstream the bluff body for the SM1 flame but not for the SMA2 flame. The time-averaged velocity and temperature comparison also shows reasonable agreement. The study shows that the sensitivity of the flame structure to inflow conditions can be captured in the present LES without requiring any model changes.     

Large eddy simulation of soot evolution in turbulent nonpremixed bluff body flames

Large Eddy Simulation (LES) is utilized to investigate soot evolution in a series of turbulent nonpremixed bluff body flames featuring different bluff body diameters. The modeling framework relies on recent development in the soot subfilter Probability Density Function (PDF) model that can correctly account for the distribution of soot with respect to mixture fraction, correcting errors in previous soot subfilter PDF models that significantly overpredict soot oxidation. With the previous soot subfilter PDF model, no soot was observed outside of the recirculation zone in past studies on similar bluff body flames. Results obtained with the current LES modeling approach compare favorably with the experimental measurements of the flow field and the soot volume fraction. Notably, the current LES modeling approach correctly predicts large soot volume fractions in the recirculation zone, a decrease in the soot volume fraction through the high-strain neck region, and then an increase again in the downstream jet-like region. Consistent with the experimental measurements, the larger bluff body diameter, with its larger recirculation zone with longer residence times, generates more soot in the recirculation zone and also more soot in the high-strain neck region. Analysis of the soot volume fraction source terms lead to mechanistic understanding of soot evolution in the entirety of the bluff body flames. Most of the soot generated in the recirculation zone is oxidized but some escapes unoxidized and is passively transported through the neck region. Virtually no new soot forms in the downstream jet-like region, and the increase in the soot volume fraction in the jet-like region is due to acetylene-based surface growth of the soot transported through the neck region. This mechanism could not be predicted with the previous soot subfilter PDF model, with the recent soot subfilter PDF model being critical in the understanding of this basic mechanism.     

LES of the Sydney swirl flame series: A study of vortex breakdown in isothermal and reacting flows

This work examines the flow and mixing in selected non-premixed cases of the Sydney Swirl Flame series by Large Eddy Simulation. A mixture fraction approach with a steady flamelet model, based on a detailed chemical mechanism, is applied to determine the chemical state in the flame. The isothermal case N29S054 is simulated to provide insights into the flow field, the resolution requirements for the simulation, and to allow for various measures of validation and verification. For the reactive case, the high-speed hydrogen/methane flame SMH1 is chosen for its similarity to the non-swirling Sydney flames and its good stability. In experiments carried out previously at Sydney University, vortex breakdown has been observed and in the isothermal case, this is clearly predicted by the LES. However, no vortex breakdown is observed in the simulations of flame SMH1, which necessitates further studies on this and similar flames to investigate this phenomenon. Studies of the low-velocity flames SM1 and SM2 show that reactive vortex breakdown can be predicted successfully. This difficulty in the prediction of vortex breakdown is another indication that the Sydney Swirl Flame series, especially at high velocities of the central jet, is an interesting and challenging test-case for the development of combustion LES.     

Subgrid combustion modeling of 3-D premixed flames in the thin-reaction-zone regime

Large eddy simulation (LES) is used to investigate three-dimensional (3D) lean premixed turbulent methane–air flames in the thin-reaction-zone regime. In this regime, the Kolmogorov scale is smaller than the preheat zone thickness, but larger than the reaction zone thickness. Past numerical studies of similar flames were primarily direct numerical simulation either in two-dimensions or using the artificially thickened flame approach in 3D. For an LES the effect of small (unresolved) scales on the scalar field must be, modeled accurately to capture the correct flame structure. A subgrid combustion model based on the linear-eddy-mixing (LEM) model is used within an LES framework (called LEM–LES hereafter) to capture the 3D flame-structure of the highly stretched premixed flames. A finite-rate, one-step methane–air chemistry with a non-unity Lewis number formulation is used in this study. The simulated flame structure resembles flames experimentally studied in the thin-reaction-zone regime. Even though the preheat zone is broadened by the penetration of small eddies, the chemical reaction zone remains thin and localized. This feature is captured properly in the current LEM–LES approach. The flame structure and other statistics such as the flame area evolution, curvature, and strain-rate statistics computed using the LEM–LES are also in good agreement with the past DNS studies.     

Large eddy simulation of piloting effects on turbulent swirling flames

Pilot flames, created by additional injectors of pure fuel, are often used in turbulent burners to enhance flame stabilization and reduce combustion instabilities. The exact mechanisms through which these additional rich zones modify the flame anchoring location and the combustion dynamics are often difficult to identify, especially when they include unsteady hydrodynamic motion. This study presents Large Eddy Simulations (LES) of the reacting flow within a large-scale gas turbine burner for two different cases of piloting, where either 2 or 6% of the total methane used in the burner is injected through additional pilot flame lines. For each case, LES shows how the pilot fuel injection affects both flame stabilization and flame stability. The 6% case leads to a stable flame and limited hydrodynamic perturbations in the initial flame zone. The 2% case is less stable, with a small-lift-off of the flame and a Precessing Vortex Core (PVC) in the cold stabilization zone. This PVC traps some of the lean cold gases issuing from the pilot passage stream, changes the flame stabilization point and induces instability.     

Self-oscillations in a jet flow and gaseous flame with strong swirl

Investigation results on unsteady flow dynamics in a gaseous jet flame with strong swirl, vortex breakdown, and precession of a vortex core obtained by panoramic optical methods are presented, as well as the results of theoretical analysis of the fastest growing modes of hydrodynamic instability. Characteristics of the most unstable self-oscillating mode in the initial region of the turbulent strongly swirling propane-air jet burning in the atmospheric air in the form of a lifted flame are determined. Analysis of data by principal component analysis and linear stability analysis revealed that evolution of the dominant self-oscillating mode corresponds to quasi-solid rotation with constant angular velocity of the spatial coherent structure consisting of a jet spiral vortex core and two spiral secondary vortices.     

Influence of pressure on near nozzle flow field and soot formation in laminar co-flow diffusion flames

Soot formation from combustion devices, which tend to operate at high pressure, is a health and environmental concern, thus investigating the effect of pressure on soot formation is important. While most fundamental studies have utilised the co-flow laminar diffusion flame configuration to study the effect of pressure on soot, there is a lack of investigations into the effect of pressure on the flow field of diffusion flames and the resultant influence on soot formation. A recent work has displayed that recirculation zones can form along the centreline of atmospheric pressure diffusion flames. This present work seeks to investigate whether these zones can form due to higher pressure as well, which has never been explored experimentally or numerically. The CoFlame code, which models co-flow laminar, sooting, diffusion flames, is validated for the prediction of recirculation zones using experimental flow field data for a set of atmospheric pressure flames. The code is subsequently utilised to model ethane-air diffusion flames from 2 to 33 atm. Above 10 atm, recirculation zones are predicted to form. The reason for the formation of the zones is determined to be due to increasing shear between the air and fuel steams, with the air stream having higher velocities in the vicinity of the fuel tube tip than the fuel stream. This increase in shear is shown to be the cause of the recirculation zones formed in previously investigated atmospheric flames as well. Finally, the recirculation zone is determined as a probable cause of the experimentally observed formation of a large mass of soot covering the entire fuel tube exit for an ethane diffusion flame at 36.5 atm. Previously, no adequate explanation for the formation of the large mass of soot existed.     

Effect of fuel mixture fraction and velocity perturbations on the flame transfer function of swirl stabilized flames

A novel methodology is developed to decompose the classic Flame Transfer Function (FTF) used in the thermo-acoustic stability analysis of lean premix combustors into contributions of different types. The approach is applied, in the context of Large Eddy Simulation (LES), to partially-premixed and fully-premixed flames, which are stabilized via a central recirculation zone as a result of the vortex breakdown phenomenon. The first type of decomposition is into contributions driven by fuel mixture fraction and dynamic velocity fluctuations. Each of these two contributions is further split into the components of turbulent flame speed and flame surface area. The flame surface area component, driven by the pure dynamic velocity fluctuation, which is shown to be a dominant contribution to the overall FTF, is also additionally decomposed over the coherent flow structures using proper orthogonal decomposition. Using a simplified model for the dynamic response of premixed flames, it is shown that the distribution of the FTF, as obtained from LES, is closely related to the characteristics of the velocity field frequency response to the inlet perturbation. Initially, the proposed method is tested and validated with a well characterized laboratory burner geometry. Subsequently, the method is applied to an industrial gas turbine burner.     

Soot studies of laminar diffusion flames with recirculation zones

This paper describes the unusual sooting structure of three flames established by the laminar recirculation zones of a centerbody burner. The vertically mounted burner consists of an annular air jet and a central fuel jet separated by a bluff-body. The three ethylene fueled flames are identified as: fully sooting, donut-shape, and ring-shape sooting flames. Different shapes of the soot structures are obtained by varying the N₂ dilution in the fuel and air jets while maintaining a constant air and fuel velocity of 1.2 m/s. All three flames have the unusual characteristic that the soot, entrained into the recirculation zone, follows discrete spiral trajectories that terminate at the center of the vortex. The questions are what cause: (1) the unusual sooting structures and (2) the spiral trajectories of the soot? Flame photographs, laser sheet visualizations, and calculations with a 2D CFD-based code (UNICORN) are used to answer these questions. The different sooting structures are related to the spiral transport of the soot, the spatial location of the stoichiometric flame surface with respect to the vortex center, and the burnout of the soot particles. Computations indicate that the spiral trajectories of the soot particles are due to thermophoresis.     

An experimental study of partially premixed flame sound

The occurrence of oscillating combustion and combustion instability has led to resurgence of interest in the causes, mechanisms, suppression, and control of combustion noise. Noise generated by enclosed flames is of greater practical interest but is more complicated than that by open flames, which itself is not clearly understood. Studies have shown that different modes of combustion, premixed and non-premixed, differ in their sound generation characteristics. However, there is lack of understanding of the region bridging these two combustion modes. This study investigates sound generation by partially premixed flames. Starting from a non-premixed flame, air was gradually added to achieve partial premixing while maintaining the fuel flow rate constant. Methane, ethylene, and ethane partially premixed flames were studied with hydrogen added for flame stabilization. The sound pressure generated by methane partially premixed flames scales with M⁵ compared to M³ for turbulent non-premixed methane flames. Also, the sound pressure generated by partially premixed flames of ethane and ethylene scales as M^4.5. With progressive partial premixing, spectra level increases at all frequencies with a greater increase in the high-frequency region compared to the low-frequency region; flames develop a peak and later a constant level plateau in the low frequency region. The partially premixed flames of methane, ethylene, and ethane generate a similar SPL as a function of equivalence ratio when the fuel volume flow rate is matched. However, when fuel mass flow rate is matched, the ethane and ethylene flames produce a similar SPL, which is lower than that produced by the methane flame.     

LES of oxy-fuel jet flames using the Eulerian Stochastic Fields method with differential diffusion

The Eulerian Stochastic Fields (ESF) Monte Carlo method to solve the transported PDF (TPDF) equation is extended to account for differential diffusion effects by incorporating species individual molecular diffusivities. The method has been applied in Large Eddy simulation (LES) to non-piloted oxy-fuel jet flames at different Reynolds numbers experimentally investigated by Sevault et al. [1]. Due to the high H₂ content in the fuel stream and CO₂ in the oxidizer these flames pose new challenges to combustion modeling as the flame structures are different compared to CH₄/air flames. The simulations show very good agreement with the experiments in terms of mixture fraction conditional mean values for temperature and mayor species on the fuel lean side and the reaction zone, deviations on the fuel rich side are discussed. The trend and location of localized extinction is reproduced well in the simulations, as well as differential diffusion effects in the near field. Additionally, it is shown that a neglect of differential diffusion in the combustion model leads to a lifted flame.     

Effect of swirl on combustion dynamics in a lean-premixed swirl-stabilized combustor

The effect of inlet swirl on the flow development and combustion dynamics in a lean-premixed swirl-stabilized combustor has been numerically investigated using a large-eddy-simulation (LES) technique along with a level-set flamelet library approach. Results indicate that when the inlet swirl number exceeds a critical value, a vortex-breakdown-induced central toroidal recirculation zone is established in the downstream region. As the swirl number increases further, the recirculation zone moves upstream and merges with the wake recirculation zone behind the centerbody. Excessive swirl may cause the central recirculating flow to penetrate into the inlet annulus and lead to the occurrence of flame flashback. A higher swirl number tends to increase the turbulence intensity, and consequently the flame speed. As a result, the flame surface area is reduced. The net heat release, however, remains almost unchanged because of the enhanced flame speed. Transverse acoustic oscillations often prevail under the effects of strong swirling flows, whereas longitudinal modes dominate the wave motions in cases with weak swirl. The ensuing effect on the flow/flame interactions in the chamber is substantial.     

Impact of co- and counter-swirl on flow recirculation and liftoff of non-premixed oxy-flames above coaxial injectors

Controlling the flame shape and its liftoff height is one of the main issues for oxy-flames to limit heat transfer to the solid components of the injector. An extensive experimental study is carried out to analyze the effects of co- and counter-swirl on the flow and flame patterns of non-premixed oxy-flames stabilized above a coaxial injector when both the inner fuel and the annular oxidizer streams are swirled. A swirl level greater than 0.6 in the annular oxidizer stream is shown to yield compact oxy-flames with a strong central recirculation zone that are attached to the rim of central fuel tube in absence of inner swirl. It is shown that counter-swirl in the fuel tube weakens this recirculation zone leading to more elongated flames, while co-swirl enhances it with more compact flames. These results obtained for high annular swirl levels contrast with previous observations made on gas turbine injectors operated at lower annular swirl levels in which central recirculation of the flow is mainly achieved with counter-rotating swirlers. Imparting a high inner swirl to the central fuel stream leads to lifted flames due to the partial blockage of the flow at the injector outlet by the central recirculation zone that causes high strain rates in the wake of the injector rim. This partial flow blockage is more influenced by the level of the inner swirl than its rotation direction. A global swirl number is then introduced to analyze the structure of the flow far from the burner outlet where swirl dissipation takes place when the jets mix. A model is derived for the global swirl number which well reproduces the evolution of the mass flow rate of recirculating gases measured in non-reacting conditions and the flame liftoff height when the inner and outer swirl levels and the momentum flux ratio between the two streams are varied.     

预混气体燃烧火焰闪烁现象分析

在低速射流的预混火焰和扩散火焰中都存在火焰闪烁现象。对扩散火焰,其机理已比较明确,是由于浮力诱导引起的一种水力学不稳定性。而对预混火焰闪烁现象则存在水力学不稳定性和热驱动不稳定性两种观点。本文根据水力学不不稳定性观点,把预混火焰的闪烁现象看成是包围火焰锋面的已燃混气层中内、外区间在垂直方向上的相对脉动,应用Ｋｅｌｖｉｎ－Ｈｅｌｍｈｏｌｔｚ不稳定性机理进行了分析,获得了火焰闪烁频率与重力和压力的关系式,并与已有的结果作了对比。     

Characteristics and structure of inverse flames of natural gas

Characteristics and structure of nominally non-premixed flames of natural gas are investigated using a burner that employs simultaneously two distinct features: fuel and oxidiser direct injection, and inverse fuel and oxidiser delivery. At low exit velocities, the result is an inverse diffusion flame that has been noted in the past for its low NO_x emissions, soot luminosity, and narrow stability limits. The present study aimed at extending the burner operating range, and it demonstrated that the inverse flame exhibits a varying degree of partial premixing dependent on the discharge nozzle conditions and the ratio of inner air jet and outer fuel jet velocities. These two variables affect the flame length, temperature distributions, and stability limits. Temperature measurements and Schlieren visualisation show areas of enhanced turbulent mixing in the shear region and the presence of a well-mixed reaction zone on the flame centreline. This reaction zone is enveloped by an outer diffusion flame, yielding a unique double-flame structure. As the fuel–air equivalence ratio is decreasing with an increase in the inner jet velocity, the well-mixed reaction zone extends considerably. These findings suggest a method for establishing a flame of uniform high temperature by optimising the coaxial nozzle geometry and flow conditions. The normalised flame length is decreasing exponentially with the air/fuel velocity ratio. Measurements demonstrate that the inverse flame stability limits change qualitatively with varying degree of partial premixing. At the low premixing level, the flame blow-out is a function of the inner and outer jet velocities and the nozzle conditions. The flame blow-out at high degree of partial premixing occurs abruptly at a single value of the inner air jet velocity, regardless of the fuel jet velocity and almost independent of the discharge nozzle conditions.     

重力对扩散射流火焰动态特性的影响

本文探讨重力对扩散射流火焰动态特性的影响规律。结果表明,火焰闪烁现象是一种浮力诱导不稳定性,在浮力消失或反向重力场中,不存在这种不稳定性现象,闪烁频率与燃料射流速度无直接关系,但涡的大小随燃料射流速度的增大而增大。存在触发火焰闪烁的临界高度,闪烁频率与重力成平方根关系式。反向重力情况下,也存在浮力稳定型平面火焰,它反映了浮力与火焰的耦合作用。     

Stability characteristics and flowfields of turbulent non-premixed swirling flames

A simple, yet representative, burner geometry is used for the investigation of highly swirling turbulent unconfined, non-premixed, flames of natural gas. The burner configuration comprises a ceramic faced bluff-body with a central fuel jet. The bluff-body is surrounded by an annulus that delivers a swirling primary flow of air. The entire burner assembly is housed in a wind tunnel providing a secondary co-flowing stream of air. This hybrid bluff-body/swirl burner configuration stabilizes complex turbulent flames not unlike those found in practical combustors, yet is amenable to modelling because of its well-defined boundary conditions. Full stability characteristics including blow-off limits and comprehensive maps of flame shapes are presented for swirling flames of three different fuel mixtures: compressed natural gas (CNG), CNG–air (1:2 by volume) and CNG–H₂ (1:1 by volume).

It is found that with increased fuel flow, flame blow-off mode may change with swirl number, S_g. At low swirl, the flame remains stable at the base but blows off in the neck region further downstream. At higher swirl numbers, the flames peel off completely from the burner's base. Swirling CNG–air flames are distinct in that they only undergo base blow-off. In the low range of swirl number, increasing S_g causes limited improvement in the blow-off limits of the flames investigated and (for a few cases) can even lead to some deterioration over a small intermediate range of S_g. It is only above a certain threshold of swirl that significant improvements in blow-off limits appear. Six flames are selected for further detailed flowfield and composition measurements and these differ in the combination of swirl number, primary axial velocity through the annulus, U_s, and bulk fuel jet velocity, U_j. Only velocity field measurements are presented in this paper. A number of flow features are resolved in these flames, which resemble those already associated with non-reacting swirling flows of equivalent swirl obtained with the present burner configuration. Additionally, asymmetric flowfields inherent to some flames are revealed where the fluidic centreline of the flow (defined in the two-dimensional (U–W velocity pair) velocity field by the ?ω? = 0 tangential velocity contour), meanders strongly on either side of the geometric centreline downstream by about one bluff-body diameter. Flow structures revealed by the velocity data are correlated to flame shapes to yield a better understanding of how the velocity field influences the flames physical characteristics.     

LES-CMC simulations of a turbulent bluff-body flame

The large Eddy simulations (LES)-conditional moment closure (CMC) method with detailed chemistry is applied to a bluff-body stabilized flame. Computations of the velocity and mixture fraction fields show good agreement with the experiments. Temperature and major species are well-predicted throughout the flame with the exception of the flow regions in the outer shear layer close to the nozzle where the pure mixing between hot recirculating products and fresh oxidizer cannot be captured. LES-CMC generally improves on results obtained with RANS-CMC and on LES that uses one representative flamelet to model the dependence of reactive species on mixture fraction. Simulated CO mass fractions are generally in good agreement with the experimental data although a 10% overprediction can be found at downstream positions. NO predictions show a distinct improvement over the flamelet approach, however, simulations overpredict NO mass fractions at all downstream locations due to an overprediction of temperature close to the nozzle. The potential of LES-CMC to predict unsteady finite rate effects is demonstrated by the prediction of endothermic—or “flame cooling”—regions close to the neck of the recirculation zone that favours ethylene production via the methane fuel decomposition channel.     

PDF calculations of piloted premixed jet flames

A series of piloted premixed jet flames with strong finite-rate chemistry effects is studied using the joint velocity-turbulence frequency-composition PDF method. The numerical accuracy of the calculations is demonstrated, and the calculations are compared to experimental data. It is found that all calculations show good agreement with the measurements of mean and rms mixture fraction fields, while the reaction progress is overpredicted to varying degrees depending on the jet velocity. In the calculations of the flame with the lowest jet velocity, the species and temperature show reasonable agreement with the measurements, with the exception of a small region near the centerline where products and temperature are overpredicted and fuel and oxidizer are underpredicted. In the calculations of the flame with the highest jet velocity, however, the overprediction of products and temperature and underprediction of fuel and oxidizer is far more severe. An extensive set of sensitivity studies on inlet boundary conditions, turbulence model constants, mixing models and constants, radiation treatment, and chemical mechanisms is conducted to show that any parameter variation offers little improvement from the base case. To shed light on these discrepancies, diagnostic calculations are performed in which the chemical reactions are artificially slowed. These diagnostic calculations serve to validate the experimental data and to quantify the amount by which the base case calculations overpredict reaction progress. Improved calculations of this flame are achieved only through artificially slowing down the chemical reaction by a factor of about 10. The mixing model behavior in this combustion regime is identified as a likely cause for the observed discrepancy in reaction progress.