發(fā)動機(jī)部件防火試驗(yàn)的數(shù)值模擬研究

尹莉萍，方鵬亞，單晨晨

(鄭州航空工業(yè)管理學(xué)院 航空工程學(xué)院，河南 鄭州 450046)

0 引言

《運(yùn)輸類飛機(jī)適航規(guī)定》(CCAR-25)[1]和《航空發(fā)動機(jī)適航規(guī)定》(CCAR-33)[2]都對航空發(fā)動機(jī)部件防火提出了明確要求，而且航空發(fā)動機(jī)部件一旦失火，必然會造成嚴(yán)重的安全事故，因此為了通過發(fā)動機(jī)適航驗(yàn)證，保證發(fā)動機(jī)部件及整個(gè)發(fā)動機(jī)的安全性，必須開展發(fā)動機(jī)部件防火符合性驗(yàn)證，而防火試驗(yàn)通常是最常用的、最有效的發(fā)動機(jī)部件防火符合性驗(yàn)證方法。為此，美國聯(lián)邦航空局(Federal Aviation Administration,FAA)、歐洲航空安全局(European Aviation Safety Agency, EASA)和中國民用航空局分別發(fā)布了《動力裝置防火的試驗(yàn)方法、標(biāo)準(zhǔn)和準(zhǔn)則》(AC 20-135)[3]、《民用飛機(jī)機(jī)載設(shè)備環(huán)境條件和試驗(yàn)方法-指定火區(qū)的防火試驗(yàn)》(ISO2685)[4]和《民用飛機(jī)機(jī)載設(shè)備環(huán)境條件和試驗(yàn)方法-指定火區(qū)的防火試驗(yàn)》(HB 6167.14-92)[5]等標(biāo)準(zhǔn)和文件，以指導(dǎo)防火試驗(yàn)的開展。根據(jù)AC 20-135可知，目前在生產(chǎn)的能用于防火試驗(yàn)的燃燒器中只有NexGen燃燒器，因此為了進(jìn)一步指導(dǎo)防火試驗(yàn)，需要開展NexGen燃燒器的火焰特性(即NexGen燃燒器在指定位置產(chǎn)生的標(biāo)準(zhǔn)火焰溫度和熱流密度分布)和典型發(fā)動機(jī)部件對火焰響應(yīng)的研究。由于NexGen燃燒器價(jià)格昂貴(燃燒器及其配套設(shè)備價(jià)格高達(dá)數(shù)百萬)、對實(shí)驗(yàn)條件要求較高且存在一定的安全隱患，因此越來越多的人采用數(shù)值計(jì)算的方法研究NexGen燃燒器的火焰特性，得到NexGen燃燒器產(chǎn)生的標(biāo)準(zhǔn)火焰分布情況，并以此為基礎(chǔ)研究發(fā)動機(jī)部件對標(biāo)準(zhǔn)火焰的響應(yīng)情況。

本文旨在通過數(shù)值計(jì)算的方法研究發(fā)動機(jī)部件防火試驗(yàn)，首先通過與實(shí)驗(yàn)數(shù)據(jù)比較建立NexGen燃燒器自由火焰的數(shù)值計(jì)算方案，由此獲得NexGen燃燒器出口指定位置處的標(biāo)準(zhǔn)火焰溫度和熱流密度分布，然后以此為初始邊界條件研究發(fā)動機(jī)部件在標(biāo)準(zhǔn)火焰中的響應(yīng)情況，從而建立發(fā)動機(jī)部件防火試驗(yàn)的數(shù)值計(jì)算方法，為發(fā)動機(jī)部件防火試驗(yàn)及申請人開展防火符合性驗(yàn)證提供一定的參考和指導(dǎo)。

1 NexGen燃燒器自由火焰的數(shù)值計(jì)算方法

首先通過查閱相關(guān)資料建立NexGen燃燒器的幾何模型[6]和計(jì)算域模型，以觀察出口火焰分布，然后對建立的幾何模型實(shí)施網(wǎng)格劃分，并通過網(wǎng)格無關(guān)性驗(yàn)證確認(rèn)網(wǎng)格模型的有效性，最后通過將計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果對比，驗(yàn)證計(jì)算方法的有效性[7]。NexGen燃燒器的計(jì)算網(wǎng)格模型如圖1所示。

圖1 NexGen燃燒器的計(jì)算網(wǎng)格模型

在ANSYS Fluent中設(shè)置相應(yīng)的邊界條件進(jìn)行數(shù)值計(jì)算，為獲得適用的仿真計(jì)算方案，根據(jù)NexGen燃燒器內(nèi)旋轉(zhuǎn)擴(kuò)張流動、出口射流特點(diǎn)以及燃油和空氣混合燃燒情況，結(jié)合同類型燃燒流場的研究結(jié)果，選擇RNGk-ε模型、Realizablek-ε模型、SSTk-ω模型三種湍流模型[8]和渦耗散模型(EDM)、混合分?jǐn)?shù)PDF模型兩種燃燒模型[9，10]進(jìn)行組合研究。

計(jì)算時(shí)具體邊界條件設(shè)置如表1所示。表1中，v為速度，T為溫度，I為湍流強(qiáng)度，d為水力直徑，f為質(zhì)量分?jǐn)?shù)，m為質(zhì)量流量，r為半徑。對于空氣采用速度進(jìn)口和壓力出口邊界條件；由于兩種燃燒模型采用的化學(xué)反應(yīng)機(jī)理相差較大，因此邊界條件的設(shè)置也有一些不同，基于氣相燃燒的渦耗散燃燒模型對燃油采用質(zhì)量進(jìn)口，而對于混合分?jǐn)?shù)PDF模型則采用離散相模型模擬燃油顆粒與空氣之間的相互作用。

表1 邊界條件設(shè)置

計(jì)算得到了RNGk-ε模型、Realizablek-ε模型、SSTk-ω模型三種湍流模型和渦耗散模型(EDM)、混合分?jǐn)?shù)PDF模型兩種燃燒模型6種組合方式下的NexGen燃燒器產(chǎn)生的火焰分布，下面通過指定位置處的溫度分布和平均熱流密度兩個(gè)參數(shù)對計(jì)算結(jié)果進(jìn)行比較分析。

6種湍流模型和燃燒模型組合下計(jì)算得到的NexGen燃燒器出口指定位置處的溫度分布如圖2所示。

由圖2可知，不同湍流模型和燃燒模型組合下燃燒器出口指定位置處的溫度分布有很大差別，但也有一些共同之處。共同之處有：①在高溫氣體浮升力的作用下，高溫區(qū)域(火焰)都上移了一定距離；②高溫區(qū)域形狀基本一致，都是寬度大于高度的類長方形結(jié)構(gòu)，且形狀與燃燒器出口形狀一致，說明湍流模型和燃燒模型對火焰的縱切面形狀幾乎沒有影響；③中間部分高溫區(qū)域溫度變化較緩慢，而在火焰邊緣溫度迅速下降，直至室溫。

差別之處是:有些模型是中心溫度高，隨著距中間距離的增加溫度逐漸減?。挥行┠Ｐ褪侵虚g溫度相對較低，但隨著距中間距離的增加，溫度先升至最高又迅速降低；而且溫度分布的對稱特性相差較大。此外，通過分析圖2中的6個(gè)溫度分布圖可知，采用渦耗散燃燒模型計(jì)算得到的最高溫度明顯高于混合分?jǐn)?shù)PDF模型，這與它們所采用的燃燒反應(yīng)機(jī)理有很大關(guān)系，渦耗散燃燒模型采用的是一步燃燒的化學(xué)反應(yīng)機(jī)理，未考慮中間的化學(xué)反應(yīng)，而這些多是吸熱反應(yīng)，因此計(jì)算的最高溫度偏高。

圖2 6種湍流模型和燃燒模型組合下的燃燒器出口指定位置處溫度分布

辛辛那提大學(xué)防火試驗(yàn)中心的研究人員通過實(shí)驗(yàn)得到了該面上的溫度分布矩陣[11]。我們在此基礎(chǔ)上通過MATLAB軟件得到了溫度分布云圖，如圖3所示。

圖3 實(shí)驗(yàn)測量的溫度分布云圖

通過比較圖3和圖2可知，Realizablek-ε湍流模型疊加混合分?jǐn)?shù)PDF燃燒模型計(jì)算得到的溫度分布與實(shí)驗(yàn)結(jié)果吻合最好，都是中間溫度較高，邊緣處溫度較低。

同時(shí)，辛辛那提大學(xué)也通過實(shí)驗(yàn)得到了指定位置處的平均熱流密度值為105.6 kW/m2，數(shù)值計(jì)算邊界條件中的空氣流量和燃油流量與實(shí)驗(yàn)一致，數(shù)值計(jì)算得到的不同湍流模型和燃燒模型組合下的平均熱流密度值如表2所示。

表2 6種湍流模型和燃燒模型計(jì)算的平均熱流密度 kW/m2

由表2可知，不同湍流模型和燃燒模型組合下產(chǎn)生的指定位置處的熱流密度值相差較大，其中Realizablek-ε湍流模型結(jié)合混合分?jǐn)?shù)PDF燃燒模型計(jì)算的平均熱流密度值與實(shí)驗(yàn)測量值最接近。

通過以上分析可知，Realizablek-ε湍流模型結(jié)合混合分?jǐn)?shù)PDF燃燒模型的計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合最好，最適用于NexGen燃燒器出口火焰分布的數(shù)值模擬。

2 發(fā)動機(jī)部件表面溫度響應(yīng)

為了進(jìn)一步研究發(fā)動機(jī)部件的防火能力，將Realizablek-ε湍流模型和混合分?jǐn)?shù)PDF燃燒模型計(jì)算得到的指定位置及其附近的熱氣流速度、溫度和熱流密度等參數(shù)導(dǎo)出，作為發(fā)動機(jī)部件的初始火焰場，然后將相應(yīng)的發(fā)動機(jī)部件置入該流場，研究該發(fā)動機(jī)部件上的溫度變化。

以2024鋁板材為例，研究在標(biāo)準(zhǔn)火焰下鋁板上的溫度分布[12]。將150 mm×150 mm的鋁平板部件置入標(biāo)準(zhǔn)火焰中，計(jì)算鋁平板部件表面的溫度分布，并以2024鋁板材的燒熔失效溫度作為鋁板失效的判定標(biāo)準(zhǔn)。計(jì)算得到的鋁平板部件上的溫度分布如圖4所示。

圖4 鋁平板表面溫度分布

由圖4可知，鋁板中心偏下位置處的溫度最高，達(dá)到了1 257.06 K，最低溫度為1 189 K，且隨著距板材中心位置距離的增加，溫度逐漸降低。根據(jù)董利穎[13]的實(shí)驗(yàn)結(jié)果可知，當(dāng)鋁板材背面溫度大于熱成像儀所能顯示的1 133 K 時(shí)，板材表面將形成可見空洞，屆時(shí)會有火焰穿過空洞流到板材背面，鋁板失效，無法完成其預(yù)定功能。說明在此邊界條件下鋁板將失效，但具體失效時(shí)間和失效模式仍需進(jìn)行進(jìn)一步研究。

3 結(jié)論

(1) 計(jì)算采用的湍流模型和燃燒模型對NexGen燃燒器出口火焰的數(shù)值計(jì)算結(jié)果有顯著影響。采用的湍流和燃燒模型不同，計(jì)算得到的燃燒器出口指定位置的溫度分布和熱流密度分布有很大差別。

(2) 建立了NexGen燃燒器的數(shù)值計(jì)算模型，即計(jì)算時(shí)采用Realizablek-ε湍流模型和混合分?jǐn)?shù)PDF燃燒模型，得到的NexGen燃燒器出口火焰分布與實(shí)驗(yàn)值吻合最好。

(3) 通過計(jì)算得到了標(biāo)準(zhǔn)火焰中鋁平板部件的溫度分布，鋁板中心偏下位置處溫度最高，最易被燒蝕。

(4) 建立了發(fā)動機(jī)部件防火試驗(yàn)的數(shù)值計(jì)算方法，即以計(jì)算得到的NexGen燃燒器出口火焰溫度分布和熱流密度分布及速度分布等作為初始邊界條件，建立初始火焰場，然后將需要開展防火審定的部件置入該火焰場中，研究發(fā)動機(jī)部件的溫度分布，當(dāng)最高溫度超過該部件材料的燒蝕溫度時(shí)，認(rèn)為該材料有燒蝕的風(fēng)險(xiǎn)，進(jìn)而可通過瞬態(tài)計(jì)算分析其達(dá)到燒蝕溫度的時(shí)間。