邱中華,劉志杰,朱文杰,黃立鈉,陳志堅(jiān)

(1.上?？臻g推進(jìn)研究所,上海 201112;2.上?？臻g發(fā)動(dòng)機(jī)工程技術(shù)研究中心,上海 201112)

0 引言

表面張力貯箱依靠推進(jìn)劑管理裝置Propellant Manage Devices(PMD)蓄留和傳輸推進(jìn)劑,PMD由表面張力元件組成。按工作方式,表面張力元件可以分為基于隔板形成楔角的開放式結(jié)構(gòu)(包含導(dǎo)流板[1]、楔形海綿體[2]、貯液槽裝置[3]等)和基于毛細(xì)網(wǎng)表面張力壓差的封閉式結(jié)構(gòu)(包含氣泡陷阱[3-4]、啟動(dòng)籃[5]、復(fù)合通道[6]、篩網(wǎng)收集器[7]等)。板式表面張力貯箱(以下簡(jiǎn)稱“板式貯箱”)的PMD構(gòu)成以開放式表面張力元件為主,通常包括導(dǎo)流板與出口附近的蓄液器。

板式貯箱因?yàn)榻Y(jié)構(gòu)簡(jiǎn)單、零組件少、工藝性好、造價(jià)低、擴(kuò)展性好等優(yōu)點(diǎn),在國(guó)外得到了大量應(yīng)用。然而,由于地面環(huán)境不是微重力環(huán)境,板式貯箱不方便直接進(jìn)行地面驗(yàn)證,其設(shè)計(jì)主要依賴于理論計(jì)算或仿真計(jì)算。

本文首先研究了微重力環(huán)境下導(dǎo)流板進(jìn)行液體驅(qū)動(dòng)的原理,然后據(jù)此原理進(jìn)行了導(dǎo)流板的理論計(jì)算,之后進(jìn)行了中性浮力試驗(yàn),并將試驗(yàn)結(jié)果與理論計(jì)算的結(jié)果進(jìn)行了對(duì)比。

1 導(dǎo)流板液體傳輸原理

在失重狀態(tài)下,液體推進(jìn)劑之間的主要作用力是分子引力,對(duì)于液位表面,則表現(xiàn)為表面張力。導(dǎo)流板就是依據(jù)導(dǎo)流板金屬分子與液體分子的引力原理設(shè)計(jì)的。導(dǎo)流板是一種靠近貯箱內(nèi)壁的結(jié)構(gòu),區(qū)別于傳統(tǒng)篩網(wǎng)式PMD建立封閉的流動(dòng)通道,導(dǎo)流板建立開放的流動(dòng)通道,推進(jìn)劑沿著通道流動(dòng)。傳統(tǒng)液體推進(jìn)劑具有較好的浸潤(rùn)性,趨向于吸附在金屬結(jié)構(gòu)的夾角處,在導(dǎo)流板與貯箱壁交線周圍形成液帶。圖1為沿著斜坡上導(dǎo)流板內(nèi)液體流動(dòng)的示意圖[1]。

根據(jù)Laplace方程,如果Rup>Rdown,會(huì)驅(qū)動(dòng)液體向下游流動(dòng),驅(qū)動(dòng)壓差

(1)

式中：σ為液體的表面張力系數(shù);Rup為上游液帶表面的主曲率半徑;Rdown為下游液帶表面的主曲率半徑,另一主曲率半徑均可認(rèn)為是無窮大。

考慮液體沿導(dǎo)流板的定常流動(dòng),將導(dǎo)流板簡(jiǎn)化為一維流動(dòng)區(qū)域,液體應(yīng)滿足如下的控制方程：

1)連續(xù)性方程

(2)

2)動(dòng)量方程

-Adp-τwsdx-ρAadz=ρAudu

(3)

3)氣液界面狀態(tài)方程

(4)

式中:x為沿導(dǎo)流板流動(dòng)的路徑長(zhǎng)度；u為液體的流動(dòng)速度；A為流動(dòng)截面積；Q為體積流量；p為液體的靜壓力；τw為液體沿壁面流動(dòng)的黏性剪阻力；s為潤(rùn)濕周邊長(zhǎng)度；a為加速度大??；z為沿加速度方向的液位高度坐標(biāo)；R為液體表面的主曲率半徑(另一曲率半徑為無窮大)；ρ為液體的密度。

液體沿導(dǎo)流板的流動(dòng)一般為層流,可認(rèn)為

(5)

將式(2)～式(5)進(jìn)行整理可得

(6)

2 導(dǎo)流板液體傳輸能力理論計(jì)算

導(dǎo)流板的截面可取為很多形狀,比如單豎片式、雙豎片式、三折橫片式、復(fù)雜式等,目前應(yīng)用最為廣泛的就是單豎片式[8-9],因?yàn)閹追N結(jié)構(gòu)的液體供應(yīng)能力差別不大,而單豎片式的結(jié)構(gòu)最為簡(jiǎn)單,工藝性最好[10-11]。單豎片式的截面形狀如圖2所示。

圖2 單豎片式導(dǎo)流板截面形狀Fig.2 Simple vane cross section

根據(jù)圖2可以得到s，A與R的關(guān)系,針對(duì)貯箱內(nèi)導(dǎo)流板的設(shè)置,可以得到dz/dx的值,應(yīng)用四階Runge-Kutta法,可對(duì)式(6)進(jìn)行求解。

下面使用該方法對(duì)一種經(jīng)典原理型板式貯箱內(nèi)的導(dǎo)流板[1]液體流動(dòng)特性進(jìn)行計(jì)算。該原理型貯箱為球柱形結(jié)構(gòu),PMD主要為導(dǎo)流板結(jié)構(gòu),帶中心柱,貯箱形狀如圖3所示。貯箱的內(nèi)徑為φ762 mm,柱段長(zhǎng)度為381 mm,導(dǎo)流板寬度為51 mm?？紤]側(cè)向加速度大小為a=0.000 5g,貯箱總流量為Q=5 ml/s,主要由側(cè)面的一根導(dǎo)流板與中心柱供應(yīng),液體介質(zhì)為肼推進(jìn)劑,這與衛(wèi)星在進(jìn)行位保動(dòng)作時(shí)的小推力工作狀態(tài)相對(duì)應(yīng)。主液池見圖中左側(cè)。因?yàn)閹в兄行闹Y(jié)構(gòu),所以液體供應(yīng)的通道可認(rèn)為有2條:一條是沿圖中A→B→F到達(dá)液口;另一條是沿圖中A→C→D→E到達(dá)液口。

圖3 原理型板式貯箱示意圖Fig.3 A schematic diagram of vane type tank

計(jì)算時(shí),取液口處E點(diǎn)的x坐標(biāo)為xE=0,x坐標(biāo)的取向沿中心柱進(jìn)行,即中心柱頂端xD=1 143 mm,貯箱柱段中間部位xA=1 932 mm(381×3+381π/2+381/2=1 932),最終計(jì)算的終點(diǎn)F點(diǎn)(與E點(diǎn)重合)坐標(biāo)為xF=2 721 mm(381×3+381π+381=2 721)。

計(jì)算結(jié)果如圖4所示。圖中不同的曲線代表不同的初始主曲率半徑。初始的主曲率半徑指流動(dòng)源頭(主液池)處的主曲率半徑。對(duì)圖3而言,即為A處的主曲率半徑。

圖4 原理型貯箱導(dǎo)流板理論計(jì)算結(jié)果 Fig.4 Theoretical calculation result of a vane ina principle vane type tank

可見,在側(cè)向加速度的條件下,液池內(nèi)液體沿2條路徑流向液口時(shí),液面處的主曲率半徑均逐漸減小。隨著貯箱內(nèi)推進(jìn)劑剩余量的減少,可以提供的初始主曲率半徑減小,導(dǎo)致液體流通路徑的主曲率半徑均減小,最終達(dá)到液體供應(yīng)的臨界點(diǎn)。這與實(shí)際的液體供應(yīng)情況是一致的。不同的液體剩余量對(duì)應(yīng)不同的R初始值,當(dāng)R初始值達(dá)到56 mm時(shí),達(dá)到了臨界供應(yīng)的情況,當(dāng)初始主曲率半徑小于該值時(shí),無法保證液體供應(yīng)到液口處。

理論計(jì)算方法能否應(yīng)用到實(shí)際的工程設(shè)計(jì)中,還要考慮理論計(jì)算方法與實(shí)際試驗(yàn)結(jié)果的一致性。

3 針對(duì)導(dǎo)流板的中性浮力試驗(yàn)

地面環(huán)境不是微重力環(huán)境,重力是決定液體流動(dòng)的主要因素,分子引力對(duì)液體作用的影響顯現(xiàn)不出來。因此,采用導(dǎo)流板結(jié)構(gòu)的板式貯箱在地面進(jìn)行性能試驗(yàn)比較困難[12]。目前常用的試驗(yàn)有落塔試驗(yàn)、拋物線飛行、空間搭載試驗(yàn)、中性浮力試驗(yàn)等[13-14]。各種方法均有其適用性與限制條件。考慮到中性浮力試驗(yàn)觀察時(shí)間不受限制的優(yōu)點(diǎn),采用中性浮力試驗(yàn)[15-16]進(jìn)行導(dǎo)流板在微重力條件下液體蓄留特性的試驗(yàn)驗(yàn)證。

3.1 應(yīng)用于表面張力貯箱的中性浮力試驗(yàn)原理

中性浮力試驗(yàn)的基本原理是用浮力抵消重力,廣泛應(yīng)用于航天員微重力環(huán)境的適應(yīng)性訓(xùn)練。

中性浮力試驗(yàn)設(shè)備應(yīng)用于表面張力貯箱時(shí),利用相似準(zhǔn)則達(dá)到對(duì)微重力工況的模擬[17-18]。通常利用一種液體模擬推進(jìn)劑液體,另一種液體模擬增壓氣體。為了模擬微重力工況,需要兩種液體的密度接近,且不相溶。另外,還要考慮液體推進(jìn)劑與PMD材料的浸潤(rùn)性,通常液體推進(jìn)劑與PMD具有較好的浸潤(rùn)性,而增壓氣體與PMD不浸潤(rùn)[19-20]?？紤]到液體密度是溫度的函數(shù),因此中性浮力試驗(yàn)作為表面張力貯箱模擬試驗(yàn)時(shí),對(duì)溫度要求較高。

中性浮力試驗(yàn)應(yīng)用于液體微重力條件下的液位分布時(shí),目前主要進(jìn)行靜平衡狀態(tài)的穩(wěn)態(tài)試驗(yàn),需要考慮重力與表面張力的相似性,還要考慮幾何相似。幾何相似除了指定幾何比尺外,還要考慮接觸角的一致性。重力與表面張力的相似用密度差

(7)

式中l(wèi)為幾何尺寸。應(yīng)用相似準(zhǔn)則時(shí),主要涉及到4種介質(zhì)：實(shí)際推進(jìn)劑、實(shí)際增壓氣體、模擬推進(jìn)劑的液體及模擬增壓氣體的液體。

在地面環(huán)境下,a′=g,g=9.8 m/s2,根據(jù)BOND數(shù)相等,可得

(8)

如果幾何尺寸沒有縮比,則上式可簡(jiǎn)化為

(9)

式中：下標(biāo)l為推進(jìn)劑液體;下標(biāo)g為增壓氣體;上標(biāo)′為模型參數(shù)。由式(8)和式 (9)可知,模擬的微重力加速度主要取決于如下幾個(gè)量：2種模擬介質(zhì)的密度、推進(jìn)劑的密度、模擬推進(jìn)劑的液體的表面張力系數(shù)及推進(jìn)劑的表面張力系數(shù)。考慮到介質(zhì)參數(shù)隨溫度的變化,要求中性浮力試驗(yàn)過程中對(duì)溫度的控制較為精確。實(shí)際推進(jìn)劑與貯箱壁面是浸潤(rùn)的,實(shí)際增壓氣體與壁面是不浸潤(rùn)的,為了模擬該接觸角特性,需要采用與模擬增壓氣體的介質(zhì)不浸潤(rùn)的材料制作貯箱殼體。本試驗(yàn)用水模擬增壓氣體氦氣,采用超薄F46貯箱殼體,用二苯甲烷模擬推進(jìn)劑肼,兩者的密度隨溫度變化關(guān)系如圖5所示。

圖5 兩種介質(zhì)的密度關(guān)系Fig.5 Density relationship of two medium

現(xiàn)將模擬條件設(shè)定為：推進(jìn)劑肼,增壓氣體氦氣,推進(jìn)劑溫度20 ℃?？傻媚M的重力加速度大小與溫度的對(duì)應(yīng)關(guān)系如表1所示。

表1 模擬的加速度值與溫度的關(guān)系

3.2 中性浮力試驗(yàn)結(jié)果

中性浮力試驗(yàn)用貯箱為直徑φ200系列的透明貯箱,分為球形與球柱形2種,2種貯箱對(duì)應(yīng)的導(dǎo)流板結(jié)構(gòu)如圖6所示。

分別進(jìn)行了多種溫度工況的中性浮力靜平衡試驗(yàn),試驗(yàn)結(jié)果如下所示。

圖6 試驗(yàn)用導(dǎo)流板結(jié)構(gòu)形式Fig.6 Vane structure for the experiment

1)球形貯箱,軸向加速度0.027 5 m/s2。

試驗(yàn)平衡溫度33.085 ℃,不考慮幾何縮比,根據(jù)式(9)計(jì)算可以得到對(duì)應(yīng)的加速度值0.027 5 m/s2。試驗(yàn)實(shí)物如圖7所示(為了便于觀察,二苯甲烷添加了微量的染色劑)。

圖7 中性浮力試驗(yàn)貯箱實(shí)物圖(0.027 5 m/s2)Fig.7 The tank used in buoyancy experiment(0.027 5 m/s2)

根據(jù)液位寬度b(指在導(dǎo)流板外側(cè)測(cè)量的液位橫向?qū)挾?如圖8所示)的實(shí)測(cè)值與導(dǎo)流板高度,計(jì)算得到了對(duì)應(yīng)位置(對(duì)應(yīng)于計(jì)算位置坐標(biāo)x處)的液位的主曲率半徑R,如表2所示。

表2 試驗(yàn)數(shù)據(jù)表

根據(jù)已知的x坐標(biāo)與主曲率半徑R值,可得到R的擬合曲線,同時(shí)根據(jù)理論計(jì)算方法可以得到R的變化規(guī)律(理論計(jì)算時(shí)各量的關(guān)系如圖8所示),兩者的對(duì)比如圖9所示。圖9中給出了一種試驗(yàn)曲線與兩種計(jì)算曲線,2種計(jì)算曲線對(duì)應(yīng)不同的加速度值。

圖8 理論計(jì)算各物理量關(guān)系示意圖 Fig.8 The relationship of parameters for thecalculation

圖9 試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比(0.027 5 m/s2) Fig.9 Contrast of experiment and calculationresults (0.027 5 m/s2)

2)球形貯箱,軸向加速度0.007 m/s2。

試驗(yàn)平衡溫度35.8 ℃,根據(jù)式(9)可求得對(duì)應(yīng)的加速度值0.007 m/s2。試驗(yàn)實(shí)物如圖10所示。

試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比如圖11所示。圖11給出了一種試驗(yàn)曲線與兩種計(jì)算曲線,兩種計(jì)算曲線對(duì)應(yīng)不同的加速度值。

圖10 中性浮力試驗(yàn)貯箱實(shí)物圖(0.007 m/s2)Fig.10 The tank used in buoyancy experiment (0.007 m/s2)

圖11 試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比(0.007 m/s2) Fig.11 Contrast of experiment and calculationresults (0.007 m/s2)

3)球形貯箱,軸向加速度0.036 7 m/s2。

試驗(yàn)平衡溫度31.855 ℃,根據(jù)式(9)可求得對(duì)應(yīng)的加速度值0.036 7 m/s2。試驗(yàn)實(shí)物如圖12所示。

圖12 中性浮力試驗(yàn)貯箱實(shí)物圖(0.036 7 m/s2)Fig.12 The tank used in buoyancy experiment (0.036 7 m/s2)

試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比如圖13所示。圖13給出了一種試驗(yàn)曲線與2種計(jì)算曲線,2種計(jì)算曲線對(duì)應(yīng)不同的加速度值。

圖13 試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比(0.036 7 m/s2) Fig.13 Contrast of experiment and calculationresults (0.036 7 m/s2)

4)球柱形貯箱,軸向加速度0.039 4 m/s2。

試驗(yàn)平衡溫度31.49 ℃,根據(jù)式(9)可求得對(duì)應(yīng)的加速度值0.039 4 m/s2。試驗(yàn)實(shí)物如圖14所示。

圖14 中性浮力試驗(yàn)貯箱實(shí)物圖(0.039 4 m/s2) Fig.14 The tank used in buoyancy experiment (0.039 4 m/s2)

試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比如圖15所示。圖15中給出了一種試驗(yàn)曲線與2種計(jì)算曲線,2種計(jì)算曲線對(duì)應(yīng)不同的加速度值。

5)球柱形貯箱,側(cè)向加速度0.040 4 m/s2。

試驗(yàn)平衡溫度31.35 ℃,根據(jù)式(9)可求得對(duì)應(yīng)的加速度值0.040 4 m/s2。試驗(yàn)實(shí)物如圖16所示。

圖15 試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比(0.039 4 m/s2) Fig.15 Contrast of experiment and calculationresults (0.039 4 m/s2)

圖16 中性浮力試驗(yàn)貯箱實(shí)物圖(0.040 4 m/s2側(cè)向)Fig.16 The tank used in buoyancy experiment (0.040 4 m/s2 lateral)

試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比如圖17所示。圖中給出了一種試驗(yàn)曲線與2種計(jì)算曲線,2種計(jì)算曲線對(duì)應(yīng)不同的加速度值。

圖17 試驗(yàn)結(jié)果與計(jì)算結(jié)果對(duì)比(0.040 4 m/s2側(cè)向) Fig.17 Contrast of experiment and calculationresults (0.040 4 m/s2 lateral)

3.3 中性浮力試驗(yàn)分析

通過多種微重力工況的靜平衡中性浮力試驗(yàn),得到了試驗(yàn)曲線與計(jì)算曲線的對(duì)比關(guān)系。根據(jù)上述各條曲線的對(duì)比可以發(fā)現(xiàn),中性浮力試驗(yàn)的結(jié)果與計(jì)算程序的結(jié)果具有一定的偏差,但對(duì)微重力加速度取一定的系數(shù)后,試驗(yàn)結(jié)果與計(jì)算程度結(jié)果的一致性較好。文中所列曲線的系數(shù)統(tǒng)計(jì)如表3所示。

表3 試驗(yàn)結(jié)果統(tǒng)計(jì)表

板式表面張力貯箱中的導(dǎo)流板結(jié)構(gòu)通常應(yīng)用于10-2～10-4m/s2微重力加速度下,而上表所列加速度已經(jīng)達(dá)到了10-1m/s2,是因?yàn)楸敬斡?jì)算針對(duì)直徑200 mm的貯箱。根據(jù)相似準(zhǔn)則,中性浮力試驗(yàn)基于Bond數(shù)相等,根據(jù)式(7),200 mm貯箱直徑下10-1m/s2時(shí)的導(dǎo)流板蓄留能力相當(dāng)于1 000 mm貯箱直徑下4×10-3m/s2時(shí)的導(dǎo)流板蓄留能力。

可見,對(duì)試驗(yàn)結(jié)果對(duì)應(yīng)的加速度取3～4倍系數(shù)后,試驗(yàn)結(jié)果與計(jì)算結(jié)果的一致較好。關(guān)于3～4倍系數(shù),可以從如下2個(gè)方面理解：①理論計(jì)算有一定的簡(jiǎn)化;②中性浮力試驗(yàn)是相似試驗(yàn),主要考慮了Bond數(shù)相似,無法考慮其它的相似條件?；谏鲜鰞煞矫娴脑?造成了該系數(shù)的存在。

4 工程可用的導(dǎo)流板計(jì)算方法

本文針對(duì)板式貯箱導(dǎo)流板不便進(jìn)行地面驗(yàn)證的問題,一方面建立了一種理論計(jì)算方法,另一方面利用中性浮力試驗(yàn)進(jìn)行了地面驗(yàn)證。通過計(jì)算結(jié)果與試驗(yàn)結(jié)果的對(duì)比,發(fā)現(xiàn)設(shè)置一定的系數(shù)后,2種結(jié)果的一致性較好?？紤]到中性浮力試驗(yàn)的特點(diǎn),本文所述3～4倍系數(shù)是偏保守的。從工程保守的角度,可以認(rèn)為,對(duì)理論加速度值取3倍系數(shù)后,可用于導(dǎo)流板的工程設(shè)計(jì)與計(jì)算。

基于上述認(rèn)識(shí),編制了一套適用于導(dǎo)流板計(jì)算的工程軟件,可直接用于導(dǎo)流板的設(shè)計(jì),軟件界面如圖18所示。

圖18 軟件界面Fig.18 The software interface

5 結(jié)語

1)針對(duì)板式貯箱導(dǎo)流板不便進(jìn)行地面驗(yàn)證的問題,進(jìn)行了理論計(jì)算研究與中性浮力試驗(yàn)研究,通過計(jì)算結(jié)果與試驗(yàn)結(jié)果的對(duì)比迭代,形成了一種經(jīng)過修正的導(dǎo)流板液體傳輸能力計(jì)算方法。

2)將導(dǎo)流板液體傳輸能力計(jì)算方法進(jìn)行了程序化,形成了一種工程可用的便捷的導(dǎo)流板計(jì)算軟件。

3)通過本文研究,促進(jìn)了板式表面張力貯箱設(shè)計(jì)技術(shù)的發(fā)展,加深了從業(yè)人員對(duì)于導(dǎo)流板性能的認(rèn)識(shí)。