高宇鵬，李俊杰 ，劉巍，李敬軒

1.北京動(dòng)力機(jī)械研究所，北京 100074 2.北京航空航天大學(xué) 宇航學(xué)院，北京 100191

0 引 言

燃燒不穩(wěn)定現(xiàn)象廣泛存在于火箭發(fā)動(dòng)機(jī)、航空發(fā)動(dòng)機(jī)以及大型空氣加熱器等裝置中，會(huì)使燃燒器內(nèi)產(chǎn)生不可控的巨大壓力脈動(dòng)，造成發(fā)動(dòng)機(jī)回火、熄火，燃燒器壁面過(guò)熱，導(dǎo)致發(fā)動(dòng)機(jī)壽命縮短甚至損毀等一系列問(wèn)題。20世紀(jì)60年代，在對(duì)F-1火箭發(fā)動(dòng)機(jī)初始方案的44次試車中，超過(guò)半數(shù)的試車都發(fā)生了燃燒不穩(wěn)定現(xiàn)象；其振幅甚至達(dá)到了燃燒室平均壓力的100%。因此，為盡可能減少不穩(wěn)定燃燒帶來(lái)的巨大危害，通常采用實(shí)驗(yàn)和數(shù)值仿真等方法對(duì)燃燒不穩(wěn)定性機(jī)理進(jìn)行研究。

王迪等對(duì)單噴嘴模型發(fā)動(dòng)機(jī)進(jìn)行了實(shí)驗(yàn)研究，測(cè)量了不同燃燒室和噴嘴縮進(jìn)長(zhǎng)度對(duì)燃燒穩(wěn)定性的影響，結(jié)果發(fā)現(xiàn)增加縮進(jìn)長(zhǎng)度會(huì)對(duì)縱向高頻燃燒不穩(wěn)定性產(chǎn)生阻尼作用，燃燒室長(zhǎng)度對(duì)縱向高頻燃燒不穩(wěn)定性的影響比縮進(jìn)長(zhǎng)度更為明顯。Richecoeur等對(duì)液氧/氣態(tài)甲烷多噴嘴模型燃燒室進(jìn)行了試驗(yàn)研究，結(jié)果發(fā)現(xiàn)橫向壓力震蕩與燃燒和火焰之間存在強(qiáng)烈的耦合作用，可增大火焰膨脹率以及火焰向壁面的傳熱率。Bibik等拍攝了燃燒室內(nèi)部氣流的旋轉(zhuǎn)過(guò)程，發(fā)現(xiàn)氣流旋轉(zhuǎn)方向與中心噴嘴旋轉(zhuǎn)相同，認(rèn)為中心噴嘴是激發(fā)切向不穩(wěn)定的重要因素。Méry等實(shí)驗(yàn)研究了液氧/甲烷燃燒室在存在擾動(dòng)的條件下，聲學(xué)擾動(dòng)與火焰之間的相互作用。結(jié)果表明：在中等擾動(dòng)條件下，液體射流存在橫向周期性振蕩；大幅度擾動(dòng)時(shí)，火焰隨聲學(xué)擾動(dòng)而周期性振蕩。Sliphorst等通過(guò)實(shí)驗(yàn)研究發(fā)現(xiàn)，燃燒不穩(wěn)定激勵(lì)主要來(lái)自于火焰與聲學(xué)振蕩的耦合作用。

Urbano等通過(guò)大渦模擬（Large Eddy Simulation, LES）對(duì)功率為80 MW的全尺寸燃燒器進(jìn)行仿真，結(jié)果發(fā)現(xiàn)一階切向和一階縱向不穩(wěn)定對(duì)燃燒室的不穩(wěn)定性起重要作用，并且與混合燃料的噴注過(guò)程強(qiáng)烈耦合。Shimizu等等對(duì)液體火箭發(fā)動(dòng)機(jī)（Liquid Rocket Engine, LRE）進(jìn)行了切向數(shù)值仿真，認(rèn)為平均流動(dòng)分布控制一階切向不穩(wěn)定的發(fā)生。Staffelbach等對(duì)自激振蕩的環(huán)形燃燒室進(jìn)行了大渦模擬，成功捕捉了燃燒室自激振蕩頻率，發(fā)現(xiàn)其特征是由兩種不同振幅的旋轉(zhuǎn)模式疊加而成。王方、張志浩和夏朝陽(yáng)等均采用LES對(duì)發(fā)動(dòng)機(jī)內(nèi)燃燒不穩(wěn)定性進(jìn)行詳細(xì)的研究。Li等對(duì)ORACLESL長(zhǎng)火焰燃燒室進(jìn)行了預(yù)測(cè)。采用大渦模擬獲得系統(tǒng)的火焰?zhèn)鬟f函數(shù)，并將此傳遞函數(shù)帶入到Li和Morgans等開(kāi)發(fā)的燃燒不穩(wěn)定預(yù)測(cè)和模擬程序—OSCILOS中，最終得到了精準(zhǔn)的熱–聲振蕩頻率和幅值的預(yù)測(cè)值。

綜上所述，前人已通過(guò)實(shí)驗(yàn)和數(shù)值的方法對(duì)熱聲耦合不穩(wěn)定燃燒進(jìn)行了研究，但是大尺寸燃燒器燃燒不穩(wěn)定性的預(yù)測(cè)一直存在計(jì)算量太大的問(wèn)題。本文將通過(guò)熱–聲解耦的方法，對(duì)某型空氣加熱器燃燒不穩(wěn)定性進(jìn)行預(yù)測(cè)，為加熱器熱試提供技術(shù)支撐。

1 燃燒-聲學(xué)解耦方法

本文基于Li等的方法，將燃燒對(duì)擾動(dòng)的響應(yīng)計(jì)算與聲學(xué)計(jì)算解耦，即分別計(jì)算火焰對(duì)聲學(xué)擾動(dòng)的響應(yīng)（通常表征為火焰?zhèn)鬟f函數(shù)）和聲學(xué)系統(tǒng)在熱源擾動(dòng)下的響應(yīng)。因此，對(duì)于燃燒的數(shù)值模擬只需考慮燃燒室內(nèi)的火焰對(duì)來(lái)流擾動(dòng)的響應(yīng)，不需要考慮包括集氣環(huán)和噴管在內(nèi)的整個(gè)發(fā)動(dòng)機(jī)系統(tǒng)，顯著降低了計(jì)算量，提高了計(jì)算速度。

1.1 噴管聲學(xué)邊界條件計(jì)算

在聲學(xué)邊界條件計(jì)算中，除了噴管外，都可處理為壅塞或聲學(xué)閉口。由于噴管復(fù)雜的頻率特性，需要考慮其聲學(xué)邊界條件隨頻率的變化，通常表征為傳遞函數(shù)的形式。一般拉瓦爾噴管的聲學(xué)系統(tǒng)可以表示成如圖1所示模型。當(dāng)噴管上游來(lái)流的物性參數(shù)發(fā)生擾動(dòng)時(shí)，下游氣流會(huì)產(chǎn)生壓力與速度的振蕩，產(chǎn)

圖1 噴管的聲學(xué)模型Fig.1 The acoustic model of nozzle

對(duì)于亞聲速流噴管（圖1），將向下游傳播的聲波施加在入口處，向上游傳播的聲波施加在出口處。假設(shè)噴管中的氣體是理想氣體，建立一維線性化歐拉方程如下：

式中，u為氣體速度，p為氣體壓力，為比熱比，s為氣體熵，為 氣體定壓比熱容，為當(dāng)?shù)芈曀伲蠘?biāo)“”為擾動(dòng)量。

將上述控制方程進(jìn)行無(wú)量綱化處理，當(dāng)馬赫數(shù)不為1時(shí)，進(jìn)行傅里葉變換可以得到：

式中，=/L，為 噴管長(zhǎng)度，x為距噴管入口的軸向距離；矩陣是由不變量構(gòu)成的矢量。矩陣和的形式分別如下：

利用Magnus展開(kāi)法可以求得式（2）的解，并求得最終的聲學(xué)反射系數(shù)。

1.2 數(shù)值計(jì)算方法

1.2.1 燃燒室三維建模

加熱器整體模型如圖2所示，包含噴注面板、燃燒室、噴管等部件。其中，噴注面板上包含48個(gè)燃料噴嘴和36個(gè)冷卻空氣噴嘴。燃燒室直徑為350 mm，長(zhǎng)度為800 mm。在對(duì)比了考慮和不考慮集氣環(huán)的燃燒室聲學(xué)特性后，發(fā)現(xiàn)集氣環(huán)對(duì)燃燒室不穩(wěn)定模態(tài)的穩(wěn)定性影響基本可以忽略。且由于各種連接件的壅塞特性，因此，在燃燒不穩(wěn)定性計(jì)算中忽略了集氣環(huán)這一部件。同時(shí)，由于噴注面板上噴嘴分布的高度對(duì)稱性，將采用1/6模型進(jìn)行計(jì)算，如圖3、4所示。

圖2 加熱器各部件示意圖Fig.2 Schematic diagram of the heater components

圖3 燃燒室模型示意圖Fig.3 Schematic diagram of combustion chamber model

為確?；鹧嫦嚓P(guān)參數(shù)計(jì)算準(zhǔn)確，對(duì)噴注面板下游200 mm范圍內(nèi)的網(wǎng)格進(jìn)行加密，其各向尺寸小于1 mm。經(jīng)驗(yàn)證計(jì)算，在網(wǎng)格量達(dá)到265萬(wàn)時(shí)，滿足網(wǎng)格無(wú)關(guān)性要求；瞬態(tài)計(jì)算時(shí)間步長(zhǎng)為10s時(shí)，滿足時(shí)間無(wú)關(guān)性要求，此時(shí)全局最大Courant數(shù)小于5。

1.2.2 模型選擇和參數(shù)設(shè)置

由于燃燒室內(nèi)總溫為1 743 K，總壓達(dá)8.2 MPa，當(dāng)預(yù)混燃料采用乙醇/液氧時(shí)，二者均處于超臨界態(tài)。超臨界流體流出噴嘴后，其氣液界面較為模糊，密度為液態(tài)的1/3，而黏度只有液態(tài)的1/12～1/4，因此其流動(dòng)行為更偏向于氣體。在數(shù)值仿真中，燃料按連續(xù)相進(jìn)行處理。

圖4 燃燒室網(wǎng)格劃分示意圖Fig.4 Divided mesh of combustion chamber model

采用ANSYS Fluent對(duì)燃燒室內(nèi)的燃燒過(guò)程進(jìn)行數(shù)值仿真，其中內(nèi)部湍流流動(dòng)采用Realizable模型求解；燃燒采用Eddy-Disspation模型求解。該燃燒模型假定燃燒速率遠(yuǎn)大于燃料混合速率，即燃料一經(jīng)混合便完全燃燒，忽略復(fù)雜的化學(xué)反應(yīng)速率，采用乙醇一步總包反應(yīng)機(jī)理：

為確保在高溫高壓的燃燒室環(huán)境下，燃料物性參數(shù)（特別是對(duì)當(dāng)?shù)芈曀儆绊憳O大的密度）數(shù)值準(zhǔn)確，采用Peng-Robinson狀態(tài)方程計(jì)算。燃料入口參數(shù)如表1所示，各噴嘴入口邊界條件采用質(zhì)量流量入口，燃燒室出口邊界條件則采用壓力出口。

表1 預(yù)混燃料入口工況Table 1 The working conditions of premixed fuel at the inlet

得到燃燒室穩(wěn)態(tài)場(chǎng)后，在燃料入口處施加不同頻率、振幅的正弦擾動(dòng)（式（6）），即可獲得各噴嘴的熱釋放率隨時(shí)間的變化關(guān)系，進(jìn)而得到各噴嘴的火焰?zhèn)鬟f函數(shù)：

1.2.3 火焰?zhèn)鬟f函數(shù)的定義

通過(guò)對(duì)每個(gè)噴嘴出口的火焰熱釋放率進(jìn)行處理，提取激振頻率上的擾動(dòng)幅值和相位，并帶入火焰?zhèn)鬟f函數(shù)的定義式（式（7）），即可得到各火焰在不同頻率和不同擾動(dòng)幅值下的火焰?zhèn)鬟f函數(shù)。

1.3 聲學(xué)模型

1.3.1 聲學(xué)模擬模型

對(duì)整個(gè)燃燒室建立Helmholtz方程并耦合火焰?zhèn)鬟f函數(shù)，采用噴管聲學(xué)邊界條件，即可求得整個(gè)系統(tǒng)的聲學(xué)模態(tài)。

對(duì)于一個(gè)有熱源的區(qū)域，流體的密度隨壓力和 熵的變化而變化，即：

結(jié)合火焰?zhèn)鬟f函數(shù)的定義式（7），最終可推導(dǎo)得到單極域源的表達(dá)式為：

式中，下標(biāo)b和u分別代表已燃?xì)怏w和未燃?xì)怏w，S為噴嘴截面積，V為熱源體積。

1.3.2 聲學(xué)邊界阻抗模型

在實(shí)際工況中，聲波在出口處存在一定的反射，即為聲阻抗R。如圖5所示，設(shè)入射波強(qiáng)度為A，反射波強(qiáng)度為A。

圖5 聲波反射與聲學(xué)邊界阻抗示意圖Fig.5 Schematic diagram of sound reflection and acoustic boundary impedance

當(dāng)集氣環(huán)與噴注器和燃燒室解耦后，噴注器進(jìn)口可認(rèn)為是壓力節(jié)點(diǎn)，空氣燃料混合噴嘴入口可看作是硬聲場(chǎng)邊界，空氣噴嘴入口為軟聲場(chǎng)邊界。燃燒室出口與噴管相連，其聲學(xué)邊界條件采用1.1小節(jié)中推導(dǎo)得到的聲學(xué)反射系數(shù)模型。

聲阻抗定義為：

2 燃燒室不穩(wěn)定模態(tài)估計(jì)

由于空氣加熱器燃燒室的噴嘴面積總和與燃燒室橫截面積之比遠(yuǎn)遠(yuǎn)小于1，且噴嘴內(nèi)與燃燒室的溫度相差較大（約1 500 K），因此噴嘴前的聲學(xué)系統(tǒng)可與燃燒室內(nèi)的聲學(xué)系統(tǒng)解耦。燃燒室出口為一壅塞的噴管，燃燒室內(nèi)的聲學(xué)系統(tǒng)也可以與噴管下游的聲學(xué)系統(tǒng)解耦，故只需計(jì)算燃燒室內(nèi)的聲學(xué)模態(tài)。在粗估燃燒室聲學(xué)模態(tài)的頻率時(shí)，燃燒室進(jìn)口可認(rèn)為是聲學(xué)閉口條件，出口由于噴管前馬赫數(shù)小于0.1,也可近似認(rèn)為是聲學(xué)閉口。

燃燒室的各個(gè)聲學(xué)模態(tài)的頻率f可以根據(jù)經(jīng)典公式進(jìn)行粗估：

式中，、、1、、、、和分別為縱向模態(tài)數(shù)、切向模態(tài)數(shù)、徑向模態(tài)數(shù)、燃燒室聲速、徑向波數(shù)、縱向波數(shù)、垂直于縱向方向的波數(shù)和燃燒室的直徑。為 貝塞爾函數(shù) dJ(k)/d=0的根，可通過(guò)查表2得到。

表2 不同模態(tài)下 k m,j 的值Table 2 The value of k m,j in different modes

結(jié)合加熱器燃燒室結(jié)構(gòu)參數(shù)以及燃燒室內(nèi)部聲速（約為837 m/s），可以粗估燃燒室的各聲學(xué)不穩(wěn)定模態(tài)頻率，如表3所示。

表3 燃燒室各不穩(wěn)定模態(tài)頻率估計(jì)值Table 3 Estimated frequency of each unstable mode of the combustion chamber

考慮到實(shí)際加熱器燃燒室最易出現(xiàn)（n=1，m=1）的不穩(wěn)定模態(tài)，因此最可能出現(xiàn)的不穩(wěn)定頻率應(yīng)為（n=1，m=1，j=1）模態(tài)的頻率1 406 Hz。這一粗估頻率將作為后續(xù)計(jì)算的初始參考值。

3 結(jié)果與討論

3.1 聲學(xué)反射系數(shù)計(jì)算結(jié)果

對(duì)1.1小節(jié)所得的數(shù)學(xué)方法進(jìn)行求解，可得0～3 000 Hz頻率范圍內(nèi)噴管聲學(xué)反射系數(shù)的幅頻、相頻特性（圖6藍(lán)色圈）。對(duì)計(jì)算結(jié)果進(jìn)行五階多項(xiàng)式擬合，可求得噴管聲學(xué)反射系數(shù)的近似解析表達(dá)式，其頻響特性曲線如圖6中紅色實(shí)線所示。此結(jié)果可作為燃燒室聲學(xué)計(jì)算的出口邊界條件，用于計(jì)算燃燒室的各聲學(xué)模態(tài)。

圖6 聲學(xué)反射系數(shù)的幅頻、相頻特性曲線Fig.6 Amplitude frequency and phase frequency characteristic curve of acoustic reflection coefficient

3.2 燃燒室內(nèi)數(shù)值計(jì)算結(jié)果

3.2.1 溫度場(chǎng)分析

由于粗估燃燒室最不穩(wěn)定頻率為1 406 Hz，因此取擾動(dòng)頻率f =1 400 Hz、擾動(dòng)振幅=04的工況進(jìn)行分析。圖7為t=0.021 161 s時(shí)刻，燃燒室距噴注面板50、100、200、400 mm處截面的總溫云圖。

圖7 不同截面處的總溫云圖Fig.7 Total temperature contour map at different cross-sections

從圖7中可以看到，燃料噴嘴出口近場(chǎng)截面的溫度分布并不均勻，中心噴嘴受周圍冷卻氣流的影響更為明顯，存在明顯的低溫環(huán)（圖7中黑色區(qū)域）；隨著軸向距離的增加，冷熱空氣逐步混合，在距噴注面板約400 mm處幾乎混合均勻，燃燒室內(nèi)部溫度分布趨于穩(wěn)定。

3.2.2 熱釋放率分析

圖8是t=0.021 161 s、f=1 400 Hz、3種不同擾動(dòng)振幅下，燃燒室中心面上的火焰熱釋放率鋒面與流線圖。其中，熱釋放率鋒面可近似認(rèn)為是火焰鋒面。在施加擾動(dòng)后，燃燒室的中心與邊緣靠近壁面處會(huì)分別產(chǎn)生兩個(gè)較為強(qiáng)烈的渦旋。處于中心位置的渦旋對(duì)火焰鋒面的拉伸作用明顯。在擾動(dòng)振幅相對(duì)較小時(shí)，渦旋強(qiáng)度較低（流線相對(duì)稀疏平整），火焰鋒面平整，形狀較為舒展；隨擾動(dòng)振幅增加，中心渦旋強(qiáng)度增加（流線更為緊湊卷曲），對(duì)火焰鋒面的影響愈發(fā)明顯，火焰產(chǎn)生明顯褶皺，且長(zhǎng)度變短。邊緣渦旋對(duì)于流量脈動(dòng)幅值的響應(yīng)并不是很明顯。在高頻擾動(dòng)下，中心渦旋造成的速度擾動(dòng)直接對(duì)火焰產(chǎn)生擾動(dòng)，造成火焰面周期性卷曲和拉伸，從而產(chǎn)生火焰?zhèn)鞑ニ俣葦_動(dòng)，引起火焰熱釋放率的波動(dòng)。

本文分別檢測(cè)了每一個(gè)噴嘴的熱釋放率隨時(shí)間的變化，以圖8中標(biāo)注的噴嘴為例，特別對(duì)受到冷空氣擠壓的火焰進(jìn)行討論。圖9（a）、（b）分別為該噴嘴在擾動(dòng)振幅=04、不同擾動(dòng)頻率f下和擾動(dòng)頻率f1 400 Hz、不同振幅下，火焰熱釋放率隨時(shí)間變化的曲線。

圖8 不同振幅下的燃燒室內(nèi)熱釋放率云圖及流線圖Fig.8 Contour of heat release rate and streamline inside the combustion chamber under different amplitudes

從圖9（a）可以看到，當(dāng)振幅相同時(shí)，該噴嘴出口火焰熱釋放率對(duì)頻率的響應(yīng)不同—不同擾動(dòng)頻率對(duì)應(yīng)不同的熱釋放率脈動(dòng)幅值。理論上，當(dāng)對(duì)入口流量給定相同的擾動(dòng)振幅時(shí)，擾動(dòng)頻率越高，熱釋放脈動(dòng)幅值應(yīng)更低。但對(duì)于該噴嘴，在擾動(dòng)頻率f=1 400 Hz時(shí)，其幅值相對(duì)f=1 250 Hz并未出現(xiàn)明顯的下降，即燃燒室內(nèi)部的不穩(wěn)定頻率可能出現(xiàn)在f=1 400 Hz附近，這一結(jié)果也證明了第2節(jié)的粗估結(jié)果是可靠的。

從圖9（b）可以看到，當(dāng)擾動(dòng)頻率f=1 400 Hz時(shí)，該噴嘴出口火焰對(duì)于不同振幅的響應(yīng)同樣存在差異。隨擾動(dòng)振幅的增加，熱釋放脈動(dòng)幅值增加；但熱釋放率脈動(dòng)幅值并不隨流量擾動(dòng)幅值的增加而持續(xù)增加。當(dāng)=06 時(shí)，熱釋放脈動(dòng)幅值與=04 時(shí)并無(wú)明顯差異，這也說(shuō)明了該噴嘴出口火焰對(duì)于大流量脈動(dòng)的響應(yīng)并不敏感。

圖9 火焰熱釋放率隨時(shí)間變化曲線Fig.9 Dimensionless flame heat release rate vs.time curve

3.3 燃燒室內(nèi)聲學(xué)模態(tài)預(yù)測(cè)結(jié)果

圖10為示例噴嘴（圖8）對(duì)應(yīng)的火焰?zhèn)鬟f函數(shù)幅值和相位圖。從圖中可以看到，在任意擾動(dòng)頻率和擾動(dòng)振幅下，該噴嘴的火焰?zhèn)鬟f函數(shù)幅值均小于1，即火焰熱釋放率相對(duì)于入口流量擾動(dòng)的響應(yīng)呈抑制作用。并且隨著擾動(dòng)流量的增加，傳遞函數(shù)幅值在各頻率下均有明顯的下降趨勢(shì)。

圖10 火焰?zhèn)鬟f函數(shù)的幅值和相位Fig.10 The amplitude and phase of the flame transfer function

將各噴嘴的火焰?zhèn)鬟f函數(shù)結(jié)果代入1.3小節(jié)的聲學(xué)模型中，即可對(duì)存在流量擾動(dòng)的燃燒室的穩(wěn)定性進(jìn)行預(yù)測(cè)，預(yù)測(cè)結(jié)果如圖11所示。該圖為燃燒室無(wú)量綱（n=1，m=1，j=1）模態(tài)的壓力分布。壓力分布模式主要存在2種，分布在沿燃燒室直徑方向上互相垂直。經(jīng)1.3小節(jié)中的模型計(jì)算，當(dāng)擾動(dòng)振幅分別為 02、04、06時(shí)，（n=1，m=1，j=1）模態(tài)的特征值分別為1 389.9+5.975 i、1 389.9+6.060 3 i和1 389.9+6.142 4 i。該特征值為復(fù)數(shù)，其實(shí)部為模態(tài)頻率，虛部為增長(zhǎng)率。當(dāng)虛部大于0時(shí)，代表燃燒室在該聲學(xué)模態(tài)上是穩(wěn)定的；當(dāng)虛部小于0時(shí)，代表燃燒室在該聲學(xué)模態(tài)上是不穩(wěn)定的。

圖11 不同擾動(dòng)振幅下，燃燒室最不穩(wěn)定模態(tài)的總聲壓場(chǎng)Fig.11 The total sound pressure field of the most unstable mode in the combustion chamber

結(jié)果表明，在不同擾動(dòng)幅值下，燃燒室的（n=1，m=1，j=1）模態(tài)均穩(wěn)定。這與該型加熱器的試驗(yàn)結(jié)果一致。該型加熱器在后期試驗(yàn)中并未出現(xiàn)燃燒不穩(wěn)定現(xiàn)象。另一型加熱器（燃燒室直徑為460 mm，長(zhǎng)度為800 mm）在試驗(yàn)過(guò)程中出現(xiàn)燃燒不穩(wěn)定現(xiàn)象，振蕩頻率為1 100 Hz。基于第2節(jié)中的方法對(duì)該燃燒室的不穩(wěn)定模態(tài)進(jìn)行預(yù)測(cè)，得到振蕩頻率為1 135 Hz，與試驗(yàn)測(cè)量頻率符合很好，進(jìn)一步驗(yàn)證了本文方法良好的預(yù)測(cè)能力。

4 結(jié) 論

1）采用Magnus展開(kāi)法獲得了噴管的聲學(xué)邊界條件。通過(guò)五階多項(xiàng)式擬合獲得了噴管聲反射系數(shù)的近似解析表達(dá)式。

2）采用ANSYS Fluent對(duì)全尺寸燃燒器進(jìn)行了非穩(wěn)態(tài)數(shù)值仿真，得到了燃燒室內(nèi)部存在擾動(dòng)時(shí)的物理場(chǎng)信息，發(fā)現(xiàn)位于噴注面板中心處的噴嘴受到周圍冷卻空氣的影響較大。提取了每一個(gè)噴嘴的熱釋放率信息，獲得了噴嘴熱釋放率脈動(dòng)與速度擾動(dòng)之間的火焰?zhèn)鬟f函數(shù)。

3）將所得聲學(xué)邊界條件與火焰?zhèn)鬟f函數(shù)代入燃燒室聲學(xué)模型中，成功預(yù)測(cè)該燃燒室最易發(fā)生不穩(wěn)定燃燒的頻率為1 389.9 Hz，且不同擾動(dòng)振幅下，燃燒室的（n=1，m=1，j=1）模態(tài)均穩(wěn)定，不會(huì)發(fā)生不穩(wěn)定燃燒現(xiàn)象，與試驗(yàn)結(jié)果一致。