預(yù)混湍流火焰面褶皺結(jié)構(gòu)網(wǎng)絡(luò)拓?fù)溲芯?/h1>
2018-03-16 06:26:35王金華聶要輝黃佐華
實(shí)驗(yàn)流體力學(xué) 2018年1期
王金華, 聶要輝, 常 敏, 張 猛, 黃佐華
(西安交通大學(xué) 能源與動(dòng)力工程學(xué)院, 西安 710065)
0 引 言
湍流火焰廣泛存在于燃?xì)廨啓C(jī)、航空發(fā)動(dòng)機(jī)、內(nèi)燃機(jī)、工業(yè)燃燒器等各種燃燒設(shè)備中，是寬廣時(shí)間空間尺度上湍流與火焰耦合作用的復(fù)雜物理化學(xué)問(wèn)題，其研究具有重要的工程和科學(xué)意義。根據(jù)Borghi[1]等人的湍流火焰模式分區(qū)理論，實(shí)際燃燒室湍流火焰主要位于火焰片模式及薄反應(yīng)區(qū)模式下，此時(shí)，湍流火焰的研究主要集中在湍流與局部層流火焰的相互作用。目前主要是從2個(gè)方面對(duì)此開(kāi)展實(shí)驗(yàn)研究:
一是從實(shí)驗(yàn)出發(fā)，研究湍流特性，包括湍流強(qiáng)度、積分尺度等對(duì)火焰的影響[2-4]，以及火焰化學(xué)反應(yīng)和熱釋放對(duì)湍流場(chǎng)的影響[5-6]。二是在數(shù)據(jù)處理方面，通過(guò)一定的數(shù)學(xué)統(tǒng)計(jì)分析方法對(duì)湍流火焰數(shù)據(jù)進(jìn)行提取和挖掘，找到湍流火焰中的特征參數(shù)，比如湍流火焰速度[7]、湍流火焰分形維數(shù)[8-9]等。湍流火焰的數(shù)據(jù)處理一般分為2種，一種是基于歐拉思想的火焰面平均，研究火焰面整體的變化，包括火焰面密度、火焰面曲率的PDF分布[3,10]、火焰面分形維數(shù)等。二是基于拉格朗日思想，把火焰面看做一系列包含溫度、速度等信息的點(diǎn)，研究離散點(diǎn)隨時(shí)間空間的變化[11-13]。
在湍流火焰中，火焰自身特性與不同湍流尺度復(fù)雜的耦合作用使得湍流火焰面發(fā)生褶皺，其中有一部分褶皺結(jié)構(gòu)很深入火焰已燃區(qū)或未燃區(qū)，這部分褶皺區(qū)域被稱作湍流火焰面上的關(guān)鍵褶皺結(jié)構(gòu)。它們是湍流火焰刷厚度增加的主要原因，也控制著湍流火焰速度等關(guān)鍵參數(shù)，是湍流與火焰相互作用的具體體現(xiàn)。根據(jù)Chen[14]等人的模擬結(jié)果，湍流火焰面“孤島”結(jié)構(gòu)的形成機(jī)理也與這些結(jié)構(gòu)相關(guān)。但是，傳統(tǒng)的曲率PDF分析方法并不能區(qū)分出火焰上這些關(guān)鍵褶皺區(qū)域。這是因?yàn)椋@些結(jié)構(gòu)很深入火焰的已燃區(qū)或未燃區(qū)，但由于其面積較大，導(dǎo)致單位長(zhǎng)度內(nèi)的弧度增量較小，故曲率并不一定大。另外，在前述2種思想的湍流火焰數(shù)據(jù)處理過(guò)程中，湍流火焰面上的空間相關(guān)信息也無(wú)法體現(xiàn)出來(lái)。
近幾年來(lái)，網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)已經(jīng)成功運(yùn)用于生物、物理等領(lǐng)域，分析復(fù)雜系統(tǒng)中節(jié)點(diǎn)的聯(lián)系規(guī)律和機(jī)理，其優(yōu)點(diǎn)是可以標(biāo)記出系統(tǒng)中具有某些特征的關(guān)鍵節(jié)點(diǎn)或結(jié)構(gòu)。在湍流方面，Liu[15]等人利用時(shí)間信息構(gòu)建的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)揭示了充分發(fā)展的三維湍流中存在的冪指規(guī)律，M. Meenatchidevi[16]等人利用時(shí)間網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)研究了熱聲震蕩的機(jī)理，分析了火焰由不震蕩到震蕩的關(guān)鍵臨界狀態(tài)，找到了熱聲震蕩的臨界點(diǎn)。
本文利用網(wǎng)絡(luò)結(jié)構(gòu)方法構(gòu)建湍流火焰面的拓?fù)浣Y(jié)構(gòu)，標(biāo)記對(duì)湍流火焰面有較大影響的關(guān)鍵褶皺結(jié)構(gòu)，分析湍流與火焰的相互作用規(guī)律。
1 實(shí)驗(yàn)裝置和方法
1.1 湍流實(shí)驗(yàn)裝置
實(shí)驗(yàn)裝置由供氣系統(tǒng)，OH-PLIF系統(tǒng)以及湍流本生燈系統(tǒng)組成。OH-PLIF系統(tǒng)包括Nd:YAG激光器(型號(hào)為Quanta-Ray Pro 190，頻率為10Hz，脈沖時(shí)間10ns，脈沖能量300mJ)、染料激光器以及片光透鏡組組成。供氣系統(tǒng)由氣瓶、混合室以及MKS質(zhì)量流量計(jì)(1179A，量程：0～50 SLM和1159A，量程：0～5 SLM)組成。關(guān)于OH-PLIF系統(tǒng)的詳細(xì)信息見(jiàn)文獻(xiàn)[17-18]。湍流本生燈剖面圖如圖1(a)所示，出口直徑D=20mm。A、B、C、D代表孔板與本生燈出口間的不同距離。湍流產(chǎn)生板的形狀如圖1(b)所示，S和P分別代表不同形狀的湍流產(chǎn)生板，利用不同的孔板放置位置及形狀，可以產(chǎn)生不同強(qiáng)度的湍流。
圖1 湍流燃燒實(shí)驗(yàn)裝置
湍流流場(chǎng)特征參數(shù)利用熱線風(fēng)速儀(Dantec，Streamline 90N)測(cè)量，并采用各向同性湍流假設(shè)計(jì)算獲得。在本生燈出口10mm處，沿著本生燈中心線各取4和8mm為測(cè)量點(diǎn)，分別測(cè)量這5個(gè)點(diǎn)當(dāng)?shù)氐耐牧鲝?qiáng)度、積分尺度、泰勒尺度以及柯氏尺度。然后進(jìn)行平均，從而得到冷態(tài)下整個(gè)流場(chǎng)的信息，如圖2所示。
圖2 本生燈湍流場(chǎng)測(cè)量示意圖
本文實(shí)驗(yàn)中，采用甲烷作為燃料，這是因?yàn)榧淄榈挠行ewis數(shù)接近于1，可以排除熱擴(kuò)散不穩(wěn)定性的影響。表1所示為不同工況下實(shí)驗(yàn)參數(shù)。實(shí)驗(yàn)中采用3種不同的孔板，S1、P2、P3，D表示孔板放置高度距離本生燈出口73.6mm，φ為當(dāng)量比，Uave為本生燈出口速度，m/s;u′為湍流強(qiáng)度，m/s；l0、l、h分別為積分尺度、泰勒尺度和柯氏尺度,mm。Li表示火焰自身不穩(wěn)定性引起的空間波長(zhǎng)[19]，計(jì)算過(guò)程如圖3(a)所示。
表1 甲烷實(shí)驗(yàn)工況Table 1 Experiment condition of methan flames
圖3(b)給出了實(shí)驗(yàn)工況在湍流火焰模式分區(qū)中的位置，可以看到，設(shè)計(jì)工況處于3個(gè)湍流火焰模式中：褶皺層流火焰面模式、波紋板式火焰面模式和薄反應(yīng)區(qū)。
(a) 不穩(wěn)定性波長(zhǎng)的計(jì)算
(b) 湍流火焰分布圖
Fig.3CalculationofDLinstabilitywavelengthandexperimentalconditionsinthecombustiondiagram
1.2 湍流火焰面提取
圖4所示為火焰面邊界提取過(guò)程。利用OH-PLIF裝置可以得到湍流火焰OH基的分布，其邊界為未燃區(qū)與已燃區(qū)邊界，被看做為湍流火焰面。把ICCD相機(jī)得到的原始OH-PLIF圖片剪裁后灰度化，再利用本課題組發(fā)展的自適應(yīng)閾值法可以精確地得到火焰面的位置信息[20]。經(jīng)過(guò)剪裁后的OH圖片大小為21.44mm×70.02mm，像素點(diǎn)為120×490。
圖4 火焰前鋒面結(jié)構(gòu)提取過(guò)程
2 網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建方法
2.1 網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建理論
對(duì)湍流火焰的關(guān)鍵褶皺結(jié)構(gòu)進(jìn)行仔細(xì)分析可以發(fā)現(xiàn)，相比湍流火焰面的非褶皺部分，它們可以看到更多的火焰面。因此，我們利用“visible[21]”方法構(gòu)建湍流火焰的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)。
要得到1個(gè)系統(tǒng)的拓?fù)鋱D，需要有節(jié)點(diǎn)、鏈路及映射函數(shù)3個(gè)信息。其中，湍流火焰拓?fù)浣Y(jié)構(gòu)的節(jié)點(diǎn)就是提取的火焰前鋒面上所有像素點(diǎn)。有了節(jié)點(diǎn)之后，對(duì)于1個(gè)節(jié)點(diǎn)，可以和其它節(jié)點(diǎn)之間連成直線。由于計(jì)算機(jī)中直線的計(jì)算是通過(guò)像素來(lái)表示的，因此，我們采用Bresenham算法[22]來(lái)計(jì)算2個(gè)節(jié)點(diǎn)之間的連線。
在湍流火焰中，湍流穿過(guò)高溫火焰面之后，湍流強(qiáng)度以及湍流渦會(huì)發(fā)生較大的變化。因此，為了表征不同空間尺寸的渦對(duì)臨近火焰面的影響，我們定義節(jié)點(diǎn)的關(guān)聯(lián)如下。同一個(gè)空間尺度下，在火焰面的未燃區(qū)內(nèi)，如果2個(gè)節(jié)點(diǎn)之間的連線和火焰面不相交，則這2個(gè)節(jié)點(diǎn)是“可見(jiàn)”的，在物理上，說(shuō)明未燃區(qū)內(nèi)這2個(gè)節(jié)點(diǎn)可能受到同一湍流渦的影響，如果2個(gè)節(jié)點(diǎn)之間是“可見(jiàn)”的，那么這2個(gè)節(jié)點(diǎn)存在空間上的聯(lián)系，這個(gè)就是火焰的映射函數(shù)。如圖5所示的節(jié)點(diǎn)N1、N2、N3、N4，在規(guī)定的空間尺度下，我們認(rèn)為N1和N4是“不可見(jiàn)”的，N1和N2、N3是“可見(jiàn)”的。
定義了湍流火焰拓?fù)浣Y(jié)構(gòu)的關(guān)系之后，對(duì)于任意湍流火焰面，可以得到其拓?fù)浣Y(jié)構(gòu)信息。此時(shí)，就可以得到整個(gè)湍流火焰節(jié)點(diǎn)的鄰接矩陣。拓?fù)浣Y(jié)構(gòu)的鄰接矩陣定義如下：
(1)
為了對(duì)得到的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)進(jìn)行分析，計(jì)算出每個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)度[23]k：
(2)
其中kn表示第n個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)度。節(jié)點(diǎn)度的概率密度為P(k)，其變化描述出網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)的動(dòng)態(tài)特性。由于Bresenham算法在相鄰節(jié)點(diǎn)上算法的不確定性，在后續(xù)分析中，節(jié)點(diǎn)度小于5的節(jié)點(diǎn)被舍去。
圖5 湍流火焰拓?fù)錁?gòu)建示意圖
2.2 一般曲線網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)
為了驗(yàn)證構(gòu)建方法的正確性，我們首先對(duì)一般的正弦曲線和直線進(jìn)行了網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建。其中，直線和正弦曲線的像素點(diǎn)及圖片大小和剪裁后的OH-PLIF圖片大小一致。圖6所示為一般曲線及其網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)在Gephi軟件中可視化結(jié)果。在Gephi可視化圖中，每個(gè)圓點(diǎn)對(duì)應(yīng)曲線中每個(gè)節(jié)點(diǎn)，其中圓的大小和顏色深度反應(yīng)原圖上每個(gè)像素節(jié)點(diǎn)度的大小。可以看出，在圖6(a)中，只有直線2個(gè)端點(diǎn)的節(jié)點(diǎn)度比較小，其它節(jié)點(diǎn)度一樣。根據(jù)圖6(b)和(c)可以看到，利用前述方法對(duì)曲線進(jìn)行網(wǎng)絡(luò)拓?fù)錁?gòu)建，節(jié)點(diǎn)度大的點(diǎn)可以標(biāo)記出曲線褶皺部分。
3 湍流火焰網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建
3.1 網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)與傳統(tǒng)PDF分析比較
在湍流火焰中，湍流與火焰的相互作用使得湍流火焰面有很多褶皺結(jié)構(gòu)。它們主要可以分為以下2種：第一種褶皺結(jié)構(gòu)具有較大的面積，很深入已燃?xì)饣蛭慈細(xì)猓鐖D7(a)中紅色實(shí)線方框所示，第二種褶皺結(jié)構(gòu)是湍流火焰面上微小突起，如圖7(b)中紅色虛線方框所示。但是，這2種褶皺結(jié)構(gòu)對(duì)整個(gè)湍流火焰的整體貢獻(xiàn)是不一樣的，第一種褶皺結(jié)構(gòu)對(duì)湍流火焰具有較大的影響，控制著湍流火焰刷厚度、湍流火焰體積、湍流火焰速度等參數(shù)，是湍流與火焰相互作用的具體體現(xiàn)[24-25]，是湍流火焰面上的關(guān)鍵褶皺區(qū)域。
(a) 直線以及網(wǎng)絡(luò)結(jié)構(gòu)節(jié)點(diǎn)度可視化
(b) 正弦曲線及網(wǎng)絡(luò)結(jié)構(gòu)節(jié)點(diǎn)度可視化
(c) 利用節(jié)點(diǎn)度重構(gòu)的正弦曲線
圖7 不同工況下湍流火焰OH-PLIF圖片
利用傳統(tǒng)的曲率PDF(Probability density function)方法無(wú)法準(zhǔn)確表征出不同工況湍流火焰面上這種關(guān)鍵褶皺區(qū)域的不同。這是因?yàn)?，這些關(guān)鍵的褶皺區(qū)域雖然很深入已燃?xì)怏w，但是，由于面積較大，導(dǎo)致單位長(zhǎng)度內(nèi)的弧度增量較小，曲率較小，比如圖7(b)中紅色虛線內(nèi)火焰褶皺結(jié)構(gòu)的曲率比圖7(a)中紅色實(shí)線的大。因此，利用曲率PDF分布無(wú)法區(qū)分出來(lái)圖7(a)，(b)所示不同工況下湍流火焰結(jié)構(gòu)的不同。
圖8所示為不同工況下湍流火焰曲率PDF分布。
圖8 不同工況下湍流火焰曲率PDF分布
對(duì)比圖7(a)～(f)的3組圖片，可以看出在不同的當(dāng)量比下，湍流火焰結(jié)構(gòu)有比較明顯的變化，但是曲率的PDF分布卻不能表征出這種變化。
由于曲率PDF無(wú)法反映不同火焰關(guān)鍵褶皺區(qū)域，因此，我們利用網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)來(lái)標(biāo)記湍流火焰面的關(guān)鍵褶皺區(qū)域。采用2.1節(jié)所述方法，可以對(duì)湍流火焰前鋒面進(jìn)行網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建。圖9所示為P2_D工況下，對(duì)不同時(shí)刻3張OH-PLIF圖片進(jìn)行網(wǎng)絡(luò)拓?fù)錁?gòu)建后的結(jié)果。圖9(a)～(c)為利用每個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)度重構(gòu)的湍流火焰面，(d)～(f)為對(duì)應(yīng)的網(wǎng)絡(luò)結(jié)構(gòu)在Gephi軟件中利用Fruchterman reingold[26]算法可視化結(jié)果。由圖9(a)～(c)可知，利用網(wǎng)絡(luò)結(jié)構(gòu)可以對(duì)湍流火焰每個(gè)區(qū)域內(nèi)節(jié)點(diǎn)進(jìn)行標(biāo)記，其中標(biāo)記出來(lái)具有較大節(jié)點(diǎn)度的節(jié)點(diǎn)位于湍流火焰前鋒面的褶皺區(qū)域，即圖中顏色較深區(qū)域的火焰結(jié)構(gòu)。仔細(xì)觀察可知，這些褶皺區(qū)域很深入地進(jìn)入到火焰的已燃區(qū)或者未燃區(qū)，也就是湍流火焰面的關(guān)鍵褶皺區(qū)域。圖9(d)～(f)中，每個(gè)圓圈代表其對(duì)應(yīng)火焰面圖片上的1個(gè)像素點(diǎn)，其中顏色較深、半徑較大的節(jié)點(diǎn)，對(duì)應(yīng)了(a)～(c)中顏色較深火焰結(jié)構(gòu)處的像素點(diǎn)，圖中的黑色細(xì)實(shí)線代表了各個(gè)像素點(diǎn)在網(wǎng)絡(luò)結(jié)構(gòu)間的聯(lián)系，即2個(gè)節(jié)點(diǎn)是“可見(jiàn)”的。
圖9 湍流火焰拓?fù)浣Y(jié)構(gòu)及可視化
3.2 不同湍流強(qiáng)度下火焰關(guān)鍵褶皺區(qū)域的形成原因
在網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)中，每個(gè)節(jié)點(diǎn)節(jié)點(diǎn)度的分布可以反映出系統(tǒng)的關(guān)鍵特性[27]，相似系統(tǒng)節(jié)點(diǎn)度的分布會(huì)展現(xiàn)出一定的相似性。湍流火焰的每1張OH-PLIF圖片都可以得到其對(duì)應(yīng)的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)及圖片上每個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)度，然后對(duì)50張OH-PLIF圖片進(jìn)行統(tǒng)計(jì)，得到所有節(jié)點(diǎn)度的概率分布，如圖10所示?？梢钥吹?，低湍流強(qiáng)度下，不同當(dāng)量比湍流火焰拓?fù)浣Y(jié)構(gòu)節(jié)點(diǎn)度的概率分布不同，當(dāng)量比f(wàn)=0.7的湍流火焰大節(jié)點(diǎn)度的概率比f(wàn)=1.0的高。這是因?yàn)楣?jié)點(diǎn)度大的點(diǎn)屬于湍流火焰面中關(guān)鍵的褶皺區(qū)域，而根據(jù)圖7的OH-PLIF圖片可以看到，低湍流強(qiáng)度下，相同湍流強(qiáng)度時(shí)，f=0.7的湍流火焰面上關(guān)鍵褶皺區(qū)域更多。對(duì)比同一當(dāng)量比下，P2_D和P3_D這2種不同湍流強(qiáng)度的湍流火焰，可以發(fā)現(xiàn)，其網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)節(jié)點(diǎn)度的分布基本相同。這說(shuō)明了在低湍流強(qiáng)度下，湍流火焰的關(guān)鍵褶皺區(qū)域是由于火焰的自身特性引起的，由于實(shí)驗(yàn)燃料甲烷的熱擴(kuò)散不穩(wěn)定性很小，那么此時(shí)引起湍流火焰面褶皺的主要因素是火焰自身的DL不穩(wěn)定性。
當(dāng)湍流強(qiáng)度增大到S1_D孔板時(shí)，可以看到，其網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)大節(jié)點(diǎn)度的概率更大，表明此時(shí)湍流火焰面關(guān)鍵褶皺區(qū)域更多，這個(gè)現(xiàn)象也可以從圖7中看出。在S1_D孔板產(chǎn)生的湍流中可以看到，當(dāng)量比為0.7和1.0的湍流火焰網(wǎng)絡(luò)結(jié)構(gòu)節(jié)點(diǎn)度概率分布基本一致。說(shuō)明在高湍流強(qiáng)度下，湍流火焰前鋒面關(guān)鍵褶皺區(qū)域的形成跟火焰自身特性無(wú)關(guān)，主要是受到湍流強(qiáng)度增加從而湍流渦尺度減小的影響，此時(shí)火焰自身不穩(wěn)定性對(duì)火焰褶皺影響較小。
(a) 低湍流強(qiáng)度下拓?fù)浣Y(jié)構(gòu)節(jié)點(diǎn)度概率分布
(b) 高湍流強(qiáng)度下拓?fù)浣Y(jié)構(gòu)節(jié)點(diǎn)度概率分布
3.3 火焰自身不穩(wěn)定性在湍流火焰空間發(fā)展
本生燈湍流火焰屬于附著在燃燒器出口的湍流火焰，其與球形火焰等自由傳播火焰最大的區(qū)別就是來(lái)流氣體在火焰刷的切向有分量[28]，從而使火焰面的褶皺部分向下游傳播。也就是說(shuō)，由于來(lái)流氣體的作用，下游火焰的火焰刷以及火焰面密度會(huì)受到上游火焰褶皺結(jié)構(gòu)影響，從而影響湍流局部消耗速率[7]，這也是湍流火焰的非局部效應(yīng)[29]。現(xiàn)階段，對(duì)于湍流燃燒局部消耗速率的大部分模型都是基于自由傳播的湍流火焰[30-31]，很少關(guān)注湍流在本生燈火焰的非局部效應(yīng)。對(duì)于DL不穩(wěn)定性在本生燈湍流火焰中發(fā)展的研究則更少。本文在低湍流強(qiáng)度下，利用標(biāo)記出來(lái)的關(guān)鍵褶皺區(qū)域可以對(duì)DL不穩(wěn)定性在本生燈湍流火焰中的發(fā)展規(guī)律進(jìn)行研究。由圖9(a)～(c)可以看到，利用網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)標(biāo)記出來(lái)的關(guān)鍵褶皺區(qū)域位于靠近本生燈的頂端處，而在底部，這種關(guān)鍵褶皺區(qū)域很少。對(duì)P3_D湍流下的湍流火焰進(jìn)行分析，我們可以得到類似的結(jié)論。這說(shuō)明在本生燈湍流火焰中，火焰由于自身不穩(wěn)定性因素引起的火焰面褶皺以及火焰刷厚度的增加是隨著空間發(fā)展的，越靠近本生燈火焰的頂端，這種影響越明顯。在本生燈湍流火焰中，對(duì)500張得到的OH-PLIF圖片進(jìn)行提取火焰前鋒面后疊加，可以得到火焰面進(jìn)展變量的分布，為了得到火焰前鋒面上關(guān)鍵褶皺區(qū)域?qū)ν牧骰鹧娴挠绊?，我們定義2個(gè)火焰刷厚度：水平火焰刷厚度dT,h以及中心線火焰刷厚度dT,0。根據(jù)不同的定義方法可以得到不同的火焰刷厚度，但是在相同的定義方式下，其變化規(guī)律相同[32]。本文定義中心線火焰刷厚度dT,0為本生燈中心線上進(jìn)展變量=0.1和0.5之間的厚度，水平火焰刷厚度dT,h為在本生燈出口高度一定距離h處，進(jìn)展變量=0.1和0.5火焰面之間水平方向的距離，如圖11(a)所示，圖11(b)和(c)所示為P2_D湍流強(qiáng)度下，f=0.7和1.0的湍流火焰刷，其中心線火焰刷厚度dT,0分別為20.50和13.63mm。這是因?yàn)?，該湍流?qiáng)度下，當(dāng)量比f(wàn)=0.7的湍流火焰的不穩(wěn)定性比1.0大，湍流火焰面上關(guān)鍵褶皺區(qū)域多，所以其頂端處中心線火焰刷厚度大。
圖11 不同工況下火焰刷厚度及定義
Fig.11Schematicofdefinitionofflamebrushthicknessandtwoflamebrushunderdifferentconditions
圖12所示為低湍流強(qiáng)度下，不同當(dāng)量比湍流火焰水平火焰刷厚度(dT,0/D)隨本生燈出口高度(h/D)無(wú)量綱后的變化。可以看到，在距離本生燈出口高度一定范圍內(nèi)，不同當(dāng)量比湍流火焰的水平火焰刷厚度基本一致，在無(wú)量綱出口距離達(dá)到一定高度處才表現(xiàn)出不同。說(shuō)明在本生燈湍流火焰中，火焰自身不穩(wěn)定性對(duì)湍流火焰的影響跟火焰自身的發(fā)展時(shí)間有關(guān)，在距離出口一定距離內(nèi)，火焰自身不穩(wěn)定對(duì)湍流火焰的影響體現(xiàn)不出。
圖12 水平火焰刷厚度隨出口高度的變化
Fig.12ThenormalisedhorizotalflamebrushthicknesswithrespecttothenormalisedaxialdistancefromBunsenburnerexit
4 結(jié) 論
(1) 利用網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)方法可以標(biāo)記出本生燈湍流火焰面的關(guān)鍵褶皺結(jié)構(gòu)。
(2) 在低湍流強(qiáng)度下，湍流火焰關(guān)鍵褶皺結(jié)構(gòu)主要由自身不穩(wěn)定性引起，在高湍流強(qiáng)度下，火焰關(guān)鍵褶皺結(jié)構(gòu)主要受湍流尺度影響。
(3) 在本生燈湍流火焰這種附著火焰中，火焰自身不穩(wěn)定性引起的火焰褶皺與火焰發(fā)展距離有關(guān)，在距本生燈出口一定距離處，火焰自身不穩(wěn)定性對(duì)火焰褶皺影響較小。隨著火焰向下游發(fā)展，其對(duì)火焰面影響逐漸增大，火焰褶皺程度增加。
[1]Peters N. Turbulent combustion[M]. Cambridge University Press, 2000.
[2]Tamadonfar Parsa, Gülder ?mer L. Flame brush characteristics and burning velocities of premixed turbulent methane/air Bunsen flames[J]. Combustion and Flame, 2014, 161(12): 3154-3165.
[3]Fragner R, Halter F, Mazellier N, et al. Investigation of pressure effects on the small scale wrinkling of turbulent premixed Bunsen flames[J]. Proceedings of the Combustion Institute, 2015, 35(2): 1527-1535.
[4]Poludnenko A Y, Oran E S. The interaction of high-speed turbulence with flames: Global properties and internal flame structure[J]. Combustion and Flame, 2010, 157(5): 995-1011.
[5]Lipatnikov A N, Chomiak J. Effects of premixed flames on turbulence and turbulent scalar transport[J]. Progress in Energy and Combustion Science, 2010, 36(1): 1-102.
[6]Nishiki S. Modeling of flame-generated turbulence based on direct numerical simulation databases[J]. Proceedings of the Combustion Institute, 2002, 29: 2017-2022.
[7]Lipatnikov A N, Chomiak J. Turbulent flame speed and thickness: phennomenology, evaluation, and application in multi-dimensional simulations[J]. Progress in Energy & Combustion Science, 2002, 28: 1-74.
[8]Fureby C. A fractal flame-wrinkling large eddy simulation model for premixed turbulent combustion[J]. Proceedings of the Combustion Institute, 2005, 30(1): 593-601.
[9]Cintosun Esen, Smallwood Gregory J, Gülder ?mer L. Flame surface fractal characteristics in premixed turbulent combustion at high turbulence intensities[J]. AIAA Journal, 2007, 45(11): 2785-2789.
[10]Bradley D. Application of a reynolds stress, stretched flamelet, mathematical model to computations of turbulent burning velocities andcomparison with experiments[J]. Combustion & Flame, 1994, 96: 221-248.
[11]Yeung P K. Lagrangian statistics from direct numerical simulations of isotropic turbulence[J]. Journal of Fluid Mechanics, 1989, 207: 531-586.
[12]Pope S. Lagrangian PDF methods for turbulent flows[J]. Annu Rev Fluid Mech, 1994, 26: 23-63.
[13]Chaudhuri S. Life of flame particles embedded in premixed flames interacting with near isotropic turbulence[J]. Proceedings of the Combustion Institute, 2015, 35(2): 1305-1312.
[14]Chen H J. The mechanism of two-dimensional pocket formation in lean premixed methane-air flames with implications to turbulent combustion[J]. Combustion & Flame, 1999: 15-48.
[15]Liu C, Zhou W X, Yuan W K. Statistical properties of visibility graph of energy dissipation rates in three-dimensional fully developed turbulence[J]. Physica A: Statistical Mechanics and its Applications, 2010, 389(13): 2675-2681.
[16]Murugesan Meenatchidevi, Sujith R I. Combustion noise is scale-free: transition from scale-free to order at the onset of thermoacoustic instability[J]. Journal of Fluid Mechanics, 2015, 772: 225-245.
[17]張猛, 王金華, 謝永亮, 等. 利用OH_PLIF測(cè)量CH4/H2/空氣混合氣湍流燃燒速率[J]. 燃燒科學(xué)與技術(shù), 2013, 19(6): 512-516.
Zhang M, Wang J H, Xie Y L, et al. Measurement of turbulent burning velocity of CH4/H2/Air mixtures using OH-PLIF[J]. Journal of Combustion Science and Technology, 2013, 19(6): 512-516.
[18]Zhang M, Wang J H, Wu J, et al. Flame front structure of turbulent premixed flames of syngas oxyfuel mixtures[J]. International Journal of Hydrogen Energy, 2014, 39(10): 5176-5185.
[19]Kobayashi H, Tamura T, Maruto K, et al. Burning velocity of turbulent premixed flames in a high pressure environment[J]. Proceedings of the Combustion Institute, 1996, 26(1): 389-396.
[20]張猛, 王金華, 俞森彬, 等. 自適應(yīng)閾值二值法提取湍流火焰前鋒面結(jié)構(gòu)[J]. 燃燒科學(xué)與技術(shù), 2016, 22(3): 212-217.
Zhang M, Wang J H, Yu S B, et al. Flame front tracking of turbulent premixed flames using adaptive threshold binarization[J]. Journal of Combustion Science and Technology, 2016, 22(3): 212-217.
[21]Luque B, Lacasa L, Ballesteros F, et al. Horizontal visibility graphs: Exact results for random time series[J]. Physical Review E Statistical Nonlinear & Soft Matter Physics, 2009, 80(2): 046103.
[22]Bresenham J E. Algorithm for computer control of a digital plotter[J]. IBM Systems Journal, 1965, 4(1):25-30.
[23]Barabasi A L, Albert R. Emergence of scaling in random networks[J]. Science, 1999, 286(5439): 509.
[24]Hamlington P E, Poludnenko A Y, Oran E S. Interactions between turbulence and flames in premixed reacting flows[J]. Physics of Fluids, 2011, 23(12): 125111.
[25]Chakraborty N, Klein M, Swaminathan N. Effects of lewis number on the reactive scalar gradient alignment with local strain rate in turbulent premixed flames[J]. Proceedings of the Combustion Institute, 2009, 32(1): 1409-1417.
[26]Fruchterman T M J, Reingold E M. Graph drawing by force-directed placement[J]. Software Practice & Experience, 2010, 21(11): 1129-1164.
[27]Scholz M. Node similarity as a basic principle behind connectivity in complex networks[J]. Computer Science, 2015: 1-7.
[28]Boyer L, Quinard J. On the dynamics of anchored flames[J]. Combustion & Flame, 1990, 82(1): 51-65.
[29]Lieuwen T. Local consumption speed of turbulent premixed flames-An analysis of “memory effect”[J]. Combustion & Flame, 2010, 157: 955-965.
[30]Clavin P, Williams F A. Theory of premixed-flame propagation in large-scale turbulence[J]. Journal of Fluid Mechanics, 2006, 90(3): 589-604.
[31]Aldredge R C, Williams F A. Influence of wrinkled premixed-flame dynamics on large-scale, low- intensity turbulent flow[J]. Journal of Fluid Mechanics, 2006, 228(228): 487-511.
[32]Tamadonfar Parsa, Gülder ?mer L. Effects of mixture composition and turbulence intensity on flame front structure and burning velocities of premixed turbulent hydrocarbon/air Bunsen flames[J]. Combustion and Flame, 2015, 162(12): 4417-4441.
猜你喜歡
基于通聯(lián)關(guān)系的通信網(wǎng)絡(luò)拓?fù)浒l(fā)現(xiàn)方法網(wǎng)絡(luò)安全與數(shù)據(jù)管理(2022年2期)2022-05-23 13:25:46 動(dòng)漫人物衣服褶皺的畫(huà)法（1）學(xué)生導(dǎo)報(bào)·中職周刊(2019年7期)2019-06-11 11:06:10 電子制作(2018年23期)2018-12-26 01:01:16 可壓縮Navier-Stokes方程平面Couette-Poiseuille流的線性不穩(wěn)定性數(shù)學(xué)物理學(xué)報(bào)(2018年2期)2018-05-14 07:32:05 一點(diǎn)褶皺Coco薇(2017年6期)2017-06-24 23:08:56 勞斯萊斯古斯特與魅影網(wǎng)絡(luò)拓?fù)鋱D汽車(chē)維修技師(2017年10期)2017-03-17 02:25:01 電測(cè)與儀表(2016年5期)2016-04-22 01:13:46 增強(qiáng)型體外反搏聯(lián)合中醫(yī)辯證治療不穩(wěn)定性心絞痛療效觀察中國(guó)繼續(xù)醫(yī)學(xué)教育(2015年1期)2016-01-06 01:36:15 復(fù)雜褶皺構(gòu)造分析與找煤江西煤炭科技(2015年2期)2015-11-07 03:10:25 夢(mèng)幻褶皺Coco薇(2015年7期)2015-08-13 22:33:38
王金華, 聶要輝, 常 敏, 張 猛, 黃佐華
(西安交通大學(xué) 能源與動(dòng)力工程學(xué)院, 西安 710065)
0 引 言
湍流火焰廣泛存在于燃?xì)廨啓C(jī)、航空發(fā)動(dòng)機(jī)、內(nèi)燃機(jī)、工業(yè)燃燒器等各種燃燒設(shè)備中，是寬廣時(shí)間空間尺度上湍流與火焰耦合作用的復(fù)雜物理化學(xué)問(wèn)題，其研究具有重要的工程和科學(xué)意義。根據(jù)Borghi[1]等人的湍流火焰模式分區(qū)理論，實(shí)際燃燒室湍流火焰主要位于火焰片模式及薄反應(yīng)區(qū)模式下，此時(shí)，湍流火焰的研究主要集中在湍流與局部層流火焰的相互作用。目前主要是從2個(gè)方面對(duì)此開(kāi)展實(shí)驗(yàn)研究:
一是從實(shí)驗(yàn)出發(fā)，研究湍流特性，包括湍流強(qiáng)度、積分尺度等對(duì)火焰的影響[2-4]，以及火焰化學(xué)反應(yīng)和熱釋放對(duì)湍流場(chǎng)的影響[5-6]。二是在數(shù)據(jù)處理方面，通過(guò)一定的數(shù)學(xué)統(tǒng)計(jì)分析方法對(duì)湍流火焰數(shù)據(jù)進(jìn)行提取和挖掘，找到湍流火焰中的特征參數(shù)，比如湍流火焰速度[7]、湍流火焰分形維數(shù)[8-9]等。湍流火焰的數(shù)據(jù)處理一般分為2種，一種是基于歐拉思想的火焰面平均，研究火焰面整體的變化，包括火焰面密度、火焰面曲率的PDF分布[3,10]、火焰面分形維數(shù)等。二是基于拉格朗日思想，把火焰面看做一系列包含溫度、速度等信息的點(diǎn)，研究離散點(diǎn)隨時(shí)間空間的變化[11-13]。
在湍流火焰中，火焰自身特性與不同湍流尺度復(fù)雜的耦合作用使得湍流火焰面發(fā)生褶皺，其中有一部分褶皺結(jié)構(gòu)很深入火焰已燃區(qū)或未燃區(qū)，這部分褶皺區(qū)域被稱作湍流火焰面上的關(guān)鍵褶皺結(jié)構(gòu)。它們是湍流火焰刷厚度增加的主要原因，也控制著湍流火焰速度等關(guān)鍵參數(shù)，是湍流與火焰相互作用的具體體現(xiàn)。根據(jù)Chen[14]等人的模擬結(jié)果，湍流火焰面“孤島”結(jié)構(gòu)的形成機(jī)理也與這些結(jié)構(gòu)相關(guān)。但是，傳統(tǒng)的曲率PDF分析方法并不能區(qū)分出火焰上這些關(guān)鍵褶皺區(qū)域。這是因?yàn)椋@些結(jié)構(gòu)很深入火焰的已燃區(qū)或未燃區(qū)，但由于其面積較大，導(dǎo)致單位長(zhǎng)度內(nèi)的弧度增量較小，故曲率并不一定大。另外，在前述2種思想的湍流火焰數(shù)據(jù)處理過(guò)程中，湍流火焰面上的空間相關(guān)信息也無(wú)法體現(xiàn)出來(lái)。
近幾年來(lái)，網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)已經(jīng)成功運(yùn)用于生物、物理等領(lǐng)域，分析復(fù)雜系統(tǒng)中節(jié)點(diǎn)的聯(lián)系規(guī)律和機(jī)理，其優(yōu)點(diǎn)是可以標(biāo)記出系統(tǒng)中具有某些特征的關(guān)鍵節(jié)點(diǎn)或結(jié)構(gòu)。在湍流方面，Liu[15]等人利用時(shí)間信息構(gòu)建的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)揭示了充分發(fā)展的三維湍流中存在的冪指規(guī)律，M. Meenatchidevi[16]等人利用時(shí)間網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)研究了熱聲震蕩的機(jī)理，分析了火焰由不震蕩到震蕩的關(guān)鍵臨界狀態(tài)，找到了熱聲震蕩的臨界點(diǎn)。
本文利用網(wǎng)絡(luò)結(jié)構(gòu)方法構(gòu)建湍流火焰面的拓?fù)浣Y(jié)構(gòu)，標(biāo)記對(duì)湍流火焰面有較大影響的關(guān)鍵褶皺結(jié)構(gòu)，分析湍流與火焰的相互作用規(guī)律。
1 實(shí)驗(yàn)裝置和方法
1.1 湍流實(shí)驗(yàn)裝置
實(shí)驗(yàn)裝置由供氣系統(tǒng)，OH-PLIF系統(tǒng)以及湍流本生燈系統(tǒng)組成。OH-PLIF系統(tǒng)包括Nd:YAG激光器(型號(hào)為Quanta-Ray Pro 190，頻率為10Hz，脈沖時(shí)間10ns，脈沖能量300mJ)、染料激光器以及片光透鏡組組成。供氣系統(tǒng)由氣瓶、混合室以及MKS質(zhì)量流量計(jì)(1179A，量程：0～50 SLM和1159A，量程：0～5 SLM)組成。關(guān)于OH-PLIF系統(tǒng)的詳細(xì)信息見(jiàn)文獻(xiàn)[17-18]。湍流本生燈剖面圖如圖1(a)所示，出口直徑D=20mm。A、B、C、D代表孔板與本生燈出口間的不同距離。湍流產(chǎn)生板的形狀如圖1(b)所示，S和P分別代表不同形狀的湍流產(chǎn)生板，利用不同的孔板放置位置及形狀，可以產(chǎn)生不同強(qiáng)度的湍流。
圖1 湍流燃燒實(shí)驗(yàn)裝置
湍流流場(chǎng)特征參數(shù)利用熱線風(fēng)速儀(Dantec，Streamline 90N)測(cè)量，并采用各向同性湍流假設(shè)計(jì)算獲得。在本生燈出口10mm處，沿著本生燈中心線各取4和8mm為測(cè)量點(diǎn)，分別測(cè)量這5個(gè)點(diǎn)當(dāng)?shù)氐耐牧鲝?qiáng)度、積分尺度、泰勒尺度以及柯氏尺度。然后進(jìn)行平均，從而得到冷態(tài)下整個(gè)流場(chǎng)的信息，如圖2所示。
圖2 本生燈湍流場(chǎng)測(cè)量示意圖
本文實(shí)驗(yàn)中，采用甲烷作為燃料，這是因?yàn)榧淄榈挠行ewis數(shù)接近于1，可以排除熱擴(kuò)散不穩(wěn)定性的影響。表1所示為不同工況下實(shí)驗(yàn)參數(shù)。實(shí)驗(yàn)中采用3種不同的孔板，S1、P2、P3，D表示孔板放置高度距離本生燈出口73.6mm，φ為當(dāng)量比，Uave為本生燈出口速度，m/s;u′為湍流強(qiáng)度，m/s；l0、l、h分別為積分尺度、泰勒尺度和柯氏尺度,mm。Li表示火焰自身不穩(wěn)定性引起的空間波長(zhǎng)[19]，計(jì)算過(guò)程如圖3(a)所示。
表1 甲烷實(shí)驗(yàn)工況Table 1 Experiment condition of methan flames
圖3(b)給出了實(shí)驗(yàn)工況在湍流火焰模式分區(qū)中的位置，可以看到，設(shè)計(jì)工況處于3個(gè)湍流火焰模式中：褶皺層流火焰面模式、波紋板式火焰面模式和薄反應(yīng)區(qū)。
(a) 不穩(wěn)定性波長(zhǎng)的計(jì)算
(b) 湍流火焰分布圖
Fig.3CalculationofDLinstabilitywavelengthandexperimentalconditionsinthecombustiondiagram
1.2 湍流火焰面提取
圖4所示為火焰面邊界提取過(guò)程。利用OH-PLIF裝置可以得到湍流火焰OH基的分布，其邊界為未燃區(qū)與已燃區(qū)邊界，被看做為湍流火焰面。把ICCD相機(jī)得到的原始OH-PLIF圖片剪裁后灰度化，再利用本課題組發(fā)展的自適應(yīng)閾值法可以精確地得到火焰面的位置信息[20]。經(jīng)過(guò)剪裁后的OH圖片大小為21.44mm×70.02mm，像素點(diǎn)為120×490。
圖4 火焰前鋒面結(jié)構(gòu)提取過(guò)程
2 網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建方法
2.1 網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建理論
對(duì)湍流火焰的關(guān)鍵褶皺結(jié)構(gòu)進(jìn)行仔細(xì)分析可以發(fā)現(xiàn)，相比湍流火焰面的非褶皺部分，它們可以看到更多的火焰面。因此，我們利用“visible[21]”方法構(gòu)建湍流火焰的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)。
要得到1個(gè)系統(tǒng)的拓?fù)鋱D，需要有節(jié)點(diǎn)、鏈路及映射函數(shù)3個(gè)信息。其中，湍流火焰拓?fù)浣Y(jié)構(gòu)的節(jié)點(diǎn)就是提取的火焰前鋒面上所有像素點(diǎn)。有了節(jié)點(diǎn)之后，對(duì)于1個(gè)節(jié)點(diǎn)，可以和其它節(jié)點(diǎn)之間連成直線。由于計(jì)算機(jī)中直線的計(jì)算是通過(guò)像素來(lái)表示的，因此，我們采用Bresenham算法[22]來(lái)計(jì)算2個(gè)節(jié)點(diǎn)之間的連線。
在湍流火焰中，湍流穿過(guò)高溫火焰面之后，湍流強(qiáng)度以及湍流渦會(huì)發(fā)生較大的變化。因此，為了表征不同空間尺寸的渦對(duì)臨近火焰面的影響，我們定義節(jié)點(diǎn)的關(guān)聯(lián)如下。同一個(gè)空間尺度下，在火焰面的未燃區(qū)內(nèi)，如果2個(gè)節(jié)點(diǎn)之間的連線和火焰面不相交，則這2個(gè)節(jié)點(diǎn)是“可見(jiàn)”的，在物理上，說(shuō)明未燃區(qū)內(nèi)這2個(gè)節(jié)點(diǎn)可能受到同一湍流渦的影響，如果2個(gè)節(jié)點(diǎn)之間是“可見(jiàn)”的，那么這2個(gè)節(jié)點(diǎn)存在空間上的聯(lián)系，這個(gè)就是火焰的映射函數(shù)。如圖5所示的節(jié)點(diǎn)N1、N2、N3、N4，在規(guī)定的空間尺度下，我們認(rèn)為N1和N4是“不可見(jiàn)”的，N1和N2、N3是“可見(jiàn)”的。
定義了湍流火焰拓?fù)浣Y(jié)構(gòu)的關(guān)系之后，對(duì)于任意湍流火焰面，可以得到其拓?fù)浣Y(jié)構(gòu)信息。此時(shí)，就可以得到整個(gè)湍流火焰節(jié)點(diǎn)的鄰接矩陣。拓?fù)浣Y(jié)構(gòu)的鄰接矩陣定義如下：
(1)
為了對(duì)得到的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)進(jìn)行分析，計(jì)算出每個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)度[23]k：
(2)
其中kn表示第n個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)度。節(jié)點(diǎn)度的概率密度為P(k)，其變化描述出網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)的動(dòng)態(tài)特性。由于Bresenham算法在相鄰節(jié)點(diǎn)上算法的不確定性，在后續(xù)分析中，節(jié)點(diǎn)度小于5的節(jié)點(diǎn)被舍去。
圖5 湍流火焰拓?fù)錁?gòu)建示意圖
2.2 一般曲線網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)
為了驗(yàn)證構(gòu)建方法的正確性，我們首先對(duì)一般的正弦曲線和直線進(jìn)行了網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建。其中，直線和正弦曲線的像素點(diǎn)及圖片大小和剪裁后的OH-PLIF圖片大小一致。圖6所示為一般曲線及其網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)在Gephi軟件中可視化結(jié)果。在Gephi可視化圖中，每個(gè)圓點(diǎn)對(duì)應(yīng)曲線中每個(gè)節(jié)點(diǎn)，其中圓的大小和顏色深度反應(yīng)原圖上每個(gè)像素節(jié)點(diǎn)度的大小。可以看出，在圖6(a)中，只有直線2個(gè)端點(diǎn)的節(jié)點(diǎn)度比較小，其它節(jié)點(diǎn)度一樣。根據(jù)圖6(b)和(c)可以看到，利用前述方法對(duì)曲線進(jìn)行網(wǎng)絡(luò)拓?fù)錁?gòu)建，節(jié)點(diǎn)度大的點(diǎn)可以標(biāo)記出曲線褶皺部分。
3 湍流火焰網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建
3.1 網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)與傳統(tǒng)PDF分析比較
在湍流火焰中，湍流與火焰的相互作用使得湍流火焰面有很多褶皺結(jié)構(gòu)。它們主要可以分為以下2種：第一種褶皺結(jié)構(gòu)具有較大的面積，很深入已燃?xì)饣蛭慈細(xì)猓鐖D7(a)中紅色實(shí)線方框所示，第二種褶皺結(jié)構(gòu)是湍流火焰面上微小突起，如圖7(b)中紅色虛線方框所示。但是，這2種褶皺結(jié)構(gòu)對(duì)整個(gè)湍流火焰的整體貢獻(xiàn)是不一樣的，第一種褶皺結(jié)構(gòu)對(duì)湍流火焰具有較大的影響，控制著湍流火焰刷厚度、湍流火焰體積、湍流火焰速度等參數(shù)，是湍流與火焰相互作用的具體體現(xiàn)[24-25]，是湍流火焰面上的關(guān)鍵褶皺區(qū)域。
(a) 直線以及網(wǎng)絡(luò)結(jié)構(gòu)節(jié)點(diǎn)度可視化
(b) 正弦曲線及網(wǎng)絡(luò)結(jié)構(gòu)節(jié)點(diǎn)度可視化
(c) 利用節(jié)點(diǎn)度重構(gòu)的正弦曲線
圖7 不同工況下湍流火焰OH-PLIF圖片
利用傳統(tǒng)的曲率PDF(Probability density function)方法無(wú)法準(zhǔn)確表征出不同工況湍流火焰面上這種關(guān)鍵褶皺區(qū)域的不同。這是因?yàn)?，這些關(guān)鍵的褶皺區(qū)域雖然很深入已燃?xì)怏w，但是，由于面積較大，導(dǎo)致單位長(zhǎng)度內(nèi)的弧度增量較小，曲率較小，比如圖7(b)中紅色虛線內(nèi)火焰褶皺結(jié)構(gòu)的曲率比圖7(a)中紅色實(shí)線的大。因此，利用曲率PDF分布無(wú)法區(qū)分出來(lái)圖7(a)，(b)所示不同工況下湍流火焰結(jié)構(gòu)的不同。
圖8所示為不同工況下湍流火焰曲率PDF分布。
圖8 不同工況下湍流火焰曲率PDF分布
對(duì)比圖7(a)～(f)的3組圖片，可以看出在不同的當(dāng)量比下，湍流火焰結(jié)構(gòu)有比較明顯的變化，但是曲率的PDF分布卻不能表征出這種變化。
由于曲率PDF無(wú)法反映不同火焰關(guān)鍵褶皺區(qū)域，因此，我們利用網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)來(lái)標(biāo)記湍流火焰面的關(guān)鍵褶皺區(qū)域。采用2.1節(jié)所述方法，可以對(duì)湍流火焰前鋒面進(jìn)行網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)構(gòu)建。圖9所示為P2_D工況下，對(duì)不同時(shí)刻3張OH-PLIF圖片進(jìn)行網(wǎng)絡(luò)拓?fù)錁?gòu)建后的結(jié)果。圖9(a)～(c)為利用每個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)度重構(gòu)的湍流火焰面，(d)～(f)為對(duì)應(yīng)的網(wǎng)絡(luò)結(jié)構(gòu)在Gephi軟件中利用Fruchterman reingold[26]算法可視化結(jié)果。由圖9(a)～(c)可知，利用網(wǎng)絡(luò)結(jié)構(gòu)可以對(duì)湍流火焰每個(gè)區(qū)域內(nèi)節(jié)點(diǎn)進(jìn)行標(biāo)記，其中標(biāo)記出來(lái)具有較大節(jié)點(diǎn)度的節(jié)點(diǎn)位于湍流火焰前鋒面的褶皺區(qū)域，即圖中顏色較深區(qū)域的火焰結(jié)構(gòu)。仔細(xì)觀察可知，這些褶皺區(qū)域很深入地進(jìn)入到火焰的已燃區(qū)或者未燃區(qū)，也就是湍流火焰面的關(guān)鍵褶皺區(qū)域。圖9(d)～(f)中，每個(gè)圓圈代表其對(duì)應(yīng)火焰面圖片上的1個(gè)像素點(diǎn)，其中顏色較深、半徑較大的節(jié)點(diǎn)，對(duì)應(yīng)了(a)～(c)中顏色較深火焰結(jié)構(gòu)處的像素點(diǎn)，圖中的黑色細(xì)實(shí)線代表了各個(gè)像素點(diǎn)在網(wǎng)絡(luò)結(jié)構(gòu)間的聯(lián)系，即2個(gè)節(jié)點(diǎn)是“可見(jiàn)”的。
圖9 湍流火焰拓?fù)浣Y(jié)構(gòu)及可視化
3.2 不同湍流強(qiáng)度下火焰關(guān)鍵褶皺區(qū)域的形成原因
在網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)中，每個(gè)節(jié)點(diǎn)節(jié)點(diǎn)度的分布可以反映出系統(tǒng)的關(guān)鍵特性[27]，相似系統(tǒng)節(jié)點(diǎn)度的分布會(huì)展現(xiàn)出一定的相似性。湍流火焰的每1張OH-PLIF圖片都可以得到其對(duì)應(yīng)的網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)及圖片上每個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)度，然后對(duì)50張OH-PLIF圖片進(jìn)行統(tǒng)計(jì)，得到所有節(jié)點(diǎn)度的概率分布，如圖10所示?？梢钥吹?，低湍流強(qiáng)度下，不同當(dāng)量比湍流火焰拓?fù)浣Y(jié)構(gòu)節(jié)點(diǎn)度的概率分布不同，當(dāng)量比f(wàn)=0.7的湍流火焰大節(jié)點(diǎn)度的概率比f(wàn)=1.0的高。這是因?yàn)楣?jié)點(diǎn)度大的點(diǎn)屬于湍流火焰面中關(guān)鍵的褶皺區(qū)域，而根據(jù)圖7的OH-PLIF圖片可以看到，低湍流強(qiáng)度下，相同湍流強(qiáng)度時(shí)，f=0.7的湍流火焰面上關(guān)鍵褶皺區(qū)域更多。對(duì)比同一當(dāng)量比下，P2_D和P3_D這2種不同湍流強(qiáng)度的湍流火焰，可以發(fā)現(xiàn)，其網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)節(jié)點(diǎn)度的分布基本相同。這說(shuō)明了在低湍流強(qiáng)度下，湍流火焰的關(guān)鍵褶皺區(qū)域是由于火焰的自身特性引起的，由于實(shí)驗(yàn)燃料甲烷的熱擴(kuò)散不穩(wěn)定性很小，那么此時(shí)引起湍流火焰面褶皺的主要因素是火焰自身的DL不穩(wěn)定性。
當(dāng)湍流強(qiáng)度增大到S1_D孔板時(shí)，可以看到，其網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)大節(jié)點(diǎn)度的概率更大，表明此時(shí)湍流火焰面關(guān)鍵褶皺區(qū)域更多，這個(gè)現(xiàn)象也可以從圖7中看出。在S1_D孔板產(chǎn)生的湍流中可以看到，當(dāng)量比為0.7和1.0的湍流火焰網(wǎng)絡(luò)結(jié)構(gòu)節(jié)點(diǎn)度概率分布基本一致。說(shuō)明在高湍流強(qiáng)度下，湍流火焰前鋒面關(guān)鍵褶皺區(qū)域的形成跟火焰自身特性無(wú)關(guān)，主要是受到湍流強(qiáng)度增加從而湍流渦尺度減小的影響，此時(shí)火焰自身不穩(wěn)定性對(duì)火焰褶皺影響較小。
(a) 低湍流強(qiáng)度下拓?fù)浣Y(jié)構(gòu)節(jié)點(diǎn)度概率分布
(b) 高湍流強(qiáng)度下拓?fù)浣Y(jié)構(gòu)節(jié)點(diǎn)度概率分布
3.3 火焰自身不穩(wěn)定性在湍流火焰空間發(fā)展
本生燈湍流火焰屬于附著在燃燒器出口的湍流火焰，其與球形火焰等自由傳播火焰最大的區(qū)別就是來(lái)流氣體在火焰刷的切向有分量[28]，從而使火焰面的褶皺部分向下游傳播。也就是說(shuō)，由于來(lái)流氣體的作用，下游火焰的火焰刷以及火焰面密度會(huì)受到上游火焰褶皺結(jié)構(gòu)影響，從而影響湍流局部消耗速率[7]，這也是湍流火焰的非局部效應(yīng)[29]。現(xiàn)階段，對(duì)于湍流燃燒局部消耗速率的大部分模型都是基于自由傳播的湍流火焰[30-31]，很少關(guān)注湍流在本生燈火焰的非局部效應(yīng)。對(duì)于DL不穩(wěn)定性在本生燈湍流火焰中發(fā)展的研究則更少。本文在低湍流強(qiáng)度下，利用標(biāo)記出來(lái)的關(guān)鍵褶皺區(qū)域可以對(duì)DL不穩(wěn)定性在本生燈湍流火焰中的發(fā)展規(guī)律進(jìn)行研究。由圖9(a)～(c)可以看到，利用網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)標(biāo)記出來(lái)的關(guān)鍵褶皺區(qū)域位于靠近本生燈的頂端處，而在底部，這種關(guān)鍵褶皺區(qū)域很少。對(duì)P3_D湍流下的湍流火焰進(jìn)行分析，我們可以得到類似的結(jié)論。這說(shuō)明在本生燈湍流火焰中，火焰由于自身不穩(wěn)定性因素引起的火焰面褶皺以及火焰刷厚度的增加是隨著空間發(fā)展的，越靠近本生燈火焰的頂端，這種影響越明顯。在本生燈湍流火焰中，對(duì)500張得到的OH-PLIF圖片進(jìn)行提取火焰前鋒面后疊加，可以得到火焰面進(jìn)展變量
圖11 不同工況下火焰刷厚度及定義
Fig.11Schematicofdefinitionofflamebrushthicknessandtwoflamebrushunderdifferentconditions
圖12所示為低湍流強(qiáng)度下，不同當(dāng)量比湍流火焰水平火焰刷厚度(dT,0/D)隨本生燈出口高度(h/D)無(wú)量綱后的變化。可以看到，在距離本生燈出口高度一定范圍內(nèi)，不同當(dāng)量比湍流火焰的水平火焰刷厚度基本一致，在無(wú)量綱出口距離達(dá)到一定高度處才表現(xiàn)出不同。說(shuō)明在本生燈湍流火焰中，火焰自身不穩(wěn)定性對(duì)湍流火焰的影響跟火焰自身的發(fā)展時(shí)間有關(guān)，在距離出口一定距離內(nèi)，火焰自身不穩(wěn)定對(duì)湍流火焰的影響體現(xiàn)不出。
圖12 水平火焰刷厚度隨出口高度的變化
Fig.12ThenormalisedhorizotalflamebrushthicknesswithrespecttothenormalisedaxialdistancefromBunsenburnerexit
4 結(jié) 論
(1) 利用網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)方法可以標(biāo)記出本生燈湍流火焰面的關(guān)鍵褶皺結(jié)構(gòu)。
(2) 在低湍流強(qiáng)度下，湍流火焰關(guān)鍵褶皺結(jié)構(gòu)主要由自身不穩(wěn)定性引起，在高湍流強(qiáng)度下，火焰關(guān)鍵褶皺結(jié)構(gòu)主要受湍流尺度影響。
(3) 在本生燈湍流火焰這種附著火焰中，火焰自身不穩(wěn)定性引起的火焰褶皺與火焰發(fā)展距離有關(guān)，在距本生燈出口一定距離處，火焰自身不穩(wěn)定性對(duì)火焰褶皺影響較小。隨著火焰向下游發(fā)展，其對(duì)火焰面影響逐漸增大，火焰褶皺程度增加。
[1]Peters N. Turbulent combustion[M]. Cambridge University Press, 2000.
[2]Tamadonfar Parsa, Gülder ?mer L. Flame brush characteristics and burning velocities of premixed turbulent methane/air Bunsen flames[J]. Combustion and Flame, 2014, 161(12): 3154-3165.
[3]Fragner R, Halter F, Mazellier N, et al. Investigation of pressure effects on the small scale wrinkling of turbulent premixed Bunsen flames[J]. Proceedings of the Combustion Institute, 2015, 35(2): 1527-1535.
[4]Poludnenko A Y, Oran E S. The interaction of high-speed turbulence with flames: Global properties and internal flame structure[J]. Combustion and Flame, 2010, 157(5): 995-1011.
[5]Lipatnikov A N, Chomiak J. Effects of premixed flames on turbulence and turbulent scalar transport[J]. Progress in Energy and Combustion Science, 2010, 36(1): 1-102.
[6]Nishiki S. Modeling of flame-generated turbulence based on direct numerical simulation databases[J]. Proceedings of the Combustion Institute, 2002, 29: 2017-2022.
[7]Lipatnikov A N, Chomiak J. Turbulent flame speed and thickness: phennomenology, evaluation, and application in multi-dimensional simulations[J]. Progress in Energy & Combustion Science, 2002, 28: 1-74.
[8]Fureby C. A fractal flame-wrinkling large eddy simulation model for premixed turbulent combustion[J]. Proceedings of the Combustion Institute, 2005, 30(1): 593-601.
[9]Cintosun Esen, Smallwood Gregory J, Gülder ?mer L. Flame surface fractal characteristics in premixed turbulent combustion at high turbulence intensities[J]. AIAA Journal, 2007, 45(11): 2785-2789.
[10]Bradley D. Application of a reynolds stress, stretched flamelet, mathematical model to computations of turbulent burning velocities andcomparison with experiments[J]. Combustion & Flame, 1994, 96: 221-248.
[11]Yeung P K. Lagrangian statistics from direct numerical simulations of isotropic turbulence[J]. Journal of Fluid Mechanics, 1989, 207: 531-586.
[12]Pope S. Lagrangian PDF methods for turbulent flows[J]. Annu Rev Fluid Mech, 1994, 26: 23-63.
[13]Chaudhuri S. Life of flame particles embedded in premixed flames interacting with near isotropic turbulence[J]. Proceedings of the Combustion Institute, 2015, 35(2): 1305-1312.
[14]Chen H J. The mechanism of two-dimensional pocket formation in lean premixed methane-air flames with implications to turbulent combustion[J]. Combustion & Flame, 1999: 15-48.
[15]Liu C, Zhou W X, Yuan W K. Statistical properties of visibility graph of energy dissipation rates in three-dimensional fully developed turbulence[J]. Physica A: Statistical Mechanics and its Applications, 2010, 389(13): 2675-2681.
[16]Murugesan Meenatchidevi, Sujith R I. Combustion noise is scale-free: transition from scale-free to order at the onset of thermoacoustic instability[J]. Journal of Fluid Mechanics, 2015, 772: 225-245.
[17]張猛, 王金華, 謝永亮, 等. 利用OH_PLIF測(cè)量CH4/H2/空氣混合氣湍流燃燒速率[J]. 燃燒科學(xué)與技術(shù), 2013, 19(6): 512-516.
Zhang M, Wang J H, Xie Y L, et al. Measurement of turbulent burning velocity of CH4/H2/Air mixtures using OH-PLIF[J]. Journal of Combustion Science and Technology, 2013, 19(6): 512-516.
[18]Zhang M, Wang J H, Wu J, et al. Flame front structure of turbulent premixed flames of syngas oxyfuel mixtures[J]. International Journal of Hydrogen Energy, 2014, 39(10): 5176-5185.
[19]Kobayashi H, Tamura T, Maruto K, et al. Burning velocity of turbulent premixed flames in a high pressure environment[J]. Proceedings of the Combustion Institute, 1996, 26(1): 389-396.
[20]張猛, 王金華, 俞森彬, 等. 自適應(yīng)閾值二值法提取湍流火焰前鋒面結(jié)構(gòu)[J]. 燃燒科學(xué)與技術(shù), 2016, 22(3): 212-217.
Zhang M, Wang J H, Yu S B, et al. Flame front tracking of turbulent premixed flames using adaptive threshold binarization[J]. Journal of Combustion Science and Technology, 2016, 22(3): 212-217.
[21]Luque B, Lacasa L, Ballesteros F, et al. Horizontal visibility graphs: Exact results for random time series[J]. Physical Review E Statistical Nonlinear & Soft Matter Physics, 2009, 80(2): 046103.
[22]Bresenham J E. Algorithm for computer control of a digital plotter[J]. IBM Systems Journal, 1965, 4(1):25-30.
[23]Barabasi A L, Albert R. Emergence of scaling in random networks[J]. Science, 1999, 286(5439): 509.
[24]Hamlington P E, Poludnenko A Y, Oran E S. Interactions between turbulence and flames in premixed reacting flows[J]. Physics of Fluids, 2011, 23(12): 125111.
[25]Chakraborty N, Klein M, Swaminathan N. Effects of lewis number on the reactive scalar gradient alignment with local strain rate in turbulent premixed flames[J]. Proceedings of the Combustion Institute, 2009, 32(1): 1409-1417.
[26]Fruchterman T M J, Reingold E M. Graph drawing by force-directed placement[J]. Software Practice & Experience, 2010, 21(11): 1129-1164.
[27]Scholz M. Node similarity as a basic principle behind connectivity in complex networks[J]. Computer Science, 2015: 1-7.
[28]Boyer L, Quinard J. On the dynamics of anchored flames[J]. Combustion & Flame, 1990, 82(1): 51-65.
[29]Lieuwen T. Local consumption speed of turbulent premixed flames-An analysis of “memory effect”[J]. Combustion & Flame, 2010, 157: 955-965.
[30]Clavin P, Williams F A. Theory of premixed-flame propagation in large-scale turbulence[J]. Journal of Fluid Mechanics, 2006, 90(3): 589-604.
[31]Aldredge R C, Williams F A. Influence of wrinkled premixed-flame dynamics on large-scale, low- intensity turbulent flow[J]. Journal of Fluid Mechanics, 2006, 228(228): 487-511.
[32]Tamadonfar Parsa, Gülder ?mer L. Effects of mixture composition and turbulence intensity on flame front structure and burning velocities of premixed turbulent hydrocarbon/air Bunsen flames[J]. Combustion and Flame, 2015, 162(12): 4417-4441.
猜你喜歡
網(wǎng)絡(luò)安全與數(shù)據(jù)管理(2022年2期)2022-05-23 13:25:46
學(xué)生導(dǎo)報(bào)·中職周刊(2019年7期)2019-06-11 11:06:10
數(shù)學(xué)物理學(xué)報(bào)(2018年2期)2018-05-14 07:32:05
Coco薇(2017年6期)2017-06-24 23:08:56
汽車(chē)維修技師(2017年10期)2017-03-17 02:25:01
中國(guó)繼續(xù)醫(yī)學(xué)教育(2015年1期)2016-01-06 01:36:15
江西煤炭科技(2015年2期)2015-11-07 03:10:25
Coco薇(2015年7期)2015-08-13 22:33:38
