合成气预混层流火焰结构的实验和数值研究

第48卷第7期西安交通大学学报VoL 48 No. 72014年7月JOURNAL OF XI'AN JIAOTONG UNIVERSITYJul.2014DoI:10.7652/ xjtuxb201407007合成气预混层流火焰结构的实验和数值研究卫之龙,王金华,舒新建,谢永亮,王锡斌,黄佐华(西安交通大学动力工程多相流国家重点实验室,710049,西安)摘要:利用 OH-PLIF方法获得了当量比分别为0.6、0.8、1.0、1.2,CO2或N2稀释比分别为3%5%时,合成气/空气/稀释气本生灯预混层流火焰中OH基的分布,结合 STAR-CD模拟计算所得火焰中的流场和组分分布进一步分析了火焰结构。研究结果表明:随着混合气当量比的增加,OH基髙浓度分布区域由火焰前锋面附近转移到火焰边缘;混合气较稀时,火焰前锋面附近OH基浓度最高且沿已燃区方向逐步递减,火焰顶端处OH基浓度减小,模拟计算结果显示火焰顶端并未发生燃料泄漏;化学当量比下,火焰前锋面附近和火焰边缘区域OH基浓度较高,火焰前锋面附近出现了预混燃烧区和扩散燃烧区,该区域中OH基呈现“W”型分布;受N2和CO2稀释的影响,混合气层流燃烧速度降低,火焰前锋面拉长,CO2对火焰结构的影响比N2更显著;火焰前锋面附近OH基浓度减小,扩散燃烧区OH基浓度增大,说明火焰的预混燃烧有所减弱,扩散燃烧有所加强。关键词:合成气;OH-PLIF方法; STAR-CD模拟;OH基;稀释气中图分类号:TK411文献标志码:A文章编号:0253-987X(2014)07-003407Experimental and Numerical Study on Structure ofLaminar Premixed Syngas-Air-Diluents FlamesWEI Zhilong, WANG Jinhua, SHU Xinjian, XIE Yongliang, WANG Xibin, HUANG Zuohua(State Key Laboratory of Multiphase Flow in Power Engineering, Xi'an Jiaotong University, Xi'an 710049, China)Abstract: The structures of laminar premixed syngas-air-diluents flames were obtaineddifferent equivalence ratios(0.6,0.8,1.0,1. 2), diluents(N2/cO2), and dilution ratios(3%and 5%)by using OH-PLIF technique. Flame structure was analyzed on the basis of the flowfield and species distribution which were calculated with STAR-CD. The results show that highoh distribution shifts towards the outer flame boundary with the increase of equivalence ratio;for lean mixture, the region of high OH radical concentration locates along the conic primaryreaction boundary OH radical concentration at the flame cone tip decreases, and the numericasimulation results show that no leakage occurs at the flame tip: for mixture at stoichiometricratio, high OH radical concentration appears along the flame front and on the flame surfacemeaning the occurrence of a premixed combustion zone and a diffusion combustion zone, and theoh distribution appears to resemble a "W"shape; due to the effects of diluents n2 and CO2, thelaminar burning velocity decreases, leading to the extension of the flame front, and COz has agreater impact on the flame structure compared to n2: Oh radical concentration decreases alongthe flame front while increases in the diffusion zone, which means the premixed combustion isweakened while the diffusion combustion is intensified收稿日期:2013-1219.作者简介:卫之龙(1990-),男,硕生生;王金华(通H中国煤化工科学基金资助项目(51376004,51006080)。CNMHG项目,国家自然网络出版时间:201404-16网络出版地址:http∥www.cnkinet/kcms/detail/61.1069.t.20140416.1746.006.html第7期卫之龙,等合成气预混层流火焰结构的实验和数值研究Keywords: syngas; OH-PLIF; STAR-CD; OH distribution; diluent随着石油资源的日益减少和排放法规的日趋严预混层流本生灯火焰中间产物OH基的分布,并结格,寻找切实可行的石油替代燃料具有重要意义。合 STAR-CD模拟计算结果,分析研究了预混层流甲醇由于来源丰富、可规模化生产成本较低廉,近年火焰结构特性和火焰自身不稳定特性。考虑到发动来便成为替代燃料研究中的一个热点。甲醇可以直机通常采用废气再循环(EGR)来降低燃烧温度和接用于发动机,但是热值较汽油低,在发动机运行过抑制NO2的生成,本文分别采用CO2和N2作为稀程中消耗量较大1。因此,提高甲醇在发动机上的释气,测量和分析了不同当量比和不同稀释比条件热效率、改善发动机性能,成为甲醇燃料发动机研究下合成气火焰的结构特性。的重点内容。利用发动机废气余热可以将甲醇裂解,从而得到合成气(H2和CO的体积比为2:1),1实验装置和实验过程将其作为发动机燃料,可以提高燃料能量品位和发图1为实验装置示意图,由 OH-PLIF系统、供气动机效率2。系统和本生灯组成。 OH-PLIF系统包括Nd:YAG目前,有关甲醇裂解所得合成气的研究多集中激光器(型号为 Quanta-Ray Pro190,频率为10Hz在发动机实验方面35。H2特殊的强反应性和输运脉冲时间为10ns,脉冲能量为300mJ)染料激光特性,使得大比例H2的合成气在火焰面结构和火器( Sirah Prsc-c300)、能量检测器、BBO倍频焰自身不稳定性上表现出与传统碳氢燃料的显著差器、透镜组、带有光线增强器( Lavision vc08-0094)异,进而影响到燃烧室内流动与火焰相互作用,最终的CCD相机( La vision Image Prox)及加装OH带影响发动机缸内燃烧过程。本生灯作为一种常见的通滤波器( La vision vz08-0222的紫外线滤光镜头燃烧器,能够提供稳定的层流火焰结构,是研究火焰( Nikon Rayfact PF10545MFUV)。供气系统由气结构和火焰自身不稳定性的有效手段。Jin等利用瓶、预混室(容积为400mL)和流量计(MKs流量计,平面激光诱导荧光研究了不同掺氢比及不同当量比型号分别为1179A和1559A,精度为1/1000L/min)对合成气本生灯预混层流火焰结构的影响。Law组成。为确保层流燃烧,混合气的Re<2000。本生等利用相机图片研究了不同碳氢燃料的预混层流本灯出口直径为5mm,出口长度为30cm(大于50倍生灯火焰,分析了燃料优先扩散和火焰拉伸的耦合管径),底部装有防止回火的单向阀。作用对火焰顶端结构的影响8。 Kozlovsky等利用数值模拟计算研究了路易斯数对本生灯层流火焰0=R结构的影响9。 Bouvet等利用化学发光方法研究Dye激光冷却水出口了CO/H2在不同比例下的预混层流本生灯火焰的YAG激光冷却水进口灯混合器顶端开口现象[10。由上述研究可以看出,预混层流火焰结构的研究不仅可以帮助认识火焰传播、局部数招驾菜器熄火、着火极限等基本的火焰现象,而且有利于理解优先扩散、路易斯数、火焰拉伸等因素对火焰不稳定性的影响,进而解释湍流燃烧等复杂燃烧现象。此图1实验装置示意图外,作为燃料燃烧的一种常见的中间产物,OH基广合成气由纯度(质量分数)分别为9.99%的泛存在于反应区和已燃区,通常被用来界定燃烧反H2和9.9的CO严格按照流量法以2:1的质量应区1,OH基分布在一定程度上能够反映火焰的比配置而成。稀释比φ定义为稀释气在混合气中局部燃烧状况,在燃烧过程中大量存在且易于检测所占的体积分数,即因此OH基通常作为燃烧诊断以及燃烧分析的标定物。 Makmool等利用相机图片、化学发光方法和Dilution +Vair +v平面激光诱导荧光研究了液化石油气预混火焰的火实验过程中,YAG激光器产生355mm的激光焰结构,表明平面激光诱导荧光是研究火焰燃烧特性经染料激光器中国煤化工769mm,该的最可靠、准确的方法12波长可用于激CNMHG(1,0)跃迁本文利用平面激光诱导荧光( OH-PLIF)测量的Q1(8)线,每个脉冲激光能量约为8mJ。脉冲激http∥www.jdxb.cnhttp∥zkxb.xjtu.edu.cn西安交通大学学报第48卷光在通过片光镜头组后产生高约50mm的片光,片光通过本生灯出口的中心轴线激发火焰中的OH生灯喷管基产生波长为308nm左右的荧光。在与片光垂直的方向上安装了ICCD相机,通过调整拍摄延迟使相机与激光同步,从而保证ICCD相机可以完整捕压力边界捉到激发的荧光信号。焰反应区混合气入口本文研究了不同当量比(中=0.6,0.8,11.2),CO2或N2稀释比中分别为3%、5%时的合(a)整体图成气预混层流火焰结构。每个工况下,ICCD相机连续拍摄100张OH基分布照片,并将100张照片的平均结果用于OH基分布分析。2 STAR-CD模拟计算利用 OH-PLIF获得的火焰结构中OH基的分布是实验中重要的信息,火焰中的OH基分布在反应区和已燃区,单一的OH基分布信息不足以解释火焰结构。CFD模拟计算可以提供较多与火焰结构相关的信息,如速度场温度场组分场等,利用这些丰富的数据信息、结合OH-PLIF实验图像,有助(b)反应区局部放大图于分析火焰中某些现象的成因以及各因素对于火焰图2 STAR-CD模型网格的整体图和反应区局部放大图结构的影响和作用本文采用 SATR-CD软件来模拟某些工况下合成气的预混层流火焰,本生灯模型按11比例建立。由于本生灯火焰是轴对称的三维火焰,所以可将三维火焰计算简化为二维火焰计算,计算网格为单层的六面体网格。为获得较好的模拟结果,对火焰反应区网格进行加密,最小网格尺寸为0.2mm,网格总数为16800。图2为模型网格的整体图和反应区局部放大图。模型中燃烧化学反应的机理采用Scott13的CO和H2简化机理(12种组分,38步基元反应)。计算时间为0.07s,可保证计算火焰达到图3Re=1800时合成气火焰数码照片稳定状态时间步长为106s。初始坏境条件为P0.1MPa,T。=298K。着当量比的增加,火焰前锋面高度降低。这是由于随着当量比的增加,混合气的层流燃烧加速,而混合3结果和分析气来流速度基本不变,使得火焰前锋面向着混合气图3为Re=1800时合成气火焰的数码照片。来流方向移动,由此导致了火焰前锋面高度降低,面由图3可以看出随着当量比的增大,火焰前锋面发积减小。图4为火焰前锋面高度随当量比的变化,光强度增强且发光的颜色随之变化。这是由于随着其也表现出了图3的变化规律。当量比的增加,混合气中燃料比例的变化使得火焰图5为Re=1800时合成气火焰的OH-PLIF前锋面处的反应更加激烈,锋面厚度增加,发光强度照片。图中火焰区域内颜色较深的部分代表OH增大。φ=0.6时,火焰前锋面发出的蓝光是由火焰基浓度高,中心轴附近OH浓度梯度最大处是火焰前锋面上激发态的CH活性基发出的;=1.2时,前锋面中国煤化工当量比的增加,混合气较浓,火焰前锋面发出的光略带黄光,这是含OH基高沐CNMH〔峰面附近向火焰碳原子团受激发发出的1。由图3还可以看出,随边缘转移,呈现出了不同的分布形态。这些变化在http∥www.jdxb.cnhttp∥zkxb.xjtu.edu.cn第7期卫之龙,等:合成气预混层流火焰结构的实验和数值研究w(OH)0.4841×1020.4150×100.2767×10-2102075×102(a)OH的质量分数图4火焰锋面高度随当量比的变化温度kK1893速度ms7.34916654.2053.157(b)火焰前锋面速度和温度图5Re=1800时合成气火焰 OH-PLIF照片w(CO)w(H2)0.7463×1011066×101普通数码照片中是观察不到的,反映了 OH-PLIF0.6397×100.9137×100.5331×1010.7614×102技术表征火焰锋面结构的优势。在混合气较稀的条件下,OH基主要分布于火焰前锋面附近;在混合气0.4265×100.6091X10较浓的条件下,OH基不仅分布于火焰前锋面附近,0.3198×1004569×10-2而且在火焰边缘也有较高浓度的分布,且随着当量0.2132×10.3046×10比增加浓度有所提高,这说明此时火焰形成了预混0.1066×10-10.523×10-2燃烧区和扩散燃烧区。此外,混合气较稀时,OH基0.6599×100.6303×10的分布出现了顶端开口,而在化学当量比条件下,(c)CO、H2质量分数图6∮=0.6时合成气火焰 OH-PLIF照片及STAROH基出现了“W”型的分布。CD模拟结果顶端开口是层流本生灯火焰的一个经典现象,该现象的研究对于更好地理解火焰燃烧特性具有重响:①本生灯层流火焰在火焰前锋面顶端处拉伸率要意义。本实验中,φ=0.6条件下数码照片并未表达到最大,导致此处混合气中优先扩散作用加强,使现出明显的顶端开口,但是在OH-PIF图片中OH得在火焰前锋面顶端的燃烧强度降低;②火焰前锋基分布出现了顶端开口,这说明此时火焰虽然没有面速度随着火焰前锋面高度的升高而增大,且在锋出现顶端开口,但是顶端燃烧强度已经有所降低,所面顶端达到最大值,速度方向是由预混燃烧区指向以OH基分布出现了顶端开口。图6为φ=0.6时扩散燃烧区,所以锋面顶端的OH基受到的输运影合成气火焰OH-PILF照片及 STAR-CD模拟结果。响较强。火焰温度可以直接反映燃烧的剧烈程由图6a可以看出,OH基表现出了相同的分布趋度,因此可以根据已燃区温度分布来判断燃烧强势,并且顶端开口存在,表明STAR-CD模拟计算结度。由图6b可以看出,火焰前锋面顶端的最高温度果可较为准确、定性地反映实际燃烧情况。由于较低这说明Ⅵ中国煤化工度减小。综OH基的分布是化学反应与输运共同作用的结果,合以上两方面CNMHG大的拉伸率,因此解释OH基在顶端开口需要考虑两方面的影由此加强了优先扩散的影响,导致该处燃烧强度降http∥www.jdxb.cnhttp∥zkxb.xjtu.edu.cn38西安交通大学学报第48卷低,OH基的生成量较其他反应区域有所减少,这是反应率。由图8可以看出,在化学当量比条件下,造成OH基顶端开口的主要原因;受速度场输运的H2总反应率及其峰值低于CO,说明在该条件下影响,OH基的浓度减小,这进一步加剧了顶端开相较于CO,H2在锋面附近可以更快地与混合气中口。由图6c中CO及H2的分布还可以看出,顶端的O2迅速反应且完全消耗掉,而CO与混合气中开口并未伴随锋面顶端处的燃料泄漏,这与之前的O2反应较慢,不能在锋面附近快速、完全消耗掉,从研究结果是一致的)。而向外扩散且与环境中的O2继续进行反应,形成在φ=0.8的条件下,顶端开口消失。这是由于扩散燃烧。因此,火焰反应区分为预混燃烧区和扩随着当量比的增加,火焰前锋面高度降低,火焰顶端散燃烧区,OH基呈现出“W”型分布。前锋面的曲率半径增大,使得锋面顶端附近拉伸率减小,由此削弱了顶端结构对质量扩散的促进作用混合气中燃料比例增大,导致火焰顶端的优先扩散日-0.02作用削弱。因此,顶端开口在该情况下消失。在化学当量比条件下,OH基出现了“W”型分布,了解“W”型分布的成因可以更好地理解燃烧微观过程,实现调整参数对火焰的控制。图7为合成气火焰OH-PLIF照片与 STAR-CD模拟计算的火焰传播方向CO和H2分布。由图7可以看出:在=0.8条件图84=.0时合成气火焰中H与CO的总反应率下,CO和H2可以在火焰前锋面附近基本消耗掉,图9为Re=1800、中=1.0时不同N2稀释比下因此该条件下OH基在火焰前锋面附近浓度较高;在∮=1.0条件下,火焰前锋面附近H2可以完全消消的合成气火焰 OH-PLIF照片。由图9可以看出,随着N2稀释比的增加,火焰前锋面拉长,前锋面附耗掉,但是火焰前锋面附近CO未能完全消耗,出现了明显的扩散现象。利用 CHEMKIN2.0和GRI近及顶端聚集的OH基浓度降低,火焰边缘的OHmech3.0机理可以对H2和CO总反应率进行基浓度明显增加。这是由于混合气中的N2本身不计算,图8为d=1.0时合成气火焰中H2和CO总参与反应,但吸收热量,而N2稀释比增加会导致反w(Co)应物的放热量减少,燃烧产物的吸热量增多,从而使0.923×100.130×10-1火焰温度降低,H2及CO反应速率随之降低。因0.791×100.111×10此,火焰前锋面附近的反应强度削弱,火焰层流燃烧0.659×100.929×1速度降低,导致火焰前锋面拉长。火焰温度降低,在0.527×100sy×102火焰前锋面附近H和CO因反应强度降低而不能完0.264×100.371×10-2全消耗掉,且扩散到火焰边缘,导致扩散燃烧加剧。0.132×100.293×100946×10-9(a)中=0.8w(CO)0.1180.68×10-10.1010.144×100.840×10-10.20×10-10.672×10-10.960×1020.504×100.336×100.480×10-20.168×100.240×10-20.346X100.238×10-7图9Re=18000时不同N2稀释比下的合成(b)=1.0气火凵中国煤化工图7合成气火焰 OH-PLIF照片与 STAR-CD模拟CNMHG计算的CO和H2分布图10为ke15、p-1,0、,=5%时不同稀http∥www.jdxb.cnhttp∥zkxb.xjtu.edu.cn卫之龙,等:合成气预混层流火焰结构的实验和数值研究释气下的火焰 OH-PLIF照片。由图10可以看出,稀释气本生灯预混层流火焰OH基分布,利用混合气在掺混CO2之后,火焰前锋面进一步拉长, STAR-CD对火焰中的流场和组分分布进行了模拟火焰前锋面附近的OH基浓度进一步降低,火焰边计算,结合 OH-PLIF实验图像研究了火焰前锋面缘OH基浓度有所增加。这是由于掺混CO2对火结构,得到的主要结论如下。焰结构的影响与掺混N2的影响基本相同,CO2的(1)混合气在较稀条件下,OH基主要分布在火吸热能力比N2更强,因此对火焰温度及反应速率焰前锋面附近且向已燃区方向发展,浓度呈现逐步的影响更大。CO2掺混比例的增加会抑制CO的氧递减的趋势。因火焰前锋面顶端拉伸率最大,此处化反应,导致预混燃烧减弱,扩散燃烧增强。混合气中优先扩散作用加强,火焰前锋面顶端的燃烧强度降低,OH基的生成量减少,加之锋面顶端输运的影响,OH基顶端出现开口。模拟计算结果显示,此时并未伴随顶端燃料泄漏。(2)化学当量比条件下,混合气OH基在火焰前锋面附近,火焰边缘区域浓度较高。H2的净消耗率远大于CO,H2在火焰前锋面附近很快完全消耗掉,而CO不能在火焰前锋面被完全消耗而发生扩散,所以火焰出现预混燃烧区和扩散燃烧区,OH基呈现“W”型分布。(3)掺混N2和CO2会导致火焰前锋面拉长,同未稀释N2稀样CO2稀释时OH基浓度在火焰锋面附近减小,在扩散燃烧区图10Re=1800、中=1.0、中=5%时不同稀释气下的火焰 OH-PLIF照片增大,说明火焰的预混燃烧有所减弱,扩散燃烧有所加剧。在相同稀释比下,CO2由于具有更强烈的吸图11为φ=1.0时火焰前锋面高度及层流火焰热能力和对CO氧化反应的抑制作用所以对火焰速度随φ的变化。由层流本生灯计算层流火焰速度结构的影响比N2更显著。公式SL=sin(6/2)可知,出口气体速度保持不变时,层流火焰速度越小,θ越小,从而火焰前锋面越高。参考文献:CO2,高度[1] AOLAH G, GOEPPERT A, PRAKASH SG K→CO2,速度中1N2,速度Weinheim, Germany Wiley, 2006: 179-230140[2]冯春晃,高孝洪.甲醇裂解燃料点燃式发动机的研究].内燃机学报,1989,7(2):151-157128绥FENG Chunhuang, GAo Xiaohong. Researches on124spark ignition engine fueled with dissociated methanol[J]. Transactions of CSICE, 1989, 7(2): 151-157.中116112[3]姚春德,徐元利,杨建军,等.甲醇裂解气对点燃式中%发动机性能影响研究[].工程热物理学报,2009,30图11φ=1.0时火焰前锋面高度及层流火焰速度随(2):353-356.中的变化YAO Chunde, XU Yuanli, YANG Jianjun, et al.Effects of dissociated methanol on perfo由图11可以看出,由于掺混CO2和N2,层流[J]. Je火焰速度下降,火焰前锋面高度上升。在相同稀释比下,相比于掺混N2,掺混CO2对层流火焰速度的[4]杨成,王炳辉,励金祥,等,柴油机进气预混微量甲影响更显著,火焰前锋面高度增加更显著。醇裂解气的试验研究[J].车用发动机,2009,3:774结论YANG C中国煤化工 Jinxian,era.E:CNMH Geled with diesel利用 OH-PLIF激光测量技术获得了合成气/blended with dissociated methanol gas [J].Vehiclehttp∥www.jdxb.cnhttp∥zkxb.xjtu.edu.cr西安交通大学学报第48卷Engine,2009,3:77-81density equation in nonpremixed turbulent combustion[5]方葛文,高孝洪.点燃式发动机掺烧甲醇裂解气的研[]. Proceedings of the Combustion Institute, 1996究[冂].内燃机学报,1993,11(2):131-136.26(1):35-42.FANG Gewen, GAO Xiaohong. Study on spark- [16] SMITH G P, GOLDEN D M, FRENKLACH M, etgnition engine fueled with gasoline blended with dissoal.GriMech-3.0[eb/Ol].[2013-11-20].http:/ciated methanol gas [J]. Transactions of CSICEwww.mw.berkeleyedu/grimech,1993,11(2):131-136.L6FUJ, TANG C L,JNW,etal. Study on laminar[本刊相关文献链接]lame speed and flame structure of syngas with varied陈朝阳,张春化,陈昊,等,当量比对二甲醚氧化分解过程的compositions using OH-PLIF and spectrograph [JJ. 8. 2014, 48(5):32-36+77. [doi: 10. 7652/xjtuxb201405International Journal of Hydrogen Energy, 2013, 38 006(3):1636-1643.魏若男,周蓉芳,周竹杰,等.双电极点火时间间隔对天然气[7] MIZOMOTO M, ASAKA Y, IKAI S,etal. Effects空气预混合气层流燃烧的影响,2014,48(1):7-12.[doi:of preferential diffusion on the burning intensity of DOI: 10. 7652/xjtuxb201401002]curved flames[门. Proceedings of the Combustion In-康婵,杨星,刘杰,等.负电场下点电极和网状电极对预混稀stitute,1985,20(1):1933-1939燃火焰的影响.2014,48(1):31-36.[doi:10.7652/ xjtuxb[8] LAW CK, ISHIZUKA S, CHO P. On the opening of 201401006]premixed Bunsen flame tips[门. Combustion Science张猛,王金华,谢永亮,等.CH4/H2/air预混潴流火焰前锋面and Technology, 1982, 28(3/4):89-96结构探测,2013,47(11):14-18.[doi:10.7652/ xjtuxb2013[9 KOZLOVSKY G, SIVASHINSKY G I. On open and11003]losed tips of Bunsen burner flames[J. Theoretical王金华,张猛,谢永亮,等,合成气纯氧高压预混湍流火焰结and Computational Fluid Dynamics,1994, 6(2 /3): #%f *. 2013, 47(10): 87-90. doi: 10. 7652/xjtuxb201310l81-192015[10] BOUVET N, CHAUVEAU C, GOKALP I,etal.孟禅文,杨星,廉婵,等.直流电场对预混CH,/O2/N2火焰Characterization of syngas laminar flames using the传播特性影响的试验研究.2013,47(7):13-17.[doi:10.Bunsen burner configuration [J]. International Journal 7652/ xjtuxb201307003]of Hydrogen Energy, 2011, 36(1):992-1005张进华,李婷,王孙安,等。可变视场下的火灾探测算法[I1] QUANG V N, PHILLIP H P. The time evolution of a 2012, 46(10): 29-35. [doi: 10. 7652/ xtuxb201210006]vortex-flame interaction observed via planar imaging of唐安东,孟祥文,周蓉芳,等.非均匀电场对火焰传播速率的CH and OH [J]. Proceedings of the Combustion Insti- t# 2012, 46(9): 16-20. [ doi: 10. 7652/ xjtuxb201209004]tute,l996,26(1):357-364周蓉芳,魏若男,周竹杰,等,点火能量对天然气空气预混合12] MAKMOOL U, JUGJAI S, TIA S. Structures and气层流燃烧的影响.2012,46(7):21-25.[doi:10.7652/performances of laminar impinging multiple premixed xjtuxb201207005LPG-air flames []. Fuel, 2013, 112(10):254-262刘德新,刘斌.二次喷油优化直噴汽油机冷启动排放特性的[13] DAVIS S G, JOSHI A V, WANG H, et al. An opti- #* 2012, 46(1): 13-18. [doi: 10.7652/xjtuxb201201003]mized kinetic model of h2/ CO combustion[冂].Pro刘亦夫,刘兵,刘亮,等.天然气缸内直喷发动机在不同喷射s of the Combustion Institute,2005,30(1):和点火时刻下的排放与燃烧特性,2011,45(5):1216.[doi1283-1292.10.7652/ xjtuxb20105003][14]蒋德明,陈长佑,杨嘉林,等.高等车用内燃机原理高辉,吴筱敏,高忠权,等.盲源分离法在火花塞高子电洸信M].西安:西安交通大学出版社,2006号分离中的应用.2010,44(1):22-26.[doi:10.7652/ xjtuxb[15] KALMTHOUT E V, VEYNANTE D, CANDEL S.201001006Direct numerical simulation analysis of flame surface(编辑苗凌)中国煤化工CNMHGhttp∥www.jdxb.cnhttp∥zkxb.xjtu.edu.cn