孫曉盟，王洪志，林國慶，曲國平，郝 明，時 黛

?

U型管式換熱器進口截面流場數(shù)值模擬

孫曉盟1，王洪志2，林國慶3*，曲國平1，郝 明4，時 黛5

（1.長春特種設(shè)備檢測研究院，吉林 長春 130011；2.中國石油管道公司沈陽龍昌管道檢測中心，遼寧 沈陽 110000 ；3.吉林化工學院機電工程學院，吉林 吉林 132022；4.北方華錦化學工業(yè)集團有限責任公司，遼寧 盤錦 124000；5. 內(nèi)蒙古民族大學機械工程學院，內(nèi)蒙古 通遼 028000）

建立了U型管式換熱器進口截面的三維穩(wěn)態(tài)流動數(shù)學模型，求得了U型管式換熱器內(nèi)部的壓力場、速度場和溫度場分布；在此基礎(chǔ)上對U型管式換熱器內(nèi)部的溫度場、壓力場和速度場進行了討論。研究結(jié)果表明，一方面，增大U型管式換熱器進口熱流體速度，可以增加U型管式換熱器的換熱量，增大出口截面的速度，增大內(nèi)部壓強，提高內(nèi)部溫度；另一方面，運行時間越長，U型管式換熱器內(nèi)部溫度越低、出口截面速度越大，總傳熱率越低，壓力損失先減小后趨于穩(wěn)定；同時，離U型管換熱器越近的外導流筒沖刷腐蝕越嚴重。

U型管換熱器；溫度場；速度場；Fluent

管殼式換熱器在泵站、煉油廠等化工生產(chǎn)中被廣泛使用，它的結(jié)構(gòu)簡單、堅固、制造較容易，處理能力大，適應(yīng)性能，操作彈性較大，尤其在高壓、高溫和大型裝置中使用更為普遍[1]。U型管換熱器是一種典型的管殼式換熱器，其管子彎成U形，由于其結(jié)構(gòu)簡單，造價比其他換熱器便宜，承壓能力強，適用于管、殼壁溫差較大或殼程介質(zhì)易結(jié)垢需要清洗，又不宜采用浮頭式和固定管板式場合。特別適用于管內(nèi)走清潔而不易結(jié)垢的高溫、高壓、腐蝕性大的物料[2]?；趂luent軟件在管殼式換熱器中具有一定的可行性[3,4]，而且U型管換熱器在fluent數(shù)值模擬方面研究較少，本文建立了U型管換熱器的模型，研究了U型管換熱器進口截面溫度場與速度場之間的關(guān)系，通過熱流體溫度與速度數(shù)值方面的不同，用來在U型管換熱器設(shè)備腐蝕方面提供理論依據(jù)。

1 數(shù)學模型建立及求解

1.1 基本假設(shè)

圖1 U型管換熱器模型

建立了U型管換熱器進口截面模型，采用前處理軟件GAMBIT進行幾何建模和網(wǎng)格劃分，采用FLUENT進行數(shù)學模型求解計算，計算結(jié)果采用TECPLOT軟件把數(shù)值模擬得到的結(jié)果進行圖像擬合。U型管換熱器通過進口熱流進料，通過與U型管中冷流換熱以達到降低溫度的目的。在GAMBIT建立的U型管換熱器如圖1所示。

U型管換熱器中的熱流體與冷流體在進口端主要呈垂直交叉流動，熱流體的流動狀況主要受到污垢、換熱管排管根數(shù)結(jié)構(gòu)、換熱管溫度等得影響，本文對U型管換熱器熱流體做如下假設(shè)：

（1）不考慮熱流體密度體積變化，將熱流體視為不可壓縮流體；

（2）非定常流動；

（3）U型換熱管中冷流體視為恒溫；

（4）不考慮冷流體流動；

（5）不考慮污垢系數(shù)影響。

1.2 控制方程

連續(xù)方程[5]：

質(zhì)量守恒方程：

能量方程：

初始條件

式中：u、u、u—、、三個方向的速度分量,m/s;

—垂直方向的體積力，即重力，考慮是不可壓縮流體， 所以為常數(shù);

—時間,s。

1.3 網(wǎng)格劃分與計算方法

運用GAMBIT對U型管換熱器進口截面進行建模和網(wǎng)格劃分，網(wǎng)格劃分方法為六面體網(wǎng)格。U型管換熱器殼體直徑為1 300 mm，進口直徑為650 mm，進口頂部與殼體距離為150 mm，U型管換熱管直徑為20 mm，管間距為32 mm，U形換熱管根數(shù)為564根，外倒流筒數(shù)共41個，由于對稱關(guān)系，換熱管數(shù)和外倒流筒數(shù)各取一半。截面寬700 mm。U型管材質(zhì)為10#鋼，熱流體為HF，密度為0.881 4 kg/m3，熱導率為0.045 4 W/(m·K)，動力粘度為1.72×10-5 kg/(m·s)，溫度為313K；冷流體為水，溫度為307 K；殼體材質(zhì)為10#鋼，密度為7.86×103 kg/m3，熱導率為44.2 W/(m·K)，比熱容477 J/(kg·K)；環(huán)境壓強為101 325 Pa，重力加速度為9.8 m/s2，環(huán)境溫度為288 K。劃分網(wǎng)格如圖2所示，網(wǎng)格節(jié)點數(shù)量約為2.41×106個，采用SIMPLE算法。連續(xù)方程和能量方程收斂判據(jù)殘差分別為10-3和10-6。在計算機配置處理器為Inter(R) Core(TM)2 Duo CPU 2.20 GHz、主板Lenovo和內(nèi)存2 GB。

圖2 U型管換熱器的三維網(wǎng)格模型

1.4 邊界條件

2 模擬結(jié)果與分析

2.1 速度場分布

從圖3(a) U型管換熱器進口截面速度矢量圖中可以看出，由于受到重力的影響，HF在流經(jīng)最上部外導流筒分流流量最大并且流量隨著高度的降低依次遞減，達到U型管換熱器最下端的外導流筒分流出的HF流量越少。由于HF的腐蝕因素與其濃度、溫度以及HF流量有關(guān)[6, 7]，因此，U型管換熱器上部外導流筒進口附近所受腐蝕與其他區(qū)域相比較為嚴重，受損的幾率大大增加。

圖3(b)為通過計算后得到的U型管換熱器截面出口處的速度平均值與時間的曲線。由于出口截面的直徑增大，與進口速度相比，出口截面的速度要小很多。隨著時間的增加，出口截面速度隨之增大。同時，進口初始速度越大，出口截面的速度也越大。

2.2 溫度場分布

從圖4中可以看出，從U型管換熱器[8]上部到下部，溫度呈現(xiàn)逐漸升高的趨勢并且有一定的溫度梯度，這是因為外導流筒位置分布的不同而造成的結(jié)果。隨著速度的增大，低溫區(qū)域逐漸減小，高溫區(qū)域不斷擴大。這說明U型管換熱器進口速度越大，U型管表面溫度越高。

圖4 U型管換熱器運行3 h對稱面溫度分布

從圖5中分析，U型管換熱器殼程上部中心區(qū)域溫度最低，溫度約為309 K。這是因為HF在U型管換熱器分流過程中流經(jīng)U型管換熱器殼程上部中心區(qū)域的流量較少，這減少了HF與U型管的直接接觸，換熱效果較其他區(qū)域較差。因此，U型管換熱器上部中心區(qū)域的溫度最低。隨著U型管換熱器進口速度的增大，U型管換熱器截面出口等溫線向著U型管換熱器殼程上部中心移動，309 K的區(qū)域逐漸減小，311.5 K的區(qū)域增大較為明顯，這說明熱流體HF的流量增加，會導致U型管換熱器內(nèi)部溫度的上升。

2.3 壓力場分布

從圖6中可以看出，U型管換熱器進口速度的增加，U型管換熱器內(nèi)部的壓強越大。在U型管換熱器內(nèi)部壓強區(qū)域從下到上逐漸降低，這是因為在U型管換熱器對稱面區(qū)域，由于U型管換熱器對稱面殼體下部存在外導流筒。從外導流筒分流出的HF量對U型管換熱器中心對稱區(qū)域起著增大U型管換熱器內(nèi)部壓強的作用。由于重力的作用，HF流向U型管換熱器上部區(qū)域的流量逐漸減少，U型管換熱器內(nèi)部壓強也隨之降低。

圖5 U型管換熱器運行3 h出口截面溫度分布

圖6 U型管換熱器運行3h對稱面壓力分布

從圖7中分析，主要壓強分布與圖6一致，但是在U型管換熱器殼程上部存在著一定的負壓區(qū)域，其進口速度越大，該區(qū)域的壓強越小。

圖7 U型管換熱器運行3 h截面出口壓強分布

2.4 流場關(guān)系

從圖8(a) U型管換熱器壓力損失隨時間的關(guān)系曲線中分析可知，U型管換熱器進口速度越大，壓力損失也就越大。U型管換熱器剛啟動不久，壓力損失最大。在運行到0.5~1.5 h區(qū)間內(nèi)，換熱器內(nèi)部壓力損失急劇降低。在運行到1.5 h之后，壓力損失趨于穩(wěn)定。

圖8 U型管換熱器壓力損失、出口總傳熱率及平均溫度隨時間的關(guān)系曲線

從圖8(b) U型管換熱器出口總傳熱率隨時間的關(guān)系曲線中可知，U型管換熱器進口速度越大，總傳熱量越大。隨著U型管換熱器運行時間的增加，總傳熱量隨之降低，但降低的幅度不大。從圖8(c) U型管換熱器出口平均溫度隨時間的關(guān)系曲線可知，U型管換熱器進口速度越小，出口截面的溫度也就越小。在運行1~1.5 h和2.5~3 h區(qū)間內(nèi)，溫度降低較快。在運行3 h之內(nèi)，溫度降幅維持在0.6 K以內(nèi)。

3 結(jié)論

（1）對U型管換熱器進行CFD數(shù)值模擬，提高U型管換熱器進口速度，出口截面速度越大，內(nèi)部壓強越大，U型管換熱器內(nèi)部溫度將上升，對換熱器的換熱效果起到一定的積極作用。

（2）U型管換熱器進口速度越大，運行時間越長，總傳熱率越低，出口截面溫度越低，出口截面速度越大，壓力損失先減小后趨于穩(wěn)定。

（3）外導流筒離U型管換熱器進口處越近，所受的沖刷腐蝕越嚴重。

［1］王志魁, 化工原理(第三版)[M]. 北京: 化學工業(yè)出版社， 2004-10: 187.

［2］王爭昪, 王學生, 戰(zhàn)洪仁, 等. 淺談U型管換熱器[J]. 廣東化工, 2009, 36(10):172-188.

［3］王艷云, 李志安, 劉紅禹, 等. FLUENT 軟件對管殼式換熱器殼程流體數(shù)值模擬方法可行性的驗證[J]. 管道技術(shù)與設(shè)備, 2007(6): 46-48.

［4］簡棄非,肖愷.基于FLUENT軟件的強化傳熱管特性三維數(shù)值模擬研究[J].低溫與特氣,2006,24 (6): 6-9.

［5］朱紅鈞,林元華,謝龍漢. FLUENT流體分析及仿真實用教程[M].北京:人民郵電出版社,2010:13-17．

［6］中國腐蝕與防護學會《金屬腐蝕手冊》編輯委員會．金屬腐蝕手冊[M]．上海：上?？茖W技術(shù)出版社，1987．

［7］王洪志, 胡傳順, 梁平,等. 氫氟酸濃度對Ni-Cu-P鍍層耐蝕性的影響[J]. 石油化工高等學校學報, 2011, 24(4):84-87．

［8］林國慶,王茂廷,王飛等.肋片式換熱器有限元熱應(yīng)力分析[J].中國安全生產(chǎn)科學技術(shù),2012,8(7):148-151．

Numerical Simulation of U- tube Heat Exchanger Inlet Section Flow Field

1,2,3,1,4,5

（1.Changchun Special Equipment Inspection and Research Institute, Jilin Changchun 130011, China; 2.Shenyang Longchang pipeline inspection center, Liaoning Shenyang 110000,China; 3. School of Mechanical and Electrical Engineering ,Jllin Institute of Chemical Technology,Jilin Jilin132022, China; 4.North Huajin Chemical Industries Group Corporation, Liaoning Panjin 124000, China; 5. College of Mechanical Engineering, Inner Mongolia University for the nationalities, Inner Mongolia Tongliao 028000, China）

The three-dimensional steady flow mathematical model about inlet section of U tube type heat exchanger was established, and then the pressure field distribution, temperature field distribution and velocity field distribution of the U type heat exchanger were obtained by the model.Based on this, the pressure field, temperature field and velocity field of the U type heat exchanger were discussed. The results show that, increasing inlet velocity of hot fluid in the U type heat exchanger, for one thing ,can increase heat exchange; for another, the longer the running time, the lower the internal temperature of the heat exchanger,the higher the outlet velocity, the lower the total heat transfer rate, the pressure loss first decreases and then tends to stability; corrosion erosion effect in all regions of U type heat exchanger is enhanced, especially the area between the U type heat exchanger inlet and the casing wall.

U-tube heat exchange；Temperature field；Velocity field；Fluent

TQ 051.7

A

1671-0460（2014）06-1117-04

2013-11-02

孫曉盟（1985-），男，吉林長春人，助理工程師，2009年畢業(yè)于吉林化工學院過程裝備與控制工程專業(yè)，研究方向為特種設(shè)備安全評定技術(shù)：從事特種設(shè)備檢驗檢測工作。E-mail：15943090361@139.com。

林國慶（1986-），男，助教，碩士研究生學歷，研究方向為過程裝備安全評定技術(shù)及能源工程。E-mail：lgq0726@126.com。