喬宏斌，王 順，胡申華

(南京工程學院 電力仿真與控制工程中心,南京 211167)

水-酒精混合蒸氣在板式換熱器中凝結(jié)時的壓降特性

喬宏斌，王 順，胡申華

(南京工程學院 電力仿真與控制工程中心,南京 211167)

通過設(shè)計搭建板式換熱器凝結(jié)換熱實驗臺，研究了不同酒精質(zhì)量分數(shù)(0%、1%、10%和50%)、蒸氣壓力(70 kPa、80 kPa和90 kPa)和蒸氣流速(15 m/s和25 m/s)下，水-酒精混合蒸氣在板式換熱器中凝結(jié)時的熱側(cè)壓降隨出口蒸氣干度的變化特性.在實驗數(shù)據(jù)的基礎(chǔ)上，采用Lockhart-Martinelli模型計算壓降，并對其進行修正.結(jié)果表明：壓降隨著酒精質(zhì)量分數(shù)的增大而升高，與純水蒸氣相比，在同樣工況下，50%酒精質(zhì)量分數(shù)的混合蒸氣壓降增大近1倍；在相同的酒精質(zhì)量分數(shù)、蒸氣流速和出口蒸氣干度下，蒸氣壓力對壓降影響明顯，蒸氣壓力越大壓降越?。粚嶒炛蹬c擬合值的對比誤差在±20%.

混合蒸氣；板式換熱器；凝結(jié)；壓降；傳熱

符號說明：

A——板片換熱面積，m2

X——Lockhart-Martinelli模型參數(shù)

AC——PHE單通道流通截面積，m2

Δp——總壓降，Pa

cp——比定壓熱容，J/(kg·K)

Δx——干度變化

de——PHE流通等效直徑，m

Δh——蒸氣凝結(jié)潛熱，kJ/kg

f——摩擦因數(shù)

δ——板片厚度、凝結(jié)液膜的厚度，m

L——PHE板片的寬度，m

ρ——密度，kg/m3

p——壓力，kPa

φ——兩相摩擦分液相表觀系數(shù)

q——熱流密度，W/m2

qm—質(zhì)量流量，kg/s

n——流道數(shù)目

T——溫度，K

U——流速，m/s

w——混合蒸氣酒精質(zhì)量分數(shù)，%

x——PHE出口蒸氣干度

xm——PHE內(nèi)部平均蒸氣干度

下標

CW——冷卻水

L——液相

G—— 氣相

in——進口

m——平均值

mix——混合項

out——出口

板式換熱器(Plate Heat Exchanger,PHE)是由一系列具有一定波紋形狀的金屬片疊裝而成的新型高效換熱器.國內(nèi)外學者對板式換熱器的研究涉及內(nèi)容很多，主要研究熱點離不開如何提高其換熱效果和降低壓降.低濃度水-酒精混合蒸氣在物體表面凝結(jié)時，由于Marangoni效應(yīng)，其傳熱系數(shù)相較于純水蒸氣能提高2～8倍[1-2].研究表明這種低濃度提高傳熱系數(shù)的現(xiàn)象在板式換熱器中也存在，因此研究水-酒精混合蒸氣的壓降特性顯得非常重要.在壓降方面，很多學者對氣-液冷凝換熱進行了定性定量地分析.

通過以上文獻可以看出，氣液兩相流在板式換熱器中壓降的研究以制冷空調(diào)行業(yè)中制冷劑的蒸發(fā)和冷凝研究為主.筆者通過設(shè)計和搭建能在正負壓下運行的具有高氣密性的PHE實驗系統(tǒng)，選取能產(chǎn)生Marangoni珠狀凝結(jié)換熱的水-酒精混合物作為混合蒸氣，在不同蒸氣壓力、流速和酒精質(zhì)量分數(shù)下進行蒸氣部分凝結(jié)換熱壓降實驗，分析其變化規(guī)律.并利用實驗結(jié)果建立壓降理論模型，從而指導板式換熱器的設(shè)計，促進中低溫能源的利用.

1 實驗系統(tǒng)

凝結(jié)換熱實驗系統(tǒng)如圖1所示.整個實驗系統(tǒng)從功能上可以分為5部分：混合工質(zhì)循環(huán)系統(tǒng)、冷卻水循環(huán)系統(tǒng)、輔助凝氣器冷卻水循環(huán)系統(tǒng)、抽真空-冷卻系統(tǒng)和數(shù)據(jù)采集系統(tǒng).混合蒸氣從蒸氣發(fā)生器進入板式換熱器中進行凝結(jié)，未完全凝結(jié)的蒸氣進入輔助凝氣器繼續(xù)冷卻.同時，輔助凝氣器也作為整個蒸氣循環(huán)系統(tǒng)的調(diào)壓設(shè)備.空氣作為不凝結(jié)氣體對蒸氣凝結(jié)換熱的影響明顯，為避免在實驗過程中不凝結(jié)氣體進入混合工質(zhì)循環(huán)系統(tǒng)并保證混合蒸氣質(zhì)量分數(shù)不變，實驗臺設(shè)計了抽真空部分，以盡量減少不凝結(jié)氣體對換熱的影響，實現(xiàn)了對系統(tǒng)氣密性的更高要求.

1-安全閥；2-壓力變送器；3-渦街流量計；4-蒸氣側(cè)進熱電偶；5-水側(cè)出熱電偶；6-蒸氣側(cè)差壓變送器；7-蒸氣側(cè)出熱電偶；8-水側(cè)進熱電偶；9-板式換熱器；10-渦輪流量計；11、14-調(diào)節(jié)閥；12-管殼式換熱器；13-抽真空閥；15-數(shù)據(jù)采集模塊；16-計算機；17-冷卻水箱；18-循環(huán)泵；19-管道增壓泵；20-電加熱絲；21-蒸氣發(fā)生器.

圖1 Marangoni凝結(jié)換熱實驗系統(tǒng)圖

Fig.1 Schematic view of the Marangoni condensation test rig

實驗開始時，先運行輔助凝汽器冷卻水循環(huán)系統(tǒng)和抽真空-冷卻系統(tǒng)，抽至系統(tǒng)的最低壓力后通電加熱，當系統(tǒng)壓力逐漸上升至接近工況壓力時，根據(jù)流速需求，調(diào)節(jié)蒸氣發(fā)生器加熱絲的加熱功率，并開啟冷卻水循環(huán)系統(tǒng)，調(diào)節(jié)管殼式換熱器的冷卻水流量，將未凝結(jié)的蒸氣完全凝結(jié)，并維持整個系統(tǒng)的壓力.壓力穩(wěn)定后，開始采集實驗數(shù)據(jù).當板式換熱器的冷卻水溫升高到混合蒸氣幾乎不在板式換熱器中凝結(jié)時，實驗結(jié)束.

實驗過程中，通過流量計測出蒸氣的質(zhì)量流量qm,G.由于混合蒸氣循環(huán)基本上為一個封閉系統(tǒng)，通過下式可計算得到單通道蒸氣的實際流速：

(1)

水-酒精混合蒸氣質(zhì)量分數(shù)通過使用廣泛的UNIFAC基團貢獻法來計算確定[18].

通過測量冷卻水進、出口溫度和流量，可以計算出冷卻水的熱流密度：

qCW=qm,CWcp(TCW,out-TCW,in)/A

(2)

混合蒸氣在板式換熱器中釋放熱量后，其干度變化為

(3)

進入板式換熱器的混合蒸氣干度可視為1，在板式換熱器中干度變化為Δx，取出口蒸氣干度作為依據(jù)，研究壓降隨出口蒸氣干度變化的規(guī)律.

文中參照文獻選取混合蒸氣中的酒精質(zhì)量分數(shù)分別為0%、1%、10%和50%.不同質(zhì)量分數(shù)和壓力下的蒸氣飽和溫度如表1所示.實驗壓力條件為70 kPa、80 kPa和90 kPa，流速取15 m/s和25 m/s.

表1 不同質(zhì)量分數(shù)和壓力下對應(yīng)的飽和溫度

Tab.1 Saturated temperature at different vapor pressures and ethanol mass fraction

酒精質(zhì)量分數(shù)/%壓力/kPa飽和溫度/℃09096.88093.67090.119096.68093.47089.9109095.48092.27088.6509089.08086.07082.6

在換熱過程中，冷卻水和蒸氣為逆流換熱方式.其中板片尺寸為212 mm×465 mm，接口進、出口直徑均為50 mm，其他參數(shù)見表2.不確定度分析結(jié)果見表3.

表2 板式換熱器幾何參數(shù)

表3 不確定度分析

在PHE的蒸氣凝結(jié)過程中，壓降是其重要的特性參數(shù).蒸氣凝結(jié)總壓降由實驗臺的進、出口差壓變送器獲得.蒸氣凝結(jié)總壓降由兩相摩擦壓降Δpf、重力壓降Δpg、加速壓降Δpa和附加壓降Δpadd構(gòu)成，其表達式如下：

Δp=Δpf+Δpg+Δpa+Δpadd

(4)

Wang等[3]研究發(fā)現(xiàn)重力壓降、加速壓降和附加壓降分別占總壓降的0.2%、4%和2.1%，則在板式換熱器中，兩相摩擦壓降是影響蒸氣凝結(jié)總壓降的主要部分.因此主要對兩相摩擦壓降進行分析研究，并在計算結(jié)果上進行修正，以對蒸氣側(cè)總壓降進行預(yù)測.

2 實驗數(shù)據(jù)分析

2.1 酒精質(zhì)量分數(shù)對熱側(cè)壓降的影響

圖2(a)、圖2(b)分別為70 kPa和90 kPa壓力下板式換熱器蒸氣側(cè)壓降隨出口蒸氣干度的變化曲線.由圖2可知，在不同酒精質(zhì)量分數(shù)下壓降變化規(guī)律相同，都隨著出口蒸氣干度的增大而升高.在不同壓力下，酒精質(zhì)量分數(shù)對板式換熱器熱側(cè)壓降影響較大，隨著酒精質(zhì)量分數(shù)的不斷增大，板式換熱器熱側(cè)壓降也不斷升高.在70 kPa時，酒精質(zhì)量分數(shù)為0%、1%和10%的壓降曲線差別較明顯，可以比較清楚地看出其變化趨勢.當壓力達到90 kPa時，各酒精質(zhì)量分數(shù)間的壓降變化十分接近.

(a)

(b)

從圖2可以發(fā)現(xiàn)，相同出口蒸氣干度下，酒精質(zhì)量分數(shù)越大壓降越高.高酒精質(zhì)量分數(shù)時壓降升高1倍左右.隨著酒精質(zhì)量分數(shù)的增大，混合蒸氣的凝結(jié)潛熱隨之減小，在同樣的冷卻強度下，混合蒸氣干度變化更加劇烈，氣液界面也更加紊亂，造成了更大的流動阻力.從微觀上說，酒精質(zhì)量分數(shù)的增大意味著混合蒸氣中酒精分子的數(shù)目增多，酒精分子的體積是水分子的3.7倍，遠大于水分子，導致酒精分子與水分子以及水分子之間的碰撞幾率增大，流體各層之間的切向力會增大，從而增大了流動阻力.

2.2 蒸氣壓力對熱側(cè)壓降的影響

圖3給出了不同蒸氣壓力下板式換熱器熱側(cè)壓降的變化.從圖3可以看出，不同蒸氣壓力下，蒸氣側(cè)壓降都隨出口蒸氣干度的增大而升高.出口蒸氣干度越大，說明混合蒸氣中氣相所占比例越大，混合蒸氣的體積流量要比低出口蒸氣干度時大，則氣液兩相摩擦壓降升高，壓損隨之升高.同樣的實驗條件下，蒸氣壓力越高，壓降越低.

(a)

(b)

2.3 混合蒸氣流速對熱側(cè)壓降的影響

圖4(a)和圖4(b)為壓力70 kPa，酒精質(zhì)量分數(shù)分別為1%和50%時，不同混合蒸氣流速下板式換熱器熱側(cè)壓降隨出口蒸氣干度的變化曲線.從圖4可以看出，在相同酒精質(zhì)量分數(shù)下，混合蒸氣流速對壓降影響較為明顯，25 m/s時的壓降明顯高于15 m/s時.這是由于混合蒸氣流速越大，摩擦阻力越大，從而導致壓降升高.同時，由于板式換熱器是由3塊板片組成的單流程單通道，流通截面積較小，這將進一步使壓降升高.

(a)

(b)

3 實驗數(shù)據(jù)與預(yù)測關(guān)聯(lián)式的對比分析

由于阻力過大將引起板式換熱器的壓降過高，這將影響板式換熱器的運行情況.Lockhart-Martinelli模型已成功應(yīng)用于預(yù)測板式換熱器的凝結(jié)壓降.利用該模型預(yù)測實驗兩相壓降，需先分別計算液相和氣相的壓降.蒸氣側(cè)進、出口液相和氣相的質(zhì)量流量可以通過傳熱量計算得到，因此單相的摩擦因數(shù)可由下式計算得到：

(5)

則液相和氣相的壓降分別由式(6)和式(7)計算得出：

(6)

(7)

Lockhart-Martinelli模型參數(shù)X和兩相摩擦分液相表觀系數(shù)的計算式如下:

(8)

(9)

式(9)中的常數(shù)C可由表4確定.

表4 Lockhart-Martinelli模型常數(shù)C確定值

Tab.4Cvalue determination for the Lockhart-Martinelli model

液體-氣體流動狀態(tài)C湍流-湍流20層流-湍流12湍流-層流10層流-層流5

考慮到板式換熱器氣液兩相流動的復(fù)雜性，將式(9)中的C取為16.則兩相摩擦壓降為

Δpf=φ2ΔpL

(10)

引入修正系數(shù)K對式(10)進行修正，即實驗中的整體總壓降為

Δp=KΔpf

(11)

對于式(11)的待求修正系數(shù)，通過Matlab編寫程序，進行實驗數(shù)據(jù)線性回歸擬合，計算獲得修正系數(shù)K=1.316.其中上述各物性參數(shù)按國際標準水和水蒸氣熱力性質(zhì)IAPWS-IF97公式編寫的Matlab子程序調(diào)用計算獲得.

圖5給出了壓降實驗值和擬合值的對比曲線.從圖5可以看出，統(tǒng)計結(jié)果誤差范圍處于±20%，擬合效果較好.

圖5 實驗值與擬合值的誤差

4 結(jié) 論

(1) 隨著出口蒸氣干度的增大，換熱量減小，混合蒸氣凝結(jié)量隨之減少，壓降升高.

(2) 在不同壓力下，壓降隨著酒精質(zhì)量分數(shù)的增大而升高.

(3) 在相同酒精質(zhì)量分數(shù)、蒸氣流速和出口蒸氣干度的情況下，高壓下的壓降比低壓下的壓降高.與純水相比，在低酒精質(zhì)量分數(shù)(1%)時，壓降幾乎相同；在中酒精質(zhì)量分數(shù)(10%)和高酒精質(zhì)量分數(shù)(50%)時，壓降明顯升高.

(4) 隨著混合蒸氣流速的增大，摩擦阻力也隨之增大,從而板式換熱器的熱側(cè)壓降升高.

(5) 凝結(jié)壓降采用應(yīng)用廣泛的Lockhart-Martinelli模型進行修正，結(jié)果表明實驗值與擬合值的誤差范圍處于±20%.

[1] UTAKA Y,WANG Shixue.Characteristic curves and the promotion effect of ethanol addition on steam condensation heat transfer[J].International Journal of Heat and Mass Transfer,2004,47(21):4507-4516.

[2] 胡申華,嚴俊杰,王進仕.非等溫表面Marangoni凝結(jié)特性的研究[J].動力工程學報,2012,32(1):27-30.

HU Shenhua,YAN Junjie,WANG Jinshi.Study of Marangoni condensation characteristics on non-isothermal surfaces[J].Journal of Chinese Society of Power Engineering,2012,32(1):27-30.

[3] WANG L K,SUNDEN B,YANG Q S.Pressure drop analysis of steam condensation in a plate heat exchanger[J].Heat Transfer Engineering,1999,20(1):71-77.

[4] 倪曉華,夏清,蕭淵.板式換熱器的換熱與壓降計算[J].流體機械,2002,30(3):22-25,32.

NI Xiaohua,XIA Qing,XIAO Yuan.Calculations of heat transfer and pressure drop in plate heat exchangers[J].Fluid Machinery,2002,30(3):22-25,32.

[5] YAN Y Y,LIO H C,LIN T F.Condensation heat transfer and pressure drop of refrigerant R-134a in a plate heat exchanger[J].International Journal of Heat and Mass Transfer,1999,42(6):993-1006.

[6] KUO W S,LIE Y M,HSIEH Y Y,etal.Condensation heat transfer and pressure drop of refrigerant R-410A flow in a vertical plate heat exchanger[J].International Journal of Heat and Mass Transfer,2005,48(25/26):5205-5220.

[8] LONGO G A,GASPARELLA A.Heat transfer and pressure drop during HFC refrigerant vaporisation inside a brazed plate heat exchanger[J].International Journal of Heat and Mass Transfer,2007,50(25/26):5194-5203.

[9] LONGO G A.R410A condensation inside a commercial brazed plate heat exchanger[J].Experimental Thermal and Fluid Science,2009,33(2):284-291.

[10] LONGO G A.Heat transfer and pressure drop during hydrocarbon refrigerant condensation inside a brazed plate heat exchanger[J].International Journal of Refrigeration,2010,33(5):944-953.

[11] HAO Peng,XIANG Ling.Neural networks analysis of thermal characteristics on plate-fin heat exchangers with limited experimental data[J].Applied Thermal Engineering,2009,29(11/12):2251-2256.

[12] NILPUENG K,WONGWISESS.Two-phase gas-liquid flow characteristics inside a plate heat exchanger[J].Experimental Thermal and Fluid Science,2010,34(8):1217-1229.

[15] LONGO G A,ZILIO C.Condensation of the low GWP refrigerant HFC1234yf inside a brazed plate heat exchanger[J].International Journal of Refrigeration,2013,36(2):612-621.

[16] LONGO G A,ZILIO C,RIGHETTI G,etal.Condensation of the low GWP refrigerant HFO1234ze(E) inside a brazed plate heat exchanger[J].International Journal of Refrigeration,2014,38(1):250-259.

[17] LEE E,KANG H,KIM Y.Flow boiling heat transfer and pressure drop of water in a plate heat exchanger with corrugated channels at low mass flux conditions[J].International Journal of Heat and Mass Transfer,2014,77(4):37-45.

[18] WANG Jinshi,YAN Junjie,HU Shenhua,etal.Marangoni condensation heat transfer of water-ethanol mixtures on a vertical surface with temperature gradients[J].International Journal of Heat and Mass Transfer,2009,52(9/10):2324-2334.

Pressure Drop in Plate Heat Exchanger During Condensation of Ethanol-Water Mixed Vapor

QIAOHongbin,WANGShun,HUShenhua

(Electrical Power Simulation and Control Engineering Center,Nanjing Institute of Technology,Nanjing 211167,China)

Experimental tests were carried out in a plate heat exchanger to study the effects of outlet vapor quality on the changes of hot side pressure drop during condensation of ethanol water mixed vapor at different ethanol mass fractions (0%,1%,10% and 50%),different vapor pressures (70 kPa,80 kPa and 90 kPa),and different vapor flow rates (15 m/s and 25 m/s).Based on the experimental data,the Lockhart-Martinelli model was used to calculate and correct the pressure drop.Results show that the pressure drop increases with the rise of ethanol mass fraction,which would be doubled for the mixed vapor with ethanol mass fraction of 50% as compared to pure water vapor,under same experimental conditions.Under conditions with the same ethanol concentration,vapor flow rate and outlet vapor quality,the vapor pressure has obvious effects on the pressure drop; the higher the vapor pressure is,the smaller the pressure drop will be.The error between experimental data and fitted results is within 20%.

mixed vapor; plate heat exchanger; condensation; pressure drop; heat transfer

2015-12-22

2016-02-16

國家自然科學基金資助項目(51166013)；南京工程學院校級科研基金資助項目(YKJ201533)

喬宏斌(1969-)，男，江蘇寶應(yīng)人，講師，碩士，主要從事火電廠的仿真開發(fā)及凝結(jié)換熱方面的研究.電話(Tel.)：17705188894；E-mail：123697618@qq.com.

1674-7607(2017)01-0033-06

TK121

A 學科分類號：470.10