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.     

Critical conditions at liftoff limit of a laminar n-butane-air jet diffusion flame

The stability mechanism of laminar coflow jet diffusion flames in normal gravity has been studied computationally and experimentally. N-butane, the heaviest alkane in a gaseous state at ambient temperature and pressure, is used as the fuel since the reaction mechanism is similar to that of higher (liquid) hydrocarbons. The critical mean n-butane jet and coflowing air velocities at flame stability limits are measured using a small fuel tube burner (0.8 mm inner diameter). The time-dependent, axisymmetric numerical code with a detailed reaction mechanism (58 species and 540 reactions), molecular diffusive transport, and a radiation model, reveals a flame structure. A fuel-lean peak reactivity spot (i.e., reaction kernel), possessing the hybrid nature of diffusion-premixed flame structure at a constant temperature of ≈1560 K, is formed at the flame base and controls the flame stability. In a near-quiescent environment, the flame base resides below the fuel tube exit plane and thereby premixing is limited. As the coflowing air velocity is increased incrementally under a fixed fuel jet velocity, the flame base moves slightly above (≈1 mm) the burner exit and vigorous premixed combustion becomes prevailing. The local heat-release rate at the reaction kernel nearly doubles due to the increased convective oxygen flux (i.e., a blowing effect). The local Damköhler number, newly defined as a ratio of the square root of the local heat-release rate and the local velocity, decreases gradually first and drops abruptly at a critical threshold value and the flame base lifts off from the burner rim. The calculated coflow air velocity at liftoff is ≈0.38 m/s at the fuel jet velocity of 2 m/s, which is consistent with an extrapolated measured value of 0.41 m/s. This work has determined the critical Damköhler number at the stability limit quantitatively, for the first time, for laminar jet diffusion flames.     

Bifurcation characteristics and flame dynamics of a ducted non-premixed flame with finite rate chemistry

The influence of system parameters such as the flame location, Peclet number and Damköhler number on the bifurcation characteristics and flame dynamics of a ducted non-premixed flame with finite rate chemistry is presented in this paper. In the bifurcation plot with flame location as the bifurcation parameter, subcritical Hopf bifurcation is found for lower values of flame location and supercritical Hopf bifurcation for higher values of flame location, for all the Damköhler numbers used in this study. The flame shapes are captured at eight different phases of a cycle of time series data of acoustic velocity at both the fold and Hopf points for bifurcation with flame location as the parameter. We find that the range of flame height variations at the Hopf point is more than the range of flame height variations obtained at the fold point. We also find that the flame oscillates in the same phase as pressure fluctuation but in a phase different from both velocity and heat release rate fluctuations in the region of hysteresis for bifurcation with flame location. The non-dimensional hysteresis width is plotted as a function of Damköhler number for variation of flame location in the subcritical region. An inverse power law relation is found between the non-dimensional hysteresis width and the Damköhler number. The bifurcation plot with Peclet number as parameter shows a subcritical Hopf bifurcation.     

双组分燃料的层流火焰速度模型

以往关于层流火焰速度的理论分析均只考虑单组分燃料,本文对双组分燃料的平面火焰进行了大活化能渐近理论分析。在理论分析中,将火焰结构分为预热区、化学反应区和平衡区,并在大活化能假设下对各个区域分别求解了关于温度与燃料质量分数的微分方程。根据每两个区域分界面上满足的结合条件,本文推导出了双组分燃料的层流火焰速度模型。该模型表明双组分燃料层流火焰速度的平方为各个单组分燃料层流火焰速度平方的加权平均。     

贫燃料预混燃烧的回火特性研究

回火问题是贫燃料预混燃烧面临的主要问题之一。本文采用计算和实验相结合的方法研究甲烷与富氢合成气贫预混燃烧的回火现象,得到不同燃料、不同稳定方式之间的回火特性。研究结果表明,回火极限可以关联为丕雷数模型,环形稳定器的回火稳定性最好,其次为杆稳定器,旋流稳定器的稳定性最差;环形稳定的甲烷预混火焰的回火过程为边缘稳定,适当加入边缘空气同轴射流后变为中心回火,且同轴射流速度存在最佳范围可以提高回火稳定性。     

Measurement of unstable burning velocities of iso-octane–air mixtures at high pressure and the derivation of laminar burning velocities

A new technique is reported for measuring burning velocities at high pressures in the final stages of two inwardly propagating flame kernels in an explosion bomb. The flames were initiated at diametrically opposite spark electrodes, close to the wall, in quiescent mixtures. Measurements of pressure and flame kernel propagation speeds by high-speed photography showed the burning velocities to be elevated above the corresponding laminar burning velocities as a result of the developing flame instabilities. The enhancement increased with increase in pressure and decreased with increase in Markstein number. When the Markstein number was negative, instabilities could be appreciable, as could the enhancement. For the iso-octane–air mixtures investigated, where the mixtures had well-characterised Markstein numbers or critical Peclet numbers at the relevant pressures and temperatures, it was possible to explain the enhancement quantitatively by the spherical explosion flame instability theory of Bechtold and Matalon, provided the critical Peclet number was that observed experimentally, and allowance was made for the changing pressure. With this theoretical procedure, it was possible to derive values of laminar burning velocity from the measured values of burning velocity over a wide range of equivalence ratios, pressures, and temperatures. The values became less reliable at the higher temperatures and pressures as the data on Markstein and critical Peclet numbers became less certain. It was found that with iso-octane as the fuel the laminar burning velocity decreased during isentropic compression.     

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.     

Triple flame: Inherent asymmetries and pentasectional character

A two-dimensional triple-flame numerical model of a laminar combustion process reflects flame asymmetric structural features that other analytical models do not generate. It reveals the pentasectional character of the triple flame, composed of the central pure diffusion-flame branch and the fuel-rich and fuel-lean branches, each of which is divided into two sections: a near-stoichiometric section and a previously unreported near-flammability-limits section with combined diffusion and premixed character. Results include propagation velocity, fuel and oxidiser mass fractions, temperature and reaction rates. Realistic stoichiometric ratios and reaction orders match experimental planar flame characteristics. Constant density, a one-step reaction, and a mixture fraction gradient at the inlet as the simulation parameter are imposed. The upstream equivalence ratio or the upstream reactant mass fractions are linear or hyperbolic functions of the transverse coordinate. The use here of experimental kinetics data differs from previous analytical works and results in flame asymmetry and different flammability limits. Upstream mixture composition gradient affects propagation velocity, flame curvature, diffusion flame reaction rate, and flammability limits. Flammability limits extend beyond those of a planar flame due to transverse heat and mass diffusion causing the pentasectional character.     

微尺度扩散火焰特性的数值解析

本文以均匀空气流中圆管形成的甲烷射流扩散火焰为对象,用数值解析的方法研究了微尺度扩散火焰的火焰结构和燃烧特性。燃烧反应采用甲烷/空气一步总包反应,喷管壁面采用绝热条件。在Re一定情况下,改变喷口尺寸和喷口流速考察了微扩散火焰的结构和火焰熄灭的尺度效应。数值结果表明,随着喷口直径的增大,微火焰的上方出现回流; Re=12条件下,在喷口直径=0.07 mm时存在熄灭极限;稳定燃烧区的最小发热率约为0.5 W;微尺度条件下,Da数对火焰结构和火焰的熄灭有一定的影响。     

On the existence of steady-state gaseous microgravity spherical diffusion flames in the presence of radiation heat loss

Whether steady-state gaseous microgravity spherical diffusion exist in the presence of radiation heat loss is an important fundamental question and has important implications for spacecraft fire safety. In this work, experiments aboard the International Space Station and a transient numerical model are used to investigate the existence of steady-state microgravity spherical diffusion flames. Gaseous spherical diffusion flames stabilized on a porous spherical burner are employed in normal (i.e., fuel flowing into an ambient oxidizer) and inverse (i.e., oxidizer flowing into an ambient fuel) flame configurations. The fuel is ethylene and the oxidizer oxygen, both diluted with nitrogen. The flow rate of the reactant gas from the burner is held constant. It is found that steady-state gaseous microgravity spherical diffusion flames can exist in the presence of radiation heat loss, provided that the steady-state flame size is less than the flame size for radiative extinction, and the flame develops fast enough that radiation heat loss does not drop the flame temperature below the critical temperature for radiative extinction (1130 K). A simple model is provided that allows for the identification of initial conditions that can lead to steady-state spherical diffusion flames. In the spherical, infinite domain configuration, the characteristic time for the diffusion-controlled system to effectively reach steady-state is found to be on the order of 100,000 s. Despite a narrow range of attainable conditions, flames that exhibit steady-state behavior are observed aboard the ISS for up to 870 s, even with the constraint of a finite boundary. Steady-state flames are simulated using the numerical model for over 100,000 s.     

Influence of Burner Geometry on Heat Transfer Characteristics of Methane/Air Flame Impinging on Flat Surface

An experimental study has been conducted to find the heat transfer characteristics of methane/air flames impinging normally to a flat surface using different burner geometries. The burners used were of nozzle, tube, and orifice type each with a diameter of 10 mm. Due to different exit velocity profiles, the flame structures were different in each case. Because of nearly flat velocity profile, the flame spread was more in case of orifice and nozzle burners as compared to tube burner. Effects of varying the value of Reynolds number (600–2500), equivalence ratio (0.8–1.5) and dimensionless separation distance (0.7–8) on heat transfer characteristics on the flat plate have been investigated for the tube burner. Different flame shapes were observed for different impingement conditions. It has been observed that the heat transfer characteristics were intimately related to flame shapes. Heat transfer characteristics were discussed for the cases when the flame inner reaction cone was far away, just touched, and was intercepted by the plate. Negative heat fluxes at the stagnation point were observed when the inner reaction cone was intercepted by the plate due to impingement of cool un-burnt mixture directly on the surface. Different heat transfer characteristics were observed for different burner geometries with similar operating conditions. In case of tube burner, the maximum heat flux is around the stagnation point and decay is faster in the radial direction. In case of nozzle and orifice burner, the heat transfer distribution is more uniform over the surface.     

Laminar flame propagation on a horizontal fuel surface: Verification of classical Emmons solution

This work analyses the classical Emmons (1956) solution of flat plate laminar flame combustion on a film of liquid fuel. A two-dimensional (2D) numerical model developed for this purpose has been benchmarked with experimental results available in the literature for methanol. In the parametric study, numerical predictions have been compared with Emmons classical solution. The study shows that the Emmons solution is valid in a range of Reynolds numbers where flame anchors near the leading edge of the methanol pool and the combustion zone is confined around the hydrodynamic and thermal boundary layers. However, in cases of low free stream velocities the combustion zone is beyond the boundary layer zone and the Emmons solution deviates. In cases of very high free stream velocities, the flame moves away from the leading edge and anchors at a location downstream. The Emmons solution is not applicable in this case as well. For the fuel considered in this study (methanol), accounting for thermal radiation, employing an optically thin radiation model, allows better agreement between experimental and numerical temperature profiles but does not affect the mass burning rates.     

Dilution effects on laminar jet diffusion flame lengths

Many studies have examined the stoichiometric lengths of laminar gas jet diffusion flames. However, these have emphasized normal flames of undiluted fuel burning in air. Many questions remain about the effects of fuel dilution, oxygen-enhanced combustion, and inverse flames. Thus, the stoichiometric lengths of 287 normal and inverse gas jet flames are measured for a broad range of nitrogen dilution. The fuels are methane and propane and the ambient pressure is atmospheric. Nitrogen addition to the fuel and/or oxidizer is found to increase the stoichiometric lengths of both normal and inverse diffusion flames, but this effect is small at high reactant mole fraction. This counters previous assertions that inert addition to the fuel stream has a negligible effect on the lengths of normal diffusion flames. The analytical model of Roper is extended to these conditions by specifying the characteristic diffusivity to be the mean diffusivity of the fuel and oxidizer into stoichiometric products and a characteristic temperature that scales with the adiabatic flame temperature and the ambient temperature. The extended model correlates the measured lengths of normal and inverse flames with coefficients of determination of 0.87 for methane and 0.97 for propane.     

Experimental low-stretch gaseous diffusion flames in buoyancy-induced flowfields

The present study experimentally investigates the structure and instabilities associated with extremely low-stretch (1 s⁻¹) gaseous diffusion flames. Ultra-low-stretch flames are established in normal gravity by bottom burning of a methane/nitrogen mixture discharged from a porous spherically symmetric burner of large radius of curvature. OH-PLIF and IR imaging techniques are used to characterize the reaction zone and the burner surface temperature, respectively. A flame stability diagram mapping the response of the ultra-low-stretch diffusion flame to varying fuel injection rate and nitrogen dilution is explored. In this diagram, two main boundaries are identified. These boundaries separate the stability diagram into three regions: sooting flame, non-sooting flame, and extinction. Two distinct extinction mechanisms are noted. For low fuel injection rates, flame extinction is caused by heat loss to the burner surface. For relatively high injection rates, at which the heat loss to burner surface is negligible, flame radiative heat loss is the dominant extinction mechanism. There also exists a critical inert dilution level beyond which the flame cannot be sustained. The existence of multi-dimensional flame phenomena near the extinction limits is also identified. Various multi-dimensional flame patterns are observed, and their evolutions are studied using direct chemiluminescence and OH-PLIF imaging. The results demonstrate the usefulness of the present burner configuration for the study of low-stretch gaseous diffusion flames.     

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.     

Combustion dynamics of inverted conical flames

An inverted conical flame anchored on a central bluff-body in an unconfined burner configuration features a distinctive acoustic response. This configuration typifies more complex situations in which the thermo-acoustic instability is driven by the interaction of a flame with a convective vorticity mode. The axisymmetric geometry investigated in this article features a shear region between the reactive jet and the surrounding atmosphere. It exhibits self-sustained oscillations for certain operating conditions involving a powerful flame collapse phenomenon with sudden annihilation of flame surface area. This is caused by a strong interaction between the flame and vortices created in the outer jet shear layer, a process which determines the amplitude of heat release fluctuation and its time delay with respect to incident velocity perturbations. This process also generates an acoustic field that excites the burner and synchronizes the vortex shedding mechanism. The transfer functions between the velocity signal at the burner outlet and heat release are obtained experimentally for a set of flow velocities fluctuations levels. It is found that heat release fluctuations are a strong function of the incoming velocity perturbation amplitude and that the time delay between these two quantities is mainly determined by the convection of the large scale vortices formed in the jet shear layer. A model is formulated, which suitably describes the observed instabilities.     

Measurements of local mixture fraction of reacting mixture in swirl-stabilised natural gas-fuelled burners

Local, time-dependent measurements of mixture fraction of the reacting mixture were obtained in a swirl-stabilised natural gas-fuelled, nominally non-premixed burner using the intensity of chemiluminescence from OH^∗ and CH^∗ radicals. The measurements quantified the mean, rms of fluctuations and probability density functions of local mixture fraction at the stabilisation region of the flame. In addition, the probability of flame presence and the degree of lean or rich versus stoichiometric reaction is reported. The burner was operated for three air flow Reynolds numbers (Re=18970, 29100 and 57600), at an overall equivalence ratio of 0.32, without and with imposed oscillations to the air flow of the burner at the resonance frequency of 350 Hz. Results show that combustion occurred in a partially premixed mode for all flow conditions, although fuel and air were injected separately in the reaction zone. The mean local mixture fraction was nearly stoichiometric at the base of the flame without imposed air oscillations, but with large fluctuations leading to around 80% of lean or rich reaction. The degree of non-stoichiometric reaction increased with axial distance from the burner exit and Reynolds number and lean reaction dominated. Imposed air oscillations led to lifted flames and increased the degree of non-stoichiometric reaction for Re=18970 and 29100, whereas the flame remained attached onto the injector for Re=57600 and little modification of the mixture fraction was observed.