胡 齊，陳 健，李 永

（北京控制工程研究所，北京100190）

微重力板式貯箱推進劑重定位的數(shù)值仿真

胡 齊，陳 健，李 永

（北京控制工程研究所，北京100190）

板式表面張力貯箱內(nèi)推進劑重定位對確定推進劑分布情況、研究晃動影響、提高控制精度等具有重要意義.為研究板式貯箱內(nèi)推進劑重定位的規(guī)律，對微重力下板式貯箱內(nèi)液體重定位問題進行數(shù)值仿真.計算時使用三維非定常兩相流動流體（VOF）模型，對某一板式貯箱壽命末期在不同微重力加速度情況下各種重定位過程進行數(shù)值仿真，得到各種工況下重定位的全過程，以及定位后推進劑的分布情況.數(shù)值仿真結(jié)果為板式貯箱的設計提供有利依據(jù).

微重力；板式表面張力貯箱；推進劑；重定位；數(shù)值仿真

表面張力貯箱是應用最為廣泛的衛(wèi)星推進劑貯箱，它利用表面張力進行液體輸送和氣液分離，為發(fā)動機或推力器提供不夾氣的推進劑［1］.從其推進劑管理裝置（PMD，propellantmanagement device）類型來分，表面張力貯箱可以分為篩網(wǎng)式和板式兩種.第一代表面張力貯箱都是篩網(wǎng)式的結(jié)構(gòu)，推進劑管理裝置使用許多毛細篩網(wǎng)，利用液體在篩網(wǎng)上的表面張力，通過篩網(wǎng)來收集推進劑，同時阻止氣體進入管路通道.第二代表面張力貯箱推進劑管理裝置則以板式結(jié)構(gòu)為主，盡量減少篩網(wǎng)的使用，利用流體在導流板上的表面張力驅(qū)動推進劑，對貯箱內(nèi)的推進劑實行全管理，可以滿足各種流量要求，并適應于不同的微重力環(huán)境，尤其是微重力相對較低的大型衛(wèi)星平臺［2］.

板式表面張力貯箱具有結(jié)構(gòu)簡單、質(zhì)量輕、可靠性高等一系列優(yōu)勢［3］.但是板式貯箱的性能不能在地面＋g環(huán)境下進行充分驗證，對其內(nèi)部流動規(guī)律的驗證主要依賴于數(shù)值仿真和微重力試驗驗證［4］.由于微重力試驗驗證費用昂貴而且需要通過數(shù)值仿真提供預估的結(jié)果，因此數(shù)值仿真在板式貯箱設計中具有極其重要的意義.在板式貯箱的設計中，對微重力環(huán)境下貯箱內(nèi)推進劑重定位的研究至關重要，該研究對確定推進劑分布情況、計算質(zhì)心位置、研究晃動影響、提高控制精度等都具有重要意義.由于條件所限，以往對推進劑重定位的研究偏重于最終定位結(jié)果，隨著計算流體力學的發(fā)展和計算機性能的提高，求解推進劑重定位全過程成為可能［5－6］.本文使用三維非定常流體的兩相流動模型，對某一板式貯箱壽命末期在不同微重力加速度情況下各種重定位過程進行數(shù)值仿真，得到各種工況下液體重定位的全過程，以及重新定位后推進劑的分布情況.

1 VOF模型基本原理

VOF模型是一種求解氣液兩相流動和自由表面流動的算法.最早由美國Los Alamos科學實驗室開發(fā)用于處理含有自由面的二維粘性流體的數(shù)值計算方法［7－8］.該方法通過引入流體體積組分函數(shù)及其控制方程來表示混合流體的密度并跟蹤自由面的位置.若設α為液相的體積組分，當網(wǎng)格中的體積組分α＝1時，表示該網(wǎng)格內(nèi)完全是液體，當網(wǎng)格中的體積組分α＝0時，表示網(wǎng)格中完全是氣體，當0＜α＜1表示網(wǎng)格中含有氣液兩相組分.

建立VOF模型時，有兩個關鍵問題需要解決：一個是對氣液交界面的重構(gòu)問題，另一個是交界面隨著時間的運移問題.以兩相流動為例，當在某一網(wǎng)格內(nèi)流體體積組分0＜α＜1時，說明這個網(wǎng)格是氣液界面，界面重構(gòu)就是指如何確定界面在網(wǎng)格中的位置.一般來說，VOF模型中界面重構(gòu)的方法有兩類：分段定常界面構(gòu)造法（PCIC，piecewise constant interface construction）和分段線性界面構(gòu)造法（PLIC，piecewise linear interface construction）.其中PCIC方法構(gòu)造精度較低；而PLIC方法是后來發(fā)展的，構(gòu)造精度較高.交界面隨時間的運移算法有兩種，即算子分裂算法和算子不分裂算法.算子分裂算法將界面從一個時刻到另一個時刻的運移過程看作界面分別在不同方向運動的合成；算子不分裂算法是直接根據(jù)網(wǎng)格中流體的當前位置和速度，判斷下一個時刻流體到達的位置，直接得到下一個時刻各個網(wǎng)格內(nèi)的流體分布.

（1）物性方程

流體的物性是由流體中不同相的體積組分確定的，物性方程就是不同體積組分流體物性的表達式.貯箱內(nèi)部流動只有氣液兩相，假設第二相的體積組分被跟蹤，那么流體的密度屬性方程為：

式中，ρ為混合流體的密度，α1、α2分別為第一相和第二相的體積組分；ρ1、ρ2分別為第一相和第二相的密度，它們?yōu)榻o定值.

（2）連續(xù)方程

流體連續(xù)方程的基本表達式如下：

式中，u為混合流體速度，R為源項.

劉雁衡睜開眼，說了一句令石警官很不滿的話。劉雁衡把凍僵的左手手指蜷起來，和右手一起把簫抱到胸前，淡淡對四小姐說：“對不起，你擋了我的風?！?/p>

根據(jù)物性方程（1）和連續(xù)方程（2）可以得到各相體積分數(shù)的運移方程：

式中，αq為第q相的體積組分，R是由凈相變引起的，R＝Re－Rc，其中Re表示蒸氣生成率，Rc表示蒸氣凝結(jié)率.在本模型中，由于R很小，故在求解時將此項可忽略不計.

（3）動量方程

混合流體的動量方程為：

式中，p為貯箱壓力，uj，ui分別為二相混合流的液相和氣相速度，xj，xi分別為二相混合流的液相和氣相位置，t為時間，gj為微重力加速度，μ為粘性系數(shù).由于表面張力的作用，在界面處有：

式中，σ為表面張力系數(shù)，r為液帶曲率半徑，聯(lián)立式（4）和式（5），得：

這就是混合流體的動量方程.

2 計算模型

計算選用如圖1所示的貯箱結(jié)構(gòu)，液體位于貯箱下半球蓄液器區(qū)域.采用結(jié)構(gòu)化網(wǎng)格對蓄液器區(qū)域進行網(wǎng)格劃分，如圖2和圖3所示.整個蓄液器內(nèi)計算區(qū)域網(wǎng)格數(shù)為203 866，整個蓄液器與貯箱計算區(qū)域網(wǎng)格數(shù)為178 368.

圖1 貯箱結(jié)構(gòu)圖

計算時，采用液態(tài)酒精作為模擬介質(zhì)，表面張力系數(shù)為0.022 03 N／m.假設整個蓄液器內(nèi)以及蓄液器與貯箱壁區(qū)域是處于密封狀態(tài)，將各個壁面都設為固壁邊界條件，即采用無穿透、無滑移邊界條件.

在初始情況下，假設模擬推進劑處于沉底狀態(tài).計算域包括液體和氣體兩部分，在進行計算時，先設定計算域內(nèi)的初始體積組分分布作為計算域內(nèi)氣液兩相分布的初始條件.初始化時設定液面高度，在該區(qū)域內(nèi)模擬推進劑的體積組分為1；其余區(qū)域模擬推進劑的體積組分為0.

圖2 蓄液器網(wǎng)格圖

圖3 蓄液器與貯箱壁區(qū)域網(wǎng)格圖

3 計算結(jié)果及分析

在壽命末期貯箱內(nèi)推進劑應該定位在蓄液器區(qū)域即貯箱出液口處，從而貯箱可以長時間為推力器提供不夾氣的推進劑.根據(jù)衛(wèi)星在軌工作情況，考慮零重力、10－3g微重力（－X、－Y方向），10－4g（－X方向）等微重力環(huán)境.

本文采用VOF兩相流動模型，模擬各種微重力情況下蓄液器內(nèi)和蓄液器與貯箱壁區(qū)域內(nèi)的推進劑的重定位過程，具體結(jié)果如下.

3.1 蓄液器內(nèi)推進劑重定位過程

（1）各種微重力加速度下液面的分布

圖4描述在各種微重力加速度下蓄液器內(nèi)重定位過程中液面的最終分布狀態(tài)，圖中深色部分代表模擬介質(zhì).

（2）零重力情況下重定位過程

在零重力情況下，推進劑將從沉底狀態(tài)沿著蓄液器弧壁向上緩慢爬升，并能夠到達蓄液器頂部.由于蓄液器頂部與蓄液器弧壁夾角較小，在表面張力的作用下有利于推進劑的蓄留，之后推進劑繼續(xù)沿著頂部先前緩慢爬升，但是由于頂部是平面而不是弧面，表面張力不足以支撐推進劑繼續(xù)先前爬升，最終重定位的狀態(tài)如圖4（a）所示.

（3）10－3g微重力（－Y方向）情況下重定位過程

在10－3g微重力（－Y方向）情況下，推進劑重定位過程存在振蕩的現(xiàn)象，推進劑在表面張力的作用下沿著蓄液器壁面向上爬升，在達到一定高度之后向下跌落，然后又沿著壁面向上爬升，如此振蕩多次，逐漸達到平衡狀態(tài)，形成如圖4（b）所示的凹面.推進劑重定位過程所引起的振蕩會造成整個貯箱質(zhì)心移動，這對衛(wèi)星的控制精度會產(chǎn)生一定影響，因此對整個振蕩過程的準確模擬能夠為提高衛(wèi)星控制精度提供有利參考價值.

（4）10－3g微重力（－X方向）情況下重定位過程

在10－3g微重力（－X方向）情況下，在微重力和表面張力作用下，推進劑從沉底狀態(tài)向－X方向運動，迅速充滿－X方向的半個蓄液器，并將其頂部淹沒，若此時頂部沒有密封，推進劑將從該位置進入蓄液器外部.然后液面形狀也開始振蕩，并最終趨于平衡，達到圖4（c）所示的狀態(tài).

圖4 不同微重力下蓄液器內(nèi)液面最終分布

3.2 蓄液器與貯箱壁區(qū)域內(nèi)推進劑重定位過程

（1）各種微重力加速度下液面的分布

圖5描述各種微重力加速度下蓄液器與貯箱壁區(qū)域重定位過程中液面的最終分布狀態(tài)，圖中深色部分代表模擬介質(zhì).

圖5 不同微重力下蓄液器與貯箱區(qū)域液面最終分布

（2）零重力情況下重定位過程

在零重力情況下推進劑在表面張力作用下將從沉底狀態(tài)沿著貯箱壁和蓄液器外壁弧面向上緩慢爬升，期間出現(xiàn)液面不斷的振蕩，直到表面張力不足以支撐推進劑繼續(xù)向上爬升，逐漸形成較規(guī)則的液面，最終重定位的狀態(tài)如圖5（a）所示.

（3）10－3g微重力（－Y方向）情況下重定位過程

在10－3g微重力（－Y方向）情況下，推進劑重定位過程仍然存在振蕩的現(xiàn)象，開始階段，推進劑在表面張力的作用下沿著貯箱壁面和蓄液器壁面向上爬升，在達到一定高度之后向下跌落，然后又沿著壁面向上爬升，如此振蕩多次，逐漸達到平衡狀態(tài)形成較規(guī)則液面，最終重定位狀態(tài)如圖5（b）所示.

（4）10－4g微重力（－X方向）情況下重定位過程

在10－4g微重力（－X方向）情況下，開始階段，由于微重力相對表面張力較小，推進劑在表面張力的作用下仍然沿著貯箱壁面和蓄液器壁面向上爬升，期間液面也仍然存在振蕩的現(xiàn)象，在達到一定高度之后，在微重力作用下，推進劑緩慢地向－X方向運動，液面同時存在振蕩，但由于微重力的值較小，振蕩幅度大幅減小，最終趨于平衡，達到圖5（c）所示的狀態(tài).

4 結(jié) 論

（1）通過數(shù)值仿真得到不同工況下推進劑重定位的全過程以及定位后推進劑的分布情況，結(jié)果表明這種板式PMD可以有效地蓄留推進劑，為推力器提供不夾氣的推進劑，驗證其對推進劑管理的良好性能.

（2）對不同工況下推進劑重定位過程的研究，可以為整個貯箱在衛(wèi)星不同工作狀態(tài)下推進劑的分布提供參考，進而為板式貯箱的設計提供有利依據(jù).

（3）推進劑重定位過程是個動態(tài)過程，存在反復振蕩現(xiàn)象.對推進劑重定位的過程數(shù)值仿真得到推進劑液面振蕩的全過程，進而為提高衛(wèi)星控制精度提供有益的參考.

［1］ Rollins JR，Grove R K，Jaekle Jr D E.Twenty－three years of surface tension propellant management system design，development，manufacture，test，and operation［R］.AIAA 85－1199

［2］ 李永，潘海林，魏延明.第二代表面張力貯箱的研究與應用進展［J］.宇航學報，2007，28（2）：503－507

［3］ Debreceni M J.W illiam D L，Jaekle Jr D E，Graffer A C.Design and development of the AXAF－IPS PMD＆PMD integration［R］.AIAA 97－2812

［4］ Netter G，Renner U，Dreyer M.Design and verification of a standard surface tension propellant tank［C］.The 35thJoint Propulsion Conference and Exhibit，Los Angeles，California，June 20－24，1999

［5］ 潘海林，洪儒珍.擋板對液體推進劑晃動力的影響［J］.中國空間科學技術，2000（1）：6－12

［6］ 潘海林，魏延明，洪儒珍.微重力環(huán)境中重力梯度加速度引起的流體晃動力和力矩的模擬［J］.宇航學報，1998，19（2）：66－72

［7］ Torrey M D，Cloutman L D，M jolsness R C.NASAVOF2D：a computer program for incompressible flows with free surfaces［M］.Los Alamos National Lab，LA－10612－MS，1985

［8］ Hirt CW，Nichols B D.Volume of fluid（VOF）method for the dynamics of free boundaries［J］.Journal of Computational Physics，1981（39）：201－225

Numerical Simulation for Reorientation of Propellant in a Vane Tank under Microgravity

HU Qi，CHEN Jian，LIYong
（Beijing Institute of Control Engineering，Beijing 100190，China）

Reorientation of propellant in a vane surface tension tank is of great significance for the ascertainment of propellant distribution，and research on liquid shaking mechanism，and improvement of control accuracy and so on.In order to investigate a law of propellant reorientation，propellant reorientation processes in vane type tank under different microgravity conditions are numerically simulated.The paper introduces a twophase flow model of three-dimensional unsteady VOF（volume of fluid），and based on the model numerical-ly simulates all kinds of reorientation processes in a vane type tank during near depletion phase under different microgravity acceleration conditions，then all processes of propellant reorientation under different working conditions and propellant distributions after reorientation are obtained.The simulation results provide advantageous foundation for design of the vane tank.

microgravity；vane surface tension tank；propellant；reorientation；numerical simulation

V43

A

1674-1579（2009）06-0054-04

2009-09-22

胡 齊（1985—），男，江西人，碩士研究生，研究方向為航天推進技術（e-mail：hqbuaa03＠sohu.com）.