国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

微重力環(huán)境下低溫流體的數(shù)值模擬與特性研究

2016-10-09 06:14潘雁頻冶文蓮王麗紅曲家闖蘭州空間技術(shù)物理研究所真空技術(shù)與物理重點(diǎn)實(shí)驗(yàn)室蘭州730000
真空與低溫 2016年2期
關(guān)鍵詞:表面張力液面氣液

張 鎧,潘雁頻,冶文蓮,王麗紅,曲家闖(蘭州空間技術(shù)物理研究所 真空技術(shù)與物理重點(diǎn)實(shí)驗(yàn)室,蘭州 730000)

微重力環(huán)境下低溫流體的數(shù)值模擬與特性研究

張 鎧,潘雁頻,冶文蓮,王麗紅,曲家闖
(蘭州空間技術(shù)物理研究所 真空技術(shù)與物理重點(diǎn)實(shí)驗(yàn)室,蘭州 730000)

在重力場(chǎng)下,貯箱中流體的狀態(tài)是已知的,流體狀態(tài)是按規(guī)律分層分布的,而在微重力環(huán)境下,貯箱中流體狀態(tài)是未知的,氣相和液相的分布是不規(guī)則的,其界面的幾何形狀也是不確定的。通過采用VOF方法,對(duì)貯箱中的液氫在微重力條件下進(jìn)行了數(shù)值模擬研究,主要考慮了不同重力加速度,接觸角和不同充液量等因素的影響,從而得到在微重力條件下影響氣液分布的主要因素。

微重力;液氫;重力加速度;接觸角;充液量

0 引言

為了促進(jìn)人類向更寬更廣的深空領(lǐng)域推進(jìn)和發(fā)展,實(shí)現(xiàn)人類在月球的居住以及火星載人探測(cè)等任務(wù)是未來航天發(fā)展的趨勢(shì),其中低溫液體在航天技術(shù)中的應(yīng)用越來越廣泛,從作為推進(jìn)劑的燃料、宇航員呼吸用的氧和氮,以及其他用途的氫、氙、甲烷等,都可以以低溫液體的形式貯存,需要在微重力環(huán)境下對(duì)低溫流體的數(shù)值與特性進(jìn)行模擬研究。

1 低溫流體貯箱仿真模擬計(jì)算

在重力場(chǎng)下,貯箱中流體的狀態(tài)是已知的,流體狀態(tài)是按規(guī)律分層分布的,而在微重力環(huán)境下,貯箱中流體狀態(tài)是未知的,氣相和液相的分布是不規(guī)則的,其界面的幾何形狀也是不確定的,兩種環(huán)境下的流體狀態(tài)如圖1所示。

圖1 不同環(huán)境下的流體狀態(tài)圖

從圖1可以看出,重力場(chǎng)下的流體狀態(tài)是按規(guī)律分層分布的,而在微重力下的流體狀態(tài)是不規(guī)則的,對(duì)幾種不同工況下的低溫貯箱的氣液界面進(jìn)行了數(shù)值模擬,分析了在微重力條件下影響貯箱內(nèi)氣液分布的主要因素。

1.1 VOF模型

CFD(Computational Fluid Dynamics)[1]軟件,基于有限體積法,通過計(jì)算機(jī)數(shù)值計(jì)算和圖像顯示,對(duì)包含有流體流動(dòng)和熱傳導(dǎo)等相關(guān)物理現(xiàn)象的系統(tǒng)做出分析??梢钥醋魇窃诹鲃?dòng)基本方程(質(zhì)量守恒方程、動(dòng)量守恒方程、能量守恒方程)的控制下對(duì)流動(dòng)的數(shù)值模擬。通過分析,可以得到極其復(fù)雜問題的流場(chǎng)內(nèi)各個(gè)位置上的基本物理量(如速度、壓力、溫度、濃度等)的分布,以及這些物理量隨時(shí)間的變化情況。CFD技術(shù)對(duì)低溫氣液兩相流研究具有很大的可行性

模擬氣液兩相流動(dòng)與傳熱的數(shù)值方法主要分為顆粒軌道模型、雙流體模型和運(yùn)動(dòng)界面追蹤技術(shù),其中運(yùn)動(dòng)界面追蹤技術(shù)能準(zhǔn)確、全面地反映氣液兩相間的關(guān)系。目前,F(xiàn)LUENT是廣泛使用的數(shù)值模擬軟件,采用流體體積(VOF)方法追蹤相間的運(yùn)動(dòng)界面,能較準(zhǔn)確的模擬相間的蒸發(fā)和冷凝現(xiàn)象。

VOF模型[6]依靠的是兩種或多種流體(或相)沒有互相穿插的事實(shí)。對(duì)增加到模型里的每一附加相,就引進(jìn)一個(gè)變量,即計(jì)算單元里相的容積比率。在每個(gè)控制容積內(nèi),所有相的volume fraction的和為1。所有變量及其屬性的區(qū)域被各相共享并且代表了容積平均值,只要每一相的容積比率在每一位置是可知的。在任何給定單元內(nèi)的變量及其屬性或者純粹代表了一相或者代表了相的混合,這取決于容積比率值。在單元中,如果第q相流體的容積比率記為αq,那么就有三個(gè)條件是可能的:當(dāng)αq=0時(shí),第q相流體在單元中是空的;當(dāng)αq=1時(shí),第q相流體在單元中是充滿的;當(dāng)0<αq<1時(shí),單元中包含了第q相流體和一相或者其他多相流體的界面。對(duì)于氣液兩相有αl+αg=1;其中αl為液相體積,αg為氣相體積。

1.2 計(jì)算模型

由于計(jì)算的復(fù)雜性,采用二維模型。模型包括一段柱段和兩個(gè)橢圓形的封頭,其中貯箱內(nèi)徑為1.5 m,圓柱體部分高為1.5 m,箱頂部以及底部橢圓形封頭高均為0.4 m。計(jì)算中選取液氫和氫氣分別作為液相與氣相的模擬工質(zhì),取環(huán)境溫度為-259℃,此時(shí)表面張力系數(shù)為2.99×10-3N/m,初始填充液體量為箱體體積的一半,初始?xì)庖航缑嬖O(shè)為水平界面,然后對(duì)于不同的填充液體量,不同的重力加速度環(huán)境,不同的壁面浸潤(rùn)性條件即取不同的接觸角下非定常運(yùn)動(dòng)的情況進(jìn)行數(shù)值模擬,得到最終穩(wěn)定的液面形狀。

1.3 控制方程

由于流體流速很低,可以當(dāng)不可壓縮流體考慮,并且采用VOF方法求解流體的運(yùn)動(dòng)和自由界面的變形,因此描述其流體運(yùn)動(dòng)和流體體積函數(shù)的控制方程可以用方程式(1)~(3)表示:

上述三個(gè)方程分別表示為流體的連續(xù)性方程、動(dòng)量方程和流體體積函數(shù)的守恒方程。其中:v為流體的運(yùn)動(dòng)速度;g為重力加速度;ρ為流體的密度;p為流體的壓力;Fσ為表面張力;τ為粘性應(yīng)力張量,對(duì)于牛頓流體有[2-3]:τ=2μS。

應(yīng)變率張量S由式(4)給出:

對(duì)于上述控制方程,還需要加上相應(yīng)的邊界條件和初始條件,才能使問題封閉。流體運(yùn)動(dòng)在各個(gè)壁面均設(shè)置為固壁面條件,采用了無穿透、無滑移的邊界條件,即用二階迎風(fēng)格式對(duì)方程進(jìn)行離散,速度和壓力耦合采用SIMPLE方法。

2 計(jì)算結(jié)果與分析

2.1 微重力條件下氣液界面隨時(shí)間變化情況

在微重力10-5g時(shí),接觸角10°,液氫填充率在50%情況下,初始?xì)庖航缑嬖O(shè)為水平界面,對(duì)貯箱內(nèi)的氣液分布隨時(shí)間變化進(jìn)行數(shù)值模擬,如圖2所示。

從圖2數(shù)值模擬可以看出,在微重力環(huán)境下,壁面有良好的浸潤(rùn)性時(shí),液體在表面張力的作用下,沿著壁面爬升,傾向于形成大球形氣泡。

2.2 不同重力加速度對(duì)液面變化的影響

在接觸角10°,液氫填充率在50%情況下,初始?xì)庖航缑嬖O(shè)為水平界面,通過改變重力加速度對(duì)貯箱內(nèi)的氣液分布進(jìn)行數(shù)值模擬,得到穩(wěn)定的氣液界面如圖3所示。

從圖3模擬可以看出,隨著重力加速度的增加,即Bo數(shù)的增加,重力對(duì)液面變化的影響逐漸提高。在重力加速度為g時(shí),此時(shí)重力對(duì)液面變化起主要作用,表面張力因素可以忽略,液面形狀為地球上常規(guī)的氣液界面即水平面的方向發(fā)展。在重力加速度很小,即處于微重力環(huán)境下,此時(shí)表面張力起主要作用,氣液界面與地球上常規(guī)的氣液界面有很大的差異,當(dāng)液體與壁面呈浸潤(rùn)性時(shí),液體會(huì)沿著貯箱的壁面爬升,直到勢(shì)能達(dá)到最小。

不同重力加速度對(duì)液面變化可以由無量綱數(shù)Bond數(shù)表示,液面狀態(tài)主要是受到表面張力和重力的影響,無量綱數(shù)Bond數(shù)即表示表面張力和重力的相對(duì)影響,其表達(dá)如式(5):

式中:σ為表面張力;d為特征尺寸;?ρ為兩種流體相的密度差;g為重力加速度。

2.3 不同接觸角對(duì)液面變化的影響

接觸角即在水平放置的均勻平整的固體表面上滴一滴液滴,當(dāng)液滴達(dá)到如圖4所示的穩(wěn)定形態(tài)[4]。以固、液、氣三相相交點(diǎn)作液氣界面的切線,與固液界面的夾角,用θ表示。

圖2 液面隨時(shí)間變化過程圖

圖3 不同重力加速度穩(wěn)態(tài)液面情況圖

圖4 液滴在固體表面的形態(tài)圖

當(dāng)固液氣三相接觸并處于平衡狀態(tài)時(shí),接觸角滿足Yong氏方程:

式中:σsv為氣固表面張力系數(shù);σsl為液固表面張力系數(shù);σlv為氣液表面張力系數(shù);當(dāng)溫度和氣體壓力不變時(shí),σsv、σsl、σlv均為常數(shù)。在這種條件下,式(6)等價(jià)于θ=const。

如果0<θ<90°,表示附著力大于內(nèi)聚力,液體浸潤(rùn)固體;

如果90°<θ<180°,則附著力小于內(nèi)聚力,呈現(xiàn)非浸潤(rùn)狀態(tài)。

在微重力10-5g下,液氫填充率在50%情況下,初始?xì)庖航缑嬖O(shè)為水平界面,改變液體相對(duì)固壁的浸潤(rùn)性,即改變接觸角的大小,對(duì)貯箱內(nèi)的氣液分布進(jìn)行數(shù)值模擬,結(jié)果如圖5所示。

由數(shù)值模擬結(jié)果可以看出,在微重力環(huán)境下,液體的接觸角對(duì)于氣液界面的影響占主導(dǎo)作用。隨著接觸角的增大,即液體相對(duì)于固壁浸潤(rùn)性的減弱,液體沿固壁向上爬升的能力降低。且在接觸角大于90°時(shí),液體不會(huì)向上爬升,反而會(huì)向下凹陷,使得整個(gè)液面呈向下凹陷的錐形。

2.4 不同液氫填充率對(duì)液面變化的影響

在微重力環(huán)境下,接觸角取10°,改變液體填充率,通過數(shù)值模擬液面從初始時(shí)刻的平面自由釋放到趨于穩(wěn)定的過程,最終趨于穩(wěn)定后的氣液界面如圖6所示。

由數(shù)值模擬結(jié)果可以看出,在微重力環(huán)境下,接觸角取10°,不同液體的填充率對(duì)于氣液界面的影響作用不大,都是由于表面張力的作用,液體沿著壁面爬升,傾向于形成大球形氣泡。

圖6 不同液氫填充率穩(wěn)態(tài)液氫液面情況圖

3 結(jié)論

通過采用VOF方法數(shù)值模擬在微重力條件下貯箱內(nèi)氣液界面的變化情況,對(duì)影響“氣枕”形變的接觸角、重力加速度、表面張力等特性參數(shù)進(jìn)行了分析,獲得了其對(duì)氣液界面變化的影響規(guī)律:在微重力環(huán)境下,壁面有良好的浸潤(rùn)性時(shí),液體在表面張力的作用下,沿著壁面爬升,傾向于形成大球形氣泡;在重力加速度為g時(shí),此時(shí)重力對(duì)液面變化起主要作用,表面張力因素可以忽略,液面形狀為地球上常規(guī)的氣液界面即水平面的方向發(fā)展;在重力加速度很小,處于微重力環(huán)境下,此時(shí)表面張力起主要作用,微重力環(huán)境下,液體的接觸角對(duì)于氣液界面的影響占主導(dǎo)作用,并且隨著接觸角的增大,液體相對(duì)于固壁浸潤(rùn)性的減弱,液體沿固壁向上爬升的能力降低;不同液體的填充率對(duì)于氣液界面的影響作用不大,都是由于表面張力的作用,液體沿著壁面爬升,傾向于形成大球形氣泡。這些規(guī)律為下一步根據(jù)已確定氣液界面和氣枕位置來分析貯箱內(nèi)溫度場(chǎng)和壓力場(chǎng)變化起到指導(dǎo)作用。

[1]丁義鋒,汪洋,王小軍,等.微尺度高效換熱器的熱流計(jì)算與設(shè)計(jì)[J].真空與低溫,2015,21(2):99-102.

[2]劉趙淼,趙婷婷,申峰.重力和接觸角對(duì)表面張力貯箱內(nèi)液體流動(dòng)的影響[J].力學(xué)學(xué)報(bào),2015,47(3):430-440.

[3]李章國(guó),劉秋生,紀(jì)巖,等.航天器貯箱氣液自由界面追蹤數(shù)值模擬[J].空間科學(xué)學(xué)報(bào),2008,28(1):69-73.

[4]張世舉,程延海,邢方方,等.接觸角與表面自由能的研究現(xiàn)狀與展望[J].煤礦機(jī)械,2011,32(10):8-10.

[5]Salzman J A.Low-gravity reorientation in a scale-model Centaurliquid-hydrogentank[R].NASA,1981.

[6]IrtCW,nicholsBD.Volumeoffluid(VOF)methodfordynamics of free boundary[J].Journal of Computational Physics,1981,39(1):201-225.

NUMERICAL SIMULATION AND CHARACTERISTIC RESEARCH OF CRYOGENIC FLUID UNDER MICROGRAVITY ENVIRONMENY

ZHANG Kai,PAN Yan-ping,YE Wen-lian,WANG Li-hong,QU Jia-chuang
(Science and Technology on Vacuum Technology and Physics Laboratory,Lanzhou Institute of Physics,Lanzhou 730000,China)

In the gravity field,the state of fluid in the tank is known,which is based on the law of hierarchical distribution.While under the circumstance of microgravity,the state of fluid in the tank is unknown.And the distribution of gaseous phase and liquid phase is irregular.So the geometrical shape of the interface is unsure.This paper uses VOF method to simulate the liquid hydrogen storage tank under the circumstance of microgravity,which mainly considers the influence of different factors,such as gravitational acceleration,contact angle and different working fluid filling.Thus the numerical simulation can summarizes the principal factor on the gas-liquid interface in microgravity.

microgravity;liquid hydrogen;gravitational acceleration;contact angle;working fluid filling

TB61

A

1006-7086(2016)02-0085-05

10.3969/j.issn.1006-7086.2016.02.005

2015-12-25

張鎧(1990-),男,安徽省淮北市人,碩士研究生,主要從事制冷低溫方面的研究。E-mail:zhang117607@126.com。

猜你喜歡
表面張力液面氣液
雙輥薄帶連鑄結(jié)晶輥面對(duì)液面波動(dòng)的影響
運(yùn)載火箭氣液組合連接器動(dòng)態(tài)自動(dòng)對(duì)接技術(shù)
微重力下兩相控溫型儲(chǔ)液器內(nèi)氣液界面仿真分析
氣液分離罐液位計(jì)接管泄漏分析
吸管“喝”水的秘密
白金板法和白金環(huán)法測(cè)定橡膠膠乳表面張力的對(duì)比
神奇的表面張力
GY-JLY200數(shù)據(jù)記錄儀測(cè)試動(dòng)液面各類情況研究
基于新型C4D的小管道氣液兩相流流型辨識(shí)方法
表面張力
基隆市| 福泉市| 津南区| 乌海市| 应用必备| 南雄市| 左云县| 平度市| 东乡| 临沂市| 巫山县| 文山县| 左云县| 平顺县| 康保县| 宝山区| 治县。| 长丰县| 浙江省| 肇源县| 张北县| 铜梁县| 福鼎市| 龙井市| 拜泉县| 兴仁县| 成武县| 启东市| 吴忠市| 永新县| 卓尼县| 尤溪县| 易门县| 宁国市| 三门峡市| 乐都县| 敦化市| 抚远县| 和顺县| 六盘水市| 从江县|