譚錕,王云,張飛,張洋洋,李俊杰
(南昌航空大學(xué) 飛行器工程學(xué)院,南昌 330063)
噴油環(huán)是微型渦噴發(fā)動機(jī)的重要組成部分,噴油環(huán)包括噴油環(huán)體及噴油針兩部分,其中噴油針噴油量的均勻性決定了燃燒室溫度場的均勻性。國外較多研究大型渦噴發(fā)動機(jī)噴油環(huán)及單個噴嘴等核心部件,或研究燃油速度變化[1-3],目前國內(nèi)主要是針對大型渦噴發(fā)動機(jī)噴油環(huán)及噴嘴動態(tài)特性及其內(nèi)氣液兩相流動進(jìn)行數(shù)值模擬研究[4],而對微渦噴發(fā)動機(jī)噴油環(huán)流量分布及均勻特性相關(guān)的研究較少。
本文采用數(shù)值模擬的方法對微渦噴發(fā)動機(jī)噴油環(huán)噴油針出口流量特性進(jìn)行分析改進(jìn)。首先建立噴油環(huán)內(nèi)部流場數(shù)學(xué)模型并劃分網(wǎng)格,根據(jù)所研究噴油環(huán)的結(jié)構(gòu)特點(diǎn)和實(shí)際測試數(shù)據(jù)確定合理的邊界條件,應(yīng)用CFX軟件通過數(shù)值模擬的方法對噴油環(huán)內(nèi)部流場分布、噴油針出口流量及速度特性進(jìn)行初步數(shù)值模擬計算,并與實(shí)驗(yàn)結(jié)果進(jìn)行比較,分析原因并給出相應(yīng)解決方案,以期為研究噴油環(huán)內(nèi)部流場并做出分析奠定基礎(chǔ)。
某微渦噴發(fā)動機(jī)噴油環(huán)如圖1所示,噴油環(huán)由進(jìn)油管和分布在環(huán)體上的16根噴油針組合而成。通過測繪可得噴油環(huán)具體數(shù)據(jù)如下:進(jìn)油管外徑3.2 mm,壁厚0.1 mm;環(huán)體管外徑5 mm,壁厚0.5 mm;噴油針管外徑1 mm,壁厚0.1 mm。
圖1 微渦噴發(fā)動機(jī)噴油環(huán)Fig.1 Micro-turbojet fuel injection ring
燃油從噴油針出口高速噴出,之后在高溫蒸發(fā)管內(nèi)霧化,使燃油能夠在燃燒室內(nèi)充分燃燒,霧化的效果也會直接影響燃燒效率。若各噴油針出口的流量不同則進(jìn)入蒸發(fā)管內(nèi)的燃油量不均勻,最終導(dǎo)致燃燒室內(nèi)燃油燃燒不均勻,而未能在燃燒室內(nèi)充分燃燒的霧化燃油會流至下一級渦輪,直接造成渦輪損傷,影響渦輪的工作壽命,因此應(yīng)保證從噴油針噴入蒸發(fā)管內(nèi)的燃油流量基本相等,從而使燃燒室內(nèi)的燃油充分燃燒。而噴油針出口流量(kg/s)與噴油針在環(huán)體上分布位置及孔徑有關(guān)。
對噴油環(huán)做流量測試實(shí)驗(yàn)所需儀器:壓力計、油泵、調(diào)壓計、電源、計重器、油箱、三通管、計時器、記號筆、導(dǎo)線若干、內(nèi)徑4 mm和1 mm橡膠管若干;500 ml容器若干。組裝好所有實(shí)驗(yàn)設(shè)備及儀器(如圖2所示),通電等待測試實(shí)驗(yàn)。實(shí)驗(yàn)原理:接通電源啟動油泵時開始計時,將油泵快速調(diào)至設(shè)定值,油泵將油箱內(nèi)的航空煤油高壓壓出;一路流到壓力計,監(jiān)測壓力情況;另一路流到噴油環(huán)進(jìn)油管再到環(huán)體,從環(huán)體上均勻分布的噴油針高速流出,經(jīng)延長橡膠管流入容器內(nèi),試驗(yàn)結(jié)束,記錄時間,并用計重器稱重各容器內(nèi)的燃油質(zhì)量,對比各容器內(nèi)的質(zhì)量流量。
圖2 噴油環(huán)流量測試實(shí)驗(yàn)Fig.2 Fuel injection flow test
通過多次噴油環(huán)流量測試實(shí)驗(yàn)分析對比,最終取一組實(shí)驗(yàn)結(jié)果。設(shè)定入口總壓為3.3 bar,作為實(shí)驗(yàn)壓力;為了使實(shí)驗(yàn)結(jié)果更為準(zhǔn)確,時間設(shè)定為3 min;待試驗(yàn)結(jié)束稱重各容器內(nèi)航空煤油質(zhì)量,可得各噴油針單位時間內(nèi)質(zhì)量流量(如表1所示),總出口流量為0.024 22 kg/s,換算成體積還可計算出各噴油針出口速度。
表1 各噴油針出口質(zhì)量流量Table 1 The injector outlet mass flow
2.2.1 物理模型及數(shù)學(xué)模型的建立
建立物理模型(如圖3所示);燃油在管道內(nèi)高速流動過程中情況比較復(fù)雜,模型按照三維流動來分析;考慮到燃油粘度變化、可壓縮性以及燃油和噴油環(huán)壁面間的摩擦,在數(shù)值模擬計算過程中需要運(yùn)用CFX軟件求解質(zhì)量和動量守恒方程。對于流動是湍流時,需要添加k-Epsilon運(yùn)輸方程。
圖3 噴油環(huán)結(jié)構(gòu)模型Fig.3 Fuel injection ring structure model
2.2.2 CFX計算網(wǎng)格的劃分及邊界條件設(shè)定
噴油環(huán)內(nèi)部流道網(wǎng)格局部模型如圖4所示,通過對比分析可以得出:噴油針流道的幾何尺寸相對于整個噴油環(huán)環(huán)體和進(jìn)油管流道具有尺寸小且復(fù)雜存在一定彎度??紤]到燃油在噴油針流道內(nèi)可能會存在劇烈的湍流和旋流運(yùn)動情況,對噴油環(huán)環(huán)體和噴油針流道需要不同精度網(wǎng)格劃分可以提高精算結(jié)果的精度。噴油環(huán)流道采用四面體非結(jié)構(gòu)網(wǎng)格劃分,網(wǎng)格數(shù)量大約150萬,其四面體網(wǎng)格占95%以上;在環(huán)體至噴油針段流道截面積驟然減小造成壓力急劇變化,對此情況需要對網(wǎng)格做進(jìn)一步的細(xì)分,使網(wǎng)格的尺寸足夠小,從而保證計算結(jié)果的精確。
圖4 噴油環(huán)局部網(wǎng)格模型Fig.4 Injector ring partial grid model
在CFX預(yù)處理中,對噴油環(huán)模型設(shè)定邊界條件,設(shè)定進(jìn)油管的入口截面為Inlet,各噴油針出口截面為Outlet并分別定義,其余面均設(shè)為Wall。材料定義欄中,給定燃油摩爾質(zhì)量為220 kg/kmol,密度為770 kg/m3,工作介質(zhì)(航空煤油)動力粘度μ=2.21。
2.2.3 數(shù)值模擬分析結(jié)果
等徑排布噴油環(huán)速度流線圖如圖5所示,可以看出:距離進(jìn)油管越遠(yuǎn)的噴油針出口速度流線密集度越低,反之越高。在數(shù)值模擬過程中發(fā)現(xiàn):燃油在流動過程中存在劇烈的湍流和漩流。
圖5 等徑排布噴油環(huán)速度流線圖Fig.5 Isometric distribution of injection ring velocity flow chart
等徑排布噴油環(huán)壓力圖如圖6所示,從進(jìn)油口至環(huán)體最后到噴油針存在明顯的壓力損失,這是由于流體轉(zhuǎn)彎太急,壓力損失較大,容易在拐角處形成低壓區(qū);通過數(shù)值模擬計算發(fā)現(xiàn)各噴油針出口流量不一致。
圖6 等徑排布噴油環(huán)壓力圖Fig.6 Equal arrangement of spray ring spray pressure map
通過CFX數(shù)值模擬,當(dāng)各噴油針出口流量趨于穩(wěn)定不再發(fā)生變化時,計算單位時間各噴油針出口流量,如表2所示。
表2 各噴油針出口質(zhì)量流量Table 2 The injector outlet mass flow
從表2可以看出:距離進(jìn)油管最近的1、2、15和16號噴油針出口流量基本相等;離進(jìn)油管最遠(yuǎn)的7、8、9、10號噴油針出口流量減少非常明顯;總出口流量為0.024 178 4 kg/s。對比分析可知,噴油針距出口流量與距離進(jìn)油管的遠(yuǎn)近及孔徑有關(guān)。
通過數(shù)值模擬的方法分析噴油環(huán)與實(shí)驗(yàn)出口流量結(jié)果對比(如表1和表2所示),可以看出:數(shù)值模擬得到的結(jié)果與實(shí)驗(yàn)結(jié)果非常接近,由于誤差很小,說明該方案切實(shí)可行。
為了進(jìn)一步提高改善燃燒室燃燒性能,有必要對噴油環(huán)噴油結(jié)構(gòu)裝置進(jìn)行優(yōu)化改進(jìn),具體方案有兩種:一種是采用噴油針不等距分布,以改善沿噴油環(huán)環(huán)向空間的噴油均勻性,但該方法會導(dǎo)致后面燃油蒸發(fā)管的不均勻布置帶來的結(jié)構(gòu)問題;一種是采用不等徑噴油針方案,通過調(diào)整噴油針噴口孔徑來保持各個位置的噴油量盡量均勻。第二種方案對原有結(jié)構(gòu)沒有影響,簡單可行。
噴油針管徑改進(jìn)前后方案的尺寸如圖7所示。
圖7 噴油針管徑改進(jìn)后前尺寸Fig.7 Injection needle diameter improved after the former size
通過數(shù)值模擬分析距離進(jìn)油管最近的噴油針(1號)和最遠(yuǎn)的噴油針(7號)可知:1號噴油針與環(huán)體的連接處附近的旋流、湍流較7號噴油針相更為劇烈一些,如圖8和圖9所示。進(jìn)油管至噴油環(huán)環(huán)體部分及環(huán)體至噴油針壓降都十分明顯(如圖10所示),這是由于管道截面積急劇減小而造成的。
(a) 1號
(b) 7號圖8 不等徑噴油針距1號和7號噴油針與環(huán)體連接處附近旋流情況圖Fig.8 Unequal diameter injection needle from the 1st and 7st injection needle and the ring near the junction of the cyclone situation
(a) 1號
(b) 7號圖9 不等徑噴油針距1號和7號噴油針與環(huán)體連接處附近湍流情況圖Fig.9 Unequal diameter injection needle from the 1st and 7st injection needle and the ring near the junction of the turbulent situation
圖10 不等徑噴油針部分壓力圖Fig.10 Particle pressure map of unequaldiameter injection needle
在數(shù)值模擬中,當(dāng)各噴油針流量趨于穩(wěn)定后對各噴油針出口流量進(jìn)行計算(如表3所示),可以看出:在改進(jìn)后不等徑噴油針方案中距離最遠(yuǎn)的7~10號噴油針出口流量得到較好的改善;1~3、14~16號噴油針未改變管徑,所以流量變化不大;4~6、11~13號噴油針其質(zhì)量流量也大致趨近于17.0 g;總出口流量為0.027 350 25 kg/s。
表3 各噴油針出口質(zhì)量流量Table 3 The injector outlet mass flow
出口質(zhì)量流量效果對比如圖11所示,綜合考慮通過采用不等徑噴油針方案,可較好地控制各噴油針的出口質(zhì)量流量,且使各噴油針出口質(zhì)量流量差值最小,以保證各噴油針出口質(zhì)量流量均勻。
圖11 出口質(zhì)量流量效果對比Fig.11 Comparison of export mass flow effect
通過數(shù)值模擬一定時間后,當(dāng)各噴油針出口流量趨于穩(wěn)定且不再發(fā)生劇烈變化,根據(jù)各噴油針出口質(zhì)量流量來計算各噴油針出口平均速度。改進(jìn)前后噴油針出口速度對比如圖12所示。
圖12 改進(jìn)前后噴油針出口速度對比Fig.12 Comparison of improvement before and after the injection needle outlet speed
從圖12可以看出:改進(jìn)前各噴油針出口速度為13.0~18.0 m/s,計算其速度方差并記為D(X1)=2.294 740;改進(jìn)后各噴油針出口速度為15.5~18.0 m/s,計算其速度方差并記為D(X2)=0.748 507;改進(jìn)前后各噴油針出口速度方差D(X2)遠(yuǎn)小于D(X1),表明不等徑噴油針結(jié)構(gòu)方案取得的效果更好。
(1) 各等徑排布噴油針單位時間內(nèi)出口質(zhì)量流量不均勻,距離進(jìn)油管越遠(yuǎn)的噴油針出口流量越小,反之越大;從進(jìn)油口至環(huán)體和噴油針的壓力損失十分明顯,其原因是由于管道截面積急劇減小。
(2) 數(shù)值模擬出的結(jié)果與實(shí)驗(yàn)結(jié)果非常接近;基于該模擬方法,提出不等徑噴油針結(jié)構(gòu)方案,以改進(jìn)之前存在的噴油環(huán)噴油量不均勻的問題;通過數(shù)值模擬分析對比發(fā)現(xiàn),不等徑噴油針結(jié)構(gòu)方案可以較好地改善各噴油針出油量不均勻問題,通過調(diào)整噴油針口徑來使得各噴油針質(zhì)量流量趨于均勻。
(3) 數(shù)值模擬結(jié)果的兩組出口速度方差D(X1)=2.294 740、D(X2)=0.748 507,可以看出不等徑噴油針結(jié)構(gòu)方案的方差最小,表明通過調(diào)整噴油針口徑,還可以較好地保證各噴油針燃油噴出的速度,從而保證進(jìn)入蒸發(fā)管內(nèi)的燃油速度都較為均勻。