劉 暢，肖澤寧，胡聲超，陳 益，韓 露

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

0 引 言

在空間中，微重力效應(yīng)會(huì)對(duì)推進(jìn)劑貯箱中氣液界面產(chǎn)生影響，表面張力在決定界面位置的因素中占據(jù)主導(dǎo)地位，液體將貼附在貯箱壁面上，氣體則處于液體中間。推進(jìn)劑傳輸過(guò)程中如果出現(xiàn)夾氣現(xiàn)象，輕則會(huì)造成發(fā)動(dòng)機(jī)工作不穩(wěn)定，重則將導(dǎo)致發(fā)動(dòng)機(jī)完全失效。為了滿足空間發(fā)動(dòng)機(jī)及在軌推進(jìn)劑流體傳輸中液體不夾氣的要求，貯箱中需要設(shè)置推進(jìn)劑管理裝置，實(shí)現(xiàn)氣液分離及液體獲取。

葉片、海綿、篩網(wǎng)通道3種傳統(tǒng)液體獲取方式在常規(guī)推進(jìn)劑領(lǐng)域得到了成功應(yīng)用[1,2]，但應(yīng)用于低溫推進(jìn)劑時(shí)存在諸多困難。葉片無(wú)法滿足大流量供應(yīng)要求；海綿在橫向及負(fù)向過(guò)載下極易失效，且對(duì)于大型貯箱，海綿質(zhì)量過(guò)大；篩網(wǎng)通道液體獲取裝置（Liquid Acquisition Device，LAD）是目前具有較好應(yīng)用前景的低溫推進(jìn)劑液體獲取方式，但其機(jī)理復(fù)雜，存在泡破點(diǎn)低、過(guò)網(wǎng)損失大、篩網(wǎng)后液體易蒸發(fā)等技術(shù)難題，目前國(guó)內(nèi)外尚無(wú)可用的飛行產(chǎn)品，技術(shù)成熟度低[3,4]。

本文以空間發(fā)動(dòng)機(jī)及流體傳輸為應(yīng)用背景，開(kāi)展推進(jìn)劑液體獲取技術(shù)研究，根據(jù)篩網(wǎng)通道式LAD壓降計(jì)算模型，進(jìn)行了不同參數(shù)下篩網(wǎng)通道壓力損失的計(jì)算，制定了地面出流試驗(yàn)方案，對(duì)模型計(jì)算的準(zhǔn)確性進(jìn)行了驗(yàn)證。

1 地面過(guò)載環(huán)境下篩網(wǎng)通道LAD壓降模型

1.1 篩網(wǎng)通道LAD壓降

無(wú)論在地面環(huán)境下還是在空間微重力環(huán)境下，只要篩網(wǎng)兩側(cè)的壓差不超過(guò)篩網(wǎng)的泡破壓力，氣體就無(wú)法通過(guò)篩網(wǎng)，就可以實(shí)現(xiàn)氣液分離。因此，對(duì)于篩網(wǎng)通道LAD，保證純液相供應(yīng)的條件為篩網(wǎng)兩側(cè)氣液壓力差應(yīng)小于篩網(wǎng)的泡破點(diǎn)[5～7]，即：

造成篩網(wǎng)兩側(cè)氣體和液體產(chǎn)生壓力差的因素包括5個(gè)方面：液柱壓力、過(guò)網(wǎng)局部損失、通道內(nèi)的沿程損失、動(dòng)壓損失以及其他壓力損失，其他壓力損失主要包括由于振動(dòng)、推進(jìn)劑晃動(dòng)、初始時(shí)刻液體沖擊造成篩網(wǎng)孔變形等因素造成的壓力損失，即：

式中ΔPH為通道內(nèi)液柱壓力；ΔPA為篩網(wǎng)過(guò)網(wǎng)損失；ΔPF為通道流動(dòng)摩擦沿程損失；ΔPD為通道流動(dòng)速度頭損失；ΔPO為其他壓力損失。

考慮地面1g過(guò)載條件下，通道垂直放置在貯箱中的情形，液體出流方向向上，如圖1所示。通道截面為矩形。初始時(shí)刻，貯箱中充滿液體，通道完全浸沒(méi)在液體中，液體通過(guò)篩網(wǎng)進(jìn)入通道并從通道上部流出。當(dāng)通道完全浸沒(méi)在液體中時(shí)，篩網(wǎng)被完全濕潤(rùn)，通道中的流動(dòng)為單相液體流動(dòng)。隨著液體不斷出流，貯箱中液面不斷下降，篩網(wǎng)逐漸露出液面，篩網(wǎng)兩側(cè)氣體與液體的壓力差不斷增加，由于表面張力作用，篩網(wǎng)阻止氣體進(jìn)入通道內(nèi)部，通道中仍然充滿液體。當(dāng)液面下降至一定程度，篩網(wǎng)兩側(cè)壓差大于篩網(wǎng)的泡破壓后，氣體進(jìn)入通道，篩網(wǎng)通道LAD失效。

圖1 篩網(wǎng)通道壓力損失示意Fig.1 Pressure Drop of Screen Channel

1.2 液柱壓力

液柱壓力ΔPH與液相、氣相的密度和重力加速度g有關(guān)，通過(guò)式（3）計(jì)算：

式中ρL為液體的密度；ρG為氣體的密度；g為重力加速度。

1.3 篩網(wǎng)過(guò)網(wǎng)損失

篩網(wǎng)過(guò)網(wǎng)損失ΔPA為

式中A為篩網(wǎng)粘性修正系數(shù)；μ為設(shè)計(jì)溫度下液體的動(dòng)力粘度；ve為推進(jìn)劑過(guò)網(wǎng)速度；B為篩網(wǎng)慣性修正系數(shù)。

式（4）中，過(guò)網(wǎng)速度ve是指液體在篩網(wǎng)孔中的流動(dòng)速度，數(shù)值上和通道內(nèi)液體流速線性相關(guān)。等號(hào)右邊，第1項(xiàng)值遠(yuǎn)遠(yuǎn)大于第2項(xiàng)值，可忽略第2項(xiàng)損失，因此可簡(jiǎn)化為

式中C為過(guò)網(wǎng)損失修正系數(shù)；v為通道內(nèi)液體的流速。

式（5）中，C與液體過(guò)網(wǎng)損失相關(guān)，由篩網(wǎng)的固有結(jié)構(gòu)所決定。可通過(guò)流阻試驗(yàn)，獲得一組流量-流阻試驗(yàn)數(shù)據(jù)，采取曲線擬合的方法，將系數(shù)C求出。對(duì)于不同的試驗(yàn)介質(zhì)，帶入相應(yīng)的μ值即可求解過(guò)網(wǎng)流阻損失。

1.4 通道流動(dòng)摩擦沿程損失

通道流動(dòng)摩擦沿程損失ΔPF與通道內(nèi)液體的雷諾數(shù)、液體流速、液體密度以及通道的幾何尺寸有關(guān)，通過(guò)式（6）計(jì)算：

式中Re為通道內(nèi)液體的雷諾數(shù)；v為通道內(nèi)液體的流速；ν為液體的運(yùn)動(dòng)粘度；d為通道寬度。

1.5 通道流動(dòng)速度損失

通道流動(dòng)速度損失ΔPD通過(guò)式（7）計(jì)算：

式中為通道內(nèi)液體的質(zhì)量流量；W為通道截面的寬度；H為通道截面的長(zhǎng)度。

2 計(jì)算結(jié)果及分析

根據(jù)上述壓降模型，對(duì)篩網(wǎng)通道LAD在地面過(guò)載環(huán)境下的壓降進(jìn)行計(jì)算。計(jì)算中以水作為介質(zhì)，加速度取為重力加速度。篩網(wǎng)規(guī)格為325×2300荷蘭斜紋網(wǎng)，以此確定篩網(wǎng)過(guò)網(wǎng)損失修正系數(shù)。分別計(jì)算了不同的篩網(wǎng)露出高度下的各項(xiàng)壓力損失變化情況，圖2為液柱高度不斷增加時(shí)各項(xiàng)壓力損失及篩網(wǎng)通道LAD壓降變化情況，表1為不同流量、液柱高度下篩網(wǎng)通道LAD壓降計(jì)算結(jié)果，從通道完全浸沒(méi)到通道逐漸露出液面，通道露出高度逐漸增大。

圖2 篩網(wǎng)通道LAD壓降隨通道內(nèi)液柱高度變化Fig.2 Pressure Drop of Screen Channel LAD Versus Liquid Column Heights

表1 不同流量、液柱高度下LAD壓降計(jì)算結(jié)果Tab.1 Pressure Drop in Different Flow Rates and Column Heights

由圖2和表1可知，篩網(wǎng)通道LAD壓降主要受篩網(wǎng)過(guò)網(wǎng)損失和液柱壓力損失影響，速度損失和沿程損失的影響較小，從表1計(jì)算結(jié)果來(lái)看，在不同流量、不同液柱壓力下，速度損失和沿程損失均遠(yuǎn)小于篩網(wǎng)過(guò)網(wǎng)損失和液柱壓力損失。篩網(wǎng)過(guò)網(wǎng)損失和液柱壓力損失均隨篩網(wǎng)露出高度的增大而增大，液柱壓力與篩網(wǎng)露出高度成正比，而篩網(wǎng)過(guò)網(wǎng)損失的增長(zhǎng)則隨篩網(wǎng)露出高度的增加而逐漸加快。當(dāng)篩網(wǎng)露出高度接近篩網(wǎng)總高度時(shí)，篩網(wǎng)過(guò)網(wǎng)損失將大幅增加，導(dǎo)致篩網(wǎng)通道LAD壓降快速增大，這是由于篩網(wǎng)過(guò)網(wǎng)損失與過(guò)網(wǎng)速度近似成正比，隨著篩網(wǎng)露出高度不斷增大，介質(zhì)能夠流通的面積不斷減小，而介質(zhì)流量維持不變，因此過(guò)網(wǎng)速度不斷增大，特別是介質(zhì)中浸沒(méi)的篩網(wǎng)高度較小時(shí)，過(guò)網(wǎng)速度將快速增大，導(dǎo)致篩網(wǎng)過(guò)網(wǎng)損失大幅增加。

不同通道內(nèi)流量下液柱壓力、篩網(wǎng)過(guò)網(wǎng)損失計(jì)算結(jié)果如圖3所示，不同通道內(nèi)流量下篩網(wǎng)通道LAD壓降計(jì)算結(jié)果如圖4所示。過(guò)網(wǎng)損失與介質(zhì)過(guò)網(wǎng)速度有關(guān)，在篩網(wǎng)完全浸沒(méi)的情況下，過(guò)網(wǎng)速度和通道內(nèi)流量成正比，因此隨著流量增大，過(guò)網(wǎng)損失逐漸增大。液柱壓力不受通道內(nèi)流量影響。

圖3 通道內(nèi)流量對(duì)液柱壓力、過(guò)網(wǎng)損失的影響Fig.3 Influence of Flow Rates on Pressure Drop of Liquid Column Heights and Screen

圖4 通道內(nèi)流量對(duì)篩網(wǎng)通道LAD壓降的影響Fig.4 Influence of Flow Rates on Pressure Drop of Screen Channel LAD

由圖3、4可知，通道內(nèi)流量對(duì)LAD壓降影響顯著，實(shí)際工程應(yīng)用中，通道內(nèi)流量取決于發(fā)動(dòng)機(jī)流量或液體傳輸流量、貯箱數(shù)量以及貯箱內(nèi)通道的個(gè)數(shù)，為了保證LAD可靠實(shí)現(xiàn)氣液分離，需要合理地選取通道內(nèi)流量。降低流量可降低LAD壓降，提高氣液分離可靠性，降低系統(tǒng)流阻，但特定的發(fā)動(dòng)機(jī)流量或液體傳輸流量下，降低通道流量需要通過(guò)增加通道數(shù)量和增大篩網(wǎng)面積實(shí)現(xiàn)，這會(huì)增加系統(tǒng)質(zhì)量，因此需要進(jìn)行合理的設(shè)計(jì)。

3 試驗(yàn)結(jié)果及分析

3.1 試驗(yàn)系統(tǒng)方案

為驗(yàn)證地面過(guò)載環(huán)境下篩網(wǎng)通道LAD壓降模型的正確性，通過(guò)地面環(huán)境下的出流試驗(yàn)對(duì)模型進(jìn)行驗(yàn)證。出流試驗(yàn)采用倒置出流方案，通道出液口位于通道頂部，液體出流方向向上，試驗(yàn)原理如圖5所示。

圖5 液體獲取裝置出流試驗(yàn)原理Fig.5 Schematic Diagram of Outflow Test for LAD

采用倒置式出流方案有利于考察通道內(nèi)各項(xiàng)水力壓力損失，為驗(yàn)證及修正模型提供數(shù)據(jù)支撐[3]。LAD縱向放置在模擬貯箱中，試驗(yàn)前先潤(rùn)濕篩網(wǎng)，并使液面沒(méi)過(guò)篩網(wǎng)，介質(zhì)充填LAD內(nèi)部。通過(guò)增壓氣體控制裝置向模擬貯箱增壓，使介質(zhì)以一定流量通過(guò)LAD及出流管路流出模擬貯箱，隨著模擬貯箱中液面不斷下降，篩網(wǎng)逐漸露出液面。在模擬貯箱及通道內(nèi)部設(shè)置壓力測(cè)點(diǎn)，模擬貯箱壓力測(cè)點(diǎn)位于模擬貯箱頂部，通道內(nèi)壓力測(cè)點(diǎn)位于通道頂部并通過(guò)測(cè)壓管路引至模擬貯箱外部，試驗(yàn)前測(cè)壓管路需要充滿介質(zhì)?？紤]到LAD的壓降特性與介質(zhì)的密度和粘性有關(guān)，因此可用水作為介質(zhì)開(kāi)展研究，其規(guī)律可拓展至低溫推進(jìn)劑。

3.2 篩網(wǎng)通道LAD方案

篩網(wǎng)為標(biāo)準(zhǔn)件產(chǎn)品，目前常用的有200×1400和325×2300兩種規(guī)格的不銹鋼荷蘭斜紋網(wǎng)。325×2300規(guī)格篩網(wǎng)等效（最大）孔徑約15 μm，對(duì)應(yīng)的泡破壓（無(wú)水乙醇）不小于6000 Pa。水的表面張力約為酒精的3.27倍，因此325×2300規(guī)格篩網(wǎng)用水作為工質(zhì)的泡破壓應(yīng)不小于19636 Pa。本文試驗(yàn)中LAD篩網(wǎng)選取325×2300規(guī)格的不銹鋼荷蘭斜紋網(wǎng)，圖6為325×2300規(guī)格篩網(wǎng)在掃描電鏡下的形貌。

圖6 325×2300規(guī)格篩網(wǎng)在掃描電鏡下的形貌Fig.6 Scanning Electron Microscopy Image of a 325×2300 Dutch Twill Screen

篩網(wǎng)通道LAD設(shè)計(jì)為長(zhǎng)方體形式，篩網(wǎng)固定在其中一面，通道截面為長(zhǎng)方形，尺寸為50 mm×40 mm，通道壁厚3 mm，材質(zhì)為鋁合金。篩網(wǎng)安裝采用可互換設(shè)計(jì)，可實(shí)現(xiàn)不同規(guī)格篩網(wǎng)試驗(yàn)。篩網(wǎng)與通道安裝處進(jìn)行密封設(shè)計(jì)，保證篩網(wǎng)泡破失效前不泄漏氣體。篩網(wǎng)通道LAD及在模擬貯箱中的安裝照片如圖7所示。由于篩網(wǎng)易受多余物堵塞，影響其流阻，因此試驗(yàn)前使用優(yōu)于篩網(wǎng)過(guò)濾精度的過(guò)濾器對(duì)試驗(yàn)介質(zhì)進(jìn)行多次過(guò)濾，保證試驗(yàn)結(jié)果的正確性。

圖7 篩網(wǎng)通道LAD及安裝Fig.7 Screen Channel LAD and Installation

3.3 試驗(yàn)結(jié)果分析

圖8為篩網(wǎng)通道LAD在地面1g過(guò)載下出流試驗(yàn)結(jié)果。橫坐標(biāo)為篩網(wǎng)露出高度，也就是通道內(nèi)部液柱高度。試驗(yàn)過(guò)程中，隨著液體不斷出流，實(shí)時(shí)記錄模擬貯箱氣枕壓力與通道內(nèi)頂部介質(zhì)的總壓，二者之差即為篩網(wǎng)通道LAD的壓降。由圖8可知，隨著篩網(wǎng)露出高度不斷增大，篩網(wǎng)通道LAD的壓降隨之不斷增大。篩網(wǎng)通道LAD壓降變化速率與篩網(wǎng)露出高度相關(guān)，以流量0.4 kg/s工況為例，篩網(wǎng)露出的初期，篩網(wǎng)通道LAD壓降變化速率較小，篩網(wǎng)露出高度超過(guò)400 mm后，篩網(wǎng)通道LAD壓降變化速率顯著增加。

圖8 325×2300篩網(wǎng)不同流量下的試驗(yàn)結(jié)果Fig.8 Test Result of 325×2300 Screen for Different Flow Rates

圖8中同時(shí)給出了本文理論模型的計(jì)算結(jié)果，可以看出流量較小時(shí)，理論模型計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合的較好。隨著流量增加，理論模型計(jì)算結(jié)果與試驗(yàn)結(jié)果之間出現(xiàn)一定偏差，這是通道內(nèi)速度分布的變化造成的。在通道內(nèi)部，沿著出流方向介質(zhì)的流動(dòng)速度是有一定分布的，在通道底部，也就是通道內(nèi)流動(dòng)的上游，介質(zhì)流動(dòng)速度小，而在通道的頂部，也就是流動(dòng)的下游，隨著液體不斷從篩網(wǎng)進(jìn)入，由于流量需要守恒，通道截面上介質(zhì)平均流動(dòng)速度不斷增大。這說(shuō)明介質(zhì)流動(dòng)速度在流動(dòng)方向上存在特定的分布，也就是介質(zhì)的動(dòng)壓存在沿出流方向的分布，這導(dǎo)致篩網(wǎng)的過(guò)網(wǎng)速度在沿出流方向上是不同的，進(jìn)而過(guò)網(wǎng)損失在沿出流方向上也是不同的。隨著篩網(wǎng)露出高度不斷變化，介質(zhì)的動(dòng)壓分布隨之產(chǎn)生變化，因此過(guò)網(wǎng)損失也隨之變化，本文的模型中假設(shè)過(guò)網(wǎng)速度沿出流方向分布是均勻的，在流量較小時(shí)，該假設(shè)引入的誤差較小，而當(dāng)流量增大時(shí)，雷諾數(shù)隨之增大，過(guò)網(wǎng)損失和速度損失的計(jì)算偏差逐漸增大，導(dǎo)致了試驗(yàn)結(jié)果與理論模型計(jì)算結(jié)果之間的偏差。這也是本文模型的后續(xù)改進(jìn)方向。

4 結(jié) 論

a）對(duì)于篩網(wǎng)通道LAD，保證純液相供應(yīng)的條件為篩網(wǎng)兩側(cè)氣液壓力差應(yīng)小于篩網(wǎng)的泡破點(diǎn)。篩網(wǎng)兩側(cè)氣液壓力差即篩網(wǎng)通道LAD壓降，造成壓力損失的因素包括液柱壓力損失、篩網(wǎng)過(guò)網(wǎng)損失、通道內(nèi)的沿程損失、動(dòng)壓頭損失以及其他壓力損失，其中主要影響因素為液柱壓力損失和篩網(wǎng)過(guò)網(wǎng)損失。

b）在地面環(huán)境下，篩網(wǎng)露出液面后，液柱壓力與篩網(wǎng)露出高度成正比，篩網(wǎng)過(guò)網(wǎng)損失主要受過(guò)網(wǎng)速度影響，隨著液面降低及浸沒(méi)在介質(zhì)中的篩網(wǎng)面積逐漸變小，過(guò)網(wǎng)損失顯著增加。

c）試驗(yàn)結(jié)果與模型計(jì)算結(jié)果在小流量下吻合較好，大流量下出現(xiàn)偏差，原因?yàn)橥ǖ纼?nèi)部速度分布的影響變得更為顯著。