樓波，許建紅，林振冠，徐毅

?

直線移動(dòng)火源擴(kuò)散火焰加速運(yùn)動(dòng)下的傾角特性

樓波，許建紅，林振冠，徐毅

(華南理工大學(xué)電力學(xué)院，廣東廣州，510640)

利用高速攝像實(shí)時(shí)獲得勻加速直線運(yùn)動(dòng)下移動(dòng)火源擴(kuò)散火焰的圖像序列，結(jié)合數(shù)字圖像處理和力場(chǎng)分析研究火焰傾角受速度、加速度影響的變化規(guī)律。研究結(jié)果表明：處于非慣性系的火焰傾角增加規(guī)律經(jīng)歷3個(gè)不同速率的變化階段；處于非慣性系的火焰傾角與火焰合力傾角并不是在同一方向上，火焰合力傾角要先于火焰傾角發(fā)生偏移，二者之間存在夾角，且大加速度工況下夾角大，隨著速度增加，空氣阻力增加，夾角減小，火焰傾角與火焰合力傾角同時(shí)趨近90°。

移動(dòng)火源；火焰傾角；火焰力場(chǎng)；慣性力

圖像攝像儀結(jié)合數(shù)字圖像處理技術(shù)是研究燃燒火焰特征的有效方法[1?7]。Fan等[1]利用高速攝像儀捕獲了火焰在熱中尺度石英玻璃制成的徑向通道的燃燒動(dòng)態(tài)過(guò)程，分析了非對(duì)稱火焰的燃燒特征。白衛(wèi)東等[4?5]通過(guò)對(duì)火焰圖像處理提取了火焰亮度、質(zhì)心偏移距離和圓形度等特征量，借用狀態(tài)判斷算法對(duì)火焰進(jìn)行識(shí)別和檢測(cè)。吳晉湘等[6]利用實(shí)時(shí)攝像的方法分析了不同轉(zhuǎn)速和氣流射流角產(chǎn)生離心力和科氏力對(duì)火焰結(jié)構(gòu)的實(shí)驗(yàn)影響。蔣曉剛等[7]通過(guò)實(shí)驗(yàn)對(duì)人工擾流條件下的丁烷射流燃燒形成的火焰傾角進(jìn)行實(shí)例分析，探討空氣擾流速度變化對(duì)火焰傾角的影響。但目前國(guó)內(nèi)外學(xué)者對(duì)移動(dòng)火源擴(kuò)散火焰[8?10]的基礎(chǔ)研究較少，仍未很清楚了解移動(dòng)火源擴(kuò)散火焰的內(nèi)在本質(zhì)和燃燒機(jī)理，本文作者通過(guò)分析運(yùn)動(dòng)速度和力場(chǎng)對(duì)火焰形態(tài)的影響規(guī)律，為進(jìn)一步認(rèn)識(shí)移動(dòng)火源下火焰特性提供一定依據(jù)。與靜止的火焰不同，移動(dòng)火源是指著火時(shí)火源點(diǎn)處于運(yùn)動(dòng)狀態(tài)。隨著各種交通工具的增多及其速度的加快，表面起火(移動(dòng)火源擴(kuò)散火焰)發(fā)生的機(jī)率也在增加。Lou等[8?10]在移動(dòng)火源擴(kuò)散火焰前期研究工作中，對(duì)圓周運(yùn)動(dòng)下和直線勻速運(yùn)動(dòng)下移動(dòng)火源進(jìn)行描述并通過(guò)實(shí)驗(yàn)揭示周圍氣流變化對(duì)蠟燭擴(kuò)散火焰影響的規(guī)律。本文作者利用無(wú)限大空間勻加速直線運(yùn)動(dòng)火源燃燒實(shí)驗(yàn)平臺(tái)和美國(guó)柯達(dá)公司Motion Corder Analyzer高速攝像儀對(duì)移動(dòng)火源擴(kuò)散火焰進(jìn)行拍攝和記錄，研究其加速運(yùn)動(dòng)下的傾角特性。

1 實(shí)驗(yàn)

1.1 實(shí)驗(yàn)系統(tǒng)和方法

圖1所示為勻加速直線移動(dòng)火源擴(kuò)散火焰實(shí)驗(yàn)系統(tǒng)圖，燈芯固定在滑道小車上，小車由繩子牽引，繩子的另一端經(jīng)過(guò)定滑輪和滑輪組，通過(guò)砝碼提供動(dòng)力(圖上重物)，使小車做勻加速直線運(yùn)動(dòng)。通過(guò)調(diào)節(jié)高速攝像儀與火源距離，確保能拍攝火源清晰圖像。拍攝火源從靜止開(kāi)始移動(dòng)的整個(gè)過(guò)程。試驗(yàn)中，通過(guò)改變砝碼的質(zhì)量來(lái)改變加速度。在滑道上確定一段拍攝范圍，設(shè)置攝像儀拍攝頻率，本實(shí)驗(yàn)圖像拍攝幀數(shù)設(shè)置為500幀/s，根據(jù)圖像的張數(shù)和位置確定火源的速度和加速度。本文研究火源速度范圍為0.012~1.560 m/s，加速度分別為2.5，3.6，5.4和6.5 m/s2，探討了火源運(yùn)動(dòng)狀態(tài)火焰傾角的變化。

1—高速攝像機(jī)；2—計(jì)算機(jī)主機(jī)系統(tǒng)；3—圖像采集卡；4—圖像監(jiān)視器；5—移動(dòng)火源擴(kuò)散火焰；6—柴油燃燒器；7—小車；8—滑輪組；9—砝碼；10—定滑輪

1.2 實(shí)驗(yàn)圖像

實(shí)驗(yàn)中獲取了加速度分別為2.5，3.6，5.4和6.5 m/s2的移動(dòng)火源擴(kuò)散火焰圖像。圖2所示為加速度為5.4 m/s2，速度范圍為0.054~1.350 m/s的移動(dòng)火源擴(kuò)散火焰勻加速直線運(yùn)動(dòng)圖像序列。該序列是從拍攝圖像的原始序列中每隔4張取出1張，共取25張，相鄰兩張的間隔時(shí)間為0.01 s，便于計(jì)算火焰傾角。

由圖2可見(jiàn)：速度從0.054 m/s增加到0.162 m/s時(shí)，火焰面基本呈圓柱狀形，火焰逐漸發(fā)生偏移，但較為平緩；速度從0.216 m/s增加到1.080 m/s時(shí)，火焰外形變化明顯，偏移程度加快；而速度從1.134 m/s增加到1.350 m/s時(shí)，火焰傾斜形態(tài)變化較為平緩，基本與水平方向平行。

v/(m?s?1): (a) 0.054; (b) 0.108; (c) 0.162; (d) 0.216; (e) 0.270; (f) 0.324; (g) 0.378; (h) 0.432; (i) 0.486; (j) 0.540; (k) 0.594; (l) 0.648; (m) 0.702; (n) 0.756; (o) 0.810; (p) 0.864; (q) 0.918; (r) 0.972; (s) 1.026; (t) 1.080; (u) 1.134; (v) 1.188; (w) 1.242; (x) 1.296; (y) 1.350

2 實(shí)驗(yàn)結(jié)果分析

2.1 火焰傾角

火焰在力場(chǎng)的作用下會(huì)產(chǎn)生一定的傾角，樓波等[8]曾用火焰頂點(diǎn)和底邊中點(diǎn)連線直觀表示圓周運(yùn)動(dòng)的火焰傾斜角度。但直線加速運(yùn)動(dòng)時(shí)，火焰受慣性力作用，火焰彎曲較大，邊緣形狀不規(guī)則，無(wú)法用直線法確定火焰偏轉(zhuǎn)角度。本文火焰傾角定義為：火焰區(qū)域圖像質(zhì)心與火焰根部中點(diǎn)′的連線，該連線和垂直方向的夾角即為火焰傾角，其中火焰根部中點(diǎn)′是火焰根部點(diǎn)與點(diǎn)之間線段的中點(diǎn)，傾角定義示意圖3(b)所示。文中通過(guò)matlab編程對(duì)拍攝火焰圖片用最大類間方差法[11?12]迭代出二值化圖像并進(jìn)行邊緣提取，確定和標(biāo)記火焰質(zhì)心和火焰根部中點(diǎn)，后使用E-Ruler電子測(cè)算角度工具測(cè)量并記錄火焰傾角。圖像質(zhì)心坐標(biāo)[13?14]求解公式如下：

其中：和為火焰圖片任意點(diǎn)坐標(biāo)；c和c為火焰質(zhì)心坐標(biāo)；(,)是二值函數(shù)，(,)=1特征區(qū)域，(,)=0背景?；鹧尜|(zhì)心和傾角示意圖如圖3所示。

(a) 火焰質(zhì)心圖(“+”為火焰質(zhì)心標(biāo)記)；(b) 火焰傾角定義示意圖

圖4所示為移動(dòng)火源在加速度2.5，3.6，5.4和6.5 m/s2時(shí)火焰傾角隨速度的變化曲線。由圖4可知：火焰傾角變化規(guī)律經(jīng)歷3個(gè)階段，特別在加速度高時(shí)(=6.5 m/s2)更加明顯：1) 第1階段，當(dāng)火源運(yùn)動(dòng)速度較低，火焰傾角變化趨勢(shì)較快，如加速度=6.5 m/s2的工況，火源速度從0.06 m/s增加到0.65 m/s，火焰傾角則從3°增大到61°。2) 第2階段，火焰傾角變化趨勢(shì)較第1階段慢，增加速率降低，火源速度從0.71m/s增加到1.1 m/s，火焰傾角則從65°增大到86°。3) 第3階段，此時(shí)空氣阻力進(jìn)一步增大，但火焰傾角并沒(méi)有繼續(xù)快速增加，而是趨于平緩，火源速度從1.23 m/s增加到1.56 m/s，火焰傾角接近90°。

a/(m?s?2): 1—2.5; 2—3.6; 3—5.4; 4—6.5

2.2 火焰力場(chǎng)對(duì)傾角的影響

在一個(gè)非慣性系中，當(dāng)系統(tǒng)存在一加速度時(shí)，則慣性力G遵循公式：，負(fù)號(hào)表示與系統(tǒng)運(yùn)動(dòng)方向相反。處于加速度運(yùn)動(dòng)狀態(tài)的火源主要受到3個(gè)力的作用，分別是火焰自身重力g，火焰浮力f，火焰移動(dòng)產(chǎn)生的與空氣之間的空氣阻力Z，由加速度引起的慣性力G。圖5所示是加速度為3.6 m/s2下速度分別為0.14，0.54和1.26 m/s的火焰受力分析示意圖。

v/(m?s?1): (a) 0.14; (b) 0.54; (c) 1.26

當(dāng)速度為0.14 m/s時(shí)，火焰剛從靜止開(kāi)始，因其速度小，阻力較小，相對(duì)而言，作用在火焰的浮力和慣性力較大，火焰受到合力指向左上方，火焰傾角較小(傾角18°，見(jiàn)圖4，因火源加速運(yùn)動(dòng)速度時(shí)刻變化，火焰受到合力傾角往往大于火焰傾角)。當(dāng)速度為0.46 m/s時(shí)，速度增大，阻力增大，火焰合力逆時(shí)針偏移，火焰傾角增大。速度繼續(xù)增加到1.26 m/s，火焰傾角在合力作用下繼續(xù)增大，這時(shí)空氣阻力在火焰的力場(chǎng)中占據(jù)主導(dǎo)地位。浮力和慣性力退居其次，即在火焰高速運(yùn)動(dòng)工況下，火焰的合力可近似看作是火焰所受到的空氣阻力，故使火焰往水平方向偏移，偏角接近90°。

2.3 慣性力對(duì)火焰傾角的影響

火源載體以加速度開(kāi)始運(yùn)動(dòng)的一瞬間，由于速度近乎為零，火焰在水平面上空氣阻力非常小，主要受慣性力作用，方向與火源載體運(yùn)動(dòng)方向相反。圖6所示是=0.04 s，加速度分別為2.5，3.6，5.4和6.5 m/s2的火焰開(kāi)始運(yùn)動(dòng)的一瞬間傾斜形態(tài)，在同一速度下，加速度越小的火焰在開(kāi)始運(yùn)動(dòng)的一瞬間，所受的慣性力則越小，火焰上下形態(tài)整體偏移均勻，偏移角度較大。加速度大的火焰在開(kāi)始運(yùn)動(dòng)的一瞬間，所受的水平向左的慣性力也越大，導(dǎo)致火焰下面部分移動(dòng)，而中上部區(qū)域未開(kāi)始發(fā)生偏移，整個(gè)火焰質(zhì)心偏離不大，即大加速度工況下的火焰傾角發(fā)生偏移較遲滯。圖4中火焰傾角隨速度和加速度變化曲線的初始運(yùn)動(dòng)階段正是反映了這種規(guī)律。

a/(m?s?2): (a) 2.5; (b) 3.6; (c) 5.4; (d) 6.5

2.4 火焰合力傾角

火焰載體做勻速直線運(yùn)動(dòng)，火焰所受的合力與火焰傾角在同一方向上，這時(shí)火焰合力傾角和火焰傾角保持一致。而本文研究的非慣性系的勻加速直線運(yùn)動(dòng)中，火焰傾角與火焰合力傾角并不是在同一方向上，火焰合力傾角要先于火焰傾角發(fā)生偏移，兩者之間的夾角記為，如圖7所示。在圖4中，例如對(duì)加速度=6.5 m/s2的傾角特性曲線，大加速度下水平方向慣性力大，合力傾角較大，但這時(shí)火焰下部移動(dòng)較快，中上部基本不動(dòng)，整個(gè)火焰質(zhì)心偏離不大，火焰傾角較小，兩者之間的夾角較大。隨著運(yùn)動(dòng)進(jìn)行，到第2階段，在加速度大的工況下，慣性力使火焰圖像傾角快速增加，很快趕上并超過(guò)加速度小的工況下的火焰傾角。第3階段由于速度較快，空氣阻力起決定作用，慣性力的影響減弱，合力與火焰傾角之間的夾角減小并趨近于0°，最后的火焰偏角接近并穩(wěn)定在90°位置。

圖7 火焰合力傾角α、火焰傾角β和合力與火焰傾角之間的夾角γ的示意圖

3 誤差分析

本研究的火焰傾角是通過(guò)實(shí)驗(yàn)平臺(tái)和高速攝像儀得到火焰圖像，然后在MATLAB軟件編程獲取火焰圖像的質(zhì)心和傾角。但是在火焰圖像的采集和處理的過(guò)程中，由于火焰載體的速度、拍攝角度、現(xiàn)場(chǎng)氣流干擾、火焰圖像預(yù)處理等因素的影響，實(shí)測(cè)數(shù)據(jù)存在誤差。

3.1 系統(tǒng)誤差

方法誤差主要包括：1) 拍攝的圖像是二維且豎直于水平面的，高速攝像機(jī)固定于某一位置，火焰運(yùn)動(dòng)發(fā)生位移時(shí)，鏡頭軸線并不總是垂直于火焰，因此拍攝到的圖像是火焰面在拍攝平面上的投影，比實(shí)際火焰傾斜形態(tài)小些。2) 最大類間方差法迭代出二值化圖像分割閾值會(huì)影響獲取的火焰二值化圖像，MATLAB軟件環(huán)境下計(jì)算出的火焰傾角，計(jì)算出的火焰傾角與真實(shí)火焰傾角相比，會(huì)存在誤差。

3.2 隨機(jī)誤差

攝取火焰運(yùn)動(dòng)圖像序列的現(xiàn)場(chǎng)是大空間封閉的環(huán)境，但也可能出現(xiàn)氣流干擾，影響到火焰燃燒的正常形態(tài)，導(dǎo)致高速攝像機(jī)輸出的火焰圖像可能出現(xiàn)隨機(jī)誤差。因此，在實(shí)驗(yàn)中，一方面盡量減少干擾氣流的存在；另一方面，通過(guò)增加實(shí)驗(yàn)次數(shù)取平均值消除隨機(jī)誤差。

4 結(jié)論

1) 火焰傾角變化規(guī)律經(jīng)歷3個(gè)階段：第1階段，火源速度從0.06 m/s增加到0.65 m/s，火焰傾角則從3°增大到61°，火焰傾角變化趨勢(shì)較快；第2階段，火源速度從0.71 m/s增加到1.1 m/s，火焰傾角則從65°增大到86°，火焰傾角變化趨勢(shì)變?yōu)槠骄?；?階段，當(dāng)火焰速度達(dá)到高速度1.56 m/s以后，火焰傾角逐漸接近或等于90°。在加速度大的工況下，火焰傾角第1階段滯后，第2和第3階段增速明顯。

2) 處于非慣性系的火焰傾角與火焰合力傾角并不是在同一方向上，火焰合力傾角大于火焰傾角，兩者之間存在夾角。速度小時(shí)慣性力作用顯著，夾角較大，而速度大時(shí)阻力起決定作用；夾角減小，火焰傾角與火焰合力傾角同時(shí)趨近90°。

參考文獻(xiàn)：

[1] Fan A W, Wan J L, Maruta K, et al. Flame dynamics in a heated meso-scale radial channel[J]. Proceedings of Combustion Institute, 2013, 34: 3351?3359.

[2] Fan A W, Minaev S, Kumar S, et al. Experimental investigation on flame pattern formations of DME–air mixtures in a radial microchannel[J]. Combustion and Flame, 2010, 157: 1637?1642.

[3] Fan A W, Minaev S, Kumar S, et al. Regime diagrams and characteristics of flame patterns in radial microchannels with temperature gradients[J]. Combustion and Flame, 2008, 153: 479?489.

[4] 白衛(wèi)東, 嚴(yán)建華, 馬增益, 等. 基于支持向量機(jī)的火焰狀態(tài)識(shí)別方法[J]. 動(dòng)力工程, 2004, 24(4): 548?551. BAI Weidong, YAN Jianhua, MA Zengyi, et al. Method of flame identification based on support vector machine[J]. Power Engineering, 2004, 24(4): 548?551

[5] 段悅, 袁昌明. 火災(zāi)探測(cè)中動(dòng)態(tài)火焰的數(shù)字圖像處理[J]. 中國(guó)計(jì)量學(xué)院學(xué)報(bào), 2009, 20(1): 55?58. DUAN Yue, YUAN Changming. Digital image processing of dynamical flame in fire detection[J]. Journal of China University of Metrology, 2009, 20(1): 55?88.

[6] 吳晉湘, 王恩宇, 劉聯(lián)勝, 等. 體積力場(chǎng)對(duì)預(yù)混火焰面形狀的影響[J]. 燃燒科學(xué)與技術(shù), 2005, 11(6): 535?538. WU Jinxiang, WANG Enyu, LIU Liansheng, et al .Influence of volume force field on the premixed flame shapes[J]. Journal of Combustion Science and Technology, 2005, 11(6): 535?538.

[7] 蔣曉剛, 金良安, 遲衛(wèi), 等. 人工擾流條件下噴射火焰的傾角計(jì)算[J]. 消防科學(xué)與技術(shù), 2011, 30(11): 987?990. JING Xiaogang, JIN Liangan, CHI Wei, et al. Calculative method for oblique angle of ejecting flame interfered by man-made jet flow[J]. Fire Science and Technology, 2011, 30(11): 987?990.

[8] 樓波, 徐毅, 林振冠. 圓周運(yùn)動(dòng)下移動(dòng)火源擴(kuò)散火焰的圖像特性分析[J]. 燃燒科學(xué)與技術(shù), 2013, 19(1): 60?66. LOU Bo, XU Yi, LIN Zhenguan. Images characteristic analysis of moving fire diffusion flame in circular motion[J]. Journal of Combustion Science and Technology, 2013, 19(1): 60?66.

[9] 樓波, 丁利, 龍新峰, 等. 直線運(yùn)動(dòng)下移動(dòng)火源燃燒實(shí)驗(yàn)和數(shù)值模擬[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 44(3): 1240?1245. LOU Bo, DING Li, LONG Xinfeng, et al. Mobile fire combustion experiment and numerical simulation during linear motion[J]. Journal of Central South University (Science and Technology), 2013, 44(3): 1240?1245.

[10] 樓波, 許建紅, 林振冠, 等. 直線運(yùn)動(dòng)下移動(dòng)火源擴(kuò)散火焰的分形結(jié)構(gòu)特征[J]. 中國(guó)科學(xué)技術(shù)大學(xué)學(xué)報(bào), 2013, 43(8): 679?686. LOU Bo, XU Jianhong, LIN Zhenguan, et al. Fractal structure characteristic of the diffusion flame of moving fire on linear motion[J]. Journal of University of Science and Technology of China, 2013, 43(8): 679?686.

[11] 齊麗娜, 張博, 王戰(zhàn)凱. 最大類間方差法在圖像處理中的應(yīng)用[J]. 無(wú)線電工程, 2006, 36(7): 25?26. QI Lina, ZHANG Bo, WANG Zhankai.Application of the OTSU method in image processing[J]. Radio Engineering of China, 2006, 36(7): 25?26.

[12] 王磊, 段會(huì)川. Otsu方法在多閾值圖像分割中的應(yīng)用[J]. 計(jì)算機(jī)工程與設(shè)計(jì), 2008, 29(11): 2844?2845. WANG Lei, DUAN Huichuan. Application of Otsu’ method in multi-threshold image segmentation[J]. Computer Engineering and Design, 2008, 29(11): 2844?2845.

[13] 白衛(wèi)東, 嚴(yán)建華, 池涌, 等. PCA和SVM在火焰監(jiān)測(cè)中的應(yīng)用研究[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2004, 24(2): 186?191. BAI Weidong, YAN Jianhua, CHI Yong, et al. A research on application of PCA and SVM to flame monitoring[J]. Proceedings of the CSEE, 2004, 24(2): 186?191.

[14] 方舟, 張君安. 基于小波分析實(shí)現(xiàn)圖像重心位移的檢測(cè)[J]. 空軍工程大學(xué)學(xué)報(bào), 2004, 5(6): 1009?3516. FANG Zhou, ZHANG Junan. Based on wavelet analysis to realize the image center of gravity displacement detection[J]. Journal of Air Force Engineering University, 2004, 5(6): 1009?3516.

Inclination characteristic of diffusion flame of moving fire on linear acceleration motion

LOU Bo, XU Jianhong, LIN Zhenguan, XU Yi

(School of Electric Power, South China University of Technology, Guangzhou 510640, China)

The image series of the diffusion flame of moving fire on acceleration motion were acquired by using high-speed video camera and the flame inclination characteristic was studied. The results show that flame inclination increase undergoes three stages of changing velocity in non-inertial system. It is not in same direction between flame resultant and flame inclination. In non-inertial system, the flame resultant angle shifts first to the flame inclination and forms an intersection angle. Moreover, the higher acceleration flame is, the more sluggish flame inclination will be. With the increase of flame speed, the air resistance increases and the intersection angle decreases, and the flame resultant angle and the flame inclination finally approach to 90°.

moving fire; flame inclination; flame force field; inertial force

10.11817/j.issn.1672-7207.2015.07.044

TK16

A

1672?7207(2015)07?2707?06

2014?07?19；

2014?10?13

廣東省自然科學(xué)基金資助項(xiàng)目(S2013010016748)；廣東省能源高效清潔利用重點(diǎn)實(shí)驗(yàn)室項(xiàng)目(KLB10004) (Project(2008A060301002) supported by the Natural Science Foundation of Guangdong Province, China; Project(KLB10004) supported by Key Laboratory of Efficient and Clean Energy Utilization of Guangdong Province, China)

樓波，博士，副教授，從事燃燒理論與應(yīng)用研究：E-mail: loubo@scut.edu.cn

(編輯 楊幼平)