甲烷/氧气层流同轴射流扩散火焰OH*自由基的数值研究

OH*自由基是火焰中主要的激发态自由基之一,它所产生的化学发光可用于描述火焰的结构、拉伸率、氧燃当量比和热释放速率等特征信息,因此被广泛应用于火焰燃烧状态的在线诊断。以甲烷/氧气层流同轴射流扩散火焰作为研究对象,采用GRI-Mech 3.0机理结合OH*自由基生成和淬灭反应进行数值计算,对OH*自由基的二维分布特性进行研究,分析不同区域内OH*自由基的生成路径,并探讨不同氧燃当量比例和不同喷嘴出口尺寸对OH*自由基强度和分布特性的影响。模拟结果与实验研究基本吻合,表明计算模型能够准确描述火焰中OH*自由基的二维分布。结果表明：在甲烷/氧气层流同轴射流扩散火焰中,OH*自由基存在两种不同形态的分布区域,分别由反应CH+O₂=OH*+CO和H+O+M=OH*+M生成;随着氧燃当量比提高,OH*自由基的分布区域逐渐向火焰下游扩张,根据其分布形态的变化可以对火焰燃烧状况进行判断;如果OH*自由基仅分布于火焰的上游区域且呈断开形态,则说明火焰处于贫氧燃烧状态。如果OH*分布呈环状形态,则说明火焰处于富氧燃烧状态;相同氧气流量条件下,缩小喷嘴出口的环隙尺寸有助于加强燃料和氧气的化学反应程度,从而使火焰中OH*自由基的摩尔分数显著提高,增强OH*化学发光的辐射强度,提高火焰光谱诊断的准确性。     

热氧喷嘴水煤浆扩散火焰辐射光谱特性研究

火焰的辐射光谱可为燃烧诊断提供诸多信息,因此目前对简单的气态火焰自由基辐射特性已进行了大量研究,而关于非均相火焰的辐射光谱特性研究则相对较少。采用改进的热氧喷嘴技术在敞开空间下直接点燃水煤浆,并利用光纤光谱仪和紫外成像系统,着重对甲烷和水煤浆火焰的辐射光谱及OH*的二维分布特性进行研究。结果表明：与甲烷火焰的光谱辐射相比,水煤浆火焰不仅存在OH*,CH*和C₂*特征辐射,还产生了Na*,Li*,K*和H*的发射谱线,并出现了连续的黑体辐射,这些光谱辐射特征可作为水煤浆气化或燃烧的标志,也可作为水煤浆是否点燃的判据;通入水煤浆后,OH*强度明显下降,而CH*和C₂*强度增大。对比甲烷火焰OH*二维分布,水煤浆火焰OH*峰值强度明显下降,化学反应区域面积显著减小;沿着火焰传播方向,甲烷和水煤浆火焰轴向的OH*强度均呈先增大后减小的趋势;甲烷火焰径向的OH*在反应核心区出现了双峰形态分布,而水煤浆火焰OH*径向始终呈单峰分布。随着氧碳当量比增大,水煤浆火焰OH*的存在范围扩大,说明氧气的增加促进了OH*的产生;随水煤浆流量提高,OH*的反应核心区域缩小,峰值强度明显下降,CH*,C₂*,Na*,Li*,K*和H*的强度显著增强,连续的黑体辐射强度也明显增大,这些辐射光谱的变化可用于表征操作负荷的变化。     

基于高光谱技术的甲烷层流预混火焰自由基特性研究

高光谱技术提供了空间和光谱维度的信息,同时基于传统黑体模型的实验技术和计算方法不适用于甲烷火焰的辐射特性,而火焰中自由基的高光谱信息反映了火焰结构、组分浓度分布等燃烧的多方面特征,能够为燃烧模型的完善提供依据。利用高光谱技术在不同当量比和不同流量下研究了甲烷预混火焰中自由基的空间和光谱特性。对不同当量比的研究表明,随着当量比的增加,火焰中心处的CH*和C*₂自由基的辐射强度先增加后降低,而燃烧区域内二者的平均辐射强度一直增加,火焰中心处的点可以表征局部的燃烧状态,而燃烧区域内辐射均值表征热释率等整体燃烧状态,定量给出了两种方法的不同趋势。火焰中心处的CH*自由基辐射强度在当量比为1.01时达到峰值,而C*₂自由基辐射强度在当量比为1.12时达到峰值,两种自由基的辐射峰值可以分别作为燃烧中反应强度和稳定性的判据。当量比可以由C*₂和CH*辐射强度之比来表征,修正了C*₂/CH*和当量比的线性关系,提出应使用燃烧区域内C*₂和CH*的平均辐射强度之比,并提出了该比值与当量比的二次关系。利用高光谱技术生成了燃烧区域内C*₂/CH*的云图,得到了详细的空间信息,当量比大于1时,首次在火焰面附近发现了明显的过渡区,体现了高光谱技术的优势。对当量比保持不变情况下的不同流量的研究表明,随着流量的增加,火焰高度增加,而火焰顶部和火焰中心的自由基的浓度分布基本不发生变化,揭示了实验工况下流动的特征时间远小于化学反应特征时间,化学反应过程未受到明显影响。应用高光谱较好的识别出了火焰中的多种自由基,研究了甲烷层流预混火焰中自由基辐射特性及其随着不同当量比和流量变化的趋势,对燃烧现象和机理的认识具有重要意义。     

甲烷-空气预混火焰中OH∗,CH∗发光的定量测量及模拟

基于燃烧化学自发光的诊断技术对发动机诊断、监控有重要意义.针对碳氢燃料燃烧中OH*，CH*激发态物质的生成机理，及其与释热率、当量比的关系进行了实验与模拟探究.首先，利用提出的辐射标定手段对当量比0.7~1.33范围内甲烷-空气预混火焰进行了化学发光量化测量，通过波长分辨的光学收集系统，同时获得各发光组分的浓度，具有很强的便利性.然后采用一维燃烧反应模拟，对与实验工况相同条件下的发光辐射进行定量计算，并对比了释热分布与激发态物质（OH*，CH*，C₂*，CO₂*）的相互关系，计算结果表明，在甲烷-空气层流火焰中，OH*，CH*最合适标识释热率，C₂*次之，CO₂*与释热率分布几乎无相关性.通过实验与计算的对比结果，分析了现有OH*，CH*的各反应通道和常数的准确度，并评估了两自发光组分的主要生成反应路径.      

甲烷/空气同轴射流火焰反应热与OH~*/CH~*定量表征特性研究

激发态粒子是化学反应过程的天然示踪剂,当前激发态粒子研究主要都是在一维简单火焰下进行。为了考察复杂情况下激发态分布对化学反应热的定量表征规律,进行了甲烷/空气同轴射流扩散火焰的试验,开展了甲烷/空气详细燃烧机理与激发态燃烧机理的同轴射流预混火焰和扩散火焰的数值模拟,分析了OH~*和CH~*的分布特性,研究了激发态对反应热的表征关系。结果表明:通过试验ICCD相机和相应滤光片获取的OH~*和CH~*的化学发光图像和数值模拟中OH~*和CH~*摩尔分数分布的模拟结果吻合, OH~*分布主要分为三个燃烧区域, CH~*分布主要为两个燃烧区域。扩散火焰中OH~*和CH~*分布呈现单峰,反应热呈现双峰现象,反应热与激发态变化趋势相似,达到第一峰值后激发态逐渐减小为零,而反应热达到第一峰值后先降低再缓慢上升到第二峰值,最后减小至零。扩散火焰中沿轴向方向,当C_2H+O=CH~*+CO(R12)反应速率达到峰值时,反应热达到第一波峰; H+O+M=OH~*+M(R2)反应速率达到峰值时,反应热达到第二波峰。预混火焰中随着局部当量比的增加, OH~*和CH~*的摩尔分数明显增大,分布区域更广;反应热和激发态粒子OH~*和CH~*的分布趋势一样,在激发态粒子OH~*和CH~*质量分数增大时,反应热也增大,当激发态粒子OH~*和CH~*的质量分数达到峰值的时候,反应热也达到峰值;沿轴向方向, OH~*和CH~*的四个生成反应速率的峰值都在同一位置,反应热也达到最大值, CH+O_2=OH~*+CO(R1)和C_2H+O=CH~*+CO(R12)反应速率相比H+O+M=OH~*+M(R2)和C_2H+O_2=CH~*+CO_2(R11)反应速率更快。     

掺氢甲烷射流扩散火焰燃烧特性分析

采用数值计算与实验相结合的方法研究了掺氢甲烷射流扩散火焰的燃烧特性。结果表明，热量的传递主要是通过热气流对流进行，上游高温气流快速沿轴向流动，径向热量传递较弱；而下游轴向速度降低，热量径向传递增强。喷嘴附近伴流气边界较为稳定，而下游在涡旋作用下出现显著的扰动。射流速度对火焰特性有较大影响，增大射流速度后，火焰高度、辐射强度以及CO、NO、CO₂、H₂O浓度皆显著增加，且辐射强度峰值向下游移动。掺氢量对火焰特性也有重要影响，随着H₂含量增加，燃料向下游传播距离缩短，CO、NO、CO₂浓度降低，H₂O浓度增加。     

N2稀释对双旋流合成气非预混燃烧特性的影响

利用平面激光诱导荧光（PLIF）、高温细丝热电偶及红外气体分析仪对不同N₂稀释量下的双旋流合成气非燃烧流场进行了实验研究。实验结果表明,随着N₂稀释量的增加:双旋流合成气火焰的基本结构发生变化,火焰会出现典型的推举火焰特征;火焰锋面被拉伸,主要反应区域轴向长度和径向宽度增大,喷嘴出口附近火焰锋面由M型分布逐渐转变为W型分布,轴线上未反应区域不断扩大;火焰张角与穿透深度增大,火焰推举高度减小;燃烧排气温度略有减小,NO排放明显降低。     

He/Ar/Kr/F2混合气体中准分子KrF*和Kr2F*的动力学研究

本文用时间积分谱的方法研究了快放电激励下He/Ar/Kr/F₂混合气体中Kr₂F^*的动力学,探讨了Kr₂F^*形成的通道特性,理论分析和实验结果表明,在快放电激励下,三体碰撞过程KrF^*+Kr+M→Kr₂F^*+M是形成Kr₂F^*的主要通道,而置换反应ArF^*+Kr→KrF^*+Ar又能有效地产生KrF^*。实验测量了Kr₂F^*形成的有关速率常数。 关键词：     

一维湍流模型对He平面羽流和CH4/H2/N2射流火焰的数值模拟

本文采用一维湍流模型（ODT）对氦气平面羽流和CH₄/H₂/N₂射流火焰进行数值模拟,和前人的实验结果进行定量地对比。结果表明,ODT模型能够准确地预测平面羽流基本特征,湍流涡团的分布同流场拉伸率之间具有密切的关系,涡团强度的分布能够直观地表明当地的湍流强度。ODT模型埘CH₄/H₂/N₂瞬态火焰的模拟定性反应了火焰特性及其与湍流作用的规律,对温度-混合分数的预测值和实测值进行比较,发现甲烷火焰燃烧在富燃料侧并未达到平衡状态,因而基元反应对火焰特性的预测具有重要作用。     

基于欧拉喷雾与火焰面/反应进度变量湍流燃烧模型的PRF燃料喷雾燃烧模拟研究

在开源计算流体软件OpenFOAM环境下,将基于欧拉方法的Σ-Y喷雾模型与非稳态火焰面/反应进度变量湍流燃烧模型相耦合,发展用于高温高压环境中液体燃料喷雾湍流燃烧高精度计算模型,分别对非燃烧和燃烧工况下的五种典型参比燃料(Primary Reference Fuel,PRF)的燃油喷射雾化与湍流燃烧过程开展数值研究。结果表明：该新型耦合模型能够准确的预测PRF燃料的喷雾和着火燃烧特性;所开发的重构数值喷雾纹影图像和燃烧OH^*图像处理方法能够很好地捕捉到试验的滞燃期、火焰浮起长度及喷雾火焰结构;研究揭示了不同比例的PRF燃料对喷雾及着火燃烧过程的影响特性,为替代燃料在发动机上的高效应用提供了理论指导。     

Scaling relations for the length of coaxial oxy-flames with and without swirl

An extensive experimental study is carried out to analyze scaling laws for the length of methane oxy-flames stabilized on a coaxial injector. The central methane fuel stream is diluted with N₂, CO₂ or He. The annular air stream is enriched with oxygen and can be impregnated with swirl. Former studies have shown that the stoichiometric mixing length of relatively short flames is controlled by the mixing process taking place in the vicinity of the injector outlet. This property has been used to derive scaling laws at large values of the stoichiometric mixture fraction. It is shown here that the same relation can be extended to methane oxy-flames characterized by small values of the stoichiometric mixture fraction. Flame lengths are determined with OH* chemiluminescence measurements over more than 1000 combinations of momentum ratio, annular swirl level and composition of the inner and outer streams of the coaxial injector. It is found that the lengths of all the flames investigated without swirl collapse on a single line, whose coefficients correspond to within 15% of flame lengths obtained for fuel and oxidizer streams at much larger stoichiometric mixture fractions. This relation is then extended to the case of swirling flames by including the contribution of the tangential velocity in the flow entrainment rate and is found to well reproduce the mixing degree of the two co-axial streams as long as the flow does not exhibit a vortex breakdown bubble. At higher swirl levels, when the flow features a central recirculation region, the flame length is found to also directly depend on the oxygen enrichment in the oxidizer stream.     

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

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

Decreasing liftoff height behavior in diluted laminar lifted methane jet flames

Stabilization of laminar lifted coflow jet flames of nitrogen-diluted methane was investigated experimentally and numerically. As the fuel jet velocity was increased, two distinct behaviors in liftoff height were observed depending on the initial fuel mole fraction; a monotonically increasing trend and a decreasing and then increasing trend (U-shaped behavior). The former was observed in the jet-developing region and the latter in the jet-developed region. Because the decreasing behavior of liftoff height with jet velocity has not been observed at ambient temperature, the present study focuses on decreasing liftoff height behavior. To elucidate the physical mechanism underlying the U-shaped behavior, numerical simulations of reacting jets were conducted by adopting a skeletal mechanism. The U-shaped behavior was related to the buoyancy. At small jet velocities, the relative importance of the buoyancy over convection was strong and the flow field was accelerated in the downstream region to stabilize the lifted flame. As the jet velocity increased, the relative importance of buoyancy decreased and the liftoff height decreased. As the jet velocity further increased, the flame stabilization was controlled by jet momentum and the liftoff height increased.     

Stabilized flames in hybrid aluminum-methane-air mixtures

A premixed methane–air bunsen-type flame is seeded with micron-sized (d₃₂ = 5.6 μm) atomized aluminum powder over a wide range of solid fuel concentrations. The burning velocities of the resulting two-phase hybrid flame are determined using the total surface area of the inner flame cone and the known volumetric flow rate, and spatially resolved flame spectra are obtained with a spectral scanning system. Flame temperatures are derived through polychromatic fitting of Planck’s law to the continuous part of the spectrum. It is found that an increase in the solid fuel concentration changes the aluminum combustion regime from low temperature oxidation to full-fledged flame front propagation. For stoichiometric methane–air mixtures, the transition occurs in the aluminum concentration range of 140–220 g/m³ and is manifested by the appearance of AlO sub-oxide bands and an increase in the flame temperature to 2500 K. The flame burning velocity is found to decrease only slightly with an increase in aluminum concentration, in contrast to the rapid decrease in flame speed, followed by quenching, that is observed for flames seeded with inert SiC particles. The observed behavior of the burning velocity and flame temperature leads to the conclusion that intense aluminum combustion in a hybrid flame only occurs when the flame front propagating through the aluminum suspension is coupled to the methane–air flame.     

Strain rate and fuel composition dependence of chemiluminescent species profiles in non-premixed counterflow flames: comparison with model results

A detailed comparison has been conducted between chemiluminescence (CL) species profiles of OH^?, CH^?, and C₂ ^?, obtained experimentally and from detailed flame kinetics modeling, respectively, of atmospheric pressure non-premixed flames formed in the forward stagnation region of a fuel flow ejected from a porous cylinder and an air counterflow. Both pure methane and mixtures of methane with hydrogen (between 10 and 30 % by volume) were used as fuels. By varying the air-flow velocities methane flames were operated at strain rates between 100 and 350 s^?1, while for methane/hydrogen flames the strain rate was fixed at 200 s^?1. Spatial profiles perpendicular to the flame front were extracted from spectrograms recorded with a spectrometer/CCD camera system and evaluating each spectral band individually. Flame kinetics modeling was accomplished with an in-house chemical mechanism including C₁–C₄ chemistry, as well as elementary steps for the formation, removal, and electronic quenching of all measured active species. In the CH₄/air flames, experiments and model results agree with respect to trends in profile peak intensity and position. For the CH₄/H₂/air flames, with increasing H₂ content in the fuel the experimental CL peak intensities decrease slightly and their peak positions shift towards the fuel side, while for the model the drop in mole fraction is much stronger and the peak positions move closer to the fuel side. For both fuel compositions the modeled profiles peak closer to the fuel side than in the experiments. The discrepancies can only partly be attributed to the limited attainable spatial resolution but may also necessitate revised reaction mechanisms for predicting CL species in this type of flame.     

Dynamic behavior of a freely-propagating turbulent premixed flame under global stretch-rate oscillations

Dynamic features of a freely propagating turbulent premixed flame under global stretch rate oscillations were investigated by utilizing a jet-type low-swirl burner equipped with a high-speed valve on the swirl jet line. The bulk flow velocity, equivalence ratio and the nominal mean swirl number were 5 m/s, 0.80 and 1.23, respectively. Seven velocity forcing amplitudes, from 0.09 to 0.55, were examined with a single forcing frequency of 50 Hz. Three kinds of optical measurements, OH-PLIF, OH^* chemiluminescence and PIV, were conducted. All the data were measured or post-processed in a phase-locked manner to obtain phase-resolved information. The global transverse stretch rate showed in-phase oscillations centering around 60 (1/s). The oscillation amplitude of the stretch rate grew with the increment of the forcing amplitude. The turbulent flame structure in the core flow region varied largely in axial direction in response to the flowfield oscillations. The flame brush thickness and the flame surface area oscillated with a phase shift to the stretch rate oscillations. These two properties showed a maximum and minimum values in the increasing and decreasing stretch periods, respectively, for all the forcing amplitudes. Despite large variations in flame brush thickness at different phase angles, the normalized profiles collapse onto a consistent curve. This suggests that the self-similarity sustains in this dynamic flame. The global OH^* fluctuation response (i.e. response of global heat-release rate fluctuation) showed a linear dependency to the forcing velocity oscillation amplitudes. The flame surface area fluctuation response showed a linear tendency as well with a slope similar to that of the global OH^* fluctuation. This indicated that the flame surface area variations play a critical role in the global flame response.     

Experimental study of air dilution in oxy-liquid fuel flames

The influence of oxidizer dilution in oxy-liquid ethanol flames is experimentally investigated by using a coaxial air-assisted injector positioned in a vertical combustion chamber. This study accounts for the influence of a two-phase mode since two different injector geometries are used: for the first configuration, a vaporization mode is observed at nominal power in oxy conditions, while for the second one, a brush mode is observed. Dilution with air is applied by keeping oxidizer velocity constant. Flame structure is observed through CH emission: dilution leads to an increase in the flame diameter, and collective effects of two-phase combustion are encouraged. The effect of dilution on oxy flame stability is also studied: for a given oxygen mass fraction in the oxidizer, the oxidizer flow rate is increased until extinction occurs. Dilution leads to a less stable flame, which may be essentially explained by the decrease in laminar flame speed with dilution. For high oxidizer dilution levels, the change in flame structure might be another parameter to consider. Finally, species concentrations are measured using a standard gas sampling technique. NO and CO evolutions with dilution are different between both two-phase combustion regimes. An empirical approach based on thermal NO mechanism and CO oxidation reaction enables one to explain the evolutions for brush mode. For vaporization mode, the residence time in burned gases is also to be considered.     

氮气稀释下丙烷层流扩散火焰的抬升特性研究

本文实验研究了氮气稀释条件下,丙烷层流扩散火焰的抬升特性。主要研究了混合物组分及射流出口速度对丙烷抬升火焰抬升高度的影响,得到了氮气稀释下丙烷抬升火焰的再附着滞后现象,并采用计算方法得出了理论火焰抬升高度,与实验数据进行了对比分析,二者符合很好。     

Experimental comparison of opposed and concurrent flame spread in a forced convective microgravity environment

Flame spread experiments in both concurrent and opposed flow have been carried out in a 5.18-s drop tower with a thin cellulose fuel. Flame spread rate and flame length have been measured over a range of 0–30 cm/s forced flow (in both directions), 3.6–14.7 psia, and oxygen mole fractions 0.24–0.85 in nitrogen. Results are presented for each of the three variables independently to elucidate their individual effects, with special emphasis on pressure/oxygen combinations that result in earth-equivalent oxygen partial pressures (normoxic conditions). Correlations using all three variables combined into a single parameter to predict flame spread rate are presented. The correlations are used to demonstrate that opposed flow flames in typical spacecraft ventilation flows (5–20 cm/s) spread faster than concurrent flow flames under otherwise similar conditions (pressure, oxygen concentration) in nearly all spacecraft atmospheres. This indicates that in the event of an actual fire aboard a spacecraft, the fire is likely to grow most quickly in the opposed mode as the upstream flame spreads faster and the downstream flame is inhibited by the vitiated atmosphere produced by the upstream flame. Additionally, an interesting phenomenon was observed at intermediate values of concurrent forced flow velocity where flow/flame interactions produced a recirculation downstream of the flame, which allowed an opposed flow leading edge to form there.     

Stabilization, propagation and instability of tribrachial triple flames

A tribrachial (or triple) flame is one kind of edge flame that can be encountered in nonpremixed mixing layers, consisting of a lean and a rich premixed flame wing together with a trailing diffusion flame all extending from a single point. The flame could play an important role on the characteristics of various flame behaviors including lifted flames in jets, flame propagation in two-dimensional mixing layers, and autoignition fronts. The structure of tribrachial flame suggests that the edge is located along the stoichiometric contour in a mixing layer due to the coexistence of all three different types of flames. Since the edge has a premixed nature, it has unique propagation characteristics. In this review, the propagation speed of tribrachial flames will be discussed for flames propagating in mixing layers, including the effects of concentration gradient, velocity gradient, and burnt gas expansion. Based on the tribrachial edge structure observed experimentally in laminar lifted flames in jets, the flame stabilization characteristics including liftoff height, reattachment, and blowout behaviors and their buoyancy-induced instability will be explained. Various effects on liftoff heights in both free and coflow jets including jet velocity, the Schmidt number of fuel, nozzle diameter, partial premixing of air to fuel, and inert dilution to fuel are discussed. Implications of edge flames in the modeling of turbulent nonpremixed flames and the stabilization of turbulent lifted flames in jets are covered.