微重力下空間循環離心泵空化性能研究

劉錦濤，李 永，史紀軍，梁紅義，羅先武

（1.北京控制工程研究所北京市高效能及綠色宇航推進工程技術研究中心，北京100094；2.北京衛星環境工程研究所，北京100094；3.清華大學熱能工程系，北京100084）

微重力下空間循環離心泵空化性能研究

劉錦濤1，李 永1，史紀軍2，梁紅義1，羅先武3

（1.北京控制工程研究所北京市高效能及綠色宇航推進工程技術研究中心，北京100094；2.北京衛星環境工程研究所，北京100094；3.清華大學熱能工程系，北京100084）

針對現有空化模型對微重力下復雜空化流動預測精度不高導致衛星在軌加注用循環離心泵的工作壓力低易發生空化的問題，建立了基于表面張力系數的空化兩相流計算模型預測微重力環境下循環離心泵空化流動，使用該模型計算了泵在空化時的性能，預測了泵內部發生空化時的空化形態以及空化體積，對比分析了常重力下的泵空化時的外特性以及空化流動形態。研究表明，微重力下表面張力的作用導致泵內部空化體積降低以及揚程的提高，微重力下表面張力在一定程度上可以抑制空化的發展。

微重力；離心泵；空化；性能

1 引言

我國從20世紀80年代末開始對采用空間循環泵流體回路實現主動熱控制進行研究［1］。神舟飛船載人任務的成功表明，單相流體回路系統能夠滿足熱控設計要求［2］。美國在火星探測中大量應用了單相流體回路熱控系統，其中以循環驅動泵熱控形式為主［3?4］。航天器中使用的熱控泵空蝕破壞現象突出，空化性能的好壞直接決定了熱空泵的使用壽命。衛星使用的熱控泵由于無法更換，其壽命大于載人飛船中使用的熱控泵，因此熱控泵的空化性能備受關注。

國際上已有很多采用空間循環泵對航天器進行主動熱控制的實例，在載人航天器上的應用較多，如：NASA的雙子星飛船、天空實驗室空間站、航天飛機，蘇聯的聯盟飛船、禮炮空間站，以及各國合作的國際空間站等［5?7］。在非載人航天器上運用空間循環泵取得成功的典型例子是NASA從20世紀90年代開始的一系列火星探測計劃：包括1997年發射的MPF（Mars Path?finder）、2003年發射的MER（Mars Exploration Rovers）以及2011年發射的MSL（Mars Science Laboratory），這些無人探測飛船上使用了噴氣推進實驗室（Jet Propulsion Laboratory，JPL）研制的空間循環泵［8?10］。國際上針對熱控泵的研究較多，但是針對微重力環境下熱控泵的空化性能研究較少。國內，熱控泵的研究仍處于探索階段，現在還未發現針對熱控泵微重力環境下空化性能研究的報道。

本文針對循環離心泵的流體驅動過程，建立微重力環境下循環離心泵空化計算模型，進而使用該模型探索泵內部空化流動特性，對比分析了微重力和常重力環境對循環離心泵內部流動的影響。

2 空化計算模型

基于均質混合流假設［11］，將液體與蒸汽混合成的空化流動混合物簡化為一種密度與粘性隨時間與空間變化的混合流體。

空化流動的基本控制方程如式（1）～（2）［12］：

式中，ρ、U（ui，uj）、p分別是流體的密度、速度和壓力，μm、μt分別是流體的分子粘性和湍流粘性，?為Level?set函數，δ為Dirac Delta函數，ρm為混合流體的密度。法向向量n＝Δ?／｜Δ?｜，界面曲率k＝Δ（Δ?｜Δ?｜）。方程式（2）中，σ為表面張力系數，而方程式等式右邊的最后一項為表面張力影響源項，表征了表面張力對空化流場的作用。

為了捕捉多種介質之間的界面，研究中采用了Level?set方法。其中，液體與蒸汽之間的界面記為S1，氣體與液體之間的界面記為S2。S1和S2界面以體積分數進行區分，在蒸汽與液體界面中，當液體體積分數分0時認定為S1界面，在液體與氣體界面，當蒸汽體積分數為0時認定為S2界面。Level set函數?則通過定義一個帶符號的距離來確定。對界面S1，Level?set函數可表示為式（3）［13］：

對界面S2，Level?set函數可表示為式（4）［13］：

式中d1指與蒸汽與液體分界面的距離，d2指與氣體與液體交界面的距離。

方程式（2）中表面張力影響源項可表示為式（5）：

在方程式（1）和（2）中，混合流體的物性參數可由方程式（6）和（7）確定［13］：

式中，下標“v”、“l”和“a”分別表示蒸汽、液體和氮氣（或其它氣體）。ε為小量規整函數，表征網格大小。下標“m”表示混合物。H（?）為光順化的Heaviside函數，定義為式（8）［13］：

對于蒸汽與氮氣（或其它氣體）之間的相變，可由式（9）～（10）所示質量輸運方程表示［13］：

方程式（9）中的相變源項可表示為式（11）～（12）［13］：

式中，Ce、Cc為基于實驗確定的經驗系數，α為體積分數，下標“nuc”表示非凝聚性氣體，Rb＝5e–5m，為初始狀態下氣核半徑。

方程式（1）～（12）組成求解本任務中循環離心泵內空化流動的基本方程。該方法考慮了表面張力和相間界面對空化發展的影響，適合于微重力條件下空化流動的求解。

3 空化流動計算方法

選取一循環離心泵模型，其設計轉速（8000 r／min）、設計流量（8 L／min）下對應的葉輪揚程為10 m，采用RANS方法來進行數值模擬。為了獲得更精確的結果，采用了k?ωSST湍流模型以及上述改進的空化模型，同時考慮了表面張力和接觸角的影響。，在微重力情況下，接觸角設置為80°；在常重力情況下，接觸角設置為90°［14］。

為了減少空化模擬對計算資源的影響，采用如圖1所示的單流道計算域。所有網格均為結構化網格。為了滿足k?ωSST模型的要求，y plus必須控制在100之內。

在空化流動方法研究中，以一種循環離心泵的設計方案為例，進行了水力性能與空化性能的比較。該方案中，葉輪出口直徑為34 mm，葉片數為5，出口角為30°，葉片包角為130°。

邊界條件的設置如下：

1）計算域的進口設置為全壓；

2）出口設置為質量流量；

3）在流道的兩側設置周期性條件；

4）其他為無滑移邊界條件。

循環離心泵的運行條件給定為：

1）轉速為8000 r／min，流量為8 L／min（0.48 m3／h）；

2）計算中，逐漸降低計算域進口的NPSH，直至揚程下降3%以上；

3）在常重力下，重力加速度設置為g0＝9.807 m／s2；而在微重力下，重力加速度設置為g＝0.000 01g0。

4 空化性能分析

為了探究表面張力對循環離心泵空化性能的影響，本文分別采用Level Set方法和原空化計算模型對循環離心泵空化性能進行了模擬。表1為用兩種不同方法在不同NPSH情況下數值模擬得出的空泡體積大小。從表1中可以發現，采用Level Set方法得到的空泡體積小于采用常規空化模擬方法得到的體積，表面張力項的引入提高了空化邊緣區域的壓力，減小了空化區域，從而可以推斷出表面張力抑制了空化的發展。這一結論與空泡動力學中得到的理論分析結果一致。

圖2顯示了在常重力和微重力情況下葉輪揚程和效率與流量的關系曲線。可以發現微重力情況下的揚程略微高于常重力下的揚程。然而，兩種情況下的效率曲線卻幾乎沒有區別。在設計轉速（8000 r／min）、設計流量（8 L／min）下對應的葉輪揚程為10.7 m，效率為65%。因此，重力加速度對循環離心泵的水力性能影響基本可以忽略不計。

表1 不同模擬方法得出的空泡體積比較Table 1 Cavitation volumeof different simulation methods

圖3 為循環離心泵在微重力條件下設計點即Q＝8 L／min時葉輪中截面的壓強分布，而圖4為在常重力條件下的壓強分布。可以發現，在不同的重力條件下，葉輪內部的總壓和靜態的分布是十分相似的。然而仍然存在一些差異，在微重力情況下葉片出口處的壓強要高于常重力下的壓強，導致微重力下的揚程要稍微高于常重力情況下。除此之外，兩者的壓強分布非常類似，因此兩者的效率也幾乎相同。

圖5為離心泵揚程隨有效空化余量變化的曲線。在微重力情況下，NPSH小于7 m時，揚程就呈現出下降的趨勢，而在常重力情況下，只有NPSH小于3 m才呈現出揚程下降趨勢。可以從圖中計算出兩種重力環境下揚程隨NPSH下降的臨界值NPSHr。在微重力情況下的NPSHr為0.84 m，而常重力情況下的NPSHr為0.94 m。表2為不同重力情況下單流道內空泡體積總量。可以發現相同NPSH時，在微重力情況下的空泡總量總是略小于常重力情況下的空泡總體積。因此，可以推斷微重力情況下循環離心泵的空化性能將優于常重力情況下的空化性能。

表2 微重力與常重力下的空泡體積Tab le 2 Cavitation volume in m icrogravity and norm al gravity

圖6 、7分別為微重力與常重力情況下不同NPSH下葉輪空泡體積分布，空泡表面定義為空泡體積率α為10%的等值面。根據觀察可知，兩種重力加速度下的空化發展具有基本相同的規律。

當離心泵進口壓力逐漸降低時，葉輪內壓力隨之下降，葉輪內出現空化，此時空泡主要集中在葉片吸力面的頭部。隨NPSH值的進一步降低，葉片吸力面空化區域增大，空泡在流道內靠近葉片吸力面一側堆積并向流道中間擴展。隨空化的加劇，流道進口段已聚集較多空泡，吸力面空化區已延伸至葉片中部，此時葉輪流道內堆積了大量空泡，造成嚴重堵塞并導致過流面積受限，影響流體正常流動，使離心泵做功能力變差，揚程明顯下降。

通過比較可知，在相同的NPSH時，常重力環境下泵葉輪內的空泡體積稍大；而當NPSH＝1 m時，常重力環境下空泡朝葉輪流道中間發展的趨勢更明顯，從而導致更加顯著的揚程下降。

5 結論

1）新建立的空化模擬方法適合于微重力環境下循環離心泵空化性能預測。

2）在微重力與常重力加速度下，空化演變具有相同的趨勢：葉輪初生空化均發生在葉片吸力面頭部。隨著空化加劇，空泡在流道內靠近葉片吸力面一側堆積并向流道中間擴展。

3）從葉輪內空泡體積和NPSHr可以得出，循環離心泵在微重力情況的空化性能將優于常重力情況下的空化性能。

（References）

［1］ 張立，姜軍，趙啟偉.單相流體回路在衛星熱控制中的應用研究［J］.中國空間科學技術，2008，28（3）：65?71. Zhang L，Jiang J，Zhao QW.Application of single?phase flu?id loop system onthermal control of satellite［J］.Chinese Space Science and Technology，2008，28（3）：65?71.（in Chi?nese）

［2］ 姜軍，張寧莉，徐濟萬.載人飛船主動流體循環回路離心泵可靠性試驗研究［J］.中國空間科學技術，2004，24（2）：33?37. Jiang J，Zhang N L，Xu JW.Reliability test study of eccen?tric pump in the active fluid cooled loop system of manned spacecraft［J］.Chinese Space Science and Technology，2004，24（2）：33?37.（in Chinese）

［3］ Bhandari P，Birur G，Karlmann P，et al.Mars science labo?ratorymechanically pumped fluid loop for thermal control?de?sign，implementation，and testing［J］.SAE International Journal of Aerospace，2009，4（1）：299?310.

［4］ Birur G C，Bhandari P，Gram M B，et al.Integrated pump assembly?an active cooling system for Mars pathfinder thermal control［R］.SAE Technical Paper961489，1996.

［5］ Birur G C，O’Donnell T P.Advanced thermal control tech?nologies for space sciencemissions at Jet Propulsion Laborato?ry［C］／／Space Technology and Applications International Fo?rum?2001.AIPPublishing，2001，552（1）：263?270.

［6］ ESA?technology?thermal control：single?phase mechanically?pumped fluid loops［EB／OL］.［2007?03?27］.http：／／www. esa.int／TEC／Thermal＿control／SEM4CLBE8YE＿0.html.

［7］ Van Benthem R C，de Grave W，van Es J，et al.Develop?ment of amechanically pumped fluid loop for 3 to 6 kW pay?load cooling［R］.SAE Technical Paper2009?01?2350，2009.

［8］ 張立，姜軍，趙啟偉.衛星用兩種單相流體回路方案對比研究［J］.航天器工程，2010，19（4）：35?40. Zhang L，Jiang J，Zhao Q W.Comparison between two schemes of single?phase fluid loop system for satellite［J］. Spacecraft Engineering，2010，19（4）：35?40.（in Chinese）

［9］ 聞建龍，沙毅，李傳君，等.無過載超低比轉速離心泵水力設計［J］.江蘇大學學報：自然科學版，2003，24（6）：21?24. Wen JL，Sha Y，LiC J，etal.Hydraulic design of non?over?load，super?low specific speed Centrifugal Pump［J］.Jour?nalofJiangsuUniversity（NaturalScienceEdition），2003，24（6）：21?24.（in Chinese）

［10］ 朱祖超.小流量高揚程泵的特點與應用［J］.水泵技術，1998（5）：10?12. Zhu Z C.Characteristic and application of small flow rate and high head pump［J］.Pump Technology，1998（5）：10?12.（in Chinese）

［11］ Echouchene F，Belmabrouk H，Le Penven，L，etal.Numer?ical simulation ofwall roughnesseffects in cavitating flow［J］. International journal of heat and fluid flow，2011，32（5）：1068?1075.

［12］ Ji B，Luo XW，Wu Y L，Peng X X and Duan Y L.Numeri?cal analysis of unsteady cavitating turbulent flow and shedding horseshoe vortex structure around a twisted hydrofoil［J］.In?ternational Journal of Multiphase Flow，2013，51：33?43.

［13］ Sussman M，Smereka P，Osher S.A level set approach for computing solutions to incompressible two?phase flow［J］. Journal of Computational physics，1994，114（1）：146?159.

［14］ Wang T，Li H X，Zhao J F，Three?Dimensional Numerical Simulation of Bubble Dynamics in Microgravity under the In?fluence of Nonuniform Electric Fields［J］.Microgravity sci?ence and technology，2016，28（2）：133?142.

Research on Performance of Space Cycling Centrifugal Pum p in M icrogravity

LIU Jintao1，LIYong1，SHIJijun2，LIANG Hongyi1，LUO Xianwu2
（1.Beijing Institute of Control Engineering，Beijing Engineering Research Center of Efficient and Green Aerospace Propulsion Technology，Beijing 100094，China；2.Beijing Institute of Spacecraft Environment Engineering，Beijing 100094，China；3.Department Of Thermal Engineering，Tsinghua University，Beijing 100084，China）

Most commercial cavitation models are not accurate enough to predict the cavitating flow in microgravity.However，the centrifugal pump used in the refueling system of satellitemay cause cavitation when running under low pressure conditions.To study the cavitating flow in the centrifugal pump in space，a cavitationmodel considering the surface tension inmicrogravity was founded.The performance of the pump with cavitation was calculated and the cavitating flow in space centrifugal pump was investigated in detail.The cavitation performance of the pump inmicrogravity and in nor?mal gravity on the ground were compared.The results showed that the surface tension inmicrogravity reduced the volume of the cavitating flow and increased the head of the pump.The surface tension may restrain the development of the cavitation in a pump in microgravity.

microgravity；centrifugal pump；cavitation；performance

TQ051.7

A

1674?5825（2017）01?0033?06

2016?01?04；

2017?01?09

國家自然科學基金（51406010）

劉錦濤，男，博士，工程師，研究方向為微重力流體機械流動穩定性以及多相流。E?mail：liujintao86＠hotmail.com