王根娟, 楊臧健, 孟 晟, 王明曉, 鐘英杰

(浙江工業(yè)大學(xué) 能源與動力工程研究所, 杭州 310014)

背景紋影定量化在層流軸對稱火焰溫度場測量中的應(yīng)用研究

王根娟, 楊臧健, 孟 晟, 王明曉, 鐘英杰*

(浙江工業(yè)大學(xué) 能源與動力工程研究所, 杭州 310014)

本文以本生型甲烷/空氣層流預(yù)混火焰為研究對象，研究了背景紋影技術(shù)在層流軸對稱火焰溫度場測量中的應(yīng)用?？紤]到背景尺度對窗口和相機參數(shù)的限制問題，采用了多尺度小波噪點背景。比較各類運動圖像處理技術(shù)的特點，選用變分光流法獲取光線穿過火焰后的偏轉(zhuǎn)角。搭建實驗臺并進行背景紋影火焰測溫實驗，實驗中發(fā)現(xiàn)，在選用多尺度小波噪點背景的情況下，由變分光流算法獲得的像素位移分布圖的噪聲小于同等條件下由互相關(guān)算法得到的結(jié)果。最后，假設(shè)火焰呈軸對稱分布，結(jié)合Gladstone-Dale公式與理想氣體狀態(tài)方程分別獲得了甲烷火焰當(dāng)量比為1.06和0.83這2種實驗條件下的溫度場，所獲得的溫度分布與Raman-LIF法的測溫結(jié)果相比，趨勢基本一致。

背景紋影；火焰溫度；變分光流算法；互相關(guān)算法；多尺度小波噪點背景

0 引 言

傳統(tǒng)紋影是通過一系列的光學(xué)器件組合，利用穿過流場光線的偏轉(zhuǎn)，將人眼無法觀察到的流場內(nèi)密度變化展現(xiàn)為圖像的技術(shù)。傳統(tǒng)紋影雖然在可視化流場的觀測中已有著廣泛的應(yīng)用，但它還存在著諸多局限性，如傳統(tǒng)紋影有著相對復(fù)雜的系統(tǒng)、使用不靈活、成本相對較高、只能獲得流場的某一個方向上的偏轉(zhuǎn)角信息，而且鑒于其易受光源亮度、流場均勻程度、鏡片反射率、流場對光線的透射率、刀口切入深度、相機光圈大小等因素的影響，傳統(tǒng)紋影很難實現(xiàn)定量化。

21世紀(jì)初，Richard和Meier等[1-2]在傳統(tǒng)紋影的基礎(chǔ)上發(fā)展了背景紋影技術(shù)。背景紋影系統(tǒng)設(shè)置簡單、使用靈活、成本較低、可以獲得流場2個方向上的偏轉(zhuǎn)角信息，且容易定量化，在國內(nèi)外研究中已有廣泛的應(yīng)用。如Klinge等[3]聯(lián)用背景紋影和PIV技術(shù)獲得了風(fēng)洞內(nèi)流場的密度和速度分布。Goldhahn等[4]詳細地分析了背景紋影在流場測量中靈敏度和精度等問題。Elsinga等[5]將彩色紋影和背景紋影技術(shù)相結(jié)合獲得了流場中的剪切層和流體的膨脹與壓縮信息。Leopold等[6-7]用改進后的彩色背景紋影技術(shù)獲得了超音速風(fēng)洞內(nèi)流場和直升機翼片周圍的流場。Ramanah等[8]將背景紋影技術(shù)應(yīng)用于超聲速沖擊波流場的顯示。Kotchourkoa等[9]利用背景紋影獲得了氫氣射流的密度。周昊等[10]研究了背景紋影技術(shù)在火焰溫度場測量中的應(yīng)用，發(fā)現(xiàn)所得的溫度場符合實際情況。

背景紋影算法方面，多數(shù)學(xué)者采用PIV技術(shù)中的互相關(guān)算法來分析像素偏移情況。少數(shù)學(xué)者引入了光流算法，如Atcheson B[11]等詳細地介紹了光流算法，并比較分析了各種光流算法在背景紋影中應(yīng)用情況，發(fā)現(xiàn)由于光流算法在背景計算中的取樣點遠遠多于互相關(guān)算法，可大大改進所得結(jié)果。

國內(nèi)外關(guān)于背景紋影的研究表明，該技術(shù)是一種可靠的獲得流場密度分布信息的有效檢測手段，與PIV等技術(shù)的聯(lián)用，可實現(xiàn)流場的可視化測試。但是，目前此類研究的重點主要關(guān)注流場信息定性的分布及變化情況，對參數(shù)定量化的研究相對欠缺，仍在逐步開展。例如，Klinge[3]等學(xué)者已對該技術(shù)的定量化研究進行了探索。S. Tokgoz等[12]開始研究用背景紋影技術(shù)進行溫度測量，在他們的薄液層溫度測量研究中發(fā)現(xiàn)，背景紋影定量化技術(shù)不僅可測得薄液層穩(wěn)定狀態(tài)下的溫度，也可以測得其流動時的溫度。但總的來看，多數(shù)研究中的測量對象僅限于無化學(xué)反應(yīng)、無明顯溫度變化的流體中，在火焰等復(fù)雜流體中的應(yīng)用研究則相對較少，難以滿足應(yīng)用需要。

基于以上分析，有必要分析和探討背景紋影定量化技術(shù)在火焰研究中應(yīng)用的可行性及準(zhǔn)確性，嘗試從紋影測量的流場密度分布出發(fā)，基于Gladstone-Dale公式和理想氣體方程，獲得流場溫度信息。本文比較了各類常用背景的特點，采用了多尺度小波噪點背景，接著在應(yīng)用小波噪點背景的基礎(chǔ)上詳細分析了互相關(guān)算法和變分光流算法在獲得圖像像素偏移值上的優(yōu)劣，并搭建了背景紋影技術(shù)測火焰溫度的實驗系統(tǒng)，最終通過背景紋影的定量化計算方法獲得了火焰的溫度分布，與Q.V.Nguyen[13]等的實驗結(jié)果相比較具有一定的準(zhǔn)確性。

1 背景紋影原理與算法

區(qū)別于傳統(tǒng)的紋影，背景紋影的偏轉(zhuǎn)角獲得是通過對比具有一定時間間隔前后拍攝的2張圖片上像素差別得到的，其原理如圖1所示，實線表示無火焰時從背景發(fā)出到相機成像平面上的光線，虛線表示火焰存在時，因流場折射率變化發(fā)生偏轉(zhuǎn)后進入到相機成像平面上的光線。通過比較前后2圖像上位置的變化，可以得到偏轉(zhuǎn)角分布。

圖1 背景紋影原理示意圖

光線經(jīng)過一個折射率有變化的流場，就會發(fā)生偏轉(zhuǎn)，其偏轉(zhuǎn)角θ和折射率n的關(guān)系可用下式表述[14]：

(1)

其中，x方向為光線入射方向。紋影正是利用這一原理工作的，它反映了流場折射率沿空間一階導(dǎo)數(shù)的變化。而本文所用的是層流本生型火焰，可近似認為其流場是軸對稱分布的，那么，式(1)中的偏轉(zhuǎn)角θ和折射率n的關(guān)系可改寫成：

(2)

式中:rn指火焰截面最外圈半徑，ri指θi所對應(yīng)的半徑，y等于ri。

通過背景紋影獲得流場折射率信息，進一步可處理獲得溫度分布，實現(xiàn)過程首先需要假設(shè)條件并計算獲得定量化的流場密度分布，具體方法是通過計算前后2張背景圖片的像素偏移來獲得光線偏轉(zhuǎn)角。通過式(1)可由光線偏轉(zhuǎn)角獲得折射率。在得到流場的折射率后，可以通過Gladstone-Dale公式由折射率獲得密度。在獲得密度的基礎(chǔ)上，結(jié)合理想氣體狀態(tài)方程可以將溫度和密度聯(lián)系起來，完成整個由背景紋影定量化獲得流場溫度的過程[10]。

1.1 多尺度小波噪點背景

背景的選擇在背景紋影中占有很重要的地位。幾個常用的背景如圖2所示。大多數(shù)背景紋影會選擇隨機噪點[11]。隨機噪點和有規(guī)則噪點的最大問題是它們具有尺度限制。特別是應(yīng)用互相關(guān)算法時，一個窗口的粒子數(shù)不宜太多也不宜太少，那么為了適應(yīng)不同尺度的噪點，每次取窗口大小時都需要預(yù)先判定，具有一定的隨機性。尺度限制同樣影響相機的設(shè)定，如果背景離相機太遠，相機鏡頭的分辨率和焦距就需要相應(yīng)提高，否則特定尺度的隨機噪點或有規(guī)則噪點背景就會失去作用。

(a) 隨機噪點 (b)有規(guī)則噪點 (c)多尺度小波噪點

圖2 幾個常用的背景

Fig.2 Common backgrounds

由于多尺度小波噪點背景更接近自然背景，連續(xù)性更好，因此它的引入可以有效地解決隨機噪點背景和有規(guī)則噪點背景的尺度限制問題，適用性更廣。此外，多尺度小波噪點背景可以使得應(yīng)用變分光流算法時獲得的像素流場更加光滑。

二維多尺度小波噪點背景圖像生成步驟如下，其流程圖如圖3所示：

(1) 首先確定背景圖像尺寸，設(shè)為2n×2n。

(2) 取某一分辨率下的尺寸2k×2k，為k

(3) 將該圖像進行小波分析，取其低通信號，即對圖片進行降采樣，記為Idown(k-1)。

(4) 再將降采樣后的圖片進行小波重構(gòu)，升采樣為原始圖片尺寸，記為Idownup(k)。

(5) 最后將所有n-1個尺度的圖片都升采樣至尺度n，取平均權(quán)重疊加在一起就形成了多尺度小波噪點圖片。

本文中，為了方便圖像尺寸的縮放采用“Haar”小波作為變換小波。

圖3 多尺度小波噪點背景的生成流程

1.2 互相關(guān)算法

在利用前后2張背景圖片像素偏移值的計算中，需要利用圖像運動分析算法，目前常用的算法是互相關(guān)算法，同時也有少數(shù)學(xué)者引用了光流算法。互相關(guān)算法是PIV粒子圖像處理技術(shù)中常用的粒子測速方法。它是通過計算2個查詢區(qū)內(nèi)粒子速度的互相關(guān)性來獲得粒子速度的。二維互相關(guān)函數(shù)定義如下[15]：

(3)

式中：I1、I2為2張圖像查詢區(qū)內(nèi)的像素灰度實函數(shù)。

由于互相關(guān)算法每迭代一次都需要查問整個圖像空間域，計算量大。因此，Willert和Gharib[16]提出了在互相關(guān)算法中加入快速傅里葉變化(FFT)，較大程度地提高了互相關(guān)算法的運算速度，其流程圖如圖4所示。首先在2張圖像上取一樣的查詢區(qū)，分別對查詢區(qū)內(nèi)的像素灰度函數(shù)作快速傅里葉變換，然后通過互相關(guān)函數(shù)確定查詢區(qū)間的最佳匹配，最后作傅里葉逆變換得到相應(yīng)查詢區(qū)內(nèi)的速度。

圖4 基于快速傅里葉變換的互相關(guān)算法流程圖

1.3 變分光流算法

由Horn和Schunck[17]最先引進的變分光流法是運動分析里最基本的方法。該算法假設(shè)圖像序列在短時間內(nèi)亮度恒定不變，即

(4)

式中：I(x+δx,y+δy,t+δt)是t+δt時刻的圖片像素灰度，I(x,y,t)是t時刻的圖片像素灰度。將式(4)右邊進行泰勒級數(shù)展開，并舍去高階項后，可得到光流梯度守恒方程:

(5)

1.3.1 Horn-Schunck算法

在式(5)的基礎(chǔ)上，Horn和Schunck提出了一個全局光滑約束項，為

(6)

(7)

于是，光流場的計算問題就變成了求解式(7)的最小值問題。

1.3.2 Lucas-Kanade算法

Lucas和Kanade則認為光流梯度守恒方程在局部圖像域Ω上成立，即

(8)

同樣地，該方法也存在著未知數(shù)大于方程數(shù)的問題。該問題可通過由最小二乘法得到一個近似解來解決[21]。于是，目標(biāo)函數(shù)就變?yōu)椋?/p>

(9)

最小化式(9)后，就可以得到關(guān)于灰度梯度的線性方程組，為

(10)

互相關(guān)算法與變分光流算法相比，一方面，由于互相關(guān)算法計算過程中需要取查詢區(qū)，而查詢區(qū)內(nèi)的樣本數(shù)目有限，因此由互相關(guān)算法獲得的結(jié)果在一定程度上受查詢區(qū)選取情況的影響。而變分光流算法所取的區(qū)域為整個背景的像素點，相比于互相關(guān)算法，樣本數(shù)目較龐大。因此，基于變分光流的Horn-Schunck算法和Lucas-Kanade算法能獲得更加致密的速度場，更適用于灰度變化較連續(xù)的背景。另一方面，由于Horn-Schunck算法是基于全局平滑項來控制方程收斂的，當(dāng)圖像局部粒子密度差別較大時，單一變量控制無法得到局部平滑收斂解，此時誤差會相對較大，該種情況下更適合采用Lucas-Kanade算法和互相關(guān)算法。

2 實驗設(shè)置

2.1 實驗系統(tǒng)

為了實現(xiàn)用背景紋影定量化技術(shù)獲得火焰溫度場，本文建立了如圖5所示的背景紋影實驗系統(tǒng)，其實物圖如圖6所示。該實驗系統(tǒng)可以分為燃燒實驗臺和背景紋影2部分。

圖5 背景紋影實驗系統(tǒng)圖

圖6 背景紋影實驗實物圖

燃燒實驗臺用于產(chǎn)生本生型甲烷/空氣層流預(yù)混火焰，其本生燈管口是外徑為20mm、壁厚2mm、長度為100mm的不銹鋼圓管。

背景紋影系統(tǒng)主要由相機和背景2部分組成。相機型號為佳能EOS-550D，背景圖片為多尺度小波噪點背景。為使相機的快門和感光度選擇有利于拍出低噪圖像，特在背景后加設(shè)了1個250W的鹵素?zé)簟嶒灂r，相機、火焰、背景需設(shè)置在同一直線上。由于本生燈火焰位置是固定的，因此相機和背景的位置依據(jù)本生燈火焰位置而定，相機的高度通過三腳架調(diào)整，背景的位置通過三維坐標(biāo)架配合步進電機完成，其定位精度為0.1mm。

為了滿足變分光流算法中亮度不變的假設(shè)，同時提高實驗的精度，本實驗特在以下幾個方面做了相應(yīng)的改進：(1)選擇本生型層流預(yù)混火焰為測量對象，流場相對透明，密度梯度相對較小，可減少反射和散射的影響；(2)背景后設(shè)置鹵素?zé)?，可保證光源穩(wěn)定；(3)選用的背景圖案接近漫反射，可減少從不同角度觀察時光線發(fā)生的變化；(4)為減少圖像體積折射的再聚焦效應(yīng)，選擇較小的相機光圈。

實驗中，背景與火焰的距離ZD=0.4m，背景與相機鏡頭的距離ZB=0.8m，相機鏡頭焦距fc=50mm。設(shè)背景圖像上實際像素位移為Δy′，相機拍攝的圖像像素位移為Δy，可以得到背景紋影光線偏轉(zhuǎn)角的計算公式為：

(11)

2.2 實驗結(jié)果

圖7為實驗獲得的火焰背景紋影圖像。從圖中可以看出，由互相關(guān)算法獲得的紋影圖像噪聲較大，由Horn-Schunk算法和Lucas-Kanade算法獲得的紋影圖像噪聲較小。這是因為多尺度小波噪點背景更接近自然背景，像素分布較為連續(xù)，不利于提高互相關(guān)算法的信噪比。

為了更好地比較在多尺度小波噪點背景下互相關(guān)算法、Horn-Schunk算法和Lucas-Kanade算法的計算結(jié)果，我們?nèi)【嚯x本生燈管口上方200個像素處為研究對象，如圖7(a)中的紅線所示。實驗中，互相關(guān)算法中查詢區(qū)的大小為8×8，重疊區(qū)域為查詢區(qū)大小的50%。Horn-Schunck算法中規(guī)范化因子alpha取25，迭代步數(shù)為1000，在進行迭代前先對原始圖片進行高斯濾波。Lucas-Kanade算法中鄰域選擇為20×20，同樣在進行迭代前先對原始圖片進行高斯濾波。這3種算法的火焰背景紋影像素偏移結(jié)果如圖8所示。

(a)

(b)

(c)

(d)

圖7 (a)多尺度小波噪點火焰背景紋影圖 (b)互相關(guān)算法 (c)Horn-Schunck算法 (d)Lucas-Kanade算法

Fig.7 (a)BOS picture of flame (b)cross-correlation algorithm (c)Horn-Schunck algorithm (d)Lucas-Kanade algorithm

圖8 3種算法的背景紋影像素偏移計算結(jié)果比較

從圖8中可以看出，用互相關(guān)算法計算的像素偏移結(jié)果曲線波動明顯較大，說明噪聲比較嚴重。而由Horn-Schunck算法和Lucas-Kanade算法產(chǎn)生的像素偏移結(jié)果曲線則較為平滑，說明噪聲較小，且2種算法結(jié)果基本一致。這是因為互相關(guān)算法得到的是查詢區(qū)內(nèi)粒子的平均速度，它的數(shù)據(jù)點遠小于2種變分光流算法所用的數(shù)據(jù)點。這就進一步加大了噪聲對互相關(guān)算法的影響。

從整個實驗中，我們可以發(fā)現(xiàn)，多尺度小波噪點背景與2種變分光流算法的組合，可以得到相對平滑準(zhǔn)確的像素偏移值。它不僅可以增強實驗的自適應(yīng)性，而且可以減少調(diào)試和選擇的影響。這為后續(xù)溫度場的計算提供了良好的基礎(chǔ)條件。

在獲得火焰背景紋影像素偏移值后，利用式(11)就可得到光線偏轉(zhuǎn)角。接著，利用軸對稱斷層投影算法可以得到火焰折射率分布，本文利用的是直接積分法，該方法假設(shè)在半徑ri和半徑ri+1之間，折射率成線性變化,從而折射率梯度可以表達為：

(12)

利用式(2)對每一個(ri,ri+1)區(qū)間分段直接積分，即可得到：

(13)

本實驗中，采用了2組當(dāng)量比分別為1.06和0.83的甲烷層流預(yù)混火焰為測量對象，實驗結(jié)果如圖9和10所示。

從圖9中可以看到，當(dāng)量比為1.06的甲烷層流預(yù)混火焰最高溫度約為1800K，溫度分布從管口軸線向外逐漸升高，然后又降到環(huán)境溫度。這是因為管口軸線處的來流預(yù)混氣溫度接近環(huán)境溫度，預(yù)混氣經(jīng)預(yù)熱后在火焰鋒面處發(fā)生化學(xué)反應(yīng)，溫度逐漸升高。當(dāng)溫度到達最高點后，隨著甲烷氣體被消耗殆盡，又逐漸回到環(huán)境溫度。

(a) 由Lucas-Kanade算法獲得的像素偏移值

(b) 火焰溫度垂直截面云圖

Fig.9 The temperature field of methane/air laminar premixed flame(φ=1.06)

(a) 由Lucas-Kanade算法獲得的像素偏移值

(b) 火焰溫度垂直截面云圖

Fig.10 The temperature field of methane/air laminar premixed flame (φ=0.83)

將圖10與9相比，可以發(fā)現(xiàn)當(dāng)量比為0.83和1.06的甲烷層流預(yù)混火焰的溫度分布大致相同，都是從管口軸線向外逐漸升高，然后又降到環(huán)境溫度。不同的是，當(dāng)量比為0.83的甲烷火焰最高溫度只有1400K左右，低于當(dāng)量比為1.06的甲烷火焰。

將本文所獲得的溫度場結(jié)果與Q.V.Nguyen[13]等人用Raman-LIF測定的甲烷本生型火焰溫度場結(jié)果相比，趨勢基本一致。

3 結(jié) 論

在火焰溫度的測量中，背景紋影定量化技術(shù)可以相對簡單地獲得軸對稱火焰的溫度場。本文在背景紋影實驗臺的基礎(chǔ)上，利用多尺度小波噪點背景比較了互相關(guān)算法、Horn-Schunck算法和Lucas-Kanade算法這3種算法在火焰背景紋影圖片像素偏移計算中的優(yōu)劣。最后基于軸對稱假設(shè)，再結(jié)合Gladstone-Dale公式和理想氣體假設(shè)得到了預(yù)混火焰溫度。將多尺度小波噪點背景和變分光流算法結(jié)合起來，一方面使得實驗設(shè)置更加簡單；另一方面由于多尺度小波噪點背景接近于實際自然背景的特點，也將進一步擴展背景紋影定量化技術(shù)在工程實踐中的應(yīng)用。本文獲得的結(jié)論主要如下：

(1) 隨機噪點背景具有尺度限制問題，而采用多尺度小波噪點背景更適用于背景紋影實驗，它可以使得實驗設(shè)置和調(diào)試更加簡單。

(2) 在利用多尺度小波噪點背景時，基于變分光流的Horn-Schunck算法和Lucas-Kanade算法相比于互相關(guān)算法，所獲得的圖像像素偏移值噪聲更小，數(shù)據(jù)更平滑。文中得出的這個結(jié)論，在一定程度上為大流場下背景像素偏移算法的選擇提供了新的參考。

(3) 當(dāng)量比為1.06的甲烷空氣預(yù)混火焰最高溫度為1800K，大于當(dāng)量比為0.83的甲烷預(yù)混火焰溫度。不過2者溫度分布規(guī)律基本一致，都是從管口軸線向外逐漸升高，在達到最高點后又降到環(huán)境溫度，符合本生燈火焰的特點。

[1] Richard H, Raffel M. Principle and applications of the background oriented schlieren (BOS) method[J]. Measurement Science & Technology, 2001, 12(9): 1576-1585.

[2] Richard H, Raffel M, Rein M, et al. Demonstration of the applicability of a Background Oriented Schlieren (BOS) method[C]. 10th International Symposium on Applications of Laser Techniques in Fluid Mechanics, 2002: 145-156.

[3] Klinge F, Kirmse T, Kompenhans J. Application of quantitative Background Oriented Schlieren (BOS): investigation of a wing tip vortex in a transonic wind tunnel[C]. The 4th Pacific Symposium on Flow Visualization and Image Processing, 2003: 1-11.

[4] Goldhahn E, Seume J.The background oriented schlieren technique: sensitivity, accuracy, resolution and application to a three-dimensional density field[J]. Experiments in Fluids, 2007, 43(2-3): 241-249.

[5] Elsinga G E, Oudheusden B W, Scarano F, et al. Assessment and application of quantitative schlieren methods: calibrated color schlieren and background oriented schlieren[J]. Experiments in Fluids, 2004, 36(2): 309-325.

[6] Leopold F. The application of the Colored Background Oriented Schlieren Technique (CBOS) to free-flight and in-flight measurements[J]. Journal of Flow Visualization and Image Processing, 2009, 16(4): 279-293.

[7] Sourgen F, Leopold F, Klatt D. Reconstruction of the density field using the Colored Background Oriented Schlieren Technique (CBOS)[J]. Optics and Lasers in Engineering, 2012, 50(1): 29-38.

[8] Ramanah D, Raghunath S, Mee D J, et al. Background oriented schlieren for flow visualisation in hypersonic impulse facilities[J]. Shock Waves, 2007, 17(1): 65-70.

[9] Kotchourkoa N, Kuznetsovb M, Kotchourkob A, et al. Concentration measurements in a round hydrogen jet using Background Oriented Schlieren (BOS) technique[J]. International Journal of Hydrogen Energy, 2014, 39(11): 6201-6209.

[10] 周昊, 呂小亮, 李清毅, 等. 應(yīng)用背景紋影技術(shù)的溫度場測量[J]. 中國電機工程學(xué)報, 2011, 31(5): 63-67.

Zhou H, Lyu X L, Li Q Y, et al. Temperature measurements using the background oriented schlieren technique[J]. Proceedings of the CSEE, 2011, 31(5): 63-67.

[11] Atcheson B, Heidrich W, Ihrke I. An evaluation of optical flow algorithms for background oriented schlieren imaging[J]. Experiments in Fluids, 2009, 46(3): 467-476.

[12] Tokgoz S, Geisler R, Van Bokhoven L J A, et al. Temperature and velocity measurements in a fluid layer using background-oriented schlieren and PIV methods[J]. Measurement Science and Technology, 2012, 23: 1-10.

[13] Nguyen Q V, Dibble R W, Carter C D, et al. Raman-LIF measurements of temperature, major species, OH, and NO in a methane-air Bunsen flame[J]. Combustion and Flame, 1996, 105: 499-510.

[14] Settles G. Schlieren and shadowgraph techniques[M]. Berlin: Springer, 2001.

[15] Barron J L, Fleet D J, Beauchemin S S. Performance of optical flow techniques[J]. International Journal of Computer Vision, 1994, 12(1): 43-77.

[16] Willert C E, Gharib M. Digital particle image velocimetry[J]. Experiments in Fluids, 1991, 10(4): 181-193.

[17] Horn B K, Schunck B. Determining optical flow[J]. Artificial Intelligence, 1981, 17(1): 185-203.

[18] Nagel H H. Displacement vectors derived from second-order intensity variations in image sequences[J]. Computer Vision, Graphics and Image Processing, 1983, 21(1): 85-117.

[19] Lucas B D, Kanade T. Iterative image registration technique with an application to stereo vision[C]//Proceedings of the 7th International Joint Conference on Artificial Intelligence, 1981, 2: 674-679.

[20] Kearney J K, Thompson W B, Boley D L. Optical flow estimation: an error analysis of gradient-based methods with local optimization[J]. IEEE, 1987, PAMI-9(2): 229-244.

[21] Lucas B D. Generalized image matching by the method of differences[D]. Pittsburgh: Carnegie-Mellon University, 1984.

[22] Qin Xiao, Xiao Xudong, Ishwar K P, et al. Effect of varying composition on temperature reconstructions obtained from refractive index measurements in flames[J]. Combustion and Flame, 2002, 128(1-2): 121-132.

(編輯：張巧蕓)

Application of quantitative background oriented schlieren in laminar axisymmetric flame temperature measurement

Wang Genjuan, Yang Zangjian, Meng Sheng, Wang Mingxiao, Zhong Yingjie*

(Institute of Energy and Power Engineering, Zhejiang University of Technology, Hangzhou 310014, China)

Background Oriented Schlieren (BOS) is a new type of schlieren with the advantages of being simple, flexible and cheap. The application of BOS in the measurement of the laminar axisymmetric flame temperature field is studied with the methane/air laminar premixed flame. Random or regular noise is usually used as the background of BOS, but this would limit the interrogation area size and camera parameters. Therefore the multi-scale wavelet noise background is adopted in this paper. The cross-correlation method and the variational optical flow method under the constraint condition of Horn-Schunck and Lucas-Kanade are compared. A BOS experiment platform which can generate a methane/air laminar premixed flame is built, and BOS images are obtained. The pixel shift distribution is calculated from the BOS image and the result shows that the noise of the variational optical flow method is much smaller than that of the cross-correlation method under the same condition when the multi-scale wavelet noise background is used. Finally, under the axisymmetric assumption, the temperature distribution of the flame is obtained from deflection angle using the Gladstone-Dale relationship and the ideal gas state equation. In the experiment, the flame temperature fields are measured at the equivalent ratios ofφ=1.06 andφ=0.83, respectively. It is basically consistent with the result acquired by Raman-LIF.

background oriented schlieren;flame temperature;variational optical flow algorithms;cross-correlation algorithms;multi-scale wavelet noise background

1672-9897(2016)02-0103-08

10.11729/syltlx20150083

2015-05-26;

2015-09-20

浙江省科技廳資助項目(2014C31034)

WangGJ,YangZJ,MengS,etal.Applicationofquantitativebackgroundorientedschliereninlaminaraxisymmetricflametemperaturemeasurement.JournalofExperimentsinFluidMechanics, 2016, 30(2): 103-110. 王根娟, 楊臧健, 孟 晟, 等. 背景紋影定量化在層流軸對稱火焰溫度場測量中的應(yīng)用研究. 實驗流體力學(xué), 2016, 30(2): 103-110.

TK311

A

王根娟(1990-)，女，浙江紹興人，碩士研究生。研究方向：火焰流場參數(shù)測試、燃燒診斷。通信地址：浙江省杭州市下城區(qū)潮王路18號浙江工業(yè)大學(xué)能源與動力工程研究所(310014)。E-mail：wanggenjuan0520@163.com

*通信作者 E-mail: zhong_yingjie@zjut.edu.cn