Jet-A表征燃料_乙醇掺混燃料层流燃烧特性研究_何旭.pdf

Jet-A 表征燃料/乙醇掺混燃料层流燃烧特性研究何旭1,2，江正晖1，桑正1，刘泽昌1，冯光远1，杨青1,2，蒋厚实1（1.北京理工大学 机械与车辆学院，北京100081；2.北京电动车辆协同创新中心，北京100081）摘 要：为了研究初始温度 470 K、初始压力 0.1 MPa 下乙醇掺混对航空煤油层流火焰燃烧特性的影响，文中利用定容燃烧弹设备结合仿真进行对比分析。结果表明，乙醇的加入显著提升了航空煤油的层流火焰燃烧速度，但仿真结果与试验数据之间存在较大的偏差，在低当量比时尤为严重，需进一步优化机理。在此基础上通过敏感性分析，发现对层流火焰燃烧速度影响较大的 4 个基元反应，通过对 4 个基元反应的参数进行调整，修改得到新模型，该模型与试验值存在较好的吻合度。关键词：航空煤油；乙醇；层流火焰燃烧速度；敏感性分析中图分类号：TK421.2 文献标志码：A 文章编号：1001-0645(2023)07-0712-07DOI：10.15918/j.tbit1001-0645.2022.174Laminar Combustion Characteristics of Jet-ASurrogate Fuel/EthanolHE Xu1,2，JIANG Zhenghui1，SANG Zheng1，LIU Zechang1，FENG Guangyuan1，YANG Qing1,2，JIANG Houshi1（1.School of Mechanical Engineering,Beijing Institute of Technology,Beijing 100081,China;2.Collaborative Innovation Center of Electric Vehicles in Beijing,Beijing 100081,China）Abstract：To study the effect of ethanol blending on the combustion characteristics of aviation kerosene laminarflame at an initial temperature of 470 K and an initial pressure of 0.1 MPa,a constant volume incendiary bombwas used in this paper to conduct a comparative analysis combined with simulation.The results show that the ad-dition of ethanol significantly improves the laminar burning velocity of aviation kerosene,but there is a large de-viation between the simulation results and the experimental data,especially at low equivalence ratios,and themechanism needs to be further optimized.On this basis,a sensitivity analysis was carried out,finding a great in-fluence of four elementary reactions on the laminar burning speed.Adjusting the parameters of the four element-ary reactions,a new model was obtained,being in good agreement with the experimental values.Key words：aviation kerosene；ethanol；laminar flame burning speed；sensitivity analysis 随着我国能源需求的不断上涨、能源进口依赖程度的不断加深以及双碳行动方案的印发1 2，对传统燃料清洁化、可再生化已是迫在眉睫.航空煤油作为传统的碳氢燃料，燃烧会产生大量的 CO、NOx、碳烟等污染物3 4.乙醇作为传统的生物质含氧燃料可通过生物质发酵等方式获得具有可再生性，具有较高的火焰速度和较低的未燃 HC、NOx等污染物排放，但其较低的能量密度和较高的水溶性限制了其完全替代航空煤油的潜力，因此在航空煤油中掺混乙醇是一种较为实际有效的方案5.由于航空煤油通常为 C7-C16烃类混合物6，不同产地、不同批次的煤油组分都会存在很大差异；另一方面，实际航空发动机在工作过程中，涉及到几百种组分、成千上万种基元反应，同时还伴随剧烈的 收稿日期：2022 10 24基金项目：国家自然科学基金面上项目（21961122007）作者简介：何旭（1976），男，博士，副教授，E-mail：.通信作者：蒋厚实（1990），男，博士，助理研究员，E-mail：.第 43 卷第 7 期北 京 理 工 大 学 学 报Vol.43No.72023 年 7 月Transactions of Beijing Institute of TechnologyJul.2023流动、传热等物理和化学变化，直接模拟其真实燃烧过程、构建全面的反应动力学是非常复杂且不切实际的.对于上述问题，使用简单可靠的表征燃料来替代真实复杂的航空煤油来开展相关基础研究已经得到广泛的认可.另外，作为湍流燃烧基础的层流燃烧特性，不仅体现了燃料的放热能力等基本燃烧性能；而且可以用于验证燃料的化学反应机理，为湍流燃烧的深入研究提供重要数据.因此，燃料的层流燃烧特性是目前进行燃烧基础研究的一个重要课题.近年来国内外众多学者对航空煤油的表征燃料开展了大量的试验研究.肖保国等7在张若凌等8对航空煤油冷却模拟的基础上提出了摩尔分数为 79%正癸烷、13%三甲基环己烷和 8%乙基苯的三组分RP-3 表征燃料.有研究提出以摩尔分数为 14%正癸烷、10%正十二烷、10%甲苯、30%异十六烷、36%甲基环己烷五组分 RP-3 表征燃料并研究了其与乙醇掺混对预混火焰的影响，发现乙醇的化学效应对层流燃烧速度的影响最大9 12.DAGAUT 等13 14使用正癸烷、正丙苯、正丙基环己烷来表征航空煤油，提出了 209 种组分与 1 673 种反应的氧化机理成功还原了煤油-氧气-氮气的火焰结构.DOOLEY 等15 16基于实际航空煤油的碳氢比和辛烷值，提出了一种以摩尔分数为42.67%正癸烷、33.02%异辛烷、24.34%甲苯的三组分表征燃料，并通过试验对该三组分替代燃料的熄灭极限和点火特性进行了验证.由此可见，针对航空煤油的表征燃料层流燃烧特性已有大量的研究基础，但针对乙醇与其表征燃料掺混的研究尚较缺乏.因此，文中的目标是开展乙醇掺混航空煤油表征燃料的层流燃烧特性研究.1 试验系统与分析方法定容燃烧弹是目前最广泛使用的层流燃烧速度测量装置17 19，和其它装置相比仅需少量的反应物，且能承受更高的试验压力.通过高速纹影摄像系统，可以获取预混层流火焰半径随时间的变化规律，并且更容易捕捉球形火焰受拉伸效应的影响.文中将选用定容燃烧弹进行相关试验.1.1 试验装置文中使用的测试系统如图 1 所示，由容弹主体、温度控制系统、配气系统、点火系统、高速纹影摄像系统组成.容弹主体由不锈钢制成，呈球形，密封性能良好，容弹设计可承受4 MPa 压力、可耐受650 K 温度，容弹的内部容积为 33.5 L.两根直径为 1.5 mm 的钨钢电极间隔 1.5 mm，对称布置在容弹上下两侧.数据采集所使用的高速相机的帧率设置为 10 000 fps，采集图像分辨率为 1 0241 024.光源狭缝真空泵燃料注射口环境大气排气口反射镜压力变送器压力表主反射镜温度变送器温度表刀口高速摄影相机氮气同步信号触发操作及数据采集系统压力传感器合成空气图 1 层流燃烧测试系统示意图Fig.1 Schematic diagram of laminar combustion test 1.2 试验工况为了便于分析验证，文中采用 Dooley 的表征燃料配比（42.67%正癸烷、33.02%异辛烷、24.34%甲苯）与乙醇掺混来开展混合燃料的层流燃烧特性研究15 16.液体燃料在实际试验研究中多使用体积参数，第 7 期何旭等：Jet-A 表征燃料/乙醇掺混燃料层流燃烧特性研究713便于燃料的配制和使用，因此本研究中混合燃料乙醇的添加含量以体积参数为基准，乙醇掺混比 的定义：=VethanolVethanol+Vkerosene（1）VethanolVkerosene式中：和分别为配置燃料时乙醇和三组分表征燃料的体积.文中将用 E0、E30、E60 分别表示掺混比为 0%、30%、60%.试验工况参数如表 1所示.表 1 试验工况表Tab.1 Experimental initial condition参数名称数值试验温度T/K470试验压力pi/MPa0.1掺混比E0，E30，E60当量比0.8，0.9，1.0，1.1，1.2，1.3，1.4 1.3 试验过程试验中所用到的合成空气为实验级（21%N2；79%O2）纯度大于 99.99%，液体燃料（乙醇；甲苯；正癸烷；异辛烷）纯度大于 99.5%.试验操作如下：首先利用真空泵将容弹内部抽至 6 kPa 以下.使用微量进液器抽取该工况下所需体积的燃料，开启容弹上的燃料注入阀门，将燃料通过阀门注入容弹内腔，然后迅速关闭阀门.静置约 5 min，待液体燃料充分汽化后打开空气进气阀门，向容弹内腔缓慢地充入空气至预设压力值.静置 10 min，待混合物混合均匀后，使用点火系统同时触发电极点火和相机拍摄，记录预混火焰发展图像.拍摄结束后，打开排气阀门，用新鲜空气冲洗容弹，重复 3 次以保证容弹内部无燃料与废气残留.检查测试结果有效后存储数据，每个工况点重复 3 次.1.4 数据处理利用高速相机拍摄得到的火焰纹影图像，通过将火焰原始图像与背景图像相减，得到火焰轮廓，利用 canny 算子20对其进行边界扫描并通过球形火焰面积拟合得到火焰半径 Rf.如式(2)所示，通过求得的 Rf求导即可得到拉伸火焰传播速度 Sb.Sb=dRfdt（2）式中，t 为火焰传播时间.火焰的拉伸率 K 定义为K=1AdAdt=2RfdRfdt=2RfSb（3）对应的，火焰的曲率 为=2Rf（4）S0b拉伸火焰传播速度和火焰拉伸率之间的关系可以通过曲率模型21和准稳态模型22的非线性拟合公式外推得到无拉伸火焰传播速度.相关公式如下：Sb=S0b(1Lb)（5）(SbS0b)2ln(SbS0b)2=2LbKS0b（6）式中：Lb为马克斯坦长度，马克斯坦长度反映的是火焰速度对拉伸的敏感性.文中通过拟合曲线与试验数据的吻合度选取拟合方法中较为合适的一种.根据火焰前缘的质量守恒定律，层流燃烧速度SL可以表示为SL=buS0b（7）式中：b为易燃气体密度；u为未燃气体密度.文中使用 ChemKin17.0 中 PREMIX 模块，分别结合 LI 等23提出的一种四组分燃料（正庚烷、异辛烷、甲苯、乙醇）的简化机理，KELLEY 等24提出的正癸烷骨架机理来实现试验预测.通过对上述两种机理的整合，初步形成了包含 73 种组分和 289 步基元反应的正癸烷、异辛烷、甲苯、乙醇的四组分燃烧反应模型，以匹配层流燃烧速度.计算中考虑 Soret 效应，允许最大网格点数为 1 000，自适应网格数量设为 50，求解梯度与曲率均设为 0.03.1.5 误差分析在本研究的层流火焰速度测量过程中，存在以下不确定性来源，包括：燃料与空气纯度、注射器、压力表精度、温控系统精度、点火能量、热辐射与压力波、火焰数据后处理等因素，这些因素对层流火焰速度测量的不确定性的贡献度可使用一维模型得到评估25.本研究中层流燃烧速度的不确定性由 Moffat 理论26计算得到：USL=(GSL)2+(tM1SSLM)2（8）SSLSLGSL式中：M 为试验的总次数；t(M1)为自由度为 M1 的 t分布在置信区间为 95%所对应的值；为的标准差；系统不确定度可以通过下式计算：GSL=ni=1(SL(xi)xiOi)2（9）714北 京 理 工 大 学 学 报第 43 卷式中：xi为任何一个影响 SL不确定度的变量；Oi为变量 xi的误差（包含当量比，初始温度和压力27）.对此，燃料与空气的纯度均在 99%以上，压力表精度为0.1 kPa，热电偶精度为1 K，燃料当量比误差范围为1.02.4%；火焰图像最小像素点为 1，后处理程序对每帧图像的半径计算值误差约为 0.11 mm，同时考虑到重复试验过程中的随机误差与系统误差，层流火焰速度的试验测量的不确定度5%.T=1 ms2 ms3 ms4 ms5 ms6 ms7 msE0E30E60图 2 不同掺混比下火焰传播图像（p=0.1 MPa，T=470 K，=1.0）Fig.2 Flame propagation images under different blending ratios(p=0.1 MPa,T=470 K,=1.0)2 结果与讨论 2.1 不同掺混比对火焰特性的影响图2 显示了初始温度为470 K、初始压力为0.1 MPa、当量比为 1.0 时，乙醇的掺混比分别为 E0、E30、E60 时的层流火焰发展图像.对比掺混乙醇时的燃料火焰图像和未掺混乙醇时的火焰图像，可以明显地观察到添加乙醇时同一时刻的火焰半径更大，反映出乙醇这种含氧燃料的添加提升了航空煤油表征燃料的燃烧效果，使火焰传播速度加快.同时，在火焰的发展初期，乙醇的添加不会使火核发生显著变化，随着燃烧过程的不断进行，乙醇对燃烧的促进作用更加明显.图 3 为在初始温度为 470 K、初始压力为 0.1 MPa时，掺混比分别为 E0、E30、E60 下的火焰传播速度与火焰拉伸率的关系.图中扇形面积（两条虚线中间所夹的区域）表示所选择的火焰半径范围为 825 毫米.将规定范围内的数据通过拟合度最高的外推法得到无拉伸火焰传播速度.由图可知，随着乙醇掺混比的提升，未拉伸火焰传播速度总体呈现上升趋势.E30 与 E60 的未拉伸火焰传播速度相较 E0 分别平均提升了 5.59%与 10.63%.在当量比为 1.1 时，拉伸率为 0 的情况下，E30 与 E60 的未拉伸火焰传播速度分别为 4.5 m/s 与 4.67 m/s，相较 E0 分别提高了 4.9%与 8.9%.由此可见，乙醇的加入使火焰传播速度提升显著.02004006008001 0000.01.02.03.04.05.06.0E0R=8 mmR=25 mmSb/(ms1)K/s1(a)E0(b)E30(c)E6002004006008001 0000.01.02.03.04.05.06.0E30Sb/(ms1)K/s1 02004006008001 0000.01.02.03.04.05.06.0E60Sb/(ms1)K/s1 T=470 K p=0.1 MPaR=8 mmR=25 mmT=470 K p=0.1 MPa=0.8=0.9=1.0=1.1=1.2=1.3=1.4 Fitting=0.8=0.9=1.0=1.1=1.2=1.3=1.4 FittingR=8 mmR=25 mmT=470 K p=0.1 MPa=0.8=0.9=1.0=1.1=1.2=1.3=1.4 Fitting图 3 不同掺混比时的火焰传播速度与火焰拉伸率的关系Fig.3 The relationship between stretched flame propagation speed and stretch rate 通过式（7）即可将火焰传播速度转化为层流火焰燃烧速度，如图 4 所示，在试验误差范围内，试验结果与仿真结果的发展趋势基本相同，即在当量比为 0.71.4 之间时，层流火焰燃烧速度先增大后减第 7 期何旭等：Jet-A 表征燃料/乙醇掺混燃料层流燃烧特性研究715=1.1小，层流火焰燃烧速度最大值出现在略微高于化学计量数一侧，本试验中的速度最大处为时，不同的乙醇掺混比不影响这一变化趋势.随着乙醇掺混比例的提高，层流火焰燃烧速度增大.然而仿真结果与试验结果在数值上存在较大偏差，模型计算值整体偏低，低当量比的数据相差最大.为进一步分析原因，文中将对层流火焰燃烧速度进行敏感性分析.由于模型计算值与试验值之间的差异主要反映在不同的当量比上，因此文中需要确定不同当量比下影响多组分燃料的层流火焰燃烧速度的关键基元反应.0.70.80.91.01.11.21.31.41.5405060708090T=470 Kp=0.1 MPaSL/(cms1)当量比 E0 E30 E60点:实验数据线:仿真值图 4 T=470 K、p=0.1 MPa 时不同掺混比下的层流火焰燃烧速度随当量比的变化Fig.4 Laminar burning velocities with equivalent ratio of various blendingratios when T=470 K、p=0.1 MPa 文中需要确定不同当量比下影响多组分燃料的层流火焰燃烧速度的关键基元反应.2.2 敏感性分析在仿真软件中通过对层流燃烧速度的敏感性分析，（即为通过给某一初始参数以微小变量x，研究目标参量（本研究为层流燃烧速度）的变化量y），可以评价不同的基元反应对温度、压力、掺混比等参数的敏感性，从而发现主导层流燃烧速度的关键基元反应，并对主要基元反应进行修订.无量纲质量流量的敏感系数如下：Si=kimomoki（10）式中：ki为第 i 个元素反应的化学反应速率；mo为质量燃烧流量.文中使用敏感性分析方法，来确认不同当量比下影响 E0 的层流火焰燃烧速度的关键基元反应，结果如下图 5 所示.其中：反应 R67、R86、R162 的正敏感性系数最大，表明这 3 个反应对燃烧过程的影响作用最为显著，这主要是因为 R67 和 R162 是非常重要的链式分支反应，对反应的起始阶段有着极大的影响.相反，R84 是一个典型的链终止反应，该反应通过消耗反应过程中的 H 自由基，降低了反应中 H 自由基的浓度，对反应起抑制作用.R203、R32 同样会消耗 H 自由基进而抑制反应的进程.R89 这类反应通过将自由基转化为稳定的产物如 CH4，抑制了层流火焰燃烧速度.随着当量比的增加，R67、R190、R36 的正敏感系数逐渐增大，说明其对层流火焰燃烧速度的影响作用在增强，高当量比下对层流火焰燃烧速度能产生更大的影响；相反，R86 基元反应在低当量比下对层流火焰燃烧速度的影响作用更加显著.对于负敏感系数下的反应，R84 能抑制低当量比下层流燃烧速度，高燃料浓度时的抑制效果明显减弱.R203、R89、R32 这三个反应在高当量比时对层流火焰燃烧速度的抑制作用最为显著，在低当量比时几乎不对火焰速度产生影响.R68 和 R242 对燃料的浓度不敏感，不同当量比下的负敏感性系数基本相同，对层流火焰燃烧速度的抑制作用无明显变化.2.3 机理修正机理模型的修改参数主要是基于化学动力学中用来定量地描述基元反应过程速率的阿伦尼乌斯公式三参数阿伦尼乌斯公式为K=ATneEa/RT（11）式中：K 为速率常数；R 为摩尔气体常量；T 为热力学温度；Ea为表观活化能；A 为指前因子（也称频率因子）.综合上述敏感性分析，结合不同当量比下试验值与仿真的拟合结果，可以得知需要提升部分正敏 0.20.10.00.10.20.30.4R67 O2+H=O+OHR86 CO+OH=CO2+HR162 HCO+M=H+CO+MR190 C2H2+H(+M)=C2H3(+M)R36 C6H5+O2=C6H5O+OR84 H+O2(+M)=HO2(+M)R203 C2H3+H=C2H2+H2R163 HCO+O2=CO+HO2R32 C6H5+H(+M)=C6H6(+M)R89 CH3+H(+M)=CH4(+M)R68 H+OH+M=H2O+M敏感度=1.3=1.1=0.9R242 C3H6+H=C3H5-A+H2E0T=470 Kp=0.1 MPa图 5 初始机理的敏感性分析Fig.5 Sensitivity analysis of the initial mechanism716北 京 理 工 大 学 学 报第 43 卷感系数的反应，使得层流火焰燃烧速度得到提升.同时，需要对部分负敏感系数的基元反应进行修正，使得整体层流火焰燃烧速度控制在合理范围.对此，参照多位学者28 32的相关研究，对原始整合机理中的部分反应的系数（A，n，Ea）进行了修正，表 2 显示了主要基元反应的参数.如图 6 所示，修正后的模型整体上有效提升了各工况下的层流火焰燃烧速度，低当量比时的火焰速度提升幅度最大，但随着乙醇掺混比的增高吻合度有所下降，高当量比下的速度提升微弱，在当量比处于化学计量数附近时，火焰速度的仿真结果略微高于试验值.总的来说，修正后的新模型计算得到的层流火焰燃烧速度与试验数据的拟合度更高，这与前文中有关基元反应的敏感性分析一致，符合预期结果.表 2 修正前后主要基元反应系数的对比Tab.2 Comparison of Modified Primitive Response Coefficients基元反应原模型参数修正后的模型参数参考文献AnEaAnEaCO+OH=CO2+H7.021042.1355.71.401051.92 347.727HCO+O2=CO+HO27.5810120.0410.03.0010120.00.028C3H6+H=C3H5-A+H23.641052.54 361.25.0010120.01 505.729H+O2(+M)=HO2(+M)4.6510120.440.05.8010120.40.030 0.70.80.91.01.11.21.31.41.55060708090SL/(cms1)当量比 E0 E30 E60点:实验值线:仿真值T=470 Kp=0.1 MPa图 6 修正后的机理与试验数据的对比Fig.6 Comparison of the revised mechanism with experimental data 3 结论文中讨论了采用定容燃烧弹测量初始温度 470 K，初始压力0.1 MP，掺混比为E0、E30、E60，当量比0.81.4 下航空煤油表征燃料与乙醇掺混的层流火焰.通过分析发现：乙醇的加入显著提升了航空煤油的层流火焰燃烧速度.通过现有的整合的机理发现仿真值与试验数据的吻合度较差，通过对当量比的敏感性分析，对4 个基元反应进行参数调整，修改得到新机理模型，新机理模型计算得到的层流火焰燃烧速度与试验数据的拟合度更高.参考文献：肖献法.国务院2030年前碳达峰行动方案给交通运输、1 汽车等行业定调J.商用汽车,2021(11):16 22.XIAO Xianfa.The state councils action plan for carbondioxide peak before 2030 sets the tone for transportation,automobile and other industriesJ.Commercial Vehicle,2021(11):16 22.(in Chinese)赵琳,刘文峰,李斌,等.政策导向下的新能源公交车节能减排效果评估J.北京理工大学学报,2016,36(增刊2):99 102.ZHAO Lin,LIU Wenfeng,LI Bin,et al.Evaluation on energy-saving and emission-reduction results of new energy busesbased on policy guidanceJ.Transactions of Beijing instituteof Technology,2016,36(suppl 2):99 102.(in Chinese)2 何旭,张志鹏,吴昊,等.基于二维激光诱导炽光法的棉籽油扩散火焰碳烟生成特性J.北京理工大学学报,2019,39(3):235 240+326.HE Xu,ZHANG Zhipeng,WU Hao,et al.Investigation on thesoot formation of cottonseed oil diffusion flame by two-dimensional laser induced incandescenceJ.Transactions ofBeijing institute of Technology,2019,39(3):235 240+326.(in Chinese)3 FOONG T M,MORGANTI K J,BREAR M J,et al.The effectof charge cooling on the ron of ethanol/gasoline blendsJ.SaeInternational Journal of Fuels&Lubricants,2013,6(1):34 43.4 YANG Q,LIU Z,HOU X,et al.Measurements of laminarflame speeds and flame instability analysis of E30-air premixedflames at elevated temperatures and pressuresJ.Fuel,2020,259:116223.5 范学军,俞刚.大庆RP-3航空煤油热物性分析J.推进技术,2006,27(2):187 192.FANG Xuejun,YU Gang.Analysis of the mophysicalproperties of Daqing RP-3 aviation keroseneJ.Journal ofPropulsion TechnologY,2006,27(2):187 192.(in Chinese)6 肖保国,杨顺华,赵慧勇,等.RP-3航空煤油燃烧的详细和简 7 第 7 期何旭等：Jet-A 表征燃料/乙醇掺混燃料层流燃烧特性研究717化化学动力学模型J.航空动力学报,2010,25(9):1948 1955.XIAO Baoguo,YANG Shunhua,ZHAO Huiyong,et al.Detailed and reduced chemical kinetic mechanisms for RP-3aviation kerosene combustionJ.Journal of AeronauticalDynamics,2010,25(9):1948 1955.(in Chinese)ZHANG R L,JIN J,LE J L.The simulation of endothermicfuel flow in cooling channels of Scramjet J.8 曾文,刘靖,张治博,等.一种新的RP-3航空煤油模拟替代燃料J.航空动力学报,2017,32(10):2314 2320.ZENG Weng,LIU Jing,ZHANG Zhibo,et al.A new surrogatefuel of RP-3 keroseneJ.Journal of Aeronautical Dynamics,2017,32(10):2314 2320.(in Chinese)9 曾文,陈欣,马洪安,等.RP-3航空煤油层流燃烧特性的实验J.航空动力学报,2015,30(12):2888 2896.ZENG Weng,CHEN Xin,MA Hongan,et al.Experiment onlaminar combustion characteristics of RP-3 keroseneJ.Journal of Aeronautical Dynamics,2015,30(12):2888 2896.(in Chinese)10 MA H,XIE M,ZENG W,et al.Experimental study oncombustion characteristics of Chinese RP-3 keroseneJ.Chinese Journal of Aeronautics,2016,29(2):375 385.11 LIU J,HU E,ZENG W,et al.A new surrogate fuel foremulating the physical and chemical properties of RP-3keroseneJ.Fuel,2020,259:116210.12 DAGAUT P.On the kinetics of hydrocarbons oxidation fromnatural gas to kerosene and diesel fuelJ.Physical ChemistryChemical Physics,2002,4(11):2079 2094.13 DAGAUT P,EL BAKALI A,RISTORI A.The combustion ofkerosene:Experimental results and kinetic modelling using 1-to 3-component surrogate model fuelsJ.Fuel,2006,85(7-8):944 956.14 DOOLEY S,WON S H,CHAOS M,et al.A jet fuel surrogateformulated by real fuel propertiesJ.Combustion and Flame,2010,157(12):2333 2339.15 DOOLEY S,WON S H,HEYNE J,et al.The experimentalevaluation of a methodology for surrogate fuel formulation toemulate gas phase combustion kinetic phenomenaJ.Combustion and Flame,2012,159(4):1444 1466.16 XIOURIS C,YE T,JAYACHANDRAN J,et al.Laminarflame speeds under engine-relevant conditions:Uncertaintyquantification and minimization in spherically expandingflame experimentsJ.Combustion&Flame,2016,163:270 283.17 GIANNAKOPOULOS G K,GATZOULIS A,FROUZAKISC E,et al.Consistent definitions of“flame displacementspeed”and“markstein length”for premixed flamepropagationJ.Combustion&Flame,2015,162(4):1249 1264.18 CLAVIN P.Dynamic behavior of premixed flame fronts inlaminar and turbulent flowsJ.Progress in Energy and 19 Combustion Science,1985,11(1):1 59.张艳群,孟凡荣.MATLAB在图像边缘检测中的应用J.计算机应用研究,2004(6):144 146.ZHANG Yanqun,MENG Fanrong.Application of MATLABin I mage edge detectionJ.Application Research ofComputers,2004(6):144 146.20 MARKSTEIN G H.Experimental and Theoretical Studies ofFlame-Front Stability M.Dynamics of Curved Fronts,S.l.:Cornell Aeronautical Laboratory,Inc.,1988:413 423.21 KELLEY A P,LAW C K.Nonlinear effects in the extractionof laminar flame speeds from expanding spherical flamesJ.Combustion&Flame,2009,156(9):1844 1851.22 LI Y,ALFAZAZI A,MOHAN B,et al.Development of areduced four-component(toluene/n-heptane/iso-octane/ethanol)gasoline surrogate modelJ.Fuel,2019,247:164 178.23 YU W,YANG W,TAY K,et al.Development of a newskeletal mechanism for decalin oxidation under enginerelevant conditionsJ.Fuel,2018,212:41 48.24 ZHENG C.On the accuracy of laminar flame speedsmeasured from outwardly propagating spherical flames:Methane/air at normal temperature and pressureJ.Combustion&Flame,2015,162(6):2442 2453.25 MOFFAT R J.Describing the uncertainties in experimentalresultsJ.Experimental Thermal&Fluid Science,1988,1(1):3 17.26 LI Q,HU E,YU C,et al.Measurements of laminar flamespeeds and flame instability analysis of 2-methyl-1-butanolair mixturesJ.Fuel,2013,112(3):263 271.27 FISHER E M,PITZ W J,CURRAN H J,et al.Detailedchemical kinetic mechanisms for combustion of oxygenatedfuelsJ.Proceedings of the Combustion Institute,2000,28(2):1579 1586.28 RANZI E,FRASSOLDATI A,STAGNI A,et al.Reducedkinetic schemes of complex reaction systems:fossil andbiomass-derived transportation fuelsJ.International Journalof Chemical Kinetics,2014,46(9):512 542.29 钟北京,姚文斐.基于特征值分析的正癸烷骨架和总包简化机理J.物理化学学报,2014,46(9):512 542.ZHONG Beijing,YAO W