陳士強(qiáng)，范瑞祥，張曙輝，黃 兵，黃 輝

（1. 北京宇航系統(tǒng)工程研究所，北京，100076；2. 中國(guó)運(yùn)載火箭技術(shù)研究院，北京，100076）

引射循環(huán)預(yù)冷回流低溫貯箱流動(dòng)與傳熱特性研究

陳士強(qiáng)1，范瑞祥2，張曙輝1，黃 兵1，黃 輝1

（1. 北京宇航系統(tǒng)工程研究所，北京，100076；2. 中國(guó)運(yùn)載火箭技術(shù)研究院，北京，100076）

運(yùn)載火箭貯箱內(nèi)部流動(dòng)與傳熱特性與火箭息息相關(guān)，中國(guó)新一代運(yùn)載火箭氧系統(tǒng)普遍采用引射循環(huán)預(yù)冷技術(shù)。介紹了低溫動(dòng)力系統(tǒng)引射循環(huán)預(yù)冷過程中低溫貯箱內(nèi)涉及到主要物理過程，采用基于自由界面追蹤的VOF模型，研究了引射循環(huán)預(yù)冷回流低溫貯箱內(nèi)氣液兩相流動(dòng)與傳熱特性。結(jié)果表明：地面狀態(tài)循環(huán)回流氣液兩相流對(duì)貯箱內(nèi)相分布、溫度、壓力和速度具有顯著影響，并進(jìn)一步討論了引射循環(huán)預(yù)冷工程優(yōu)化原則，為相關(guān)工程應(yīng)用提供指導(dǎo)。

引射循環(huán)預(yù)冷回流；低溫貯箱；流動(dòng)；傳熱

0 引 言

貯箱是運(yùn)載火箭重要結(jié)構(gòu)件和推進(jìn)劑貯存裝置，其內(nèi)部流動(dòng)與傳熱特性與運(yùn)載火箭息息相關(guān)。中國(guó)新一代運(yùn)載火箭氧系統(tǒng)普遍采用了引射循環(huán)預(yù)冷技術(shù)。在低溫動(dòng)力系統(tǒng)引射循環(huán)預(yù)冷過程中，貯箱通過輸送管向發(fā)動(dòng)機(jī)系統(tǒng)輸入低溫推進(jìn)劑，經(jīng)過高溫管路、泵等組件后，低溫推進(jìn)劑（包含蒸發(fā)形成的氣態(tài)推進(jìn)劑和引射氣體）回流至貯箱內(nèi)，并與原有推進(jìn)劑相互作用，影響貯箱內(nèi)流動(dòng)與傳熱。在引射循環(huán)預(yù)冷進(jìn)行過程中，貯箱內(nèi)涉及到低溫推進(jìn)劑復(fù)雜氣液兩相流動(dòng)與傳熱過程，主要包括低溫推進(jìn)劑溫度分層、貯箱內(nèi)引射氣體運(yùn)動(dòng)、回流氣液混合物與貯箱內(nèi)推進(jìn)劑傳熱傳質(zhì)作用、低溫推進(jìn)劑與氣枕接觸面的相變與換熱、氣枕壓力變化與溫度分層等。

Han[1]、Panzarella[2]等對(duì)氦引射循環(huán)預(yù)冷過程低溫貯箱內(nèi)的氣液兩相流動(dòng)和傳熱特性進(jìn)行了分析，為后續(xù)工作奠定了基礎(chǔ)；Cho等[3,4]采用基于一維、平衡傳質(zhì)非平衡傳熱模型研究了引射氦氣進(jìn)入液氧貯箱后系統(tǒng)的溫度特性，為引射循環(huán)預(yù)冷工程應(yīng)用提供了有益的指導(dǎo)，然而該模型在流動(dòng)與傳熱細(xì)節(jié)分析方面仍存在不足；文獻(xiàn)[5]介紹了在低溫動(dòng)力系統(tǒng)引射循環(huán)預(yù)冷全系統(tǒng)特性研究方面開展的一些工作。本文構(gòu)建了一種二維引射循環(huán)預(yù)冷回流低溫貯箱內(nèi)氣液兩相流動(dòng)與傳熱特性模型，分析了地面狀態(tài)循環(huán)回流氣液兩相流對(duì)貯箱內(nèi)相分布、溫度、壓力和速度的影響，力求從更多細(xì)節(jié)加強(qiáng)對(duì)引射循環(huán)預(yù)冷的規(guī)律性認(rèn)識(shí)，為相關(guān)工程應(yīng)用提供參考。

1 模型構(gòu)建

1.1 幾何模型與網(wǎng)格劃分

貯箱計(jì)算模型采用二維模型[6]，軸向長(zhǎng)度設(shè)定為12.8 m，筒段內(nèi)徑為Φ2.3 m，上、下橢球底高0.7 m。推進(jìn)劑輸送管位于下底中心位置，內(nèi)徑為Φ0.22 m；預(yù)冷回流管內(nèi)徑為Φ0.042 m，與貯箱接口位于下底弧段中心位置并平行于法線方向。圖1為貯箱二維幾何模型及網(wǎng)格劃分。

1.2 計(jì)算模型與初始邊值條件

引射氣體為氦氣，貯箱內(nèi)低溫推進(jìn)劑為液氧，計(jì)算模型為基于自由界面追蹤的VOF方法，包含液氧與氧氣之間的相變傳質(zhì)過程，采用FLUENT 6.3.26計(jì)算平臺(tái)。

1.2.1 基本假定

a）工質(zhì)為牛頓流體，氣態(tài)推進(jìn)劑可壓縮；

b）氣液之間存在熱相變，忽略不可冷凝氣體（如氦氣）在低溫推進(jìn)劑中的溶解；

c）貯箱內(nèi)氣液兩相流動(dòng)過程中存在較為明顯的連續(xù)相界面。

1.2.2 控制方程

a）基本控制方程。

整體質(zhì)量守恒方程：

式中 ρ為密度；u為速度矢量；Sm為源項(xiàng)，包括從分散的二級(jí)相中加入到連續(xù)相的質(zhì)量（如液滴的蒸發(fā)）和自定義源項(xiàng)。

j相質(zhì)量守恒方程為

j=1ijm˙ji分別為i相向j相和j相向i相的質(zhì)量傳遞。

動(dòng)量守恒方程在慣性（非加速）坐標(biāo)系中i方向上的表達(dá)式為

式中 μ為粘性系數(shù)；I為單位張量。

能量守恒方程為

式中 e為比內(nèi)能；keff為有效熱傳導(dǎo)系數(shù)；T為溫度；Sh為源項(xiàng)。

b）附加方程。

附加方程包括狀態(tài)方程和比內(nèi)能方程：

式中 R為理想氣體常數(shù)；H為比焓。

c）表面張力方程（微重力工況需特別考慮）。

表面張力選擇Brackbill提出的連續(xù)表面張力（Continuum Surface Force，CSF）模型，作為動(dòng)量方程的一個(gè)源項(xiàng)。

d）不可冷凝氣體方程。

對(duì)于存在3種流體的系統(tǒng)，如液氧、氧氣和氦氣，需添加混合氣中不可冷凝氣體體積含量的控制方程：

式中ncρ為不可冷凝氣體密度；Dnc為不可冷凝氣體擴(kuò)散系數(shù)。同時(shí)，不可冷凝氣體對(duì)蒸汽的分壓、傳熱和混合氣體的物性存在一定影響。

e）界面追蹤方程。

在整個(gè)流場(chǎng)中定義函數(shù)F，在每個(gè)劃分的網(wǎng)格中，F(xiàn)為流體體積與網(wǎng)格體積之比，滿足如下方程：

式中 當(dāng)F=1時(shí)，網(wǎng)格中充滿某種流體；F=0時(shí)，網(wǎng)格內(nèi)不含流體；在流體界面處0＜F＜1。模擬中需要求出每個(gè)網(wǎng)格中的F值，然后構(gòu)造出運(yùn)動(dòng)界面形狀。

f）湍流控制方程。

湍流控制方程可由-kε標(biāo)準(zhǔn)模型計(jì)算得出：

1.2.3 初始邊值條件

貯箱內(nèi)液氧初始高度為11.0 m，氣枕組分為氧氣；推進(jìn)劑出流口為質(zhì)量流量邊界，預(yù)冷回流口為速度邊界，貯箱頂端氣枕排氣口為壓力邊界；壁面為自然對(duì)流邊界，換熱系數(shù)為1.0 W/（m2?K），是結(jié)合工程經(jīng)驗(yàn)及保溫層后的折算系數(shù)；氣枕初始?jí)毫?.12 MPa。

氦引射循環(huán)預(yù)冷回流貯箱模型典型參數(shù)設(shè)定為：推進(jìn)劑出流口邊界為1.0 kg/s的純液氧，溫度91.8 K；預(yù)冷回流口邊界速度為2.8 m/s，方向沿邊界法向，溫度96.5 K，液氧體積分?jǐn)?shù)25%，氦氣體積分?jǐn)?shù)20%。

2 分析與討論

2.1 典型特性分析

圖2～6分別給出了氦引射循環(huán)預(yù)冷回流貯箱內(nèi)不同時(shí)刻氧氣體積分?jǐn)?shù)、氦氣體積分?jǐn)?shù)、速度、溫度和壓力分布。

t=15 s時(shí)，氦引射循環(huán)預(yù)冷回流貯箱內(nèi)的流線和溫度場(chǎng)分布情況如圖7所示。

氦引射循環(huán)預(yù)冷回流貯箱內(nèi)主要流動(dòng)與傳熱具有如下特點(diǎn)：

a）貯箱推進(jìn)劑出流口預(yù)冷流速對(duì)貯箱內(nèi)傳熱和流動(dòng)過程影響非常微弱；氦引射循環(huán)預(yù)冷系統(tǒng)起動(dòng)后，氣液兩相流逐漸從回流口管路進(jìn)入貯箱，貯箱內(nèi)部各參數(shù)初始處于微弱擾動(dòng)狀態(tài)。

b）隨著預(yù)冷的進(jìn)行，高速回流氣液混合物進(jìn)入貯箱并與液氧進(jìn)行劇烈的動(dòng)量交互而降速；受到自身重力影響及上層液氧、貯箱壁面的阻礙作用，回流推進(jìn)劑帶動(dòng)周圍液氧一方面沿回流入射方向流動(dòng)，另一方面向箱底下沉形成逆時(shí)針環(huán)流；環(huán)流帶來貯箱底部液氧壓力的波動(dòng)并影響回流推進(jìn)劑的入射方向，使其逐漸靠近回流口附近的壁面形成沿壁面向上的回流氣液兩相射流；而入射方向的改變促進(jìn)了環(huán)流的擴(kuò)大，使環(huán)流涉及的范圍越來越廣，同時(shí)破碎的氣泡受到浮力作用開始上浮，進(jìn)一步促進(jìn)了環(huán)流的發(fā)展，最終形成環(huán)繞整個(gè)貯箱內(nèi)部貼近箱壁的大環(huán)流，且頂部速度高于底部，大環(huán)流內(nèi)部受到氣液兩相混合作用和湍流擾動(dòng)形成多個(gè)復(fù)雜的無規(guī)則小環(huán)流。

c）在組分特性方面，預(yù)冷回流的氣體速度較高，與液氧進(jìn)行劇烈摻混后破碎為大量小氣泡，分散在液氧中的大量氧氣無法迅速完成換熱和冷凝過程，仍以氣態(tài)形式存在，引起貯箱液面一定程度的升高和波動(dòng)，氦氣的進(jìn)入進(jìn)一步加劇了貯箱內(nèi)組分分布的混亂程度（對(duì)于不同的預(yù)冷回流口邊界條件，貯箱內(nèi)組分情況會(huì)有明顯差異）。

d）在溫度特性方面，高溫回流氣液兩相推進(jìn)劑會(huì)引起周圍液氧溫度的升高，壁面自然對(duì)流換熱和氣枕壓力的升高會(huì)進(jìn)一步惡化貯箱內(nèi)液氧溫度品質(zhì)，而液氧的蒸發(fā)和氦氣入射對(duì)蒸發(fā)吸熱的促進(jìn)作用[1,2]有利于體系溫度的降低，綜合表現(xiàn)為回流口附近的核心區(qū)和近壁面區(qū)保持較高的溫度，貯箱內(nèi)大部分區(qū)域溫度處于92.5 K左右（見圖7），大環(huán)流和小渦流的存在有利于貯箱內(nèi)推進(jìn)劑溫度的均勻分布，破壞了可能存在的局部過熱和嚴(yán)重的熱分層效應(yīng)。

e）15 s后箱內(nèi)流動(dòng)基本穩(wěn)定，氣體的不斷逸出引發(fā)液面的劇烈波動(dòng)，同時(shí)氣枕和箱內(nèi)推進(jìn)劑壓力也會(huì)受到影響而有所增加。

與自然循環(huán)預(yù)冷回流貯箱[6,7]相比，氦引射循環(huán)預(yù)冷回流貯箱具有不同特性，主要表現(xiàn)在如下幾個(gè)方面：

a）由于氣相組分中氦氣的存在，回流混合物密度較小，受浮力影響較大，貯箱內(nèi)原有液氧對(duì)回流混合物沿徑向流動(dòng)的阻礙作用更為明顯，使得回流混合物在箱內(nèi)未形成環(huán)流時(shí)便呈現(xiàn)近似沿軸向向上流動(dòng)，并且受環(huán)流影響更為顯著，回流開始2 s后已在回流口附近形成貼壁面流動(dòng)。

b）氦氣與氧氣以大氣團(tuán)和小氣泡的形式混合在一起，受到氦氣低密度影響，混合氣上浮速度較自然循環(huán)預(yù)冷中純氧氣更快，氣泡在10 s左右便可以到達(dá)氣枕，系統(tǒng)內(nèi)大環(huán)流也隨之基本形成。

c）氦氣泡沿壁面上浮后大部分逸出到氣枕中，小部分隨環(huán)流重新回流到貯箱中部，受浮力作用未能進(jìn)入貯箱底層。

2.2 工程優(yōu)化討論

在工程實(shí)際中，氦引射循環(huán)預(yù)冷開始前系統(tǒng)一般處于自然循環(huán)預(yù)冷狀態(tài)[5]（引射氦氣未注入），氦氣注入后由于回流組分中氦氣的存在，密度降低，強(qiáng)化了貯箱內(nèi)原有液氧對(duì)回流混合物沿徑向流動(dòng)的阻礙作用，易于近似沿軸向向上流動(dòng)，環(huán)境的自然對(duì)流換熱和箭體內(nèi)部的主動(dòng)熱吹除對(duì)貯箱近壁面處推進(jìn)劑加熱持續(xù)作用，氦引射循環(huán)預(yù)冷快速回流特性對(duì)箱內(nèi)推進(jìn)劑溫度品質(zhì)的提高具有明顯的促進(jìn)作用；同時(shí)氦氣受浮力和環(huán)流影響基本不會(huì)出現(xiàn)在貯箱底部出流口附近，并且對(duì)氧氣上浮起到牽引促進(jìn)作用，可以較好地規(guī)避可能出現(xiàn)的出流推進(jìn)劑夾氣風(fēng)險(xiǎn)。

氦氣快速通過貯箱內(nèi)液氧而逸出進(jìn)入氣枕，如果長(zhǎng)時(shí)間引射，受貯箱排氣能力限制，氣枕壓力會(huì)受到較大影響，對(duì)系統(tǒng)穩(wěn)定性不利，因此需要較大的初始?xì)庹砣莘e。

射前的大量過冷液氧補(bǔ)加可以進(jìn)一步提高貯箱推進(jìn)劑溫度品質(zhì)，降低體系含氣率?；亓骰旌衔餃囟容^高，受箱內(nèi)低溫液氧冷卻后溫度會(huì)有所降低，同時(shí)質(zhì)量流量很小，對(duì)貯箱內(nèi)推進(jìn)劑溫度升高影響不明顯。箱內(nèi)液氧升溫汽化后會(huì)吸收汽化潛熱，氦氣的存在會(huì)降低氧分壓，進(jìn)一步促進(jìn)液氧汽化，引起推進(jìn)劑溫度降低，利于體系溫度穩(wěn)定。

氦氣質(zhì)量流量是影響貯箱內(nèi)流動(dòng)與傳熱特性的重要參數(shù)：流量偏低會(huì)降低引射回流對(duì)貯箱內(nèi)推進(jìn)劑的摻混作用，不利于推進(jìn)劑溫度品質(zhì)提高；流量偏高則增大貯箱排氣能力需求，同時(shí)增加出流推進(jìn)劑加氣的風(fēng)險(xiǎn)；對(duì)于基礎(chǔ)級(jí)以上的動(dòng)力系統(tǒng)，過高的氦氣需求增加了結(jié)構(gòu)質(zhì)量和氦氣用量，不利于運(yùn)載能力的提高。

引射循環(huán)預(yù)冷回流氦氣泡的破裂和上浮、貯箱內(nèi)大環(huán)流和局部繞流引起貯箱液面較大晃動(dòng)，可能進(jìn)一步引發(fā)流固耦合振動(dòng)，不利于系統(tǒng)穩(wěn)定，工程研制中需特別考慮。

3 結(jié)束語

本文構(gòu)建了一種二維引射循環(huán)預(yù)冷回流低溫貯箱內(nèi)氣液兩相流動(dòng)與傳熱特性模型，分析了地面狀態(tài)循環(huán)回流氣液兩相流對(duì)貯箱內(nèi)相分布、溫度、壓力和速度的影響，初步比對(duì)引射循環(huán)預(yù)冷與自然循環(huán)預(yù)冷回流貯箱特性的不同，結(jié)合工程實(shí)際分析影響引射循環(huán)預(yù)冷回流貯箱特性的關(guān)鍵因素和優(yōu)化原則，為工程應(yīng)用提供參考。

[1] Han V N. Ground pre-pressurization by helium bubbling for cryogenic upper stages[R]. AIAA 2001-3833, 2001.

[2] Panzarella C H, Kassemi M. On the validity of purely thermodynamic descriptions of two-phase cryogenic fluid storage[J]. Journal of Fluid Mechanics, 2003(484): 41-68.

[3] Cho N Y, Kwon O, Kim Y G, Jeong S K. Investigation of helium injection cooling to liquid oxygen propellant chamber[J]. Cryogenics, 2006(46): 132-142.

[4] 陳士強(qiáng), 范瑞祥, 黃兵, 黃輝. 運(yùn)載火箭低溫動(dòng)力系統(tǒng)注氣式循環(huán)預(yù)冷過程的AMESim仿真研究[J]. 載人航天, 2014, 20(5): 413-420.

[5] 程向華, 厲彥忠, 陳二鋒. 不同回流位置液體火箭發(fā)動(dòng)機(jī)循環(huán)預(yù)冷回路特性[J]. 推進(jìn)技術(shù), 2008, 29(6): 646-650.

[6] 陳二鋒, 厲彥忠, 程向華. 液體火箭發(fā)動(dòng)機(jī)自然循環(huán)預(yù)冷回路的數(shù)值研究[J]. 航空動(dòng)力學(xué)報(bào), 2009, 24(3): 698-704.

Study on Flow and Heat Transfer Characteristics of Cryogenic Tank with Gas Injection Circulation Precooling Backflow

Chen Shi-qiang1, Fan Rui-xiang2, Zhang Shu-hui1,Huang Bing1, Huang Hui1
(1. Beijing Institute of Astronautical Systems Engineering, Beijing, 100076; 2. China Academy of Launch Vehicle Technology, Beijing, 100076)

The characteristics of flow and heat transfer in the tank of rocket are closely related. The oxygen system of the new generation rocket is widely used the technology of injection circulation precooling backflow. Several main physical processes of cryogenic tank with gas injection circulation precooling for launch vehicle propulsion system are introduced, and the flow and heat transfer characteristics are investigated based on the free surface tracing VOF model. The influences of circulation precooling gas-liquid two-phase backflow to cryogenic tank are analyzed, including phase distribution, temperature variation, pressure fluctuation, and circulation flow velocity. The optimization principles of gas injection circulation precooling engineering application are discussed based on the results of research.

Injection circulation precooling backflow; Cryogenic tank; Flow; Heat transfer

V421.33

A

1004-7182(2017)03-0036-05

10.7654/j.issn.1004-7182.20170308

2016-03-24；

2016-08-22

陳士強(qiáng)（1986-），男，博士，工程師，主要研究方向?yàn)檫\(yùn)載火箭動(dòng)力系統(tǒng)總體設(shè)計(jì)