立式三分螺旋折流板冷凝器的數值計算

林 麗， 陳亞平， 吳嘉峰， 齊雋楠

（1.東南大學能源與環境學院，南京 210096；2.煙臺市建筑設計研究股份有限公司，煙臺 264600）

立式三分螺旋折流板冷凝器的數值計算

林 麗1，2， 陳亞平1， 吳嘉峰1， 齊雋楠1

（1.東南大學能源與環境學院，南京 210096；2.煙臺市建筑設計研究股份有限公司，煙臺 264600）

利用立式三分螺旋折流板冷凝器之傾斜折流板的疏液強化凝結傳熱是一種創新的改進.以數值模擬的方式，研究立式三分螺旋折流板冷凝器模型中傾斜折流板的疏液和螺旋通道內的流動汽液分離對凝結換熱強化的影響.模擬結果表明，凝結換熱效果隨折流板傾斜角度減小而強化，且在含不凝性氣體時，傾斜角為15°，25°和35°的3種三分螺旋折流板冷凝器方案和弓形折流板換熱器方案比無折流板光管束Nusselt凝結模型方案的平均凝結換熱系數分別高出89.7%，75.7%，58.2%和33.1%.

立式冷凝器；三分螺旋折流板；數值模擬；凝結換熱；液膜

Key words：vertical condenser；trisection helical baffle heat exchangers；numerical simulation；condensation heat transfer；liquid film

立式冷凝器相比臥式冷凝器，具有占地面積小、抽管束方便和結構簡便等優點［1-2］.然而，傳統弓形折流板結構的立式冷凝器由于其殼側的凝結換熱系數明顯低于臥式冷凝器，限制了立式冷凝器的應用.例如，電廠高、低壓加熱器采用臥式換熱器結構時，不僅本身占地面積較大，還要留有與管束長度相當的抽管束用的空間，且水平方向抽取管束的操作也不方便；另一方面，高大的汽輪機廠房有著現成的起吊設備，如果能夠通過強化立式換熱器的管外凝結換熱，推廣使用立式高、低壓加熱器，不僅能有效減小占地面積，還能使其結構更加簡單，無需設置滾輪導向軌等部件，并減少維修時間，這對于發電行業極有吸引力；同樣，立式螺旋折流板冷凝器對于制冷、化工等生產過程中廣泛應用的冷凝設備的結構改進也有參考或借鑒意義.

強化凝結傳熱，一直是科研的熱點課題.目前，工業上一般僅能實現穩定的膜狀凝結.減薄液膜層的厚度，及時將傳熱表面上產生的凝結液排除是提高豎壁凝結換熱系數的關鍵手段.已有采用縱向金屬掛線［3］、縱槽管、螺紋槽管等強化方案［4］，使凝結液在表面張力的作用下流向溝槽的谷底，減薄波峰和凸出部分的液膜.彭曉峰等［5］曾提出一種分段冷凝、中間排液的強化管內凝結技術思路，利用凝結的初始段液膜較薄、換熱系數較高的原理，構造分段的短管傳熱面.陳亞平等［6-7］根據上述一些研究思路，在三分螺旋折流板換熱器基礎上，又提出了立式三分螺旋折流板冷凝器方案.其傾斜折流板不僅具有固定管束，抑制管束振動的作用，還具有將管束分隔成短管段，形成“排液盤”的功能.每塊扇形折流板下游處折邊形成擋液堰，引導凝結液流向外圍殼體內壁，直接流向底部液囊而使管束表面液膜減薄，從而強化凝結傳熱.

本文通過數值模擬計算的方法研究立式三分螺旋折流板冷凝器的性能，計算管束表面液膜厚度、蒸汽凝結量和凝結換熱系數等的沿程變化以及折流板傾斜角對凝結性能強化的影響，有利于更好地表述采用三分螺旋折流板強化立式冷凝器凝結傳熱的機理，顯示可能達到的強化傳熱效果，指導不同工況下的立式冷凝器設計，是進行性能實測試驗的補充.管束壁面的凝結傳熱模型基于Nusselt［8］理論提出，并綜合文獻［9-12］的研究成果，考慮了氣液交界面剪切力等影響因素.由于不凝性氣體對工業換熱器中的凝結過程具有很大影響，忽略不凝性氣體影響的計算結果與實際情況相差甚遠，因此盡管目前該領域的理論還很不完善，本文仍采用了Martin等［13］推薦的不凝性氣體質量分數減弱因子的方法，考慮了不凝性氣體的存在對凝結傳熱的弱化.將三分螺旋折流板的立式冷凝器與無折流板和弓形折流板冷凝器的傳熱數據相比較，以此模擬和揭示三分螺旋折流板強化凝結傳熱的作用機理.

1 立式三分螺旋折流板冷凝器模型

1.1 計算模型

立式三分螺旋折流板及弓形折流板冷凝器的結構模型如圖1所示，其結構參數如表1所示，螺旋折流板及弓形折流板參數如表2所示，該模型參數與擬進行的性能試驗是一致的.

圖1 立式螺旋折流板與弓形折流板冷凝器結構示意圖Fig.1 Structure diagrams of vertical trisection helical baffle and segmental baffle condensers

表1 立式三分螺旋折流板與弓形折流板冷凝器結構參數Tab.1 Structure parameters of vertical trisection helical baffle condenser and segmental baffle condenser

表2 立式三分螺旋折流板與弓形折流板冷凝器折流板參數Tab.2 Baffles structure parameters of vertical trisection helical baffle condenser and segmental baffle condenser

在立式三分螺旋折流板冷凝器中，蒸汽由三分螺旋折流板導向在管內呈螺旋流動，傾斜掠過冷凝管并在管外凝結，凝結液沿管壁流下，遇折流板的阻擋后沿傾斜折流板流動，折流板上的液膜流到殼體內壁，然后沿殼體內壁流至換熱器底部.而在立式弓形折流板冷凝器中，蒸汽由弓形折流板引導在管內呈“Z”字形流動，蒸汽在管外凝結，凝結液沿管壁流下，在弓形折流板處堆積直至從折流板缺口處溢出.蒸汽在立式冷凝器的管子上汽、液兩相流動和凝結過程的物理模型如圖2所示.參考蒸汽膜狀凝結的Nusselt理論，并考慮蒸汽流動對液膜的剪切作用，假設：a.管子表面覆蓋一層液膜，在有不凝性氣體存在時，液膜表面還覆蓋一層擴散層，液膜和擴散層的厚度相對于換熱管管徑都足夠小，故可不考慮管子的曲率；b.在動量方程中忽略慣性力、壓力梯度和對流項；c.忽略表面張力的效果；d.凝結液膜的流動是穩定的層流，忽略粘性擴散；e.凝結液膜的物理性質恒定；f.不考慮蒸汽螺旋流動離心力的作用，但考慮了蒸汽隨凝結過程進行流速的變化.取沿管長的豎直方向為x方向，與換熱管垂直的水平方向為y方向建立直角坐標系.

實際工作中，錄音錄像可能會因為一些客觀原因而中斷，例如設備故障、損壞，天氣情況惡劣或者電量、存儲空間不足等。原則上，錄音錄像中斷的應當暫停訊問等取證工作，待錄音錄像恢復之后繼續進行。如果出現緊急情情況不能中斷取證工作，而專業錄音錄像又無法使用時，可以允許使用手機等其他設備替代，但必須立即向監察機關領導報告，事后出具書面說明，但是該錄音錄像的證明力可能會受到影響。

對于純蒸汽，忽略蒸汽流動條件下可得Nusselt經典理論解，管壁液膜傳熱可以看成是凝結液的熱傳導，液膜局部換熱系數hNu（x）為

式中，λ為液膜導熱系數；δ（x）為液膜厚度沿管長x分布；r為蒸汽凝結潛熱；ρc為凝結液密度；μc為凝結液動力粘度；ts為蒸汽飽和溫度；tw為冷凝管壁溫.

由于不凝性氣體對凝結傳熱的影響不可忽略，故采用圖2所示的雙膜模型，按照Martin等［13］推薦的方法計算，其總凝結換熱系數hz及液膜和擴散層的換熱系數的計算公式分別為

圖2 立式螺旋折流板冷凝器單根管上的凝結液膜物理模型示意圖Fig.2 Condensation diagram on a tube of vertical helical baffle condenser

式中，hc和hv分別為含不凝性氣體條件下液膜側和擴散層側凝結換熱系數；ti為氣液交界面溫度；tv為蒸汽溫度；Yg為不凝性氣體的質量分數；f（Yg）為不凝性氣體減弱因子.在本文中將不凝性氣體質量分數取為3%.

液膜厚度可按照δ（x）=λ／hc（x）計算.考慮蒸汽流動對凝結液膜的剪切作用時，液膜動量方程及其邊界條件為

式中，u為液膜沿x方向的流動速度；τix為氣液交界面x方向剪切力.利用式（7）及邊界條件（7a）整理得速度u（y）的分布

式中，τix=τisinθ.采用Wallis［14］的計算方法τi=fρv／2.并根據Wallis管內剪切力的摩擦因子計算公式構造管外剪切力的摩擦因子計算式

式中，R為傳熱管外半徑；uv為蒸汽流速wv的x方向分量.單位管長管外一周液膜質量為

單位時間單位面積管外蒸汽微分凝結量dm為

dm，τix和δ（x）均為未知量，τix亦與δ（x）相關，可利用能量平衡計算dm.根據液膜導熱量、擴散層導熱量與蒸汽凝結放出潛熱熱量的能量平衡

由此可求得dm，迭代可得δ（x）.

弓形折流板凝結模型與立式三分螺旋折流板冷凝器有所不同.弓形折流板背流側的低速區為蒸汽流動死區，占總面積20%，同樣采用含不凝性氣體的凝結模型計算；流速較高區域采用上述考慮交界面剪切力的凝結模型算法.由于弓形折流板立式冷凝器中蒸汽流速不均勻的特點，導致不同區域的管外液膜厚度不同，為方便計算弓形折流板立式冷凝器液膜厚度與凝結換熱系數，本文采用不同區域管外液膜厚度與凝結換熱系數的平均值.

1.2 數值模擬

以立式三分螺旋折流板凝結冷凝器為例，利用Matlab軟件對其換熱情況及凝結液膜的分布進行數值模擬，主要步驟如下：

a.先假定交界面溫度ti=（tw+tv）／2，液膜初始厚度δ，蒸汽凝結量mc0=0；

b.計算帶不凝性氣體條件下的hv和hc；

c.由式（13）和（12）更新ti和dm；

d.更新δ（x）；

e.循環迭代直至最近兩次蒸汽凝結量之差小于設定誤差輸出數據；

f.更新速度等計算參數，計算下一管長直至計算至冷凝器長度.

2 數值模擬結果與分析

圖3顯示了蒸汽含3%不凝性氣體，蒸汽初始流速為10 m／s，蒸汽溫度與壁面溫差為15 K條件下，3種不同傾斜角的立式三分螺旋折流板冷凝器δ，dm，mc和hz沿管長分布情況.作為比較，圖中同時給出了弓形折流板和無折流板的光管束方案以及在純蒸汽凝結條件下Nusselt經典解的數據.

圖3（a）中，含不凝性氣體的無折流板方案（Nu -gas）以及在純蒸汽的Nusselt凝結模型（Nu）的液膜厚度都沿管長逐漸增加，而三分螺旋折流板和弓形折流板凝結模型（seg）中液膜厚度均呈周期性分布.螺旋折流板方案每一螺旋周期內液膜沿管長逐漸增厚，遇折流板沿折流板排泄到殼體內壁流至殼體底部，從而減薄了傳熱管表面液膜厚度.弓形折流板液膜厚度變化亦呈現周期分布，但其折流板的疏液能力差，且在缺口處的折流板支撐距離是中間部分的2倍，管外液膜逐漸增厚，使得其折流板處平均液膜較螺旋折流板方案的厚.圖3（a）表明，純蒸汽Nusselt凝結模型液膜最厚，但是由于此時只有一層熱阻，其平均換熱系數是最高的.在含不凝性氣體的條件下，無折流板方案的平均液膜最厚，弓形折流板方案其次，而螺旋折流板方案的平均液膜厚度隨傾斜角減小而遞減.

圖3（b）中純蒸汽Nusselt凝結模型和含不凝性氣體時無折流板方案的微分凝結量曲線均為連續遞減型，但前者的數值比后者大得多，純蒸汽Nusselt凝結模型的平均微分凝結量最大，含不凝性氣體時無折流板方案的平均微分凝結量最小.其余方案在凝結的第一個周期的前段曲線是重合的；經過折流板的疏液作用呈現周期性從大到小變化，即初期因液膜較薄，微分凝結量較大，隨液膜厚度增加微分凝結量逐漸減小.弓形折流板方案的平均微分凝結量低于螺旋折流板方案；螺旋折流板方案的平均微分凝結量隨傾斜角增加而減小.

圖3（c）顯示了各方案的積分凝結量沿管長的分布.可見在3／4管長內純蒸汽Nusselt凝結模型的凝結量是最大的，在后1／4管長段被傾斜角15°螺旋折流板方案超出；傾斜角25°和35°的螺旋折流板方案分列第三和第四；弓形折流板方案和含不凝性氣體時無折流板方案分別為倒數第二和第一.在含不凝結性氣體條件下，采用螺旋折流板方案能夠接近甚至超過純蒸汽凝結的凝液量和傳熱系數，可見螺旋折流板疏液對強化凝結傳熱具有重要的意義.

圖3 不同方案的立式冷凝器內δ，dm，mc和hz分布Fig.3 Distributions ofδ，dm，mcand hzin different vertical condenser schemes

圖3（d）顯示了各方案的局部換熱系數沿管長的分布.與圖3（b）的曲線類似，純蒸汽Nusselt凝結模型和含不凝性氣體時無折流板方案的微分凝結量曲線均為連續遞減型，純蒸汽Nusselt凝結模型的平均換熱系數最大，而含不凝性氣體時無折流板方案的平均換熱系數遠小于其它凝結模型.對于有折流板的方案，凝結換熱系數在每個周期入口處達到最大值，隨管長快速減小.螺旋折流板凝結模型平均局部換熱系數隨傾斜角增加而減小.弓形折流板方案的局部凝結換熱系數低于螺旋折流板方案.傾斜角15°，25°和35°的立式三分螺旋折流板方案和弓形折流板方案凝結換熱系數相對于含不凝性氣體時無折流板方案的比值分別為1.897，1.757，1.582和1.331.

3 結 論

本文通過數值模擬研究立式三分螺旋折流板冷凝模型中傾斜折流板的疏液和流動對凝結換熱的強化作用，建立了含不凝性氣體條件下的不同傾斜角的立式三分螺旋折流板冷凝器方案、弓形折流板方案和無折流板光管束Nusselt凝結模型，并對各方案的凝結液膜厚度、蒸汽凝結量和凝結換熱系數沿管程的分布情況進行了分析比較.其主要結論如下：

a.相同條件下，純蒸汽Nusselt凝結模型的平均液膜最厚，但其只有一層液膜熱阻；含不凝性氣體無折流板光管束Nusselt凝結模型其次；弓形折流板方案第三；不同傾斜角螺旋折流板方案的平均液膜厚度隨傾斜角減小而遞減.有折流板方案的液膜厚度均呈現周期性變化.

b.純蒸汽Nusselt凝結模型和含不凝性氣體無折流板光管束Nusselt凝結模型的微分凝結量和局部換熱系數均呈現先快后慢連續下降的趨勢，但前者的數值比后者大得多；其余有折流板方案的微分凝結量和局部換熱系數均呈現周期性變化.

c.相對于含不凝性氣體時無折流板光管束Nusselt凝結模型方案，傾斜角15°，25°和35°的立式三分螺旋折流板方案和弓形折流板方案的凝結換熱系數的相對比值分別為1.897，1.757，1.582和1.331.

［1］ 胡雷鳴.垂直管外凝結的雙側強化傳熱研究及其在發電廠低壓加熱器上的應用［D］.重慶：重慶大學，2001.

［2］ 梁平，龍新峰，樓波，等.電站立式給水加熱器傳熱強化及無銅化運行研究［J］.上海理工大學學報，2001，23（3）：274-278.

［3］ Thomas D G.Enhancement of film condensation heat transfer rates on vertical tubes by vertical wires［J］. Industrial and Engineering Chemistry Fundamentals，1967，6（1）：97-103.

［4］ 顧維藻，神家銳，馬重芳，等.強化傳熱［M］.北京：科學出版社，1990.

［5］ 彭曉峰，吳迪，張揚.高性能冷凝器技術原理與實踐［J］.化工進展，2007，26（1）：97-104.

［6］ 陳亞平.適合于正三角形排列布管的螺旋折流板換熱器［J］.石油化工設備，2008，37（6）：1-5.

［7］ 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.

［8］ Nusselt W.Dieoberfl?chen kondensation des wasser dampfes［J］.VDI，1916，60：541-546.

［9］ Fujii T，Honda H，Oda K.Condensation of steam on a horizontal tube—the influence of oncoming velocity and thermal condition at the wall［C］∥Proceeding of the 18thNational Heat Transfer Conference，1979.

［10］ Shekriladze I，Gomelauri G.Theoretical study of laminar film condensation of flowing vapor［J］.Heat Mass Transfer，1966，9：581-591.

［11］ Rohsenow W M，Webber JH，Ling A T.Effect of vapor velocity on laminar and turbulent film condensation［J］.Trans ASME，1956，78：1637-1643.

［12］ Sparrow E，Gregg M.A boundary-layer treatment of laminar film condensation［J］.Trans ASME J Heat Transfer，1959，81：13-17.

［13］ Martin-V J M，Jimenez M A，Martin-F F，et al. Comparison of film condensation models in presence of non-condensable gases implemented in a CFD code［J］.Heat Mass Transfer，2005，41（11）：961-976.

［14］ Wallis G B.One-dimensional two-phase flow［M］.New York：McGraw-Hill，1969.

（編輯：董 偉）

Numerical Simulation of Vertical Trisection Helical Baffle Condensers

LINLi1，2， CHENYa-ping1， WUJia-feng1， QIJun-nan1
（1.School of Energy and Environment，Southeast University，Nanjing 210096，China；2.Yantai Architectural Design and Research Co.，Ltd.，Yantai 264600，China）

A novel structure of vertical trisection helical baffle feed water heater can enhance the heat transfer by effectively draining the condensate with inclined baffles and separating water droplets from vapor due to centrifugal force in helix channels.The draining film effect of inclined baffles and the impact of spiral flow on condensation heat transfer in a vertical trisection helical baffle condenser were numerically simulated.These simulation schemes include the Nusselt condensation models under the conditions of pure steam or in presence of non-condensable gas，the vertical trisection helical baffle condenser schemes with baffle inclined angles of 15°，25°and 35° and the segmental baffle scheme in the presence of non-condensable gas.The results demonstrate that the condensation heat transfer increases with the decrease of inclined angle of helical baffles，and in the presence of non-condensable gas，the average condensation heat transfer coefficients of three trisection helical baffle schemes with inclined angles of 15°，25°and 35°and the segmental baffle scheme are 89.7%，75.7%，58.2%and 33.1%higher than that of the bare tube bundle of Nusselt condensation model.

TK 172

A

2013-08-31

國家自然科學基金資助項目（51276035）；江蘇省科技創新與成果轉化專項引導資金項目（BY2011155）第一作者：林 麗（1987-），女，碩士研究生.研究方向：強化傳熱.E-mail：804620928＠qq.com

陳亞平（1956-），男，教授.研究方向：強化傳熱.E-mail：ypgchen＠sina.com

1007-6735（2014）04-0333-05

10.13255／j.cnki.jusst.2014.04.006