豎直翅片間液橋體積計算模型

莊大偉，楊藝菲，胡海濤，丁國良

?

豎直翅片間液橋體積計算模型

莊大偉，楊藝菲，胡海濤，丁國良

(上海交通大學制冷與低溫工程研究所，上海200240)

管翅式換熱器在析濕工況下會在翅片間形成大量液橋，液橋體積會影響換熱器中冷凝水的排出。提出一種將液面彎曲線沿液橋三相接觸線積分的方法計算豎直翅片間液橋體積，其中液橋三相接觸線采用橢圓方程描述，液面彎曲線模型基于Young-Laplace方程的計算結(jié)果擬合得到。通過液橋三相接觸線與液面彎曲線觀測實驗對提出的液橋體積計算模型進行了驗證，結(jié)果表明，95%的計算結(jié)果與實驗的誤差范圍在±15%內(nèi)，平均誤差為7.12%。

液橋體積；豎直翅片；模型；數(shù)值分析；實驗驗證

引 言

翅片管式換熱器廣泛應用于空調(diào)系統(tǒng)[1-3]，而且經(jīng)常運行在濕工況下[4-6]。當翅片溫度低于濕空氣的露點溫度時，濕空氣中的水蒸氣會冷凝析出，并在翅片表面形成小液滴[7-8]，小液滴會隨冷凝過程的不斷進行進一步生長并互相靠攏，從而在翅片間形成液橋。豎直翅片間的液橋會阻塞空氣流道，惡化換熱以及增大壓降[9-11]。因此，研究豎直翅片間的液橋很重要。

液橋的研究包括3個方面：接觸線、接觸角和體積。現(xiàn)有研究表明豎直翅片間液橋的接觸線為橢圓形[11-12]；接觸角可以通過與方位角有關(guān)的三次多項式表達[12]。液橋的體積決定液橋受到的重力，目前已有將液橋假設(shè)為雙曲線體的體積模型[13]。但觀察發(fā)現(xiàn)豎直翅片間液橋不能簡單地采用雙曲線體計算液橋的體積，因此需要建立一個能夠計算豎直翅片間液橋體積的模型，用于評價液橋的排水特性。

液橋體積模型可以通過對液固接觸面和氣液接觸面的積分求得，其中液固接觸面是由接觸線包圍而成，氣液接觸面是由彎曲液面沿接觸線旋轉(zhuǎn)而成，所以建立液橋體積模型的關(guān)鍵是彎曲液面曲線和接觸線方程。彎曲液面和橢圓接觸面會受到表面張力和重力的影響，導致液橋接觸角在不同方位下的不同，因此在計算液橋體積時需要考慮重力和表面張力對液橋形狀的影響。目前液橋接觸線方程已有研究成果[12-16]，包括水平液橋的圓形接觸線[13-14]和豎直的橢圓形接觸線[15-16]。對于本工作研究的豎直液橋，接觸線為橢圓形，橢圓方程與液橋尺寸、翅片間距和表面特性有關(guān)。目前接觸角的研究是針對水平板間液橋的[17-19]，當液橋形成在豎直平板間時，在重力的作用下液橋?qū)l(fā)生形變，其接觸角會隨方位角變化而變化，因此不能用水平液橋接觸角模型計算豎直液橋的接觸角。

本工作考慮到在重力和表面張力作用下液橋的接觸角隨方位角的變化，建立了能夠計算豎直翅片間液橋體積的模型。

1 數(shù)學模型

1.1 建模對象

豎直翅片間形成的液橋如圖1所示。液橋由空氣-水交界面和水-翅片交界面兩部分包裹而成。通過已知的翅片表面潤濕性、翅片間距和液橋尺寸可以得到空氣-水彎曲線和水-翅片接觸線，然后分別得到空氣-水交界面和水-翅片交界面，最后通過積分的方法求得液橋的體積。

1.2 水-翅片交界面的接觸線方程

水-翅片交界面的形狀如圖2所示。其形狀可以表示為橢圓方程

式中，表示接觸線橢圓的半長軸，可以反映液橋的尺寸；表示接觸線橢圓的長短比[12]，可以通過式(2)求得

式中，表示Bond數(shù)，表示重力加速度，表示水的密度，表示翅片間距，表示水的表面張力系數(shù)。

1.3 空氣-水交界面的彎曲線方程

豎直平板間的彎曲線輪廓如圖3所示。彎曲線方程可以由Young-Laplace基本方程[式(4)]推導得到。方程左項是考慮重力作用下彎曲線液面的內(nèi)外壓差，可由式(5)求解得到；方程右項可以由曲率計算公式[式(6)]求解得到[20]。將式(5)和式(6)代入式(4)可以得到液橋彎曲線的曲率公式[式(7)]。

式中，Δ是交界面壓差，1和2是主曲率半徑，Δ是水和濕空氣密度差，是重力加速度，是彎曲線上任意一點的縱坐標，是表面張力系數(shù)。

對式(7)左右兩側(cè)同時進行積分，如式(8)所示，即可以得到彎曲線的斜率公式[式(9)]，式(9)中的系數(shù)可以通過式(10)計算。

式(9)是彎曲線的斜率方程，但該方程并沒有封閉解析解。因此，為了更方便地求解，本工作引入四次方程[式(11)]擬合彎曲線，并通過式(12)和式(13)所示的邊界條件確定四次方程的待定系數(shù)。

式中，1、2和3是待定系數(shù)，0表示彎曲線最低點0的縱坐標。

因此，最終的彎曲線形狀的四次方程可以表示如下

圖4是解析方程與四次擬合方程關(guān)于彎曲線在軸右半部分的對比，由圖可知四次曲線與解析曲線吻合較好，因此可以用四次曲線表示液橋的彎曲線的形狀，從而為進一步積分求解液橋體積提供方便。

1.4 液橋體積模型

式(14)中的參數(shù)和是與方位角有關(guān)的變量，具體關(guān)系如下所示[12]

式中，()表示隨方位角變化的橢圓接觸線弦長；0表示方位角為0時的橢圓弦長，即橢圓長半軸；()表示隨方位角變化的液橋接觸角。

液橋體積可以通過積分求得。對于的積分從–/2到/2，對于的積分從0到2π，如圖5所示。液橋的體積可以用積分的形式表達為

2 模型驗證

2.1 實驗目的

通過可視化實驗觀測液橋輪廓，并對比實驗中液橋體積與模型計算得到的液橋體積結(jié)果，進而達到驗證模型的目的。

2.2 實驗原理與裝置

本實驗的工作原理是在豎直平板間滴注液體形成液橋，通過CCD相機拍攝液橋輪廓，并采用圖像處理技術(shù)進行分析。

實驗需要測量的量包括平板的前進接觸角A、液橋長度和平板間距。通過這3個實驗量作為輸入量代入計算模型中計算液橋體積，并對比液橋體積的模型計算結(jié)果與實驗滴注的液橋體積值，從而得到模型預測值與實驗值的誤差。

液橋體積的實驗驗證裝置包括液體注射系統(tǒng)、翅片樣件系統(tǒng)、拍攝系統(tǒng)和旋轉(zhuǎn)系統(tǒng)，如圖6所示。

2.3 實驗步驟

本實驗包括兩個步驟。

（1）在豎直翅片間注射蒸餾水，形成液橋。

豎直翅片由兩個平板組成，通過調(diào)節(jié)墊片厚度可以調(diào)節(jié)翅片間距。翅片樣件通過螺絲與接觸角測量儀測試平臺進行固定。通過注射器向翅片間注射蒸餾水，從而形成液橋。

（2）拍攝液橋輪廓，并獲得液橋體積數(shù)據(jù)。

通過旋轉(zhuǎn)系統(tǒng)可以使拍攝角度達到0°～90°。CCD照相機可以在0°時拍攝液橋的正視圖，在90°時拍攝液橋的側(cè)視圖，從而獲得液橋的基本參數(shù)。通過接觸角測量儀自帶的液體注射系統(tǒng)可以控制滴入液橋的體積，從而與模型預測值進行對比。

2.4 驗證結(jié)果與討論

圖7顯示了模型預測與實驗數(shù)據(jù)關(guān)于液橋體積的誤差，其中符號×表示無重力計算模型與實驗的誤差，符號○表示有重力計算模型與實驗的誤差。無重力計算模型可以在±40%的誤差范圍內(nèi)描述100%的點，平均誤差為20.15%；有重力計算模型可以在±15%的誤差范圍內(nèi)描述95%的點，平均誤差為7.12%。

對于考慮了重力作用的計算模型，預測結(jié)果與實驗結(jié)果的誤差是由于對固體表面的前進接觸角和后退接觸角描述的不準確引起的。在實際的實驗樣件中，固體平板是表面特性不均勻的，在平板上的各點其前進接觸角和后退接觸角不一致，通過實驗對平板上、下、左、右、中5個不同位置點進行測量，取平均值作為實驗的測量值，見表1。在模型預測中，固體表面特性是均勻的，各點的前進接觸角和后退接觸角一致，直接采用實驗得到的平均值作為輸入?yún)?shù)。因此，模型預測的結(jié)果與實驗結(jié)果存在誤差。

表1 平板不同位置點的前進接觸角和后退接觸角

豎直翅片間液橋的最小接觸角會在前進接觸角和后退接觸角之間變化。在工程計算中，為了簡化液橋體積的計算過程，最小接觸角通常被認為近似等于后退接觸角。這種近似很可能會對液橋體積的計算造成較大的誤差，因此需要了解液橋最小接觸角的取值對液橋體積計算的影響。

圖8反映了最小接觸角取值對液橋體積計算結(jié)果的影響。由圖可知，液橋的最小接觸角會對液橋體積的計算結(jié)果造成很大的影響。采用min=R的近似計算液橋體積會造成最大17.3%的計算誤差。因此，在計算液橋體積時，不能用后退接觸角直接代替液橋的最小接觸角進行計算，而應該采用本文使用的式(16)進行求解。

3 結(jié) 論

（1）豎直翅片間液橋體積可以通過對彎曲線沿著接觸線進行積分求得。

（2）水-翅片交界面的接觸線可以用與長短比有關(guān)的橢圓方程表示；空氣-水交界面的彎曲線可以由四次方程對彎曲線方程進行描述。

（3）液橋形狀的模型預測結(jié)果和可視化圖像吻合較好；液橋體積的模型計算結(jié)果可以在±15%的誤差范圍內(nèi)描述95%的點，平均誤差為7.12%。

（4）液橋體積的計算不能將后退接觸角近似等于最小接觸角，這種近似會造成最大17.3%的計算誤差。

符 號 說 明

Bo——Bond數(shù) c——待定系數(shù) d——平板間距，mm g——重力加速度，9.81 m·s-2 k——待定系數(shù) L——長軸長度，mm p——壓力，Pa R——主曲率半徑 V——體積，μl β——伸長比 γ——方位角，(°) θ——接觸角，(°) ρ——密度，kg·m-3 σ——表面張力，N·m-1 下角標 A——前進 R——后退

References

[1] MA X K, DING G L, ZHANG Y M,Airside heat transfer and friction characteristics for enhanced fin-and-tube heat exchanger with hydrophilic coating under wet conditions[J]. International Journal of Refrigeration, 2007, 30(7): 1153-1167.

[2] WANG C C, LIN Y T, LEE C J. An airside correlation for plain fin-and-tube heat exchangers in wet conditions[J]. International Journal of Heat and Mass Transfer, 2000, 43(10): 1869-1872.

[3] 馬有福, 袁益超, 陳昱, 等. 翅片螺距對鋸齒螺旋翅片換熱管特性的影響[J]. 化工學報, 2011, 62(9): 2484-2489. MA Y F, YUAN Y C, CHEN Y,Effects of fin pitch on heat transfer and flow resistance of serrated spiral-finned-tube banks[J]. CIESC Journal, 2011, 62(9): 2484-2489.

[4] WANG Q, HAN X H, SOMMERS A,A review on application of carbonaceous materials and carbon matrix composites for heat exchangers and heat sinks[J]. International Journal of Refrigeration, 2012, 35(1): 7-26.

[5] 閔敬春, 彭曉峰, 王曉東. 豎壁上液滴的脫落直徑[J]. 應用基礎(chǔ)與工程科學學報, 2002, 10(1): 57-62. MIN J C, PENG X F, WANG X D. Departure diameter of a drop on a vertical plate[J]. Journal of Basic Science and Engineering, 2002, 10(1): 57-62.

[6] 楊藝菲, 莊大偉, 胡海濤, 等. 濕工況下平翅片平面凝水形成及運動過程的數(shù)值模擬與實驗驗證[J]. 化工學報, 2014, 65(S2): 140-147. YANG Y F, ZHUANG D W, HU H T,Numerical simulation and experimental validation of water condensing and moving on plain-fin surface under dehumidifying conditions[J]. CIESC Journal, 2014, 65(S2): 140-147.

[7] LIU Z L, HUANG L Y, GOU Y J,Experimental investigations of frost release by hydrophilic surfaces[J]. Frontiers of Energy and Power Engineering in China, 2010, 4(4): 475-487.

[8] ZHUANG D W, HU H T, DING G L,. Numerical model for liquid droplet motion on vertical plain-fin surface[J]. HVAC and R Research, 2014, 20(3): 332-343.

[9] KORTE C, JACOBI A M. Condensate retention effects on the performance of plain-fin-and-tube heat exchangers: retention data and modeling[J]. Journal of Heat Transfer, 2001, 123(5): 926-936.

[10] WANG C C, LEE W S, SHEU W J,Investigation of the airside performance of the slit fin-and-tube heat exchangers[J]. International Journal of Refrigeration, 1999, 22(8): 595-603.

[11] QI Z G. Water retention and drainage on air side of heat exchangers—a review[J]. Renewable and Sustainable Energy Reviews, 2013, 28:1-10.

[12] YANG Y F, ZHUANG D W, HU H T,Observation and description of the shape of water bridge retaining between vertical plane surfaces[J]. International Journal of Refrigeration, 2016, 64:20-31.

[13] HARRIS, J W, STOCKER H. Handbook of Mathematics and Computational Science[M]. New York: Springer-Verlag, 1998.

[14] CHEN Y C, ZHAO Y Z, GAO H L,Liquid bridge force between two unequal-sized spheres or a sphere and a plane[J]. Particuology, 2011, 9(4): 374-380.

[15] AKBARI A, HILL R J, VAN DE VEN T G M. Liquid-bridge breakup in contact-drop dispensing: liquid-bridge stability with a free contact line[J]. Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 2015, 92(2): 022404.

[16] MESEGUER J, SLOBOZHANIN L A, PERALES J M. A review on the stability of liquid bridges[J]. Advances in Space Research, 1995, 16(7): 5-14.

[17] VOGEL T I. Stability of a liquid drop trapped between two parallel planesⅡ[J]. SIAM Journal on Applied Mathematics, 1989, 49: 1009-1028.

[18] PETKOV P V, RADOEV B P. Statics and dynamics of capillary bridges[J]. Colloid Surface A, 2014, 460: 18-27.

[19] CHEN H, TANG T, AMIRFAZLI A. Modeling liquid bridge between surfaces with contact angle hysteresis[J]. Langmuir, 2013, 29(10): 3310-3319.

[20] BATEMAN H. Partial Differential Equations of Mathematical Physics[M]. UK: Cambridge University Press, 1959: 170.

Model for calculating water bridge volume retained between vertical fins

ZHUANG Dawei, YANG Yifei, HU Haitao, DING Guoliang

(Institute of Refrigeration & Cryogenics Engineering, Shanghai Jiao Tong University, Shanghai 200240, China)

When a fin-and-tube heat exchanger worked as dehumidifier, lots of water bridges would form between vertical fins and the water volume of these bridges could affect discharge of water condensate. A model for the calculation of water bridge volume was established by integration of meniscus curve along triple contact line of the water bridge. An ellipse equation was used to simulate the triple contact line while a model of meniscus curve at the water-air interface was obtained from calculation results of the Young-Laplace equation. The model for calculating water bridge volume was validated by experimental study on triple contact lines and meniscus curves. The model predicted water volumes of these bridges 95% times aligned with the experimental data within the deviation limit of ±15%, which had a mean deviation of 7.12%.

water bridge volume; vertical fin; model; numerical analysis; experimental validation

2016-04-01.

DING Guoliang, glding@sjtu.edu.cn

10.11949/j.issn.0438-1157.20160411

TK 124

A

0438—1157（2016）10—4080—06

國家自然科學基金項目（51576122）；國家自然科學基金創(chuàng)新研究群體科學基金項目（51521004）；上海市優(yōu)秀學術(shù)帶頭人計劃項目（16XD1401500）。

2016-04-01收到初稿，2016-06-03收到修改稿。

聯(lián)系人：丁國良。第一作者：莊大偉（1987—），男，博士后。

supported by the National Natural Science Foundation of China (51576122), the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (51521004) and the Program of Shanghai Academic Research Leader (16XD1401500).