郭　亞，林　麗，陳亞平，吳嘉峰

（東南大學(xué)能源與環(huán)境學(xué)院，江蘇南京 210096）

?

立式螺旋折流板冷凝器內(nèi)制冷劑流型和換熱特性的數(shù)值模擬

郭亞*，林麗，陳亞平，吳嘉峰

（東南大學(xué)能源與環(huán)境學(xué)院，江蘇南京 210096）

[摘 要]通過數(shù)值模擬對比研究單、雙頭變角度三分螺旋折流板與變間距弓形折流板立式冷凝器的凝結(jié)換熱性能。建立了設(shè)有擋液堰和疏液缺口的三分螺旋折流板立式冷凝器和弓形折流板立式冷凝器的模型，利用Fluent軟件的混合模型模擬了蒸汽在上述不同立式冷凝器內(nèi)的速度、壓力及凝結(jié)液的分布情況。模擬結(jié)果表明，雙頭變角度三分螺旋折流板立式冷凝器能夠有效提高凝結(jié)換熱系數(shù)，其平均換熱系數(shù)可比弓形折流板方案提高35.7％，且擋液堰與疏液缺口能夠有效排液，防止凝結(jié)液被蒸汽流夾帶，克服了傳統(tǒng)立式冷凝器的水平弓形折流板不能有效疏液的局限。

[關(guān)鍵詞]立式冷凝器；三分螺旋折流板；擋液堰；疏液缺口；冷凝

*郭亞（1989-），男，碩士。研究方向：傳熱傳質(zhì)。聯(lián)系地址：南京市玄武區(qū)四牌樓2號東南大學(xué)南高院103，郵編：210096。聯(lián)系電話：18795856362。E-mail：1017451713@qq.com；陳亞平（1956-），男，教授。聯(lián)系電話：13851729402。E-mail：ypgchen@sina.com。

0 引言

立式冷凝器相比臥式冷凝器，具有占地面積小、結(jié)構(gòu)簡便、抽管束和維修方便等優(yōu)點。然而，傳統(tǒng)弓形折流板結(jié)構(gòu)的立式冷凝器由于管外凝結(jié)液在重力作用下聚集，管束表面的凝結(jié)傳熱熱阻很快沿程增大，導(dǎo)致其殼側(cè)的凝結(jié)換熱系數(shù)明顯低于臥式冷凝器，從而限制了立式冷凝器的應(yīng)用。因此強(qiáng)化立式冷凝器的管外凝結(jié)換熱，對于在發(fā)電行業(yè)推廣立式給水加熱器極具價值。

提高冷凝器傳熱系數(shù)的關(guān)鍵在于減薄凝結(jié)液膜。彭曉峰等[1]曾提出利用凝結(jié)初始段液膜較薄、換熱系數(shù)較高的原理，構(gòu)造分段短管傳熱面加中間排液的思路，并根據(jù)該原理設(shè)計了一種板翅式冷凝器。本課題組借鑒上述思路[2-9]在三分螺旋折流板換熱器基礎(chǔ)上，提出了立式三分螺旋折流板冷凝器方案[10]。該方案主要針對高、低壓給水加熱器等殼側(cè)凝結(jié)的管殼式換熱器，在垂直管束的殼側(cè)設(shè)置周向重疊三分螺旋折流板，將管束分隔成短管段，使傾斜折流板具有“排液盤”的功能，不斷將管束上的冷凝液刮下。蒸汽則將在折流板的作用下沿螺旋通道流動。同時，在折流板下游邊設(shè)置擋液堰和在折流板弧段上設(shè)置疏水缺口，可及時將從管束表面流下的凝結(jié)液沿殼體內(nèi)壁排泄到換熱器底部，有效減薄傳熱管壁外液膜厚度，增加凝結(jié)量，提高換熱系數(shù)。

由于螺旋折流板冷凝器幾何結(jié)構(gòu)發(fā)生變化從而導(dǎo)致殼程流體流動狀態(tài)發(fā)生改變，進(jìn)而造成整體傳熱特性和壓力分布的改變，使得傳統(tǒng)的分析方法不能很好地滿足換熱器的設(shè)計要求。隨著計算數(shù)學(xué)的發(fā)展和眾多成熟商業(yè)軟件的開發(fā)，以及計算機(jī)硬件技術(shù)的不斷進(jìn)步，采用試驗和數(shù)值模擬相結(jié)合的研究方法已逐漸成為許多工程研究領(lǐng)域的主流研究方法。本文采用數(shù)值模擬方法研究立式三分螺旋折流板冷凝器的換熱特性，采用Gambit軟件建立了三分螺旋折流板立式冷凝器結(jié)構(gòu)模型以及帶有擋液堰和疏液缺口的折流板結(jié)構(gòu)模型，并利用Fluent軟件的混合模型模擬蒸汽在上述不同立式冷凝器內(nèi)的速度、壓力及凝結(jié)液的分布情況。本文還將三分螺旋折流板冷凝器與現(xiàn)有的弓形折流板立式冷凝器的換熱特性進(jìn)行了分析比較。

1 數(shù)值計算方法

1.1 幾何模型與網(wǎng)格劃分

三分螺折流板和弓形折流板立式冷凝器結(jié)構(gòu)如圖1所示，冷凝器主要由圓筒形殼體、換熱管束、折流板和管殼程流體進(jìn)、出口接管組成?？紤]到蒸汽的體積流量隨著凝結(jié)過程而逐漸變小，螺旋折流板冷凝器和弓形折流板冷凝器均采用截面積分段逐漸變小的結(jié)構(gòu)。螺旋折流板采用三種傾斜角度，弓形折流板采用三種間距。螺旋折流板采用單頭和雙頭兩種方案，單頭螺旋折流板冷凝器中三段的傾斜角依次為35°、25°和15°；雙頭螺旋折流板換熱器是在單頭螺旋折流板換熱器的基礎(chǔ)上，將螺距較大的35°和25°傾斜角折流板段布置為雙頭螺旋折流板，而螺距較小的15°傾斜角折流板段仍為單頭。折流板結(jié)構(gòu)如圖2所示，在折流板下游邊設(shè)置折彎擋液堰，在折流板的弧段設(shè)置疏液缺口，可有效減薄傳熱管壁外液膜厚度，提高換熱系數(shù)，并及時將從管束表面流下的凝結(jié)液沿殼體內(nèi)壁流到換熱器底部。表1和表2分別為三分螺旋折流板冷凝器和弓形折流板冷凝器參數(shù)表和折流板布置方案表。

圖1　立式冷凝器結(jié)構(gòu)示意圖

圖2　折流板結(jié)構(gòu)圖

表1　三分螺旋折流板冷凝器參數(shù)表

表2　變角度三分螺旋折流板方案冷凝器的折流板布置方案

采用Gambit軟件建立立式三分螺旋折流板與弓形折流板三維模型并劃分網(wǎng)格。圖3為變角度傾斜角三分螺旋折流板方案的網(wǎng)格圖，換熱管處采用加密處理，經(jīng)過網(wǎng)格獨立性檢驗，確定網(wǎng)格數(shù)量約為250萬。

圖3　三分螺旋折流板網(wǎng)格圖

1.2 控制方程

考慮到本模型蒸汽相與液相速度不同且速度分布差異較大的特點，選用Fluent多相流模型中的混合(Mixture)模型進(jìn)行模擬分析?；旌夏Ｐ褪且环N簡化的多相流模型，用于模擬各相有速度不同的多相流，通過假定多相流在短空間尺度上局部平衡的方法，來增強(qiáng)相間的耦合性?；旌夏Ｐ涂梢酝ㄟ^建立兩相和混合相的連續(xù)性、動量和能量方程，相對速度方程，以及第二相(液相)的體積分?jǐn)?shù)方程等關(guān)系式來模擬兩相流動。其連續(xù)性方程、動量方程和能量方程等守恒方程[11-12]如下：

式中：

αk——第k相的體積分?jǐn)?shù)；

n——相數(shù)；

S→Mk、SEm——質(zhì)量源項和能量源項；

Ek——相間傳遞能量；

ρk——第k相密度；

ρm——混合密度；

μm——混合粘性；

λ——有效熱傳導(dǎo)率；

keff——有效→熱傳導(dǎo)率。

ρm、μm、分別由體積分量加權(quán)平均獲得：

不可壓縮相相間傳遞能量Ek= hk，這里hk是顯熱焓值?？蓧嚎s相的相間傳遞能量計算公式為：相對速度（也指滑流速度）被定義為液相c的速度相對于汽相v的速度：

液相的體積分?jǐn)?shù)方程：

式中：

凝結(jié)模型中的特點是汽相質(zhì)量隨凝結(jié)而減小，放出凝結(jié)熱后，產(chǎn)生相應(yīng)質(zhì)量的凝結(jié)液。本文采用質(zhì)量與能量源項來描述傳熱傳質(zhì)過程，液相、汽相的質(zhì)量源項和混合相的能量源項分別為：

其中ΔH為凝結(jié)焓。

1.3 邊界條件和計算設(shè)置

邊界條件的設(shè)置主要有冷熱流體的入口、出口、外殼壁面、折流板以及換熱管壁面。定義冷、熱流體的入口為速度入口條件，給定入口速度以及相應(yīng)的溫度條件和湍流值；定義出口處為自由出流條件；外殼壁面和折流板定義為不可滲透、絕熱邊界條件；換熱管壁面定義為與管內(nèi)冷水換熱的耦合邊界條件；蒸汽定義為主相，水為輔相，并在蒸汽相中導(dǎo)入蒸汽質(zhì)量源項，水相中導(dǎo)入液體質(zhì)量源項，蒸汽混合相中導(dǎo)入能量源項。由于湍流效應(yīng)對流動與傳熱的影響造成計算結(jié)果的收斂不穩(wěn)定。因此，本文模擬采用對旋轉(zhuǎn)流具有針對性的RNG κ-ε湍流模型。本模型采用分離求解器及非穩(wěn)態(tài)隱式格式的計算方法；壓力和速度耦合采用基于交錯網(wǎng)格的SIMPLE算法，動量、能量以及湍流參量的求解采用二階迎風(fēng)格式。定義動量方程迭代收斂的條件為方程中各變量平均殘差絕對值小于10-4；能量方程迭代收斂的條件為方程中各變量平均殘差絕對值小于10-7。

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

圖4顯示了蒸汽在變角度單、雙頭三分螺旋折流板與弓形折流板立式冷凝器內(nèi)由第一段流至第二段的速度矢量疊加圖。由圖4(a)、圖4(b)顯示的主視圖可見在三分螺旋折流板冷凝器內(nèi)蒸汽在螺旋折流板的作用下呈螺旋流動。在折流板背面單頭方案有明顯的低速區(qū)域，而雙頭方案由于折流板布置較密集使周期內(nèi)大部分蒸汽呈現(xiàn)較好的螺旋流動。由圖4(c)顯示的子午面切片圖中可見在弓形折流板冷凝器內(nèi)蒸汽呈近似“Z”字形流動，且在折流板缺口處速度矢量分布較密集，而在折流板背流區(qū)較稀疏。蒸汽主流穿過缺口處后，在慣性作用下以較高速度撞擊下一塊折流板，然后貼著折流板上表面轉(zhuǎn)為橫向沖刷，主流向下穿過下一個缺口，小股流體向上形成回流，在折流板的背面形成流動死區(qū)，如此呈周期性分布。

圖5顯示了蒸汽在變角度單、雙頭三分螺旋折流板與弓形折流板立式冷凝器內(nèi)的壓力分布圖(橫切片壓力分布截取距頂部0.31 m處的數(shù)據(jù))。由圖可見，在三分螺旋折流板方案與弓形折流板方案中蒸汽壓力均沿豎直方向逐漸降低，但弓形折流板方案壓降比三分螺旋折流板方案大得多。雙頭方案由于折流板較多，其壓降比單頭方案略大。在三分螺旋折流板方案中，由于流體在折流板形成的螺旋通道內(nèi)流動時會受到離心力的作用，中心部分的流體在離心力的驅(qū)動下不斷的向外圍流動，因此會造成壓力分布呈現(xiàn)外圍高、中心低。

圖4　蒸汽速度云圖和速度矢量疊加圖

圖5　壓力分布圖（橫切片位置：0.31 m）

圖6顯示了單、雙頭三分螺旋折流板與弓形折流板立式冷凝器內(nèi)凝結(jié)液體積分?jǐn)?shù)分布沿程主視圖和中間段局部子午面切片圖。總體可見，單、雙頭三分螺旋折流板方案下部的凝結(jié)液體積分?jǐn)?shù)比弓形折流板方案的高得多，反映其優(yōu)良的凝結(jié)強(qiáng)化性能。由中間段的局部放大圖可見，三分螺旋折流板方案中管束上被折流板刮下的凝結(jié)液沿折流板傾斜流下，堆積在折流板擋液堰處，沿?fù)跻貉吡髦翚んw內(nèi)壁，再穿過疏液缺口沿殼體內(nèi)壁流下，從而使凝結(jié)液排泄過程連續(xù)快速進(jìn)行。由于雙頭螺旋折流板方案的折流板螺距更短，凝結(jié)液膜較薄，其凝結(jié)效果更好。圖中亦可看出三分螺旋折流板方案兩折流板之間的殼體壁面上有凝結(jié)液區(qū)，可見，設(shè)有疏液缺口與擋液堰的傾斜折流板能夠有效排液，促進(jìn)強(qiáng)化傳熱。相反弓形折流板方案由于水平布置的弓形折流板不能將凝結(jié)液及時排泄，使凝結(jié)液在折流板表面堆積，并被高速流過的主流蒸汽夾帶或沿傳熱管子流下，而沿殼體內(nèi)壁幾乎沒有凝結(jié)液流下。

圖6　凝結(jié)液體積分?jǐn)?shù)分布圖（放大圖位置：0.5 m～0.6 m）

圖7顯示了蒸汽在變角度單、雙頭三分螺旋折流板與弓形折流板立式冷凝器內(nèi)子午面上5根管子表面局部換熱系數(shù)分布圖（螺旋折流板方案的截面位于0.3 m～0.5 m處，弓形折流板方案位于0.3 m～0.5 m處）。由圖可見，雙頭三分螺旋折流板方案換熱系數(shù)的平均值較高，而單頭方案的一個螺旋周期中的上半部分的局部換熱系數(shù)與雙頭方案雖然相近，但其下半部分的數(shù)值則較低。主要由于雙頭折流板使兩折流板之間的換熱管變短，改善了單頭方案周期下部換熱較差的缺陷。在弓形折流板立式冷凝器中也可看出折流板背面由于大部分凝結(jié)液被折流板刮除，使其折流板下部局部換熱系數(shù)較高。但總體上其局部換熱系數(shù)還是比螺旋折流板方案低得多。雙頭和單頭三分螺旋折流板方案的平均換熱系數(shù)可比弓形折流板方案分別提高35.7％和22.4％。

圖7　子午面上五根管子表面的局部換熱系數(shù)分布圖（截面位置：0.3 m～0.5 m）

3 結(jié)論

采用Fluent軟件中Mixture多相模型模擬了三分螺旋折流板立式冷凝器變傾斜角單、雙頭方案的殼側(cè)蒸汽速度、壓力、凝結(jié)水體積分?jǐn)?shù)和局部換熱系數(shù)分布情況，并將其與弓形折流板方案進(jìn)行對比分析，得出以下結(jié)論。

1）三分螺旋折流板立式冷凝器內(nèi)蒸汽沿折流板通道呈螺旋流動，能夠加強(qiáng)流體沖刷管壁表面，減小液膜厚度，傾斜角度分段逐漸減小的折流板能夠有效彌補(bǔ)因凝結(jié)過程蒸汽體積流量的減小而造成的流速下降；將折流板的下游邊折彎成擋液堰和在折流板弧段設(shè)置疏液缺口能夠引導(dǎo)凝結(jié)液流至殼體內(nèi)壁，避免被蒸汽夾帶，有效排液，從而增強(qiáng)殼側(cè)管束上的凝結(jié)換熱。雙頭和單頭三分螺旋折流板立式冷凝器方案的平均換熱系數(shù)可比弓形折流板方案分別提高35.7％和22.4％。

2）三分螺旋折流板凝結(jié)立式冷凝器由于蒸汽流動均勻，流動阻力較小，殼側(cè)壓力降明顯低于弓形折流板立式冷凝器方案。其中雙頭螺旋折流板方案因增加了折流板的摩擦較單頭螺旋折流板方案的流阻略大。

三分螺旋折流板立式冷凝器的綜合換熱性能明顯優(yōu)于傳統(tǒng)的弓形折流板立式冷凝器，為在發(fā)電廠中推廣立式給水加熱器開拓了新的思路。

參考文獻(xiàn)：

[1]彭曉峰,吳迪,張揚(yáng).高性能冷凝器技術(shù)原理與實踐[J].化工進(jìn)展,2007,26(1):97-104.

[2]陳亞平.適合于正三角形排列布管的螺旋折流板換熱器[J].石油化工設(shè)備,2008,37(6):1-5.

[3]陳亞平,王偉晗,李彥晴.三分螺旋折流板換熱器殼側(cè)傳熱特性研究[J].工程熱物理學(xué)報,2010,31(11):1905-1908.

[4]陳亞平,董聰,操瑞兵,等.周向重疊三分螺旋折流板換熱器殼側(cè)性能試驗研究[J].化工學(xué)報,2012,33(7):1194-1196.

[5]CHEN Y P,SHENG Y J,DONG C,et al.Numerical simulation on flow field in circumferential overlap trisection helical baffle heat exchanger[J].Applied Thermal Engineering,2013,50 (1):1035-1043.

[6]黃慶楠,陳杰,楊文剛,等.螺旋折流板管殼式換熱器內(nèi)LNG混合介質(zhì)流動冷凝特性仿真平臺開發(fā)及影響因素分析[J].制冷技術(shù),2015,35(6):18-24,66.

[7]李志文,賴學(xué)江,楊麗明.整體針翅管和螺旋扭曲管在制冷機(jī)換熱器上運用的探討[J].制冷技術(shù),2002,22(2):16-18.

[8]余江海,陸震,范林.蒸發(fā)式冷凝器應(yīng)用現(xiàn)狀及存在問題探討[J].制冷技術(shù),2001,21(2):33-36.

[9]王聰,翟曉強(qiáng),李卉,等.相變蓄冷換熱器性能的實驗研究[J].制冷技術(shù),2014,34(2):11-15.

[10]林麗,陳亞平,吳嘉峰,等.立式三分螺旋折流板冷凝器的數(shù)值計算[C].中國工程熱物理學(xué)會傳熱傳質(zhì)學(xué)學(xué)術(shù)會議論文,重慶,2013.

[11]ALIZADEHDAKHEL A,RAHIMI M.CFD modeling of flow and heat transfer in a thermosyphon[J].International Communications in Heat and Mass Transfer,2010,37(3):312-318.

[12]GHORAI S,NIGAM K D P.CFD modeling of flow profiles and interfacial phenomena in two-phase flow in pipes[J].Chemical Engineering and Processing Journal,2006,45(1):55-65.

Numerical Simulation of Flow Pattern and Heat Transfer Characteristics of Refrigerant in Vertical Condenser with Helical Baffles

GUO Ya*,LIN Li,CHEN Ya-ping,WU Jia-feng
(School of Energy and Environment,Southeast University,Nanjing,Jiangsu 210096)

[Abstract]The condensation heat transfer characteristics of the vertical condensers with variable angled trisection helical baffles with single or double thread and with a variable spanned segmental baffle were numerically studied and compared.The physical models of the trisection helical baffle condensers and the segmental baffle one were developed.The distributions of the velocity and pressure of vapor and the volumes of fraction of condensate in these heat exchangers were simulated with the mixing module of Fluent software.The simulation results show that,the variable angled trisection helical baffle vertical condenser with double-threads could significantly improve the condensation heat transfer coefficient,and the average shell-side heat transfer coefficient of the vertical condenser with dual-thread variable angled trisection helical baffles is 35.7％ higher than that of the variable spanned segmental baffled one.The liquid dam plus drainage gaps on each inclined baffle plate could drain condensate effectively and prevent liquid film entrainment into vapor,thus it overcomes the limitation of the conventional vertical condensers which could not drain condensate effectively with horizontal segmental baffles.

[Keywords]Vertical condensers; Trisection helical baffles; Liquid dam; Drainage gaps; Condensation

基金項目：國家自然科學(xué)基金項目（No.51276035、51206022）；江蘇省科技創(chuàng)新與成果轉(zhuǎn)化專項引導(dǎo)資金項目（No.BY2011155）。

doi：10.3969/j.issn.2095-4468.2016.01.108