閆 博, 李 猛, 陳 力, 陳 爽, 吳運(yùn)剛, 楊富榮, 母金河

(中國(guó)空氣動(dòng)力研究與發(fā)展中心, 四川 綿陽 621000)

0 引 言

發(fā)動(dòng)機(jī)燃燒室的壓強(qiáng)通常在0.2MPa以上，溫度約為1000～3000K，且伴隨有激波邊界層干擾、混合邊界層等現(xiàn)象。同時(shí)發(fā)動(dòng)機(jī)燃燒試驗(yàn)條件極為苛刻，包括受限空間試驗(yàn)布局，受限空間引入的強(qiáng)雜散光干擾，超聲速燃燒流場(chǎng)引起的高度湍流特性，燃燒室溫

度、壓強(qiáng)分布的不均勻性等。這些因素導(dǎo)致發(fā)動(dòng)機(jī)燃燒場(chǎng)的溫度測(cè)量非常困難，亟需響應(yīng)速度快、不干擾流場(chǎng)、高時(shí)空分辨率、高靈敏度和可消除雜散光干擾的燃燒場(chǎng)溫度測(cè)量技術(shù)[2- 3]。

目前，能夠?qū)崿F(xiàn)燃燒場(chǎng)二維溫度測(cè)量的技術(shù)主要有平面激光誘導(dǎo)熒光技術(shù)和瑞利散射技術(shù)。平面激光誘導(dǎo)熒光技術(shù)通常采用雙波長(zhǎng)激發(fā)來獲取燃燒場(chǎng)的溫度信息[4- 5]，但在帶壓燃燒場(chǎng)環(huán)境下，分子間碰撞加劇，受碰撞淬滅的影響，熒光信號(hào)急劇下降，這給溫度測(cè)量帶來了極大的挑戰(zhàn)。瑞利散射技術(shù)利用激光照射燃燒場(chǎng)獲取瑞利散射信號(hào)來測(cè)量流場(chǎng)的溫度，與平面激光誘導(dǎo)熒光技術(shù)不同，瑞利散射信號(hào)強(qiáng)度與氣體的數(shù)密度成正比，在帶壓情況下，壓強(qiáng)越高，瑞利散射信號(hào)越強(qiáng)，越容易實(shí)現(xiàn)高精度溫度測(cè)量。然而，復(fù)雜的燃燒場(chǎng)測(cè)量環(huán)境會(huì)引入很強(qiáng)的雜散光(主要來自于米散射和背景散射)，這會(huì)極大影響瑞利散射信號(hào)的信噪比，降低溫度測(cè)量的精度[6- 7]。

為消除雜散光干擾，通常采用窄帶寬濾波片對(duì)信號(hào)光進(jìn)行濾波處理，但由于瑞利散射光頻率與雜散光頻率基本相同，導(dǎo)致普通的窄帶寬濾波片無法濾除大量的雜散光干擾。為此，20世紀(jì)90年代，Miles和Lempert[8]在傳統(tǒng)瑞利散射技術(shù)的基礎(chǔ)上，發(fā)展了一種用于流場(chǎng)診斷的濾波瑞利散射技術(shù)(FRS)，通過在待測(cè)高溫流場(chǎng)和探測(cè)器之間放置碘分子濾波池，實(shí)現(xiàn)了對(duì)雜散光的有效抑制，從而極大提高了瑞利散射信號(hào)的信噪比。因此，濾波瑞利散射技術(shù)一經(jīng)提出，便受到國(guó)內(nèi)外眾多研究者的廣泛關(guān)注。 Forkey等[9]校正了碘分子濾波池的吸收模型；Hoffman等[10]在甲烷/空氣火焰(含有輕微碳煙)中開展了二維溫度測(cè)量；德國(guó)DLR研究所的Doll[11- 12]和Schroll等[13]針對(duì)常壓燃燒場(chǎng)和發(fā)動(dòng)機(jī)燃燒場(chǎng)進(jìn)行了二維溫度定量測(cè)量研究；西北核技術(shù)研究所的王晟等[14]和中國(guó)空氣動(dòng)力研究與發(fā)展中心的鄭堯邦等[15]也開展了常壓燃燒場(chǎng)溫度和密度測(cè)量的初步研究。綜上所述，國(guó)內(nèi)外關(guān)于FRS技術(shù)的研究主要是針對(duì)常壓燃燒場(chǎng)溫度的測(cè)量，關(guān)于帶壓燃燒場(chǎng)的溫度測(cè)量工作較少。

因此，本文將濾波瑞利散射技術(shù)應(yīng)用于高壓火焰溫度測(cè)量，針對(duì)高壓平面層流火焰爐產(chǎn)生的帶壓燃燒場(chǎng)試驗(yàn)環(huán)境(1000～2000K，0.1～0.5MPa)，開展了甲烷/空氣預(yù)混火焰的二維溫度分布測(cè)量實(shí)驗(yàn)研究。

1 基本原理

當(dāng)線寬很窄的一束激光投射到燃燒流場(chǎng)時(shí)，將伴隨激光散射現(xiàn)象。通常，散射信號(hào)既包含氣體分子引起的瑞利散射信號(hào)，也包含大量雜散光(來自米散射和背景散射)。如圖1所示，燃燒流場(chǎng)氣體分子的瑞利散射，因其內(nèi)部大量分子無規(guī)則的熱運(yùn)動(dòng)，導(dǎo)致其譜線較入射激光有著明顯的展寬。同時(shí)，對(duì)于不同的氣體溫度和壓強(qiáng)，瑞利散射光譜線型也不同。但是對(duì)于雜散光，其頻率和線寬均與入射激光相同，且信號(hào)較強(qiáng)，在燃燒場(chǎng)溫度測(cè)量過程中屬于干擾噪聲。為減小雜散光的干擾，提高瑞利散射信號(hào)的信噪比，可通過在燃燒流場(chǎng)與探測(cè)器ICCD之間放置具有窄吸收峰的分子濾波池，并通過波長(zhǎng)調(diào)諧技術(shù)，將窄線寬激光光源的輸出波長(zhǎng)調(diào)至濾波分子的吸收峰中心(如圖1所示)，雜散光因與入射激光具有相同頻率和線寬，并正好位于濾波分子的強(qiáng)吸收峰中心而被過濾掉；而燃燒場(chǎng)的瑞利散射光將部分透過吸收峰被ICCD所探測(cè)。不同的燃燒狀態(tài)，所引起的瑞利散射譜線展寬不同，因此通過濾波池后，ICCD相機(jī)探測(cè)到的光強(qiáng)也不一樣，根據(jù)濾波瑞利散射信號(hào)的強(qiáng)弱，可以反解出燃燒場(chǎng)的溫度信息。

圖1 碘分子對(duì)激光散射信號(hào)吸收示意圖

Fig.1SchematicofiodinetransmissionprofileconvolutedwiththemolecularRayleighscatteringspectralprofile

考慮激光能量以及實(shí)際燃燒場(chǎng)的其他物理量，ICCD探測(cè)到的濾波瑞利散射信號(hào)可以表示為：

常規(guī)超聲可以顯示病灶的大小、位置、數(shù)目、形態(tài)及毗鄰情況，清晰顯示管系內(nèi)有無癌栓。但是，受到消融區(qū)域內(nèi)氣泡形成的強(qiáng)回聲的干擾，常規(guī)超聲的臨床作用大打折扣，消融形成的不規(guī)則高回聲區(qū)還會(huì)遮掩目標(biāo)腫瘤，難以形成令人滿意的消融后邊界。常規(guī)超聲還易受肋骨、肥胖、肝硬化背景等影響，使圖像顯示不清晰。

(1)

式中，C為校準(zhǔn)常數(shù)，I0為入射激光強(qiáng)度，N為燃燒場(chǎng)的數(shù)密度，ω為激光的頻率，χk和Mk為燃燒場(chǎng)中第k種分子的摩爾組分和質(zhì)量， (?σ/?Ω)為瑞利散射截面， Rk為瑞利散射光譜線型，τ(ω)為濾波分子吸收峰透過率曲線。式(1)中的數(shù)密度N可表示為N=p/KT，其中p為燃燒場(chǎng)壓強(qiáng)，K為波爾茲曼常數(shù)，T為燃燒場(chǎng)溫度。為消除校準(zhǔn)常數(shù)C以及入射激光強(qiáng)度I0等參數(shù)的影響，引入室溫T0下空氣的濾波瑞利散射信號(hào)S0為參考，并定義FRS相對(duì)強(qiáng)度S*為：

(2)

由式(2)，根據(jù)燃燒條件及濾波分子透過率曲線的條件，可建立FRS相對(duì)強(qiáng)度與溫度的對(duì)應(yīng)關(guān)系S*-T曲線，將實(shí)驗(yàn)測(cè)得的FRS相對(duì)強(qiáng)度S*與理論計(jì)算的S*-T曲線比對(duì)，即可獲得流場(chǎng)的溫度信息。

2 實(shí)驗(yàn)裝置簡(jiǎn)介

2.1 裝置光路簡(jiǎn)介

濾波瑞利散射溫度測(cè)量裝置主要由激光器、片光系統(tǒng)、能量波長(zhǎng)監(jiān)測(cè)系統(tǒng)、分子濾波池、同步時(shí)序控制器及ICCD相機(jī)組成。如圖2所示，激光器為種子注入Nd∶YAG激光器，提供在532nm波長(zhǎng)附近可調(diào)的窄線寬脈沖光源，激光器輸出能量為400mJ，激光脈寬約10ns，波長(zhǎng)選擇為532.218 5nm。激光器輸出的激光經(jīng)分束比約為24∶1的分束鏡分成2束，能量低的一束進(jìn)入能量波長(zhǎng)監(jiān)測(cè)系統(tǒng)，其中光電倍增管PMT(英國(guó)Sens- Tech公司生產(chǎn)，P30A- 03)用于對(duì)激光輸出能量監(jiān)測(cè)，F(xiàn)- P干涉儀用于對(duì)激光波長(zhǎng)漂移進(jìn)行監(jiān)測(cè)，實(shí)驗(yàn)時(shí)，通過對(duì)激光能量和波長(zhǎng)的監(jiān)測(cè)，可對(duì)探測(cè)到的濾波瑞利散射圖像進(jìn)行必要的修正；分束后，能量高的一束激光經(jīng)過由柱面鏡組組成的片光系統(tǒng)后，形成一片激光通過高壓火焰爐觀察窗照射高壓燃燒火焰。

圖2 分子濾波瑞利散射溫度測(cè)量光路示意圖

Fig.2SchematicofopticalarrangementforFRStemperaturemeasurements

高壓火焰爐如圖3(a)所示。該高壓火焰爐可以實(shí)現(xiàn)壓強(qiáng)0.1～1.0MPa之間連續(xù)可調(diào)，且壓強(qiáng)控制精度為0.01MPa。同時(shí)，火焰爐為標(biāo)準(zhǔn)McKenna平面火焰爐(如圖3(b)所示)，爐面尺寸為25mm，可通過調(diào)節(jié)爐內(nèi)壓力和CH4/空氣比例，產(chǎn)生穩(wěn)定的層流火焰。在火焰中心區(qū)域懸置熱電偶用于火焰溫度點(diǎn)測(cè)量，圖3(b)中的標(biāo)記點(diǎn)A、B和C為實(shí)驗(yàn)中所用到的熱電偶探測(cè)點(diǎn)，距離爐面的高度分別為hA=12.5mm，hB=17mm，hC=21.5mm。瑞利散射光收集成像裝置如圖4所示，火焰的散射光經(jīng)過一個(gè)碘分子濾波池(碘分子池的溫度控制在52.3℃)，實(shí)現(xiàn)對(duì)雜散光的過濾吸收后，被ICCD相機(jī)所探測(cè)。整個(gè)測(cè)量裝置的同步時(shí)序控制由一臺(tái)DG535實(shí)現(xiàn)。

(a)

(b)

圖3 (a)高溫高壓火焰爐實(shí)物圖；(b)McKenna火焰爐及其上方實(shí)驗(yàn)布局尺寸圖

Fig.3(a)Schematicofhightemperatureandpressurecombustionburner; (b)ThesizeofMcKennaburnerandtheexperimentlayout

圖4 瑞利散射光收集成像裝置實(shí)物圖

2.2 瑞利- 布里淵散射光譜理論計(jì)算

為建立FRS相對(duì)強(qiáng)度與溫度的對(duì)應(yīng)關(guān)系S*-T曲線，必須先獲得不同溫度和壓強(qiáng)條件下瑞利- 布里淵散射的光譜線型，進(jìn)而獲得不同溫度和壓強(qiáng)條件下瑞利散射信號(hào)經(jīng)過濾波池后的光強(qiáng)透過率。

對(duì)于氣體分子的瑞利- 布里淵散射光譜線型，國(guó)外學(xué)者已開展了大量的深入研究，建立了許多物理模型和計(jì)算方法。其中，Tenti等[16]建立的S6模型因精度較高而被廣泛使用，近年來Pan等[17]建立了S7模型，進(jìn)一步提高了計(jì)算精度，其計(jì)算結(jié)果也得到了大量實(shí)驗(yàn)驗(yàn)證。本文采用S7模型，計(jì)算了不同溫度、壓強(qiáng)條件下，氮?dú)夥肿拥娜鹄? 布里淵散射光譜分布，依據(jù)式(2)，建立了FRS相對(duì)強(qiáng)度與溫度的對(duì)應(yīng)關(guān)系S*-T曲線，如圖5所示。通過ICCD探測(cè)燃燒和未燃燒條件下的流場(chǎng)濾波瑞利散射信號(hào)，并依據(jù)S*-T曲線，即可得到燃燒火焰的二維溫度分布結(jié)果。

圖5 不同壓強(qiáng)下FRS相對(duì)強(qiáng)度與溫度的對(duì)應(yīng)關(guān)系

Fig.5NormalizedFRSsignalversusflowtemperaturecurveunderdifferentpressures

3 結(jié)果與分析

為驗(yàn)證濾波瑞利散射測(cè)溫裝置對(duì)米散射和背景雜散光的抑制能力，按照?qǐng)D2所示的光路，開展了濾波和未濾波條件下空氣散射流動(dòng)顯示對(duì)比實(shí)驗(yàn)。在未添加碘分子濾波池的條件下，得到空氣對(duì)入射片激光的散射圖像，如圖6(a)所示。從圖中可觀察到很強(qiáng)的背景雜散光，如熱電偶探針的反射光和平面火焰爐爐面的反射光，雜散光已經(jīng)將瑞利散射信號(hào)全部掩蓋，極大影響了探測(cè)到的瑞利散射信號(hào)的信噪比。在散射光收集光路中添加碘分子濾波池，并將入射激光波長(zhǎng)調(diào)諧至圖1所示的吸收峰中心，獲得濾波后的空氣瑞利散射圖像，如圖6(b)所示?？梢钥闯?，由平面火焰爐爐面和熱電偶探針引起的背景雜散光幾乎全部被過濾吸收掉，即系統(tǒng)有效實(shí)現(xiàn)了對(duì)雜散光的抑制，明顯提高了瑞利散射信號(hào)圖像的信噪比。

為進(jìn)一步說明設(shè)計(jì)的濾波瑞利散射測(cè)溫裝置對(duì)雜散光的抑制能力，采用圖2所示光路得到了未燃條件、不同腔壓下的瑞利散射信號(hào)值(圖6(b)中紅色方框內(nèi)的瑞利散射信號(hào)平均值)，如圖7所示。由式(1)可知，瑞利散射信號(hào)強(qiáng)度值S與氣體分子密度N成正比，而N=pV/RT，其中V為氣體腔室體積，p為腔壓，R為理想氣體常數(shù)，T為腔內(nèi)氣體溫度，故S∝p。由圖7中的線性擬合公式S=657 704.9p+1662.9可知，由雜散光干擾引起的截距項(xiàng)僅為1662.9，小于總信號(hào)強(qiáng)度的3%，說明設(shè)計(jì)的濾波瑞利散射測(cè)溫裝置具有較好的背景雜散光抑制能力。其次，利用濾波瑞利散射測(cè)溫裝置在高壓火焰爐上開展了預(yù)混燃燒火焰溫度測(cè)量實(shí)驗(yàn)。帶壓燃燒實(shí)驗(yàn)參數(shù)如表1所示，測(cè)量區(qū)域片激光的厚度約0.5mm，寬度約25mm，激光下邊緣距爐面約8mm(如圖3(b)所示)。

(a) 無濾波池條件下的散射圖像

(b) 有碘濾波池條件下的空氣瑞利散射圖像

Fig.6Scatteringimageswithoutandwiththeiodinefilteringcell

圖7 ICCD相機(jī)測(cè)得的瑞利散射信號(hào)強(qiáng)度S(20幅累加)與腔體壓強(qiáng)p的關(guān)系

Fig.7TherelationbetweentheRayleighscatteringsignalSandthechamberpressurep

表1 帶壓燃燒實(shí)驗(yàn)參數(shù)(壓強(qiáng)p，甲烷和空氣流量QCH4、QAir，當(dāng)量比φ及其預(yù)混燃?xì)饪偭魉賤premixed)Table 1 Parameters for combustion experiments (Chamber pressure p, flow rates for air through burner QAir and for fuel CH4QCH4, equivalence ratios φ, premixed gas velocity vpremixed)

圖8(a)顯示了瞬態(tài)火焰和50幅火焰圖像累加的溫度場(chǎng)分布平均結(jié)果。可以看出，隨著高壓火焰爐腔內(nèi)壓強(qiáng)的增加，火焰形狀變得尖銳，火焰截面變窄。這是因?yàn)楫?dāng)甲烷和空氣的預(yù)混氣體流量不變時(shí)，隨著腔壓的升高，單位時(shí)間內(nèi)通過的預(yù)混氣體體積減小(V=NRT/p)，預(yù)混氣體的出射速度v減小(v=V/Θ，V為單位時(shí)間內(nèi)通過的氣流體積，Θ為截面積)，進(jìn)而火

(a)

(b)

(c)

圖8 (a) 經(jīng)濾波瑞利散射技術(shù)得到的火焰溫度分布圖像；(b) 不同腔壓下，火焰爐上方1.5cm處的火焰溫度分布圖; (c)不同探測(cè)位置處，熱電偶和FRS (50幅平均)溫度測(cè)量結(jié)果

Fig.8(a)ImagesofcombustiontemperaturedistributionbytheFRStechnique; (b)Distributionofcombustiontemperatureattheheight(1.5cm)abovetheburner; (c)ComparisonofthetemperatureresultsbetweentheFRStechniqueandthethermocouple

焰形狀變得尖銳，火焰面變窄。同時(shí)，為了對(duì)比不同壓強(qiáng)下的濾波瑞利散射技術(shù)溫度測(cè)量的相對(duì)不確定度，圖8(b)對(duì)比了火焰爐上方1.5cm且平行爐面處的瞬態(tài)瑞利散射溫度分布曲線。由圖可知：0.11MPa時(shí)的溫度測(cè)量相對(duì)不確定度約在15%以內(nèi)，0.30MPa時(shí)約在10%以內(nèi)，而0.50MPa時(shí)約在7%以內(nèi)。這是因?yàn)殡S著腔壓的升高，瑞利散射信號(hào)強(qiáng)度增大，瑞利散射圖像的信噪比增強(qiáng)，造成溫度測(cè)量相對(duì)不確定度減小。最后，為了驗(yàn)證濾波瑞利散射技術(shù)的溫度測(cè)量精度，分別在A、B和C點(diǎn)(如圖3(b)和8(a)所示)進(jìn)行了熱電偶溫度點(diǎn)測(cè)量，并將熱電偶測(cè)量結(jié)果與濾波瑞利散射測(cè)溫裝置得到的結(jié)果(50幅平均)進(jìn)行對(duì)比，如圖8(c)所示。由圖可知，2種測(cè)量方法均可得到：火焰距離爐面的高度越高，溫度越低，這是因?yàn)榛鹧嬖谙蛏蟼鞑r(shí)會(huì)有熱損失，從而導(dǎo)致溫度變低；壓強(qiáng)越高，相同位置處的溫度越低，這是因?yàn)楫?dāng)預(yù)混燃?xì)?甲烷/空氣)的總量不變時(shí)，壓強(qiáng)越高，預(yù)混燃?xì)庀蛏系乃俣仍降停瑥亩鴮?dǎo)致層流火焰面整體向下偏移。此外，對(duì)于距離爐面較近的探測(cè)點(diǎn)A來說，濾波瑞利散射溫度測(cè)量結(jié)果和熱電偶測(cè)量結(jié)果偏差在10%以內(nèi)(不同壓強(qiáng))。而隨著探測(cè)點(diǎn)距離爐面高度的增加，2種測(cè)量方法之間的偏差逐漸增大，且腔壓越高，偏差越大，尤其是對(duì)于腔壓為0.50MPa中的探測(cè)點(diǎn)C來說，2種方法得到的溫度值最大偏差達(dá)到18%。這是因?yàn)殡S著探測(cè)點(diǎn)距離爐面的高度的增加，燃?xì)鈿饬鞯乃俣葴p小，受熱傳遞的影響，上方火焰出現(xiàn)擾動(dòng)(由圖9可以看出，不同時(shí)間的火焰形狀不同)，進(jìn)而導(dǎo)致2種測(cè)量方法得到的溫度偏差增大。

圖9 腔壓p=0.50MPa時(shí)，不同時(shí)間下的火焰瑞利散射信號(hào)分布圖像

Fig.9ImagesofcombustionRayleighscatteringsignaldistributionwiththechamberpressurep=0.50MPa,atdifferenttimes

4 結(jié) 論

本文基于濾波瑞利散射原理技術(shù)，以大功率窄線寬Nd∶YAG脈沖激光器為光源，以碘分子濾波池為背景光抑制器件，以ICCD相機(jī)為圖像采集設(shè)備，設(shè)計(jì)了一套濾波瑞利散射溫度測(cè)量裝置，并利用該裝置在高壓火焰爐上開展了不同壓強(qiáng)條件下甲烷/空氣預(yù)混火焰溫度場(chǎng)診斷實(shí)驗(yàn)。結(jié)果表明：

(1) 濾波瑞利散射測(cè)溫技術(shù)能有效抑制受限空間引入的背景雜散光干擾；

(2) 濾波瑞利散射測(cè)溫技術(shù)能夠獲得帶壓條件下的瞬態(tài)燃燒場(chǎng)溫度分布結(jié)果，且溫度測(cè)量的相對(duì)不確定度優(yōu)于15%；

(3) 通過與熱電偶溫度測(cè)量實(shí)驗(yàn)的結(jié)果進(jìn)行對(duì)比，兩者的偏差約為10%，驗(yàn)證了濾波瑞利散射技術(shù)具備帶壓燃燒場(chǎng)溫度的非接觸測(cè)量能力。

本文工作對(duì)實(shí)現(xiàn)受限空間、帶壓條件下的燃燒場(chǎng)溫度分布定量測(cè)量具有參考意義，下一步計(jì)劃將FRS技術(shù)應(yīng)用于內(nèi)燃機(jī)燃燒場(chǎng)和超燃沖壓燃燒場(chǎng)溫度分布定量測(cè)量實(shí)驗(yàn)。