2. 中国空气动力研究与发展中心超高速空气动力研究所高超声速冲压发动机技术重点实验室, 四川绵阳 621000
2. Science and Technology on Scramjet Laboratory, Hypervelocity Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China
燃烧过程是各类动力推进系统工作的核心要素之一, 因此燃烧研究至关重要[1]. 要想实现高效、稳定的燃烧组织, 就必须对燃烧过程中涉及的物理、化学反应过程有清晰的认识, 对组分浓度、温度、压力、速度等燃烧流场信息进行诊断和测量[2]. 但是燃烧过程往往涉及复杂的化学反应及其与流动之间的相互耦合, 传统的燃烧诊断技术, 包括Pitot管、热电偶、测压传感器等, 都属于侵入式诊断技术, 对于燃烧流场干扰较大, 具有很大的局限性, 已经不能满足燃烧学研究对燃烧诊断和流场可视化的要求. 随着激光技术的发展, 激光燃烧诊断技术得以快速发展, 由于具有对流场无干扰、时空分辨率高等优势, 已经成为燃烧实验研究的主要手段[3-4].
平面激光诱导荧光技术(planar laser-induced fluorescence, PLIF)是一种非接触式的燃烧诊断技术[5], 能够探测分子、自由基等粒子的浓度[6-7], 测量燃烧场温度[8-10], 表征复杂燃烧场, 实现火焰结构的可视化[11-13]. 该技术具有高灵敏度、高时空分辨率、高选择性以及原位无干扰等特性, 已经成为研究燃烧和火焰结构的重要工具, 在基础燃烧研究[14-15]、超声速燃烧[16-17]、航空发动机[18]等领域得到了广泛应用, 并发挥着越来越重要的作用. 在碳氢燃料燃烧过程中, OH自由基是一种关键的中间产物, 参与大多数的基元反应, 通常用于表征燃烧反应区, 对于反映火焰锋面位置、火焰燃烧反应区等具有重要意义. 而且, OH具有浓度较高、吸收较强、辐射谱线位于紫外波段等特点, 能够较为容易地被检测. 因此, OH是PLIF技术中最重要也是最普遍的一个测量对象.
实际燃烧器, 不管是燃气轮机, 还是航空发动机、超燃冲压发动机, 都是在高压条件下工作的. 虽然OH-PLIF技术在这些实际燃烧器的研究中发挥着重要作用, 但是在高压条件下开展OH-PLIF测量非常具有挑战性[19]. 一个显著的问题是, 随着压强的增加, OH荧光猝灭效应显著增强, 导致荧光信号迅速减弱. 由于OH-PLIF物理机理比较复杂, 而且高压燃烧实验比较困难, 因此量化压强对OH-PLIF的影响仍是十分困难的.
目前, 开展高压燃烧研究的实验装置主要有两大类. 一类是利用高压定容燃烧弹[20-21], 这类燃烧装置的燃烧时间极短, 火焰迅速变化, 不适合开展光学诊断测量研究. 二是采用连续燃烧炉, 如本生灯[22]、旋流燃烧器[23-24]、平面火焰炉[25-27]等. 这类高压燃烧装置可以实现连续稳定的燃烧, 适合采用光学诊断技术开展研究. 为了开展光学诊断技术研究和参数校正, 本课题组自行设计了一个高压火焰炉, 并采用滤波Rayleigh散射技术, 进行了不同压强条件下甲烷/空气预混火焰温度场测量研究[25, 27].
本文将在高压火焰炉上开展甲烷/空气预混平面火焰OH-PLIF测量研究, 具体内容安排如下: 第1节介绍实验装置, 分别对高压火焰炉和OH-PLIF测量系统进行详细描述, 第2节对实验结果进行分析讨论, 第3节对实验工况对应的荧光产额进行仿真, 并与实验结果进行对比分析, 第4节对本文内容进行总结.
1 实验装置 1.1 高压火焰炉课题组自主设计建设了一套高压火焰炉系统, 该系统由设备本体、气源系统、冷却系统、测控系统4部分组成. 设备能够承受的极限压强可达1.5 MPa, 压强在0.1~1 MPa之间连续可调, 压强控制精度优于0.01 MPa. 高压火焰炉的结构示意图如图 1(a)所示, 高压腔室的四面均安装有石英玻璃光学窗口, 便于进行光学测量研究. 高压火焰炉采用氮气作为工作气体来调节压强, 同时氮气也用作火焰保护气体. 火焰炉采用McKenna平面火焰炉, 炉面直径为25 mm. 燃烧气体采用甲烷/空气预混气体, 通过调节炉内压力可以实现不同压强(0.1~1.0 MPa)条件下甲烷/空气层流预混火焰的稳定燃烧, 其中常压，0.34，0.64，0.98 MPa条件下稳定燃烧的平面火焰照片如图 1(b)~(e)所示.
|
| 图 1 高压火焰炉结构示意图及不同压强条件下稳定燃烧的平面火焰图像 Fig.1 Schematic diagram of the high-pressure burner and photos of the flat flames at different pressures |
在本文的研究中, 甲烷/空气的燃烧当量比设置为1.0, 具体实验工况如表 1所示. 高压火焰炉的压强分别设置为0.1~1.0 MPa, 间隔是0.1 MPa. 为了保持平面火焰炉的出口速度不变, 甲烷和空气的流量随着压强的增加而成比例地增加. McKenna平面火焰炉表面的多孔板中通有冷却水管, 一方面防止燃烧器过热, 保证火焰炉的安全, 另一方面保证预混气体的出口温度稳定在约300 K.
|
下载CSV 表 1 实验工况 Tab.1 Experimental conditions |
为了研究压强对OH-PLIF的影响, 本文搭建了OH-PLIF测量实验系统, 研究对象是层流预混火焰, 如图 2所示. 激光器系统由Nd∶YAG固体激光器、染料激光器(Dye)和紫外倍频器(SHG)这3部分组成. Nd∶YAG固体激光器输出的光经倍频后产生532 nm脉冲激光, 进入染料激光器用作泵浦光, 染料激光器输出的光经SHG倍频后输出紫外激光. 紫外激光的脉宽约为10 ns, 重复频率为10 Hz, 能量约为5 mJ, 线宽约为0.3 cm-1. 根据LIFBase数据库[28], OH的激发谱线选用A2Σ+-X2Π(1, 0)电子跃迁带的P1(1)谱线, 对应波长为282.254 nm. 染料激光器输出的光采用一个波长计来监测波长, 以保证紫外激光能够恰好落在OH的激发谱线上. 实验采用McKenna平面火焰炉, 放置在高压火焰炉中. 紫外激光通过透镜组后, 形成约15 mm高的激光片光, 片光沿着平面火焰炉中心线竖直地穿过火焰. OH荧光信号由ICCD相机接收, 在镜头前装有滤光片(Semrock FF02-320/40), 以减少火焰辐射及其他杂散光的干扰. 该滤波片对应的通带波段为300~340 nm, 覆盖了OH的A-X(0, 0)和A-X(1, 1) 电子跃迁带的荧光. 采用DG535型数字脉冲发生器来控制Nd∶YAG激光器的触发信号、Q开关和ICCD相机的时序. 激光器和ICCD快门的时序控制如图 3所示, 激光器和ICCD快门是同一时刻触发的; ICCD的快门时间设置为200 ns, 以消除火焰自发辐射的影响.
|
| 图 2 高压层流预混火焰OH-PLIF测量系统示意图 Fig.2 Experimental setup of OH-PLIF measurement for the flat flames at high pressures |
|
| 图 3 激光器、ICCD时序控制示意图 Fig.3 Schematic diagram of laser and ICCD timing control |
实验按照表 1所列的每个实验工况, 设置高压火焰炉的压强和平面火焰的甲烷、空气流量. 对于每个实验工况，待高压火焰炉内压强和平面火焰稳定后, 采用PLIF技术连续拍摄50幅OH荧光瞬态图像, 将每幅图像左上角没有信号的100×100区域内平均像素值作为基底, 对每幅OH荧光图像进行背景扣除, 然后对这50幅扣除背景的图像进行数据平均处理, 得到不同压强条件下甲烷/空气预混层流火焰OH-PLIF平均图像, 如图 4所示. 需要指出的是, 图 4中每个子图的色标都是按照0.1 MPa对应的OH-PLIF图像最大值进行归一化的. 从图 3中可以看出, 随着压强的增加, OH的荧光信号越来越弱, 这是由于OH的荧光猝灭效应随压强增加而增强.
|
| 图 4 不同压强条件下甲烷/空气预混层流火焰OH-PLIF平均测量结果 Fig.4 Averaged OH-PLIF measurement results of laminar premixed CH4/air flames at different |
从图 4所示的OH-PLIF平均图像沿着火焰炉的中心线提取得到OH荧光强度, 得到不同压强条件OH荧光强度随着高度的变化曲线, 如图 5所示. 由于平面火焰可以近似看成是一维火焰, 而一维火焰的火焰面比较薄, 因此靠近炉面附近, OH荧光信号随OH浓度的增加而迅速上升. 在本实验中火焰并不是绝热的, 火焰与外围的氮气存在热交换, OH浓度在到达峰值后将会逐渐减低直至零, 相应地, OH荧光信号逐渐下降, 直至消失. 因为压强的增加会提高甲烷/空气的燃烧反应速率, 所以火焰锋面的厚度会越来越薄. 于是, OH荧光信号的峰值位置会随压强的增加而越来越靠近炉面. 同样由于猝灭效应, OH荧光信号随压强的增加而减弱.
|
| 图 5 火焰炉中心线上OH荧光强度沿着高度方向的分布图 Fig.5 Fluorescence intensity of OH vs height along the centerline of the burner |
首先, 采用开源软件Cantera[29]和GRI Mech3.0化学反应机理[30], 在绝热近似条件下, 计算表 1中所列工况对应的甲烷/空气火焰达到平衡状态时的气体组分和温度. 再采用LASKIN代码[31-32], 针对本实验的OH-PLIF探测采集方案(即激发谱线采用A-X(1, 0)频带的P1(1)线, 荧光信号包含OH A-X(0, 0)和A-X(1, 1)频带的总辐射), 结合仿真得到的甲烷/空气火焰气体组分和温度, 对荧光产额进行仿真, 得到不同压强下甲烷/空气火焰的荧光产额, 见图 6中实心正方形标记的曲线. 在图 5所示的OH荧光强度分布曲线上提取最大值, 得到OH峰值荧光信号随着压强的变化关系, 见图 6中实心圆标记的曲线.
|
| 图 6 不同压强条件下OH荧光峰值信号与荧光产额对比图 Fig.6 Comparison of peak fluorescence intensity and simulated fluorescence yield of OH at different pressures |
从荧光产额与压强的关系曲线看, 荧光产额随着压强的增加而迅速降低, 两者近似呈指数衰减关系. 1 MPa对应的荧光产额仅为0.1 MPa的1/10, 增大压力还会进一步降低荧光产额. 这也是高压条件下OH-LIF测量比较困难的原因.
从图 6中两条曲线对比来看, 实测的甲烷/空气平面火焰OH荧光强度随压强的变化趋势, 大体与仿真的荧光产额随压强的变化趋势一致. 这表明, 虽然OH荧光强度与重叠因子密切相关, 进而与压强造成的OH谱线展宽、频移等因素相关, 但是荧光淬灭效应(荧光产额背后的物理机制)在OH荧光强度-压强关系中起着决定性的作用.
4 结论本文在自行设计的高压火焰炉上, 搭建了一套OH-PLIF测量系统, 采用OH A-X(1, 0)电子跃迁带的P1(1)谱线, 对不同压强条件下甲烷/空气预混平面火焰进行OH-PLIF测量研究, 并对实验工况对应的荧光产额进行了仿真, 将仿真结果与实验结果进行了对比分析. 研究结果表明: (1)随着压强的增加, OH的荧光信号逐渐减弱; (2)平面预混火焰的锋面厚度会随压强增加而逐渐变薄, OH荧光信号的峰值位置会越来越靠近炉面; (3)荧光淬灭效应在OH荧光强度-压强关系中起着决定性的作用.
| [1] |
甯红波, 李泽荣, 李象远. 燃烧反应动力学研究进展[J]. 物理化学学报, 2016, 32(1): 131-153. Ning H B, Li Z R, Li X Y. Progress in combustion kinetics[J]. Acta Physico-Chimica Sinica, 2016, 32(1): 131-153. (in Chinese) |
| [2] |
Aldén M, Bood J, Li Z S, et al. Visualization and understanding of combustion processes using spatially and temporally resolved laser diagnostic techniques[J]. Proceedings of the Combustion Institute, 2011, 33(1): 69-97. DOI:10.1016/j.proci.2010.09.004 |
| [3] |
Egon P H, Linow S. Laser diagnostics for studies of turbulent combustion[J]. Measurement Science and Technology, 2000, 11(2): R37-R57. DOI:10.1088/0957-0233/11/2/201 |
| [4] |
Cheskis S, Goldman A. Laser diagnostics of trace species in low-pressure flat flame[J]. Progress in Energy and Combustion Science, 2009, 35(4): 365-382. DOI:10.1016/j.pecs.2009.02.001 |
| [5] |
Grisch F, Orain M. Role of planar laser-induced fluorescence in combustion research[J]. Journal of Aerospace Lab, 2009(1): AL01-11. |
| [6] |
Yu X, Yang Z, Peng J B, et al. Quantitative measurements of one-dimensional OH absolute concentration profiles in a methane/air flat flame by bi-directional laser-induced fluorescence[J]. Chinese Physics B, 2015, 24(11): 114204. DOI:10.1088/1674-1056/24/11/114204 |
| [7] |
张蕾, 孙锐, 杨振, 等. 一维平面火焰OH基绝对浓度值的测量与分析[J]. 燃烧科学与技术, 2017, 23(2): 161-165. Zhang L, Sun R, Yang Z, et al. Measurement and analysis of 1-D OH radicals absolute concentration in flat flame[J]. Journal of Combustion Science and Technology, 2017, 23(2): 161-165. (in Chinese) |
| [8] |
Chen S, Su T, Yang F R, et al. Calibration method for 2D instantaneous OH-PLIF temperature measurements in flame[J]. Chinese Optics Letters, 2013, 11(5): 053001. DOI:10.3788/COL201311.053001 |
| [9] |
苏铁, 陈爽, 杨富荣, 等. 双色平面激光诱导荧光瞬态燃烧场测温实验[J]. 红外与激光工程, 2014, 43(6): 1750-1754. Su T, Chen S, Yang F R, et al. Investigation of temperature of transient combustion using two-line PLIF[J]. Infrared and Laser Engineering, 2014, 43(6): 1750-1754. DOI:10.3969/j.issn.1007-2276.2014.06.010 (in Chinese) |
| [10] |
Yang Z, Peng J B, Yu X, et al. Two-dimensional temperature characteristics of syngas flames with nitrogen, carbon dioxide, and steam diluents by planar laser-induced fluorescence thermometry[J]. Spectroscopy Letters, 2016, 49(7): 482-490. DOI:10.1080/00387010.2016.1208249 |
| [11] |
朱家健, 赵国焱, 龙铁汉, 等. OH和CH2O平面激光诱导荧光同时成像火焰结构[J]. 实验流体力学, 2016, 30(5): 55-60, 87. Zhu J J, Zhao G Y, Long T H, et al. Simultaneous OH and CH2O PLIF imaging of flame structures[J]. Journal of Experiments in Fluid Mechanics, 2016, 30(5): 55-60, 87. (in Chinese) |
| [12] |
Zhou B, Brackmann C, Li Z S, et al. Simultaneous multi-species and temperature visualization of premixed flames in the distributed reaction zone regime[J]. Proceedings of the Combustion Institute, 2015, 35(2): 1409-1416. DOI:10.1016/j.proci.2014.06.107 |
| [13] |
Butz D, Hartl S, Popp S, et al. Local flame structure analysis in turbulent CH4/air flames with multi-regime characteristics[J]. Combustion and Flame, 2019, 210: 426-438. DOI:10.1016/j.combustflame.2019.08.032 |
| [14] |
Skiba A W, Carter C D, Hammack S D, et al. A simplified approach to simultaneous multi-scalar imaging in turbulent flames[J]. Combustion and Flame, 2018, 189: 207-211. DOI:10.1016/j.combustflame.2017.10.035 |
| [15] |
张猛, 俞森彬, 金武, 等. 预混湍流3D火焰面密度和燃料消耗率的估计[J]. 工程热物理学报, 2015, 36(7): 1578-1583. Zhang M, Yu S B, Jin W, et al. Experimental study on turbulent premixed flames of syngas oxyfuel mixtures at high pressure with OH-PLIF[J]. Journal of Engineering Thermophysics, 2015, 36(7): 1578-1583. (in Chinese) |
| [16] |
McRae C D, Johansen C T, Danehy P M, et al. Image analysis of hydroxyl-radical planar laser-induced fluorescence in turbulent supersonic combustion[J]. Journal of Propulsion and Power, 2016, 32(3): 542-559. DOI:10.2514/1.B35611 |
| [17] |
Johansen C T, Mcrae C D, Danehy P M, et al. OH PLIF visualization of the UVa supersonic combustion experi-ment: configuration A[J]. Journal of Visualization, 2014, 17(2): 131-141. DOI:10.1007/s12650-014-0197-2 |
| [18] |
胡志云, 叶景峰, 张振荣, 等. 航空发动机地面试验激光燃烧诊断技术研究进展[J]. 实验流体力学, 2018, 32(1): 33-42. Hu Z Y, Ye J F, Zhang Z R, et al. Development of laser combustion diagnostic techniques for ground aero-engine testing[J]. Journal of Experiments in Fluid Mechanics, 2018, 32(1): 33-42. (in Chinese) |
| [19] |
Matynia A, Molet J, Roche C, et al. Measurement of OH concentration profiles by laser diagnostics and modeling in high-pressure counterflow premixed methane/air and biogas/air flames[J]. Combustion and Flame, 2012, 159(11): 3300-3311. DOI:10.1016/j.combustflame.2012.06.013 |
| [20] |
Yang Q, Liu Z C, Hou X H, et al. Measurements of laminar flame speeds and flame instability analysis of E30-air premixed flames at elevated temperatures and pres-sures[J]. Fuel, 2020, 259: 116223. DOI:10.1016/j.fuel.2019.116223 |
| [21] |
Duva B C, Wang Y C, Chance L E, et al. Correlations for the laminar burning velocity and burned gas Markstein length of methane-air mixtures diluted with flue gases at high temperatures and pressures[J]. Fuel, 2020, 281: 118721. DOI:10.1016/j.fuel.2020.118721 |
| [22] |
王金华, 聂要辉, 余芊芊, 等. 高压下湍流火焰速度测量与拟合[J]. 工程热物理学报, 2016, 37(4): 902-906. Wang J H, Nie Y H, Yu Q Q, et al. Measurement and correlation of turbulent burning velocity under high pressure[J]. Journal of Engineering Thermophysics, 2016, 37(4): 902-906. (in Chinese) |
| [23] |
Boxx I, Slabaugh C, Kutne P, et al. 3 kHz PIV/OH-PLIF measurements in a gas turbine combustor at eleva-ted pressure[J]. Proceedings of the Combustion Insti-tute, 2015, 35(3): 3793-3802. |
| [24] |
Sadanandan R, Kutne P, Steinberg A, et al. Investiga-tion of the syngas flame characteristics at elevated pressures using optical and laser diagnostic methods[J]. Flow, Turbulence and Combustion, 2011, 89(2): 275-294. |
| [25] |
闫博, 李猛, 陈力, 等. 基于滤波瑞利散射技术的带压燃烧场温度测量实验研究[J]. 实验流体力学, 2019, 33(4): 27-32. Yan B, Li M, Chen L, et al. Experimental study on temperature measurement of high pressure combustion based on filtered Rayleigh scattering technology[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(4): 27-32. (in Chinese) |
| [26] |
Wang S X, Wang Z H, Han X L, et al. Experimental and numerical study of the effect of elevated pressure on laminar burning velocity of lean H2/CO/O2/diluents flames[J]. Fuel, 2020, 273: 117753. |
| [27] |
Yan B, Chen L, Li M, et al. Quantitative temperature imaging at elevated pressures and in a confined space with CH4/air laminar flames by filtered Rayleigh scatte-ring[J]. Chinese Physics B, 2020, 29(2): 024701. |
| [28] |
Luque J, Crosley D. LIFBASE: Database and spectral simulation[R]. SRI International Report MP 99-009.1, 1999.
|
| [29] |
Goodwin D G, Speth R L, Moffat H K, et al. Cantera: An object-oriented software toolkit for chemical kinetics, thermodynamics, and transport processes[EB/OL]. https://www.cantera.org, Version 2.4.0, 2018.
|
| [30] |
Smith G, Golden D M, Frenklach M, et al. GRI-Mech 3.0[EB/OL]. https://www.cantera.org.
|
| [31] |
Rahmann U, Bulter A, Lenhard U, et al. LASKIN-A Simulation Program for Time-Resolved LIF-Spectra[EB/OL]. University of Bielefeld, Faculty of Chemistry, Physical Chemistry I, http://pc1.uni-bielefeld.de/?laskin.
|
| [32] |
Chen S, Su T, Zheng Y B, et al. Numerical analysis of quantitative measurement of hydroxyl radical concentra-tion using laser-induced fluorescence in flame[J]. Chi-nese Physics B, 2016, 25(6): 060703. |