水平異型管外加裝排液板降膜流動(dòng)數(shù)值分析

李慧君, 賈怡瓊 , 王 炯

(1.華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003；2.中材節(jié)能股份有限公司，天津 300400)

水平異型管外加裝排液板降膜流動(dòng)數(shù)值分析

李慧君1, 賈怡瓊1, 王 炯2

(1.華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003；2.中材節(jié)能股份有限公司，天津 300400)

液膜熱阻是影響水平管外降膜流動(dòng)傳熱的主要因素之一，采用加速排泄管外液體及減薄液膜厚度的方法可以使其換熱加強(qiáng)。通過(guò)建立底部加裝排液板的異型管物理模型，采用fluent軟件中的VOF模型對(duì)其管外降膜流動(dòng)進(jìn)行了二維數(shù)值模擬，來(lái)研究排液板對(duì)水平異型管外的降膜流動(dòng)與傳熱的影響。并比較了模擬得到結(jié)果與文獻(xiàn)中的實(shí)驗(yàn)數(shù)據(jù)，變化趨勢(shì)吻合較好。結(jié)果表明：加裝排液板的異型管較其光管的平均液膜厚度減薄，管壁的局部Nu數(shù)增大，表明排液板可以加速排液，有助于強(qiáng)化傳熱；蛋形管具有比滴形管更好的排液性能，加裝排液板后排液效果更加顯著。

強(qiáng)化傳熱；數(shù)值模擬；降膜流動(dòng)傳熱；排液板；液膜厚度

0 引言

水平管外降膜蒸發(fā)技術(shù)是一種高效節(jié)能的薄膜蒸發(fā)技術(shù)，其在較低流速及小換熱溫差下具有相對(duì)較高的熱質(zhì)傳遞能力[1]。水平管式降膜蒸發(fā)器廣泛應(yīng)用于工業(yè)生產(chǎn)，其在節(jié)能環(huán)保方面有著十分重要的意義，引起了國(guó)內(nèi)外許多研究者的廣泛關(guān)注。文獻(xiàn)[2]得出了液膜厚度隨管徑、管束布置方式以及管間距等參數(shù)的變化規(guī)律。文獻(xiàn)[3]研究了傳熱情況下水平管外的降膜流動(dòng)，得到了液膜厚度的近似式。文獻(xiàn)[4]利用位移測(cè)微儀對(duì)液膜水平管外周向角范圍為15°～165°的液膜進(jìn)行了實(shí)驗(yàn)研究，確定了液膜厚度在周向角95°～120°的范圍內(nèi)存在最小值。文獻(xiàn)[5]比較了Turbo-CII管和光滑管，得出了Turbo-CII管外周向液膜厚度變化同光管一樣，最小值位于周向角95°～120°的范圍內(nèi)。文獻(xiàn)[6]研究了R-134a在3根垂直排列的水平銅制管外降膜蒸發(fā)的傳熱性能，得出隨著流體溫度的升高，光滑管的換熱系數(shù)增大，在“干斑”出前流量也略有增加。文獻(xiàn)[7]研究了真空狀態(tài)下6排水平圓管束的降膜蒸發(fā)，得出增加管子內(nèi)表面、提高加熱水溫度可以得到更好的換熱性能，并預(yù)測(cè)了傳熱系數(shù)和增強(qiáng)比例的相關(guān)性。管外液膜分布對(duì)于換熱器換熱性能起著至關(guān)重要的影響，液膜過(guò)厚換熱系數(shù)降低，液膜過(guò)薄可能會(huì)導(dǎo)致管壁局部干區(qū)的出現(xiàn)而造成換熱器損害[8]。大部分水平管降膜蒸發(fā)器的換熱管采用圓管，但國(guó)內(nèi)外許多學(xué)者研究表明圓管不是獲得最優(yōu)液膜分布的管型[9～19]。

國(guó)內(nèi)外大多數(shù)文獻(xiàn)都是以水平圓管為研究對(duì)象，少數(shù)涉及到異形管。但針對(duì)異型管外加裝排液板的降膜流動(dòng)研究鮮有報(bào)道。諸多研究者采用了排液板技術(shù)來(lái)強(qiáng)化膜狀冷凝換熱并取得了一定的成果，研究表明加裝排液板起到了加快管外凝結(jié)液的排泄，同時(shí)有效減薄了管外壁面凝結(jié)液膜的厚度，強(qiáng)化了凝結(jié)換熱[20]。因此，可以在水平異型管外采用排液板技術(shù)，使其換熱性能得到提高。為了研究排液板對(duì)降膜傳熱的影響，采用fluent軟件中的VOF模型對(duì)其進(jìn)行模擬，以期為降膜流動(dòng)傳熱提供參考。

1 物理模型

水平蛋形管外帶排液板降膜流動(dòng)的物理模型，如圖1所示。貼體坐標(biāo)系x沿著管壁面切線方向，y則是沿管壁外法線方向，其分速度分別由u、v表示。δ為液膜厚度，t及H分別表示排液板的厚度和高度，θ為圓周角。質(zhì)量流量為2Γ、溫度為T0的液體噴淋到管子頂部，形成液膜并沿管壁面向下流動(dòng)，通過(guò)排液板排離出去。壁面維持恒定溫度Tw。

圖1 水平蛋型管外帶排液板降膜流動(dòng)物理模型

為簡(jiǎn)化計(jì)算，作如下假設(shè)：(1) 工質(zhì)的物性參數(shù)為定常數(shù)；(2) 液膜流動(dòng)為穩(wěn)定的層流；(3)與換熱管管徑相比，液膜厚度是小量；(4) 流體不可壓縮；(5) 液體入口溫度等于飽和蒸汽溫度；(6) 不計(jì)流體相變，傳熱只考慮從管壁到汽液界面；(7) 不考慮排液板的傳熱，排液板起到排液作用。

2 數(shù)值計(jì)算

2.1 異型管及排液板幾何參數(shù)尺寸

采用與φ19的光圓管具有相同換熱面積對(duì)蛋形管、滴型管及不同長(zhǎng)短半軸比的橢圓管進(jìn)行尺寸優(yōu)化設(shè)計(jì)，各異型管及排液板的幾何參數(shù)如圖2所示。

圖2 3種水平異型管及排液板幾何參數(shù)

2.2 無(wú)量綱數(shù)

為了更好地比較不同參數(shù)下管外流動(dòng)及換熱情況，定義如下無(wú)量綱參數(shù)[11]713。

(1)

(2)

(3)

(4)

式中：P為單側(cè)管壁長(zhǎng)度，m；T、T0、Tw分別為液膜內(nèi)流體溫度、入口處流體溫度和壁溫，K；δ為液膜厚度，m；ν為流體運(yùn)動(dòng)粘度，m2/s；λ為熱導(dǎo)率，W/(m·K)；h為局部換熱系數(shù)，W/(m2·K)。

2.3 網(wǎng)格劃分及邊界條件

由于VOF模型對(duì)計(jì)算區(qū)域網(wǎng)格精度要求較高，對(duì)局部網(wǎng)格進(jìn)行加密，圖3為蛋形管網(wǎng)格劃分及邊界條件示意圖。設(shè)定管壁和液相水壁面的接觸角為20°[21]。

圖3 網(wǎng)格模型示意圖

2.4 網(wǎng)格無(wú)關(guān)性驗(yàn)證與模型驗(yàn)證

以蛋形管為例，Γ為0.3 kg/(m·s)時(shí)，采用網(wǎng)格數(shù)量分別為50 214、72 188及166 032。最大相對(duì)誤差小于2.0%，表明計(jì)算結(jié)果與網(wǎng)格無(wú)關(guān)。綜合考慮了計(jì)算精度和計(jì)算時(shí)間，擇優(yōu)選取網(wǎng)格數(shù)為72 188的網(wǎng)格劃分方案，其他異型管采用類似的網(wǎng)格進(jìn)行劃分。

為了驗(yàn)證本模型的正確性，將以E=1.5的橢圓管模擬結(jié)果與文獻(xiàn)[10]77中最相近工況進(jìn)行比較，如表1所示。

表1 本模型及文獻(xiàn)實(shí)驗(yàn)工況

由于本模型中所采用橢圓管幾何參數(shù)小于文獻(xiàn)[10]77，導(dǎo)致在周向角θ<135°時(shí)，模擬結(jié)果普遍高于實(shí)驗(yàn)結(jié)果，但變化趨勢(shì)是一致的，如圖4所示。在周向角θ>135°時(shí)，模擬結(jié)果小于實(shí)驗(yàn)值，主要原因是實(shí)驗(yàn)過(guò)程伴隨有蒸發(fā)相變，蒸發(fā)現(xiàn)象導(dǎo)致液膜波動(dòng)[10]77。液膜厚度的平均誤差小于30%，表明本模型基本正確。

圖4 本模型的液膜厚度與文獻(xiàn)[10]78的比較

2.5 模擬算例結(jié)果分析

取T0=333.16 K的飽和水為流動(dòng)工質(zhì)，壁面溫度為Tw=335.16 K，并設(shè)定外部空間充滿此溫度下的靜止的飽和蒸汽。流體在整個(gè)傳熱過(guò)程中是常物性，其參數(shù)由表2可知。

表2 流體物性參數(shù)

3種截面形狀的異型管的速度分布呈拋物線形，液膜內(nèi)最大速度出現(xiàn)在液膜中間某處，如圖5所示。在管壁處受到壁面的粘附作用，存在流動(dòng)邊界層，流體在貼壁處存在較大的正向速度梯度，迅速增大到主流速度，隨后受到汽液界面靜止的蒸汽剪切作用，速度略有減小。隨著無(wú)量綱長(zhǎng)度X的增大，異型管液膜內(nèi)最大速度逐漸增大，其原因是液膜在重力的作用下，液體沿管壁周向流動(dòng)加快。滴型管及蛋形管加裝排液板后，液膜流動(dòng)速度較其光管大，表明排液板對(duì)液體流動(dòng)起到加速作用。在無(wú)量綱長(zhǎng)度X=5.333附近處，帶排液的橢圓管E=1.5的最大速度較其光管略小，但在管壁周向其他位置處帶排液板的橢圓管沿管壁周向液膜內(nèi)最大速度較其光管的大，表明排液板對(duì)液體流動(dòng)起加速作用同樣適用于橢圓管。隨著X的增大，壁面處速度梯度減小，表明液體沿管壁周向流動(dòng)時(shí)，壁面處液體的流動(dòng)減弱，對(duì)流換熱減小。

圖5 異型管膜內(nèi)局部速度分布

圖6 異型管局部膜內(nèi)無(wú)量綱溫度分布

隨著無(wú)量綱長(zhǎng)度X的增大，各異型管局部膜內(nèi)溫度Θ分布隨之增大，如圖6所示。

圖6表明隨著X的增大，膜內(nèi)的溫度梯度隨之減小，熱邊界層隨之增大，換熱效果也隨之減弱。加裝排液板的異型管膜內(nèi)溫度Θ較其光管要小，表明加裝排液板后，膜內(nèi)的溫度梯度增大，熱邊界層減小，強(qiáng)化了換熱效果。在壁面處流體溫度和管壁面溫度相等，溫度Θ為0，熱量從管壁傳遞給液膜，其溫度逐漸升高，液膜溫度在汽液分界面達(dá)到飽和溫度，其溫度Θ為1。

各異型管的液膜厚度在管頂及管底部的分布較厚，如圖7所示。

圖7 異型管液膜厚度分布

圖7表明，主要是由于管頂受到下落的噴淋液體的沖擊，豎直方向的速度分量幾乎為零，液體沿管兩側(cè)鋪展，液體沿管壁周向下落的速度特別小，造成液體在管頂?shù)亩逊e；在異型管底部受到兩側(cè)流體的匯集，造成液膜急劇增厚。流體沿管壁的流動(dòng)速度相對(duì)于噴淋液體的下落速度小，液膜不能迅速排離管壁面，使液膜積聚。加裝排液板后，各異型管的液膜厚度均較其光管的要薄。對(duì)于橢圓管在X=5.333附近處，帶排液板的橢圓管較其光管略厚。但帶排液板的橢圓管平均液膜厚度較其光管要薄，且在管頂及管底部表現(xiàn)十分明顯。主要原因是加裝排液后，在管頂部受到排液板對(duì)液膜的加速引流作用，液膜厚度跟其光管相比變薄明顯；在管底部，排液板對(duì)匯集的液體起到加速排泄作用，液體能迅速排離換熱管，不至于在管底部匯集過(guò)厚，對(duì)強(qiáng)化換熱有利。

當(dāng)Г=0.45 kg/(m·s)時(shí)，橢圓管E=1.5異型管最厚液膜與最薄液膜的比值為1.52，加裝排液板后，其比值為1.34；并且加裝排液板后，平均液膜厚度減薄了5.5%。表明加裝排液板后，有效地減薄了厚液膜區(qū)，增大了薄液膜區(qū)，使得橢圓管外液膜分布較均勻，薄液膜區(qū)液膜變化趨勢(shì)平緩，可以使“干斑”的可能性減??；同時(shí)，加裝排液板后，管外平均液膜減薄，有利于強(qiáng)化換熱。

由于滴型管與蛋形管尺寸設(shè)計(jì)時(shí)，長(zhǎng)軸的長(zhǎng)度約相同，重力對(duì)液膜的作用可認(rèn)為是一致的。因此可以對(duì)比管型截面形狀對(duì)液膜流動(dòng)的影響。蛋形管在管頂部分液膜厚度稍厚于滴形管，由圖8可知。

圖8 滴型管與蛋形管液膜厚度分布

由于曲率半徑稍大的滴形管液膜在水平方向上鋪展速度較蛋形管的快，使得其厚度較蛋形管的薄。但沿著管壁周向向下流動(dòng)，蛋形管的液膜厚度減薄的速度比滴形管大，主要原因是液體沿蛋形管壁周向向下流動(dòng)時(shí)，液體的流動(dòng)方向與重力方向的夾角小于滴形管，液體受重力加速作用明顯。Γ為0.45 kg/(m·s)時(shí)，蛋形管液膜的平均厚度較滴形管的減小了0.5%，加裝排液板后，蛋形管的平均液膜厚度較滴形管減小4.3%，表明加裝排液板后，液膜明顯變薄，增大了傳熱量。

圖9 異型管管壁局部Nu數(shù)分布

各異型管壁局部Nu數(shù)分布在管頂及管最底部處較大，如圖9所示。主要原因是由于在管頂部受到下落流體的沖擊作用，擾動(dòng)加劇，使得流體與管壁的對(duì)流換熱加強(qiáng)；在管底部由于兩側(cè)流體的匯集及流體脫離換熱管時(shí)的撕裂作用，對(duì)換熱起到一定的強(qiáng)化作用。由圖5、圖6及圖7所示，流體隨著無(wú)量綱長(zhǎng)度X的增大，管壁面處的速度梯度減小，溫度梯度隨之逐漸減小，熱邊界層厚度隨之逐漸增大，使得流體在換熱管壁面處的換熱沿著管壁周向向下逐漸減小。各異型管加裝排液板后，管壁局部Nu數(shù)較其光管要高，主要原因是加裝排液板后，使液膜的流動(dòng)加快，壁面處的速度梯度增大，溫度梯度增大，液膜厚度減薄，增大了傳熱量。

3 結(jié)論

(1)3種截面形狀的異型管加裝排液板后，膜內(nèi)速度提高，膜內(nèi)無(wú)量綱溫度減小，熱邊界層減薄，使管外液體加快排泄，液膜厚度減小，換熱效果增強(qiáng)。

(2)帶排液板的橢圓管液膜厚度分布較其光管要薄，薄液膜區(qū)相對(duì)較大且變化趨勢(shì)平緩，可提高降膜流動(dòng)傳熱系數(shù)；減小“干斑”或者燒毀的可能性，保證了降膜蒸發(fā)設(shè)備的安全運(yùn)行。

(3)蛋形管具有比滴形管更好的排液能力，加裝排液板后，排液效果更加顯著。

[1]何茂剛, 王小飛, 張穎,等.水平管降膜蒸發(fā)器管外液體流動(dòng)研究及膜厚的模擬計(jì)算[J].熱科學(xué)與技術(shù), 2007, 6(4):319-325.

[2]邱慶剛, 陳金波.水平管降膜蒸發(fā)器管外液膜的數(shù)值模擬[J].動(dòng)力工程學(xué)報(bào)，2011, 31(5): 357-361, 374.

[3]RIFERT V G, PODBEREZNY V I, PUTILIN J V, et al.Heat transfer in thin film-type evaporator with profile tubes[J].Desalination, 1989, 74(8):363-372.

[4]HOU H, BI Q C, MA H, et al.Distribution characteristics of falling film thickness around a horizontal tube[J].Desalination, 2012, 285(3): 393-398.

[5]WANG X F, HE M G, FAN H L, et al.Measurement of falling film thickness around a horizontal tube using a laser measurement technique[J].Journal of Physics, 2009, 147(1):175-180.

[6]CHIENAL L H, CHENAL R H.An experimental study of falling film evaporation on horizontal tubes using R-134a[J].Journal of Mechanics, 2012, 28(2): 319-327.

[7]LI W, WU X Y, LUO Z, et al.Falling water film evaporation on newly-designed enhanced tube bundles[J].International Journal of Heat and Mass Transfer, 2011, 54(13): 2990-2997.

[8]WU Y, LI J, WANG N H, et al.Industrial experiments on desulfurization of flue gases by pulsed corona induced plasma chemical process[J].Journal of Electrostatics, 2003, 57(3): 233-241.

[9]齊春華, 徐克, 馮厚軍, 等.橢圓管外液膜流動(dòng)的數(shù)值模擬及傳熱實(shí)驗(yàn)研究 Ⅰ: 數(shù)值模擬研究[J].化學(xué)工業(yè)與工程, 2012, 29(4): 74-79.

[10]齊春華, 馮厚軍, 邢玉雷, 等.橢圓管外液膜流動(dòng)的數(shù)值模擬及傳熱實(shí)驗(yàn)研究 Ⅱ: 實(shí)驗(yàn)研究[J].化學(xué)工業(yè)與工程, 2012, 29(5): 43-47.

[11]羅林聰, 張冠敏, 田茂誠(chéng), 等.異形管降膜傳熱與流動(dòng)數(shù)值模擬研究[J].工程熱物理學(xué)報(bào), 2013, 34(4): 710-714.

[12]李慧君, 王炯.水平滴型管外凝結(jié)液膜分布及換熱特性研究[J].華北電力大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 42(5):94-99.

[13]羅林聰, 潘繼紅, 田茂誠(chéng), 等.管形對(duì)水平管降膜圓周膜厚和 Nusselt 數(shù)的影響[J].化工學(xué)報(bào), 2013, 64(8): 2760-2768.

[14]李慧君, 彭文平, 董楠,等.汽-氣層流凝結(jié)氣液膜內(nèi)物理量分布的研究[J].電力科學(xué)與工程, 2013, 29(8):39-45.

[15]LI Q L, MA L X, Qi S L, et al.Experimental investigation on heat transfer and flow resistance of an 1.6/1.0 elliptic cylinder in cross-flow[J].Journal of Qingdao University of Science and Technology, 2004, 25(1): 54-61.

[16]KIBBOUA R, AZZI A.Laminar film condensation on an elliptical tube embedded in porous media[J].Energy Conversion and Management, 2005, 46(15): 2359-2372.

[17]ESFAHANI J A, MODIRKHAZENI M.Entropy generation of forced convection film condensation on a horizontal elliptical tube[J].Comptes Rendus Mecanique, 2012, 340(7): 543-551.

[18]ZHANG L, YANG S, XU H.Experimental study on condensation heat transfer characteristics of steam on horizontal twisted elliptical tubes[J].Applied Energy, 2012, 97(3): 881-887.

[19]HORVAT A, LESKOVER M, MAVKO B.Comparison of heat transfer conditions in tube bundle cross-flow for different tube shapes[J].International Journal of Heat and Mass Transfer, 2006, 49(5): 1027-1038.

[20]AL-BADRI A R, GEBAUER T, LEIPERTZ A, et al.Element by ele ment prediction model of condensation heat transfer on a horizontal integral finned tube[J].International Journal of Heat and Mass Transfer, 2013, 62(1): 463-472.

[21]LEI X, LI H.Numerical simulation of the behavior of falling films on horizontal plain tubes[J].AIP Conference Proceedings, 2010, 1207(1): 998-1003.

Numerical Analysis of Falling Film Flow Outside a Horizontal Shaped Circular Tube with a Drainage Strip

LI Huijun1, JIA Yiqiong1, WANG Jiong2

(1.School of Energy Power and Mechanical Engineering, North China Electric Power University, Baoding 071003, China;2.Sinoma Energy Conservation Company Limited by Shares,Tianjin 300400,China)

Liquid film resistance is one of the main factors which influence the falling flow heat transfer of horizontal tube.It can be used to enhance the heat transfer by accelerating drainage and reducing the thickness of the liquid film.In order to study the effect of the drainage strip on falling film flow and heat transfer outside a horizontal shaped tube, a physical model of special pipe with a drainage plate at the bottom was established and the VOF model of fluent software was used to two-dimensional numerical simulate it.The numerical simulation results were in good agreement with literature experimental data.The results are as follows.Compared the shaped circular tube with a strip with it without a strip, the average film thickness is thinner and the localNunumber is greater.It is showed that the drainage plate can accelerate drainage and reduce the thickness of liquid film, which helps to enhance the heat transfer.The oval-shaped tube has better drainage performance than the drop-shaped tube, and the effect of drainage with the drainage strip is more significant than before.

heat transfer enhancement；numerical simulation；falling film flow and heat transfer； drainage strip；film thickness

10.3969/j.ISSN.1672-0792.2017.06.008

2017-04-28。

TK124

A

1672-0792(2017)06-0048-07

李慧君(1964-)，男，教授，研究方向?yàn)閺?qiáng)化換熱及數(shù)值計(jì)算，電廠熱力系統(tǒng)的節(jié)能原理及檢測(cè)診斷等。

賈怡瓊(1992-)，女，碩士研究生，主要從事強(qiáng)化換熱及數(shù)值計(jì)算方面的研究工作。

王炯(1988-)，男，碩士研究生，主要從事強(qiáng)化換熱及數(shù)值計(jì)算方面的研究工作。