Identification and analysis of instability in non-premixed swirling flames using LES

Large eddy simulations (LES) of turbulent non-premixed swirling flames based on the Sydney swirl burner experiments under different flame characteristics are used to uncover the underlying instability modes responsible for the centre jet precession and large scale recirculation zone. The selected flame series known as SMH flames have a fuel mixture of methane-hydrogen (50:50 by volume). The LES solves the governing equations on a structured Cartesian grid using a finite volume method, with turbulence and combustion modelling based on the localised dynamic Smagorinsky model and the steady laminar flamelet model respectively. The LES results are validated against experimental measurements and overall the LES yields good qualitative and quantitative agreement with the experimental observations. Analysis showed that the LES predicted two types of instability modes near fuel jet region and bluff body stabilised recirculation zone region. The mode I instability defined as cyclic precession of a centre jet is identified using the time periodicity of the centre jet in flames SMH1 and SMH2 and the mode II instability defined as cyclic expansion and collapse of the recirculation zone is identified using the time periodicity of the recirculation zone in flame SMH3. Finally frequency spectra obtained from the LES are found to be in good agreement with the experimentally observed precession frequencies.     

Analysis of kinetic mechanism performance in conditional moment closure modelling of turbulent,non-premixed methane flames

This paper presents results obtained from the application of a first-order conditional moment closure approach to the modelling of two methane flames of differing geometries. Predictions are based upon a second-moment turbulence and scalar-flux closure, and supplemented with full and reduced chemical kinetic mechanisms, ranging from a simple 12-step to a complex 1207-step mechanism. Alongside analysis of the full kinetic schemes' performance, is an appraisal of the behaviour of their derivatives obtained using mechanism-reduction techniques. The study was undertaken to analyse the practicality of incorporating kinetic models of varying complexity into calculations of turbulent non-premixed flames, and to make comparison of their performance. Despite extensive studies of the predictive ability of such schemes under laminar flame conditions, systematic evaluations have not been performed for turbulent reacting flows. This paper reflects upon the impact that selection of chemical kinetics has upon subsequent calculations and concludes that, although application of reduced schemes is more than adequate to reproduce experimental data, selection of the parent mechanism is of paramount importance to the prediction of minor species. Although widely used schemes are well documented and validated, their performances vary considerably. Thus, careful consideration must be made to their application and origins during the evaluation of combustion models.     

Edge flames stabilized in a non-premixed microcombustor

The dynamics of an edge flame confined in a non-premixed microcombustor model is studied numerically within the context of a diffusive-thermal model. Fuel and oxidizer, separated upstream by a thin plate, flow through a channel with a prescribed velocity. At the tip of the plate, the fuel and oxidizer mix and, when ignited, an edge flame is sustained at some distance from the plate. The objective in this work is to consider the effects of confinement, differential diffusion, and heat loss on the dynamics of an edge flame in a narrow channel. We consider a wide range of channel widths and allow for changing Lewis numbers, and both adiabatic conditions and heat losses along the channel walls. The results illustrate how the flame shape and standoff distance are affected by the channel width, by mixture composition through variations in Lewis numbers and by heat losses. Conditions for flame stabilization, flame oscillations and flame extinction or blowoff are predicted.     

Numerical investigation of the effects of polydisperse water sprays on extinction conditions of counterflow methane non-premixed flames

Water, sprayed in the form of tiny droplets, has emerged as a potential fire suppressant after the halon compounds such as trifluorobromomethane (CF₃Br, Halon 1301) were banned by the Montreal protocol. The size distribution of the water droplet plays a crucial role in the effectiveness of the water spray in fire suppression. A numerical investigation of the influence of size distribution of a polydisperse water spray on extinction of counterflow diffusion flames is presented in this paper. This study uses laminar finite rate model with reduced CHEMKIN chemistry for numerical simulations. The discrete phase, namely the water spray, is simulated using Lagrangian Discrete Phase Modelling approach. In this work, the polydispersity of water spray is taken into account in the numerical simulation by a suitable Rosin–Rammler distribution. Results obtained from numerical simulation are validated with the experimental results reported in the literature. This study demonstrates that the representation of the polydisperse spray by a monodisperse spray (with droplet diameter same as the SMD of the polydisperse spray) in numerical simulations is not always justified and it leads to deviation from the experimental results. The effects of number mean diameter and spread parameter on the efficacy of flame suppression are investigated for polydisperse sprays. A comprehensive comparison is done between the effectiveness of monodisperse and polydisperse water sprays. An optimum droplet diameter is obtained for monodisperse sprays for which the effectiveness of the spray is maximum. The effects of evaporation Damköhler number and Stokes number of water droplets on flame suppression have also been explained.     

Ignition of non-premixed cyclohexane and mono-alkylated cyclohexane flames

Ignition temperatures of non-premixed cyclohexane, methylcyclohexane, ethylcyclohexane, n-propylcyclohexane, and n-butylcyclohexane flames were measured in the counterflow configuration at atmospheric pressure, a free-stream fuel/N₂ mixture temperature of 373 K, a local strain rate of 120 s^?1, and fuel mole fractions ranging from 1% to 10%. Using the recently developed JetSurf 2.0 kinetic model, satisfactory predictions were found for cyclohexane, methyl-, ethyl-, and n-propyl-cyclohexane flames, but the n-butylcyclohexane data were overpredicted by 20 K. The results showed that cyclohexane flames exhibit the highest ignition propensity among all mono-alkylated cyclohexanes and n-hexane due to its higher reactivity and larger diffusivity. The size of mono-alkyl group chain was determined to have no measurable effect on ignition, which is a result of competition between fuel reactivity and diffusivity. Detailed sensitivity analyses showed that flame ignition is sensitive primarily to fuel diffusion and also to H₂/CO and C₁–C₃ hydrocarbon kinetics.     

Response of non-premixed flames to bulk flow perturbations

This paper describes the dynamics of non-premixed flames responding to bulk velocity fluctuations, and compares the dynamics of the flame sheet position and spatially integrated heat release to that of a premixed flame. The space–time dynamics of the non-premixed flame sheet in the fast chemistry limit is described by the stoichiometric mixture fraction surface, extracted from the solution of the

-equation. This procedure has some analogies to premixed flames, where the premixed flame sheet location is extracted from the G = 0 surface of the solution of the G-equation. A key difference between the premixed and non-premixed flame dynamics, however, is the fact that the non-premixed flame sheet dynamics are a function of the disturbance field everywhere, and not just at the reaction sheet, as in the premixed flame problem. A second key difference is that the non-premixed flame does not propagate and so flame wrinkles are convected downstream at the axial flow velocity, while wrinkles in premixed flames convect downstream at a vector sum of the flame speed and axial velocity. With the exception of the flame wrinkle propagation speed, however, we show that that the solutions for the space–time dynamics of the premixed and non-premixed reaction sheets in high velocity axial flows are quite similar. In contrast, there are important differences in their spatially integrated unsteady heat release dynamics. Premixed flame heat release fluctuations are dominated by area fluctuations, while non-premixed flames are dominated by mass burning rate fluctuations. At low Strouhal numbers, the resultant sensitivity of both flames to flow disturbances is the same, but the non-premixed flame response rolls off slower with frequency. Hence, this analysis suggests that non-premixed flames are more sensitive to flow perturbations than premixed flames at O(1) Strouhal numbers.     

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.     

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.     

Ignition of non-premixed counterflow flames of octane and decane isomers

Ignition temperatures of non-premixed flames of octane and decane isomers were determined in the counterflow configuration at atmospheric pressure, a free-stream fuel/N₂ mixture temperature of 401 K, a local strain rate of 130 s^?1, and fuel mole fractions ranging from 1% to 6%. The experiments were modeled using detailed chemical kinetic mechanisms for all isomers that were combined with established H₂, CO, and n-alkane models, and close agreements were found for all flames considered. The results confirmed that increasing the degree of branching lowers the ignition propensity. On the other hand, increasing the straight chain length by two carbons was found to have no measurable effect on flame ignition for symmetric branched fuel structures. Detailed sensitivity analyses showed that flame ignition is sensitive primarily to the H₂/CO and C₁–C₃ hydrocarbon kinetics for low degrees of branching, and to fuel-related reactions for the more branched molecules.     

Attachment structure of a non-premixed laminar methane flame

The stoichiometry and the flame structure of the leading edge, an anchor point, of a non-premixed methane flame were investigated. Local equivalence ratio at an anchor point was measured using local chemiluminescence spectra with a high spatial resolution of 17 × 450 μm. Spatially and spectrally resolved chemiluminescence measurements were carried out along the centerline and radius of the non-premixed laminar flame. The chemiluminescence spectra measured at the flame tip contained very strong luminous spectra, while these continuous background spectra disappeared at the blue flame tip region. The chemiluminescence spectra below the blue flame region were very similar to those measured in laminar premixed methane/air flames. Based on these results, the local equivalence ratio near the anchor point was calculated. Therefore, we measure the anchor point location, its shape, and stoichiometry using the flame spectra. At the anchor point, there was an island of lower equivalence ratio of 0.65, which can be estimated as the lower flammable limit of premixed laminar flame. The size of the anchor point was of horizontal elliptical shape less than 0.6 and 0.4 mm in vertical length, which located at 1.2 mm above the burner rim and inside of the rim.     

A comparative study of turbulence models for non-premixed swirl-stabilized flames

ABSTRACT

The accuracy of turbulent swirl-stabilized flame simulation strongly depends on the choice of turbulence model. In this study, four 3D unsteady turbulence closures, including large eddy simulation, scale-adaptive simulation, and two detached eddy simulation variants, along with four RANS models, including RNG k??, SST k?ω, transition SST, and RSM, are examined for moderate- and high-swirl case studies. It is observed that the scale-adaptive simulation provides the most accurate results for almost all variables and both swirl conditions in the reactive flow. Only the 3D unsteady models predict the vortex breakdown bubble and flame attachment state correctly. However, based on our error analysis, the flow and composition fields predicted by the RANS models are in acceptable agreement with the experimental fields, especially the ones of transition SST when higher swirl number cases or minor species concentration are of interest. Moreover, it is concluded that the viscosity ratio criterion is a better measure of the local LES quality than the turbulent kinetic energy ratio, and the accuracy of a hybrid simulation may be much more dependent on the ability of the model to operate close to the RANS mode where the grid resolution is not sufficient for a resolving simulation than the fraction of the resolved kinetic energy. Finally, the propriety of the base (RANS) model of a DES for the application of interest is important, such that DES with realizable k?? outperforms the commonly used DES with SST k?ω model.     

Assessment of subgrid-scale stress statistics in non-premixed turbulent wall-jet flames

We investigate the heat-release effects on the characteristics of the subgrid-scale (SGS) stress tensor and SGS dissipation of kinetic energy and enstrophy. Direct numerical simulation data of a non-premixed reacting turbulent wall-jet flow with and without substantial heat release is employed for the analysis. This study comprises, among others, an analysis of the eigenvalues of the resolved strain rate and SGS stress tensors, to identify the heat-release effects on their topology. An assessment of the alignment between the eigenvectors corresponding to the largest eigenvalues of these two tensors is also given to provide further information for modelling of the SGS stress tensor. To find out the heat-release effects on the dynamics of the turbulent kinetic energy and enstrophy dissipation, probability density functions (PDFs) and mean values are analysed. The mean SGS shear stress and turbulent kinetic energy both slightly increase in the buffer layer and substantially decrease further away from the wall, due to the heat-release effects. Contrary to the kinetic energy, heat release decreases the mean SGS dissipation of enstrophy in the near-wall region. Moreover, differences in the shapes of the PDFs between the isothermal and exothermic cases indicate changes in the intermittency level of both SGS dissipations. Heat release also increases the SGS stress anisotropy in the near-wall region. Although, the structure of the mean resolved strain-rate tensor only marginally differs between the isothermal and exothermic cases in the near-wall region, substantial differences are observed in the jet area, where compressibility effects are important and heat-release effects are found to promote compression states. The differences in the relative alignment between the SGS stress and resolved strain-rate tensors in the isothermal and exothermic cases are discussed in connection with the differences in the SGS dissipation of kinetic energy.     

Nano organic carbon and soot in turbulent non-premixed ethylene flames

Spectral optical techniques are combined to characterise the distribution of large-molecule soot precursors, nanoparticles of organic carbon, and soot in two turbulent non-premixed ethylene flames with differing residence times. Laser-induced fluorescence, laser-induced incandescence and light scattering are used to define distributions across the particle size distribution. From the scattering and laser-induced emission measurements it appears that two classes of particles are formed. The first ones are preferentially formed in the fuel-rich region of the flame closer to the nozzle, have sizes of the order of few nanometers but are not fully solid particles, because the constituent molecules still maintain their individual identity exhibiting strong broadband fluorescence in the UV. The second class of particles constituted by solid particles, with sizes of the order of tens of nanometers are able to absorb a sufficient number of photons to be heated to incandescent temperatures. These larger particles are formed at larger residence times in the flame since they are the result of slow growth processes such as coagulation or carbonization. The flames are also modeled in order to produce mixture fraction maps. A new discovery is that nanoparticles of organic carbon concentration, unlike soot, does correlate well with mixture fraction, independent of position in the flame. This is likely to be a significant benefit to future modelling of soot inception processes in turbulent non-premixed flames.     

Cell formation in non-premixed,axisymmetric jet flames near extinction

Systematic experiments with CO₂ diluted H₂–O₂ circular jet diffusion flames have been undertaken to study the formation of cellular flames, which occur for relatively low reactant Lewis numbers and near the extinction limit. The jet Reynolds number for all experiments was about 500, based on the centreline velocity, jet diameter and ambient fuel properties. The Lewis numbers, based on the initial mixture strength φ and ambient conditions of the investigated near-extinction mixtures, vary in the range 1.1–1.3 for oxygen and 0.25–0.29 for hydrogen (φ is defined here as the fuel-to-oxygen molar ratio normalized by the stoichiometric value). Various conditions near the extinction limit were investigated by fixing the fuel composition (H₂–CO₂ mixture), and systematically lowering the oxygen concentration in the co-flowing oxidizer stream past the point where cellular structures formed, until extinction occurred. The observed different instability states were correlated with the initial mixture strength and the proximity to the extinction limit.

The parameter space for cellularity was found to increase with decreasing initial mixture strength. For a given initial mixture strength, several cellular states were found to co-exist near the extinction limit, and the preferred number of cells (the azimuthal wave number) was observed to decrease with decreasing oxygen concentration (Damköhler number). These trends are consistent with previous theoretical work and our own stability analysis that will be reported elsewhere.     

Conditional analysis of temperature and strain rate effects on dissipation structure in turbulent non-premixed jet flames

This work presents results from simultaneous high-resolution temperature and velocity measurements in a series of turbulent non-premixed jet flames. The filtered Rayleigh scattering (FRS)-based temperature measurements demonstrate sufficient signal-to-noise (SNR) and spatial resolution to estimate the smallest scalar length scales and accurately determine dissipation rate fields. A comprehensive set of conditional statistics are used to characterize the small-scale structure, including the dependence of dissipation layer widths on Reynolds number, temperature, and dissipation magnitude. In general, the dissipation layer thickness decrease with increasing Reynolds number and increase with increasing temperature. However, dissipation layer widths show two distinct behaviors with respect to dissipation magnitude. For small dissipation values, increases in magnitude results in broadening of the dissipation layer, while for larger magnitude values of dissipation, the layer widths are thinned, highlighting the complexity of small-scale turbulent mixing. Additionally, measured ratios of the dissipation layer width to the Batchelor length scale are consistent across all Reynolds numbers and agree with previous studies in non-reacting flows. The unique aspect about the current set of measurements is the ability to examine the interaction of dissipation structure with turbulent flow parameters for the first time in turbulent non-premixed flames. Particularly, the strain rate/dissipation relationship is examined and compared to previous studies in non-reacting flows. It is found that the dissipation layers tend to align normal to the principal compressive strain axis and this tendency increases with increasing Reynolds number. For the lowest Reynolds number case, no dependence of the dissipation layer width nor dissipation rate magnitude on strain rate is found. However, for higher Reynolds numbers, a strong dependence of the dissipation layer width and dissipation rate magnitude on the principal compressive strain rate is observed. These results indicate the direct role of the compressive strain rate field on small-scale mixing structure in reacting flows.