On the oscillating flame characteristics in nonpremixed laminar coflow-jets: An experimental and numerical study

This study investigates the characteristics of oscillating lifted flames in laminar coflow-jets experimentally and numerically by varying both fuel density (by varying propane and n-butane mixtures) and coflow density (by diluting air with N₂/He mixtures). Two different lifted flame oscillation behaviors are observed depending on these parameters: oscillating tribrachial lifted flame (OTLF) and oscillating mode-change lifted flame (OMLF), where a rapid increase in flame radius is observed. The regimes of the two flames are identified from experiments, which shows that OMLF occurs only when the effect of the negative buoyancy on the flow field by the fuel heavier than air becomes significant at low fuel jet velocity. OMLFs are also identified to distinguish OTLF regime from flame extinction, which implies that an OMLF can be extinguished when the positive buoyancy becomes weak, losing its stabilizing effect, or when the negative buoyancy becomes strong, further enhancing its destabilizing effect. Transient numerical simulations of both OTLF and OMLF reveal that the OMLF occurs by a strong toroidal vortex and a subsequent counterflow-like structure induced by relatively-strong negative buoyancy. Such a drastic flow redirection significantly changes the fuel concentration gradient such that the OMLF changes its mode from a tribrachial flame mode (decreasing edge speed with fuel concentration gradient) to the premixed flame-like transition mode when the fuel concentration gradient becomes very small (increasing edge speed with fuel concentration gradient). Again, a tribrachial flame mode is recovered during a rising period of flame edge and repeats an oscillation cycle.     

Structure of propagating edge diffusion flames in hydrocarbon fuel jets

Direct numerical simulations with a C₃-chemistry model have been performed to investigate the transient behavior and internal structure of flames propagating in an axisymmetric fuel jet of methane, ethane, ethylene, acetylene, or propane in normal earth gravity (1g) and zero gravity (0g). The fuel issued from a 3-mm-i.d. tube into quasi-quiescent air for a fixed mixing time of 0.3 s before it was ignited along the centerline where the fuel–air mixture was at stoichiometry. The edge of the flame formed a vigorously burning peak reactivity spot, i.e., reaction kernel, and propagated through a flammable mixture layer, leaving behind a trailing diffusion flame. The reaction kernel broadened laterally across the flammable mixture layer and possessed characteristics of premixed flames in the direction of propagation and unique flame structure in the transverse direction. The reaction kernel grew wings on both fuel and air sides to form a triple-flame-like structure, particularly for ethylene and acetylene, whereas for alkanes, the fuel-rich wing tended to merge with the main diffusion flame zone, particularly methane. The topology of edge diffusion flames depend on the properties of fuels, particularly the rich flammability limit, and the mechanistic oxidation pathways. The transit velocity of edge diffusion flames, determined from a time series of calculated temperature field, equaled to the measured laminar flame speed of the stoichiometric fuel–air mixtures, available in the literature, independent of the gravity level.     

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.     

Effects of Soret diffusion on the laminar flame speed and Markstein length of syngas/air mixtures

The effects of Soret diffusion on premixed syngas/air flames at normal and elevated temperatures and pressures are investigated numerically including detailed chemistry and transport. The emphasis is placed on assessing and interpreting the influence of Soret diffusion on the unstretched and stretched laminar flame speed and Markstein length of syngas/air mixtures. The laminar flame speed and Markstein length are obtained by simulating the unstretched planar flame and positively-stretched spherical flame, respectively. The results indicate that at atmospheric pressure the laminar flame speed of syngas/air is mainly reduced by Soret diffusion of H radical while the influence of H₂ Soret diffusion is negligible. This is due to the facts that the main reaction zone and the Soret diffusion for H radical (H₂) are strongly (weakly) coupled, and that Soret diffusion reduces the H concentration in the reaction zone. Because of the enhancement in the Soret diffusion flux of H radical, the influence of Soret diffusion on the laminar burning flux increases with the initial temperature and pressure. Unlike the results at atmospheric pressure, at elevated pressures the laminar flame speed is shown to be affected by the Soret diffusion of H₂ as well as H radical. For stretched spherical flame, it is shown that the Soret diffusion of both H and H₂ should be included so that the stretched flame speed can be accurately predicted. Similar to the laminar flame speed, the Markstein length is also reduced by Soret diffusion. However, the reduction is found to be mainly caused by Soret diffusion of H₂ rather than that of H radical. Moreover, the influence of Soret diffusion on the Markstein length is demonstrated to decrease with the initial temperature and pressure.     

Visualization of dilute hydrogen jet flame in air flow

The structure of hydrogen jet flame diluted by CO₂ in air flow is studied by various visualization techniques, such as schlieren, direct photograph, tracer injection and reactive Mie scattering method, which allow understanding of the influence of CO₂ on the characteristics of the hydrogen jet flame. The experimental result indicates that the flame structure consists of laminar fuel jet and surrounding reaction zone near the nozzle exit. When the CO2 fraction is increased, the width of the fuel jet grows and the reaction zone is reduced in size. These observations are further confirmed by quantitative measurements of temperature and velocity fields in the flame, which are evaluated by thermocouple and particle image velocimetry (PIV), respectively. These results indicate that the flame temperature is decreased and the flow rate of the fuel jet is increased by the influence of diluents, which are due to the reduced calorific value and larger density of fuel, respectively.     

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

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

A numerical study of quasi-steady burning characteristics of a condensed fuel: effect of angular orientation of fuel surface

A numerical study of laminar diffusion flames established over a condensed fuel surface, inclined at several angular orientations in the range of –90°?θ?+90° with respect to the vertical axis, under atmospheric pressure and normal gravity environment, is presented. Methanol is employed as the fuel. A numerical model, which solves transient gas-phase, two-dimensional governing conservation equations, with a single-step global reaction for methanol–air oxidation and an optically thin radiation sub-model, has been employed in the present investigation. Numerical results have been validated against the experimental data from the present study. Thereafter, the model is used to investigate the influence of angular orientation of fuel surface on its quasi-steady burning characteristics. Results in terms of fuel mass burning rate, flame stand-off distances, temperature field, velocity profiles and oxygen contours have been presented and discussed in detail. It is observed that orientation angles in the range of –45°?θ? –30° (fuel surface facing upwards), yield the maximum mass burning rates. The flame anchoring location near the leading edge of the fuel surface, normal gradient of fuel vapor mass fraction at the surface and oxygen contours have been used to explore this unique behavior. Based on the numerical results, a theoretical correlation to predict the mass burning rate as a function of fuel surface orientation is also proposed. Furthermore, a discussion on the differences in the structure of laminar diffusion flame established over fuel surface as a function of its angular orientation is included.     

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.     

反扩散火焰温度场特性的实验研究

本文设计搭建石英管反扩散层流燃烧实验系统,研究了不同雷诺数和燃料组分条件下反扩散火焰的温度场特性。实验表明,反扩散火焰中心某一高度存在燃料和中心空气快速混合区域,体现出一定的预混燃烧特性。随中心射流雷诺数增大,火焰底部冷核区增大,高温区位置沿轴线上移,但火焰最高温度点的水平高度先升高后降低。不同燃料组分下,反扩散火焰的整体温度有所变化,但火焰温度的分布特性基本相似。     

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.     

The effect of fuel inlet turbulence intensity on H2/CH4 flame structure of MILD combustion using the LES method

Large eddy simulations (LES) are employed to investigate the effect of the inlet turbulence intensity on the H₂/CH₄ flame structure in a hot and diluted co-flow stream which emulates the (Moderate or Intense Low-oxygen Dilution) MILD combustion regime. In this regard, three fuel inlet turbulence intensity profiles with the values of 4%, 7% and 10% are superimposed on the annular mixing layer. The effects of these changes on the flame structure under the MILD condition are studied for two oxygen concentrations of 3% and 9% (by mass) in the oxidiser stream and three hot co-flow temperatures 1300, 1500 and 1750 K. The turbulence-chemistry interaction of the numerically unresolved scales is modelled using the (Partially Stirred Reactor) PaSR method, where the full mechanism of GRI-2.11 represents the chemical reactions. The influences of the turbulence intensity on the flame structure under the MILD condition are studied by using the profile of temperature, CO and OH mass fractions in both physical and mixture fraction spaces at two downstream locations. Also, the effects of this parameter are investigated by contours of OH, HCO and CH₂O radicals in an area near the nozzle exit zone. Results show that increasing the fuel inlet turbulence intensity has a profound effect on the flame structure particularly at low oxygen mass fraction. This increment weakens the combustion zone and results in a decrease in the peak values of the flame temperature and OH and CO mass fractions. Furthermore, increasing the inlet turbulence intensity decreases the flame thickness, and increases the MILD flame instability and diffusion of un-burnt fuel through the flame front. These effects are reduced by increasing the hot co-flow temperature which reinforces the reaction zone.     

Diffusion flame instabilities in a 0.75 mm non-premixed microburner

We examine the cellular instabilities of laminar non-premixed diffusion flames that arise in a polycrystalline alumina microburner with a channel wall gap of dimension 0.75 mm. Changes in the flame structure are observed as a function of the fuel type (H₂, CH₄, and C₃H₈) and diluent. The oxidizer is O₂/inert. In contrast to previous observations on laminar diffusion flame instabilities, the current instabilities occur in the direction of flow above the splitter plate, and only occur for the heavier fuel types. They are not observed in a H₂–O₂ mixture, which will only support a continuous laminar flame inside our burner, regardless of the initial mixture strength and whether or not the flame is in near-quenching conditions. The only exception is when helium is added to the H₂–O₂ mixture, raising the effective Lewis numbers of both components.     

A spray flamelet/progress variable approach combined with a transported joint PDF model for turbulent spray flames

A spray flamelet/progress variable approach is developed for use in spray combustion with partly pre-vaporised liquid fuel, where a laminar spray flamelet library accounts for evaporation within the laminar flame structures. For this purpose, the standard spray flamelet formulation for pure evaporating liquid fuel and oxidiser is extended by a chemical reaction progress variable in both the turbulent spray flame model and the laminar spray flame structures, in order to account for the effect of pre-vaporised liquid fuel for instance through use of a pilot flame. This new approach is combined with a transported joint probability density function (PDF) method for the simulation of a turbulent piloted ethanol/air spray flame, and the extension requires the formulation of a joint three-variate PDF depending on the gas phase mixture fraction, the chemical reaction progress variable, and gas enthalpy. The molecular mixing is modelled with the extended interaction-by-exchange-with-the-mean (IEM) model, where source terms account for spray evaporation and heat exchange due to evaporation as well as the chemical reaction rate for the chemical reaction progress variable. This is the first formulation using a spray flamelet model considering both evaporation and partly pre-vaporised liquid fuel within the laminar spray flamelets. Results with this new formulation show good agreement with the experimental data provided by A.R. Masri, Sydney, Australia. The analysis of the Lagrangian statistics of the gas temperature and the OH mass fraction indicates that partially premixed combustion prevails near the nozzle exit of the spray, whereas further downstream, the non-premixed flame is promoted towards the inner rich-side of the spray jet since the pilot flame heats up the premixed inner spray zone. In summary, the simulation with the new formulation considering the reaction progress variable shows good performance, greatly improving the standard formulation, and it provides new insight into the local structure of this complex spray flame.     

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.     

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.     

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.     

Control of oscillating combustion and noise based on local flame structure

To control combustion oscillations, the characteristics of an oscillating swirl injection premixed flame have been investigated, and control of oscillating combustion and noise based on local flame structure has been conducted. The r.m.s. value of pressure fluctuations and noise level show significantly large values between = 0.8 and 1.1. The beating of pressure fluctuations is observed for the large oscillating flame conditions in this combustor. Relationship between beating of pressure fluctuations and local flame structure was observed by the simultaneous measurement of CH/OH planar laser induced fluorescence and pressure fluctuations. The local flame structure and beating of pressure fluctuations are related and the most complicated flame is formed in the middle pressure fluctuating region of beating. The beating of pressure fluctuations, which plays important roles in noise generation and nitric oxide emission in this combustor, could be controlled by injecting secondary fuel into the recirculating region of oscillating flames. Injecting secondary fuel prevented lean blowout, and low NO_x combustion was also achieved even for the case of pure methane injection as a secondary fuel. By injecting secondary fuel into the recirculating region near the swirl injector, the flame lifted from the swirl injector and its reaction region became uniform and widespread, hence resulting in low nitric oxide emission. Secondary mixture injection, fuel diluted with air, is not effective for control of combustion oscillations suppression and lean blowout prevention.     

A comparison of computational and experimental lift-off heights of coflow laminar diffusion flames

As a sensitive marker of changes in flame structure, the number densities of excited-state CH (denoted CH^*), and excited-state OH (denoted OH^*) are imaged in coflow laminar diffusion flames. Measurements are made both in normal gravity and on the NASA KC-135 reduced-gravity aircraft. The spatial distribution of these radicals provides information about flame structure and lift-off heights that can be directly compared with computational predictions. Measurements and computations are compared over a range of buoyancy and fuel dilution levels. Results indicate that the lift-off heights and flame shapes predicted by the computations are in excellent agreement with measurement for both normal gravity (1g) and reduced gravity flames at low dilution levels. As the fuel mixture is increasingly diluted, however, the 1g lift-off heights become underpredicted. This trend continues until the computations predict stable flames at highly dilute fuel mixtures beyond the 1g experimental blow-off limit. To better understand this behavior, an analysis was performed, which indicates that the lift-off height is sensitive to the laminar flame speed of the corresponding premixed mixture at the flame edge. By varying the rates of two key “flame speed” controlling reactions, we were able to modify the predicted lift-off heights so as to be in closer agreement with the experiments. The results indicate that reaction sets that work well in low dilution systems may need to be modified to accommodate high dilution flames.     

Analysis of a spray diffusion flame and its extinction in a shear layer

A new semi-analytical solution for a laminar spray diffusion flame in the shear layer between fuel and oxidant streams is developed. The Stokes number is identified as a small spray droplet-related parameter to be used in a perturbation analysis of the liquid phase governing equations. Appropriate specification of an additional parameter ensures that similarity is achieved so that the concentration field of the liquid in the spray can be readily found. The coupled gas-phase equations are treated using the usual inverse of the large Zeldovitch number for the asymptotic analysis. Numerical results demonstrate how the distribution of the liquid phase in the developing shear layer between two unidirectional gas streams flowing over one another with (the possibility of) dissimilar velocities in their respective free-streams influences the flame shape, location, fuel vapour and temperature fields. An extinction analysis enables a parametric mapping of conditions for extinguishment of the spray diffusion flames to be drawn.     

Combustion of a high-velocity hydrogen microjet effluxing in air

This study is devoted to experimental investigation of hydrogen-combustion modes and the structure of a diffusion flame formed at a high-velocity efflux of hydrogen in air through round apertures of various diameters. The efflux-velocity range of the hydrogen jet and the diameters of nozzle apertures at which the flame is divided in two zones with laminar and turbulent flow are found. The zone with the laminar flow is a stabilizer of combustion of the flame as a whole, and in the zone with the turbulent flow the intense mixing of fuel with an oxidizer takes place. Combustion in these two zones can occur independently from each other, but the steadiest mode is observed only at the existence of the flame in the laminar-flow zone. The knowledge obtained makes it possible to understand more deeply the features of modes of microjet combustion of hydrogen promising for various combustion devices.