溫度場(chǎng)的三維數(shù)值分析

2013-06-20  by:廣州有限元分析、培訓(xùn)中心-1CAE.COM  來(lái)源:仿真在線

對(duì)皮托管處于高馬赫數(shù)均勻氟里昂氣體來(lái)流下不同頭部形狀、不同姿態(tài)的皮托管進(jìn)行流固耦合計(jì)算,給出了固體域溫度場(chǎng)三維數(shù)值分析結(jié)果。計(jì)算是采用ANSYS軟件的結(jié)構(gòu)模塊進(jìn)行建模,然后應(yīng)用ICEM軟件進(jìn)行網(wǎng)格的劃分,采用CFX軟件進(jìn)行數(shù)值求解和進(jìn)行數(shù)據(jù)的后處理。

核工業(yè)理化工程研究院 丁保庚 王穎 周陳龍 來(lái)源:e-works
關(guān)鍵字:皮托管 溫度場(chǎng) 三維數(shù)值分析 CFX ansys icem

1  前言

在用皮托管測(cè)量高速氣體流速時(shí),皮托管的頭部形狀和姿態(tài)對(duì)其溫度有很大的影響,過(guò)高的溫度要影響皮托管的使用,使用測(cè)試方法測(cè)量其溫度是困難的,特別是對(duì)于小型皮托管情況,所以非常有必要用數(shù)值分析的方法研究皮托管的頭部姿態(tài)對(duì)其頭部溫度的影響。

2  數(shù)學(xué)模型

    進(jìn)行計(jì)算時(shí)作了以下幾點(diǎn)假設(shè):

    a.  皮托管頭部的來(lái)流是均勻的;

    b.  忽略由于熱輻射引起的熱交換;

    c.  流體的粘性系數(shù)、熱傳導(dǎo)系數(shù)和等壓比熱系數(shù)都為常數(shù),不隨溫度變化。

由于研究的是流固耦合問(wèn)題,因此需要求解的控制方程包括固體域的溫度方程和流體域的控制方程。固體域應(yīng)用固體熱傳導(dǎo)方程,流體域的控制方程采用流體力學(xué)的N-S方程,流、固界面兩側(cè)熱通量相等。     

3  計(jì)算模型

采用ANSYS軟件^[1]的結(jié)構(gòu)模塊建立計(jì)算模型,由于幾何模型相對(duì)于中分面是對(duì)稱的,因此只需要建立一半模型。皮托管為直徑5mm,長(zhǎng)50mm的直管,頭部存在70度切角,處在半徑25mm的柱形流體域當(dāng)中。計(jì)算模型如下所示(見圖1)。

應(yīng)用ICEM CFD網(wǎng)格預(yù)處理軟件^[2]進(jìn)行模型網(wǎng)格的劃分,采用塊體結(jié)構(gòu)拓?fù)鋭澐志W(wǎng)格的方法進(jìn)行,網(wǎng)格采用八節(jié)點(diǎn)六面體網(wǎng)格,具體網(wǎng)格如下(見圖2和圖3)。

應(yīng)用CFX流體計(jì)算分析軟件^[3]設(shè)定邊界條件進(jìn)行流固耦合求解,皮托管處于來(lái)流馬赫數(shù)為5.0,來(lái)流壓力和溫度分別為600Pa和20℃的均勻氟里昂氣體來(lái)流中。皮托管的端部邊界條件為20℃。

流體計(jì)算采用湍流模型,具體采用SST(剪應(yīng)力輸運(yùn)模型)模型。

4  計(jì)算結(jié)果及分析

4.1 頭部無(wú)切角、均勻來(lái)流無(wú)攻角情況

當(dāng)皮托管頭部為平頭沒有切角來(lái)流無(wú)攻角時(shí),計(jì)算結(jié)果如下所示(見圖4~6)。

圖4 頭部無(wú)切角時(shí)皮托管的溫度分布

圖5 頭部無(wú)切角時(shí)皮托管對(duì)稱面上氣體等壓線圖

圖6 頭部無(wú)切角時(shí)皮托管對(duì)稱面上氣體等溫線圖

由計(jì)算結(jié)果可以發(fā)現(xiàn),整個(gè)皮托管呈現(xiàn)頭部溫度較高,最高溫度達(dá)到247℃,隨后沿管道方向逐漸降低的溫度分布(如圖4所示)。圖5和圖6中顯示皮托管頭部形成了非常明顯的正激波,使得皮托管頭部附近的氣體溫度急劇上升,這些溫度升高的氣體,通過(guò)熱傳導(dǎo)將熱量傳遞到皮托管的頭部,使得其頭部的溫度也急劇升高,最高溫度達(dá)342℃。激波后的高溫氣體沿流動(dòng)方向隨著熱量的散失溫度逐漸下降,因而其通過(guò)熱傳導(dǎo)傳到皮托管的熱量也沿著管道的方向逐漸降低,直至環(huán)境溫度。     

4.2  頭部有70度切角、均勻來(lái)流無(wú)攻角情況

當(dāng)流無(wú)攻角時(shí),計(jì)算結(jié)果如下所示(見圖7~9)。

圖7 來(lái)流無(wú)攻角情況皮托管的溫度分布

    圖8 來(lái)流無(wú)攻角情況皮托管對(duì)稱面上氣體等壓線圖

    圖9 來(lái)流無(wú)攻角情況皮托管對(duì)稱面上氣體等溫線圖

當(dāng)皮托管頭部被切掉20度,成70度切角形狀時(shí),皮托管仍然呈現(xiàn)頭部溫度較高,尾部溫度較低的溫度分布,如圖7所示。從流場(chǎng)的計(jì)算結(jié)果(如圖8、圖9),在皮托管頭部尖角處形成了斜激波,然后逐漸過(guò)渡為脫體激波,這是使得其頭部附近的氣體溫度以及皮托管溫度升高的主要原因。但皮托管頭部最高溫度下降到了240℃,說(shuō)明由于頭部存在切角,減弱了皮托管頭部的正激波,使得其頭部的溫度較平頭時(shí)有所下降。

4.3 頭部有70度切角、來(lái)流有5度攻角情況

當(dāng)來(lái)流存在5度的攻角時(shí),其計(jì)算結(jié)果如下所示(見圖10~12)。

圖10 來(lái)流有5度攻角時(shí)皮托管的溫度分布

圖11 來(lái)流有5度攻角皮托管對(duì)稱面上氣體等壓線圖

    圖12來(lái)流有5度攻角時(shí)皮托管對(duì)稱面上氣體等溫線圖

當(dāng)來(lái)流存在5度攻角時(shí),整個(gè)皮托管的溫度分布狀況以及頭部最高溫度與無(wú)攻角時(shí)基本相同(圖10),其流場(chǎng)結(jié)果(圖11、圖12)也相差不大,說(shuō)明來(lái)流存在小角度(5度范圍之內(nèi))攻角,對(duì)皮托管的溫度分布以及附近流場(chǎng)幾乎沒有影響。

4.4 來(lái)流有15度攻角情況

當(dāng)來(lái)流存在15度攻角時(shí),其計(jì)算結(jié)果如下所示(見圖13~15)。

圖13 來(lái)流有15度攻角時(shí)皮托管的溫度分布

圖14 來(lái)流有15度攻角皮托管對(duì)稱面上氣體等壓線圖

    圖15來(lái)流有15度攻角時(shí)皮托管對(duì)稱面上氣體等溫線圖

當(dāng)來(lái)流存在15度攻角時(shí),整個(gè)皮托管的溫度分布趨勢(shì)(圖13)以及其流場(chǎng)結(jié)果(圖14、圖15)與第二、第三種情況基本一致,說(shuō)明當(dāng)攻角增大時(shí)(15度以內(nèi)),皮托管頭部的氣體流動(dòng)狀態(tài)并沒有受到很大的影響,但皮托管頭部的溫度上升到250℃,說(shuō)明來(lái)流攻角的變化范圍增大以后,會(huì)造成皮托管頭部溫度的升高。     

5  結(jié)論

通過(guò)計(jì)算分析,獲得了以下幾點(diǎn)結(jié)論:

a)  在工作介質(zhì)為氟里昂氣體,來(lái)流馬赫數(shù)和壓力分別為5.0和600Pa的均勻來(lái)流流場(chǎng)中,皮托管頭部溫度較高,會(huì)在240℃以上;

b)  當(dāng)皮托管頭部存在切角時(shí),可以有效的減弱皮托管頭部的激波損耗,使其頭部溫度下降;

 c)  來(lái)流存在小角度(5度范圍之內(nèi))攻角時(shí),對(duì)皮托管溫度分布以及頭部附近流場(chǎng)影響不大,因此在進(jìn)行測(cè)量時(shí),對(duì)于皮托管的安裝偏差在5度之內(nèi)無(wú)明顯的影響;

d)  來(lái)流攻角增大到15度會(huì)導(dǎo)致皮托管頭部溫度的升高,但對(duì)其流場(chǎng)及溫度分布的整個(gè)趨勢(shì)沒有太大的影響。

相關(guān)標(biāo)簽搜索：溫度場(chǎng)的三維數(shù)值分析 Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識(shí) ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓(xùn)