预混火焰拉伸和曲率效率的物理分析

湍流燃烧的基本火焰结构是拉伸的曲面涡管;拉伸流场中的管形火焰模拟了湍流燃烧的细微结构。本文对平面预混火焰、拉伸预混火焰和管形拉伸预混火焰进行了质量、能量和组分的守恒分析。通过对比这几种火焰,揭示了火焰拉伸效果是通过优先扩散改变火焰温度和熄火极限;而火焰曲率通过增强或削弱这种优先扩散效果来影响火焰温度,影响的大小和火焰厚度与火焰半径的比值呈正比。     

CH4/O2/N2层流火焰瞬态响应特性数值分析

用数值分析方法研究CH₄/O₂/N₂层流扩散火焰的瞬态响应特性.采用详细的GRI-Mech 3.0机理(包含53种组分,325个基元反应)描述CH₄氧化和NO_x生成.首先比较火焰面稳态结构的计算结果和实验数据,以验证数值方法的可靠性.用台阶跃变的火焰拉伸率来模拟瞬态流场对火焰面局部结构的影响,给出了火焰面结构(温度、组分浓度)的瞬态响应曲线,分析了火焰面的响应特性.着重探讨了不同拉伸率跃变幅度对响应特性的影响,发现火焰面的响应对于拉伸率正向跃变和负向跃变并不对称,而是相反,且在小的拉伸率跃变范围内火焰面响应时间和拉伸率跃变幅度近似成反比关系.另外,温度的平均响应时间远大于一个典型湍流燃烧场的流动时间尺度,说明火焰面非稳态效应对于湍流燃烧数值模拟有重要意义.     

浮力扩散火焰对低频外加声场的响应特性

本文通过实验研究了浮力扩散火焰在外加声场作用下的响应行为和熄火机理.实验发现,声场中的气体燃料火焰熄火时,火焰处的声压与声频呈正相关.但该正比关系受火焰和环境等因素影响,无法精确判断熄火的临界条件.因此,本文基于火焰在声场中的偏移量脉动,为声场的临界灭火条件构造了数值判据.研究结果表明,该判据对目标火焰的依赖性较低,且...     

层流拉伸非预混火焰影响因素的数值研究

采用简化机理(19个基元反应,14种组分)和半详细机理(79个基元反应,32种组分)对层流拉伸非预混火焰的结构进行了数值模拟,通过与实验数据的对比发现:在由当量混合的标量耗散率表征的火焰拉伸一定时,考虑不同分子扩散系数和不同的化学反应机理均对火焰面结构的准确模拟有着重要影响.此外计算了考虑热辐射时,随着当量混合的标量耗散率的改变而变化的每单位火焰面面积放热速率和燃烧场的最高温度,并分析了层流非预混火焰的熄火现象.     

湍流扩散火焰的非稳态火焰面模拟

采用稳态的和非稳态的火焰面模型同时对一个湍流甲烷射流扩散火焰进行了数值模拟,比较了两者对湍流平均火焰结构、活性自由基和污染物(氮氧化物)排放的模拟效果。速度场采用κ-ε模型计算,守恒标量混合物分数的分布通过其概率密度函数(PDF)输运方程的求解得到。稳态的火焰面结构由查询火焰面数据库得到,而非稳态的火焰面结构由火焰面方程和流场方程耦合求解来计算。采用详细的GRI—Mech 3.0机理描述甲烷的氧化和氮氧化物的形成。数值模拟结果和实验数据作了广泛的对比,验证了火焰面模型对湍流扩散燃烧的定量模拟能力。     

稀释H2/CO混合气与空气对撞燃烧的数值研究

采用多组分混合物中扩散速度表达式并考虑Soret效应,在物理空间求解同时加入了多种稀释成分(CO2,H2O,N2)的H2/CO混合气与空气对撞燃烧形成的非预混火焰;考察了稀释成分加入与否以及稀释成分浓度的调整对于火焰结构以及相应熄火拉伸率产生的影响.结合数值计算中采用的化学反应机理,从稀释给H2/CO混合气所带来的物理...     

本生灯火焰和Smithll火焰的实验研究

一、引言 本生灯火焰顶部开口现象实际上是火焰面大曲率变化处的火焰局部熄火。除了化学反应动力学性质引起的火焰局部熄火外,对局部熄火机理的理解主要有两种。一种认为是由流场不均匀造成的气动拉伸机理;另一种认为是燃料和氧化剂分层扩散的机理。最初对本生灯火焰顶部开口现象的解释是认为二种机理都在起作用。后来的研究指     

钝体驻定湍流扩散火焰的数值研究——燃烧模型比较

分别采用标量联合的概率密度函数方法、稳态火焰面模型、Euler非稳态火焰面模型和基于有限体积/Monte Carlo混合算法的完备PDF模型对钝体驻定的Sydney湍流扩散火焰HM1进行数值模拟,以比较不同燃烧模型的性能,并比较标量联合的概率密度函数方法和Euler非稳态火焰面模型对氮氧化物排放预测的差异.计算结果和实验数据的比较表明,采用概率密度函数方法计算化学反应可以得到更好的结果但计算量较大,而用火焰面模型求解计算量较小,在接近完全燃烧的情形下,其计算结果比较合理.     

碳氢扩散火焰实验研究：燃烧气氛对火焰光谱与温度的影响

辐射是各种燃烧过程中热传递的主要方式。在不同的火焰中,辐射光谱分布十分复杂。在这项工作中,利用光谱仪测量了可见光(200～900 nm),近红外(900～1 700 nm)和中红外(2 500～5 000 nm)波段火焰的光谱强度,分析了空气和富氧气氛下扩散火焰的光谱特征。并基于光谱分析,定量得到了火焰中碳烟以及气体发射的辐射力,计算了火焰的温度分布。结果表明,空气燃烧中的火焰温度低于富氧燃烧中的火焰温度。在空气气氛下,火焰中的碳烟和气体均对中的热辐射起着重要作用。而在富氧气氛下,气体对于火焰热辐射更为重要。在可见光和近红外波段,由于在空气气氛下火焰中碳烟的大量形成,光谱曲线显示出了良好连续性。而富氧气氛下火焰的辐射光谱降低。在中红外波段,空气气氛下火焰的气体辐射明显弱于富氧气氛下火焰的气体辐射。     

PDF方法模拟钝体驻定的湍流扩散火焰

采用标量联合的概率密度函数方法,对钝体驻定的湍流射流扩散Sydney火焰HM1进行数值模拟,结合当地自适应建表方法加速化学反应计算,用修正的LRR-IP雷诺应力模型求解速度场.首次对3种不同规模的甲烷化学反应动力学机理进行研究,并与实验数据进行比较,结果表明,模型和反应机理很好地预测了速度场和标量场的变化及局部熄火现象,而考虑反应机理中的C2化学对火焰HM1的影响不大.     

Bifurcation and extinction limit of stretched premixed flames with chain-branching intermediate kinetics and radiative loss

Premixed counterflow flames with thermally sensitive intermediate kinetics and radiation heat loss are analysed within the framework of large activation energy. Unlike previous studies considering one-step global reaction, two-step chemistry consisting of a chain branching reaction and a recombination reaction is considered here. The correlation between the flame front location and stretch rate is derived. Based on this correlation, the extinction limit and bifurcation characteristics of the strained premixed flame are studied, and the effects of fuel and radical Lewis numbers as well as radiation heat loss are examined. Different flame regimes and their extinction characteristics can be predicted by the present theory. It is found that fuel Lewis number affects the flame bifurcation qualitatively and quantitatively, whereas radical Lewis number only has a quantitative influence. Stretch rates at the stretch and radiation extinction limits respectively decrease and increase with fuel Lewis number before the flammability limit is reached, while the radical Lewis number shows the opposite tendency. In addition, the relation between the standard flammability limit and the limit derived from the strained near stagnation flame is affected by the fuel Lewis number, but not by the radical Lewis number. Meanwhile, the flammability limit increases with decreased fuel Lewis number, but with increased radical Lewis number. Radical behaviours at flame front corresponding to flame bifurcation and extinction are also analysed in this work. It is shown that radical concentration at the flame front, under extinction stretch rate condition, increases with radical Lewis number but decreases with fuel Lewis number. It decreases with increased radiation loss.     

A LES-CMC formulation for premixed flames including differential diffusion

A finite volume large eddy simulation–conditional moment closure (LES-CMC) numerical framework for premixed combustion developed in a previous studyhas been extended to account for differential diffusion. The non-unity Lewis number CMC transport equation has an additional convective term in sample space proportional to the conditional diffusion of the progress variable, that in turn accounts for diffusion normal to the flame front and curvature-induced effects. Planar laminar simulations are first performed using a spatially homogeneous non-unity Lewis number CMC formulation and validated against physical-space fully resolved reference solutions. The same CMC formulation is subsequently used to numerically investigate the effects of curvature for laminar flames having different effective Lewis numbers: a lean methane–air flame with Le_eff = 0.99 and a lean hydrogen–air flame with Le_eff = 0.33. Results suggest that curvature does not affect the conditional heat release if the effective Lewis number tends to unity, so that curvature-induced transport may be neglected. Finally, the effect of turbulence on the flame structure is qualitatively analysed using LES-CMC simulations with and without differential diffusion for a turbulent premixed bluff body methane–air flame exhibiting local extinction behaviour. Overall, both the unity and the non-unity computations predict the characteristic M-shaped flame observed experimentally, although some minor differences are identified. The findings suggest that for the high Karlovitz number (from 1 to 10) flame considered, turbulent mixing within the flame weakens the differential transport contribution by reducing the conditional scalar dissipation rate and accordingly the conditional diffusion of the progress variable.     

A review of chemical diffusion: Criticism and limits of simplified methods for diffusion coefficient calculation

A review of the physics and modelling of mass diffusion involving different gaseous chemical species is firstly proposed. Both accurate and simplified models for mass diffusion involve the calculation of individual species diffusion coefficients. Since these are computationally expensive, in CFD they are commonly estimated by assuming constant Lewis or Schmidt numbers for each chemical species. The constant Lewis number assumption is particularly used. As a matter of fact, these assumptions have never been theoretically justified nor verified in practical flames. The only published information are the first observations by Smooke and Giovangigli about the Lewis number against temperature distributions in methane–air premixed and counterflow diffusion one-dimensional flames. The aim of this work is to verify these assumptions. Functional dependences of molecular properties appearing in these numbers are made explicit to show that while Sc _i depends only on composition, Le _i depends also on temperature and therefore it certainly cannot be assumed constant in a flame. Then, accurately calculating molecular properties, distributions of these characteristic numbers against temperature are obtained a posteriori from numerical simulations of different flames, premixed and non-premixed, and burning different fuels. For non-premixed flames, individual species Lewis number distributions are broad for most of the species considered in this article, whilst they are tight for premixed flames. Some attention is focused on the particular shape of Lewis distributions in non-premixed flames: they are characterized by four or five (when extinction is experienced) branches associated to precise regions in the flame (basically, lean, rich and stoichiometric combusting zones). Instead, the Schmidt distributions are always tighter, also when extinctions take place: for many species they can be approximatively assumed constant. Finally, a simplified procedure to estimate individual species diffusion coefficients is suggested, assuming the median of non-premixed flame Schmidt distributions has a constant value for each chemical species.     

Structure and dynamic response of radiative diffusion flames

Nitrogen-diluted hydrogen burning in air is modeled numerically using a constant density and one-step reaction model in a plane two-dimensional counterflow configuration. An optically thin assumption is used to investigate the effects of radiation on the dynamics, structure, and extinction of diffusion flames. While there exist dual steady-state extinction limits for the 1D radiative flame response, it is found that as the 1D radiative extinction point is approached the 1D low-stretch diffusion flame exhibits oscillatory response, even with sub-unity Lewis number fuel. These radiation-induced limit cycle oscillations are found to have increasing amplitude and decreasing frequency as the stretch rate is reduced. Flame oscillation eventually leads to permanent extinction at the stretch rate which is larger than the steady-state radiative extinction value. Along the 1D radiative response curve, the transition from 1D flame to 2D structure and the differences in the resulting 2D flame patterns are also examined using a variety of initial profiles, with special emphasis on the comparison of using the initial profiles with and without a flame edge. Similar to the previous studies on the high-stretch adiabatic edge flames using the same configuration, the high-stretch radiative flames are found to resist 1D blow-off quenching through various 2D structures, including propagating front and steady cellular flames for initial profiles with and without flame edges. For all initial profiles studied, the low-stretch radiative flames are also found to exhibit different 2D flame phenomena near the 1D radiative extinction limit, such as transient cellular structures, steady cellular structures, and pulsating ignition fronts. Although the results demonstrate the presence of low-stretch and high-stretch 2D bifurcation branches close to the corresponding 1D extinction limits irrespective of the initial profile used, particular 2D flame structures in certain stretch rate range are initial profile dependent. The existence of two-dimensional flame structures beyond the 1D steady-state radiative extinction limit suggests that the flammable range is expanded as compared to that predicted by the 1D model. Hence, multi-dimensional flame patterns need to be accounted for when determining the flammability limits for a given system.     

Diffusional-thermal instability in strained diffusion flames with unequal Lewis numbers

Diffusional–thermal instability is analysed for near-extinction counterflow diffusion flames to examine the instability characteristics of strained diffusion flamelets in turbulent flames, with the additional intention of providing a guideline to future experimental investigations. Attention is focused on the linear stability of the instability patterns appearing in the unstrained direction of two-dimensional counterflow diffusion flames, which is treated by the near-equilibrium regime of activation-energy asymptotics with Lewis numbers close to unity. The effects of unequal Lewis numbers for fuel and oxidizer are also taken into account by introducing an effective Lewis number. The resulting formulation describing linear stability of the harmonically decomposed disturbances turns out to be identical to the formulation derived previously for equal fuel and oxidizer Lewis numbers. For effective Lewis numbers less than unity, cellular instability is predicted for the entire range of the equivalence ratio, and the threshold Lewis number maintains a value slightly less than unity. On the other hand, for effective Lewis numbers sufficiently greater than unity, two types of oscillatory instabilities are found. As the effective Lewis number increases from unity, a travelling instability is first encountered for a range of finite wavelengths, and a pulsating instability emerges immediately above the travelling instability. These two types of oscillatory instabilities are predicted only for equivalence ratios sufficiently greater than unity because the threshold Lewis numbers for these instabilities are found to be infinity at unity equivalence ratio. For large values of the equivalence ratio, which is typical of most hydrocarbon flames, oscillatory instabilities are predicted only for flames burning extremely heavy hydrocarbon fuels or for flames heavily diluted by light inert gases.     

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.     

旋流燃烧器复合小火焰模型的大涡模拟

对低旋燃烧器内的流动和燃烧进行了大涡模拟,其中化学反应分别采用传统的扩散小火焰模型和所提出的复合小火焰模型描述。复合小火焰模型借助于燃烧区索引的概念区分当地的燃烧模式,若当地的燃烧为扩散燃烧,则调用扩散小火焰库,否则调用预混小火焰库。数值结果与实验数据的对比表明,采用大涡模拟方法结合小火焰模型能够很好地模拟实验室尺度的低旋燃烧器,且采用复合小火焰模型能够得到更加符合实际的结果。     

Homoclinic bifurcations in radiating diffusion flames

We analyse the dynamics of a model describing a planar diffusion flame with radiative heat losses incorporating a single step kinetic using timestepping techniques for Lewis number equal to one. We construct the full bifurcation diagram with respect to the Damköhler number including the branches of oscillating solutions. Based on this analysis we found, for the first time, homoclinic bifurcations that mark the abrupt disappearance of the nonlinear oscillations near extinction as reported in experiments.     

The effect of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame

The influence of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame was investigated by numerical simulation using three different transport property calculation methods. One simulation included preferential diffusion and the other two neglected preferential diffusion. The results show that the neglect of preferential diffusion or the use of unity Lewis number for all species results in a significant underprediction of soot volume fraction. The peak soot volume fraction is reduced from 8.0 to 2.0 ppm for the studied flame when preferential diffusion is neglected in the simulation. Detailed examination of numerical results reveals that the underprediction of soot volume fraction in the simulation neglecting preferential diffusion is due to the slower diffusion of some species from main reaction zone to PAH and soot formation layer. The slower diffusion of these species causes lower PAH formation rate and thus results in lower soot inception rate and smaller particle surface area. The smaller surface area further leads to smaller surface growth rate. In addition, the neglect of preferential diffusion also leads to higher OH concentration in the flame, which causes the higher specific soot oxidation rate. The lower inception rate, smaller surface growth rate and higher specific oxidation rate results in the lower soot volume fraction when preferential diffusion is neglected. The finding of the paper implies the importance of preferential diffusion for the modeling of not only laminar but maybe also some turbulent flames.     

Prediction of local extinction and re-ignition effects in non-premixed turbulent combustion using a flamelet/progress variable approach

The flamelet/progress variable approach (FPVA) has been proposed by Pierce and Moin as a model for turbulent non-premixed combustion in large-eddy simulation. The filtered chemical source term in this model appears in unclosed form, and is modeled by a presumed probability density function (PDF) for the joint PDF of the mixture fraction Z and a flamelet parameter λ. While the marginal PDF of Z can be reasonably approximated by a beta distribution, a model for the conditional PDF of the flamelet parameter needs to be developed. Further, the ability of FPVA to predict extinction and re-ignition has also not been assessed. In this paper, we address these aspects of the model using the DNS database of Sripakagorn et al. It is first shown that the steady flamelet assumption in the context of FPVA leads to good predictions even for high levels of local extinction. Three different models for the conditional PDF of the flamelet parameter are tested in an a priori sense. Results obtained using a delta function to model the conditional PDF of λ lead to an overprediction of the mean temperature, even with only moderate extinction levels. It is shown that if the conditional PDF of λ is modeled by a beta distribution conditioned on Z, then FPVA can predict extinction and re-ignition effects, and good agreement between the model and DNS data for the mean temperature is observed.