新型錐形孔折流板管殼式換熱器數(shù)值模擬研究

楊 鵬 陳忠海

(河北建筑工程學(xué)院，河北 張家口 075000)

0 引 言

管殼式換熱器作為重要的換熱設(shè)備，因結(jié)構(gòu)簡單、運(yùn)行可靠和易于維護(hù)等特點(diǎn)而廣泛應(yīng)用于石油化工、冶金電力、能源動力等領(lǐng)域，最高可達(dá)到整個換熱器市場的70%[1].

管殼式換熱器殼側(cè)流動結(jié)構(gòu)是影響換熱器的關(guān)鍵因素，也是管殼式換熱器的研究重點(diǎn)，以傳統(tǒng)的弓形折流板換熱器為代表的管殼式換熱器，得到了長足的發(fā)展，在弓形折流板的基礎(chǔ)上進(jìn)行了許多的改進(jìn)和優(yōu)化，并提出了很多的管側(cè)支撐結(jié)構(gòu)，如雙弓形折流板，曲面折流板，螺旋折流板，折流桿，整圓形折流板等.眾多學(xué)者在這個領(lǐng)域進(jìn)行了大量的理論分析、數(shù)值模擬和實(shí)驗(yàn)研究.華媛[2]提出一種新型弓形折流板結(jié)構(gòu)(扇形夾角折流板)，通過改變板間距和扇面夾角比較換熱和流動情況.喬智晶[3]發(fā)現(xiàn)在弓形折流板換熱器殼側(cè)工質(zhì)呈“Z”形流動，在缺口處易形成湍流，在折流板背面形成類似三角形傳熱死區(qū).劉琪[4]對大小孔折流板與給定波紋管組合而成的換熱器進(jìn)行了流體流動與傳熱研究.

基于上述背景提出種錐形孔折流板管殼式換熱器,這種換熱器的特點(diǎn)是：(1)能有效的支撐管束從而避免管束流體發(fā)生誘導(dǎo)振動；(2)折流板小孔類似于錐形形噴嘴可使流體產(chǎn)生射流，一定流速的流體沖刷管束外壁，從而增加流體湍流度，減小管壁液體邊界層，有效強(qiáng)化了傳熱.本文利用Gambit建立了錐形孔折流板管殼式換熱器的物理模型，使用Fluent對其進(jìn)行數(shù)值模擬，研究其傳熱性能以及流動性能，并將之與傳統(tǒng)的弓形折流板換熱器進(jìn)行比較.

1 數(shù)值模擬理論

1.1 湍流模型

數(shù)值模擬中湍流模型使用k-ε標(biāo)準(zhǔn)模型[5]，標(biāo)準(zhǔn)k-ε方程是在動量方程中引入了兩個新變量k(湍流動能)與ε(湍流耗散率)，其中這兩個變量的方程為：

(1)

(2)

2 錐形孔折流板換熱器傳熱模型

2.1 物理模型

表1 錐形孔折流板換熱器結(jié)構(gòu)參數(shù)

圖1 錐形孔折流板換熱器 圖2 弓形折流板換熱器

圖3 錐形孔折流板內(nèi)部結(jié)構(gòu)圖 圖4 殼側(cè)流場側(cè)面示意圖

圖5 錐形孔三維示意圖 圖6 錐形孔平面示意圖

2.2 邊界條件的設(shè)定及其簡化

表2 邊界條件

換熱器管程介質(zhì)為水，殼程介質(zhì)為空氣，簡化處理[6]：殼體壁面采用不可滲透，無滑移絕熱邊界，殼側(cè)流體穩(wěn)態(tài)不可壓縮，不考慮折流板與管道與殼體的間隙.

表3 換熱器管殼程流體物性參數(shù)

2.3 網(wǎng)格的劃分及求解設(shè)置

利用Gambit軟件對三維管殼式換熱器進(jìn)行網(wǎng)格劃分，網(wǎng)格數(shù)為2018206個單元，壓力速度耦合使用SIMPLE算法、動量、能量、湍流動能和湍耗散率離散采用二階迎風(fēng)格式[7].

3 模擬結(jié)果及其分析

3.1 速度場分析

圖7和圖8分別是弓形折流板和錐形孔折流板換熱器縱截面上的速度矢量圖，如圖7所示，弓形折流板換熱器殼側(cè)流體呈Z字型流動，且沿Z字型流體流速較其他位置高.在靠近折流板拐角處有流動死區(qū).如圖8所示錐形孔折流板換熱器殼側(cè)流體為縱向流體且經(jīng)過錐形孔時流速增大，而后逐漸減小.殼側(cè)流體充分沖刷管束，且無明顯流動死區(qū).

圖7 弓形折流板換熱器縱截面速度矢量圖 圖8 錐形孔折流板換熱器縱截面速度矢量圖

3.2 溫度場分析

圖9和圖10分別為弓形折流板換熱器和錐形孔折流板換熱器縱截面上的溫度云圖.如圖10所示，錐形孔折流板換熱器殼側(cè)流體溫度縱向分布均勻，而圖7弓形折流板換熱器殼側(cè)流體，縱向溫度分布變化大，且有換熱死區(qū)，直觀上換熱效果較弓形折流板換熱器更好.

圖9 弓形折流板換熱器縱截面溫度云圖 圖10 錐形孔換熱器縱截面溫度云圖

3.3 壓力場分析

圖11和圖12分別為弓形折流板換熱器和錐形孔折流板換熱器的壓力云圖，如圖11所示，弓形折流板換熱器經(jīng)過折流板時，由于存在換熱死區(qū)，拐角處出現(xiàn)低壓區(qū).如圖12所示，殼側(cè)流體入口壓降開始較弓形折流板大，在經(jīng)過錐形孔時壓降較大，并沿著流體方向逐漸減小.

圖11 弓形換熱器縱截面壓力云圖 圖12 錐形孔折流板換熱器縱截面壓力云圖

4 換熱性能比較分析

換熱系數(shù)根據(jù)以下公式求出：

(3)

(4)

(5)

A=NπDL

(6)

其中Cp為比熱容[J·(kg·℃)-1]，m1、m2殼側(cè)和管側(cè)的進(jìn)出口流量,A為換熱面積，T1、T4、T3、T2分別為管殼側(cè)進(jìn)口和出口工質(zhì)的溫度，N為換熱管數(shù)量,D為換熱管直徑,L為換熱管長度.

表4 不同換熱器的殼程平均換熱系數(shù)

5 結(jié) 語

1)弓形折流板換熱器殼側(cè)流體在經(jīng)過折流板時，由于受到折流板的阻撓，其流場分布呈Z字形，在折流板一側(cè)產(chǎn)生流動死區(qū)，而錐形孔折流板換熱器的流場分布程射流狀，射流狀態(tài)下的流體流速較大，且經(jīng)過錐形孔時，殼側(cè)流體噴射沖刷管束的面積較大.

2)弓形折流板換熱器由于殼側(cè)流體對管束橫向沖擊較多，管束橫向振動較大，而錐形孔折流換熱器殼測流體沿管束縱向流動，管束承受壓降小，對管束振動較少.

3)綜合分析及計(jì)算得出錐形孔折流板換熱器換熱系數(shù)略大于傳統(tǒng)弓形折流板換熱器換熱系數(shù)，換熱效果優(yōu)于弓形折流板換熱器.