王柏輝, 史寶魯, 鄧明君, 李英斌, 尚文錦

(1.宇航動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 西安 710043; 2.中國(guó)航天員科研訓(xùn)練中心人因工程重點(diǎn)實(shí)驗(yàn)室, 北京 100094;3.中國(guó)航天員科研訓(xùn)練中心, 北京 100094)

1 引言

在載人航天失重條件下,氣體和液體將處于混流狀態(tài),不能像地面重力場(chǎng)條件下自動(dòng)實(shí)現(xiàn)氣體和液體分離[1-2],但在航天器中又廣泛存在對(duì)水、尿液、冷卻液、推進(jìn)劑等各類液體工質(zhì)的處理要求。因此無論是航天員飲用水,還是諸如系統(tǒng)控溫、尿收集、電解制氧等復(fù)雜的液體輸送和處理過程,都必須以通過氣液分離獲得液態(tài)流體為基本條件[3-4]。

微重力條件下的氣液分離技術(shù)通常分為靜態(tài)水氣分離技術(shù)[5]和動(dòng)態(tài)水氣分離技術(shù)[6],其原理均是通過一種力(如毛細(xì)力、離心力等)的作用來代替重力。與靜態(tài)分離相比,動(dòng)態(tài)分離的優(yōu)勢(shì)是在完成氣液分離的條件下,還可以提供部分管路循環(huán)壓力。旋片式動(dòng)態(tài)水氣分離器(Dynamic Gas-Liquid Separator with Rotary Impeller, DGLSRI)通過葉輪旋轉(zhuǎn)產(chǎn)生的離心力使氣液混合物旋轉(zhuǎn),從而實(shí)現(xiàn)不同密度的流體分離。該裝置是微重力環(huán)境下氣液分離的常用手段,其分離能力強(qiáng)、速度快、適用范圍廣[7-8],已應(yīng)用于環(huán)境控制與生命保障系統(tǒng)各氣液處理系統(tǒng)中。

目前,各國(guó)已研制出不同類型的動(dòng)態(tài)水氣分離器。1995 年,俄羅斯對(duì)應(yīng)用于和平號(hào)空間站尿液和衛(wèi)生廢水處理系統(tǒng)中的動(dòng)態(tài)水氣分離器進(jìn)行了性能測(cè)試研究[9]。該裝置具有出色的氣液分離效率,同時(shí)得到了不同出口流量下水的背壓,但并未對(duì)該裝置的功耗進(jìn)行約束。2002 年,美國(guó)研制了用于電解制氧系統(tǒng)的RSA(Rotary Separator Accumulator)[10],利用腔體內(nèi)的旋轉(zhuǎn)盤創(chuàng)建離心力場(chǎng)將氫氣與水分離[11],同時(shí)也集成了儲(chǔ)能器的功能,其缺點(diǎn)是裝置功耗較大。2005 年,美國(guó)研制了鼓型旋轉(zhuǎn)組件內(nèi)置葉片的RDS(Rotary Drum Separator),并將其用于Sabatier 二氧化碳還原系統(tǒng)中以去除水中的甲烷。該分離器采用2400 r/min 的轉(zhuǎn)速,將液體出口壓力提升到了138 kPa[12-14]。

當(dāng)前的旋片式動(dòng)態(tài)水氣分離器存在體積大、功耗高、輸送能力不足等問題,對(duì)此本文開展旋片式動(dòng)態(tài)水氣分離器輸送能力提升研究,提出準(zhǔn)穩(wěn)態(tài)環(huán)境下輸送能力的求解方法,通過仿真計(jì)算得到葉輪結(jié)構(gòu)等不同參數(shù)變化對(duì)輸送能力的影響;同時(shí)將仿真分析與試驗(yàn)相結(jié)合,提出可靠的旋片式動(dòng)態(tài)水氣分離器輸送能力提升方法,以期為該裝置的設(shè)計(jì)優(yōu)化提供理論依據(jù)。

2 設(shè)計(jì)方法

2.1 旋片式動(dòng)態(tài)水氣分離器工作原理

旋片式動(dòng)態(tài)水氣分離器工作原理如圖1 所示。該裝置特點(diǎn)是通過旋轉(zhuǎn)的葉輪結(jié)構(gòu)使氣液混合物獲得離心加速度。由于液體密度大,受到的離心力大,在葉輪的攪動(dòng)下被甩到腔體壁面上形成水膜。開始工作時(shí),液體出口為關(guān)閉狀態(tài),氣體出口始終打開。氣液混合物從靠近轉(zhuǎn)軸的入口流入,在葉輪的作用下實(shí)現(xiàn)氣液混合物的旋轉(zhuǎn)分離。當(dāng)液體水膜達(dá)到一定厚度時(shí),液體出口打開,液體通過腔體外緣處的排液通道排出,氣體從旋軸中心的排氣通道流出。旋轉(zhuǎn)葉輪是實(shí)現(xiàn)氣液分離的核心組件,主要起到捕獲液體、傳導(dǎo)離心力的作用,實(shí)現(xiàn)機(jī)械能向液體能量的轉(zhuǎn)化。

圖1 旋片式動(dòng)態(tài)水氣分離器工作原理示意圖Fig.1 Schematic diagram of working Principle of DGLSRI

2.2 輸送能力分析方法

旋片式動(dòng)態(tài)水氣分離器輸送能力提升的關(guān)鍵是提升其液體出口處的壓力。集液過程中,液體出口始終關(guān)閉,液體在腔體壁面積累形成液環(huán),氣體和液體之間出現(xiàn)交界面,如圖2 所示。葉輪轉(zhuǎn)動(dòng)的角速度為ω,圓柱形腔體半徑為R2,葉輪中心到氣液交界面距離R1,h為水膜厚度。氣液混合物在葉輪的作用下以一定的角速度繞轉(zhuǎn)軸運(yùn)動(dòng),滿足能量守恒方程,常用伯努利方程來描述能量守恒關(guān)系。在旋轉(zhuǎn)坐標(biāo)系中,伯努利方程式可由式(1)表示。

圖2 DGLSRI 流場(chǎng)示意圖Fig.2 Schematic diagram of flow field of DGLSRI

式中,p1、p2表示圖2 中1、2 兩點(diǎn)的壓力,z1、z2表示單位質(zhì)量流體的位置勢(shì)能,u1、u2表示流道的牽連速度,w1、w2表示液體的相對(duì)速度,ρ為液體密度,g為重力加速度,hs,1-2表示能頭損失,常由摩擦等因素造成。

首先對(duì)模型進(jìn)行如下假設(shè):

1)氣液混合物在裝置腔體內(nèi)實(shí)現(xiàn)完全分離,即不考慮液體中含有的氣泡和氣體中的液滴;

2)忽略葉片與壁面之間的間隙;

3)將流體看作是理想流體;

4)引入無限葉片假設(shè),即葉片足夠多無限薄。因此腔體內(nèi)流體的流動(dòng)可視為軸對(duì)稱流動(dòng),流體的轉(zhuǎn)速與葉片轉(zhuǎn)速相同均為ω。

取液體通道的一條流線1 ～2,對(duì)于1、2 兩點(diǎn)可寫出理想流體的伯努利方程,見式(2)。

式中,ρw表示液體密度。由此可推出2 點(diǎn)處的壓力p2,見式(3)。

如果不考慮重力影響,則可簡(jiǎn)化為式(4)。

在分析時(shí)引入水膜厚度h的概念,即h=R2-R1,代入式(4)中可得式(5):

因此,由式(5)可得,增加葉輪直徑、葉輪轉(zhuǎn)速及水膜厚度均有利于提升該裝置的液體出口壓力。但此式未考慮縫隙及有限葉片的影響,因此不能作為最終結(jié)果計(jì)算。

DGLSRI 總功耗由兩部分組成,一部分是電機(jī)功耗,主要包括機(jī)械損失和電磁損失等;另一部分是腔體功耗,下文簡(jiǎn)稱功耗,這部分功耗發(fā)生在分離腔體內(nèi),包括葉輪對(duì)流體做功的功耗和流動(dòng)損失等。腔體功耗可以通過仿真計(jì)算得到,而電機(jī)功耗必須結(jié)合實(shí)驗(yàn)耦合求解。

本文將應(yīng)用計(jì)算流體力學(xué)(Computational Fluid Dynamics,CFD)仿真方法和實(shí)驗(yàn)方法研究動(dòng)態(tài)水氣分離器輸送能力提升方法。

2.3 CFD 求解方法

本文采用一種旋片式動(dòng)態(tài)水氣分離器輸送能力的CFD 求解方法[15],通過該方法可求解不同結(jié)構(gòu)、不同參數(shù)下該裝置的液體出口壓力、功耗和分離效率,該方法的核心在于模型的選擇。

DGLSRI 是處理航天器中氣液混合物的裝置,其內(nèi)部存在著復(fù)雜的氣液流動(dòng)。在計(jì)算DGLSRI 腔體內(nèi)部湍流運(yùn)動(dòng)時(shí),由于其運(yùn)動(dòng)的復(fù)雜性,因此采用建立湍流模型結(jié)合Reynolds 時(shí)均方程計(jì)算的方法。腔體內(nèi)部的流動(dòng)主要是葉輪旋轉(zhuǎn)引起的湍流,且葉輪與壁面間縫隙的大小對(duì)最終流動(dòng)狀態(tài)影響很大。常用的湍流模型有k-ε 模型和k-ω 模型,k-ε 模型對(duì)邊界層的計(jì)算和剪切流的計(jì)算具有很好的適用性,k-ω 模型對(duì)低Re流動(dòng)、自由剪切流動(dòng)以及可壓縮流動(dòng)的計(jì)算效果更很好,但均未考慮湍流切應(yīng)力的運(yùn)輸影響。因此,采用剪切應(yīng)力運(yùn)輸(Shear-Stress Transport,SST kω)模型,通過將k-ε 模型和k-ω 模型以加權(quán)平均的方式組合起來,近壁處采用k-ω 模型,遠(yuǎn)壁處采用k-ε 模型,克服了原k-ω 模型中邊界層外ω對(duì)自由流相應(yīng)值的依賴性,對(duì)近壁區(qū)域流動(dòng)的計(jì)算更加精確。

流體體積(Volume of Fluid,VOF)模型是處理2 種或以上不相融流體流動(dòng)的表面跟蹤方法,通過它可獲得不同組分間的交界面及每個(gè)計(jì)算單元內(nèi)各組分的體積分?jǐn)?shù)。該裝置處理的兩相流為互不相融的流體,且VOF 模型對(duì)分界面的追蹤有利于確定臨界排水點(diǎn),對(duì)分離效率的提高具有重要意義。因此本文選用VOF 模型作為兩相流處理模型。

對(duì)于旋轉(zhuǎn)機(jī)械多流體區(qū)域,采用多參考系(Multiple Reference Frame,MRF)模型,該模型在假定流動(dòng)為定常的情況下,根據(jù)流體區(qū)域的實(shí)際旋轉(zhuǎn)情況進(jìn)行多個(gè)區(qū)域的劃分,對(duì)不同的流體區(qū)域應(yīng)用不同參考系求解計(jì)算。該模型可有效模擬旋轉(zhuǎn)機(jī)械內(nèi)部流動(dòng)。

為了使該裝置能夠更好地應(yīng)用于空間環(huán)境,對(duì)比了在重力條件下和微重力條件下的仿真結(jié)果,可得出在高轉(zhuǎn)速情況下,離心力的影響大于重力的影響,重力對(duì)該裝置性能參數(shù)的影響可以忽略不計(jì)。

綜上,根據(jù)CFD 求解方法,本文進(jìn)一步仿真分析不同影響因素對(duì)液體出口壓力和功耗的影響規(guī)律,進(jìn)而設(shè)計(jì)試驗(yàn)裝置,通過開展輸送能力提升試驗(yàn),以驗(yàn)證仿真分析的正確性;最終通過綜合仿真和試驗(yàn)結(jié)果,提出輸送能力的提升方法。

3 輸送能力影響因素仿真

根據(jù)理論分析和工程實(shí)際,葉輪(電機(jī))轉(zhuǎn)速、葉輪結(jié)構(gòu)和葉輪尺寸等均會(huì)影響旋片式動(dòng)態(tài)水氣分離器的輸送能力。通過CFD 求解方法分別對(duì)電機(jī)轉(zhuǎn)速、葉片數(shù)量和葉輪尺寸對(duì)輸送能力的影響展開仿真分析。

3.1 電機(jī)轉(zhuǎn)速影響

根據(jù)式(5),在引入無限葉片假設(shè)的情況下,液體出口壓力與轉(zhuǎn)速的平方成正比。但在有限葉片時(shí),由于渦旋和縫隙的影響,液體出口壓力與電機(jī)轉(zhuǎn)速將不是嚴(yán)格平方關(guān)系。通過仿真計(jì)算,得到電機(jī)轉(zhuǎn)速分別為1150 r/min、1600 r/min、2000 r/min、2200 r/min 和2400 r/min 時(shí)DGLSRI輸送能力仿真結(jié)果如圖3 所示,橫軸為水膜厚度比h/R,其中h為水膜厚度,R為腔體直徑。

圖3 不同電機(jī)轉(zhuǎn)速輸送能力仿真結(jié)果Fig.3 Simulation results of conveying capacity under different motor speeds

圖3(a)結(jié)果表明:①在電機(jī)轉(zhuǎn)速不變時(shí),隨水膜厚度的增加,液體出口壓力呈拋物性型增加,增長(zhǎng)速率逐漸減小;轉(zhuǎn)速越高,液體出口壓力變化曲線的斜率越大。②水膜厚度相同時(shí),液體出口壓力會(huì)隨電機(jī)轉(zhuǎn)速的增加而增大;水膜厚度相同時(shí),液體出口壓力約與電機(jī)轉(zhuǎn)速的平方成正比,與式(5)給出的變化規(guī)律基本相同。

圖3(b)結(jié)果表明:①轉(zhuǎn)速相同時(shí),功耗隨水膜厚度的增加而增大,但增長(zhǎng)速率逐漸減小;當(dāng)水膜厚度增大到一定值時(shí),隨水膜厚度的增加,功耗基本不變,達(dá)到最大值。②水膜厚度相同時(shí),電機(jī)轉(zhuǎn)速增加約1 倍,最大功耗值增大至原功耗的8倍左右。因此,電機(jī)轉(zhuǎn)速增加對(duì)裝置功耗的貢獻(xiàn)遠(yuǎn)大于其對(duì)液體出口壓力的貢獻(xiàn)。

綜上,僅采用增加電機(jī)轉(zhuǎn)速的方式來提高DGLSRI 的輸送能力是不可行的。電機(jī)轉(zhuǎn)速增加雖然會(huì)增大該裝置的液體出口壓力,但也會(huì)造成功耗的急劇增大,而航天器中對(duì)功耗限制嚴(yán)格。

3.2 葉片數(shù)量影響

為便于研究,2.2 節(jié)中引入了無限葉片假設(shè),氣液混合物流動(dòng)受到嚴(yán)格的約束。而實(shí)際情況中,葉輪中葉片數(shù)量有限且具有一定的厚度,這會(huì)導(dǎo)致流動(dòng)滑移現(xiàn)象。流動(dòng)滑移的存在會(huì)使腔體內(nèi)流體運(yùn)動(dòng)產(chǎn)生壓力的理論值減小。

葉片數(shù)量不同會(huì)使流體受到葉片的夾持作用不同,對(duì)輸送能力的影響也不同。對(duì)電機(jī)轉(zhuǎn)速2200 r/min、水膜厚度0.34R、葉片數(shù)分別為4 片、6 片和8 片的情況進(jìn)行仿真計(jì)算,可以得到不同葉片數(shù)量情況下液體出口壓力與功耗值,如表1所示。

表1 不同葉片數(shù)量下液體出口壓力與功耗結(jié)果Table 1 Results of liquid outlet pressure and power consumption under different blades

通過分析可知,葉片數(shù)量增多時(shí)葉片對(duì)流體的夾持作用增加,流體的隨動(dòng)性更好,流動(dòng)滑移對(duì)液體出口壓力的減弱作用在下降。表1 結(jié)果表明,葉片數(shù)量增加至8 片時(shí),在增大液體出口壓力的同時(shí)可以降低功耗,因此后續(xù)研究均圍繞8 葉片的葉輪結(jié)構(gòu)展開。

3.3 葉輪尺寸影響

葉輪直徑和葉片寬度是葉輪最重要的尺寸參數(shù),由式(5)可知,當(dāng)引入無限葉片假設(shè)并忽略縫隙時(shí),液體出口壓力與葉輪直徑成正比,但與葉片寬度無關(guān)。由于縫隙的存在且葉片數(shù)量是有限的,因此葉片寬度也會(huì)對(duì)液體出口壓力產(chǎn)生影響。

仿真條件如下:葉輪直徑為115 ～140 mm,葉片寬度為6 ～30 mm,葉片數(shù)量為8,電機(jī)轉(zhuǎn)速為2200 r/min,水膜厚度為25 mm。輸送能力仿真計(jì)算結(jié)果如圖4 所示。

圖4 不同葉輪尺寸輸送能力仿真結(jié)果Fig.4 Simulation results of conveying capacity under different impeller sizes

圖4(a)結(jié)果表明: ①固定葉片寬度時(shí),隨葉輪直徑的增加,液體出口壓力隨之增大。但由于葉輪與腔體壁面間縫隙粘性底層的存在以及腔體內(nèi)流體渦旋的影響,液體出口壓力隨葉輪直徑的關(guān)系曲線會(huì)呈現(xiàn)出一定的波動(dòng)。當(dāng)葉片寬度較小,如為5 mm 時(shí),液體出口壓力與葉輪直徑的關(guān)系曲線波動(dòng)較大。②固定葉輪直徑時(shí),葉片寬度的增加同樣會(huì)使液體出口壓力增大,這一結(jié)論在有限葉片數(shù)時(shí)成立。當(dāng)葉輪直徑不變,葉片寬度從30 mm 開始減小時(shí),液體出口壓力會(huì)隨之緩慢減小;當(dāng)葉片寬度減小至10 mm 以下時(shí),隨葉片寬度減小,液體出口壓力的下降程度在逐漸增大。這一現(xiàn)象在葉輪直徑較小時(shí)更為明顯。

圖4(b)結(jié)果表明: ①固定葉片寬度時(shí),隨葉輪直徑的增加,功耗不斷增大,功耗的增長(zhǎng)速率變大。葉片寬度較大時(shí),葉輪直徑增加,功耗的增長(zhǎng)速率快;而葉片寬度較小時(shí),功耗增長(zhǎng)速率較為緩慢。因此,若要通過增大葉輪直徑增加裝置的輸送能力時(shí),應(yīng)盡量選擇小葉片寬度以免功耗過大,超過額定功耗。②固定葉輪直徑時(shí),葉片寬度的增加也會(huì)使功耗增大。同樣,當(dāng)葉輪直徑較大時(shí),葉片寬度增加,功耗增長(zhǎng)較快;當(dāng)葉輪直徑較小時(shí),隨葉片寬度增加,功耗增長(zhǎng)的較為緩慢。

通過上述分析,增大葉輪直徑和葉片寬度均會(huì)同時(shí)增加液體出口壓力和功率,且葉輪直徑對(duì)輸送能力的影響程度遠(yuǎn)大于葉片寬度的影響。

3.4 不同葉輪尺寸最大輸送能力仿真計(jì)算

DGLSRI 設(shè)計(jì)過程中有一項(xiàng)重要的技術(shù)要求,即保證分離效率、低于額定功耗的情況下可達(dá)到的最大液體出口壓力。為防止水從氣體出口流出,在水膜厚度達(dá)到一定值時(shí)需打開液路出口,輸送腔體中的液體。實(shí)際工程試驗(yàn)中發(fā)現(xiàn),當(dāng)氣液混合物中氣體流量小時(shí),液體積累至旋轉(zhuǎn)軸附近氣體出口處仍無液滴產(chǎn)生。因此,對(duì)不同葉輪尺寸下液體充滿至旋轉(zhuǎn)軸附近所能達(dá)到的液體出口壓力和功耗(即最大輸送能力)展開研究。根據(jù)3.3 節(jié)中的仿真條件,得到不同葉輪尺寸下該裝置可達(dá)到的最大輸送能力如圖5 所示。

圖5 不同葉輪尺寸最大輸送能力Fig.5 Maximum conveying capacity under different impeller sizes

圖5(a)結(jié)果表明: ①最大液體出口壓力與圖4(a)同水膜厚度下的液體出口壓力規(guī)律不同,改變?nèi)~片寬度對(duì)最大液體出口壓力的影響很小。固定葉輪直徑的情況下,當(dāng)葉片寬度大于10 mm時(shí),隨葉片寬度增加最大液體出口壓力基本不變;當(dāng)葉片寬度小于10 mm 時(shí),隨葉片寬度減小,最大液體出口壓力也會(huì)降低。②葉輪直徑的增加會(huì)大幅提升該裝置的最大液體出口壓力,且最大液體出口壓力與葉輪直徑呈線性關(guān)系。

雖然改變?nèi)~片寬度對(duì)最大液體出口壓力的影響很小,但會(huì)對(duì)功耗產(chǎn)生不可忽視的影響。圖5(b)結(jié)果表明: ①固定葉輪直徑的情況下,葉片寬度增加,最大功耗隨之增大,葉輪直徑為110 mm時(shí),葉片寬度從6 mm 增加到30 mm,最大功耗會(huì)增加24 W。這樣的增加量在實(shí)際設(shè)計(jì)過程中是不可忽略的。②葉輪直徑仍是影響最大功耗的主要因素。固定葉片寬度的情況下,最大功耗隨葉輪直徑的增加而增大,且葉輪直徑的變化對(duì)最大功耗的影響效果顯著。

因此,為盡可能降低最大功耗,需著重關(guān)注葉輪直徑尺寸的設(shè)計(jì)。與葉片寬度相比,葉輪直徑改變對(duì)輸送能力的影響更大。因此,在固定轉(zhuǎn)速的情況下,增加葉輪直徑是增大液體出口壓力的主要手段。而由此帶來的功耗升高可通過減小葉片寬度使其降低至額定功耗以下。

4 輸送能力提升試驗(yàn)

本文設(shè)計(jì)了一種DGLSRI 的裝置結(jié)構(gòu),模型及實(shí)物圖如圖6 所示,該裝置由電機(jī)、殼體、聯(lián)軸器、葉輪組件、接頭接口等部分組成。

圖6 試驗(yàn)裝置圖Fig.6 Drawing of the experimental apparatus

4.1 不同電機(jī)轉(zhuǎn)速條件下的試驗(yàn)結(jié)果

對(duì)不同轉(zhuǎn)速條件下的液體出口壓力和總功耗進(jìn)行測(cè)試,測(cè)試條件如下:葉輪直徑為115 mm,葉片寬度為10 mm,葉片數(shù)量為8 片,電機(jī)轉(zhuǎn)速為800 r/min、1200 r/min、1600 r/min、2000 r/min 和2200 r/min。

不同轉(zhuǎn)速下輸送能力隨水膜厚度的變化關(guān)系如圖7 所示。圖7(a)結(jié)果表明: ①同一電機(jī)轉(zhuǎn)速下,液體出口壓力隨水膜厚度增加呈拋物線型增加;水膜厚度越大,曲線增長(zhǎng)越緩慢。②增大電機(jī)轉(zhuǎn)速會(huì)使液體出口壓力大幅增加。圖7(b)結(jié)果表明: ①裝置空載時(shí)功耗低,向空腔內(nèi)注水會(huì)使功耗急劇增加;水膜厚度達(dá)到約0.4R以后,總功耗基本不再變化。②增大電機(jī)轉(zhuǎn)速會(huì)使功耗大幅增加。

圖7 不同轉(zhuǎn)速下輸送能力試驗(yàn)結(jié)果Fig.7 Experimental results of conveying capacity at different motor speed

在試驗(yàn)條件下,不同轉(zhuǎn)速下液體出口壓力和總功耗的增長(zhǎng)模式與仿真結(jié)果一致。在2200 r/min 轉(zhuǎn)速下,取一組仿真和試驗(yàn)數(shù)據(jù)對(duì)比見表2。仿真與試驗(yàn)結(jié)果比較接近,但存在偏差,主要原因?yàn)榉抡媾c試驗(yàn)環(huán)境的能量損失不同。出現(xiàn)偏差的原因有: ①仿真未考慮葉輪表面粗糙度影響,實(shí)際裝置存在摩擦損失; ②由于殼體內(nèi)部裝配間隙等存在誤差,液體出口壓力由于泄漏損失而偏低。

表2 仿真和試驗(yàn)數(shù)據(jù)對(duì)比Table 2 Comparison of simulation and experimental results

4.2 不同葉輪直徑條件下的試驗(yàn)結(jié)果

對(duì)不同葉輪直徑條件下的液體出口壓力和總功耗進(jìn)行測(cè)試,測(cè)試條件如下:葉輪直徑分別為115 mm、120 mm、125 mm 和130 mm;葉片寬度為10 mm;葉片數(shù)量為8 片;電機(jī)轉(zhuǎn)速為2200 r/min。

輸送能力隨水膜厚度變化如圖8 所示。圖8(a)結(jié)果表明: ①同一葉輪直徑下,液體出口壓力隨水膜厚度呈拋物線型增加,增長(zhǎng)模式與式(5)相同;水膜厚度越大,曲線增長(zhǎng)越緩慢。②增大葉輪直徑會(huì)使液體出口壓力增加。因此,可通過增加葉輪直徑來增大液體出口壓力,但提升效果不如增加電機(jī)轉(zhuǎn)速顯著。

圖8 不同葉輪直徑輸送能力試驗(yàn)結(jié)果Fig.8 Experimental results of conveying capacity under different impeller diameters

圖8(b)結(jié)果表明: ①裝置空載時(shí)功耗低,向空腔內(nèi)注水會(huì)使功耗急劇增加;水膜厚度達(dá)到約0.4R后,總功耗基本不再變化。②增大葉輪直徑會(huì)使功耗增加。當(dāng)葉輪直徑從115 mm 增大到130 mm 時(shí),使總功耗最大值增加100 W 以上,因此需選擇合適的葉輪直徑以避免總功耗大幅增加。

該試驗(yàn)條件下,不同葉輪直徑下液體出口壓力和總功耗的增長(zhǎng)模式與仿真結(jié)果相同。取不同葉輪直徑的最大液體出口壓力仿真和試驗(yàn)值如表3 所示。仿真與試驗(yàn)結(jié)果偏差較小,最大偏差仍小于15%,且偏差來源分析同4.1 節(jié)。

表3 不同葉輪直徑仿真和試驗(yàn)數(shù)據(jù)對(duì)比Table 3 Comparison of simulation and experimental results under different impeller diameters

4.3 不同葉片寬度條件下的試驗(yàn)結(jié)果

對(duì)不同葉片寬度條件下的液體出口壓力和總功耗進(jìn)行測(cè)試,測(cè)試條件如下:葉輪直徑為115 mm;葉片寬度分別為4 mm、6 mm、8 mm 和10 mm;葉片數(shù)量為8 片;電機(jī)轉(zhuǎn)速為2200 r/min。

輸送能力隨水膜厚度的變化關(guān)系如圖9 所示。圖9(a)結(jié)果表明: ①同一葉輪寬度下,液體出口壓力隨水膜厚度呈拋物線型增加;水膜厚度越大,曲線增長(zhǎng)越緩慢,但由于裝置注滿水時(shí)可忽略泄露損失影響,因此當(dāng)水膜厚度大于0.6R時(shí),曲線增速相較前一范圍變快; ②減小葉片寬度會(huì)使液體出口壓力下降;葉片寬度越小,減小葉片寬度液體出口壓力下降幅度越大。原因是葉片寬度越小,縫隙存在導(dǎo)致的泄露損失影響越大。因此,可通過增加葉片寬度的方式增大液體出口壓力,葉片寬度通常根據(jù)DGLSRI 腔體容積要求調(diào)整。

圖9 不同葉片寬度輸送能力試驗(yàn)結(jié)果Fig.9 Experimental results of conveying capacity under different blade widths

圖9(b)結(jié)果表明: ①裝置空載時(shí)功耗低,向空腔內(nèi)注水會(huì)使功耗急劇增加;水膜厚度達(dá)到約0.4R以后,總功耗基本不再變化; ②增大葉片寬度會(huì)使功耗增加。與葉輪直徑和電機(jī)轉(zhuǎn)速相比,葉片寬度對(duì)總功耗的影響相對(duì)較小。綜合考慮液體出口壓力和總功耗,當(dāng)裝置功耗不滿足系統(tǒng)要求時(shí),優(yōu)先考慮減小葉片寬度。

試驗(yàn)條件下,不同葉片寬度下液體出口壓力和總功耗的增長(zhǎng)模式與仿真結(jié)果相同。取不同葉片寬度的最大液體出口壓力仿真和試驗(yàn)值如表4所示。最大偏差低于13%,驗(yàn)證了仿真方法的正確性。

表4 不同葉片寬度仿真和試驗(yàn)數(shù)據(jù)對(duì)比Table 4 Comparison of simulation and experimental results under different blade widths

4.4 輸送能力提升方法

電機(jī)轉(zhuǎn)速、葉輪直徑和葉輪寬度是影響輸送能力的最主要因素,在設(shè)計(jì)過程中需重點(diǎn)考慮這3 個(gè)參數(shù)的設(shè)計(jì)。經(jīng)過分析,提出DGLSRI 輸送能力的提升方法如下:

1)確定滿足系統(tǒng)要求的電機(jī)轉(zhuǎn)速、葉輪直徑和葉片寬度。這3 種參數(shù)的增加會(huì)同時(shí)增大液體出口壓力和功耗,首先根據(jù)系統(tǒng)對(duì)液體出口壓力的要求和對(duì)功耗、體積的限制等要求,確定電機(jī)轉(zhuǎn)速和葉輪尺寸的大致范圍;通過仿真計(jì)算,確定合適的電機(jī)轉(zhuǎn)速、葉輪直徑和葉片寬度。優(yōu)先考慮增大電機(jī)轉(zhuǎn)速以提高液體出口壓力,減小寬度以降低功耗。

2)采用低轉(zhuǎn)速分離、高轉(zhuǎn)速輸送的工作模式。根據(jù)DGLSRI 的工作特點(diǎn),氣液兩相流進(jìn)入水氣分離器腔內(nèi)后,首先完成氣液分離,當(dāng)液體積累量達(dá)到閾值時(shí)開始液體輸出。氣液分離的液體積累過程往往較長(zhǎng),而最終的液體輸送過程則時(shí)間較短。因此,采用低轉(zhuǎn)速分離、高轉(zhuǎn)速輸送的工作模式可使裝置的功耗降低,該方法的工作模式如圖10 所示。采用該工作模式DGLSRI 將長(zhǎng)時(shí)間在低轉(zhuǎn)速模式下工作,可以大大降低裝置運(yùn)行的功耗。后續(xù)將持續(xù)開展實(shí)驗(yàn)研究,對(duì)該工作模式進(jìn)行進(jìn)一步驗(yàn)證和優(yōu)化。

圖10 旋片式動(dòng)態(tài)水氣分離器工作流程Fig.10 Workflow of the dynamic gas-liquid separator with rotary impeller

5 結(jié)論

本文圍繞旋片式動(dòng)態(tài)水氣分離器存在的體積大、功耗高、輸送能力不足等問題,研究了旋片式動(dòng)態(tài)水氣分離器輸送能力求解方法,提出了輸送能力提升方法,并開展了仿真及地面試驗(yàn)研究。結(jié)論如下:

1)根據(jù)DGLSRI 的工作特點(diǎn),基于CFD 技術(shù),提出了輸送能力的求解方法。方法核心是采用SST k-ω 湍流模型描述該裝置腔體內(nèi)部湍流運(yùn)動(dòng),通過VOF 多相流模型描述腔體內(nèi)部氣液兩相流動(dòng)。該方法對(duì)腔體內(nèi)氣液流體的界面追蹤具有良好效果。

2)通過仿真方法得出了不同參數(shù)對(duì)輸送能力的影響,增加電機(jī)轉(zhuǎn)速、葉輪直徑、葉片寬度均會(huì)使液體出口壓力和功耗同時(shí)增大,增加葉片數(shù)量在減小功耗的同時(shí)可提高液體出口壓力,通過試驗(yàn)測(cè)試證明了仿真方法的正確性。

3)輸送能力提升的關(guān)鍵是在滿足功耗和體積約束的條件下提升液體出口壓力。綜合考慮不同參數(shù)對(duì)該裝置液體出口壓力和功耗的影響程度,提出了輸送能力的提升方法。該方法可為在航天任務(wù)體積、功耗約束條件下進(jìn)行系統(tǒng)優(yōu)化提供技術(shù)途徑。