郭 土，馬貴陽，張一楠，孫 皓，曹燕龍

（遼寧石油化工大學(xué) 石油天然氣工程學(xué)院, 遼寧 撫順 113001）

換熱器在化工過程中起著至關(guān)重要的作用，投資費(fèi)用一般占全部投資費(fèi)用的 40%左右[1]。傳統(tǒng)的單弓形折流板換熱器換熱效率較低，殼程壓力損失較大，存在流動死區(qū)，容易結(jié)垢等[2]，難以滿足生產(chǎn)的要求。因此不斷有新型的殼程折流支撐結(jié)構(gòu)的換熱器出現(xiàn)，如折流桿換熱器[3,4]、螺旋折流板換熱器[5-7]、花格板換熱器[8]等。

錢才富，高宏宇[9]等提出一種新型管殼式換熱器——曲面弓形折流板換熱器，用CFD技術(shù)對曲面弓形折流板換熱器和普通弓形折流板換熱器的殼側(cè)流體流動與傳熱性能進(jìn)行了數(shù)值模擬與研究，發(fā)現(xiàn)和普通弓形折流板相比，曲面弓形折流板所引起的殼程流體速度分布在流道內(nèi)更加均勻，流線更為光滑，而且流動“死區(qū)”明顯減小。并且對曲面弓形折流板換熱器的殼側(cè)流體流動與傳熱性能進(jìn)行了實(shí)驗(yàn)驗(yàn)證，模擬和數(shù)值結(jié)果基本一致，在殼程進(jìn)口流速相等的情況下，曲面弓形折流板換熱器比普通弓形折流板換熱器的殼側(cè)壓降要降低9%～24%，殼程傳熱系數(shù)比普通弓形折流板換熱器提高3%～11%。

本文提出一種新型的管殼式換熱器——曲面弓形折流板換熱器，并且在不同進(jìn)口流速0.1～0.8 m/s（每隔0.1 m/s做一次計(jì)算）下，對曲率半徑為0.75D、1D、1.25D、1.5D（D為換熱器殼體直徑）的球面弓形折流板換熱器和普通弓形折流板換熱器的殼側(cè)流體流動與傳熱性能進(jìn)行數(shù)值模擬與研究，得到了殼程流體流場分布以及殼程壓降和傳熱系數(shù)，研究不同的球面弓形折流板曲率半徑對換熱器殼程壓降、殼程換熱系數(shù)的影響。

1 數(shù)學(xué)模型的建立

換熱器殼程流體的流動與傳熱過程需要滿足三大守恒方程。因?yàn)橥牧餍?yīng)對流動和傳熱影響非常大，所以計(jì)算模型選用標(biāo)準(zhǔn)模型， 標(biāo)準(zhǔn)的運(yùn)輸方程為：

式中：t—時(shí)間,s；

μ—流體粘度，kg/(m·s)；

ρ—流體密度，kg/m3；

Gk—均速度梯度引起的湍動能產(chǎn)生項(xiàng)；

Gb—浮力引起的湍動能產(chǎn)生項(xiàng)；

YM—可壓縮湍流脈動膨脹對總得耗散率的影響。

經(jīng)驗(yàn)常數(shù)C1ε=1.44，C3ε=1, Cμ=0.09。

2 弓形折流板三維模型的建立

2.1 幾何模型

模型分為普通弓形折流板換熱器和球面弓形折流板換熱器兩種類型，其中球面弓形折流板換熱器折流板曲率半徑分別取0.75D、1D、1.25D、1.5D，具體結(jié)構(gòu)參數(shù)見表1。0.75D球面弓形折流板換熱器結(jié)構(gòu)示意圖如圖1。換熱器殼程結(jié)構(gòu)相對復(fù)雜,流體流動和傳熱過程也十分復(fù)雜。因此在保證相關(guān)物理量準(zhǔn)確的前提下, 建立模型時(shí)作如簡化假設(shè):

表1 換熱器主要結(jié)構(gòu)參數(shù)Table 1 Main structural parameters of the heat transfer

(1)忽略折流板與換熱管、折流板與筒體間的間隙;

(2)假設(shè)流體流動和傳熱過程是均一、穩(wěn)定的;

(3)對傳熱導(dǎo)致的流體物性的變化忽略不計(jì),如密度、黏度和比熱容等;

(4)筒體外壁與外界無熱交換;

(5)流體為牛頓流體、不可壓縮、各向同性且連續(xù)。

圖1 0.75D球面弓形折流板換熱器示意圖Fig.1 The structure chart of the 0.75D spherical baffle heat exchanger

2.2 網(wǎng)格劃分及邊界條件的設(shè)置

根據(jù)表1的幾何參數(shù)分別繪制出普通弓形折流板換熱器和4種曲率半徑的球面弓形折流板換熱器的幾何體，由于弓形折流板換熱器殼程流域模型十分復(fù)雜，所以網(wǎng)格劃分采用分塊劃分的方法，間距為2.5 mm。圖2為球面弓形折流板換熱器截面網(wǎng)格劃分示意圖。

圖2 球面弓形折流板換熱器截面網(wǎng)格劃分示意圖Fig.2 The sectional mesh dividing effect of the spherical baffle heat exchanger

使用分離求解器，穩(wěn)態(tài)隱式格式求解以保證收斂的穩(wěn)定性；速度壓力耦合方式采用基于交錯網(wǎng)格的SIMPLE算法；流體為水，物性參數(shù)為等效溫度下的常數(shù)；假設(shè)入口速度均勻分布，殼體壁面和折流板采用不可滲透、無滑移絕熱邊界條件；忽略重力影響。

使用速度入口（velocity-inlet）和出流（out flow）邊界條件。入口條件的設(shè)置采用湍流強(qiáng)度和水力直徑，湍流強(qiáng)度通過（3）式得到；動量、能量以及湍流參量的求解采用二階迎風(fēng)格式。湍流模型選用標(biāo)準(zhǔn)模型。進(jìn)口流速范圍為0.1～0.8 m/s，每隔0.1 m/s做一次模擬計(jì)算。

3 數(shù)值模擬與結(jié)果分析

3.1 對換熱器壓力場的分析

殼程入口流速為0.3 m/s時(shí)，0.75D球面弓形折流板換熱器殼程壓力場中心線截面等值線圖如圖3所示：模擬結(jié)果表明，在相同入口流速下，球面弓形折流板換熱器比普通換熱器殼程壓降要小 8%～11%；且曲率半徑越小，換熱器殼程壓降越低。普通弓形折流板換熱器殼程壓降為328 Pa，曲率半徑為0.75D的球面弓形折流板換熱器壓降為294 Pa。

圖3 0.3 m/s時(shí)0.75D球面弓形折流板換熱器殼程壓力場中心線截面等值線圖Fig.3 The shell side central line sectional pressure field of the 0.75D spherical baffle heat exchanger when inlet velocity is 0.3 m/s

3.2 對換熱器溫度場的分析

不同流速下兩種換熱器的換熱系數(shù)曲線如圖4所示。在不同的殼程入口流速下，球面弓形折流板換熱器比普通換熱器換熱系數(shù)低，且曲率半徑越小，換熱器換熱系數(shù)越低。這主要是因?yàn)榍蛎婀握哿靼迩拾霃皆叫?，與殼程流體的流動曲線越吻合，這不利于造成殼程流體的擾流，在一定程度上降低了換熱系數(shù)，也正是由于這一點(diǎn)卻能大大的降低殼程壓降。

圖4 換熱器的換熱系數(shù)曲線Fig.4 Heat exchanger’s heat exchange coefficient curve

3.3 對兩種換熱器流動與傳熱綜合性能的分析

不同曲率半徑球面弓形折流板換熱器換熱性能提高效果如圖5所示。在殼程入口流速相等的條件下，球面弓形折流板換熱器曲率半徑越小換熱器殼程傳熱與流阻綜合性能指標(biāo)越大，但是增長的速度不是均勻的。如殼程進(jìn)口流速為0.8 m/s時(shí)，曲率半徑分別為0.75D、1D、1.25D、1.5D的球面弓形折流板換熱器殼程傳熱與流阻綜合性能指標(biāo)比普通弓形折流板換熱器提高的百分?jǐn)?shù)分別為9.5%、6.8%、5.5%和4.5%。隨著殼程入口流速的增大，球面弓形折流板換熱器換熱性能提高的效果越好，但是增加的幅度不大。

圖5 球面弓形折流板換熱器在不同進(jìn)口流速下綜合性能提高效果Fig.5 The increasing effect of different spherical baffle heat exchanger performance under different inlet flow rate

4 結(jié) 論

(1)在相同入口流速下，球面弓形折流板換熱器曲率半徑越小，換熱器殼程壓降越低，換熱系數(shù)也隨之降低。

(2) 在本文設(shè)計(jì)工況下，曲率半徑越小換熱器殼程傳熱與流阻綜合性能指標(biāo)pΔ?/越大，且增長的速度不是均勻的。

(3)考慮換熱器折流板加工工藝和生產(chǎn)成本問題，球面弓形折流板換熱器殼程傳熱與流阻綜合性能指標(biāo)pΔ?/比普通弓形折流板換熱器提高在10%以內(nèi)，效果并不特別理想，還有進(jìn)一步改進(jìn)提升的空間。

［1］ 孫成家，楊軍，陳保東.不同折流板換熱器的傳熱與流阻性能對比[J].節(jié)能技術(shù)，2005,23（1）:59-61.

［2］ 劉敏珊，董其伍，劉乾.折流板換熱器殼程流場數(shù)值模擬與結(jié)構(gòu)優(yōu)化[J].化工設(shè)備與管道，2006,43（2）:24-31.

［3］ 吳金星，董其伍，劉敏珊，等.折流桿換熱器殼程湍流和傳熱的數(shù)值模擬[J].高?；瘜W(xué)工程學(xué)報(bào)，2006,20（2）：213-216.

［4］ 馬雷，王英雙，楊杰，等. 變截面折流桿換熱器的流動與傳熱分析[J].工程熱物理學(xué)報(bào)，2012,33（1）：113-116.

［5］ 潘振，陳保東，商麗艷.螺旋折流板換熱器的研究與進(jìn)展[J].節(jié)能技術(shù)，2006,24（1）:81-85.

［6］ Zhang Zhengguo, Xu Tao, Fang Xiaoming. Experimental study on heat transfer enhancement of a helically baffled heat exchanger combined with three-dimensional finned tubes[J]. Applied Thermal Engineering ,2004, 24 (14/15) : 2293-2300.

［7］ Nasiruddin, Kam ran Siddiqui M H. Heat transfer augmentation in a heat exchanger tube using a baffle[J].International Journal of Heat and Fluid Flow , 2007, 28(2): 318-328.

［8］ 王英雙，張曉嶼，劉志春，等. 花格板換熱器的流動與傳熱[J].化工學(xué)報(bào)，2012,63（S1）：99-105.

［9］ 錢才富，高宏宇，孫海陽.曲面弓形折流板換熱器殼程流體流動與傳熱[J].化工學(xué)報(bào)，2011（5）：1233-1238.