熱擋板設(shè)計參數(shù)對渦聲耦合壓力振蕩影響研究*

趙小鋒，孫兵兵，李軍強

(1 海軍裝備部， 北京 100000； 2 西安近代化學(xué)研究所， 西安 710065)

0 引言

隨著軍事需求對推力要求的不斷提高，分段式裝藥設(shè)計逐漸被大型固體火箭發(fā)動機采用,較好的完成了助推任務(wù)。但其彈道曲線或多或少存在壓力振蕩問題，嚴(yán)重時甚至影響發(fā)動機的正常工作，造成任務(wù)失敗，如Ariane 5、RSRM及Titan IV等[1]。為了研究該問題，國內(nèi)外許多學(xué)者做了大量的研究工作。通過Flandro[2]、Brown[3]、Anthoine[4-5]等人多年的研究，把該現(xiàn)象稱為“渦聲耦合”，解釋了出現(xiàn)壓力振蕩的機理，并得到了一些可靠有用的抑制措施。國內(nèi)的李強[6]、秦飛[7]、張翔宇[8]、王健儒[9]、許團(tuán)委[10]、蘇萬興[11]以及張嶠[12]等人也展開了研究，從理論計算與實驗等方面研究了該問題，為工程設(shè)計人員提供了一定的理論指導(dǎo)。但無論是實驗還是理論計算均缺乏必要的數(shù)據(jù)來指導(dǎo)熱擋板參數(shù)設(shè)計，不利于大型分段式發(fā)動機的設(shè)計。因此，文中基于阿麗亞娜5型運載火箭發(fā)動機的縮比模型,利用大渦模擬方法開展了數(shù)值研究，研究了熱擋板高度和彎曲角度等設(shè)計參數(shù)對渦聲耦合壓力振蕩的影響。

1 計算模型及邊界條件

文獻(xiàn)[5-12]計算結(jié)果表明，大渦模擬能夠精確計算發(fā)動機燃燒室內(nèi)由流動引起的壓力振蕩問題，因此采用大渦模擬來進(jìn)行數(shù)值計算研究，詳細(xì)的控制方程、亞格子模型及數(shù)值方法見文獻(xiàn)[11]。

以阿麗亞娜5型運載火箭發(fā)動機縮比模型的1/30作為研究對象，為了計算方便，計算模型采用二維軸對稱模型。

參照文獻(xiàn)[11]設(shè)置邊界條件，具體如圖1所示。為了準(zhǔn)確記錄燃燒室內(nèi)壓力的振頻與振幅特性，分別對發(fā)動機頭端部監(jiān)測點1(0.011 5，0.038)、熱擋板后端監(jiān)測點2(0.327 5，0.038)以及潛入式空腔內(nèi)監(jiān)測點3(0.366，0.05)處的振頻與振幅進(jìn)行了監(jiān)測與計算。

圖1 計算模型

計算網(wǎng)格的劃分在一定程度上可以影響計算結(jié)果的精度，蘇萬興在文獻(xiàn)[11]中研究了網(wǎng)格數(shù)目對計算結(jié)果的影響，得到了能夠準(zhǔn)確計算該模型所需的最佳網(wǎng)格數(shù)目，文中采取了與之相近的網(wǎng)格劃分方法，網(wǎng)格數(shù)目約為14萬。網(wǎng)格劃分如圖2所示。

圖2 熱擋板附近的網(wǎng)格劃分

2 計算結(jié)果與分析

2.1 擋板高度對壓力振蕩特性影響

文獻(xiàn)[11]研究了熱擋板位置對發(fā)動機內(nèi)的壓力振蕩的影響，發(fā)現(xiàn)合理設(shè)計熱擋板位置可以降低壓力振蕩的幅度，較好的指導(dǎo)了熱擋板的設(shè)計，但熱擋板的高度對火箭發(fā)動機壓力振蕩影響規(guī)律不是很清楚，因此探討擋板高度對發(fā)動機壓力振蕩特性的影響很有必要。為了更加直觀的觀察熱擋板高度的改變對壓力振蕩的影響，在圖1的基礎(chǔ)上新增熱擋板高度L為1 mm、3 mm 、5 mm、7 mm、11 mm的5組算例，其他條件與圖1保持一致，并結(jié)合標(biāo)準(zhǔn)模型(圖1L為9 mm)的計算結(jié)果，共6組算例進(jìn)行對比分析。

首先分析不同算例的時均流場。圖3顯示了不同熱擋板高度模型的速度流線圖?？梢郧宄乜吹?，隨著熱擋板高度增加，回流區(qū)范圍逐漸增大，渦的中心位置橫坐標(biāo)后移，縱坐標(biāo)下移。

圖3 不同熱擋板高度模型的速度流線圖

其次分析不同算例的渦量云圖。圖4顯示了不同熱擋板高度模型渦量云圖。從云圖對比圖中可以清楚地看到，隨著熱擋板高度L的變大，平均渦量也隨之變大，這意味著熱擋板高度的增加使得障礙物旋渦脫落的強度增加；同時從渦量云圖也可以看出漩渦在噴管前部的碰撞點位置下移，渦碰撞破裂后，一部分渦量進(jìn)入潛入式噴管空腔內(nèi)。

圖4 不同熱擋板高度模型的渦量云圖

對不同熱擋板高度算例監(jiān)測點1處計算得到的壓力數(shù)據(jù)進(jìn)行處理分析，可以發(fā)現(xiàn)，隨著熱擋板高度的增大，不同工況下壓力振蕩無量綱幅度逐漸增大，由熱擋板高度L=1 mm時的0.233%增長到L=11 mm時的1.6%，幾乎增加了一個數(shù)量級，如圖5所示。因此合理的擋板高度設(shè)計對減小壓力振蕩具有非常重要的意義。

圖5 不同熱擋板高度下的壓力振蕩無量綱幅度

從圖5也可以看出，3個監(jiān)測點振幅大小順序為監(jiān)測點3>監(jiān)測點1>監(jiān)測點2。其中監(jiān)測點2中壓力振幅最小，這是由于監(jiān)測點1和監(jiān)測點3處于發(fā)動機聲壓的波腹處，對壓力振蕩敏感，監(jiān)測點2處于發(fā)動機聲壓的波節(jié)與波腹之間，對壓力振蕩沒有監(jiān)測點1和監(jiān)測點3處敏感；此外監(jiān)測點3處的振蕩幅度大于監(jiān)測點1處，這是因為監(jiān)測點3除了處于聲壓的波腹處外，還處于潛入式噴管的潛入式空腔內(nèi)，潛入式空腔的加入會使得壓力振幅增大。 Anthoine[4-5]通過理論計算及實驗測試的方法研究驗證了該結(jié)果。此外，隨著熱擋板高度的增加，監(jiān)測點3與監(jiān)測點1處的振蕩幅度的差距逐漸增大，這從側(cè)面也證實了流入潛入式噴管空腔內(nèi)的渦量逐漸增多，較好的吻合了渦量云圖的結(jié)果。

2.2 擋板彎曲角度對壓力振蕩特性影響

擋板形狀設(shè)計也是分段裝藥設(shè)計的重要組成部分，因此合理設(shè)計熱擋板非常重要，但是相關(guān)方面的研究不是很充分。以圖1所示的基本模型為基礎(chǔ)，新增熱擋板彎曲角度(熱擋板與發(fā)動機壁面的夾角)為30°、45°、60°、75°四組算例，圖1所示模型的彎曲角度為90°，其他條件與圖1保持一致，并結(jié)合標(biāo)準(zhǔn)模型的計算結(jié)果，進(jìn)行對比分析。

首先分析不同算例的時均流場。圖6顯示了不同熱擋板彎曲角度模型的時均流場。可以清楚地發(fā)現(xiàn)，隨著熱擋板彎曲角度α的增大，熱擋板后緣下游回流區(qū)渦與潛入式空腔內(nèi)部的渦結(jié)構(gòu)聯(lián)系更加緊密。熱擋板上游附近形成的回流渦的范圍逐漸變小，當(dāng)α=90°時基本消失。

圖6 不同熱擋板彎曲角度模型的速度流線圖

對不同工況下監(jiān)測點1處的壓力振蕩數(shù)據(jù)進(jìn)行處理分析，不同彎曲角度下監(jiān)測點1處的無量綱壓力振蕩頻譜圖如圖7所示。

圖7 不同熱擋板彎曲角度下的壓力頻譜分布圖

從圖7可以看出，隨著熱擋板彎曲角度的增大，不同工況下壓力振蕩最激烈幅度處所對應(yīng)的頻率從高模態(tài)向低模態(tài)轉(zhuǎn)移，且振蕩現(xiàn)象逐漸從單一頻率下劇烈的振蕩向多頻率振蕩分散，降低了振蕩幅度，這主要是因為熱擋板的彎曲改變了聲渦耦合的反饋路徑，從而影響了渦脫落與聲場之間的耦合。在彎曲角度為90°時完成了頻率的轉(zhuǎn)移過程，壓力振蕩最激烈處由三階轉(zhuǎn)移到了一階，多頻率振蕩向基頻振蕩集中，壓力振蕩幅度再次增大。該結(jié)論較好的吻合了文獻(xiàn)[13]中Jerome冷流試驗得到的結(jié)論，因此合理的設(shè)計擋板彎曲角度可以把較大的單一頻率下的壓力振蕩幅度分散為多個頻率下較小的壓力振蕩幅度，從而降低單一頻率下劇烈振蕩對發(fā)動機的危害。

3 結(jié)論

文中基于阿麗亞娜5型運載火箭發(fā)動機的縮比模型,利用大渦模擬方法開展了數(shù)值研究，研究了熱擋板高度及彎曲角度對旋渦脫落的影響，結(jié)果表明：

1)隨著熱擋板高度的增加，渦在噴管前部碰撞后流入嵌入式空腔內(nèi)的渦量增多;壓力振蕩無量綱幅度逐漸增大。

2)隨著熱擋板彎曲角度的增大，不同工況下壓力振蕩最激烈幅度處所對應(yīng)的頻率從高模態(tài)向低模態(tài)轉(zhuǎn)移，且振蕩現(xiàn)象逐漸從單一頻率下劇烈的振蕩向多頻率振蕩分散，降低了最大壓力振蕩幅度，在彎曲角度為90°時完成了頻率的轉(zhuǎn)移過程，壓力振蕩最激烈處由三階轉(zhuǎn)移到了一階，多頻率振蕩向基頻振蕩集中，壓力振蕩幅度再次增大。

3)合理的設(shè)計熱擋高度與彎曲角度可以降低壓力振蕩的無量綱幅度。