吳軍，宋豐成，潘智祥，張琛萍，郭潤(rùn)夏，張曉瑜

（1.中國(guó)民航大學(xué) 航空工程學(xué)院，天津 300300；2.中國(guó)民航大學(xué) 電子信息與自動(dòng)化學(xué)院，天津 300300）

1 引言

航空發(fā)動(dòng)機(jī)的運(yùn)行狀態(tài)監(jiān)測(cè)在飛行安全領(lǐng)域發(fā)揮著重要作用［1-2］。在航空發(fā)動(dòng)機(jī)工作時(shí)，尾氣產(chǎn)生的溫度場(chǎng)中蘊(yùn)含著大量關(guān)于發(fā)動(dòng)機(jī)燃燒狀態(tài)的信息［3-4］。這些狀態(tài)信息都可以對(duì)發(fā)動(dòng)機(jī)運(yùn)行狀態(tài)控制提供理論依據(jù)與數(shù)據(jù)支持［5-6］。以航空發(fā)動(dòng)機(jī)的尾氣為例，根據(jù)火焰燃燒特性，溫度分布按照中心對(duì)稱并沿著空氣流動(dòng)方向逐漸降低形成三維溫度場(chǎng)，采取接觸式測(cè)溫會(huì)破壞流場(chǎng)結(jié)構(gòu)，影響溫度測(cè)量，現(xiàn)階段常用光學(xué)測(cè)量法通過(guò)測(cè)量空氣折射率分布實(shí)現(xiàn)非接觸式測(cè)溫［7］。激光吸收光譜測(cè)溫法是利用光在散射后發(fā)生能量衰減來(lái)對(duì)溫度測(cè)量［8-10］，常適用于需要搭建大量測(cè)量結(jié)構(gòu)的溫度場(chǎng)截面測(cè)量，在尾氣流場(chǎng)這樣溫度場(chǎng)測(cè)量中受到限制。

現(xiàn)代定量紋影測(cè)量技術(shù)由 Moreno-Hernández［11］提出，具有精度高，可視化效果好等優(yōu)點(diǎn)。紋影系統(tǒng)結(jié)構(gòu)近似為傅里葉光學(xué)處理器，使用刀口充當(dāng)了空間頻譜的半平面掩模，對(duì)通過(guò)待測(cè)流場(chǎng)的平行光束實(shí)現(xiàn)空間濾波，可以顯示出折射率場(chǎng)中的紋影細(xì)節(jié)［12］。這一方法受到透鏡尺寸的影響，一般測(cè)量范圍在15 cm 左右，多用于實(shí)驗(yàn)室內(nèi)測(cè)量。2013 年Martinez-Gonzalez［13］采用定標(biāo)紋影法測(cè)量了加熱金屬板引起的空氣自然對(duì)流中的溫度場(chǎng)，2017 年，Alvarez-Herrera［14］通過(guò)改良紋影技術(shù)完成了發(fā)動(dòng)機(jī)內(nèi)部溫度場(chǎng)重建。

反射式紋影技術(shù)，是20 世紀(jì)中期提出的一種將高速相機(jī)與反光屏幕和紋影技術(shù)相結(jié)合的簡(jiǎn)易紋影法。該方法早期被應(yīng)用在航空彈道學(xué)和直升機(jī)旋翼的測(cè)試分析中，現(xiàn)在與高速相機(jī)攝像、反光屏幕等相結(jié)合進(jìn)行紋影測(cè)量［15］。2017 年Gihyeon Min［16］在實(shí)驗(yàn)室條件 下用此方法對(duì) 蠟燭、沸騰熱水壺、加熱器和燃?xì)饣鹁嬉鸬目諝馔牧鬟M(jìn)行了成像。經(jīng)過(guò)亮度補(bǔ)償后，干涉精度對(duì)比原始紋影精度平均提高14.15 倍。2018 年，Martinez-Gonzalez［17］利用此方法設(shè)計(jì)了一種測(cè)量流體平面溫度場(chǎng)的方法，并測(cè)量了商業(yè)熱風(fēng)槍焊接時(shí)噴嘴發(fā)出的熱空氣平面溫度場(chǎng)?？梢钥闯?，與激光診斷技術(shù)［18］相結(jié)合的光學(xué)測(cè)量方法在溫度測(cè)量方面具有更高的工程應(yīng)用價(jià)值。結(jié)合二者優(yōu)點(diǎn)的紋影測(cè)溫方法具有非接觸測(cè)量、對(duì)測(cè)量結(jié)果響應(yīng)速度快和空間分辨率高的優(yōu)勢(shì)。由于傳統(tǒng)的背景紋影法的局限性是測(cè)量范圍基于相機(jī)的成像模型，不可避免存在圖像畸變，而邊緣位置圖像偏移量本身較小，會(huì)導(dǎo)致整體測(cè)量精度的下降。本文結(jié)合紋影成像模型和激光干涉測(cè)溫等原理，提出了利用干涉條紋光程差解算折射率分布的方法。該方法在保留背景紋影法測(cè)量設(shè)備簡(jiǎn)單的基礎(chǔ)上，極大地提高了測(cè)量的穩(wěn)定性與精確性。由于激光干涉條紋對(duì)折射率場(chǎng)變化敏感且有良好的相干性，因此產(chǎn)生的干涉條紋可以放大光線偏折角，進(jìn)而有效提高紋影法溫度測(cè)量精度。

2 干涉條紋紋影軸對(duì)稱溫度場(chǎng)測(cè)溫原理

2.1 干涉條紋的形成及變形原理

條紋背景紋影法，是利用光線在通過(guò)固定狹縫后會(huì)發(fā)生的干涉現(xiàn)象對(duì)固定場(chǎng)內(nèi)影響折射率分布因素的測(cè)量。

楊氏雙縫干涉實(shí)驗(yàn)的條紋產(chǎn)生原理如圖1 所示。相干光源S1和S2由單色點(diǎn)光源S垂直射到相鄰的兩個(gè)小孔上，此時(shí)兩個(gè)小孔就相當(dāng)于兩個(gè)同相位的次級(jí)光源，光束在xOy平面上的成像區(qū)域疊加會(huì)產(chǎn)生條紋圖樣。這種干涉條紋是由強(qiáng)度相同，相位不同的光干涉產(chǎn)生的。

圖1 兩個(gè)點(diǎn)光源產(chǎn)生干涉的幾何圖示Fig.1 Geometric diagram of interference from a set of coherent light sources

設(shè)定干涉圖樣存在xOy平面上，該平面與S1S2的垂直等分線CO正交，且x軸平行于S1S2。假設(shè)兩個(gè)小孔的間隔d相對(duì)于兩孔連線到觀察面的距離a很小，對(duì)于點(diǎn)成像平面上的P（x，y）而言，存在路徑關(guān)系：

得到實(shí)際路徑差：

根據(jù)溫度場(chǎng)的變化，r1與r2路徑上的折射率值發(fā)生改變，同一位置在溫度場(chǎng)存在前后的光程差主要受到溫度影響，從而改變相位差。在不考慮折射率變化的情況下有相位差關(guān)系：

由于光在非均勻折射率場(chǎng)中傳播時(shí)會(huì)受到介質(zhì)的折射率影響改變傳播速度，相干光源在經(jīng)過(guò)特定的區(qū)域后，產(chǎn)生的相位差會(huì)發(fā)生變化。由于干涉條紋本身具有尺度的特性，在固定位置的紋影測(cè)量實(shí)驗(yàn)中，使用條紋干涉紋影測(cè)量的方法能夠得到條紋偏移的真實(shí)距離，從而解算相位差。

條紋是根據(jù)相干光源疊加產(chǎn)生的，即等波長(zhǎng)，相同傳播方向，不同相位的光線疊加，可以看作兩個(gè)振幅為a1的矢量波E1和E2的疊加：

其中：ωt由于光源波長(zhǎng)相同，頻率一致，∫krdr為光線經(jīng)過(guò)高溫區(qū)域產(chǎn)生的相位變化值。

直接通過(guò)匹配光強(qiáng)條紋的偏移求得條紋偏移量，進(jìn)而推出相位差的變化。相干光的相位變化值，可以得到光程差變化量，進(jìn)而得到光路直線上折射率值的積分。

2.2 基于干涉條紋變形的軸對(duì)稱溫度場(chǎng)重構(gòu)

在理想狀態(tài)下，火焰溫度場(chǎng)呈軸對(duì)稱均勻分布。因此當(dāng)火焰處于測(cè)量位置中心時(shí)，受到溫度影響的光線折射率分布也均勻軸對(duì)稱，從各個(gè)視角得到的紋影圖像結(jié)果相同。假設(shè)火焰滿足理想溫度分布條件，就可以利用一個(gè)角度的測(cè)量信息反解出火焰的三維溫度場(chǎng)分布。

對(duì)于軸對(duì)稱的介質(zhì)而言，Abel 反演公式可以實(shí)現(xiàn)從一維投影場(chǎng)重建二維平面結(jié)果。構(gòu)建一個(gè)高斯分布的折射率場(chǎng)n，則有光線通過(guò)的積分路徑如圖2 所示。高斯分布的折射率場(chǎng)，有距離中心距離相同處折射率相同的對(duì)稱特點(diǎn)，由于相干光線位置接近，實(shí)際產(chǎn)生的相位差等價(jià)于由不同的積分半徑造成的積分差，其中r表示被測(cè)平面任意位置的軸向半徑，R表示測(cè)量區(qū)域內(nèi)的最大軸向半徑。

圖2 光線通過(guò)軸對(duì)稱折射率場(chǎng)的相位差積分原理Fig.2 Principle of phase difference integration of light through an axisymmetric refractive index field

根據(jù)干涉條紋的產(chǎn)生性質(zhì)，相位差是光程差與光頻率的關(guān)系，當(dāng)光線穿過(guò)軸對(duì)稱折射率場(chǎng)時(shí)，由于光線的偏折角極小，在x處產(chǎn)生的相位差δ可以根據(jù)公式（8）計(jì)算：

其中的變量n為歸一化折射率，即以測(cè)試場(chǎng)周圍環(huán)境折射率為1 等比例求解的折射率分布。此歸一化的折射率n分布可以由Abel 逆變換求得：

根據(jù)此公式求解積分函數(shù)時(shí)會(huì)存在奇點(diǎn)問(wèn)題，因此求解時(shí)采用數(shù)值近似方法求解。設(shè)定兩束相干光線的距離固定，且相干光線相鄰兩段Δr=rj-ri的相位差變化量為δj，積分得到n（ri）為不同位置的折射率梯度。變量i和j表示光線發(fā)生偏折區(qū)域的順序數(shù)據(jù)點(diǎn)，即積分路徑上的順序。線性算子系數(shù)Dij與相干光線間距Δr無(wú)關(guān)系，利用反演算法可以確定：

將式（10）帶入變換公式中，利用辛普森1/3rd進(jìn)行求解：

帶入得到線性算子系數(shù)Dij：

通過(guò)Dij外推后，可以得到消除i=j處的奇點(diǎn)，得到結(jié)果：

此時(shí)得到的Dij為相位變化值與折射率之間的關(guān)系。其中，1 標(biāo)準(zhǔn)大氣壓下，波長(zhǎng)為644 nm的可見光在15 ℃且CO2含量450 ppm 的干燥空氣中折射率是1.000 276 4。空氣中光的折射率n受到傳播介質(zhì)的溫度、密度以及光線的波長(zhǎng)等條件限制，對(duì)于光線的折射率計(jì)算存在眾多經(jīng)驗(yàn)公式，符合實(shí)驗(yàn)室環(huán)境條件選取了Ciddor 公式［19］。將水蒸氣含量百分比e=0，大氣壓P=101 325 Pa，光波長(zhǎng)λ=0.54 μm，帶入可以得到空氣的折射率與溫度之間的反比關(guān)系：

2.3 條紋系統(tǒng)對(duì)溫度場(chǎng)重構(gòu)的影響

圖3 狹縫對(duì)于條紋影響示意圖Fig.3 Diagram of the effect of slits on stripes

因?yàn)闂l紋圖案是根據(jù)雙縫干涉產(chǎn)生的，條紋間距反應(yīng)了光程差的信息。在系統(tǒng)構(gòu)型固定的情況下，條紋變化就可以反映該測(cè)量空間內(nèi)空氣折射率的變化信息，因此在固定系統(tǒng)構(gòu)型時(shí)需要充分考慮干涉條紋的影響。由于激光干涉的特殊性質(zhì)，在衍射光線夾角超過(guò)6°后，得到的干涉條紋的條紋間距就已經(jīng)不能視為寬度相等，且隨著條紋亮度衰減，對(duì)于不同寬度的雙縫片的選取在實(shí)驗(yàn)中顯得更加重要。在成像平面上，激光器發(fā)出的相干光的衍射效果可以近似為夫瑯和費(fèi)衍射，利用近似公式可以得到條紋的光強(qiáng)分布。當(dāng)相干光線照射在寬度為d雙縫上，根據(jù)惠更斯-菲涅耳原理，雙縫上每個(gè)間隙都可以視為無(wú)數(shù)個(gè)新的子波源，子波源向各方向發(fā)出球面次波，這些次波在成像平面上疊加，可以得到平行于狹縫的條紋圖樣。a為雙縫中每一條縫的寬度，λ為此激光波長(zhǎng)。對(duì)于是夫瑯和費(fèi)近似，有成像平面上光強(qiáng)的公式［20］：

可以看出，紋影系統(tǒng)構(gòu)型對(duì)條紋紋影法測(cè)量精度存在影響。根據(jù)公式（17），對(duì)影響成像平面上光強(qiáng)的變量進(jìn)行仿真比較。

當(dāng)波長(zhǎng)λ=560 nm 時(shí)，對(duì)于夫瑯和費(fèi)單縫衍射的光線強(qiáng)度可以計(jì)算得出。取雙縫中心對(duì)應(yīng)成像平面為條紋中心，成像平面距離l+L=800 mm，對(duì)成像平面上條紋圖像進(jìn)行窗口為1的高斯濾波，得到不同狹縫的仿真結(jié)果如圖4所示。

假設(shè)一束相干光通過(guò)（0，0）位置的雙縫經(jīng)過(guò)溫度場(chǎng)，選取兩組成像平面，根據(jù)公式（17）可以得到光線穿過(guò)折射率場(chǎng)后產(chǎn)生的相應(yīng)的光程差，如圖5 所示。

根據(jù)折射率場(chǎng)存在前后光程差變化，根據(jù)激光波長(zhǎng)固定的原理，可以求解相位差，通過(guò)切面驗(yàn)證得到條紋偏移量結(jié)果如圖6 所示。

圖4 不同寬度雙縫衍射對(duì)應(yīng)光線強(qiáng)度對(duì)比圖Fig.4 Comparison of the light intensity corresponding to different widths of double slit diffraction

圖5 不同成像平面的干涉條紋光程差Fig.5 Interference fringe optical range difference of different imaging planes

可以看出，調(diào)節(jié)干涉系統(tǒng)的大小對(duì)偏移結(jié)果影響成放大倍數(shù)關(guān)系。當(dāng)固定雙縫和成像平面的距離，調(diào)節(jié)火焰在干涉系統(tǒng)中位置時(shí)，條紋整體寬度無(wú)影響，但隨折射率變化受到很大影響，如圖7 所示，進(jìn)而可以得到圖8 所示的干涉條紋的偏移量。

圖6 不同成像平面的干涉條紋偏移量Fig.6 Interference fringe optical range difference of different imaging planes

圖7 不同火焰位置的干涉條紋光程差Fig.7 Interference fringe optical path difference of different flame positions

圖8 不同火焰位置的干涉條紋偏移量Fig.8 Interference fringe offsets for different flame positions

當(dāng)雙縫和成像平面位置固定時(shí)，溫度場(chǎng)距離光源越近在成像平面上得到的條紋偏移越明顯。由于光線偏折角的角度很小，因此可以將光線傳播路徑近似為直線。由于激光傳播的性質(zhì)，光程差就是路徑上折射率的積分。因此可以根據(jù)條紋偏移結(jié)果解算出火焰溫度分布。

3 干涉條紋紋影測(cè)溫的數(shù)值模擬

由2.3 節(jié)可知，不同的火焰中心溫度對(duì)條紋的偏移量存在影響。設(shè)定高斯溫度場(chǎng)火焰中心位于（1 000，0）處，雙縫寬度為0.2，成像平面位于y=1 500 處（單位：mm），分別設(shè)定兩組軸對(duì)稱溫度場(chǎng)，距離中心的溫度如圖9（a）所示，并根據(jù)經(jīng)驗(yàn)公式求得兩組溫度場(chǎng)對(duì)應(yīng)的折射率分布。

圖9 數(shù)值仿真結(jié)果Fig.9 Numerical simulation

將折射率分布帶入仿真程序中，可以得到不同中心溫度火焰對(duì)應(yīng)的干涉條紋偏移效果如圖10 所示?？梢钥闯?，不同溫度的軸對(duì)稱場(chǎng)對(duì)應(yīng)的折射率場(chǎng)不同，進(jìn)而導(dǎo)致干涉條紋的最大偏移量位置，以及偏移量都不同。

圖10 不同火焰溫度的干涉條紋偏移量Fig.10 Interference fringe offsets for different flame temperatures

可以利用δ（r）采用上述Abel 逆變換，可以反演求得空間中的折射率分布n（r），由于該給定溫度場(chǎng)是標(biāo)準(zhǔn)高斯分布，選取中心角度折射率值進(jìn)行對(duì)比。將中心溫度為800 K 的條紋偏移量帶入公式（6）中，可以計(jì)算折射率分布，徑向溫度分布圖如圖所示，由于區(qū)域范圍對(duì)應(yīng)實(shí)際寬度為20 mm 之內(nèi)，選取200 組點(diǎn)進(jìn)行擬合對(duì)比，結(jié)果如圖11（b）所示。

圖11 不同中心溫度下計(jì)算結(jié)果對(duì)比Fig.11 Comparison of calculation results at different center temperatures

分別選取4 種不同的中心溫度，采用仿真方法計(jì)算條紋偏移對(duì)其產(chǎn)生的影響，并采用Abel 重建獲得溫度分布，如圖11 所示。此處產(chǎn)生最大誤差的原因是投影過(guò)程中采用了離散取值并進(jìn)行的Abel 重建，反演算法時(shí)選取了邊緣為可靠已知量，累計(jì)中會(huì)產(chǎn)生累計(jì)誤差，當(dāng)投影選取的離散間隔足夠小時(shí)，誤差可以消除。

通過(guò)圖11 的數(shù)據(jù)對(duì)比，已知條紋變化量，可以得知對(duì)應(yīng)溫度分布中心溫度。在這種溫度場(chǎng)變化趨勢(shì)固定的情況下，對(duì)中心溫度多次求偏移擬合可以得到圖12 所示的最大偏移量與溫度的關(guān)系。因此可以根據(jù)測(cè)量得到的最大偏移量推導(dǎo)出中心溫度。

圖12 條紋偏移量最大值隨中心溫度的變化Fig.12 Variation curve of maximum value of stripe deflection with center temperature

4 軸對(duì)稱溫度場(chǎng)條紋紋影重建測(cè)量實(shí)驗(yàn)

用酒精燈模擬軸對(duì)稱溫度場(chǎng)，其中由于火焰燃燒存在的不穩(wěn)定性，溫度真實(shí)結(jié)果是采用熱電偶對(duì)火焰固定區(qū)域進(jìn)行長(zhǎng)時(shí)間測(cè)量的平均值。在實(shí)驗(yàn)室條件下，采用高精度鉑銥熱電偶（WRP-130S，測(cè)量溫度范圍為0～1 600 ℃）配備XMTA-308 型數(shù)字溫度顯示器以及高精度三坐標(biāo)平移臺(tái)對(duì)酒精燈火焰進(jìn)行對(duì)應(yīng)坐標(biāo)的溫度驗(yàn)證。利用條紋紋影法和背景紋影法同時(shí)重建酒精燈的火焰截面，輔以熱電偶測(cè)得的實(shí)際溫度值。搭建干涉條紋紋影測(cè)溫系統(tǒng)如圖13。

圖13 實(shí)驗(yàn)系統(tǒng)圖Fig.13 Experimental system diagram

在紋影系統(tǒng)之中，由于溫度影響產(chǎn)生的折射率場(chǎng)的存在，激光器在成像平面上得到的干涉條紋會(huì)發(fā)生抖動(dòng)，預(yù)先標(biāo)定偏移量對(duì)激光頻率造成的影響進(jìn)行補(bǔ)償。本實(shí)驗(yàn)中使用功率為50 mW的激光模組進(jìn)行原理驗(yàn)證。如圖14 所示，先進(jìn)行小尺度（平面）驗(yàn)證，當(dāng)波長(zhǎng)為540 nm 激光時(shí)，在成像平面上得到的干涉圖樣條紋寬度在3 mm左右。

圖14 干涉條紋圖案Fig.14 Interference stripe pattern

由于50 mW 的激光經(jīng)過(guò)衍射后，中央亮條紋整體亮度較高，得到條紋的光強(qiáng)分布曲線圖15，即干擾偏移量最大值為4 pix，約單個(gè)條紋間距的10%。通過(guò)激光準(zhǔn)直，可以精確地測(cè)量固定熱電偶到火焰截面處。相機(jī)拍攝的干涉條紋由于明暗特性自身帶有寬度信息。

圖15 干涉條紋的灰度分布Fig.15 Grayscale distribution of interfering stripes

通過(guò)與無(wú)火焰狀態(tài)下條紋光強(qiáng)的波形對(duì)比，利用相位角和光強(qiáng)成正比的關(guān)系，可以求出全部拍攝區(qū)域的干涉條紋偏移量如圖16。根據(jù)推導(dǎo)的雙縫干涉的光強(qiáng)公式，可以根據(jù)條紋的光強(qiáng)衰減變化即可求得相位差。

圖16 干涉條紋的偏移量Fig.16 Offset of interference stripes

根據(jù)條紋的偏移量，可以積分得到歸一化折射率的梯度分布，進(jìn)而得到的重建溫度場(chǎng)如圖17所示。根據(jù)條紋寬度，可以確定相可以確定相機(jī)拍攝的范圍為20 cm，而由于酒精燈位于狹縫和成像平面中間位置，即重建得到的火焰邊界為10 cm 左右。

圖17 酒精燈火焰重建溫度場(chǎng)Fig.17 Reconstruction of temperature field by flame of alcohol lamp

采用背景紋影法重建火焰截面與此方法進(jìn)行比較，如圖18 所示。

實(shí)驗(yàn)將熱電偶固定在高精度平移臺(tái)上，每次沿同一方向移動(dòng)5.00 mm，根據(jù)熱電偶測(cè)量的酒精燈溫度，和重建得到的溫度場(chǎng)進(jìn)行對(duì)比。結(jié)果如圖19 所示。

圖18 背景紋影測(cè)量實(shí)驗(yàn)場(chǎng)景圖Fig.18 Background-Oriented Schlieren（BOS）system scene map

圖19 酒精燈火焰重建結(jié)果圖Fig.19 Alcohol lamp flame reconstruction results graph

考慮到熱電偶測(cè)溫存在傳熱不均勻的問(wèn)題，整體重建結(jié)果較好。對(duì)于酒精燈火焰溫度場(chǎng)的重建，背景紋影法平均誤差為22.7 ℃，對(duì)應(yīng)的測(cè)量誤差范圍為0.5 mm；相比條紋紋影法測(cè)量范圍可達(dá)到0.2 mm，對(duì)應(yīng)平均誤差為15.8 ℃。

5 結(jié)論

本文詳細(xì)地推導(dǎo)了干涉條紋紋影法軸對(duì)稱溫度場(chǎng)測(cè)溫的原理，搭建了一個(gè)干涉條紋紋影系統(tǒng)，通過(guò)對(duì)軸對(duì)稱溫度場(chǎng)的仿真和實(shí)驗(yàn)，重建出軸對(duì)稱火焰的溫度分布。

提出了基于干涉條紋紋影系統(tǒng)的軸對(duì)稱溫度測(cè)量方法。通過(guò)相機(jī)拍攝激光干涉得到的偏移條紋信號(hào)，進(jìn)而推斷出待測(cè)區(qū)域的折射率場(chǎng)分布。利用相機(jī)的成像，與傳統(tǒng)激光光譜測(cè)溫方法相比，增大了測(cè)量區(qū)域。通過(guò)引入激光干涉，簡(jiǎn)化了傳統(tǒng)紋影技術(shù)的復(fù)雜裝置的同時(shí)，保留了紋影技術(shù)高分辨率的特點(diǎn)。數(shù)值模擬和對(duì)比實(shí)驗(yàn)的結(jié)果表明，本文提出的方法可以實(shí)現(xiàn)軸對(duì)稱場(chǎng)的溫度測(cè)量。在后續(xù)研究中，將進(jìn)一步研究基于多角度層析重構(gòu)的非軸對(duì)稱溫度場(chǎng)重構(gòu)方法。