邓 凯，李 华，杨臧健，钟英杰

（浙江工业大学能源与动力工程研究所脉动技术工程研究中心，杭州 310014）

基于红外和纹影技术研究声作用下火焰锋面运动及热量演化

邓 凯＊，李 华，杨臧健，钟英杰

（浙江工业大学能源与动力工程研究所脉动技术工程研究中心，杭州 310014）

以Rijke型驻波声场为声环境，以本生型层流部分预混火焰为实验对象，研究了火焰在声场作用下的锋面运动及热量演化行为。通过火焰红外图像表征燃烧反应产生的高温产物，借助火焰纹影图像描述径向温度梯度分布，并以红外－纹影信息一体化的图像说明了高温产物与径向温度梯度的关系。结果表明，声场作用下的部分预混火焰，其高温产物附着在锋面上，表现为层分布的对称型尖头结构，下游区域较上游区域幅面更宽；轴向热量演化过程以尖头状分布方式在声周期内经历2次升降过程，向火焰下游方向输运，径向热量传递过程所形成的温度梯度轮廓表现为“手指状”，由高温产物的尖头结构及运动决定。

火焰动力学；热量演化；部分预混火焰；驻波声场；燃烧不稳定

0 引 言

研究声场与火焰的作用有着重要的工程价值和理论意义。一方面，声与火焰的耦合，将有效地组织燃烧过程，不仅强化了燃烧和传热，而且有效降低污染物的排放量，脉动燃烧技术是其典型代表；另一方面，燃烧过程中放热波动与压力脉动满足一定相位关系（瑞利准则），将诱发热声振荡，也称为燃烧不稳定，是航空发动机、燃气轮机等重要动力装置设计与运行中所遇到的技术难题，目前为止，一直没有得到有效的解决。此外，基于流动与燃烧的关系，火焰在声场速度脉动作用下的响应行为，如放热、火焰动力学等，对于认识湍流火焰燃烧特性具有重要意义。

学术界当前主要关注声场作用下贫燃预混火焰的响应行为，代表性的思想是构建描述速度脉动与放热量变化之间关系的火焰传递函数。Durox等［1］测量了不同几何形状的火焰在声场下的火焰传递函数，根据传递函数增益和相位差2方面的数据，比较了形状不同带来的非线性响应差别，分析了几何形状对火焰动力学的影响。Balachandran等［2］实验研究了速度脉动作用下湍流预混火焰的非线性响应行为，重点关注贫燃预混火焰在大振幅声场作用下的饱和机理。Lieuwen［3］以本生型锥焰为建模对象，通过求解恒定火焰传播速度的锋面动力学方程（G－方程），研究了火焰结构特征参量对声场作用下火焰动力学特性的影响，分析了火焰传递函数从线性响应发展为非线性响应的条件。

这些工作存在1个共同的特点：研究对象是化学当量比为1左右的预混火焰，其原因可能有2点，一是工程中燃气轮机燃烧室热声振荡的条件是贫燃预混燃烧，二是理论上描述火焰锋面动力学特性的G－方程是基于一维火焰的假设。值得注意的是，作为燃烧方式中同样重要的部分预混火焰（Partially premixed flame），其在声场下的响应行为则很少被人关注。最近，Kim等［4］以燃气轮机燃烧室内的部分预混旋流燃烧器火焰为对象，研究了火焰在受到声场机理下的响应行为，作者认为声场作用下部分预混火焰的响应行为目前还没有被很好地认识，原因是其间的作用关系更为复杂，涉及化学、火焰动力学等方面的难点更为突出。

本文以Rijke型驻波声场为声环境，以本生型层流部分预混火焰为关注对象，借助多种测试手段，通过分析声波振荡周期内火焰的时间序列红外图像和纹影图像，揭示火焰锋面的运动过程及由此带来的热量演化行为。

1 实验系统与方法

1.1 声场发生系统

激励火焰的声环境由Rijke管实现，管内要实现激励的声场需满足以下条件：（1）Rijke管两端开口，保证管内存在一定的平均流速；（2）距管下端1／4管长处存在热源，且热源需形成一定的温度梯度；（3）保证Rijke管为良好的声保持系统，即管侧壁需完全密封，不发生声泄漏；（4）Rijke管长径比大于10。

综合考虑上述条件，设计Rijke管为两端开口的直管，管径为150mm，管长为1.6m。如图1所示，在管内距离下端口400mm（管长的1／4）处设置一通透性较好的近似平面热源，在热源加热过程中，热对流产生的管内气体流动与热源放热过程互相作用，使声压与放热满足一定的相位关系，最终激发管内的热声振荡。这种声环境是典型的一维驻波声场［5］，图1右侧给出了压力脉动量和速度脉动量的管内分布特征，管内中点位置压力振幅最大，速度振幅最小，管口两端则与此相反，在距管下端1／4管长位置，该点的压力振幅和速度振幅的乘积最大，最易于激发热声振荡。实验中的Rijke管采用内径150mm、壁厚4mm的碳钢管，热源由螺旋状电炉丝盘绕而成。为确保管内的流动通畅，热源支架由5根宽度为10mm的硅条按72°等角度均布排列。火焰受激励位置选择压力振幅和速度振幅相对较大的上1／4位置。声压特征通过在Rijke管侧壁上设置压力传感器提取，位置与火焰基部（燃烧器出口）对应。

理论上n阶的声压和速度波动表示为：

式中：An为压力振幅；kn＝ωn／c为波数；ωn＝nπc／L为n阶声场的频率；ρ0为气体密度；c是声速；L是管长。由上述公式可知，管内声场频率由本地声速和管长决定，振幅由热源的输出功率决定。实际测得的声场参数为：火焰位置的声压级（SPL）为144dB，管内声频率为111Hz，对应周期为9ms。

图1 火焰振荡声环境示意图Fig.1 Acoustic field excited by Rijke tube

1.2 燃烧系统

燃烧系统的实验装置如图2所示，甲烷和空气分别由各自供气系统送入混合腔4，流量用热式流量计测量，预混后的气体经混气管2进入烧嘴3，烧嘴由内径13mm、厚度2.5mm的不锈钢管加工而成，在没有声场作用的情况下火焰为本生型。Rijke管上方管段设置了一对平行放置的石英视窗，尺寸为90mm× 400mm，供纹影法和红外摄像观测火焰。火焰的特征参量为甲烷流量1000m L／min，一次空气当量比为1.57，燃烧器内气流Re数为790，这种火焰在没有声场作用的条件下表现为双火焰（Double flame）结构。

图2 燃烧系统示意图Fig.2 Experimental system

1.3 纹影和红外拍摄系统

纹影系统的布置如图3所示，光源处发出的光线，经过多组光学元件后，以平行光穿过Rijke管内的火焰区域，而后经一组球面反射镜汇聚于主反射镜（二）焦点位置的光刀处，纹影图像由光刀后方的相机记录，刀口布置的方向与狭缝方向设置相同。

拍摄火焰红外图像的CCD光谱如图4所示，红外CCD的平均响应主要集中在近红外波段。拍摄中对整体火焰结构进行拍摄，相机至于视窗前方一定位置（物距约22cm），另一视窗的内侧加设黑色云母板，目的是为拍摄过程提供慢反射背景，避免由于石英玻璃表面光线发生镜面反射对图像记录构成影响。

图3 纹影摄像系统示意图Fig.3 Schlieren system

图4 红外相机CCD光谱响应曲线Fig.4 Spectral response curve of infrared CCD

1.4 相位锁定同步拍摄系统

通过压力相位锁定的记录方式，对压力信号和红外图像进行同步触发，其系统如图5所示。

图5 压力信号与红外图像同步触发试验系统Fig.5 Synchronous system for pressure signal and infrared image

信号发生器1发出连续的TTL电平触发信号，经过信号开关2后被分为2路，一路与红外相机3的触发接口连接，另一路与数据采集器4的触发通道连接，压力传感器6与数据采集器4的其它采集通道连接。相机和数据采集器设置为单次触发模式，对应的作业方式为：由TTL电平信号触发后，仪器以自身的采样频率连续采集信号。具体的试验过程为：断开触发信号开关，分别激活红外相机和数据采集器的采集状态，当二者均为等待触发信号的状态时，闭合信号开关，TTL信号分别触发红外相机和数据采集器，实现了红外时序图像和声场压力信号的同步连续采集。为确保触发时刻足够短，红外相机和数据采集器均设置为上升沿（Rising）触发，信号发生器TTL电平信号的脉冲频率为8k Hz，对应的信号周期远小于声周期。此外，红外相机的曝光时间和快门等待时间分别为0.1和10ms，数据采集器的采样频率与相机的快门时间相匹配，设置为10 k Hz。

2 结果与讨论

2.1 红外和纹影技术的优势

火焰的红外图像用于描述声周期下火焰锋面的运动特性和热量演化，具有以下优势：

（1）在红外成像的物理过程上，火焰自发光红外光谱波段的光辐射具有较高的能量。火焰红外图像是由燃烧区域能够辐射红外光线的物质在CCD上成像而得到的，对于甲烷－空气的燃烧反应而言，二氧化碳和水是主要的燃烧产物，二者主要集中在红外波段［67］，即便是缺氧条件下形成碳烟，其辐射强度最大值也处于红外区［8］，因此，形成甲烷火焰红外图像的光线具有较大的辐射强度，图像中的平均灰度值更高，火焰区域与背景区域的对比度更大。

（2）甲烷火焰的红外图像，在物理意义上，描述了燃烧反应后产生高温产物的分布及运动。这些高温产物包围在燃烧反应区域附近，近似地表征了反应区域的运动过程。顾璠等［9］利用湍流火焰的红外图像表征大涡逆序结构的做法可说明这一点。

图6对可见图像和红外图像进行了对比，对火焰整体分别采用单反和红外相机进行拍摄。其中，可见图像通过单反相机获得，快门时间为0.25ms，红外图像的快门时间为0.1ms。外轮廓上，二者都具有对称的尖头（Cusp）结构，数量和空间分布几乎相同，但红外图像比单反图像中的尖头结构尺寸更大，对应总体的火焰幅面更宽。因此，对于描述锋面运动过程中尖头结构的发展与变化而言，红外图像更能够突显特征结构的发展。

图6 火焰可见图像与红外图像对比Fig.6 Comparison of visible and infrared flame images

火焰的纹影图像则可用于描述径向的热量传递，根据Dunn－Rankin等［10］的观点，纹影信息描述了燃烧区域温度梯度较大的位置，那么火焰纹影轮廓将反映火焰前锋附近的温度梯度，由此可推导对应的热量传递。最终，火焰红外和纹影2类图像的综合信息将有效揭示甲烷火焰的热量演化过程及特征。

2.2 部分预混火焰的锋面结构

图7给出了采用相位锁定方法拍摄的不同时刻的红外火焰图像，以此描述声场作用下部分预混火焰的锋面运动过程，其中红外相机的曝光时间为0.1ms。靠近燃烧器出口附近的火焰区域为火焰的上游，远离燃烧器出口的火焰顶部区域处为火焰的下游，对不同高度的尖头结构进行编号，在45°相位图像中由上游至下游的4层尖头依次编号为1～4号，1号和2号定义为上游，3号和4号定义为下游，以此描述锋面运动过程和结构特征。

由图7可知，对于部分预混火焰锋面而言，各个相位下的火焰具有相似的空间结构特征，尖头结构仍是其在声激励下的典型动力学特征。火焰结构呈现层状分布，锋面上尖头的数量相对固定，连续的火焰区域存在约4对；火焰下游的幅面更宽，上游的幅面相对较窄，同时距离下游火焰顶部一定距离处均存在弧线状孤立区域。

对预混火焰的研究表明［3］，预混锋面在声场作用下存在2种运动行为，包括沿锋面法向的传播运动和沿锋面切向的对流运动。对流波是2种运动形式产生的主要锋面几何特征，对应着预混锋面上轴对称的褶皱和尖头结构。对流波的波长决定锋面切向相邻2个尖头结构之间的距离，由锋面对流运动形成；对流波的幅值决定尖头结构的明显程度，由锋面传播运动形成。

而根据扩散锋面的研究结果，声激励的作用下，锋面运动由单方向对流转为往复式的变速运动。Hardalupas Y等［11］提出的唯象模型（Phenomenological model）描述，这种往复式的变速运动将在扩散锋面的外围区域形成一系列的涡街，相邻2个涡团的卷吸作用产生了锋面上的尖头结构。因此，尖头结构的明显程度应主要与涡团的运动行为有关。

综上所述，对于部分预混火焰而言，其在声作用下的结构应为对流波和涡团2种行为的综合结果。根据对流波的运动特征，声周期的不同时刻，锋面形成脉动位移最大值的位置不同，在火焰的下游和上游区域分别形成火焰基部和顶部的变形。因此，将火焰锋面上有尖头出现的区域进行划分，4个区域的出现描述了沿锋面法向的传播运动。

2.3 火焰锋面的运动行为

如图7所示，上半周期：前1／4周期内（0°～90°相位），上游的1号尖头包围燃烧器外壁的面积逐渐增大，红外轮廓的幅面增加，最低点略有下降；2号尖头的尖角上翘，在空间高度略有降低；3号尖头和4号尖头之间连接部分不断减小，4号尖头在空间上绝对高度保持不变，3号尖头处于下降过程，二者之间呈现了明显的分离特征，从整体上看，火焰在该时间段内呈现下降的运动形式。随后（90°～180°相位），90°相位下的4号尖头从火焰本体上脱离，形成了空间上孤立的红外区域，如180°相位的图像所示，这一过程描述了高温产物随火焰锋面运动最终脱离火焰本体的过程；同时，火焰基部红外轮廓的幅面逐渐减小，最低点不断上升，包围燃烧器外壁的特征消失，基部轮廓的宽度在135°后几乎接近燃烧器内径，尖头结构的尖角收缩，整个火焰呈拉伸状，运动形式表现为上升。

下半周期：锋面由向上运动转为向下运动，向上拉伸过程中逐渐减小的尖头尖角在下降过程中不断扩大，在火焰下游区域，顶部逐渐展开，幅面变宽，表现为压缩的锋面变形特征。此后（270°～360°相位），火焰锋面再次上升，呈现出90°～180°相位时间段内的锋面运动行为。

图7 火焰时间序列红外图像Fig.7 Infrared images at different phases

2.4 火焰的热量演化

整个周期内火焰区域的高温产物经历了2次下降和上升的过程，同时扩散锋面上尖头结构周期性地发展。结合锋面的运动行为和尖头结构的发展过程，认为声波作用下部分预混火焰沿轴向的热量演化路径为：火焰基部处的反应热，以高温产物为载体，由上游至下游运动，沿程汇聚不同位置（各层各处的尖头结构）产生的反应热，最终从火焰顶部脱离反应区域，以较大的表面积向下游的气体传递。

图8为拍摄的火焰时间序列纹影图像，使用的纹影系统光刀沿水平方向移动，因此图中亮暗区域及其对比度反映的是火焰区域沿水平方向的温度梯度，目的是说明径向的热量演化。总体结构上，包括内锥和外轮廓2部分，外轮廓的空间分布方式与红外图像相似，以一层一层的形式存在，但红外图像的尖头结构在纹影图像上表现为“手指状”结构。周期运动过程中，纹影图像的各层结构也存在下降和上升2个阶段，下面就上述存在的异同及反映的实际物理现象进行讨论。

由于红外图像表征的是高温产物的空间分布，而纹影图像描述的是水平方向温度梯度的分布，从热量与温度梯度的关系考虑，尖头结构作为“热源”，所形成的径向温度梯度轮廓线表现“手指状”的弧线是合理的，而“热源”一层一层分布也决定了温度梯度的分布形式。为验证上述分析的正确性，拍摄了纹影和红外信息一体化的图像。具体的做法是，以红外相机作为纹影系统的记录装置，将图3中的狭缝2尽量关小，使穿过火焰区域平行光的亮度与火焰自身发出的红外光线亮度接近，同时设置红外相机的增益程度，提高CCD的感光度（ISO），获得如图9所示的图像。图9中，明亮的区域与红外信息对应，依然表现为尖头状的特征，亮暗相间的区域为纹影信息，外轮廓线为“手指状”。“手指状”的纹影轮廓与“尖状体”位于同一高度，二者在“分层”和“配合”2个方面都比较吻合。

图8 火焰时间序列纹影图像Fig.8 Schlieren images at different phases

图9 红外与纹影信息一体化图像Fig.9 Schlieren and infared integrated images

由此可知，声波作用下的部分预混火焰区域中，由于高温产物与周围二次空气之间存在温度差和流速差，在火焰锋面的运动过程中存在径向的热量传递，结构上分层的尖头结构形成分层的“手指状”温度梯度轮廓，升降式的运动方式使不同时刻的径向热量传递过程差别较大。

3 结 论

（1）声波作用下的部分预混火焰，在结构上表现为多层对称的尖头结构，尖头结构与对流波和涡团运动有关；锋面在声波振荡的周期内经历2次上升和下降的过程，热量演化由这些结构发展过程和锋面运动方式决定。

（2）轴向热量传递：火焰基部处的反应热，以高温产物为载体，由上游至下游运动，沿程汇聚不同位置（各层尖头结构）产生的反应热，最终从火焰顶部脱离反应区域，以较大的表面积向下游的气体传递。

（3）径向热量传递：形成分层的“手指状”温度梯度轮廓，发生径向传热过程的轴向位置随相位的改变周期性变化。

［1］Durox D，Schuller T，Noiray N，et al.Experimental analysis of nonlinear flame transfer functions for different flame geometries［J］.Proceedings of the Combustion Institute，2009，32（1）：1391－1398.

［2］Balachandran R，Ayoola B O，Kaminski C F，et al.Experimental investigation of the nonlinear response of turbulent premixed flames to imposed inlet velocity oscillations［J］.Combustion and Flame，2005，143（1－2）：37－55.

［3］Lieuwen T.Nonlinear kinematic response of premixed flames to harmonic velocity disturbances［J］.Proceedings of the Combustion Institute，2005，30（2）：1725－1732.

［4］Kim K T，Lee J G，Quay B D，et al.Response of partially premixed flames to acoustic velocity and equivalence ratio perturbations［J］.Combustion and Flame，2010，157（9）：1731－1744.

［5］马大猷.黎开管振荡的简单理论［J］.声学学报，2001，26（4）：289－294.Ma D Y.Simple theory of Rijke tube oscillation［J］.Acta Acustica，2001，26（4）：289－294.

［6］刘晅亚，陆守香，朱迎春.水雾作用下甲烷／空气预混火焰的光谱特性［J］.燃烧科学与技术，2008，14（1）：44－49.Liu X Y，Lu S X，Zhu Y C.Study of methane／air premixed flame emission spectrum under the influence of water mist［J］.Combustion Science and Technology，2008，14（1）：44－49.

［7］王喜世，伍小平，秦进.用红外热成像方法测量火焰温度的实验研究［J］.激光与红外，2001，31（3）：169－172.Wang X S，Wu X P，Qin J.Experimental study on flame temperature measurement by infrared thermography［J］.Laser and Infrared，2001，31（3）：169－172.

［8］Turns S R.燃烧学导论：概念与应用［M］.姚强，李水清，王宇译.北京：清华大学出版社，2009.

［9］顾璠，杨宏旻，刘勇，等.湍流射流与扩散火焰大涡拟序结构的波动特性研究［J］.工程热物理学报，1999，20（6）：773－776.Gu F，Yang H M，Liu Y，et al.Wave characteristics of large eddy coherent structures of turbulent jets and diffusion flames［J］.Journal of Engineering Thermophysics，1999，20（6）：773－776.

［10］Dunn－Rankin D，Weinberg F.Location of the schlieren image in premixed flames：axially symmetrical refractive index fields［J］.Combustion and Flame，1998，113（3）：303－311.

［11］Hardalupas Y，Selbach A.Imposed oscillations and nonpremixed flames［J］.Progress in Energy and Combustion Science，2002，28（1）：75－104.

Investigation of heat transfer and flame dynamics under acoustic excitation based on infrared and shadow method

Deng Kai＊，Li Hua，Yang Zangjian，Zhong Yingjie
（Institute of Energy and Power Engineering，Engineering Research Center of Pulse Technology，Zhejiang University of Technology，Hangzhou 310014，China）

The heat transfer of partially premixed laminar flame under Rijke acoustic standing wave was investigated.The flame dynamics and hot products could be captured by the infrared images and the radial temperature gradient distribution could be shown by schlieren imaging.The relationship between the hot products and the radial temperature is characterized with infraredschlieren information integration.The results show that the hot products attached on the flame front under acoustic excitation present a symmetric structure of tip layer distribution.And the downstream region is wider than the upstream region.In one sound cycle the axial heat evolution shows a tip shape transporting heat to the downstream of the flame through two lifting processes.The temperature gradient induced by the radial heat transfer takes a finger shape whose evolvement is decided by the development of the cusps.

flame dynamics；heat transfer；partially premixed flame；acoustic standing wave；combustion instability

TK31

A

（编辑：李金勇）

1672－9897（2016）06－0026－06

10.11729／syltlx20160069

2016－04－26；

2016－08－19

国家自然科学基金项目（51106139）

＊通信作者E－mail：dkai＠zjut.edu.cn

Deng K，Li H，Yang Z J，et al.Investigation of heat transfer and flame dynamics under acoustic excitation based on infrared and shadow method.Journal of Experiments in Fluid Mechanics，2016，30（6）：26－31.邓 凯，李 华，杨臧健，等.基于红外和纹影技术研究声作用下火焰锋面运动及热量演化.实验流体力学，2016，30（6）：26－31.

邓 凯（1984－），女，湖北天门人，博士，讲师。研究方向：火焰动力学／燃烧污染物／燃烧不稳定性。通信地址：浙江省杭州市潮王路18号浙江工业大学能源与动力工程研究所（310014）。E－mail：dkai＠zjut.edu.cn