何 超，瞿瑩，王華玲，蔣志鵬，吳宏亮

（1.淮陰工學(xué)院 機(jī)械與材料工程學(xué)院，江蘇 淮安 223003；2.淮陰工學(xué)院 江蘇省先進(jìn)制造技術(shù)重點(diǎn)實(shí)驗(yàn)室，江蘇 淮安 223003；3.淮安市鼎新電子有限公司，江蘇 淮安 223200）

管殼式換熱器廣泛應(yīng)用于電站、船舶、石油化工及制冷空調(diào)等領(lǐng)域的換熱設(shè)備中。而管殼式換熱器在這些領(lǐng)域工作過程中，往往存在著大量的相變換熱及兩相流動(dòng)問題，對(duì)換熱器本身的傳熱效率及使用壽命等有著至關(guān)重要的影響[1]。統(tǒng)計(jì)數(shù)字表明，在工業(yè)用換熱器中有近50%的換熱器均在殼程存在兩相流動(dòng)，諸如重沸器、凝結(jié)器和蒸發(fā)器等[2]。

管殼式換熱器殼程結(jié)構(gòu)復(fù)雜，流體在殼程的流動(dòng)也很復(fù)雜，使用實(shí)驗(yàn)方法很難得到換熱器內(nèi)部具體的流場(chǎng)和溫度場(chǎng)，且耗時(shí)和昂貴，傳統(tǒng)的理論計(jì)算基礎(chǔ)上進(jìn)行建模實(shí)驗(yàn)驗(yàn)證的設(shè)計(jì)方法顯得越來越低效。隨著計(jì)算機(jī)技術(shù)的不斷升級(jí)以及CFD技術(shù)的快速發(fā)展，采用數(shù)值模擬的方法來對(duì)換熱器進(jìn)行流體與傳熱的特性計(jì)算，預(yù)先得出傳熱與阻力的準(zhǔn)數(shù)關(guān)聯(lián)式，為實(shí)際尺寸模型的設(shè)計(jì)提供可靠的依據(jù)成為了一種趨勢(shì)，是繼實(shí)驗(yàn)方法后的一種重要換熱器研究手段[3]。因此，以常用的管殼式換熱器（弓形折流板換熱器）作為研究對(duì)象，該換熱器殼程流體發(fā)生冷凝相變及兩相流動(dòng)，采用CFD模擬軟件建立管殼式換熱器三維實(shí)體模型，對(duì)殼程流體的速度場(chǎng)、溫度場(chǎng)和冷凝水體積分布情況進(jìn)行了分析研究[4]。

1 換熱器殼程流場(chǎng)數(shù)值模擬

1.1 基本控制方程

流體介質(zhì)的流動(dòng)必然會(huì)受到物理學(xué)守恒定律的支配，流體介質(zhì)的流動(dòng)要滿足三個(gè)基本的物理量守恒定律：物質(zhì)的質(zhì)量守恒定律、物質(zhì)的動(dòng)量守恒定律和物質(zhì)的能量守恒定律。如果流體流動(dòng)屬于湍流流動(dòng)狀態(tài)，則還要滿足湍流輸送方程。殼程介質(zhì)為飽和水蒸氣，其換熱后發(fā)生冷凝，涉及傳質(zhì)和相變，因此還需要添加一個(gè)組分質(zhì)量守恒方程。

（1）質(zhì)量守恒方程（連續(xù)性方程）

對(duì)于不可壓縮流體，ρ為常數(shù)，故不可壓縮流體三維流動(dòng)的質(zhì)量守恒方程可表示為：

（2）動(dòng)量守恒方程

（3）能量守恒方程

（4）組分質(zhì)量守恒方程

式中：U為速度矢量，u、v、w分別為速度矢量U在x、y、z坐標(biāo)方向的分量，ρ為流體密度，t為時(shí)間，μ為流體的動(dòng)力黏度，p為壓力，Su、Sv、Sw三個(gè)動(dòng)量守恒方程的廣義源項(xiàng)，T為溫度，λ為流體的導(dǎo)熱系數(shù)，Cp為比熱容，ST為粘性耗散項(xiàng)，Cs為組分s的體積濃度，Ds為擴(kuò)散系數(shù)，Ss為系統(tǒng)內(nèi)s組分質(zhì)量的產(chǎn)生率[5]。

1.2 計(jì)算模型

進(jìn)行數(shù)值模擬的換熱器為普通管殼式換熱器，換熱器結(jié)構(gòu)為單管程、單殼程和單弓形折流板，其結(jié)構(gòu)尺寸見表1。

表1 換熱器結(jié)構(gòu)尺寸

1.3 計(jì)算模型與計(jì)算條件的簡(jiǎn)化

由于換熱器殼程結(jié)構(gòu)復(fù)雜，在保證相關(guān)物理量準(zhǔn)確的前提下，建立模型時(shí)作如下簡(jiǎn)化假設(shè)：（1）折流板與換熱管、折流板與殼體之間的間隙忽略不計(jì)；（2）拉桿、定距管、不凝氣出口管等結(jié)構(gòu)忽略；（3）流體流動(dòng)和傳熱過程穩(wěn)定，不考慮污垢熱阻；（4）流體視為不可壓縮流，各向同性且連續(xù)；（5）氣相看作純組分，忽略液膜熱阻[6]。

計(jì)算邊界條件的簡(jiǎn)化：（1）殼體光滑，無滑移邊界條件；（2）假設(shè)換熱管外壁溫度為恒溫；（3）中心面為對(duì)稱邊界；（4）速度進(jìn)口和壓力出口邊界條件，出口表壓力為0 MPa。

1.4 幾何模型建立及網(wǎng)格劃分

該軟件以有限元法為基礎(chǔ)，采用迭代方法進(jìn)行數(shù)值模擬計(jì)算，需要對(duì)模型進(jìn)行網(wǎng)格劃分。建立的模型為三維模型，網(wǎng)格劃分時(shí)運(yùn)用對(duì)稱特性，實(shí)際計(jì)算時(shí)只取一半的模型進(jìn)行，從而減小了模型尺寸，降低了計(jì)算成本。簡(jiǎn)化后的幾何模型如圖1，劃分完的網(wǎng)格的模型如圖2。

圖1 簡(jiǎn)化幾何模型

圖2 網(wǎng)格劃分

1.5 求解設(shè)置

1.5.1 物理場(chǎng)選擇

換熱器殼程和管程流體均為湍流流動(dòng)，同時(shí)熱流體和冷流體之間存在傳熱情況。對(duì)于上述情形，該軟件的流體流動(dòng)模塊下設(shè)置了非等溫流動(dòng)，再下設(shè)湍流k-ε物理場(chǎng)，該接口將流體傳熱與湍流k-ε接口相結(jié)合，自動(dòng)添加了非等溫流動(dòng)多物理場(chǎng)耦合，可用于模擬傳熱與流體流動(dòng)之間的耦合。在湍流k-ε物理場(chǎng)中，物理模型選擇了傳統(tǒng)雷諾平均模型（RANS）的k-ε二方程模型，流體為不可壓縮流動(dòng)，壁處理為壁函數(shù)。

1.5.2 定義邊界條件

邊界條件設(shè)置為速度入口和壓力出口邊界條件，開啟重力條件，重力加速度方向?yàn)?z軸方向。殼程飽和蒸汽入口速度為0.4 m/s，入口溫度為155℃，管程冷卻水入口速度為0.01 m/s，入口溫度為50℃。在流體傳熱物理場(chǎng)中，環(huán)境溫度為20℃，環(huán)境絕對(duì)壓力為1個(gè)標(biāo)準(zhǔn)大氣壓，定義換熱器所有壁為薄層。殼程飽和蒸汽相變溫度為151.85℃，冷凝潛熱2108.4 kJ/kg。湍流條件采用湍流長度和湍流強(qiáng)度，湍流長度尺度取0.07r，其中r為入口通道半徑，湍流強(qiáng)度I通過經(jīng)驗(yàn)公式（7）求得。除了對(duì)稱面，其余外部邊界都是熱絕緣壁面。在進(jìn)行數(shù)值模擬時(shí)，因?yàn)槌隹诔３?huì)發(fā)生回流，因此使用壓力出口邊界來代替速度出口邊界條件以獲得更好的收斂效果，并設(shè)定出口為法向流并抑制回流。壁面的邊界選擇無滑移的剛性壁面，除了對(duì)稱面，其余所有外部邊界都是熱絕緣壁面。

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

2.1 速度流線圖

由圖3（a）可見由于入口處截面積較小，流體流入殼程時(shí)，由于流動(dòng)截面積增大，流速快速降低。殼程內(nèi)流體由于換熱器殼壁和折流板的作用，整體呈“Z”字形向出口流動(dòng)，在折流板的缺口下方，形成了明顯的旋渦，而且流速較低，在折流板根部靠近殼壁處，則較少有流線經(jīng)過，說明在該處存在流動(dòng)死區(qū)，區(qū)域內(nèi)流體經(jīng)過量較少，傳熱影響明顯。圖3（b）顯示由于換熱管的阻擋作用，流體沿?fù)Q熱管間隙向殼程內(nèi)部流動(dòng)，當(dāng)遇到殼壁時(shí)，流體沿殼程斜下方流動(dòng)。

圖3 殼程流體的速度流線圖

2.2 溫度分布圖

圖4（a）顯示在換熱管管壁處出現(xiàn)了淺色區(qū)域，主要集中在第一折流板根部的流動(dòng)死區(qū)和第二折流板缺口下方的旋渦區(qū)，這是因?yàn)檫@兩個(gè)區(qū)域內(nèi)，新進(jìn)的流體流量少，流速低，區(qū)域內(nèi)流體與換熱管內(nèi)冷流體換熱時(shí)間相對(duì)較長，新進(jìn)流體補(bǔ)充不充分，導(dǎo)致該區(qū)域內(nèi)殼程流體溫度比附近區(qū)域要低，與管程冷流體間的相對(duì)溫差變小，根據(jù)穩(wěn)態(tài)傳熱基本公式可知，該區(qū)域的傳熱量較高溫區(qū)域要低。從圖4（b）中可看出，高溫飽和蒸汽從換熱器入口接管進(jìn)入換熱器后，溫度逐漸降低，總體的變化趨勢(shì)也呈“Z”字形，這說明溫度的分布與速度流線有著直接的關(guān)系，同時(shí)折流板根部流動(dòng)死區(qū)和缺口下發(fā)旋渦區(qū)，由于傳熱效果差，溫度比同高度的其他區(qū)域要低。

圖4 殼程流體的溫度分布圖

2.3 冷凝水體積分布圖

圖5（a）可以看出，殼程進(jìn)口段至上側(cè)的折流板區(qū)域，凝結(jié)并不明顯，這是因?yàn)楣艹塘黧w經(jīng)換熱器下部換熱后，流體溫度較高，與入口蒸汽溫差較小，傳熱系數(shù)較低，且新進(jìn)入殼程的飽和蒸汽流速快，使得冷凝形成的濕蒸汽向換熱器下方移動(dòng)，而在上側(cè)折流板至殼程出口段，蒸汽冷凝快速發(fā)生。在換熱管管壁處，冷凝水的體積分布也比附近區(qū)域要高。

圖5（b）可以看出冷凝水的體積分布也是呈“Z”字形分布，與殼程流體速度流線和溫度分布情況相吻合。由于不斷通入殼程入口的高溫蒸汽持續(xù)給前方流入換熱器內(nèi)的氣體加熱，進(jìn)口處凝結(jié)量較少，沿速度流線方向上，冷凝水體積占比從進(jìn)口到出口逐漸升高，在出口處及下管板附近，該區(qū)域體積分?jǐn)?shù)接近1，說明飽和蒸汽基本完全冷凝。在折流板背部的流動(dòng)死區(qū)及折流板缺口背部的旋渦區(qū)，冷凝水的體積占比比同高度的其他區(qū)域要高，這是因?yàn)樵搮^(qū)域新進(jìn)的高溫蒸汽流量少，蒸汽溫度較低，與管程冷流體換熱后，凝結(jié)量較大，又因?yàn)樵搮^(qū)域流體流動(dòng)緩慢，凝結(jié)形成的濕蒸汽不能快速流動(dòng)離開該區(qū)域，故區(qū)域內(nèi)冷凝水的體積分布大。

圖5 殼程冷凝水的體積分布圖

3 結(jié)論

以常用的弓形折流板管殼式換熱器為研究對(duì)象，通過對(duì)換熱器殼程三維流場(chǎng)的流動(dòng)與傳熱進(jìn)行數(shù)值模擬，得到以下結(jié)論：

（1）高溫飽和水蒸氣由入口進(jìn)入換熱器，由于換熱管的阻擋，流體沿?fù)Q熱管間隙向換熱器殼壁流動(dòng)，同時(shí)在折流板作用下，呈“Z”字形向換熱器出口流動(dòng)。換熱器殼程內(nèi)部的溫度和冷凝水體積分布也隨流體流動(dòng)情況而變化，整體也呈“Z”字形，與流體速度流向正相關(guān)。

（2）在折流板缺口背部區(qū)域出現(xiàn)旋渦區(qū)，該處流體流速較低，同時(shí)折流板根部出現(xiàn)流動(dòng)死區(qū)，這兩個(gè)區(qū)域溫度比同高度的其他區(qū)域溫度低，因此冷凝水體積分布較大，說明低流速區(qū)域冷凝現(xiàn)象較明顯，不利于傳熱進(jìn)行。

研究結(jié)果有助于對(duì)相變換熱情況下的管殼式換熱器殼程流體流動(dòng)增加認(rèn)識(shí)，并為結(jié)構(gòu)設(shè)計(jì)和優(yōu)化提供參考。