国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

基于內(nèi)窺火焰?zhèn)鞲衅骷夹g(shù)的超聲速燃燒感知實(shí)驗(yàn)研究

2022-05-21 02:19李忠朋周芮旭孟凡釗陳池李拓連歡
實(shí)驗(yàn)流體力學(xué) 2022年2期

李忠朋 ,周芮旭 ,孟凡釗 ,陳池 ,李拓 ,連歡

1.中國科學(xué)院力學(xué)研究所 高溫氣體動力學(xué)國家重點(diǎn)實(shí)驗(yàn)室,北京 100190

2.中國科學(xué)院大學(xué) 工程科學(xué)學(xué)院,北京 100190

0 引 言

未來寬域飛行包線內(nèi)提升沖壓發(fā)動機(jī)性能需要依靠先進(jìn)的燃燒組織策略及發(fā)動機(jī)系統(tǒng)控制技術(shù),其中高動態(tài)頻響傳感器及作動機(jī)構(gòu)是高性能控制系統(tǒng)(Full Authority Digital Engine Controller,F(xiàn)ADEC)的關(guān)鍵技術(shù)之一。穩(wěn)定可靠、結(jié)構(gòu)簡單、非侵入、能耗低、重量輕是沖壓發(fā)動機(jī)在線傳感器的基本要求。

壓力是表征超燃沖壓發(fā)動機(jī)內(nèi)流道流動特性及燃燒釋熱的重要物理量,且高速來流環(huán)境下壓力傳感器測量數(shù)據(jù)穩(wěn)定可靠,被廣泛應(yīng)用于沖壓發(fā)動機(jī)實(shí)驗(yàn)。壓力傳感器通常裝配在發(fā)動機(jī)壁面,激波運(yùn)動和燃燒釋熱引起的壓力變化反饋至發(fā)動機(jī)壁面具有一定延遲,因此壓力傳感器用于燃燒過程感知存在時(shí)間滯后效應(yīng)。

可調(diào)諧半導(dǎo)體激光吸收光譜技術(shù)(Tunable Diode Laser Absorption Spectroscopy,TDLAS )通過電流和溫度調(diào)諧半導(dǎo)體激光器的輸出波長,對被測物質(zhì)進(jìn)行掃描。由于被測氣體具有選擇性吸收的特性,基于Beer-Lambert定律,通過檢測吸收光譜的吸收強(qiáng)度,TDLAS可直接獲得燃燒產(chǎn)物的組分濃度、溫度等物理量,測量信息豐富,在沖壓發(fā)動機(jī)燃燒過程感知方面有獨(dú)特優(yōu)勢。然而該技術(shù)需要同時(shí)安裝發(fā)射端和接收端,增加了發(fā)動機(jī)燃燒室在線感知系統(tǒng)應(yīng)用的復(fù)雜程度。目前工業(yè)動力裝備的應(yīng)用主要包括:通過測量燃?xì)鉁囟群退魵鉂舛日{(diào)節(jié)空氣的流量優(yōu)化燃料的混合程度以提高燃燒效率;在基于甲烷的燃?xì)獍l(fā)電裝置中,通過監(jiān)測火爐內(nèi)的濕度和未燃燃?xì)鉂舛染S持系統(tǒng)的安全運(yùn)轉(zhuǎn);監(jiān)測燃?xì)廨啓C(jī)中氧氣的流量等。

此外,309 nm峰值紫外波段OH*與431 nm峰值藍(lán)光波段CH(A)*的比值OH*/CH(A)*以及天鵝波段C*與431 nm峰值藍(lán)光波段CH(A)*的比值C*/CH(A)*均與碳?xì)淙剂袭?dāng)量比φ有較好的線性正相關(guān)關(guān)系,且受湍流火焰應(yīng)力拉伸的影響不大,因此通常使用OH*/CH(A)*以及C*/CH(A)*表征碳?xì)淙剂袭?dāng)量比。對于甲烷或天然氣火焰,使用OH*/CH(A)*表征當(dāng)量比較為普遍;對于較大的碳?xì)淙剂系幕鹧?,選擇C*/CH(A)*則更為合適。此外,由于紫外波段被動自發(fā)光強(qiáng)度較弱,且紫外波段光電傳感器的性能通常低于天鵝波段光電傳感器,因此本文在基于乙烯的超燃環(huán)境下選擇C*/CH(A)*表征碳?xì)淙剂袭?dāng)量比。

光纖火焰?zhèn)鞲衅骷夹g(shù)基于火焰自發(fā)光原理,直接測量燃燒過程火焰自發(fā)光強(qiáng)度。被動燃燒化學(xué)自發(fā)光信號蘊(yùn)含了豐富的燃燒特性,通過建立特定波長化學(xué)發(fā)光強(qiáng)度與燃燒特性及發(fā)動機(jī)總體性能的關(guān)聯(lián)來感知沖壓發(fā)動機(jī)燃燒過程。該技術(shù)僅具有接收端,結(jié)構(gòu)簡單、惡劣環(huán)境適應(yīng)性強(qiáng)是該技術(shù)在發(fā)動機(jī)燃燒過程感知中的主要優(yōu)勢。

光纖火焰?zhèn)鞲衅鲝V泛應(yīng)用于內(nèi)燃機(jī)及航空發(fā)動機(jī)燃燒過程的在線監(jiān)測,高溫光纖火焰?zhèn)鞲衅髟O(shè)計(jì)技術(shù)主要由美國圣地亞國家實(shí)驗(yàn)室、奧地利AVL公司掌握。在超聲速燃燒領(lǐng)域,光纖傳感器的使用較為少見,需研究人員根據(jù)沖壓發(fā)動機(jī)工作特性進(jìn)行針對性開發(fā)。

2014年由美國空軍企業(yè)創(chuàng)新SBIR (Small Business Innovation Research)項(xiàng)目“Passive Optical Combustion Sensors for Scramjet Engine Control”支持,對沖壓發(fā)動機(jī)穩(wěn)焰模式研究做出了突出貢獻(xiàn)的Micka博士作為項(xiàng)目負(fù)責(zé)人,進(jìn)行了被動光譜光纖火焰?zhèn)鞲衅髟跊_壓發(fā)動機(jī)燃燒室的首次應(yīng)用嘗試。項(xiàng)目初步驗(yàn)證了光纖火焰?zhèn)鞲衅髟跊_壓發(fā)動機(jī)燃燒室感知方面的應(yīng)用前景。Micka通過在發(fā)動機(jī)壁面多點(diǎn)布設(shè)光纖傳感器測點(diǎn),研究了先鋒燃料關(guān)閉前后局部當(dāng)量比的變化,并基于局部當(dāng)量比預(yù)測了閃回現(xiàn)象的發(fā)生,并分析了超聲速燃燒的振蕩特性。

針對發(fā)動機(jī)的主動控制和燃燒感知需求,本文根據(jù)沖壓發(fā)動機(jī)工作特性,提出一種創(chuàng)新的光纖火焰?zhèn)鞲衅髋c沖壓發(fā)動機(jī)點(diǎn)火器一體集成化的傳感器設(shè)計(jì)方法,不需要單獨(dú)開設(shè)測量孔及光學(xué)窗口,可廣泛應(yīng)用于圓形燃燒室。下面將開展該傳感器感知沖壓發(fā)動機(jī)超聲速燃燒的實(shí)驗(yàn)研究,討論發(fā)動機(jī)超燃和亞燃工作模態(tài)的釋熱分布、局部當(dāng)量比特征和燃燒振蕩特性,基于火焰質(zhì)心的感知方法,驗(yàn)證內(nèi)窺火焰?zhèn)鞲衅骷夹g(shù)的超聲速燃燒感知價(jià)值。

1 實(shí)驗(yàn)設(shè)備與實(shí)驗(yàn)參數(shù)

1.1 直連實(shí)驗(yàn)設(shè)備及參數(shù)

中國科學(xué)院力學(xué)研究所研制的連續(xù)變馬赫數(shù)直連實(shí)驗(yàn)臺主體由流量動態(tài)調(diào)節(jié)氣源、可調(diào)加熱器、可調(diào)噴管以及實(shí)驗(yàn)段組成,如圖1所示。通過壓力調(diào)節(jié)閥實(shí)現(xiàn)空氣、氫氣、氧氣流量的動態(tài)調(diào)節(jié);采用環(huán)形多孔噴注中心燃燒的燒氫補(bǔ)氧方式對來流進(jìn)行加熱;通過氣源流量動態(tài)調(diào)節(jié)實(shí)現(xiàn)加熱器內(nèi)模擬總溫總壓實(shí)時(shí)可調(diào);可調(diào)噴管型面按照高低馬赫數(shù)雙設(shè)計(jì)點(diǎn)非對稱方法設(shè)計(jì),由機(jī)電伺服機(jī)構(gòu)驅(qū)動改變喉道面積實(shí)現(xiàn)模擬飛行馬赫數(shù)實(shí)時(shí)調(diào)節(jié)??烧{(diào)加熱器和可調(diào)噴管均采用水冷熱防護(hù)。流量動態(tài)調(diào)節(jié)氣源、可調(diào)加熱器、可調(diào)噴管通過PLC工業(yè)控制系統(tǒng)進(jìn)行精確時(shí)序控制,接收實(shí)驗(yàn)控制臺操作指令,并根據(jù)指令執(zhí)行流量調(diào)節(jié)、伺服機(jī)構(gòu)運(yùn)動等相應(yīng)動作,同時(shí)具有實(shí)驗(yàn)臺主體氣源壓力、直連臺雙螺柱安裝型推力傳感器數(shù)據(jù)采集和處理能力。實(shí)驗(yàn)臺可模擬最大氣體流量2.5 kg/s,最高總溫1 900 K,最高總壓4 Mpa,最長工作時(shí)間30 s,飛行馬赫數(shù)范圍4.5~6.5。

圖1 超燃沖壓發(fā)動機(jī)直連式實(shí)驗(yàn)臺示意圖[28]Fig.1 Schematic diagram of direct-connected experimental bench for scramjet[28]

實(shí)驗(yàn)?zāi)M燃燒室入口馬赫數(shù)2.8、飛行馬赫數(shù)5.6、總溫1 475 K、總壓1 678 kPa、加熱氣體總流量1 379 g/s、實(shí)驗(yàn)有效時(shí)間2 s。實(shí)驗(yàn)段為80 mm×40 mm的等截面隔離段和單邊擴(kuò)張雙凹腔燃燒室模型。隔離段長度為337 mm,燃燒室模型下壁面擴(kuò)張角為2°;凹腔位于隔離段入口下游447 mm處,長度為65 mm,深度為17 mm,傾角為22.5°,凹腔上游60.5 mm處安裝乙烯噴塊。凹腔底部安裝光纖火焰?zhèn)鞲衅髋c沖壓發(fā)動機(jī)點(diǎn)火器一體集成化火花塞輔助點(diǎn)火和燃燒狀態(tài)感知。

實(shí)驗(yàn)工況如表1所示,乙烯從6個直徑1.5 mm的孔中噴出,通過改變乙烯噴注壓力來實(shí)現(xiàn)當(dāng)量比的變化,進(jìn)而使發(fā)動機(jī)實(shí)現(xiàn)不同的工作狀態(tài),實(shí)驗(yàn)時(shí)間均為2 s。

表1 工況1、2的實(shí)驗(yàn)參數(shù)Table 1 Experimental parameters for condition 1 and 2

1.2 內(nèi)窺光纖傳感器系統(tǒng)

使用光纖傳感器系統(tǒng)捕捉凹腔燃燒室內(nèi)的CH*光信號與C*光信號,下文按照測點(diǎn)P1、P2、P3、P4進(jìn)行標(biāo)注和說明,安裝位置如圖2所示。

圖2 光纖傳感器的安裝位置Fig.2 Installation position of the optical fiber sensor

光纖傳感器系統(tǒng)由光纖準(zhǔn)直器、光纖和labview采集系統(tǒng)組成,如圖3所示。光纖準(zhǔn)直器經(jīng)過溫度跨度范圍20~1 600 K的系統(tǒng)熱容差分析設(shè)計(jì),考慮了光學(xué)元件曲率半徑和厚度變化以及介質(zhì)折射率受熱膨脹的影響,光學(xué)窗口直徑為2.5 mm,位于整個光路的最前端,與點(diǎn)火器集成在一起。為了使光纖準(zhǔn)直藍(lán)寶石晶體與不銹鋼殼體的氣密性、耐高溫結(jié)合,在一體化成型方面采用了高溫釬焊技術(shù)。在釬焊工藝方面,基于藍(lán)寶石晶體表面加工難度大、應(yīng)力敏感度高的物理特性,采取了焊接表面金屬化改性和真空釬焊兩步工藝措施。藍(lán)寶石晶體表面金屬化鍍鎳,通過鍍鎳層實(shí)現(xiàn)對藍(lán)寶石的致密包裹,起到應(yīng)力緩沖層的作用;真空釬焊工藝?yán)镁徛鋮s過程中焊接界面的應(yīng)力釋放和緩沖來進(jìn)行質(zhì)量控制。經(jīng)特殊設(shè)計(jì)的高溫光纖傳感器在五十余車超聲速燃燒實(shí)驗(yàn)熱考核中未見明顯損傷。

圖3 光纖傳感器系統(tǒng):Fig.3 Fiber optic sensor system

用430±10 nm的CH*帶通濾光片和514.5±10.0 nm的C*帶通濾光片對捕獲的光信號進(jìn)行濾波處理,然后通過光電倍增管對光信號進(jìn)行光電轉(zhuǎn)換。光電倍增管可以對入射的光電子進(jìn)行指數(shù)級別的放大,因此其比高速相機(jī)的CMOS傳感器具有更高的頻率響應(yīng)。本實(shí)驗(yàn)所選用光電倍增管具有300 kHz的動態(tài)頻率響應(yīng)。光電轉(zhuǎn)換后的光電壓信號通過模電轉(zhuǎn)換獲得實(shí)驗(yàn)需要的數(shù)字信號,最終實(shí)現(xiàn)1 MHz的采樣頻率,光路原理如圖4所示。

圖4 光路原理圖Fig.4 Optical path schematic

2 工作模態(tài)與穩(wěn)焰模式分析

2.1 工作模態(tài)分析

雙模態(tài)超燃沖壓發(fā)動機(jī)既能工作在超燃模態(tài)也能工作在亞燃模態(tài),這取決于燃燒室入口的馬赫數(shù)大小。由于燃燒室入口的馬赫數(shù)無法直接進(jìn)行測量,且亞燃模態(tài)與超燃模態(tài)的流場有較大差異,因此有學(xué)者通過激波串特征對2種模態(tài)做出判斷。此外,流場差異也會導(dǎo)致壓力場的不同,通過沿程布設(shè)壓力傳感器也能實(shí)現(xiàn)流場狀態(tài)的間接判斷,從而實(shí)現(xiàn)燃燒模態(tài)的識別?;谘爻虊毫Ψ植?,一維馬赫數(shù)分析方法可大致對流場的馬赫數(shù)分布進(jìn)行計(jì)算。

本文基于紋影法、沿程壓力和一維馬赫數(shù)分析方法,對燃燒模態(tài)進(jìn)行判斷。對應(yīng)的紋影圖像如圖5所示;2組固定當(dāng)量比狀態(tài)實(shí)驗(yàn)的沿程壓力如圖6所示,壓力曲線為實(shí)驗(yàn)時(shí)間內(nèi)的平均值;一維馬赫數(shù)分析如圖7所示。

圖5 工況1、2的紋影圖像Fig.5 Schlieren images of condition 1 and 2

圖6 工況1、2的沿程壓力分布Fig.6 Pressure distribution along the model of condition 1 and 2

圖7 工況1、2的沿程馬赫數(shù)分布Fig.7 Mach number distribution along the model of condition 1 and 2

當(dāng)乙烯的噴注壓力為1.5 MPa時(shí),從對應(yīng)的紋影圖像可以觀察到,乙烯橫向射流導(dǎo)致在超聲速來流中產(chǎn)生了一道略微彎曲的斜激波。在整個隔離段區(qū)間內(nèi),壓力曲線較為穩(wěn)定,但是在隔離段出口處壓力曲線出現(xiàn)一個較為尖銳的波動,幅值大約為25 kPa。由于壓力傳感器裝配在上壁面,而射流所導(dǎo)致的斜激波正好在隔離段出口的上壁面形成反射,因此這個尖銳的波動很有可能是這個激波所導(dǎo)致的。在燃燒室區(qū)域,燃燒釋熱所引起的壓力抬升對來流形成反壓,導(dǎo)致了第二條斜激波的產(chǎn)生。壓力曲線發(fā)生75 kPa的幅值抬升的位置,也與第二道斜激波的反射位置一致。

與之對應(yīng)的一維馬赫數(shù)分析也說明了在整個燃燒室內(nèi),氣流依舊保持為超聲速,這是一個非常典型的雙模態(tài)超燃沖壓發(fā)動機(jī)的超燃工作模態(tài)。

隨著乙烯噴注壓力增加到2.0 MPa,燃燒釋熱開始增強(qiáng),釋熱形成的反壓引起的邊界層分離導(dǎo)致了預(yù)燃激波串的產(chǎn)生,此時(shí)壓力曲線大幅抬升,達(dá)到200 kPa。當(dāng)量比的繼續(xù)升高使得壓力開始抬升的位置向上游移動,這是由于預(yù)燃激波串不斷向上游移動導(dǎo)致的,且壓力抬升的幅值比工況1更高。

對比工況2與工況1的的一維馬赫數(shù)曲線,可以看到,工況2氣流在燃燒室中為亞聲速狀態(tài)。結(jié)合紋影圖像、壓力曲線和一維馬赫數(shù)分析,可知工況2為亞燃工作模態(tài)。

2.2 穩(wěn)焰模式分析

在2.1小節(jié)的分析中已經(jīng)確定,工況1條件下,發(fā)動機(jī)的工作模態(tài)為超燃模態(tài);工況2條件下,發(fā)動機(jī)為亞燃模態(tài)。在2 s的實(shí)驗(yàn)時(shí)間內(nèi),計(jì)算幀率為2 000 幀/s 的CH*自發(fā)光圖片的均值與標(biāo)準(zhǔn)差,將圖片進(jìn)行偽彩化處理,得到的偽彩圖像如圖8所示。

圖8 工況1、2的CH*瞬時(shí)圖、均值圖與標(biāo)準(zhǔn)差圖Fig.8 The transient, mean and standard deviation of CH* images from condition 1 and 2

從CH*自發(fā)光的均值圖片中可以觀察到,工況1與2條件下火焰狀態(tài)最大的區(qū)別在于:工況1中,隔離段中基本沒有火焰,而工況2由于當(dāng)量比的增大,隔離段火焰的強(qiáng)度增加;前者為剪切層穩(wěn)焰模式,后者為射流尾跡穩(wěn)焰模式。值得一提的是,高速來流中的凹腔結(jié)構(gòu)會使流體在其內(nèi)部產(chǎn)生一個相對低速的回流區(qū),起到駐留火焰和維持燃燒的作用。而回流區(qū)與主流的邊界區(qū)域?yàn)榧羟袑?,剪切層?nèi)存在大量的渦結(jié)構(gòu),十分有利于燃料的混合。從流動過程上來講,高速來流壓制著噴出的乙烯,使之緊貼著壁面向下游運(yùn)動,并進(jìn)入剪切層中,其中一部分卷吸進(jìn)入凹腔中被點(diǎn)火器點(diǎn)燃,凹腔火焰同時(shí)引燃剪切層中的預(yù)混燃料。當(dāng)量比較低時(shí),燃燒釋熱所引起的壓力抬升不足以引起邊界層分離,導(dǎo)致隔離段中無法產(chǎn)生預(yù)燃激波串對來流進(jìn)行增壓減速,因此火焰無法自由傳播,由于上方的高速來流壓制,以及凹腔內(nèi)較低的局部當(dāng)量比導(dǎo)致火焰無法發(fā)展到凹腔中,最終穩(wěn)定為剪切層穩(wěn)焰模式;這是流動與燃燒共同作用的結(jié)果。工況2條件下,由于當(dāng)量比較高,燃燒釋熱更強(qiáng),巨大的反壓導(dǎo)致了隔離段中產(chǎn)生了一系列的預(yù)燃激波串,在激波串的作用下,燃燒室中的氣流速度降低為亞聲速。這時(shí)主流的壓制效應(yīng)瞬間消失,火焰迅速在亞聲速氣流中傳播,甚至上傳到隔離段中,最終形成在射流附近燃燒的尾跡穩(wěn)焰模式。這2種火焰狀態(tài)是超燃和亞燃工作模態(tài)所對應(yīng)的火焰穩(wěn)定模式,有利于科研人員從火焰狀態(tài)的角度對發(fā)動機(jī)工作模態(tài)進(jìn)行判斷。

3 內(nèi)窺光纖傳感器的燃燒診斷

3.1 釋熱分布測量

從圖2中可知,P1和P3測點(diǎn)為展向測點(diǎn),P2和P4測點(diǎn)為流向測點(diǎn),其中P4測點(diǎn)位于上游。圖9、10分別為工況1、2歸一化的釋熱分布,可以觀察到,工況2的光信號在0.6 s左右發(fā)生了一個突變,這是由于工況2為處在模態(tài)轉(zhuǎn)換邊界附近的亞燃工作模態(tài)。實(shí)驗(yàn)中展向測點(diǎn)的釋熱分布始終是不均勻的,P3測點(diǎn)的釋熱更強(qiáng),釋熱不均勻會導(dǎo)致發(fā)動機(jī)壁面受力與受熱不均勻,不利于發(fā)動機(jī)的長久運(yùn)行。展向燃燒釋熱的不均勻特征有可能是乙烯的噴注不對稱造成的。對于流向的測點(diǎn),可以看到下游測點(diǎn)的光強(qiáng)度始終高于上游測點(diǎn)的光強(qiáng)度,這是因?yàn)橄掠螠y點(diǎn)比上游測點(diǎn)處的燃料混合程度更高,燃燒更加劇烈,釋熱也更高。圖8中剪切層火焰的亮度也是從上游到下游逐漸增強(qiáng)。在工況2的射流尾跡穩(wěn)焰中,上游測點(diǎn)的光強(qiáng)度一直在衰減,而工況1剪切層穩(wěn)焰模式的光強(qiáng)度比較穩(wěn)定,這種現(xiàn)象可能與燃料的動態(tài)混合有關(guān),需要進(jìn)一步的研究和探索。

圖9 工況1條件下4個測點(diǎn)的CH*光信號Fig.9 CH* chemiluminescence of four measuring points under condition 1

3.2 振蕩特性分析

超燃沖壓發(fā)動機(jī)在正常工作時(shí),需避免因結(jié)構(gòu)設(shè)計(jì)問題或燃燒規(guī)律所帶來的劇烈壓力振蕩,這可能會對發(fā)動機(jī)的性能及結(jié)構(gòu)產(chǎn)生不利影響,包括引起熄火、降低燃燒效率、增大局部熱載荷、引起結(jié)構(gòu)的振動及破壞等。因此需要研究發(fā)動機(jī)工作時(shí)的壓力振蕩特性,以避免可能產(chǎn)生的危害。壓力振蕩一般由燃燒釋熱振蕩引起,因此,基于光纖傳感器信號的振蕩特性分析,是從本質(zhì)上去發(fā)現(xiàn)問題的創(chuàng)新性方法。本小節(jié)在定當(dāng)量比實(shí)驗(yàn)條件下,對燃燒振蕩特性進(jìn)行了測量與分析。

圖10 工況2條件下4個測點(diǎn)的CH*光信號Fig.10 CH* chemiluminescence of four measuring points under condition 2

3.2.1 壓力信號的FFT分析

在乙烯噴注壓力分別為1.5、2.0 MPa的實(shí)驗(yàn)中,對隔離段壓力測點(diǎn)CH1和凹腔內(nèi)的壓力測點(diǎn)CH3的高頻壓力信號進(jìn)行了快速傅里葉變換(FFT),其振蕩特性曲線如圖11、12所示。

圖11 工況1高頻壓力信號的頻域曲線Fig.11 FFT of pressure under condition 1

可以觀察到,2個工況的隔離段測點(diǎn)CH1和凹腔測點(diǎn)CH3均存在2 300 Hz左右的振蕩主頻。也就是說,在不同的實(shí)驗(yàn)工況下,當(dāng)量比的變化與穩(wěn)焰模式的變化并不對2 300 Hz左右的壓力振蕩主頻產(chǎn)生影響,說明該主頻的產(chǎn)生與燃燒狀態(tài)無關(guān)??紤]到無燃燒因素的的高頻振蕩機(jī)理,凹腔自激振蕩或許可以解釋這種現(xiàn)象,Rossiter綜合了聲波傳播和渦耗散的特點(diǎn),提出了計(jì)算凹腔自激振蕩的經(jīng)驗(yàn)公式:

式中:L為凹腔的長度, M、U分別為自由來流的馬赫數(shù)與速度,m為對應(yīng)耦合的模態(tài)數(shù), f為在此模態(tài)下的頻率,α和k為經(jīng)驗(yàn)常數(shù)。Heller和Bliss考慮流體的可壓縮效應(yīng)對上式進(jìn)行了修正,修正后的公式為:

圖12 工況2高頻壓力信號的頻域曲線Fig.12 FFT of pressure under condition 2

γ為氣體的比熱比,基于這2個計(jì)算凹腔自激振蕩的公式,在 L =0.1m、 m =1、 a =0.25、 k =0.57、M=2.8條件下,計(jì)算得到基于公式(1)的2 414 Hz和基于公式(2)的3 077 Hz。與前面的FFT分析結(jié)果進(jìn)行對比,發(fā)現(xiàn)Rossiter的公式更加符合當(dāng)前實(shí)際。在Heller和Bliss的實(shí)驗(yàn)研究中,發(fā)現(xiàn)凹腔的長深比 L/D對 f的預(yù)測結(jié)果影響較大,公式(2)的計(jì)算結(jié)果偏差較大可能是長深比導(dǎo)致的。值得注意的是,隔離段內(nèi)的壓力測點(diǎn)為什么會測量到凹腔自激振蕩帶來的壓力變化,或者說隔離段內(nèi)的高頻振蕩是否源于其他的振蕩機(jī)制,目前并沒有明確的解釋,有待學(xué)者進(jìn)一步的研究。

除了2 300 Hz的高頻振蕩頻率,在工況2條件下,還發(fā)現(xiàn)了500 Hz的低頻振蕩。一般燃燒室中這種百赫茲量級的振蕩主頻有可能是由熱聲耦合振蕩導(dǎo)致,具體表現(xiàn)為聲波在燃燒區(qū)和射流區(qū)之間或燃燒區(qū)與激波區(qū)之間來回振蕩所致。此外,也有可能是火焰低頻周期振蕩導(dǎo)致的,但在工況1條件下并沒有發(fā)現(xiàn)這種低頻振蕩,這是因?yàn)楣r2為亞燃模態(tài),而工況1為超燃模態(tài),可能在超燃模態(tài)下壓力振蕩無法傳導(dǎo)到壓力測點(diǎn)處,從而造成了這種差異。

3.2.2 CH自發(fā)光信號的FFT分析

基于從釋熱角度研究振蕩特性的考慮,同樣對實(shí)驗(yàn)中的火焰CH*自發(fā)光信號進(jìn)行了FFT變換,如圖13、14所示。從圖中可見,在乙烯噴注壓力為1.5 MPa的工況1條件下,光信號的FFT分析也發(fā)現(xiàn)了500 Hz左右的振蕩主頻,與高頻壓力的FFT分析結(jié)果一致。在3.2.1小節(jié)中,對低頻振蕩的產(chǎn)生做出了2種解釋,即流向熱聲振蕩和火焰周期振蕩,在超燃模態(tài)下,由于燃燒室內(nèi)為超聲速流場,因此流向的熱聲振蕩是不存在的。但從目前光信號的FFT分析結(jié)果來看,在超燃模態(tài)下依舊出現(xiàn)了500 Hz的振蕩主頻。

圖13 工況1 CH*光信號的頻域曲線Fig.13 FFT of CH* chemiluminescence under condition 1

圖14 工況2 CH*光信號的頻域曲線Fig.14 FFTof CH* chemiluminescence under condition 2

對于500 Hz的振蕩主頻,同樣也有2種解釋。一是展向的熱聲振蕩,即燃燒區(qū)的釋熱振蕩引起的壓力擾動在凹腔底部和燃燒區(qū)之間來回傳播所導(dǎo)致的。二是火焰周期振蕩,隨著當(dāng)量比的增加,這種500 Hz的低頻振蕩在工況2中慢慢發(fā)展為800 Hz左右的振蕩主頻,從展向的熱聲振蕩角度來說,這是因?yàn)楫?dāng)量比的增加伴隨著燃燒區(qū)的擴(kuò)大,同時(shí)凹腔底部與燃燒區(qū)的距離減小,于是展向的熱聲振蕩主頻逐漸增大。但基于高頻壓力的FFT分析,并沒有明顯觀察到低頻逐漸增大的現(xiàn)象。由于從釋熱振蕩到壓力傳感器的探測,需要以流體作為介質(zhì)傳導(dǎo)壓力反饋,從而涉及到跟隨性的問題,這可能是光信號和壓力信號的振蕩主頻不一致的原因。此外值得注意的是,在光信號的FFT分析中,并沒有觀察到2 300 Hz左右的高頻振蕩,這可能是凹腔自激振蕩并未與燃燒振蕩耦合導(dǎo)致的。

3.3 C2*/CH*表征局部當(dāng)量比

3.3.1 局部當(dāng)量比與釋熱的時(shí)空關(guān)聯(lián)性

L型內(nèi)燃機(jī)在工作過程中,普遍存在燃料混合不均勻的問題,而局部過度富燃或貧燃都會使燃燒不充分,導(dǎo)致有害氣體排放和燃料浪費(fèi)。燃油與空氣混合特征取決于氣缸內(nèi)的渦流和湍流,并與進(jìn)氣歧管的結(jié)構(gòu)、供油系統(tǒng)的設(shè)計(jì)和燃燒室的形狀直接相關(guān),了解燃料混合特性有助于L型內(nèi)燃機(jī)的設(shè)計(jì)。Chou等在L型內(nèi)燃機(jī)氣缸內(nèi)燃料分布的研究中提出基于不同燃料的條件下C*/CH*與當(dāng)量比的定量關(guān)系,不過這種定量關(guān)系并非普適的,而取決于發(fā)動機(jī)的工況。C*/CH*與當(dāng)量比的正相關(guān)關(guān)系在后續(xù)的層流預(yù)混火焰的研究中也得到了證實(shí)。類比于L型內(nèi)燃機(jī)的科研需要,超燃沖壓發(fā)動機(jī)也需要考慮燃料混合的均勻性問題,此外,考慮到壓力信號作為燃燒主動控制的反饋信號存在一定延遲,因此需要深入了解火焰自發(fā)光與燃燒狀態(tài)的相關(guān)性。

通過光纖傳感器,在凹腔內(nèi)同時(shí)測量了C*和CH*光信號,以C*/CH*表征當(dāng)量比,工況1、2條件下以自發(fā)光表征的釋熱和當(dāng)量比的關(guān)系分別如圖15、16所示。可以觀察到,2次實(shí)驗(yàn)的任何一個時(shí)間段,C*/CH*表征的當(dāng)量比始終是測點(diǎn)P1>P2>P4>P3,這可能是乙烯的噴注不對稱或者是一些未知的發(fā)動機(jī)模型結(jié)構(gòu)設(shè)計(jì)導(dǎo)致的。

圖15 工況1條件下CH*、C2*/CH*與時(shí)間的關(guān)系Fig.15 CH* and C2*/CH* versus time under condition 1

3.3.2 局部窗口的CH滯后時(shí)間計(jì)算

在任意毫秒級局部尺度內(nèi),對2組實(shí)驗(yàn)P3測點(diǎn)的CH*信號與C*/CH*進(jìn)行對比觀察,結(jié)果如圖17、18所示,CH*和C*/CH*都表現(xiàn)出了類周期性的振蕩,對比二者的振蕩峰值和波形的相似度,可以發(fā)現(xiàn)C*/CH*基本上是先于CH*發(fā)生變化的。為了計(jì)算CH*的延遲時(shí)間,對C*/CH*(t,t+T)的時(shí)間序列和CH*(t+τ,t+T+τ)的時(shí)間序列做相關(guān)性分析,其中T為窗口長度,t為時(shí)刻零點(diǎn),τ為延遲時(shí)間。

圖16 工況2條件下CH*、C2*/CH*與時(shí)間的關(guān)系Fig.16 CH* and C2*/CH* versus time under condition 2

固定C*/CH*的時(shí)間序列,移動CH*時(shí)間序列,如工況1的窗口長度T=6 ms,時(shí)刻零點(diǎn)t=711 ms,工況2的窗口長度T=4 ms,時(shí)刻零點(diǎn)t=269 ms,依次對應(yīng)。τ為延遲時(shí)間(0~10 ms),每一個延遲時(shí)間τ對應(yīng)一個相關(guān)系數(shù)。在微秒量級分辨率下做10 000次相關(guān)性分析,其結(jié)果如圖17、18所示??梢钥吹?組工況的相關(guān)系數(shù)大小也呈現(xiàn)出類周期性的變化。理論上,由于CH*相對于C*/CH*存在延遲,因此,相關(guān)系數(shù)曲線在初始零點(diǎn)隨著延遲時(shí)間的增大,其相關(guān)系數(shù)也會隨之增大。

圖17 工況1 CH*和C2*/CH*的局部相關(guān)性Fig.17 Local correlation between CH* and C2*/CH* under condition 1

3.3.3 全局CH滯后時(shí)間的統(tǒng)計(jì)特征

在不同測點(diǎn)選擇不同的窗口,CH*的延遲時(shí)間并不一致。為了在全局上估計(jì)延遲時(shí)間,在2 s的時(shí)間尺度內(nèi)對延遲時(shí)間進(jìn)行了統(tǒng)計(jì)分析,與單個窗口的延遲時(shí)間計(jì)算方法相同,劃分2 s的實(shí)驗(yàn)時(shí)間為500個窗口,每個窗口長度T為4 ms??紤]到CH*延遲時(shí)間為微秒級,因此選擇τ為0~1 ms。每個窗口計(jì)算1 000個相關(guān)系數(shù),取相關(guān)系數(shù)最大值對應(yīng)的延遲時(shí)間為CH*延遲時(shí)間。對500個窗口的CH*延遲時(shí)間做分布統(tǒng)計(jì),每一組實(shí)驗(yàn)不同測點(diǎn)延遲時(shí)間的統(tǒng)計(jì)特性結(jié)果如圖19、20所示。在這2個實(shí)驗(yàn)工況中,可以觀察到,絕大部分的窗口中, C*/CH*與CH*的同步性較強(qiáng),CH*延遲時(shí)間集中在0~200 us。

圖18 工況2 CH*和C2*/CH*的局部相關(guān)性Fig.18 Local correlation between CH* and C2*/CH* under condition 2

圖19 工況1 CH*延遲時(shí)間的統(tǒng)計(jì)特征Fig.19 Statistical characteristics of CH* delay time under condition 1

3.4 局部火焰質(zhì)心位置的測量

通過2個通道的光纖測點(diǎn)布設(shè)方式求解局部火焰質(zhì)心,可得到一維火焰位置的離散時(shí)間序列?;诠饫w測點(diǎn)的布設(shè)方式,沿流向的火焰位置計(jì)算公式如下所示:

式中:P1、P2、P3和P4測點(diǎn)的相對位置關(guān)系如圖2所示, x=x=11mm,I為歸一化的CH*光強(qiáng)度。對于基于CH*平面自發(fā)光成像的一維火焰質(zhì)心位置,可由相片的像素矩陣計(jì)算得到,計(jì)算表達(dá)式為:

圖20 工況2 CH*延遲時(shí)間的統(tǒng)計(jì)特征Fig.20 Statistical characteristics of CH* delay time under condition 2

式中:r和q分別為空間每個像素點(diǎn)的矩陣坐標(biāo)與像素值大小,Q為所有像素點(diǎn)值的總和。定當(dāng)量比條件下火焰質(zhì)心位置的概率分布如圖21所示?;谄矫娉上裼?jì)算的全局火焰質(zhì)心位置的概率分布顯示:在2個實(shí)驗(yàn)條件下,火焰質(zhì)心存在2個主要的穩(wěn)定位置,分別在點(diǎn)火器上游15 mm的A處和點(diǎn)火器下游25 mm的B處。當(dāng)量比的增加會導(dǎo)致B處概率峰值的降低和A處概率峰值的升高,這是因?yàn)獒専岬奶鰪?qiáng)了預(yù)燃激波串,使火焰更易于穩(wěn)定在上游區(qū)域。對火焰質(zhì)心的FFT分析未發(fā)現(xiàn)主頻的存在,因此火焰質(zhì)心在A和B之間的運(yùn)動沒有確定的周期,為不規(guī)則脈動?;诠饫w傳感器計(jì)算的局部火焰質(zhì)心為單峰形態(tài),相比于射流尾跡火焰,剪切層穩(wěn)焰模式的局部火焰質(zhì)心位置位于上游,因此光纖傳感器也可以通過質(zhì)心位置的概率分布辨別這2種穩(wěn)焰模式。

圖21 火焰質(zhì)心位置的概率分布曲線Fig.21 Probability distribution curve of flame centroid position

4 總 結(jié)

本文介紹了一種基于被動火焰自發(fā)光譜的內(nèi)窺式光纖火焰?zhèn)鞲衅?,通過發(fā)展特種高溫?zé)崦芊夤に嚱鉀Q了傳感器發(fā)動機(jī)惡劣環(huán)境適應(yīng)性問題,并驗(yàn)證了光纖火焰?zhèn)鞲衅鲾?shù)據(jù)的燃燒過程感知價(jià)值。

1)內(nèi)窺式光纖傳感器可感知燃燒室釋熱率的時(shí)空演變特性。在當(dāng)前的實(shí)驗(yàn)條件和測點(diǎn)位置上,超燃狀態(tài)相比于亞燃狀態(tài),其燃燒釋熱率明顯更弱。

2)內(nèi)窺式光纖傳感器可感知頻域燃燒振蕩特性。CH*自發(fā)光信號和壓力信號的FFT分析表明在燃燒過程中可能存在凹腔自激振蕩和展向的熱聲振蕩現(xiàn)象。燃燒釋熱振蕩引起壓力的振蕩,壓力的FFT分析證實(shí)燃燒時(shí)存在2 300 Hz的高頻振蕩,同時(shí)在亞燃條件下存在500 Hz的低頻振蕩。CH*自發(fā)光的FFT結(jié)果表明,在所有條件下均存在500~1 000 Hz的低頻振蕩,但不存在千赫茲量級的振蕩主頻,2 300 Hz的高頻振蕩可能是凹腔自激振蕩導(dǎo)致的,500 Hz的低頻振蕩則可能源于展向的熱聲振蕩。

3)內(nèi)窺式光纖傳感器C*/CH*光信號可感知局部當(dāng)量比的時(shí)空演變特性,結(jié)合CH*光信號可應(yīng)用于混合場與燃燒場關(guān)聯(lián)性的研究。C*/CH*表征的局部當(dāng)量比與CH*自發(fā)光信號在局部存在明顯的時(shí)間相關(guān)性,CH*普遍滯后于C*/CH*10~10微秒。同時(shí)由于內(nèi)窺式光纖傳感器具有流向及展向二維空間分辨能力,實(shí)驗(yàn)表明,局部當(dāng)量比與釋熱率之間不具有絕對的空間相關(guān)性,這說明混合場與燃燒場關(guān)聯(lián)具有很強(qiáng)的時(shí)間及空間尺度非線性特征,需系統(tǒng)進(jìn)行非線性定量化以獲得更豐富的燃燒場認(rèn)知。

4)流向全局和局部的火焰質(zhì)心位置統(tǒng)計(jì)特征表征了剪切層穩(wěn)焰模式和射流尾跡穩(wěn)焰模式。全局的火焰質(zhì)心位置概率分布曲線以雙峰的形式描述了不同當(dāng)量比實(shí)驗(yàn)中火焰質(zhì)心在不同位置的穩(wěn)定性變化。

龙口市| 盐山县| 宜阳县| 侯马市| 奉贤区| 左云县| 会东县| 平阳县| 班戈县| 晋江市| 唐海县| 磐安县| 开平市| 盈江县| 孝义市| 锡林浩特市| 武定县| 报价| 苍山县| 青神县| 正宁县| 南木林县| 大竹县| 恩平市| 吉水县| 文安县| 康乐县| 苍溪县| 绍兴县| 寻甸| 乌鲁木齐市| 孟村| 长沙县| 三原县| 依安县| 大冶市| 那曲县| 神木县| 盈江县| 四川省| 城市|