張灝,高忠權(quán),高煦堯,寇夢(mèng)玄,杜增暉,洪流
(1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;2.液體火箭發(fā)動(dòng)機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室,710100,西安)
隨著環(huán)境問(wèn)題的日益嚴(yán)峻,清潔能源成為各國(guó)研究的重點(diǎn),天然氣作為一種重要的清潔燃料,主要成分為甲烷(CH4),具有安全性高、排放性能好、成本低的優(yōu)點(diǎn)[1]。為了提高利用清潔燃料的精細(xì)度,燃燒檢測(cè)技術(shù)成為檢測(cè)燃料燃燒狀況、控制燃燒工況的關(guān)鍵。實(shí)際應(yīng)用中,甲烷多以湍流狀態(tài)燃燒,國(guó)內(nèi)外學(xué)者利用燃燒檢測(cè)技術(shù)對(duì)甲烷的湍流燃燒特性進(jìn)行了大量研究。
Sweeney等利用分層渦流燃燒器控制分層和渦流程度,來(lái)對(duì)甲烷/空氣的預(yù)混湍流火焰結(jié)構(gòu)進(jìn)行研究[2-3]。Askari等通過(guò)向容彈內(nèi)噴射甲烷和空氣形成局部預(yù)混的湍流混合氣,改變溫度、壓強(qiáng)、添加氫氣或稀釋氣來(lái)研究火焰燃燒速率、壓強(qiáng)等參數(shù)[4-5]。王金華等分別在本生燈和風(fēng)扇擾動(dòng)的定容燃燒彈內(nèi)對(duì)甲烷/氫氣與空氣的湍流火焰進(jìn)行了火焰結(jié)構(gòu)、傳播速度的測(cè)量,并探究了稀釋氣、Lewis數(shù)等變量或參數(shù)對(duì)火焰的影響[6-8]。
目前,多數(shù)實(shí)驗(yàn)研究是利用紋影[4-5,8]、激光[2-3,6-7]等光學(xué)檢測(cè)技術(shù)對(duì)甲烷湍流火焰進(jìn)行測(cè)定,這種檢測(cè)手段多需要精密昂貴的設(shè)備,不利于實(shí)際應(yīng)用,而離子電流法因具有結(jié)構(gòu)簡(jiǎn)單、便于應(yīng)用[9]、采集到的信號(hào)豐富且時(shí)效性強(qiáng)等優(yōu)點(diǎn)被重視,近些年來(lái)得到了國(guó)內(nèi)外廣大學(xué)者的關(guān)注。
因離子電流法具有極高的實(shí)用性,所以多數(shù)學(xué)者將其應(yīng)用在發(fā)動(dòng)機(jī)氣缸內(nèi)燃燒的研究上[10-13]。但是,發(fā)動(dòng)機(jī)實(shí)驗(yàn)在火焰直觀的可視化方面太過(guò)局限,所以有部分學(xué)者將其與光學(xué)檢測(cè)法結(jié)合,來(lái)對(duì)火焰進(jìn)行研究。李理光等研究了電場(chǎng)對(duì)甲烷層流火焰離子電流的影響,并建立了數(shù)值模型[14-15]。吳筱敏等在定容燃燒彈上,結(jié)合紋影系統(tǒng)研究了電極大小和間隙等參數(shù)對(duì)甲烷層流火焰離子電流信號(hào)的影響[16-17]。Shy等通過(guò)設(shè)計(jì)一個(gè)十字型的定容燃燒裝置,利用離子電流信號(hào)對(duì)不同當(dāng)量比、不同湍流強(qiáng)度和不同稀釋條件下的甲烷/空氣預(yù)混湍流火焰的傳播速度進(jìn)行了測(cè)量[18-20]。
綜上,可以發(fā)現(xiàn)目前利用離子電流法對(duì)甲烷預(yù)混湍流火焰檢測(cè)的研究較少,更是缺少對(duì)甲烷預(yù)混湍流火焰離子電流信號(hào)的分析,尤其是結(jié)合火焰可視化圖像進(jìn)行分析。因此,本文采用離子電流法火焰檢測(cè)技術(shù),對(duì)定容燃燒彈內(nèi)初始溫度298 K,初始?jí)簭?qiáng)100 kPa,湍流裝置轉(zhuǎn)速系數(shù)kn為0(層流)和1 650(湍流),過(guò)量空氣系數(shù)λ為1.0、1.1、1.2、1.4和1.6時(shí)的甲烷/空氣預(yù)混湍流火焰的離子電流信號(hào)進(jìn)行采集,重點(diǎn)分析了甲烷預(yù)混湍流火焰離子電流信號(hào)的成因,以及不同工況下燃燒壓強(qiáng)信號(hào)和離子電流信號(hào)的峰值大小、峰值時(shí)刻和離子電流信號(hào)上升沿?cái)?shù)量的變化規(guī)律,并進(jìn)一步探究了離子電流與湍流火焰?zhèn)鞑ニ俣戎g的關(guān)系。
本文通過(guò)對(duì)賈偉東等所采用的定容燃燒彈實(shí)驗(yàn)臺(tái)架[21]進(jìn)行改進(jìn)以實(shí)現(xiàn)湍流場(chǎng)的產(chǎn)生,其示意圖如圖1所示。實(shí)驗(yàn)裝置由點(diǎn)火系統(tǒng)、紋影系統(tǒng)、進(jìn)排氣系統(tǒng)、數(shù)據(jù)采集系統(tǒng)和湍流定容燃燒彈5部分組成,其中湍流容彈和離子電流測(cè)量系統(tǒng)的結(jié)構(gòu)簡(jiǎn)圖如圖2所示。
圖1 離子電流測(cè)量湍流火焰實(shí)驗(yàn)臺(tái)架系統(tǒng)簡(jiǎn)圖 Fig.1 Diagram of turbulent flame experimental bench system for ion current measurement
圖2 湍流容彈和離子電流測(cè)量系統(tǒng)Fig.2 Turbulent constant-volume combustor and ion current measurement system
湍流容彈彈體為鐵制正方體裝置,整體接地,內(nèi)部為直徑130 mm、長(zhǎng)180 mm的水平圓柱式腔室,整個(gè)腔室除湍流發(fā)生裝置安裝孔處有部分裸露外,其余部分均覆有聚四氟乙烯。點(diǎn)火電極位于整個(gè)腔室中心上下兩側(cè),沿徑向?qū)ΨQ布置,放電尖端相距2 mm,除放電尖端裸露之外,其余部分均作絕緣處理。湍流發(fā)生裝置和點(diǎn)火電極位于同一橫截面,是由電機(jī)、扇葉和安裝套組成的旋轉(zhuǎn)擾動(dòng)式湍流發(fā)生裝置,其軸線水平布置,且與腔室軸線垂直,并距離40 mm,扇葉旋轉(zhuǎn)中心基本與腔室圓柱面重合。湍流的發(fā)生就由對(duì)稱安裝在容彈壁內(nèi)部的4組裝置產(chǎn)生,并通過(guò)調(diào)節(jié)器來(lái)控制電機(jī)轉(zhuǎn)速,調(diào)節(jié)器所顯示的數(shù)值定義為kn值。在開(kāi)放空間利用熱敏式風(fēng)速測(cè)速儀測(cè)得湍流發(fā)生裝置在kn=1 650時(shí)所產(chǎn)生的風(fēng)速vwind基本穩(wěn)定,且4個(gè)湍流裝置產(chǎn)生的風(fēng)速基本一致,如圖3所示,所以認(rèn)為能夠產(chǎn)生可控的相對(duì)穩(wěn)定的湍流場(chǎng)。
圖3 kn=1 650時(shí)各湍流發(fā)生裝置在開(kāi)放空間測(cè)得的風(fēng)速 Fig.3 The wind speed of various turbulence generating devices measured in open space when kn=1 650
離子電流測(cè)量系統(tǒng)由測(cè)量電路、測(cè)量電極和信號(hào)采集裝置3部分組成,其中測(cè)量電極采用Φ0.8 mm的不銹鋼絲制作,兩個(gè)測(cè)量電極對(duì)稱于容彈腔水平中心面布置,軸間距為3.8 mm,并和點(diǎn)火電極在同一橫截面,斷面距點(diǎn)火電極軸心5.5 mm,另一端與腔室圓柱面重合。通過(guò)高壓模塊將12 V直流電源升壓至80 V并加載到測(cè)量電極兩端,分壓電阻為900 kΩ,濾波電容為82 nF。
本文實(shí)驗(yàn)根據(jù)甲烷/空氣的化學(xué)計(jì)量燃燒反應(yīng)式,結(jié)合道爾頓分壓定律,通過(guò)壓力變送器將甲烷和空氣按照過(guò)量空氣系數(shù)為1.0、1.1、1.2、1.4和1.6充入容彈腔內(nèi)并靜置2 min,待其混合均勻之后啟動(dòng)湍流發(fā)生裝置,調(diào)節(jié)電機(jī)轉(zhuǎn)速控制器使kn分別為0和1 650,等待2 min,使湍流場(chǎng)平穩(wěn),通過(guò)點(diǎn)火系統(tǒng)點(diǎn)燃混合氣,數(shù)據(jù)采集系統(tǒng)采集離子電流信號(hào)和燃燒壓強(qiáng)信號(hào),紋影系統(tǒng)采集容彈內(nèi)混合氣燃燒圖像。實(shí)驗(yàn)中所用到的儀器儀表如表1所示。
表1 實(shí)驗(yàn)所用儀器儀表
基于離子電流法檢測(cè)火焰的基本原理,在甲烷/空氣的預(yù)混層流火焰(kn=0)中測(cè)得過(guò)量空氣系數(shù)λ=1.0時(shí)的燃燒壓強(qiáng)和離子電流信號(hào),以及用高速攝像機(jī)記錄的離子電流信號(hào)轉(zhuǎn)折時(shí)對(duì)應(yīng)的火焰燃燒狀態(tài),如圖4所示。
(a)不同時(shí)刻的火焰燃燒狀態(tài)
從圖4可以看出,燃燒壓強(qiáng)信號(hào)僅有一個(gè)上升沿,而離子電流信號(hào)包含3個(gè)上升階段,這3個(gè)階段分別對(duì)應(yīng)火花點(diǎn)火階段(t1~t2)、火焰前鋒區(qū)離子電流(t3~t4)和焰后區(qū)離子電流(t5~t6)。
火花點(diǎn)火階段(t1~t2)產(chǎn)生的離子電流信號(hào)主要是因點(diǎn)火時(shí)混合氣被擊穿放電產(chǎn)生的感應(yīng)電勢(shì)所引起的。測(cè)量電極與點(diǎn)火電極分離能夠有效避免點(diǎn)火時(shí)離子電流信號(hào)的產(chǎn)生[22],本實(shí)驗(yàn)為更全面地測(cè)得火焰發(fā)展過(guò)程的離子電流信號(hào),測(cè)量電極尖端距點(diǎn)火電極放電尖端較近,故造成在點(diǎn)火時(shí)產(chǎn)生了一個(gè)較小的離子電流信號(hào),如圖4中放大部分所示。
火焰前鋒離子電流(t3~t4)產(chǎn)生的離子電流是火焰前鋒面接觸測(cè)量電極產(chǎn)生的。甲烷燃燒時(shí),火焰前鋒面內(nèi)發(fā)生著復(fù)雜而劇烈的化學(xué)反應(yīng),其全局反應(yīng)可以描述為
CH4+2(O2+3.76N2)?CO2+2H2O+7.52N2
(1)
反應(yīng)過(guò)程會(huì)發(fā)生化學(xué)電離,產(chǎn)生大量帶電的離子基團(tuán),其中主要的化學(xué)電離過(guò)程為
(2)
CHO+基團(tuán)隨后與H2O產(chǎn)生迅速的質(zhì)子轉(zhuǎn)移反應(yīng)
(3)
k1、k2均為反應(yīng)速率常數(shù),其中k1=5×10-14cm3/s,k2=7×10-9cm3/s,k2在數(shù)量級(jí)上遠(yuǎn)大于k1,所以式(3)的反應(yīng)速率遠(yuǎn)大于式(2),H3O+的濃度遠(yuǎn)大于CHO+的濃度,故一般認(rèn)為,H3O+的數(shù)量在正離子中居首位。之后,H3O+會(huì)與電子發(fā)生中和反應(yīng)
(4)
所以火焰前鋒面中所生成的CHO+、H3O+、e-等帶電基團(tuán)在偏置電壓所生成的電場(chǎng)中受電場(chǎng)力的作用產(chǎn)生移動(dòng),進(jìn)而形成電流信號(hào),電流的大小可表示為
(5)
式中:If為火焰前鋒面的離子電流;Nion為單位體積的電荷密度;Vrz為測(cè)量電極間反應(yīng)區(qū)域的體積;E為電場(chǎng)強(qiáng)度;ve為遷移率;r為測(cè)量電極的間距。
焰后區(qū)離子電流(t5~t6)的形成主要是火焰前鋒面接觸容彈壁未絕緣部分形成的燃燒后期的火焰前鋒面離子電流信號(hào),疊加上焰后高溫區(qū)熱電離形成的離子電流信號(hào)[23]。此時(shí)的火焰前鋒面產(chǎn)生的離子電流信號(hào)同t3~t4階段產(chǎn)生機(jī)理一樣,高溫區(qū)熱電離目前普遍認(rèn)為主要是NO發(fā)生熱電離產(chǎn)生NO+和e-,這些經(jīng)熱電離產(chǎn)生的帶電離子基團(tuán)在電場(chǎng)力的作用下產(chǎn)生遷移,形成電流[24]。至于普遍認(rèn)為焰后區(qū)的大部分自由電子是由NO產(chǎn)生的,主要是因?yàn)槠錈犭x子化所需的熱能最低,離子化率又較高,如表2所示。
電流信號(hào)中有兩段下降沿,其中t2~t3是因?yàn)辄c(diǎn)火階段產(chǎn)生的離子電流消失所致;t4~t5的下降沿是因火焰前鋒面已經(jīng)不與測(cè)量電極接觸,但尚未與容彈壁的未絕緣部分接觸,且焰后區(qū)熱電離尚未開(kāi)始或電離程度不高,所以離子電流信號(hào)產(chǎn)生下降。
本實(shí)驗(yàn)采取湍流生成裝置在kn=1 650時(shí)所形成的湍流,此時(shí)湍流生成裝置在開(kāi)放環(huán)境下測(cè)得風(fēng)速vwind為21.22 m/s左右,測(cè)得λ=1.0時(shí)的燃燒壓強(qiáng)信號(hào)和離子電流信號(hào),以及用高速攝像機(jī)記錄的離子電流信號(hào)轉(zhuǎn)折時(shí)刻對(duì)應(yīng)的湍流火焰燃燒狀態(tài),如圖5所示。
表2 已燃高溫區(qū)主要成分的熱電離能和離子化率[25]
(a)不同時(shí)刻的火焰燃燒狀態(tài)
(b)燃燒壓強(qiáng)信號(hào)和離子電流信號(hào)圖5 λ=1.0、kn=1 650時(shí)甲烷/空氣預(yù)混湍流火焰的燃燒壓強(qiáng)信號(hào)和離子電流信號(hào)Fig.5 The combustion pressure signal and the ion current signal of methane/air premixed turbulence flame when λ=1.0 and kn=1 650
從圖中可以看出,和層流火焰采集到的信號(hào)相比,燃燒壓強(qiáng)信號(hào)依舊只有一個(gè)上升沿,而離子電流信號(hào)有4段上升期,其中:t1~t2對(duì)應(yīng)點(diǎn)火階段的離子電流信號(hào);t3~t4對(duì)應(yīng)火焰前鋒面與測(cè)量電極接觸產(chǎn)生的離子電流信號(hào);t5~t6階段的上升沿是由火焰前鋒面觸壁后在湍流作用下折返與測(cè)量電極接觸產(chǎn)生的離子電流、火焰前鋒面與容彈壁未絕緣部分接觸產(chǎn)生的離子電流以及焰后區(qū)熱電離產(chǎn)生的離子電流的疊加信號(hào);t7~t8對(duì)應(yīng)火焰完全燃燒后高溫區(qū)NO熱電離產(chǎn)生的離子電流信號(hào)。
離子電流信號(hào)在t1~t7之間存在3段下降沿,分別為t2~t3、t4~t5和t6~t7。其中:t2~t3是因?yàn)辄c(diǎn)火階段產(chǎn)生的離子電流消失所致;t4~t5為火焰前鋒面與測(cè)量電極分離之后又尚未觸壁折返,且未接觸容彈壁導(dǎo)電部分,焰后區(qū)熱電離程度過(guò)低或未產(chǎn)生熱電離而產(chǎn)生的離子電流信號(hào)下降;t6~t7是因火焰燃燒完畢之后,其折返火焰前鋒面與離子電流測(cè)量電極接觸產(chǎn)生的離子電流和火焰前鋒面接觸容彈壁導(dǎo)電部分產(chǎn)生的離子電流信號(hào)消失,且消失速率大于焰后區(qū)NO熱電離產(chǎn)生離子電流信號(hào)的速率而產(chǎn)生的下降。
在層流(kn=0)和湍流(kn=1 650)的狀態(tài)下,分別測(cè)得過(guò)量空氣系數(shù)λ為1.0、1.2、1.4和1.6時(shí)的燃燒壓強(qiáng)和離子電流信號(hào),其對(duì)比圖如圖6所示。
(a)λ=1.0
(b)λ=1.2
(c)λ=1.4
(d)λ=1.6圖6 不同λ下甲烷/空氣預(yù)混層流和湍流火焰的燃燒壓強(qiáng)信號(hào)和離子電流信號(hào)對(duì)比Fig.6 ContrastFigure of methane/air premixed laminar flame’s combustion pressure signal and ion current signal and methane/air premixed turbulence flame’s when λ=1.0,1.2,1.4 and 1.6 respectively
圖6以點(diǎn)火時(shí)刻為基準(zhǔn),即點(diǎn)火時(shí)刻為0 ms開(kāi)始記錄燃燒壓強(qiáng)信號(hào)和離子電流信號(hào);離子電流信號(hào)忽略點(diǎn)火感應(yīng)電勢(shì)造成的離子電流信號(hào),即從0 μA開(kāi)始記錄因火焰造成的正向離子電流信號(hào)。
從峰值角度分析,在相同λ下,湍流的燃燒壓強(qiáng)和離子電流信號(hào)均比層流的大,而隨著λ增大,不管是湍流火焰還是層流火焰的各信號(hào)峰值均隨之減小。
從時(shí)間角度分析,在相同λ下,湍流火焰各信號(hào)到達(dá)對(duì)應(yīng)峰值的時(shí)刻均比層流火焰提前,當(dāng)λ不同時(shí),在同一傳播階段,層、湍流火焰的各測(cè)量信號(hào)到達(dá)此階段峰值的時(shí)間均隨著λ的增大而變長(zhǎng)。
從上升沿?cái)?shù)量分析,在相同λ和不同λ下,層流和湍流的燃燒壓強(qiáng)曲線均為單峰曲線,僅一個(gè)上升沿。離子電流曲線的上升沿?cái)?shù)量在相同λ下,湍流的上升沿?cái)?shù)量總是大于或等于層流的數(shù)量,而在不同λ下,層流和湍流離子電流信號(hào)的上升沿?cái)?shù)量均隨λ增大而減少,在λ=1.2時(shí),層流和湍流離子電流信號(hào)的焰后區(qū)部分均消失,但湍流還保留火焰觸壁后上升的離子電流信號(hào),在λ=1.4,1.6時(shí),此部分信號(hào)也消失,只剩火焰未觸壁前的離子電流信號(hào)。
不論從峰值大小、到達(dá)峰值時(shí)刻,還是離子電流信號(hào)的上升沿?cái)?shù)量分析,均說(shuō)明在相同λ下,一定程度的湍流對(duì)火焰的燃燒具有促進(jìn)作用。結(jié)合圖4和圖5中的火焰圖像可以發(fā)現(xiàn),相比層流火焰,湍流加強(qiáng)了火焰的不穩(wěn)定性,使火焰面出現(xiàn)大量的破裂和褶皺,形成規(guī)格小但數(shù)量多的火焰胞面,這些胞面擴(kuò)大了火焰前鋒面與未燃區(qū)混合氣的接觸面積,從而提高了燃燒速度。
火焰?zhèn)鞑ニ俣仁侨剂先紵龝r(shí)的重要參數(shù),其可以定義為未燃區(qū)混合氣沿火焰面法線方向單位時(shí)間內(nèi)向火焰面移動(dòng)的距離。層流火焰的局部火焰面可以看作一個(gè)平面,其法線方向固定,所以可以用單位時(shí)間內(nèi)局部火焰前鋒面的傳播距離定義其速度;湍流火焰因火焰面的形狀隨時(shí)都在變化,很難確定一個(gè)固定的法線方向,所以湍流火焰的傳播速度以單位時(shí)間內(nèi)火焰投影面積的變化來(lái)定義,如下式所示
(6)
式中:vcom為湍流火焰燃燒速度;St、St+Δt為t、t+Δt時(shí)刻的湍流火焰面積,火焰面積選用以兩點(diǎn)火電極中點(diǎn)為圓心、半徑為57 mm的圓形范圍內(nèi)的火焰圖像提取的面積,以避免扇葉對(duì)火焰圖像處理的影響;Δt為高速攝像機(jī)的采樣間隔0.2 ms。
火焰圖像的處理基于matlab R2020a的Image Processing Toolbox,采用形態(tài)法濾波進(jìn)行邊緣提取。將提取的火焰邊緣(黑色邊線)與火焰原圖進(jìn)行疊加,如圖7所示,可以發(fā)現(xiàn)提取的邊緣與實(shí)際邊緣重合。
圖7 湍流火焰提取邊緣與實(shí)際邊緣對(duì)比Fig.7 Contrast of turbulent flame’s extracted edge and its actual edge
為測(cè)得kn=1 650時(shí)湍流火焰離子電流和燃燒速度之間的關(guān)系,選取λ為1.0、1.1、1.2時(shí),火焰開(kāi)始接觸測(cè)量電極時(shí)刻至火焰任意徑向半徑達(dá)到57 mm時(shí)刻之間的離子電流及對(duì)應(yīng)的燃燒速度進(jìn)行分析,結(jié)果如圖8所示。
(a)λ=1.0
(b)λ=1.1
(c)λ=1.2 圖8 不同λ下湍流火焰的離子電流和燃燒速度 之間的關(guān)系 Fig.8 The relationship between turbulent flame’s ion current and its combustion velocity when kn=1 650,λ=1.0,1.1 and 1.2 respectively
擬合曲線的相關(guān)系數(shù)平方R2在λ=1.0時(shí)為0.977 6,λ=1.1時(shí)為0.938 3,λ=1.2時(shí)為0.831 9,3種工況下的R2均大于0.75,說(shuō)明在這3種工況下湍流火焰的離子電流和燃燒速度成正線性相關(guān)。
本實(shí)驗(yàn)通過(guò)對(duì)在kn=1 650下的甲烷/空氣預(yù)混湍流火焰燃燒全過(guò)程的火焰圖像、燃燒壓強(qiáng)信號(hào)和離子電流信號(hào)進(jìn)行分析,得出以下結(jié)論。
(1)當(dāng)λ=1時(shí),與層流相比,湍流預(yù)混火焰?zhèn)鞑ミ^(guò)程的離子電流信號(hào)存在更多上升沿,會(huì)出現(xiàn)一段火焰前鋒面與焰后區(qū)共同作用的離子電流信號(hào)。
(2)通過(guò)對(duì)λ分別為1.0、1.2、1.4和1.6時(shí)的層流和湍流火焰的燃燒壓強(qiáng)信號(hào)和離子電流信號(hào)進(jìn)行對(duì)比,發(fā)現(xiàn)隨著λ減小,層、湍流預(yù)混火焰的燃燒壓強(qiáng)和離子電流信號(hào)變化趨勢(shì)一樣:各傳播階段峰值變大,峰值時(shí)刻提前且離子電流信號(hào)上升沿?cái)?shù)量增加;在相同λ下,湍流對(duì)火焰的傳播有促進(jìn)作用,結(jié)合火焰圖像可以得出湍流增大了火焰前鋒面與未燃混合氣的接觸面積,從而加速燃燒,使壓強(qiáng)和離子電流信號(hào)的峰值大于層流,峰值時(shí)刻早于層流,離子電流上升沿?cái)?shù)量大于或等于層流。
(3)通過(guò)對(duì)kn=1 650,λ為1.0、1.1和1.2時(shí)的甲烷預(yù)混湍流火焰的離子電流和燃燒速度進(jìn)行相關(guān)性分析,發(fā)現(xiàn)離子電流和燃燒速度之間存在正線性相關(guān),R2分別為0.977 6、0.938 3和0.831 9。