李 易

(西北工業(yè)大學(xué)陜西省空天飛行器設(shè)計重點(diǎn)實(shí)驗(yàn)室，西安 710072)

0 引言

級間分離是運(yùn)載火箭發(fā)射任務(wù)中的關(guān)鍵環(huán)節(jié)。對于串聯(lián)式多級火箭，常用的級間分離方式可分為熱分離和冷分離兩種。熱分離是在下一級火箭燃料耗盡時利用上一級火箭發(fā)動機(jī)點(diǎn)火后產(chǎn)生的噴流沖擊力完成級間分離。采用熱分離技術(shù)的火箭可以在一級火箭燃盡后迅速獲得二級火箭發(fā)動機(jī)產(chǎn)生的推力以約束火箭飛行姿態(tài)，而無須安裝專用的輔助發(fā)動機(jī)，這有利于降低火箭整體復(fù)雜度。本文主要針對固體運(yùn)載火箭的大氣層內(nèi)級間熱分離模擬問題開展研究。

由于熱分離過程中發(fā)動機(jī)燃?xì)獾拇嬖诤头蛛x時較高的飛行馬赫數(shù)，火箭面臨極其嚴(yán)苛的氣動力、熱環(huán)境。二子級火箭發(fā)動機(jī)的燃?xì)鈺_擊一子級火箭頂部，在排出級間段后與高速來流交互作用，形成復(fù)雜的化學(xué)非平衡高溫流場。在設(shè)計火箭的級間熱分離機(jī)構(gòu)和分離程序時，分離動力學(xué)特性是重要輸入條件，因此詳細(xì)的分離動力學(xué)研究是十分必要的。然而，火箭級間熱分離的風(fēng)洞實(shí)驗(yàn)昂貴而且實(shí)施難度高，目前在實(shí)驗(yàn)?zāi)Ｐ统叨群涂芍貜?fù)性方面仍有較大困難，數(shù)值仿真便是目前研究火箭級間分離動力學(xué)的主要途徑。傳統(tǒng)的火箭級間分離仿真方法有定常彈道仿真和計算流體力學(xué)(CFD)-飛行力學(xué)耦合仿真方法。定常彈道仿真沿用傳統(tǒng)的彈道仿真方法解算剛體運(yùn)動學(xué)和動力學(xué)方程，但需要大量的氣動數(shù)據(jù)，數(shù)據(jù)可通過大量的風(fēng)洞實(shí)驗(yàn)和定常CFD計算獲得。這種方法在工業(yè)界較為流行，可在型號設(shè)計階段快速進(jìn)行參數(shù)化研究，但氣動數(shù)據(jù)生成周期長、成本高，而且無法記錄飛行器分離過程中的非定常現(xiàn)象。CFD-飛行力學(xué)耦合仿真方法是在20世紀(jì)80年代重疊網(wǎng)格技術(shù)(Chimera方法)出現(xiàn)后開始得到應(yīng)用，這種仿真可較多地計入分離過程中的非定常流動現(xiàn)象。CFD-飛行力學(xué)耦合仿真方法早期主要用于機(jī)載武器的投放模擬，近20年來隨著動網(wǎng)格技術(shù)和計算機(jī)技術(shù)的發(fā)展，這種耦合仿真方法在飛行器的各種多體分離研究中得到了廣泛應(yīng)用。

關(guān)于火箭級間熱分離，發(fā)動機(jī)燃?xì)馀c空氣混合后持續(xù)燃燒，內(nèi)外流場實(shí)質(zhì)上為化學(xué)非平衡流動。過去有關(guān)級間熱分離運(yùn)動的研究，主要針對分離初始時刻或分離過程中某一時刻的流場進(jìn)行定常模擬，較少耦合飛行力學(xué)方程進(jìn)行瞬態(tài)流場的研究，而且在定常流場的研究中主要采用單組分方法處理流場，即將火箭燃?xì)饣旌衔锱c空氣看作同一種氣體，且假設(shè)為完全氣體(不考慮熱力學(xué)效應(yīng))，較少考慮內(nèi)外流場中各氣體組分之間的非平衡效應(yīng)問題。Applebaum等曾對火箭熱分離時刻流場模擬中的羽流建模問題進(jìn)行研究，筆者曾采用非平衡流法計算了火箭級間段的流場。由于計算量龐大，較少有學(xué)者采用 CFD-飛行力學(xué)耦合方法進(jìn)行級間熱分離仿真，已發(fā)表的成果主要采用單組分完全氣體法，且較少計入分離時刻發(fā)動機(jī)的推力瞬態(tài)變化過程。筆者在過去工作中，分別采用單組分完全氣體法和雙組分完全氣體法進(jìn)行了CFD-飛行力學(xué)耦合仿真研究，并通過參考發(fā)動機(jī)試車數(shù)據(jù)實(shí)時更新燃燒室邊界條件,實(shí)現(xiàn)對火箭瞬態(tài)推力的模擬。

在運(yùn)載火箭的研制工作中，科研人員受計算條件和CFD軟件技術(shù)的限制，往往無法在CFD-飛行力學(xué)耦合仿真過程中同時考慮熱分離期間發(fā)動機(jī)的內(nèi)外非平衡流動問題。為便于大規(guī)模的參數(shù)化研究，科研人員往往需使用低熱力學(xué)等級的流動模型簡化對流場的模擬，然而對簡化后的流場仿真的保真度缺乏研究，仿真結(jié)果仍然蘊(yùn)藏著不確定性，不利于降低對昂貴的風(fēng)洞實(shí)驗(yàn)和飛行試驗(yàn)的需求。

本文使用超級計算機(jī)和化學(xué)非平衡流CFD程序，對現(xiàn)行CFD方法中可用于模擬火箭羽流的不同熱力學(xué)等級的流動模型進(jìn)行比較分析，包括多相多組分化學(xué)非平衡流、多組分化學(xué)非平衡氣體、多組分完全氣體、雙組分完全氣體和單組分完全氣體等，對各種模型的優(yōu)缺點(diǎn)和耦合仿真中適用性進(jìn)行定性總結(jié)，最后結(jié)合CFD-飛行力學(xué)耦合仿真算例進(jìn)行定量分析，供科研人員今后在類似的工程中進(jìn)行參考。

1 級間熱分離流場結(jié)構(gòu)分析

圖1為串聯(lián)式固體火箭級間熱分離時刻的流場CFD模擬結(jié)果(飛行馬赫數(shù)為2.6，高度為15 km)?；鸺七M(jìn)劑在發(fā)動機(jī)燃燒室內(nèi)通過燃燒將化學(xué)能轉(zhuǎn)換為氣體熱能，再通過拉瓦爾噴管將燃?xì)獾臒崮苻D(zhuǎn)化為噴流的動能。發(fā)動機(jī)燃?xì)馀懦鰢姽軙r為超聲速噴流，在與一子級頂端壁面相遇后形成復(fù)雜的沖擊流場。在沖擊流場中同時存在亞聲速與超聲速區(qū)域，存在復(fù)雜激波-膨脹波系，以及激波相互干涉和激波與邊界層之間的交互作用。在級間段，羽流沖擊形成的局部高壓會導(dǎo)致沿二子級箭體向上游方向的逆壓梯度的產(chǎn)生，進(jìn)而誘導(dǎo)二子級尾部壁面發(fā)生流動分離(如圖1所示)，這種分離在文獻(xiàn)中常被稱為羽流誘導(dǎo)流動分離(Plume-induced flow separation)。這種流動分離在箭體周向上往往是非對稱的，并且由于上一級火箭尾部距離箭體重心較遠(yuǎn)，會使得火箭承受較大的側(cè)向力矩。當(dāng)一子級火箭在分離力的驅(qū)動下遠(yuǎn)離二子級時，這種羽流誘導(dǎo)流動分離則會突然消失。

圖1 火箭級間段流場[4]

2 CFD-飛行力學(xué)耦合仿真方法

在過去20余年中，CFD-飛行力學(xué)耦合仿真逐漸成為一種模擬非定常流場中多物體相對運(yùn)動的有效途徑，例如模擬飛機(jī)投放、人員彈射和飛行器的動態(tài)特性研究。根據(jù)CFD方程與飛行力學(xué)方程的耦合程度，CFD-飛行力學(xué)耦合仿真可分為準(zhǔn)定常仿真和非定常仿真兩種。在準(zhǔn)定常耦合仿真中，CFD方程與飛行力學(xué)方程在每個仿真時間步內(nèi)各自獨(dú)立求解，整個仿真過程由離散在各個時間步上的定常流場計算組成；在非定常耦合仿真中，CFD方程與飛行力學(xué)方程聯(lián)立求解，物體運(yùn)動狀態(tài)發(fā)生變化時其表面局部對當(dāng)?shù)亓鲌龅臄_動可被計入，因此非定常耦合仿真可模擬飛行器的氣動阻尼現(xiàn)象。

CFD-飛行力學(xué)耦合仿真方法也被用于研究火箭的級間分離運(yùn)動，然而其應(yīng)用范圍主要是串聯(lián)式或并聯(lián)式火箭的冷分離運(yùn)動。在前期研究基礎(chǔ)上，筆者成功實(shí)現(xiàn)了火箭級間熱分離三維非定常仿真。運(yùn)用CFD-飛行力學(xué)仿真方法研究火箭熱分離動力學(xué)的關(guān)鍵是要在計算中準(zhǔn)確預(yù)測下面級火箭所受的羽流沖擊和上面級火箭尾部的流動分離，而解決這一關(guān)鍵問題的核心任務(wù)是尋找合適的發(fā)動機(jī)內(nèi)外流氣體模型。

3 計算模型

3.1 氣體模型

以鋁化復(fù)合燃料(由過氯酸銨、端羥基聚丁二烯和鋁組成)固體火箭發(fā)動機(jī)為例，在燃燒室內(nèi)其燃燒產(chǎn)物主要是氣體和液態(tài)鋁氧化物。在噴管膨脹段下游，隨著溫度降低液態(tài)鋁氧化物逐漸轉(zhuǎn)化為固態(tài)氧化鋁粒子。由NASA CEA代碼計算可知，其噴管流場中有160余種可能的氣態(tài)、液態(tài)和固態(tài)成份。高溫流場中的局部流體組分與當(dāng)?shù)販囟取毫蜁r間相關(guān)，隨著燃?xì)庀驀姽芟掠瘟鲃樱細(xì)鉁囟冉档?，其中的氣體組分會發(fā)生變化。在噴管外部，發(fā)動機(jī)燃?xì)庵械母蝗冀M分還會持續(xù)與空氣中的組分發(fā)生化學(xué)反應(yīng)。

火箭發(fā)動機(jī)燃?xì)饬魇且环N非定常多相化學(xué)非平衡流動，目前CFD技術(shù)在精確模擬發(fā)動機(jī)羽流方面無論在理論上、數(shù)值算法上還是計算技術(shù)上仍有較大困難。表1列舉了幾種不同等級的流動模型。對比計算效率時，計算網(wǎng)格點(diǎn)數(shù)量為680萬，使用192個CPU 核(2.4GHZ)。其中多組分化學(xué)非平衡流計算在技術(shù)上可以實(shí)現(xiàn)，但由于計算量過大，工業(yè)應(yīng)用中不適用于大規(guī)模3D流場計算。雙組分量熱完全氣體模型是將空氣和發(fā)動機(jī)羽流各認(rèn)為是一種量熱完全氣體(calorically perfect gas)，并且不考慮組分之間的化學(xué)反應(yīng)。單組分量熱完全氣體模型是在級間分離仿真方面常用的一種氣體模型，即將空氣和發(fā)動機(jī)羽流看作同一種氣體。該模型適用于普通計算單氣體組分的CFD代碼，計算量最小。

本節(jié)著重探討工程上可行的3種氣體模型，即多組分化學(xué)非平衡流模型，雙組分量熱完全氣體模型和單組分量熱完全氣體模型。

3.1.1 多組分化學(xué)非平衡流模型

這種氣體模型是將流動描述為一種由多種熱力學(xué)完全氣體(忽略分子的熱力學(xué)非平衡效應(yīng))組成的化學(xué)非平衡混合體(各氣體組分之間存在化學(xué)反應(yīng))。流動中各組分之間的化學(xué)反應(yīng)可由有限速率模型描述

(1)

式中，

X

表示燃燒物或反應(yīng)物，

v

和

w

分別表示反應(yīng)物和燃燒物的摩爾數(shù)，

k

和

k

分別表示前向和后向反應(yīng)速率常數(shù)。

k

=

c

T

e

-/

(2)

(3)

式中，

C

,

α

和

ε

均來自于實(shí)驗(yàn)，

K

表示平衡常數(shù)。表2為本文計算固體發(fā)動機(jī)內(nèi)外流場所使用的化學(xué)動力學(xué)模型，其中第1～8項代表H/O反應(yīng)系統(tǒng)，第9～11項代表CO/CO反應(yīng)系統(tǒng)。在火箭燃?xì)庵?，除了燃料殘余，與空氣接觸發(fā)生后續(xù)燃燒的主要是中間燃燒產(chǎn)物，例如OH和CO。表2中第12～17項代表HCl的抑制作用。

表1 流動模型比較

表2 一種固體發(fā)動機(jī)化學(xué)動力學(xué)模型

由于該氣體模型引入了多種氣體成份，CFD技術(shù)在處理每一種氣體組分時需在描述流動的Navier-Stokes(N-S)方程中增加一項輸運(yùn)方程，而且需要處理額外的化學(xué)反應(yīng)相關(guān)項，因此運(yùn)用該氣體模型會導(dǎo)致流場計算量過大，不適合工程上的大規(guī)模參數(shù)化分析。

3.1.2 雙組分量熱完全氣體模型

該模型是將發(fā)動機(jī)燃?xì)饪醋饕环N氣體組分，將大氣看作另一種氣體組分。兩種氣體之間無化學(xué)反應(yīng)，且均為量熱完全氣體，即氣體的比熱比(

γ

)保持不變。該模型僅在標(biāo)準(zhǔn)N-S方程中多引入一項輸運(yùn)方程，且無化學(xué)反應(yīng)相關(guān)項，因此計算量較低。利用大型集群計算機(jī)，可在工程任務(wù)中使用該模型。關(guān)于雙組分量熱完全氣體模型的闡述可參考文獻(xiàn)[1]。

3.1.3 單組分量熱完全氣體模型

該氣體模型將外部空氣和發(fā)動機(jī)羽流看作同一種氣體，流場計算中僅求解標(biāo)準(zhǔn)N-S方程，計算量在3種模型中最低，任務(wù)可在普通集群計算機(jī)上運(yùn)行。

關(guān)于化學(xué)非平衡流動和量熱完全氣體概念的詳細(xì)闡述可參考文獻(xiàn)[38]，固體發(fā)動機(jī)的化學(xué)動力學(xué)模型可參考文獻(xiàn)[8]。

3.2 氣體模型比較

在串聯(lián)式運(yùn)載火箭的級間分離仿真中，以模擬一、二級分離為例，須重點(diǎn)關(guān)注的是級間段的流動分離(見圖1和圖2)和作用在一級火箭上的羽流沖擊力。3種氣體模型在計算羽流沖擊力方面沒有顯著差異，但在計算流動分離時差別則較明顯。圖3為用3種氣體模型模擬級間段流動分離的計算結(jié)果(二維軸對稱計算)，其中在運(yùn)用多組分化學(xué)非平衡流模型時，引入了12種氣體組分和17種化學(xué)反應(yīng)模型。圖 3中的計算結(jié)果顯示，單組分氣體模型相對高熱化學(xué)等級的多組分氣體模型會低估二級尾部壁面上的邊界層分離長度，這種差異會影響對二級火箭級間分離過程中所受氣動力矩的預(yù)測；雙組分氣體模型則與多組分氣體模型有相似的計算結(jié)果，而雙組分氣體模型的計算量卻小得多。

圖2 二級火箭尾部非對稱流動分離(M)[29]

圖3 級間段流動分離(分別使用3種氣體模型計算)

3.3 黏性效應(yīng)

流場計算時若忽略氣體黏性效應(yīng)，即僅求解Euler方程，可以降低計算量。在計算高速飛行器的氣動力時常采用無黏(Euler)方法，在流場中無強(qiáng)烈的流動分離時使用Euler法可以得出與N-S方法近似的結(jié)果。過去在級間分離仿真中也常采用Euler法，但由于忽略黏性效應(yīng)，仿真中無法預(yù)測級間段的流動分離，從而無法詳細(xì)預(yù)測二級火箭所受氣動力矩的演化過程。

4 仿真比較

基于筆者前期研究計算結(jié)果，對上文所述計算模型做進(jìn)一步比較分析。由于多組分氣體模型不易用于大型3D計算網(wǎng)格，在仿真中僅采用雙組分氣體模型和單組分氣體模型，并結(jié)合Euler(無黏)和N-S(有黏)方法進(jìn)行綜合比較分析，即進(jìn)行了雙組分氣體有黏計算(2-species-gas, N-S)、單組分氣體有黏計算(Single-species-gas, N-S)、雙組分氣體無黏計算(2-species-gas, Euler)和單組分氣體無黏計算(Single-species-gas, Euler)。仿真初始條件如表3所示。

表3 仿真初始條件

當(dāng)火箭在分離時刻有飛行攻角時，二級火箭尾部將存在非對稱的流動分離(圖1)，致使二級火箭在兩級解鎖后承受一個初始的俯仰力矩，火箭在兩級解鎖后仿真結(jié)果顯示各種計算模型對火箭線位移的預(yù)測影響較小，如圖4和圖5所示，圖中

x

表示火箭兩級在分離坐標(biāo)系(描述相對運(yùn)動)中的軸向位移。分離坐標(biāo)系是慣性坐標(biāo)系，其原點(diǎn)位于火箭縱向軸線與分離面的交點(diǎn)，

x

沿火箭縱向軸線指向火箭頭部方向?yàn)檎?；分離坐標(biāo)系相對地面坐標(biāo)系(忽略地球自轉(zhuǎn))以火箭兩級分離初始時刻的速度保持平移。可以看出，在二級火箭發(fā)動機(jī)提供的分離力作用下，火箭兩級在分離坐標(biāo)系中的運(yùn)動方向相反；一子級火箭由于燃料耗盡，主要剩余結(jié)構(gòu)質(zhì)量，因此軸向速度的變化要明顯大于二級火箭；在地面坐標(biāo)系中，一子級火箭的飛行高度仍持續(xù)上升，但在二級火箭燃?xì)獾臎_擊作用下其飛行速度會迅速下降，而二級火箭在發(fā)動機(jī)的推力作用下速度則持續(xù)增大。因此，分離后一子級火箭不會重新“追上”二級火箭造成碰撞。

圖4 二級火箭的軸向線位移

圖5 一子級火箭的軸向線位移

由于二級火箭的軸向運(yùn)動主要取決于發(fā)動機(jī)推力，一子級火箭的軸向運(yùn)動主要取決于二級火箭發(fā)動機(jī)燃?xì)獾臎_擊，所以各種計算模型對火箭兩級分離后軸向運(yùn)動的預(yù)測影響較小。不同計算模型主要影響火箭兩級分離后的角位移，尤其是二級火箭的角位移。如圖 6和圖 7所示，無黏計算由于無法預(yù)測級間段的流動分離(圖 8)，因此明顯低估了火箭的俯仰運(yùn)動幅度；而采用單組分氣體的有黏計算，雖可模擬出火箭級間段的流動分離，但因會低估二級火箭尾部邊界層分離的長度(圖 8)，所以分離初始階段預(yù)測出的火箭俯仰力矩偏小。火箭分離過程中，隨著兩級分離距離的增加，二級火箭尾部的邊界層分離會突然消失，從而給火箭控制系統(tǒng)造成擾動，對分離初始階段流動分離峰值的準(zhǔn)確預(yù)測便有助于火箭的控制系統(tǒng)設(shè)計，改善火箭的飛行品質(zhì)。由圖 6和圖 7還可以看出，分離后二級火箭的俯仰角速度呈上升趨勢。由于運(yùn)載火箭在超聲速飛行時，其飛行攻角往往需被限制在較小的范圍內(nèi)(不超過3°)，因此二級火箭姿控系統(tǒng)在兩級分離后(0.02～0.04 s)應(yīng)及時啟動以抑制其俯仰運(yùn)動。

圖6 二級火箭的俯仰角速度

圖7 二級火箭在級間分離過程中所受的氣動力矩

圖8 火箭分離初始時刻流場

5 結(jié)論

運(yùn)用CFD-飛行力學(xué)耦合仿真方法研究火箭熱分離動力學(xué)的關(guān)鍵是尋找合適的發(fā)動機(jī)燃?xì)饽Ｐ汀Ｓ捎谟嬎阈实?，非平衡流模型還不適用于工程研究，尤其是火箭設(shè)計階段的大規(guī)模參數(shù)化仿真。雙組分氣體模型在精度上接近非平衡流模型，且計算量較小，未來在工程研究中具有較強(qiáng)的應(yīng)用潛力。單組分氣體模型會低估火箭級間段的流動分離，從而導(dǎo)致對二級火箭的俯仰或偏航運(yùn)動的預(yù)測值偏低。無黏計算則無法預(yù)測火箭級間段的流動分離。然而，由于火箭級間段的流動分離僅出現(xiàn)在分離初始階段，在工程研究中，若主要關(guān)注火箭的大范圍線位移運(yùn)動，也可采用單組分氣體模型和無黏計算，從而進(jìn)一步降低計算量。