謝涵廣,胡劍光,王 成,戴干策

(華東理工大學(xué) 化學(xué)工程聯(lián)合國家重點(diǎn)實(shí)驗(yàn)室,上海 200237)

垂直板上開窗區(qū)受限自由膜的形成

謝涵廣,胡劍光,王 成,戴干策*

(華東理工大學(xué) 化學(xué)工程聯(lián)合國家重點(diǎn)實(shí)驗(yàn)室,上海 200237)

液膜是常見的氣液接觸方式之一，依據(jù)單或雙自由面可分為壁面膜和自由膜。不同于傳統(tǒng)液膜研究，板上開窗或開孔后會同時(shí)存在這2種液膜。實(shí)驗(yàn)采用了6種液體(Ka數(shù)52～3000)和10余種不同尺寸的矩形窗(9～1152mm2)，觀察了垂直板上下降液膜通過開窗區(qū)產(chǎn)生的自由膜流動行為。結(jié)果表明，在表面張力、慣性力、重力及黏性力的共同作用下，液體通過開窗區(qū)的方式有繞流、過流、偏流和“背流”4種主要形式。隨著流量的增加，過流又可以有液滴、液柱、液膜以及它們的組合等豐富流型，同時(shí)窗內(nèi)自由膜在表面張力主導(dǎo)下易受擾動而聚并、破碎，流型呈現(xiàn)多樣性和動態(tài)性。當(dāng)流量達(dá)到臨界值時(shí)，窗口內(nèi)會形成完整穩(wěn)定的受限自由液膜，呈現(xiàn)板上壁面膜與窗口內(nèi)受限自由膜交替共存的特殊液膜——孿生液膜。本文關(guān)聯(lián)了臨界Re數(shù)、Ka數(shù)和無量綱幾何參數(shù)Nx，給出了成膜條件的經(jīng)驗(yàn)判別式。孿生液膜具有的特殊波列結(jié)構(gòu)、成膜過程的遲滯現(xiàn)象以及2種液膜和固體壁間相互作用引起的多重干擾等，都體現(xiàn)出與傳統(tǒng)的壁面液膜和自由液膜不同的流動性質(zhì)。實(shí)驗(yàn)結(jié)果將有助于工業(yè)過程中開窗和開孔的幾何優(yōu)化，強(qiáng)化傳遞過程，并進(jìn)一步豐富傳統(tǒng)的液膜研究領(lǐng)域。

孿生液膜；流型轉(zhuǎn)變；成膜；開窗板

0 引 言

液膜流動是一種常見的低速開放式流動，它蘊(yùn)涵著極其豐富的動力學(xué)行為，其中的交界面不穩(wěn)定性、自由表面波、三維時(shí)空圖像等非線性現(xiàn)象一直是流體力學(xué)工作者所研究的熱點(diǎn)[1]。在工業(yè)上也具有廣泛的應(yīng)用和研究，比如涂布流、降膜反應(yīng)器、薄膜傳熱與蒸發(fā)冷凝等等，甚至填料塔、圓盤反應(yīng)器、微流體反應(yīng)器等都可以用液膜流動理論來強(qiáng)化傳質(zhì)反應(yīng)過程。根據(jù)液膜形成特征可以分為壁面膜(一側(cè)有固體壁面支撐，僅有1個(gè)自由面)和自由膜(在流動過程中不接觸任何固體，同時(shí)擁有2個(gè)自由面)[2]。關(guān)于壁面膜，自Kapitza父子[3]開創(chuàng)傾斜板上的薄膜下落實(shí)驗(yàn)后，許多學(xué)者通過不同的測試設(shè)備和數(shù)值計(jì)算方法開展這類研究，發(fā)現(xiàn)了孤立波、三維不穩(wěn)定性、回流等特殊流動現(xiàn)象。Wasden和Dukler[4]通過數(shù)值計(jì)算發(fā)現(xiàn)液膜波動會使垂直于膜厚方向速度分量增加，有效促進(jìn)傳質(zhì)過程。Alekseenko et al.[5]、Chang和Demekhin[6]以及Kalliadasis et al.[7]對液膜流動的實(shí)驗(yàn)方法、穩(wěn)定性分析、數(shù)值計(jì)算等方面進(jìn)行了系統(tǒng)闡述。關(guān)于自由膜，由于2個(gè)自由面的存在，在改變液體流量或物性時(shí)會出現(xiàn)如液滴、液柱、液膜等一系列豐富流型。對它的研究主要集中在液膜波動的穩(wěn)定性與射流破裂，其理論在發(fā)動機(jī)燃燒和噴霧領(lǐng)域等有著重要作用。Squire[8]較早通過理論和實(shí)驗(yàn)分析了靜止氣體中自由液膜的穩(wěn)定性以及射流破裂問題。Brown[9]測定了自由下落液膜的速度分布，發(fā)現(xiàn)對于低黏液體，液膜速度接近于自由落體變化，并初步提出維持液膜穩(wěn)定的韋伯?dāng)?shù)(We)準(zhǔn)則。Lin和Robert[10]通過在自由液膜中心放置一小障礙物，發(fā)現(xiàn)自由膜受擾時(shí)波動會產(chǎn)生2種不同的形式，一類是非對稱的曲波，一類是對稱的膨脹波。隨后更多學(xué)者不斷對相關(guān)理論、模型進(jìn)行改進(jìn)，更加詳細(xì)的著作可以參見Lin[11]，Sirignano與Mehring[12]。此外也有部分學(xué)者研究了下降自由膜的流型隨流量的變化[13]。

對于2類液膜的研究，至今逾半個(gè)世紀(jì)，然由于各自不同的工業(yè)應(yīng)用背景和研究角度，對兩者的研究甚少關(guān)聯(lián)。尤其是自由膜，其流動更加復(fù)雜，研究相對較少，主要集中在穩(wěn)定性分析而缺少傳質(zhì)等方面的研究。而事實(shí)上，在現(xiàn)有化工設(shè)備中，如散堆填料和規(guī)整填料、圓盤反應(yīng)器等，往往在固體壁面上開有大小不一的孔或窗，并非完整的壁面或完全的自由空間。當(dāng)液體完全覆蓋在這樣的開窗板上時(shí)，即會形成板上壁面膜與窗口區(qū)受限自由膜交替共存的現(xiàn)象(可稱之為孿生液膜)，呈現(xiàn)出不同于兩者的新流動特性。更值得一提的是，不少工業(yè)設(shè)備如填料塔，早期的經(jīng)典理論都是在壁面膜基礎(chǔ)上建立并沿用至今，然而隨著時(shí)間的發(fā)展，壁面開孔已成一種趨勢[14]，傳統(tǒng)的理論及模型已不再完全適用。雖然部分學(xué)者發(fā)現(xiàn)開孔后填料塔內(nèi)的傳質(zhì)及流動通量有所改變[15-17]，但國內(nèi)外文獻(xiàn)均未見對窗口區(qū)的液膜流動進(jìn)行報(bào)道。

此外，現(xiàn)有的工業(yè)設(shè)備，從填料片上幾毫米的小孔到大型圓盤反應(yīng)器數(shù)米的開窗，其形狀和尺寸都相差巨大，但在工業(yè)實(shí)踐中尚無明確統(tǒng)一的指導(dǎo)理論。目前的認(rèn)知表明，當(dāng)窗口尺寸小于毛細(xì)長度時(shí)整個(gè)孿生液膜的性能更接近壁面膜(孔內(nèi)會形成液橋和死區(qū)而不再體現(xiàn)自由液膜的流動特征)；當(dāng)窗口尺寸遠(yuǎn)大于毛細(xì)長度時(shí)，在重力加速作用下，窗內(nèi)液體流動更接近自由液膜(實(shí)際上通常對自由液膜的研究中，為了防止液膜收縮，兩側(cè)都有鐵棒等固體壁面作為支撐)。因而中間尺度的開窗恰恰能最好地體現(xiàn)出孿生液膜的流動特征，也是本文的重點(diǎn)研究對象。最近，Nikolai Kolev等[18]開發(fā)了一種由橫向條片組成的高性能填料(實(shí)則為開有若干矩形窗的平板)，實(shí)驗(yàn)結(jié)果表明此種填料有效比表面積大于1，十分適用于低液體流量下的吸收過程。同時(shí)期，Hu等[19-21]開發(fā)了一種適用于粘性吸收劑的垂直板填料。與前者不同，Hu等強(qiáng)調(diào)板上要開有足夠面積的窗口，使無壁支承的自由膜與壁面膜共存，且相互作用。實(shí)驗(yàn)結(jié)果表明，在一定操作范圍內(nèi)，此種填料相比于商業(yè)波紋填料傳質(zhì)效率更高，操作壓降更低，操作通量更大。值得注意的是這兩者的開窗尺寸都在厘米級，大于常見液體的毛細(xì)長度，但又小于通常自由液膜的研究尺寸。

本文在從事液膜相關(guān)工作基礎(chǔ)上，提出孿生液膜這一新概念。即壁面膜與受限自由膜同時(shí)存在、相互作用，形成時(shí)互為前驅(qū)[19]。此處所謂受限自由膜是指自由面為封閉界面限定，區(qū)別于已有廣泛研究的無界自由膜。這一概念的提出，既豐富了傳統(tǒng)的液膜流動概念，也因其本身已經(jīng)存在于工業(yè)實(shí)際，對它的理論研究將有助于這類工業(yè)設(shè)備的開發(fā)和效率提升，亦可進(jìn)一步為傳熱、傳質(zhì)設(shè)備提供強(qiáng)化手段。王良生等[22]在對聚酯開窗圓盤反應(yīng)器的研究中，較早注意到窗口區(qū)液膜的形成機(jī)理與壁面液膜完全不同。隨后在其基礎(chǔ)上，鄧斌等[23]通過數(shù)值模擬發(fā)現(xiàn)窗口區(qū)液膜的表面更新速率要比壁面液膜的高出近40%。但是兩者均未對窗口區(qū)液膜形成及其特征進(jìn)行詳細(xì)分析。Hu等[24]最早通過VOF方法模擬分析了二維開窗區(qū)受限自由液膜的部分流動特性，但亦未論述它的形成條件。本文將在前期工作的基礎(chǔ)上，展現(xiàn)窗口區(qū)豐富多樣的流型，并系統(tǒng)地分析物性和窗口形狀對窗口區(qū)成膜性的影響。本文結(jié)構(gòu)如下：首先闡述實(shí)驗(yàn)設(shè)備和觀測方法，其次考察窗口區(qū)的各種流型，最后標(biāo)度分析流動的特征量并給出窗口區(qū)成膜條件的經(jīng)驗(yàn)判別式。

1 實(shí)驗(yàn)設(shè)備與過程

圖1是測試裝置示意圖，主要由高位槽、液體流量計(jì)、水泵、液體分布器、測試板及液體收集槽、高速攝影儀等組成。實(shí)驗(yàn)在常壓下進(jìn)行，液體由泵輸送至高位槽，經(jīng)轉(zhuǎn)子流量計(jì)調(diào)節(jié)流量，在重力作用下從高位槽流至分布器。分布器出口為寬150mm的15°傾斜板，測試板經(jīng)夾具垂直固定在金屬架上。仔細(xì)調(diào)整兩者的間隙后，液體經(jīng)斜板導(dǎo)流能夠均勻鋪展于測試板上，最終回到下方水箱。為了準(zhǔn)確記錄實(shí)驗(yàn)過程中的流動現(xiàn)象，采用美國York-Tech公司的Phantom V 9.1高速攝影機(jī)，拍攝速率為100～1000幀/s，使用配套軟件Phantom Video Player對錄像進(jìn)行后處理。

圖1 實(shí)驗(yàn)裝置圖Fig.1 Experimental setup

測試板為200mm×200mm光滑鋁板，開窗區(qū)距離液體入口約80mm(見圖2)。為了方便敘述，與液膜流動方向垂直的窗邊定義為迎流邊(長度L)，與液膜流動方向平行的窗邊定義為流向邊(長度W)。為保證液體在窗上方呈均勻液膜流動，采用細(xì)砂紙對鋁板表面輕微打磨，并用丙酮溶液仔細(xì)擦洗以改善其潤濕性，同時(shí)在每次試驗(yàn)之前進(jìn)行預(yù)潤濕。采用日本三豐SJ-201表面粗糙度測試儀，測得打磨后鋁板表面的粗糙度為Ra=0.4μm，而本文研究的液膜膜厚均大于100μm且為層流工況，因此可以判定對液膜主體流動的影響可以忽略。

圖2 測試板Fig.2 Test plate

為了考察不同物性對窗口區(qū)成膜流量的影響，測試了表1所示6種液體。其中密度使用50ml容量瓶以稱重法獲得；粘度測量采用美國TA公司(TA Instruments Inc.New Castle.DE.USA)生產(chǎn)的ARES高級流變儀測得，測量精度為±0.1mPa·s；液體表面張力由上海中晨數(shù)字技術(shù)設(shè)備有限公司(Shanghai Zhongchen Technic Apparatus Co.Ltd)生產(chǎn)的JK99C型全自動張力儀測得，測量精度為±0.01mN/m。實(shí)驗(yàn)前后重復(fù)測量，整體物性參數(shù)的誤差約為5%。

表1 實(shí)驗(yàn)流體的物理性質(zhì)Table 1 Physical properties of the liquids used

2 實(shí)驗(yàn)結(jié)果

2.1窗口區(qū)的基本流型

觀察表明液體在窗口區(qū)的流動(見圖3)可分為:(1) 繞流：液體繞過窗口邊緣，窗下方為“干區(qū)”;(2) 偏流：少量液體進(jìn)入窗區(qū)兩角，大部分液體繞過邊緣，窗下仍為干區(qū);(3) 過流：液體越過窗上邊緣進(jìn)入窗區(qū)，并潤濕窗下方壁面;(4) 背流：液體偏流、過流時(shí)會有少量液體進(jìn)入板的另一側(cè)。

上述基本流型的形成主要取決于表面張力、慣性力、重力和黏性力間的平衡。窗口上游均勻液膜向下鋪展，在達(dá)到窗口上邊緣時(shí)，會在窗口上沿受到表面張力的阻滯作用。圖中，γsv、γsl、γlv分別為氣固、固液、氣液界面張力，θ為動態(tài)接觸角(前進(jìn)角)。固體壁面上液體的鋪展，γsv=γsl+γlvcosθ(見圖4(a))，與此不同，在開窗區(qū)，窗口上邊緣界面張力明顯構(gòu)成了一個(gè)與流動方向相反的力γsl+γlvcosθ(見圖4(b))，從而液體進(jìn)入到窗口區(qū)受阻，液體堆積。Oliver[25]等人通過實(shí)驗(yàn)對銳角邊緣的單液滴移動接觸線進(jìn)行了定量化的研究。結(jié)果表明，隨著液滴體積的增大，θ會不斷增大以維持重力與界面張力的穩(wěn)定，當(dāng)θ大于臨界值時(shí)(具體量值受到物性、接觸角和銳角大小等多重參數(shù)的影響)，接觸線才會突越過銳角邊緣。本文中，液體跨過開窗區(qū)同樣需要達(dá)到臨界能量或流量，來迫使液體進(jìn)入到窗區(qū)。

圖3 窗口區(qū)的基本流型示意圖Fig.3 Basic flow patterns around the window region

圖4 壁面(a)和窗口區(qū)(b)移動接觸線示意圖Fig.4 Sketch of moving contact line on the wall (a) and up the window (b)

當(dāng)慣性力和重力克服界面張力后，液體才能進(jìn)入到窗口區(qū)。由于不存在壁面的剪切和黏附作用，此時(shí)的流動和成膜機(jī)理更多地顯現(xiàn)出自由膜的特征。但又與通常自由膜不同：上下游壁面液膜和窗口內(nèi)的自由膜(假定已經(jīng)成膜，實(shí)際上成膜前有復(fù)雜的流型演變，下詳)存在流動的連續(xù)性和相互干擾；窗上下左右壁面限定自由膜，會對自由膜流動產(chǎn)生直接作用，如液體在窗口區(qū)經(jīng)重力作用加速后會撞擊到下壁面；開窗后，壁面兩側(cè)連通，增強(qiáng)了液體在膜厚方向上的流動。

對窗區(qū)表面張力作用初步分析可知：(1) 當(dāng)液體的慣性力不足以克服表面張力時(shí)，則會類似于局部駐點(diǎn)壓力形成繞流；(2) 少量液體進(jìn)入到窗區(qū)，由于接觸到左或右側(cè)壁面(壁面液膜)，在表面張力的作用下迅速被拉離窗口區(qū)，即偏流；(3) 上游液體進(jìn)入到窗口區(qū)后，在表面張力作用下急劇收縮，在接觸到下壁面后又會重新鋪展開，即過流；(4) 由于橫向γsv的作用，以及偏流、過流時(shí)液體在板兩側(cè)的重新分配，少量液體會進(jìn)入到板的另一側(cè)，即背流。

2.2過流流型演變

圖5 不同流量下窗口區(qū)的部分流型圖(介質(zhì)為水)Fig.5 Flow patterns in the window (L=48mm,W=12mm) with different Re (water)

(1) 單滴形成：當(dāng)流量特別小時(shí)，窗口上沿形成了一段厚的液層，類似于局部駐點(diǎn)壓力，上游液體會以繞流的形式從窗口區(qū)兩側(cè)流走。隨著流量的增大，液層不斷增厚。當(dāng)Re達(dá)到13.9時(shí)，可以發(fā)現(xiàn)窗口上沿中部的厚液層開始進(jìn)入到窗口區(qū)，并在重力作用下不斷積累下落，最終帶動周圍液體以液橋的形式通過窗口區(qū)并潤濕窗下游的壁面。隨后，液橋在表面張力的拉伸下迅速斷裂，回復(fù)到液層積累的初始狀態(tài)，形成一個(gè)周期性的過流過程。

(2) 滴柱共存：隨著流量的逐漸增大，液滴下落的頻率加快，并逐漸轉(zhuǎn)變?yōu)榉€(wěn)定的液柱。與此同時(shí)迎流邊上的液滴落點(diǎn)也將會不斷增加，類似于圓環(huán)和圓柱上的液柱[26]，窗口區(qū)的落點(diǎn)之間也存在著臨界距離，此處約為10mm。隨著落點(diǎn)數(shù)量的飽和，窗口區(qū)的流動開始出現(xiàn)不規(guī)則和隨機(jī)性。由于開窗區(qū)不存在強(qiáng)的黏性剪切，表面張力成為控制液體流動的主導(dǎo)作用。液柱在輕微擾動下就會相互碰撞和聚并，靠近窗兩側(cè)的自由膜還能直接以偏流的形式被壁面液膜帶走。當(dāng)Re達(dá)到111時(shí)，窗口區(qū)的流動又重新達(dá)到一個(gè)穩(wěn)定的狀態(tài)，液體以4條液柱的形式穿過窗口區(qū)。

(3) 掛膜階段：進(jìn)一步增大流量后Re=125，液柱聚并成股，窗口區(qū)的流動重新變得不規(guī)則。值得注意的是上游液膜的前鋒不再限于壁面厚液層區(qū)，而是掛在了窗口上沿下方，此時(shí)慣性力加上重力與表面張力在窗口區(qū)達(dá)成了微妙的平衡。由于水具有極高的Ka數(shù)，在自然擾動下就能形成小的毛細(xì)波。實(shí)驗(yàn)過程中發(fā)現(xiàn)當(dāng)上游毛細(xì)波達(dá)到窗口區(qū)時(shí)，慣性力就能暫時(shí)性地克服表面張力的約束，大量的液體以片狀形式進(jìn)入到窗口區(qū)。由于此時(shí)流量尚低，液體無法鋪滿窗口區(qū)，在兩側(cè)和下方局部表面張力的作用下，不完整的液膜很快就斷裂消失。這一階段是窗口成膜的必經(jīng)階段。

(4) 成膜階段：之后隨著流量的增加，進(jìn)入到窗口的前驅(qū)液膜越來越逼近下壁面(如Re=153，窗口區(qū)的前驅(qū)液膜甚至能覆蓋窗口的大半)，并在重力、慣性力與表面張力的共同作用下快速收縮和拉伸。直到流量繼續(xù)增大，或者某一瞬間的擾動使得前驅(qū)液膜完全與下壁面貼合，從而形成一個(gè)穩(wěn)定的力平衡。如圖所示，Re=158完整的受限自由膜會展現(xiàn)出一個(gè)相對穩(wěn)定并十分規(guī)整的波列結(jié)構(gòu)。

通過上述觀察，可以發(fā)現(xiàn)窗口區(qū)液膜流動是在慣性力、表面張力、粘性力和重力共同作用的結(jié)果。只有在足夠大流量下，慣性力和重力才能克服液體的表面張力進(jìn)入到窗口區(qū)。與此同時(shí)，粘性力則能在一定程度上抑制擾動的傳播和自由液體表面的不規(guī)則運(yùn)動，促進(jìn)液體呈膜狀接觸到下壁面，最終形成完整的液膜。值得注意的是，即使沒有成膜，液體在流經(jīng)窗口區(qū)的同時(shí)也將發(fā)生自由表面的收縮與接觸到下壁面時(shí)的再鋪展，這一過程對工業(yè)上的傳質(zhì)及混合都將十分有利。

2.3標(biāo)度分析

在進(jìn)一步討論之前，先對流動的特征量進(jìn)行簡單的標(biāo)度分析。以液相特征速度Vl和特征長度Ll為基礎(chǔ)，導(dǎo)出慣性力Fli、粘性力Flv、重力Flg和表面張力Fls的表達(dá)式:

再將式(2)～(4)分別除以式(1)可以得到如下無量綱式：

從上述分析可以發(fā)現(xiàn)影響整個(gè)流動工況的主要是Re和Ka這2個(gè)無量綱數(shù)，而表面張力的作用最為明顯(實(shí)驗(yàn)范圍Ka：52.7～2.93×103，臨界成膜Re：4～140)。式(7)～(9)也在一定程度上解釋了為什么Ka，KF=Ka-3[27]以及Ga=Ka3[28]會經(jīng)常用來關(guān)聯(lián)流型轉(zhuǎn)變的Re。同時(shí)為了表征窗幾何形狀對成膜流量的影響引入無量綱數(shù)Nx(L，W)，得到了如下形式的經(jīng)驗(yàn)關(guān)聯(lián)式：

式中：C、n、b為待定常數(shù)。

2.4成膜條件關(guān)聯(lián)式擬合

在液膜研究領(lǐng)域中，大量文獻(xiàn)關(guān)注液膜的流動規(guī)律，著眼于液膜形成的很少。本文重點(diǎn)研究了受限自由膜的形成條件。圖6展示了不同矩形窗和不同物性條件下的窗口區(qū)臨界成膜流量?？梢园l(fā)現(xiàn)，隨著迎流邊長、流向邊長以及Ka數(shù)的增加，臨界成膜Re不斷增加。經(jīng)擬合得到經(jīng)驗(yàn)關(guān)聯(lián)式:

或者

圖6 不同矩形窗和物性下的成膜Re數(shù)Fig.6 Critical film formation Re with different windows and fluids

上述關(guān)聯(lián)式在一定程度上反應(yīng)了成膜過程主要是液膜慣性力和重力克服表面張力向下鋪展的過程，而窗迎流長度和流向長度的增加會延長這一鋪展過程同時(shí)降低成膜過程的穩(wěn)定性，黏性力則在一定程度上抑制液膜的不穩(wěn)定性，有利于成膜。因此成膜流量隨液體表面張力系數(shù)和窗尺寸的增大而增大，隨密度和黏度的增高而減小。當(dāng)然由于成膜過程十分復(fù)雜，迎流邊(L)與流向邊(W)發(fā)揮著各自不同的作用，上式中窗幾何參數(shù)Nx的模型過于簡化，因此還存在一定局限性。

3 結(jié) 論

本文通過實(shí)驗(yàn)觀測對垂直板開窗區(qū)的液膜流動行為及其成膜性進(jìn)行了系統(tǒng)的研究。得到結(jié)論：

(1) 下降液膜在通過開窗區(qū)時(shí)有繞流、過流、偏流和“背流”4種基本流型。對于過流，隨著流量的變化，窗口內(nèi)會出現(xiàn)單液滴、多液滴、單一或多重液柱、液膜以及它們的組合，即多種型態(tài)共存、破裂、轉(zhuǎn)變。顯然，窗口區(qū)受限自由膜的流型具有多樣性和動態(tài)性的特征。

(2) 本文研究了不同物性、不同矩形窗條件下的受限自由膜成膜流量，并經(jīng)過標(biāo)度分析擬合得到了成膜條件的經(jīng)驗(yàn)判別式。結(jié)果表明，臨界成膜Re數(shù)隨著窗迎流邊長、流向邊長和Ka數(shù)的增大而增大，Ka數(shù)更為敏感。

(3) 成膜后，開窗板上會呈現(xiàn)壁面液膜與窗內(nèi)受限自由液膜交替共存(孿生液膜)的特殊流動工況。窗口區(qū)受限自由膜特殊的波列結(jié)構(gòu)，開窗后板兩側(cè)流體的連通，流體在窗區(qū)的重力加速、對窗下壁面的撞擊，特別是壁面膜與自由膜的相互作用(即流量波動或擾動下，多種流型的收縮破裂、動態(tài)演變)等導(dǎo)致的多重干擾，將會強(qiáng)化全板上的液膜表面更新與液膜波動和湍動。

[1]胡軍,胡國輝,孫德軍.沿平板下落薄膜流動的研究綜述[J].力學(xué)進(jìn)展,2005,35(2):161-169.Hu J,Hu G H,Sun D J.A review on thin films falling along an inclined plate[J].Advances in Mechanics,2005,35(2):161-169.

[2]Batchelor G K,Moffatt H K,Worster M.Perspectives in fluid dynamics:a collective introduction to current research[M].Cambridge:Cambridge University Press,2002.

[3]Kapitza P.Dynamic stability of the pendulum with vibrating suspension point[J].Soviet Physics-JETP,1951,21(5):588-597.

[4]Wasden F K,Dukler A E.A numerical study of mass transfer in free falling wavy films[J].American Institute of Chemical Engineers Journal,1990,36(9):1379-1390.

[5]Alekseenko S,Nakoryakov V E,Pokusaev B G.Wave flow of liquid films[M].New York:Begell House,1994.

[6]Chang H-H,Demekhin E A.Complex wave dynamics on thin films[M].Elsevier,2002.

[7]Kalliadasis S,Ruyer-Quil C,Scheid B,et al.Falling liquid films[M].Springer Science &Business Media,2011.

[8]Squire H B.Investigation of the instability of a moving liquid film[J].British Journal of Applied Physics,1953,4(6):167-169.

[9]Brown D.A study of the behaviour of a thin sheet of moving liquid[J].Journal of Fluid Mechanics,1961,10(2):297-305.

[10]Lin S,Roberts G.Waves in a viscous liquid curtain[J].Journal of Fluid Mechanics,1981,112:443-458.

[11]Lin S P.Breakup of liquid sheets and jets[M].England:Cambridge University Press Cambridge,England,2003.

[12]Sirignano W A,Mehring C.Review of theory of distortion and disintegration of liquid streams[J].Progress in Energy &Combustion Science,2000,26(4):609-655.

[13]Pritchard W.Instability and chaotic behaviour in a free-surface flow[J].Journal of Fluid Mechanics,1986,165:1-60.

[14]Mackowiak J.Fluid dynamics of packed columns[M].London,New York:Springer,2010.

[15]Fair J R,Seibert A F,Behrens M,et al.Structured packing performance experimental evaluation of two predictive models[J].Industrial &Engineering Chemistry Research,2000,39(6):1788-1796.

[16]Pavlenko A,Pecherkin N,Chekhovich V,et al.Hydrodynamics in falling liquid films on surfaces with complex geometry[J].Microgravity Science and Technology,2009,21(1):207-213.

[17]Yao Y,Pavlenko A,Volodin O.Effects of layers and holes on performance of wire mesh packing[J].Journal of Engineering Thermophysics,2015,24(3):222-236.

[18]Kolev N,Kralev B,Kolev D.Gas side controlled mass transfer in a new packing with stamped horizontal lamellae operating at extremely low liquid loads[J].Chemical Engineering and Processing:Process Intensification,2013,63:44-49.

[19]戴干策,胡劍光,于建國,等.適應(yīng)粘性吸收劑的有壁與無壁液膜交替的規(guī)整填料:中國,201210234334.6[P].2013-07-10.Dai G C,Hu J G,Yu J G,et al.Wall liquid film and wall-free liquid film alternate structured filler adapting to viscous absorbent:China，201210234334.6[P].2013-07-10.

[20]Hu J,Liu J,Yu J,et al.CO2absorption into highly concentrated dea solution flowing over a vertical plate with rectangular windows[J].International Journal of Greenhouse Gas Control,2013,19:13-18.

[21]胡劍光,劉佳特,袁猛,等.新型垂直板規(guī)整填料流體力學(xué)及傳質(zhì)性能[J].化工學(xué)報(bào),2014,65(1):116-122.Hu J G,Liu J T,Yuan M,et al.Hydrodynamics and mass transfer characteristics of a novel vertical-sheet structured packing[J].Journal of Chemical Industry and Engineering,2014,65(1):116-122.

[22]王良生,戴干策.圓盤反應(yīng)器成膜性的持液量研究[J].化學(xué)反應(yīng)工程與工藝,2000,16(2):127-135.Wang L S,Dai G C.A study of film forming and hold-up in a rotating disc-ring reactor[J].Chemical Reaction Engineering and Technology,2000,16(2):127-135.

[23]鄧斌,戴干策.圓盤反應(yīng)器液膜表面更新數(shù)值模擬[J].化工學(xué)報(bào),2015,66(4):1407-1416.Deng B,Dai G C.Numerical simulation of surface renewal frequency on vertically rotating disc[J].Journal of Chemical Industry and Engineering,2015,66(4):1407-1416.

[24]Hu J,Yang X,Dai G.Numerical investigation on hydrodynamics of vertically confined free film[J].The Canadian Journal of Chemical Engineering,2016,94(2):340-348.

[25]Oliver J F,Huh C,Mason S G.Resistance to spreading of liquids by sharp edges[J].Journal of Colloid and Interface Science,1977,59(3):568-581.

[26]Brunet P,Flesselles J-M,Limat L.Dynamics of a circular array of liquid columns[J].The European Physical Journal B,2007,55(3):297-322.

[27]Schmuki P,Laso M.On the stability of rivulet flow[J].Journal of Fluid Mechanics,1990,215:125-143.

[28]Ruan B,Jacobi A M,Li L.Effects of a countercurrent gas flow on falling-film mode transitions between horizontal tubes[J].Experimental Thermal and Fluid Science,2009,33(8):1216-1225.

謝涵廣(1991-)，男，浙江臺州人，博士研究生。研究方向：液膜流動。通信地址：上海市徐匯區(qū)梅隴路130號 (200237)。E-mail:xhgshangwu@163.com

(編輯：張巧蕓)

Formationofconfinedfreefilminthewindowonaverticalperforatedplate

Xie Hanguang,Hu Jianguang,Wang Cheng,Dai Gance*

(State Key Laboratory of Chemical Engineering,East China University of Science and Technology,Shanghai 200237,China)

Liquid film is a common contact way between gas/liquid and can be divided into wall-bounded film and free film according to the film formation process and the number of free surfaces.When a plate is perforated,both kinds of films exist.This paper describes some experimental observations of free-surface flows arising when a thin liquid film flows through the window on a vertical perforated plate.Dozens of rectangular windows (9～1152mm2) and six fluids (Kafrom 52 to 3000) were used to investigate the flow mechanisms.Several typical flow patterns including pass-around flow,pass-through flow,bias flow and back-side flow are well defined under such a flow condition.Here a special focus is given to the window region,where various free-surface flow patterns composed of droplets,columns,sheets and their combinations were observed with increasing flow rate.Meantime,free film in the window is surface tension dominated and susceptible to disturbance,resulting in coalescence or break of liquid columns and films.At a critical flow rate,liquid film is able to full fill the window,forming a stable complete confined free film.Mutual influence between confined free film in the window region and wall-bounded film around shows special wavy trains,which is also called “twin liquid film”.Based on experimental data and scaling analysis,an empirical equation which relates Reynolds numberRe,Kapitza numberKaand a dimensionless lengthNxis proposed to characterize the film formation conditions.It is found that critical film formation Reynolds number increases with Kapitz number and window size.And hysteresis phenomenon is manifested by obviously different flow transitionRefor confined free film formation and breaking.The results can help window geometry optimization industrial processes to improve local heat and mass transfer.It can also enrich the traditional film flow investigations.

twin-liquid film;flow transition;film formation;perforated plate

O363.2

A

1672-9897(2017)05-0032-07

10.11729/syltlx20160197

2016-12-15;

2017-08-23

*通信作者 E-mail:gcdai@ecust.edu.cn

XieHG,HuJG,WangC,etal.Formationofconfinedfreefilminthewindowonaverticalperforatedplate.JournalofExperimentsinFluidMechanics,2017,31(5):32-38.謝涵廣,胡劍光,王 成,等.垂直板上開窗區(qū)受限自由膜的形成.實(shí)驗(yàn)流體力學(xué),2017,31(5):32-38.