楊 鵬 陳忠海
(河北建筑工程學(xué)院,河北 張家口 075000)
管殼式換熱器作為重要的換熱設(shè)備,因結(jié)構(gòu)簡(jiǎn)單、運(yùn)行可靠和易于維護(hù)等特點(diǎn)而廣泛應(yīng)用于石油化工、冶金電力、能源動(dòng)力等領(lǐng)域,最高可達(dá)到整個(gè)換熱器市場(chǎng)的70%[1].
管殼式換熱器殼側(cè)流動(dòng)結(jié)構(gòu)是影響換熱器的關(guān)鍵因素,也是管殼式換熱器的研究重點(diǎn),以傳統(tǒng)的弓形折流板換熱器為代表的管殼式換熱器,得到了長(zhǎng)足的發(fā)展,在弓形折流板的基礎(chǔ)上進(jìn)行了許多的改進(jìn)和優(yōu)化,并提出了很多的管側(cè)支撐結(jié)構(gòu),如雙弓形折流板,曲面折流板,螺旋折流板,折流桿,整圓形折流板等.眾多學(xué)者在這個(gè)領(lǐng)域進(jìn)行了大量的理論分析、數(shù)值模擬和實(shí)驗(yàn)研究.華媛[2]提出一種新型弓形折流板結(jié)構(gòu)(扇形夾角折流板),通過改變板間距和扇面夾角比較換熱和流動(dòng)情況.喬智晶[3]發(fā)現(xiàn)在弓形折流板換熱器殼側(cè)工質(zhì)呈“Z”形流動(dòng),在缺口處易形成湍流,在折流板背面形成類似三角形傳熱死區(qū).劉琪[4]對(duì)大小孔折流板與給定波紋管組合而成的換熱器進(jìn)行了流體流動(dòng)與傳熱研究.
基于上述背景提出種錐形孔折流板管殼式換熱器,這種換熱器的特點(diǎn)是:(1)能有效的支撐管束從而避免管束流體發(fā)生誘導(dǎo)振動(dòng);(2)折流板小孔類似于錐形形噴嘴可使流體產(chǎn)生射流,一定流速的流體沖刷管束外壁,從而增加流體湍流度,減小管壁液體邊界層,有效強(qiáng)化了傳熱.本文利用Gambit建立了錐形孔折流板管殼式換熱器的物理模型,使用Fluent對(duì)其進(jìn)行數(shù)值模擬,研究其傳熱性能以及流動(dòng)性能,并將之與傳統(tǒng)的弓形折流板換熱器進(jìn)行比較.
數(shù)值模擬中湍流模型使用k-ε標(biāo)準(zhǔn)模型[5],標(biāo)準(zhǔn)k-ε方程是在動(dòng)量方程中引入了兩個(gè)新變量k(湍流動(dòng)能)與ε(湍流耗散率),其中這兩個(gè)變量的方程為:
(1)
(2)
表1 錐形孔折流板換熱器結(jié)構(gòu)參數(shù)
圖1 錐形孔折流板換熱器 圖2 弓形折流板換熱器
圖3 錐形孔折流板內(nèi)部結(jié)構(gòu)圖 圖4 殼側(cè)流場(chǎng)側(cè)面示意圖
圖5 錐形孔三維示意圖 圖6 錐形孔平面示意圖
表2 邊界條件
換熱器管程介質(zhì)為水,殼程介質(zhì)為空氣,簡(jiǎn)化處理[6]:殼體壁面采用不可滲透,無滑移絕熱邊界,殼側(cè)流體穩(wěn)態(tài)不可壓縮,不考慮折流板與管道與殼體的間隙.
表3 換熱器管殼程流體物性參數(shù)
利用Gambit軟件對(duì)三維管殼式換熱器進(jìn)行網(wǎng)格劃分,網(wǎng)格數(shù)為2018206個(gè)單元,壓力速度耦合使用SIMPLE算法、動(dòng)量、能量、湍流動(dòng)能和湍耗散率離散采用二階迎風(fēng)格式[7].
圖7和圖8分別是弓形折流板和錐形孔折流板換熱器縱截面上的速度矢量圖,如圖7所示,弓形折流板換熱器殼側(cè)流體呈Z字型流動(dòng),且沿Z字型流體流速較其他位置高.在靠近折流板拐角處有流動(dòng)死區(qū).如圖8所示錐形孔折流板換熱器殼側(cè)流體為縱向流體且經(jīng)過錐形孔時(shí)流速增大,而后逐漸減小.殼側(cè)流體充分沖刷管束,且無明顯流動(dòng)死區(qū).
圖7 弓形折流板換熱器縱截面速度矢量圖 圖8 錐形孔折流板換熱器縱截面速度矢量圖
圖9和圖10分別為弓形折流板換熱器和錐形孔折流板換熱器縱截面上的溫度云圖.如圖10所示,錐形孔折流板換熱器殼側(cè)流體溫度縱向分布均勻,而圖7弓形折流板換熱器殼側(cè)流體,縱向溫度分布變化大,且有換熱死區(qū),直觀上換熱效果較弓形折流板換熱器更好.
圖9 弓形折流板換熱器縱截面溫度云圖 圖10 錐形孔換熱器縱截面溫度云圖
圖11和圖12分別為弓形折流板換熱器和錐形孔折流板換熱器的壓力云圖,如圖11所示,弓形折流板換熱器經(jīng)過折流板時(shí),由于存在換熱死區(qū),拐角處出現(xiàn)低壓區(qū).如圖12所示,殼側(cè)流體入口壓降開始較弓形折流板大,在經(jīng)過錐形孔時(shí)壓降較大,并沿著流體方向逐漸減小.
圖11 弓形換熱器縱截面壓力云圖 圖12 錐形孔折流板換熱器縱截面壓力云圖
換熱系數(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為換熱管長(zhǎng)度.
表4 不同換熱器的殼程平均換熱系數(shù)
1)弓形折流板換熱器殼側(cè)流體在經(jīng)過折流板時(shí),由于受到折流板的阻撓,其流場(chǎng)分布呈Z字形,在折流板一側(cè)產(chǎn)生流動(dòng)死區(qū),而錐形孔折流板換熱器的流場(chǎng)分布程射流狀,射流狀態(tài)下的流體流速較大,且經(jīng)過錐形孔時(shí),殼側(cè)流體噴射沖刷管束的面積較大.
2)弓形折流板換熱器由于殼側(cè)流體對(duì)管束橫向沖擊較多,管束橫向振動(dòng)較大,而錐形孔折流換熱器殼測(cè)流體沿管束縱向流動(dòng),管束承受壓降小,對(duì)管束振動(dòng)較少.
3)綜合分析及計(jì)算得出錐形孔折流板換熱器換熱系數(shù)略大于傳統(tǒng)弓形折流板換熱器換熱系數(shù),換熱效果優(yōu)于弓形折流板換熱器.