孫啟鵬

（洛陽(yáng)理工學(xué)院，洛陽(yáng) 471005）

0 引言

帶折流板的管殼式換熱器廣泛應(yīng)用在大型動(dòng)力電站、核電站、石油化工、制冷空調(diào)工程等工業(yè)部門(mén)中，而其中約一半以上涉及殼側(cè)氣液的兩相流動(dòng)和傳熱。自本世紀(jì)三十年代起，對(duì)裝有折流板的管殼式換熱器殼側(cè)單相壓降的研究一直在進(jìn)行中。在實(shí)際應(yīng)用的換熱器中，由于加工﹑裝配和管束布置的限制，不可避免地在管束與殼體內(nèi)壁之間形成一定的間隙。在這些間隙中流動(dòng)的流體對(duì)橫掠管束的主流產(chǎn)生分流作用，使殼側(cè)的總驅(qū)動(dòng)壓降減小，也使殼側(cè)實(shí)際參與換熱的流體量減少，削弱了實(shí)際換熱器的換熱能力。為了分析和研究實(shí)際換熱器殼側(cè)的流動(dòng)特性，Tinker提出了殼側(cè)流動(dòng)模型[1]，將殼側(cè)流體分為錯(cuò)流、旁路流及泄漏流等幾個(gè)流路；Palen和Taborek對(duì)Tinker模型進(jìn)行了改進(jìn)，形成目前廣為接受的殼側(cè)流動(dòng)模型。

現(xiàn)有的研究主要集中在各個(gè)流路的單相流動(dòng)特性上，并取得了一定的成果。但是工業(yè)中50%以上的管殼式換熱器涉及殼側(cè)兩相流動(dòng)和換熱，對(duì)于殼側(cè)兩相流動(dòng)特性的研究，目前則主要局限于無(wú)旁路和泄漏的理想換熱器，對(duì)于各個(gè)流路的兩相流動(dòng)特性及其對(duì)實(shí)際換熱器殼側(cè)的流型﹑壓降﹑含氣率的研究，迄今未見(jiàn)有公開(kāi)的文獻(xiàn)報(bào)道。

1 試驗(yàn)裝置與測(cè)試系統(tǒng)

本文在常溫常壓下，通過(guò)對(duì)純主流路和主流、旁流共存流路，應(yīng)用空氣—水兩相混合物研究了旁路流對(duì)TEMA—E型換熱器的壓降及流量分配特性的影響，并提出了考慮旁路的錯(cuò)流區(qū)兩相壓降通用計(jì)算關(guān)聯(lián)式。殼側(cè)流動(dòng)模型及各流路的定義見(jiàn)文獻(xiàn)[2,3]：

試驗(yàn)測(cè)試系統(tǒng)如圖1所示。試件采用有機(jī)玻璃制成。換熱器水平放置，折流板垂直切割，無(wú)泄漏流A和E。試件的幾何尺寸詳見(jiàn)表1。

表1 換熱器結(jié)構(gòu)參數(shù)

試驗(yàn)中先從近似于橫掠管束的B流路開(kāi)始進(jìn)行流動(dòng)特性試驗(yàn)，之后加入旁路流C，研究旁路對(duì)殼側(cè)壓降的影響，以及氣液流量在各流路的分配情況。

2 試驗(yàn)結(jié)果與分析

B流路為換熱器殼側(cè)橫掠管束的主流路，目前對(duì)這一流路兩相流型和壓降的研究，大多在理想方腔殼體內(nèi)進(jìn)行，且有意將錯(cuò)流區(qū)區(qū)段加長(zhǎng)，以便于觀察純橫掠管束的兩相流流型。本文在實(shí)際換熱器模型中增加了旁路擋板，以研究?jī)上嗔黧w在實(shí)際換熱器中流過(guò)純B流路的流動(dòng)特性。

兩相流動(dòng)總壓降Δp可以表示為：

其中，ΔpF為摩擦阻力壓降，Δpg為重力壓降，Δpa為加速壓降。

本文研究的管殼式換熱器為水平放置，折流板垂直切割，不存在重位壓降。對(duì)殼側(cè)錯(cuò)流區(qū)局部含氣率的測(cè)量結(jié)果表明，錯(cuò)流區(qū)內(nèi)兩相分布是不均勻的。因此，錯(cuò)流區(qū)和窗口區(qū)兩相壓降中應(yīng)包含由于含氣率變化引起的加速壓降，但是要準(zhǔn)確地得出測(cè)壓平面處含氣率的值是十分困難的。一般情況下，加速壓降的值相對(duì)于摩擦壓降的值而言是比較小的，而且，在殼側(cè)兩相流壓降的分析中，只要不發(fā)生相變，大多不考慮加速壓降[5]。對(duì)于橫掠管束的兩相摩擦壓降的計(jì)算式，本文采用Ishihara[6,7]推薦的公式，得到結(jié)果為：

圖1 試驗(yàn)測(cè)試系統(tǒng)

其中F2表示由于旁路C的存在，使錯(cuò)流區(qū)壓降減小的無(wú)量綱參數(shù)。當(dāng)ReLB＞2000時(shí)，i＝L；當(dāng)ReLB≤2000時(shí)，i＝G。于是，錯(cuò)流區(qū)兩相壓降表示為：

其中，ΔpiB表示全部氣體或液體單獨(dú)流過(guò)B流路時(shí)的壓降。

根據(jù)方程（4），對(duì)本試驗(yàn)可以將流過(guò)B＋C流路的兩相流體視為全部通過(guò)B流路，由式（2、3）求出 ，同時(shí)可求出ΔpiB，而ΔpTPBC為試驗(yàn)實(shí)測(cè)值，于是可以求出F2的大小。根據(jù)本試驗(yàn)的結(jié)果，在ReLB＞2000或ReLB≤2000條件下，F(xiàn)2的值為0.27～0.5的范圍內(nèi)，顯然，旁路的存在擴(kuò)大了錯(cuò)流區(qū)的流通截面積，同時(shí)，旁路的流動(dòng)阻力較小，從而使錯(cuò)流區(qū)整體流動(dòng)阻力減小，錯(cuò)流區(qū)兩相壓降下降，壓降的平均降幅即為F2的平均值，在40%左右，其大小取決于旁路與B流路流通截面積的相對(duì)大小。

在實(shí)際換熱器中，旁路一般總是存在的。當(dāng)旁路的寬度不同時(shí)，即使管排方式不變，仍然會(huì)使錯(cuò)流區(qū)的兩相流動(dòng)阻力特性和B、C流路的氣液分布發(fā)生變化，從而使有旁路時(shí)的錯(cuò)流區(qū)兩相摩擦壓降計(jì)算式發(fā)生變化。

對(duì)此，本文提出這樣一種思路：對(duì)于純橫掠管束的兩相流動(dòng)（B流路），在經(jīng)過(guò)大量的理想管束試驗(yàn)后，有可能得到比較通用的兩相壓降關(guān)聯(lián)式。旁路C只是繞管束的流動(dòng)，其流動(dòng)阻力特性與管排方式基本無(wú)關(guān)，而與旁路流通截面積直接相關(guān)，其截面積的相對(duì)大小決定著旁路的流動(dòng)Re數(shù)的大小，也決定著其對(duì)主流路的分流量。因此，在錯(cuò)流區(qū)有旁路時(shí)，描述這一流動(dòng)的兩相壓降計(jì)算式應(yīng)當(dāng)以B流路的壓降計(jì)算公式為基礎(chǔ)，結(jié)合C流路的結(jié)構(gòu)特點(diǎn)，對(duì)B流路的兩相壓降公式進(jìn)行合理的修正，得到一種通用的兩相壓降計(jì)算關(guān)聯(lián)式。因此，引入無(wú)量綱修正系數(shù)bcs：

由本試驗(yàn)結(jié)果統(tǒng)計(jì)得：在ReLB＞2000時(shí)，ncs的平均值為2.20，而當(dāng)ReLB≤2000時(shí)，ncs的平均值為2.0。

這樣，錯(cuò)流區(qū)的通用兩相壓降計(jì)算式為：

從式中可以看到，當(dāng)無(wú)旁路時(shí)，即SC＝0，則bcs＝1，F(xiàn)2＝1，上述公式完全恢復(fù)為純B流路的計(jì)算式，說(shuō)明無(wú)量綱修正系數(shù)bcs的物理意義是明確的。

圖2給出了應(yīng)用式（8）（9）得到的錯(cuò)流區(qū)壓降值與實(shí)測(cè)壓降值（B＋C流路）的比較，可以看到，90%的預(yù)測(cè)值與實(shí)測(cè)值的偏差在±25%以內(nèi)，說(shuō)明上述處理方法是合理的。

為計(jì)算錯(cuò)流區(qū)內(nèi)氣、液流量在B、C流路中的分布，建立了將殼側(cè)錯(cuò)流區(qū)的兩相流動(dòng)簡(jiǎn)化為分相流動(dòng)的模[2]，SLC＋SGC＝SC，SC為旁路總流通截面積，旁路中氣液之間的剪切力通過(guò)可變的SLC、SGC來(lái)體現(xiàn)。于是對(duì)殼側(cè)中部SB可直接應(yīng)用純B流路的兩相公式，并可簡(jiǎn)寫(xiě)為：

對(duì)旁路的純液部分SLC和純氣部分SGC可分別建立如下方程：其中fLC、fGC由單相旁路阻力系數(shù)公式求出。

根據(jù)質(zhì)量守恒有：

圖2 錯(cuò)流區(qū)通用壓降公式的預(yù)測(cè)值與B＋C試驗(yàn)值的比較

由于B、C流路為同方向并聯(lián)流動(dòng)，可假設(shè)ΔpTPB＝ΔpLC＝ΔpGC＝ΔpTPBC（錯(cuò)流區(qū)壓降）。于是，對(duì)式（10）～（14），已知總來(lái)流流量及已測(cè)出的錯(cuò)流區(qū)壓降，待求量為B、C流路的氣液流量及SLC，方程封閉，并可通過(guò)非線性方程組的迭代求解方法（即Newton法）解出。

由上述模型得出的計(jì)算結(jié)果表明：相對(duì)于單相流動(dòng)而言，兩相流動(dòng)由于流型的變化（氣液質(zhì)量流速的變化），使主流路B與旁路C各自氣液流量占總流量的比例范圍波動(dòng)較大。在大液量、小氣量（如泡狀流和間歇狀流）下，可以看到液相的波動(dòng)范圍較小，氣相的波動(dòng)范圍較大；而在大氣量、小液量下（如分層流和環(huán)狀流），此時(shí)液相的波動(dòng)范圍要大于氣相的波動(dòng)范圍。這說(shuō)明，殼側(cè)兩相流在B、C流路中的氣液分配比例與兩相流流型和氣液質(zhì)量流速有關(guān)，兩相流體在B、C流路中的分布也是非均勻的。

3 結(jié)論

1）應(yīng)用Ishihara模型的兩相壓降公式，得到了B流路的兩相摩擦壓降表達(dá)式。

2）本文采用無(wú)量綱結(jié)構(gòu)修正因子，以B流路的兩相壓降公式為基礎(chǔ)建立了錯(cuò)流區(qū)通用兩相壓降計(jì)算式，在一定程度上解決了錯(cuò)流區(qū)兩相壓降關(guān)聯(lián)式的任意性。

3）建立了錯(cuò)流區(qū)分相流動(dòng)模型。在已知兩相流量和錯(cuò)流區(qū)及窗口區(qū)壓降的條件下，實(shí)現(xiàn)了B、C流路氣液流量的分離，得到B、C流路氣液流量各自占總流量的份額。計(jì)算結(jié)果表明：B、C流路中，氣液流量的分布是不均勻的。主流路B和旁路C中的氣液各自占相應(yīng)總流量的比例在不同的流型下明顯不同，且比例值的波動(dòng)范圍較大。

[1] Tinker T. Shell Side Heat Transfer Characteristic of Segmentally Baffled Shell and Tube Heat Exchangers [J].ASME, 1949, 21(4): 47-130.

[2] 呂彥力, 孫啟鵬, 陶文銓, 徐斌. 旁路流對(duì)換熱器殼側(cè)氣液兩相流動(dòng)特性的影響[J]. 河南科技大學(xué)學(xué)報(bào), 2006,27(4): 22-25.

[3] 徐斌,王啟杰,馬志豪.換熱器殼側(cè)泄漏流對(duì)氣液兩相分布的影響[J]. 河南科技大學(xué)學(xué)報(bào), 2003, 24(2):48-51.

[4] Grant, I. D. R., and Chishlom, D., Two-Phase Flow on the Shell-Side of a Segmental Baffled Shell-and-Tube Heat Exchanger, [J]. Heat Transfer, 1979, Vol. 101, 38-42.

[5] Williams, C. L. Peterson, A. C. Two-Phase Flow Patterns with high-pressure Water in a Heated Four-Rod Bundle, [J].Nuclear Science and Engineering, 1978, 68, 155-169.

[6] Ishihara, K., Palen, J. W., and Taborek, J., Critical Review of Correlation for Predicting Two-Phase Flow Pressure Drop Across Tube Banks, [J].Heat Transfer Engineering,1980, 1, (3): 23-33.

[7] Briggs, A., and Sabaratnam, S., Condensation of Refrigerant on a Single Tube and a Tube Bank, [J]. Energy Res., 2003, 27: 301-314.