劉展,孫培杰,李鵬,厲彥忠,晉永華

(1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;2.上海宇航系統(tǒng)工程研究所,201108,上海)

?

升空過程中低溫液氧貯箱壓力變化及熱分層研究

劉展1,孫培杰2,李鵬2,厲彥忠1,晉永華1

(1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;2.上海宇航系統(tǒng)工程研究所,201108,上海)

針對火箭升空過程,通過編寫用戶自定義程序詳細(xì)考慮了氣動(dòng)熱以及空間輻射熱的影響,數(shù)值研究了低溫液氧箱體壓力變化及流體熱分層現(xiàn)象。在計(jì)算過程中,著重考慮了大氣物性參數(shù)隨高度的變化、飛行過程中加速度變化與氣液界面相變對箱體壓力以及箱內(nèi)流體溫度分布的影響。模擬結(jié)果表明:氣動(dòng)熱對箱體控壓頻率產(chǎn)生了較大的影響。隨著氣動(dòng)熱流的增加,箱體增壓時(shí)間變短,降壓時(shí)間變長。在飛行120 s時(shí),氣動(dòng)熱流達(dá)到最大,其對箱體壓增性能的影響也最為突出,此時(shí)箱體增壓時(shí)間最短為4 s,箱體降壓時(shí)間最長,約13 s。在無排液階段,箱體壓力呈現(xiàn)波動(dòng)變化,氣液界面處氣枕被冷凝。在該過程中,液相質(zhì)量增加了11.05 kg,氣相質(zhì)量減少了1.52 kg。在增壓排液階段,盡管增壓氣體持續(xù)注入箱內(nèi),箱體壓力仍逐漸減小,而氣相質(zhì)量則呈波動(dòng)增加。隨著時(shí)間的延長,氣液相溫度均向溫度升高的方向推進(jìn)。由于空間輻射漏熱造成了排液溫度的升高,給發(fā)動(dòng)機(jī)運(yùn)行帶來安全隱患,應(yīng)采取有效絕熱措施來減少空間漏熱。

升空過程;氣動(dòng)熱;空間輻射;壓力變化;熱分層

低溫推進(jìn)劑箱體的增壓過程及流體分層現(xiàn)象對流體的順利排出以及安全貯存有著極其重要的影響。在升空過程中,低溫液體火箭將受到劇烈的氣動(dòng)熱以及各種空間輻射熱流的影響。當(dāng)考慮這兩部分熱量對低溫箱體壓增及分層影響時(shí),計(jì)算將變得異常復(fù)雜,因此有必要對其進(jìn)行深入研究。

有關(guān)低溫箱體增壓排液過程以及熱分層現(xiàn)象,研究人員開展了不同工況計(jì)算。張超采用零維模型對低溫箱體的增壓過程進(jìn)行了計(jì)算分析[1]。文獻(xiàn)[2]在計(jì)算箱體增壓過程時(shí),考慮了氣液相間質(zhì)量傳遞的影響,經(jīng)與實(shí)驗(yàn)結(jié)果對比發(fā)現(xiàn),計(jì)算模型預(yù)測結(jié)果較好。陳陽提出了一種用于計(jì)算低溫推進(jìn)劑貯箱增壓系統(tǒng)參數(shù)分布預(yù)測的數(shù)值方法,并與實(shí)驗(yàn)結(jié)果對比,驗(yàn)證了仿真方法的可靠性[3]。文獻(xiàn)[4]采用CFD技術(shù)預(yù)測了低溫箱體的自增壓過程。文獻(xiàn)[5]對低溫箱體內(nèi)部溫度分層及壓增過程進(jìn)行了編程分析。文獻(xiàn)[6]考慮氣動(dòng)熱的影響,數(shù)值模擬了增壓排液過程中低溫箱體內(nèi)部的溫度場分布。

綜上可知,研究人員對箱體的壓增過程、溫度場分布以及增壓排液等過程分別進(jìn)行了研究。考慮到實(shí)際過程中,低溫推進(jìn)劑貯箱往往是先加氣增壓,然后再注氣排液,目前針對這一實(shí)際過程所開展的研究較少。本文基于某一實(shí)際低溫火箭推進(jìn)劑箱體,在外部氣動(dòng)熱以及空間輻射下,針對箱體加氣增壓及增壓排液過程,對箱內(nèi)壓力變化及流體熱分層現(xiàn)象進(jìn)行了分析研究。本研究可為低溫推進(jìn)劑貯箱壓增系統(tǒng)以及絕熱措施等優(yōu)化設(shè)計(jì)提供一定的技術(shù)參考。

1 研究對象

選取某一實(shí)際低溫液氧箱體作為研究對象。該箱體由筒段以及上下橢球形封頭組成,筒段內(nèi)徑2 250 mm,高度1 887 mm,封頭高度800 mm,貯箱金屬層壁厚2.4 mm。金屬層外包裹有絕熱材料,筒段和下封頭外包裹厚度為20 mm絕熱層,導(dǎo)熱系數(shù)為0.03 W/(m·K),密度取40 kg/m3,比熱容為1 470 J/(kg·K)。上封頭外部包裹變厚度絕熱材料,厚度從與筒段連接處的20 mm增加到封頭頂部80 mm,密度為75 kg/m3,比熱容為1 kJ/(kg·K),導(dǎo)熱系數(shù)與溫度的關(guān)系滿足λ=0.013 2+1.437×10-7T+1.587×10-7T2。貯箱上封頭開有直徑為200 mm的增壓口,下封頭開有直徑為120 mm的排液口,排液流率為38.744 kg/s。低溫箱體初始壓力為0.435 MPa,箱體壓力控制上下限分別為0.405 MPa、0.470 MPa。增壓氣體采用360 K高溫氣氧。

對于本文所研究工況,從發(fā)射開始,低溫推進(jìn)劑箱體將依次經(jīng)歷在大氣層內(nèi)的主動(dòng)增壓模式以及脫離大氣層后入軌段的增壓排液模式。增壓排液結(jié)束后,箭體正式在軌,具體如圖1所示。

圖1 飛行過程及增壓模式

上升過程中,箭體飛行速度逐漸增大,其加速度也逐漸增大,箭體飛行最大加速度可達(dá)37 m/s2。在上升段,增壓氣體以0.18 kg/s的質(zhì)量流率間歇性變化,而在182～487 s時(shí)間內(nèi),低溫箱體開始增壓排液,增壓氣體質(zhì)量流率始終為0.18 kg/s。

2 氣動(dòng)熱及空間輻射熱模型

2.1 氣動(dòng)熱模型

圖2反映了低溫火箭上升過程中大氣壓力、密度以及溫度隨飛行高度的變化[6-7]。隨著高度的增加,空氣變得越來越稀薄,大氣壓力及密度均隨著高度的增加而減小,具體如圖2a所示。圖2b展示了大氣溫度隨高度的變化。由于大氣壓力、溫度的變化將引起大氣各物性參數(shù)的變化,所以在計(jì)算氣動(dòng)熱時(shí)需著重考慮各物性參數(shù)的影響。

(a)大氣壓力及大氣密度

(b)環(huán)境溫度圖2 上升過程中大氣參數(shù)變化

氣動(dòng)熱計(jì)算主要采用基于邊界層理論的參考焓法[6-7],相關(guān)熱力學(xué)性質(zhì)與輸送特性如下。

Eckert參考焓h*與參考溫度T*分別為

h*=he+0.5(hw-he)+0.22(hr-he)

(1)

T*=Te+0.5(Tw-Te)+0.22(Tr-Te)

(2)

(3)

式中:Tr為恢復(fù)溫度;Te為環(huán)境溫度;當(dāng)流動(dòng)為層流時(shí),n=1/2,當(dāng)流動(dòng)為湍流時(shí),n=1/3;ue為空氣外掠速度;cp為空氣比定壓熱容;下標(biāo)e表示邊界層外緣;上標(biāo)*表示參考值。

氣動(dòng)換熱系數(shù)αx與雷諾數(shù)Rex有直接關(guān)系,當(dāng)0

(4)

當(dāng)Rex>105時(shí)

(5)

式中:hr、hw為恢復(fù)溫度以及壁溫對應(yīng)的焓值;ρ*、μ*為參考溫度對應(yīng)的密度和黏度。

氣動(dòng)加熱熱流為

qx=αx(Tr-Tw)

(6)

黏度采用Sutherland公式計(jì)算

(7)

貯箱外壁與環(huán)境的輻射換熱為

(8)

式中:ε為貯箱外壁發(fā)射率;Tw為箱體壁面溫度。

2.2 空間輻射熱模型[8-9]

當(dāng)?shù)蜏鼗鸺┘壝撾x大氣層進(jìn)入空間后,以太陽輻射、地球反照輻射、地球紅外輻射以及黑背景輻射為主的各空間輻射熱流開始占據(jù)主導(dǎo)。

2.2.1 太陽輻射 太陽光被認(rèn)為是均勻的平行光束,其輻射強(qiáng)度稱為太陽常數(shù)(S=1 353 W/m2),箱體外表面所收到的太陽輻射熱量為

q1=φ1SAt

(9)

φ1=cosβs

(10)

式中:φ1為太陽輻射角系數(shù);At為箱體在太陽光線方向上的投影面積。

2.2.2 地球反照輻射 假定地球?yàn)橐宦瓷潴w,對太陽輻射的反射遵守蘭貝特定律,反射光譜與太陽光譜相同,反照率以平均反照率ρr表示,取ρr=0.3,則地球表面對箱體外表面的反照輻射熱量為

q2=φ2ρrSAt

(11)

式中:φ2為地球反照角系數(shù)(詳見文獻(xiàn)[9])。

2.2.3 地球紅外輻射 假設(shè)地球是一個(gè)均勻輻射的熱平衡體,其表面上任一點(diǎn)紅外輻射強(qiáng)度相同。箱體外表面接收到整個(gè)地球表面的紅外輻射熱為

q3=φ3SAt(1-ρr)/4

(12)

φ3=φ2cosψ

(13)

式中:φ3為地球紅外角系數(shù);ψ為相角。

低溫箱體的總輻射熱量為

qsr=q1+q2+q3

(14)

低溫箱體向空間的輻射換熱同式(8)。

低溫箱體表面吸收外部熱量的同時(shí),也向外部輻射熱量。對于升空過程,無論是氣動(dòng)熱還是空間輻射熱,都可以看做投入輻射,則箱體表面的凈吸收熱量可表示為

(15)

環(huán)境溫度Te在上升段按照圖2b所示的規(guī)律變化,進(jìn)入空間后,黑背景溫度取為4 K,氣動(dòng)熱流qx以及輻射熱流qsr都可通過相關(guān)公式計(jì)算獲得。因此,當(dāng)貯箱表面吸收率αa、貯箱表面反射率ε給定時(shí),凈熱流可表示為貯箱外壁面穩(wěn)定溫度Tw的函數(shù)(當(dāng)投入輻射為氣動(dòng)熱時(shí),αa取值1;當(dāng)投入輻射為輻射熱時(shí),αa取值0.1～0.9)。為了準(zhǔn)確獲得該非穩(wěn)態(tài)過程各參數(shù)變化,本文根據(jù)氣動(dòng)熱及空間輻射熱公式編寫相應(yīng)的UDF,作為貯箱外壁面熱流邊界條件加載到數(shù)值模型中。

3 計(jì)算設(shè)置

3.1 模型處理與邊界設(shè)置

采用Gambit 2.4.6對所研究液氧箱體劃分二維軸對稱面網(wǎng)格,計(jì)算網(wǎng)格數(shù)取為43 000。相關(guān)的計(jì)算設(shè)置詳見文獻(xiàn)[10]。與之不同的是,增壓入口以及排液口均設(shè)為質(zhì)量入口邊界,流體進(jìn)出入箱體通過質(zhì)量流速的正負(fù)加以設(shè)置。在數(shù)值計(jì)算過程中,詳細(xì)考慮了氣液相間的質(zhì)量傳遞,相變模型詳見文獻(xiàn)[6,10]。

3.2 模型驗(yàn)證

根據(jù)之前的研究[6]可知,在低溫貯箱增壓排液過程中,采用低雷諾數(shù)k-ε模型可較好地預(yù)測箱體壓力以及流體溫度變化,計(jì)算誤差均控制在10%以內(nèi)。因此,對于本文相似研究工況,此處也采用低雷諾數(shù)k-ε模型來精確計(jì)算緊貼壁面處的流動(dòng)換熱過程。

4 結(jié)果分析

4.1 無排液階段

由于箱體初始壓力為0.435 MPa,箱內(nèi)液氧具有較大的初始過冷度,再加之氣氧過熱,因此氣相將向液相傳遞熱量,并被液相冷卻。液氧箱體壓力在前182 s內(nèi)隨時(shí)間的變化如圖3所示。在液相的冷卻下,箱體壓力從初始的0.435 MPa逐漸降低,當(dāng)降低到箱體設(shè)定的壓力下限0.405 MPa時(shí),高溫氣氧以0.18 kg/s的質(zhì)量流率注入箱體,箱體壓力得以升高。當(dāng)氣枕壓力增加到箱體所設(shè)定的壓力上限0.47 MPa時(shí),增壓氣體停止注入。相比于高溫高壓氣枕,液相仍處于大過冷度狀態(tài),因此氣液界面處仍發(fā)生氣相冷凝,氣枕壓力逐漸降低。在無排液增壓過程中,氣枕壓力將在所設(shè)定的箱體壓力上下限內(nèi)波動(dòng)變化。另外,受外部氣動(dòng)漏熱的影響,箱體壓力波動(dòng)頻率出現(xiàn)不同。上升過程中氣動(dòng)熱流的變化如圖4所示,隨著時(shí)間的增加,箱體所受氣動(dòng)熱流在約120 s時(shí)達(dá)到最大值。在飛行前72 s內(nèi),由于氣動(dòng)熱流較小,箱體壓增及降壓時(shí)間所受影響較小,分別維持在5 s與9 s。當(dāng)112 s之后,氣動(dòng)熱對箱體控壓頻率產(chǎn)生了較大影響,此時(shí)氣動(dòng)熱流逐漸增加,箱體壓增時(shí)間小于5 s,并且降壓時(shí)間大于9 s。在120 s時(shí),箱體壓增時(shí)間達(dá)到最小值4 s,降壓時(shí)間達(dá)到最大值13 s。120 s之后,氣動(dòng)漏熱減少,箱體的控壓頻率又趨于初始頻率。過高或過低的增壓降壓頻率都會(huì)對箱內(nèi)熱力過程造成影響,因此為保證低溫箱體安全運(yùn)行,控壓頻率應(yīng)在合理范圍內(nèi)。

圖3 箱體壓力在前182 s內(nèi)的變化

圖4 上升過程中的氣動(dòng)熱流變化

圖5給出了上升段氣液相質(zhì)量的變化。在該過程中,氣相被液相冷卻,氣液界面處氣相冷凝,氣相質(zhì)量mv降低,液相質(zhì)量ml增加。當(dāng)氣枕壓力降低到控壓下限時(shí),增壓氣體開始注入,氣枕質(zhì)量隨之增加。盡管高溫氣氧的注入增加了氣相的質(zhì)量,但高溫高壓的氣枕仍被過冷的液氧冷卻,因此氣相一直處于冷凝狀態(tài)。液相在冷卻氣相的同時(shí),其質(zhì)量一直增加,氣相質(zhì)量則隨著增壓氣體的間歇性注入而出現(xiàn)波動(dòng)變化。至于氣相質(zhì)量每次變化的波動(dòng)幅度不同,也主要與間歇性注入的增壓氣體有關(guān)。在整個(gè)過程中,液相質(zhì)量增加了11.05 kg,而氣相質(zhì)量變化約為1.52 kg。

圖5 無排液階段氣液相質(zhì)量變化

圖6給出了無排液增壓過程中不同時(shí)刻箱體內(nèi)部溫度場及氣液相分布。為對比觀察流體溫度分層,不同時(shí)刻溫度分布均控制在88～200 K。在外部氣動(dòng)熱作用下,緊貼壁面的流體被加熱,產(chǎn)生自然對流。在10 s時(shí),箱體內(nèi)部液相區(qū)自然對流及溫度分布均相對穩(wěn)定。受增壓氣體的影響,與氣相接觸的部分液相產(chǎn)生了渦旋,擾動(dòng)較大,溫度分布不規(guī)律。隨著時(shí)間的延長,增壓氣體影響的液相區(qū)域逐漸增大,所帶來的液相擾動(dòng)也越來越大,沿箱體高度方向的氣相向液相的熱量滲透也逐漸增加。受熱浮力的影響,高溫氣氧注入箱體后向上微翹,到達(dá)箱體壁面后再返回,在氣流的上下部形成了方向相反的渦旋,造成該區(qū)域較大的擾動(dòng)。受氣動(dòng)漏熱的影響,箱體內(nèi)部自然對流在112～180 s十分顯著。即便如此,箱體界面處氣相仍處于冷凝狀態(tài)。圖6的左半部分展示了在該過程中氣液相的分布。由于該過程沒有液相排出,不同時(shí)刻氣液相分布大致相同,但受增壓氣體擾動(dòng)的影響,界面處仍存在不明顯的波動(dòng)。

圖6 無排液階段不同時(shí)刻下流體熱分層及相分布圖

(a)氣液相中心線溫度變化

(b)液相中心線溫度變化圖7 無排液階段中心線溫度變化

圖7展示了不同時(shí)刻貯箱中心線溫度分布,可以看出,不同時(shí)刻氣液相中線溫度分布基本一致。由于氣相區(qū)中線溫度只監(jiān)測到增壓口處,該處流體溫度并沒有達(dá)到增壓氣氧的溫度,所以圖7a中所展示氣相中線溫度最高也僅在127 K左右。在10 s及60 s時(shí),氣相中線頂部溫度先升高后降低,這主要是由增壓氣體回流擾動(dòng)所致。至于在120 s時(shí),氣枕最高溫度高于180 s時(shí)的氣枕最高溫度,則主要與增壓氣體的注入有關(guān)。由于在120 s時(shí)沒有增壓氣體的注入,高溫氣體向上運(yùn)動(dòng),并在頂部累積,形成了較大的溫度梯度。在180 s時(shí),隨著高溫氣氧注入,氣相擾動(dòng)較大,溫度分布不均,向液相的傳熱也增加,最終導(dǎo)致該時(shí)刻氣枕最高溫度低于120 s時(shí)的最大值,但180 s時(shí)的氣枕區(qū)(-0.41～-0.36 m區(qū)域)主體溫度仍高于120 s時(shí)的主體溫度。圖7b展示了不同時(shí)刻液相中線溫度分布,可以看出,不同時(shí)刻液相中線溫度分布基本一致,均隨高度的增加而減小,并在-0.3～1.5 m區(qū)產(chǎn)生較大溫度梯度。至于圖中液相中線溫度隨時(shí)間的增加,并沒有呈現(xiàn)規(guī)律的遞增變化,這仍與增壓氣體的注入情況以及所帶來的液相擾動(dòng)有直接關(guān)系。

4.2 排液階段

低溫火箭末級飛離大氣層后,液氧箱體將經(jīng)歷注氣排液過程。圖8展示了排液過程中箱體壓力隨時(shí)間的變化。當(dāng)排液口打開時(shí),在高溫氣氧的活塞效應(yīng)推動(dòng)下,液氧以設(shè)定的質(zhì)量流率38.744 kg/s排出箱體。盡管在整個(gè)過程中,高溫氣氧持續(xù)注入,但由于液體的排出,氣相所占區(qū)域增大,氣枕膨脹,致使箱體壓力降低。箱體壓力在180～210 s內(nèi)出現(xiàn)了陡降,這主要與箱體形狀結(jié)構(gòu)有關(guān)。由于初始液位處于箱體上封頭處,隨著液體的排出,液位將逐漸下降到筒段。這一過程出現(xiàn)氣相體積的突然增大,最終導(dǎo)致箱體壓力的陡降。當(dāng)液位處在箱體筒段處,箱體壓力以較小的斜率近似線性下降,而當(dāng)液位處在下封頭時(shí),如在450 s以后,箱體壓力下降斜率又有所增加。在該過程中,箱體壓力由0.428 MPa降低到0.272 MPa。在整個(gè)排液過程中,增壓氣氧持續(xù)以0.18 kg/s的質(zhì)量流率注入箱體,氣氧消耗量線性增加,最終消耗量為54.54 kg。

圖8 排液過程中箱體壓力變化

圖9給出了排液過程中氣液相質(zhì)量變化曲線,隨著增壓氣氧的注入以及液氧的排出,液相的質(zhì)量近似線性下降,氣相的質(zhì)量則波動(dòng)增加。因?yàn)橐合嗳蕴幱谶^冷狀態(tài),氣相被過冷液相冷凝。由于冷凝相變量相對于液相總質(zhì)量來說特別小,因此在一定的排液速率下,液相質(zhì)量是線性降低的。至于氣相,在排液剛開始處于冷凝狀態(tài),其質(zhì)量降低;隨著氣氧注入,其質(zhì)量又逐漸增加。在冷凝與注氣的共同作用下,氣相質(zhì)量呈現(xiàn)出波動(dòng)增加的趨勢。

圖9 排液過程中氣液相質(zhì)量變化

圖10展示了排液階段箱內(nèi)流體溫度分層及氣液相分布,從中容易看出,隨著液體的排出,液位逐漸降低,氣相所占空間逐漸增大。對比不同時(shí)刻溫度分布可以發(fā)現(xiàn),排液初期液相溫度場擾動(dòng)較大。這是因?yàn)榇藭r(shí)液位距離增壓口較近,受增壓氣體影響較大,與氣相所接觸的液相區(qū)域產(chǎn)生了較大的擾動(dòng)。在240～420 s期間,液位處于筒段,在增壓氣體的活塞效應(yīng)下,流體液位下降速率較大,掩蓋了自然對流的影響,所以在遠(yuǎn)離增壓口處的氣枕區(qū)以及液相區(qū)溫度分布都比較整齊。在480 s時(shí),液位處在下封頭,增壓氣體對液相的影響變?nèi)?但排液口對剩余液體產(chǎn)生了擾動(dòng),形成了小范圍的渦流。對于氣相區(qū),在增壓口處形成了360 K小范圍的高溫區(qū)域,并在該區(qū)域產(chǎn)生擾動(dòng)。

圖10 排液階段不同時(shí)刻下流體熱分層及相分布圖

圖11 排液階段貯箱中心線溫度的變化

在低溫箱體脫離大氣層到入軌階段,實(shí)際漏入箱體的空間輻射熱流在0.6 W/m2左右。盡管空間漏熱對箱體的壓降以及氣液相冷凝過程沒產(chǎn)生明顯的影響,但受外部漏熱的影響,箱內(nèi)低溫流體溫升速率明顯增加。從圖11所展示的流體中線溫度變化可以看出,隨著時(shí)間的持續(xù),液位下降,氣液相溫度逐漸向增大的方向推進(jìn),并在界面間形成了良好的過渡。再者,在氣相區(qū)向液相的傳熱以及外部漏熱的影響下,箱內(nèi)剩余液相的溫度逐漸升高。圖12展示了排液流體溫度隨時(shí)間的變化。在整個(gè)排液過程中,隨著空間漏熱的增加,排液流體溫度由最初的88.9 K增加到91.26 K。當(dāng)進(jìn)入發(fā)動(dòng)機(jī)的低溫流體溫度要求低于91 K時(shí),在450 s后,低溫流體將不能排出,否則流體在通過泵的時(shí)候很容易造成汽蝕,并帶來安全隱患,因此需通過有效措施減小空間漏熱,以降低低溫流體的溫升。

圖12 排液流體溫度隨時(shí)間的變化

5 結(jié) 論

(1)在前182 s,低溫液氧箱體受到強(qiáng)烈的氣動(dòng)熱影響。隨著氣動(dòng)熱流的增加,箱體增壓時(shí)間逐漸縮短,降壓時(shí)間逐漸延長。當(dāng)氣動(dòng)熱在120 s達(dá)到最大值時(shí),箱體增壓時(shí)間最短為4 s,而降壓時(shí)間卻延長到13 s。由于箱內(nèi)液相過冷度較大,氣液界面處出現(xiàn)氣枕冷凝。該階段箱體壓力呈現(xiàn)波動(dòng)變化,液相質(zhì)量增加了11.05 kg,氣相質(zhì)量變化約為1.52 kg。當(dāng)增壓工況改變時(shí),與氣枕接觸的液相區(qū)域熱分層受增壓氣體影響較大。

(2)在182～487 s,隨著液相的排出,氣相體積增大,盡管增壓氣體持續(xù)注入,箱體壓力仍逐漸降低。在氣液界面處,氣枕仍被冷凝,考慮增壓氣體的注入,氣枕質(zhì)量波動(dòng)增加。受空間輻射漏熱的影響,氣相溫度隨著液位的降低向溫度升高的方向推進(jìn),并且剩余液體的溫度也逐漸升高。由于空間漏熱升高了流體的溫度,為了保證發(fā)動(dòng)機(jī)入口溫度低于設(shè)定值,需采取有效措施來減少空間漏熱。

[1] 張超, 魯雪生, 田麗亭. 火箭低溫液體推進(jìn)劑增壓系統(tǒng)數(shù)學(xué)模型 [J]. 低溫與超導(dǎo), 2005, 33(2): 35-38. ZHANG Chao, LU Xuesheng, TIAN Liting. The numerical model of rocket cryogenic liquid propellant pressurization system [J]. Cryogenic and Superconductivity, 2005, 33(2): 35-38.

[2] ZILLIAC G, KARABEYOGLU M A. Modeling of propellant tank pressurization [C]∥Proceedings of 41st AIAA/ASME/SAE/ASEE Joint Propulsion Con-

ference & Exhibit. Reston, VA, USA: AIAA, 2005: 1-25.

[3] 陳陽, 張振鵬, 楊思鋒, 等. 低溫推進(jìn)劑貯箱增壓系統(tǒng)分布參數(shù)數(shù)值仿真: Ⅱ 增壓系統(tǒng)數(shù)值模型與仿真結(jié)果 [J]. 航空動(dòng)力學(xué)報(bào), 2008, 23(2): 329-335. CHEN Yang, ZHANG Zhenpeng, YANG Sifeng, et al. Distributed parameter numerical simulation for the transients of cryogenic propellant tank pressurization system: Ⅱ Numerical model and simulation results of pressurization system [J]. Journal of Aerospace Power, 2008, 23(2): 329-335.

[4] 劉展, 厲彥忠, 王磊, 等. 在軌運(yùn)行低溫液氫箱體蒸發(fā)量計(jì)算與增壓過程研究 [J]. 西安交通大學(xué)學(xué)報(bào), 2015, 49(2): 135-140. LIU Zhan, LI Yanzhong, WANG Lei, et al. Evaporation calculation and pressurization process of on-orbit cryogenic liquid hydrogen storage tank [J]. Journal of Xi’an Jiaotong University, 2015, 49(2): 135-140.

[5] LIU Z, WANG L, JIN Y, et al. Development of thermal stratification in a rotating cryogenic liquid hydrogen tank [J]. International Journal of Hydrogen Energy, 2015, 40(43): 15067-15077.

[6] WANG L, LI Y, ZHAO Z, et al. Transient thermal and pressurization performance of LO2 tank during helium pressurization combined with outside aerodynamic heating [J]. International Journal of Heat & Mass Transfer, 2013, 62(1): 263-271.

[7] LIU Z, LI Y, XIE F, et al. Thermal performance of foam/MLI for cryogenic liquid hydrogen tank during the ascent and on orbit period [J]. Applied Thermal Engineering, 2015, 98: 430-439.

[8] 劉展, 孫培杰, 李鵬, 等. 微重力下低溫液氧貯箱熱分層研究 [J]. 低溫工程, 2016(1): 25-31. LIU Zhan, SUN Peijie, LI Peng, et al. Research on thermal stratification of cryogenic liquid oxygen tank in microgravity [J]. Cryogenics, 2016(1): 25-31.

[9] 閔桂榮. 衛(wèi)星熱控制技術(shù) [M]. 北京: 中國宇航出版社, 1991: 5-40.

[10]劉展, 孫培杰, 李鵬, 等. 地面停放低溫液氧貯箱熱物理過程研究 [J]. 西安交通大學(xué)學(xué)報(bào), 2016, 50(9): 36-42. LIU Zhan, SUN Peijie, LI Peng, et al. Research on thermal physical process of cryogenic liquid oxygen tank during the ground parking [J]. Journal of Xi’an Jiaotong University, 2016, 50(9): 36-42.

(編輯 荊樹蓉)

Research on Pressure Change and Thermal Stratification of Cryogenic Liquid Oxygen Tank during the Ascent Process

LIU Zhan1,SUN Peijie2,LI Peng2,LI Yanzhong1,JIN Yonghua1

(1. School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China; 2. Shanghai Institute of Aerospace System Engineering, Shanghai 201108, China)

A numerical simulation method was adopted to study the pressure change and fluid thermal stratification in liquid oxygen tank during the ascent process by a user-defined function (UDF) considering the aerodynamic heat and space radiations. The influences of atmosphere physical parameters, flight acceleration and phase change occurring in the liquid-vapor interface on the tank pressure and fluid temperature distribution were taken into account in the present numerical model. The calculation results showed that aerodynamic heat has great influence on the tank pressurization frequency and with the increase of aerodynamic heat flux, the tank pressurization time increases and depressurization time reduces gradually. After 120 s flight, the aerodynamic heat flux reaches its maximum value, and its influence on the tank pressurization behavior becomes most prominent. At this moment, the minimum tank pressure rise time and the maximum tank pressure reduction time are 4 s and 13 s, respectively. During the pressurization process without outflow, the tank pressure fluctuates within the set pressure range, and the phase change is condensation at the interface. The liquid mass increases about 11.05 kg while the ullage mass decreases about 1.52 kg in this process. When it is in the outflow period, the tank pressure decreases with the continuous injection of pressurized gas, and the ullage mass shows a fluctuating increase. Both the gas and liquid temperatures boost to higher temperatures with time. As the space radiation leakage leads to the temperature increase of drainage liquid, and hence resulting in safety issues for the rocket engine, some effective measures should be taken to reduce the space leakage.

ascent process; aerodynamic heat; space radiation; tank pressure change; thermal stratification

2016-01-27。 作者簡介:劉展(1988—),男,博士生;厲彥忠(通信作者),男,教授。 基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(51376142);航天低溫推進(jìn)劑技術(shù)國家重點(diǎn)實(shí)驗(yàn)室開放課題(SKLTSCP1505);上海航天核攀項(xiàng)目(ZY2015-015)。

時(shí)間:2016-07-15

10.7652/xjtuxb201611015

V511

A

0253-987X(2016)11-0097-07

網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20160715.1659.002.html