李冬冬，王 革

(哈爾濱工程大學 航天與建筑工程學院，哈爾濱 150001)

二次噴射閥門作動對側(cè)向力的影響

李冬冬，王 革

(哈爾濱工程大學 航天與建筑工程學院，哈爾濱 150001)

針對潛入式噴管燃氣二次噴射推力矢量控制系統(tǒng)，利用AUSM+格式，并結(jié)合Realizablek-ε湍流模型，對噴管內(nèi)二次噴射流場進行數(shù)值模擬。主要研究二次噴射閥門作動過程中側(cè)向力的變化規(guī)律，分析正弦波、三角波2種作動方式和5、7、10 Hz 3種頻率對側(cè)向力的影響規(guī)律。結(jié)果表明，在閥動過程中側(cè)向力變化與閥門開度變化趨勢相同，大體上側(cè)向力隨閥門開度增加而增加，隨閥門開度減小而減小，但存在一定的時間和空間滯后，側(cè)向力變化滯后于閥門位置變化，最大側(cè)向力出現(xiàn)在半個周期以后，即最大開度出現(xiàn)時刻之后，且閥門重新閉合后，很短時間內(nèi)仍有側(cè)向力殘存，相同開度下，側(cè)向力有如下大小關(guān)系：關(guān)閉過程大于準穩(wěn)態(tài)，二者又大于開啟過程中的側(cè)向力；從最大側(cè)向力出現(xiàn)時刻分析滯后時間τD，在波形上三角波下滯后時間τD要大于正弦波下的τD，在頻率上10 Hz下的τD最大，5 Hz下的τD最小;從側(cè)向力相對殘存時間τR分析時間滯后，在波形上三角波下的τR要大于正弦波下τR，在頻率上7 Hz下的τR最大，5 Hz下的τR最小。

燃氣二次噴射；推力矢量控制；數(shù)值模擬；側(cè)向力

0 引言

推力矢量控制技術(shù)是飛行器關(guān)鍵技術(shù)之一[1]，受到世界各國的普遍關(guān)注。推力矢量控制技術(shù)的發(fā)展與應用，可有效提高飛行器機動性能，減小飛行器阻力和自身重量，提高隱形能力、安全性能和生存能力，顯著提高作戰(zhàn)效能。

與機械式推力矢量控制系統(tǒng)相比，燃氣二次噴射系統(tǒng)具有頻率響應快，靈敏度高、效率高、密封性好，結(jié)構(gòu)簡單，伺服系統(tǒng)功率要求小，系統(tǒng)可靠性高等優(yōu)點，已在少數(shù)國家的裝備中使用。

美國針對燃氣二次噴射的研究，開發(fā)了PAB3D、TetrUSS、NPARC等仿真系列軟件。Waithe K A 等[2]使用PAB3D對二次噴射推力矢量控制進行了計算，并與試驗進行對比，驗證了PAB3D對復雜流場的計算能力。Lee J G 等[3]對燃氣二次噴射的氣固兩相流問題進行了數(shù)值模擬，研究了質(zhì)量流率、噴射位置等對推力矢量控制系統(tǒng)的影響。Newton J F 等[4]通過對大推力長時間工作的固體火箭發(fā)動機進行試驗研究，分析了側(cè)向控制力的產(chǎn)生機理和壓強擾動結(jié)構(gòu)。

與國外相比，國內(nèi)在此領(lǐng)域的研究雖然起步較晚，但有關(guān)專家學者已經(jīng)開展了大量有意義的研究工作。喬渭陽[5]、劉剛[6]、鄧遠擷[7]等基于二維N-S方程，模擬了二次噴射推力矢量噴管內(nèi)流場。研究表明，二次噴射系統(tǒng)可在低推力損失的情況下，有效地提供側(cè)向控制力。賴川等[8]通過PHOENICS程序，對氣體二次噴射進行了三維流場數(shù)值仿真，并對二維矩形噴管進行了等效實體和風洞試驗[9]，研究了噴射參數(shù)對弓形激波的影響。劉輝等[10]通過對三維流場的模擬，比較了不同噴射參數(shù)和不同噴管落壓比下的流場特征，分析了這些參數(shù)對矢量偏轉(zhuǎn)效率和推力系數(shù)的影響。曹熙煒等[11]利用計算流體力學(CFD)軟件，對不同溫度的二次流進行了模擬仿真分析，研究了二次流溫度對推力矢量的影響規(guī)律。吳雄等[12-14]對該問題進行了充分的實驗和理論探索，對二次噴射的流場、響應面優(yōu)化設計及閥門開啟的非定常效應等問題做了深入分析。王革、肖雪峰等[15-16]通過對雙孔燃氣二次噴射流場的研究，分析了噴射孔徑、噴射角、噴射位置、噴射孔夾角及噴射孔包抄角對側(cè)向控制力、軸向推力的影響規(guī)律與強度。

本文基于國內(nèi)研究現(xiàn)狀及發(fā)展需求，針對潛入式噴管燃氣二次噴射方案，結(jié)合目前的CFD技術(shù)，對噴管燃氣二次噴射進行更深入的數(shù)值研究，從系統(tǒng)動態(tài)作動的角度，研究二次噴射閥門作動方式及作動頻率對側(cè)向控制力的影響，探究側(cè)向控制力在閥門作動過程中的變化規(guī)律及頻率和作動方式上的差異性，為進一步的實驗研究和工程設計提供一定的參考。

1 物理問題的描述和建立

1.1 幾何模型

本文以潛入式噴管作為研究對象，采用具有一定結(jié)構(gòu)優(yōu)勢的燃氣二次噴射方案，通過連接燃燒室潛入段與噴管擴張段的直通道，將燃燒室內(nèi)的燃氣直接引入主噴管擴張段，形成二次噴射。

圖1給出了燃氣二次噴射結(jié)構(gòu)的示意圖。為了節(jié)省計算資源，計算模型中對燃燒室段的結(jié)構(gòu)不予考慮，圖2給出了計算使用模型的三維透視圖。噴射角α=65°，噴射孔與喉部直徑比值為0.277，噴管擴張比為16.90。

1.2 網(wǎng)格劃分

計算網(wǎng)格的劃分采用混合網(wǎng)格方式，在結(jié)構(gòu)形狀復雜的燃燒室內(nèi)噴射閥周圍采用非結(jié)構(gòu)化網(wǎng)格，其余部分采用結(jié)構(gòu)化網(wǎng)格，邊界層指標y+控制在30左右。

網(wǎng)格劃分總體效果如圖3所示。

1.3 湍流模型及邊界條件

考慮到流動中有射流及射流流動引起的復雜流動情況，湍流模型使用realizablek-ε模型[17]。

計算中，給定噴管進口定流量邊界條件，不考慮閥動過程中開孔后二次射流對進口邊界的反饋作用。噴管出口邊界條件為壓力出口，給定1 atm。其余為壁面邊界條件。

1.4 計算格式

本文為了得到較清晰的激波結(jié)構(gòu)，并由激波結(jié)構(gòu)判斷結(jié)果的可靠性，采用AUSM+計算格式，AUSM+具體格式見式(1)～式(4)：

(1)

馬赫數(shù)分裂函數(shù)：

(2)

壓力項分裂函數(shù)：

(3)

界面音速的確定：

(4)

2 計算結(jié)果分析

本文針對正弦波和三角波兩種作動方式，10、7、5 Hz 3種作動頻率，共計6種情況進行了計算分析,目的是研究閥動過程中側(cè)向力的變化及作動頻率和作動方式對側(cè)向力的影響。計算過程中，閥門初始開度為0 mm，前半個周期閥門由0 mm開度向最大18 mm開度運動，后半個周期閥門由18 mm開度向0 mm開度運動，而后閥門位置保持0 mm開度不變。

對于正弦波作動方式，其具體運動方程如下：

對于三角波作動方式，其具體運動方程如下：

式中v為閥門運動速度，mm/s；x為閥門開度，mm；f為閥門作動頻率，Hz；T為閥門作動周期，s；A為閥門最大開度，A=18 mm。

2.1 典型工作狀態(tài)分析

本節(jié)選取正弦波5 Hz工作狀態(tài)，對閥動過程中側(cè)向力進行分析。依據(jù)上節(jié)對閥門作動方式的描述，正弦波5 Hz工作狀態(tài)閥動過程中，閥門速度和閥門開度變化如圖4所示。

圖5是閥動過程中捕獲的典型流場馬赫數(shù)等值線圖與近壁面馬赫數(shù)分布。由馬赫數(shù)等值線圖可清晰地觀察到，二次射流與超聲速主流在主噴管擴張段相互作用產(chǎn)生的復雜波系結(jié)構(gòu)。首先，在噴射孔的上游形成了一道弓形激波，在噴射孔前端與后端各形成一道激波，且這兩道激波在下游某處相交，形成兩道新的相對較弱的激波，前端激波相交后的激波在壁面處又反射出一道激波，另一道相交激波一直延伸到噴管出口。近壁面馬赫數(shù)的變化也反映出同樣的波系結(jié)構(gòu)。

圖6是閥動過程中閥門開度及側(cè)向力隨時間的變化。側(cè)向力的變化趨勢與閥門開度的變化趨勢基本一致，側(cè)向力最大值大致出現(xiàn)在0.104 s處，滯后于閥門開度為18 mm時的0.100 s。

表1列出了動態(tài)過程中最大側(cè)向力及18 mm開度時側(cè)向力與準穩(wěn)態(tài)18 mm側(cè)向力的大小。動態(tài)過程中最大側(cè)向力與準穩(wěn)態(tài)下最大側(cè)向力僅相差0.202%，在動態(tài)過程中的18 mm開度與準穩(wěn)態(tài)18 mm時僅相差-0.405%。此外，在閥門重新關(guān)閉后，直到0.208 s，側(cè)向力才恢復到初始時刻的大小。這是因為閥門關(guān)閉后，連接通道內(nèi)氣體的壓力仍高于主噴管內(nèi)噴孔出口附近的壓力，導致仍有氣體注入主噴管，形成二次射流，產(chǎn)生側(cè)向力，隨壓力的逐漸平衡，殘存的側(cè)向力不斷減小，直至恢復到初始狀態(tài)。

圖7給出了閥動過程中閥門作動速度及側(cè)向力變化率隨時間的變化 曲線。由圖7可看出，側(cè)向力變化率并不像閥門作動速度隨時間的變化那樣是光滑連續(xù)的，它存在著強烈的非線性和不穩(wěn)定性。分析認為，這可能是由于閥動過程的每一時刻擴張段的激波是不穩(wěn)定的，即在穩(wěn)定的過程中，存在輕微的抖動現(xiàn)象。因此，導致了在閥動過程中出現(xiàn)側(cè)向力變化率的這種不穩(wěn)定性。

項目數(shù)值/kN差異率/%準穩(wěn)態(tài)最大側(cè)向力19．239標準動態(tài)最大側(cè)向力19．2780．202動態(tài)18mm開度時側(cè)向力19．161-0．405

2.2 側(cè)向力變化對比分析

2.2.1 作動方式對側(cè)向力的影響分析

分析閥動過程中作動方式對側(cè)向力的影響規(guī)律。圖8是不同頻率下閥動過程中閥門開度和側(cè)向力的對比情況。對比發(fā)現(xiàn)，側(cè)向控制力在2種作動方式下的表現(xiàn)與閥門開度的變化相似。即在閥門開啟后的一段時間和關(guān)閉前的一段時間之內(nèi)，三角波作動方式下的側(cè)向力要大于正弦波下的側(cè)向力；其余時間則相反。但與閥門開度表現(xiàn)不同的是開度相同的時刻，側(cè)向控制力并不相同，存在一定的響應差異(時間滯后)。

表2中，列出了這種有關(guān)時間滯后的部分數(shù)據(jù)。閥動中最大側(cè)向力的滯后時間約為閥門作動周期的8%，對應的工況為三角波10 Hz閥動狀態(tài)。

由表2中數(shù)據(jù)可知，在3種不同的頻率下，三角波的最大側(cè)向力滯后時間均大于正弦波下的值。同樣，閥門關(guān)閉后，側(cè)向力的相對恢復時間也表現(xiàn)出相同的規(guī)律。

從這2個角度分析表明，正弦波下側(cè)向力的響應性優(yōu)于三角波下的響應性。

2.2.2 作動頻率對側(cè)向力的影響

分析作動頻率對側(cè)向力的影響，圖9給出了閥動過程中側(cè)向力隨閥門開度的變化及準穩(wěn)態(tài)情況下側(cè)向力與開度的關(guān)系，即側(cè)向力隨空間的變化情況。表3和表4列出了幾個特殊位置的側(cè)向力大小。顯然，側(cè)向力隨空間或者閥門開度的變化規(guī)律表現(xiàn)為：側(cè)向力在閥門開啟和關(guān)閉過程中并不對稱，相同開度下，關(guān)閉過程中的側(cè)向力大于準穩(wěn)態(tài)下的值，二者又大于開啟過程中的值，且閥動過程中，側(cè)向力最大值出現(xiàn)的位置在閥門關(guān)閉過程中，稍滯后于閥門最大開度。這是由于閥動過程中二次噴射帶來的側(cè)向力變化滯后于閥門的作動。與上節(jié)所提及的側(cè)向力的時間滯后是一樣的，這是在閥門開度上的表現(xiàn)。

表2 閥動過程中最大側(cè)向力、最大側(cè)向力滯后時間及恢復時間

作動方式頻率f/Hz準穩(wěn)態(tài)側(cè)向力/kN動態(tài)側(cè)向力/kN側(cè)向力差異率/%三角波519．23918．926-1．630719．23918．826-2．1491019．23918．000-6．441正弦波519．23919．161-0．409719．23918．915-1．6831019．23918．803-2．265

在不同頻率下，相同開度時，側(cè)向力基本相同。只在三角波10 Hz工作狀態(tài)下有一些差別。排除這組數(shù)據(jù)，表3列出了不同工況下閥門作動到18 mm時的側(cè)向力大小，工況之間最大側(cè)向力的差別小于0.6%。表4列出了不同工況下閥門作動到9 mm時的側(cè)向力大小，工況之間的最大差別約2%。

由表2中關(guān)于時間滯后的數(shù)據(jù)，在最大側(cè)向力滯后時間τD方面，10 Hz下的τD大于7 Hz下的值，同時二者又大于5 Hz下的τD；相對恢復時間τR方面，7 Hz下的τR最大，5 Hz下的值次之，10 Hz下的τR最小。此外，由表2可見，最大側(cè)向力隨閥門動作頻率的變化似乎沒有明確的規(guī)律，有的頻率下，動態(tài)最大側(cè)向力大于最大閥門開度時的準穩(wěn)態(tài)值，而有的頻率下，則反之。分析認為，最大側(cè)向力的值取決于瞬態(tài)激波的位置，激波擺動頻率與閥門同步，但激波擺角可能大于準穩(wěn)態(tài)，也可能小于準穩(wěn)態(tài)。具體細節(jié)較復雜，需要進一步加密網(wǎng)格深入進行計算分析研究。

從表2可見，最大側(cè)向力滯后時間隨頻率變化很有規(guī)律，但閥門關(guān)閉后，恢復時間隨頻率變化沒有規(guī)律。分析認為，關(guān)閉后的響應時間與頻率沒有表現(xiàn)出明確關(guān)系，主要是由于氣動響應時間遠小于活塞作動周期，氣動響應頻率遠高于作動頻率。

表4 9 mm開度下側(cè)向力對比

3 結(jié)論

(1)在給定的閥門作動規(guī)律下，側(cè)向力的變化滯后于閥門的作動；閥門作動一個周期重新關(guān)閉后，由于噴管內(nèi)殘余激波的存在，短時間內(nèi)仍有側(cè)向力存在。

(2)在本文給定的2種閥門作動波形下，側(cè)向力隨時間的變化趨勢相同；就最大側(cè)向力出現(xiàn)時刻和閥門重新閉合后的流場恢復時間分析時間響應(時間滯后)，正弦波作動方式要優(yōu)于三角波。

(3)作動頻率對側(cè)向力響應影響的分析表明，在3種作動頻率下，在最大側(cè)向力滯后時間τD方面，10 Hz下的τD大于7 Hz下的值，同時二者又大于5 Hz下的τD。閥門關(guān)閉后，恢復時間隨頻率變化沒有明顯規(guī)律，分析認為，主要是由于氣動響應時間遠小于活塞作動周期，氣動響應頻率遠高于作動頻率。

本文所得到的結(jié)論是基于以下幾點假設得到的。首先，忽略發(fā)動機燃燒室部分，只截取噴管段幾何結(jié)構(gòu)進行分析；其次，在噴管進口給定了進口質(zhì)量流量，不考慮閥門開度的變化可能帶來的進口流量的改變。為了更深入地研究該問題，后期的研究將考慮這些因素對問題的影響。此外，本文只對正弦波、三角波2種波形及5、7、10 Hz 3種作動頻率進行了計算分析，后期將對其他作動波形和頻率展開研究，以期得到更普遍的規(guī)律。

[1] Kowal H J.Advances in thrust vectoring and the application of flow-control technology[R].Canadian Aeronautics and Space Joural,2002 CASI.

[2] Karen A Deere.Summary of fluidic vectoring research con-ducted at NASA Langley research center[C]//The 21st AIAA Applied Aerodynamics Conference,2003,Orlando,Florida.AIAA 2003-3800.

[3] Lee Jin-gyu,Chang Keun-shik.A 2D CD nozzle with thrust vector control under dilute gas-particle interaction[J].Computational Fluid DynamicsJournal,2003,11(4):420-431.

[4] Newton J F,Spaid F W.Interaction of secondary injectants and rocket exhaust for thrust vector control[J].AIAA 1962,32(8):1203-1211.

[5] 喬渭陽,蔡元虎,齊少軍,等.次流噴射控制推力矢量噴管的流場數(shù)值模擬[J].推進技術(shù),2000,21(16):15-20.

[6] 劉剛,楊永,李聰.二股流噴射控制推力矢量數(shù)值計算研究[J].航空計算技術(shù),2005,35(2):50-53.

[7] 鄧遠擷,鐘梓鵬,宋文艷.收斂-擴張噴管中運用次流推力矢量控制技術(shù)的計算研究[J].固體火箭技術(shù),2004,25(1):29-32.

[8] 賴川.用于推力矢量控制的氣體噴射混合流場研究[D].哈爾濱:哈爾濱工程大學,2000.

[9] 賴川,郜冶,陳步學.噴管內(nèi)二次噴射的激波測量與分析[J].推進技術(shù),2000,21(3):26-29.

[10] 劉輝,邢玉明.氣體二次噴射矢量噴管三維流場計算[J].北京航空航天大學學報,2009,35(10):1174-1178.

[11] 曹熙煒,劉宇,任軍學,等.工質(zhì)溫度對二次噴射推力矢量性能的影響[J].宇航學報,2010,31(12):2784-2787.

[12] 吳雄.固體發(fā)動機燃氣二次噴射理論與試驗研究[D].長沙：國防科技大學,2007.

[13] 吳雄,王中偉,焦紹球,等.氣體二次噴射推力矢量控制的響應面優(yōu)化設計[J].航空動力學報,2007,22(9):1569-1572.

[14] 吳雄,王中偉,焦紹球,等.燃氣二次噴射推力矢量控制系統(tǒng)啟動過程非定常效應研究[J].空氣動力學學報,2009,27(2):220-226.

[15] 王革,肖雪峰,李垚.雙孔燃氣二次噴射推力矢量性能影響因素數(shù)值分析[J].固體火箭技術(shù),2014,37(5):622-627.

[16] 肖雪峰.潛入式噴管燃氣二次噴射推力矢量控制流場計算[D].哈爾濱：哈爾濱工程大學,2013.

[17] 吳德銘,郜冶.實用計算流體力學基礎[M].哈爾濱：哈爾濱工程大學出版社,2008.

(編輯：崔賢彬)

Effect of valve actuation mode and frequency on lateral force in gas secondary injection system

LI Dong-dong,WANG Ge

(College of Aerospace and Civil Engineering,Harbin Engineering University，Harbin 150001，China )

Combining AUSM+format with Realizablek-εturbulence model,change law of lateral force in valve actuation process was explored and the effect of valve actuation mode (triangular wave,sinusoidal wave) and frequency(5,7 and 10 Hz) on lateral force was analyzed based on the simulation of flow field of nozzle in gas secondary injection system.Results show that change law of lateral force is similar with that of valve opening, in general,lateral force increases with the increasing of valve opening,and there is delay on time and space compared with valve actuation.Variations of lateral force lag behind that of the valve position.The largest lateral force occur after half period,while valve opening is less than the maximum value.There are still remnants of lateral forces in a short time after the valve is closed again.In the same valve opening,the lateral force in the closing process is greater than that of quasi-steady state,lateral force in opening process is the least among of them.Delay timeτDof triangular wave is greater than that of sinusoidal wave in the three frequencies.τDof 10 Hz is the greatest andτDof 5 Hz is the least among them in the two valve actuation mode based on the delay time of largest lateral force.Restoration timeτRof triangular wave is greater than that of sinusoidal wave in the three frequencies.τRof 7 Hz is the greatest andτRof 5 Hz is the least among them in the two valve actuation mode based on the restoration time of largest lateral force after the valve is closed again.

gas secondary injection;thrust vector control;numerical simulation;lateral force

2016-01-23；

2016-03-02。

李冬冬(1991—)，男，博士生，從事固體火箭發(fā)動機內(nèi)流場研究。E-mail：lidongdongheu@hrbeu.edu.cn

王革(1966—)，男，教授，從事火箭發(fā)動機內(nèi)彈道研究。E-mail：wangge@hrbeu.edu.cn

V435

A

1006-2793(2017)02-0151-07

10.7673/j.issn.1006-2793.2017.02.004