流動(dòng)液膜的數(shù)值模擬研究綜述

吳正人，宋朝匣，劉梅，劉秋升

?

流動(dòng)液膜的數(shù)值模擬研究綜述

吳正人1，宋朝匣1，劉梅2，劉秋升1

（1華北電力大學(xué)能源動(dòng)力與機(jī)械工程學(xué)院，河北保定 071003；2華北電力大學(xué)經(jīng)濟(jì)管理系，河北保定 071003）

液膜流動(dòng)現(xiàn)象廣泛存在于自然界中，作為一種高效傳熱傳質(zhì)技術(shù)，其在化工等領(lǐng)域有著廣泛的應(yīng)用。近幾年來(lái)，國(guó)內(nèi)外學(xué)者越來(lái)越熱衷于運(yùn)用數(shù)值模擬技術(shù)來(lái)研究液膜的流動(dòng)特性及傳熱傳質(zhì)特性。本文歸納分析了數(shù)值模擬研究中液膜自由液面的追蹤方法?？偨Y(jié)了不同壁面結(jié)構(gòu)、不同壁面傾角、液體物性、液相流量與氣相流速4個(gè)方面對(duì)液膜的流動(dòng)特性的影響規(guī)律，以及改變壁面傾角、入口雷諾數(shù)、入口添加擾動(dòng)時(shí)表面波呈現(xiàn)的波動(dòng)特性。此外，還論述了流動(dòng)液膜的傳熱傳質(zhì)特性的研究現(xiàn)狀。所得結(jié)論對(duì)流動(dòng)液膜的數(shù)值模擬研究具有一定的參考價(jià)值，最后提出了用數(shù)值模擬方法研究液膜流動(dòng)的缺陷與不足，展望了更加科學(xué)合理地研究流動(dòng)液膜的方法。

液膜；流體動(dòng)力學(xué)；兩相流；傳熱傳質(zhì)；數(shù)值模擬

在許多化工設(shè)備中，如板式蒸發(fā)冷凝器、規(guī)整填料片、膜式氣液反應(yīng)器等，都存在液膜流動(dòng)，這種膜式流動(dòng)包括不相溶的氣-液兩相流和液-液兩相流。掌握液膜的流動(dòng)特性和傳熱傳質(zhì)特性，對(duì)于分析這些設(shè)備中的工藝流程以及提高化工設(shè)備的運(yùn)行性能意義重大。

20世紀(jì)中期，Kapitza父子進(jìn)行了液體薄膜沿傾斜壁面下落的開(kāi)創(chuàng)性實(shí)驗(yàn)。至此以后的半個(gè)多世紀(jì)，對(duì)薄膜流動(dòng)的動(dòng)力學(xué)性能的研究引起了許多科研工作者的關(guān)注。此外，早期對(duì)于液膜流動(dòng)的研究主要集中在平整基底上，隨著研究的進(jìn)展和工程生產(chǎn)中的需要，對(duì)液膜流動(dòng)的研究就不僅僅局限于平整壁面。近年來(lái)，對(duì)非平整壁面或管內(nèi)、外液膜的流動(dòng)特性及傳熱傳質(zhì)特性的研究成為熱點(diǎn)。

對(duì)液膜流動(dòng)的研究主要包括理論分析、實(shí)驗(yàn)研究和數(shù)值模擬等方法，本文主要對(duì)基于數(shù)值模擬技術(shù)研究液膜流動(dòng)特性的現(xiàn)狀進(jìn)行綜述分析。

1 液膜自由界面追蹤方法

對(duì)液膜流動(dòng)進(jìn)行數(shù)值模擬時(shí)，由于液膜自由液面隨時(shí)間變化很難準(zhǔn)確追蹤捕捉到，從而使得數(shù)值模擬存在很大困難。有關(guān)學(xué)者經(jīng)過(guò)研究提出了幾種自由邊界的處理方法，包括液膜厚度方程法（height function）、線段法（line segments）、標(biāo)記粒子法（markerparticles）和VOF方法（volume of fluid）。

液膜厚度方程法是一種只能用來(lái)確定一維自由邊界面的方法，而線段法打破了這種方法的局限性，但仍然很難擴(kuò)展到三維面上。標(biāo)記粒子法雖然可以擴(kuò)展到三維問(wèn)題，但需要很大的存儲(chǔ)空間。在這4種方法中最理想、應(yīng)用最廣的是VOF方法，其不僅可以研究三維問(wèn)題，而且只用一個(gè)體積分?jǐn)?shù)函數(shù)就可以來(lái)描述流體區(qū)域中各相的復(fù)雜變化，能比較準(zhǔn)確地跟蹤捕捉到自由界面，已經(jīng)引起了越來(lái)越多計(jì)算流體研究人員的重視[1-8]。

2 液膜流動(dòng)特性

影響液膜流動(dòng)特性的因素很多，本文主要從不同壁面結(jié)構(gòu)、不同壁面傾角、液體物性、液相流量與氣相流速4個(gè)方面進(jìn)行分析總結(jié)。

2.1 不同壁面結(jié)構(gòu)的影響

谷芳等[9-11]建立了二維和三維氣-液兩相流的CFD模型，分析研究了不同板面結(jié)構(gòu)下液膜的流動(dòng)特性。李超[12]建立了二維氣-液兩相流模型，采用了谷芳等提出的表面張力動(dòng)量源項(xiàng)和氣液剪切應(yīng)力動(dòng)量源項(xiàng)，研究了板面結(jié)構(gòu)等重要參數(shù)的影響。吳治將等[13]、孫荷靜[14]分別建立了波紋板蒸發(fā)式冷凝器的二維和三維計(jì)算模型，研究了板面結(jié)構(gòu)對(duì)液膜流動(dòng)的影響，并分析了它們對(duì)波紋板蒸發(fā)式冷凝器的傳熱影響機(jī)理。劉鑫[15]分別對(duì)傾斜波形壁面和矩形結(jié)構(gòu)壁面上的液膜流動(dòng)進(jìn)行了數(shù)值模擬，探究了不同壁面結(jié)構(gòu)對(duì)液膜流動(dòng)特性的影響，并且通過(guò)改變不同結(jié)構(gòu)壁面上流體流量以及流體物性來(lái)分析液膜成膜質(zhì)量的影響因素。

由以上研究，通過(guò)對(duì)比平板和波紋板上的流動(dòng)可知，液膜在長(zhǎng)度相同的波紋板和平板上流動(dòng)時(shí)，平板的液膜流動(dòng)時(shí)間比波紋板上的短，這是由于波紋結(jié)構(gòu)使液膜的流動(dòng)路程增大，從而使液膜在板面上的停留時(shí)間變長(zhǎng)，這樣為傳熱傳質(zhì)提供了更有利的條件。

通過(guò)對(duì)不同板面結(jié)構(gòu)對(duì)液膜流動(dòng)的影響研究發(fā)現(xiàn)，波紋特征尺寸越大，液體越容易在波紋谷底形成漩渦，從而造成液體滯留，不利于板面上流體的更新。并且液膜在波紋尺寸較小的板面上容易發(fā)生斷裂，使板面上形成“干斑”，而在波紋尺寸較大的板面上容易形成連續(xù)液膜。另外，研究者還測(cè)定了不同壁面結(jié)構(gòu)下剪切應(yīng)力的分布，并分析得出近壁區(qū)漩渦的位置和大小可通過(guò)測(cè)定壁面剪切應(yīng)力來(lái)預(yù)測(cè)。

2.2 不同壁面傾角的影響

液膜厚度、壁面剪切應(yīng)力等都會(huì)隨著板面傾角變化而有一定的變化。研究結(jié)果表明，隨波紋板傾角增加，液膜厚度逐漸增加，而壁面剪切應(yīng)力呈逐漸減小趨勢(shì)，但板面傾角的變化不會(huì)改變液膜的基本流動(dòng)結(jié)構(gòu)，即液膜自由液面與板面的相位差。同時(shí)研究發(fā)現(xiàn)，隨著傾角的逐漸增大，波紋板進(jìn)出口壓降先增大后減小再增大。

2.3 液體物性的影響

液體物性對(duì)液膜的流動(dòng)狀態(tài)穩(wěn)定性以及傳熱傳質(zhì)特性有很大影響。谷芳等[10-11]對(duì)乙醇、丙酮、甘油1、甘油2和水5種液體的黏度、表面張力和接觸角等物性進(jìn)行了數(shù)值模擬。孫鳳丹[16]利用Fluent軟件對(duì)二維波紋板上的液膜進(jìn)行研究，探究了液相物性等對(duì)液膜流動(dòng)的影響。劉梅等[17]同樣對(duì)液體物性對(duì)液膜流動(dòng)特性的影響進(jìn)行了數(shù)值模擬研究。以上研究發(fā)現(xiàn)，在相同液體進(jìn)口速度下，黏度大的液體其流動(dòng)形成的液膜厚度也大，但不同流體的流動(dòng)結(jié)構(gòu)不會(huì)發(fā)生變化。另外發(fā)現(xiàn)，對(duì)于非常薄的液膜流動(dòng)，表面張力起到至關(guān)重要的作用，當(dāng)不考慮表面張力時(shí)，流體在波紋壁面上能形成連續(xù)液膜，且液膜自由表面形狀、壁面剪切應(yīng)力緊隨波紋壁面變化；當(dāng)考慮表面張力時(shí)，液膜自由表面波形與固體表面波形出現(xiàn)了很大的相位差；此外，由于水具有很大的表面張力，在固體壁面上不能形成連續(xù)液膜。壁面剪切應(yīng)力增大出現(xiàn)較大波動(dòng)，此時(shí)壁面剪切應(yīng)力符號(hào)改變，則預(yù)示著此時(shí)漩渦產(chǎn)生，從而對(duì)傳熱傳質(zhì)產(chǎn)生一定的影響。因此，在模擬計(jì)算壁面上的液膜流動(dòng)時(shí)，表面張力的影響不能忽略。

2.4 液相流量與氣相流速的影響

液相速度和氣相速度也是影響液膜流動(dòng)的重要因素。Ho等[7]采用VOF模型對(duì)降膜反應(yīng)器中的兩相流動(dòng)進(jìn)行CFD模擬；李相鵬等[18]基于精餾過(guò)程，對(duì)規(guī)整填料表面液膜的流動(dòng)特性進(jìn)行了模擬分析；吳治將、孫荷靜等[13-14]分別建立了波紋板蒸發(fā)式冷凝器的二維和三維計(jì)算模型；Xu等[19]基于Fluent軟件建立模型，研究了氣液相速度對(duì)液膜厚度分布、液相平均流速等的影響；孫鳳丹[16]同樣探究了液體流速對(duì)液膜流動(dòng)的影響。

綜合以上研究結(jié)果表明，液相流量增大時(shí)，液膜厚度增大，液相流速也隨之增大，但并不是液相流速越大越容易形成連續(xù)液膜，連續(xù)液膜的形成與液相物性有很大關(guān)系。此外，由于氣體的流動(dòng)，使液膜自由表面發(fā)生波動(dòng)，從而不能形成穩(wěn)定的膜狀流動(dòng)。氣液兩相的流向也影響液膜流動(dòng)：氣液并流時(shí)，氣相速度的增大會(huì)減薄液膜厚度；氣液逆流時(shí)，隨著氣相速度的增大，液膜被吹散成液滴，湍動(dòng)程度加劇，液體入口處液膜厚度增加，出口處液膜厚度減小，且流動(dòng)明顯變得緩慢，對(duì)傳熱傳質(zhì)產(chǎn)生更大的影響。

3 表面波特性

由于液膜表面波高度影響著液膜內(nèi)的傳熱傳質(zhì)，因此對(duì)其的關(guān)注具有深刻意義。

3.1 入口擾動(dòng)的影響

液膜表面波受入口擾動(dòng)頻率及壁面結(jié)構(gòu)的共同影響。當(dāng)入口擾動(dòng)頻率很小時(shí)，低頻率擾動(dòng)發(fā)展形成孤立波；此外，表面波的振幅和波長(zhǎng)會(huì)隨著入口擾動(dòng)頻率的增大而減小。當(dāng)入口擾動(dòng)頻率較高時(shí)，表面波會(huì)受到很強(qiáng)的干擾，甚至發(fā)生合并現(xiàn)象。當(dāng)表面波的波長(zhǎng)和基底結(jié)構(gòu)波長(zhǎng)一致時(shí)，就會(huì)發(fā)生共振現(xiàn)象。共振是沿波紋壁面下落的液膜流動(dòng)的獨(dú)特的物理現(xiàn)象。

馬學(xué)虎、薄守石等[3-4]在考慮表面張力的作用下，基于VOF方法研究了在液膜入口處添加擾動(dòng)時(shí)的液膜表面穩(wěn)態(tài)波形，分析了擾動(dòng)頻率等對(duì)液膜表面波動(dòng)演化發(fā)展的影響及波內(nèi)速度場(chǎng)和壁面剪切應(yīng)力等的變化規(guī)律。樸明日、胡國(guó)輝[20]同樣在入口處引入周期擾動(dòng)，采用 Open FOAM 軟件對(duì)沿矩形波紋底板的二維液膜流動(dòng)進(jìn)行數(shù)值模擬，研究了壁面結(jié)構(gòu)對(duì)液膜流動(dòng)結(jié)構(gòu)的影響，并對(duì)流場(chǎng)中的共振現(xiàn)象進(jìn)行了分析。Kunugi等[21-22]為研究液膜表面波特性與傳熱特性間的關(guān)系，分別在有速度擾動(dòng)和無(wú)速度擾動(dòng)的情況下對(duì)垂直壁面上的下降液膜流動(dòng)進(jìn)行了二維數(shù)值模擬，得到了在不同入口擾動(dòng)頻率下的液膜厚度、流線和傳熱系數(shù)的分布。Miyara[23]在入口處加入小振幅周期速度擾動(dòng)，分析了波頻對(duì)垂直板面及微傾斜板面的波形、流線和傳熱特性的影響并進(jìn)行數(shù)值模擬研究。Valluri等[24]用CFX軟件在≈30下對(duì)雙正弦曲線波紋壁面下的液膜的動(dòng)態(tài)發(fā)展過(guò)程進(jìn)行研究，考察板面微結(jié)構(gòu)及擾動(dòng)頻率對(duì)液膜厚度及速度場(chǎng)的影響，驗(yàn)證了他們提出的積分平衡法的正確性。Min等[6]對(duì)整個(gè)流場(chǎng)的波的傳播、運(yùn)動(dòng)和加強(qiáng)進(jìn)行了數(shù)值研究，用VOF和PLIC（分段線性截面計(jì)算方法）得到了雷諾數(shù)在200～1000范圍內(nèi)的層流波狀流動(dòng)，包括平均液膜厚度、波的振幅、頻率和流速的分布。

研究發(fā)現(xiàn)，在入口邊界處的小振幅擾動(dòng)會(huì)發(fā)展成包含一個(gè)大振幅“滾波”和小振幅毛細(xì)波的孤立波。沿著流動(dòng)方向在滾波內(nèi)產(chǎn)生渦流，渦流的波幅會(huì)發(fā)展成與滾波波幅大小相近，此外，渦流對(duì)于強(qiáng)化傳熱傳質(zhì)有重要作用。另一方面，在毛細(xì)波內(nèi)沒(méi)有渦流產(chǎn)生。高頻的擾動(dòng)相互作用，滾波的振幅減小，毛細(xì)波伴隨著滾波消失，同時(shí)在已經(jīng)發(fā)展成的波內(nèi)會(huì)有渦流產(chǎn)生，在一個(gè)特定頻率下，它會(huì)由于波的相互作用重復(fù)產(chǎn)生和消失。

3.2 雷諾數(shù)（Re）的影響

目前，有關(guān)液膜流動(dòng)的研究中涉及了有關(guān)準(zhǔn)數(shù)的影響作用，如等。其中對(duì)于波動(dòng)液膜的發(fā)展有重要影響。本文對(duì)數(shù)的影響作用進(jìn)行了總結(jié)研究。

研究表明，數(shù)增大時(shí)，平均液膜厚度變大，波幅變小，孤立波逐漸消失，小振幅的波逐漸充滿整個(gè)流場(chǎng)；另外，數(shù)增大，沿著下游方向波幅增大，上游的波合并在一起并擴(kuò)大，沒(méi)有發(fā)生衰減，最后在下游處產(chǎn)生了高振幅的孤立波。

許松林等[25]數(shù)值模擬了有無(wú)氣相并流時(shí)，隨液相數(shù)的增大，垂直下落液膜的流動(dòng)特性及實(shí)時(shí)監(jiān)測(cè)了壁面剪切力分布特點(diǎn)。劉玉峰等[26]對(duì)于高數(shù)豎壁降膜流場(chǎng)進(jìn)行了數(shù)值模擬，分析了降水膜速度分布、溫度分布、波速及平均水膜厚度等特性參數(shù)的變化規(guī)律。模擬結(jié)果發(fā)現(xiàn)，隨著流動(dòng)向下游發(fā)展，平均水膜厚度先是迅速減小，然后減小速度趨緩，最后趨于恒定值；當(dāng)初始水膜厚度一定時(shí)，隨著初始的增大，平均水膜厚度逐漸增大，水膜波動(dòng)的振幅逐漸減小，但水膜穩(wěn)定性變差。Szulczewska等[27]建立了二維兩相逆流的板狀規(guī)整填料的CFD模型，研究液相或氣相流速和流動(dòng)流體的物理化學(xué)性質(zhì)對(duì)表面波區(qū)域的影響。另外，該模型可以用來(lái)確定在填料表面觀察到斷裂液膜時(shí)的最小液相流速。于意奇等[28]對(duì)傾斜壁面上的降膜流動(dòng)進(jìn)行了數(shù)值模擬，得到了不同下二維和三維的模擬結(jié)果。結(jié)果顯示，三維模擬能更好地表現(xiàn)界面波動(dòng)。孫文倩[29]對(duì)水平管外降膜流動(dòng)進(jìn)行了數(shù)值模擬，得到管間流型隨的變化而發(fā)生變化，且液膜厚度及液膜流速隨的增大而增大。

3.3 壁面傾角的影響

當(dāng)壁面傾斜角發(fā)生變化時(shí)，表面波也會(huì)發(fā)生變化，隨著傾斜角降低，波動(dòng)發(fā)展緩慢，波峰處的循環(huán)流動(dòng)及波谷處的渦流均逐漸消失。馬學(xué)虎[3-4]、Miyara[23]等對(duì)此進(jìn)行了研究分析。

4 傳熱傳質(zhì)特性

對(duì)于氣液界面?zhèn)鳠醾髻|(zhì)的研究，最早主要集中于理論和實(shí)驗(yàn)方面。近年來(lái)，采用數(shù)值模擬方法研究液膜的傳熱傳質(zhì)特性成為熱點(diǎn)。

郭常青等[30]將板式蒸發(fā)式冷凝器氣-液兩相流問(wèn)題簡(jiǎn)化為無(wú)限大平板上的氣液兩相流動(dòng)，重點(diǎn)考察了壁面熱流密度、水的進(jìn)口溫度和空氣流速對(duì)豎直板面溫度分布、氣-液兩相界面處潛熱和顯熱傳熱熱量密度的相對(duì)關(guān)系的影響。張景衛(wèi)、朱冬生等[31-32]將數(shù)值模擬與實(shí)驗(yàn)研究相結(jié)合，對(duì)板式蒸發(fā)式冷凝器的傳熱特性進(jìn)行了研究，最后得出空氣與水并流時(shí)的傳熱效果比逆流和錯(cuò)流時(shí)要好。當(dāng)空氣與水并流時(shí)，水膜厚度變薄，傳熱熱阻減小，從而更好地強(qiáng)化傳熱。孫荷靜等[33]建立了豎直波紋管和豎直圓管兩種物理模型，并對(duì)兩種模型內(nèi)的降膜流動(dòng)進(jìn)行了數(shù)值模擬研究，分析得出，在一定范圍內(nèi)，噴淋密度越大，傳熱性能越好，且波紋管比圓管傳熱性能好。周文生等[34]采用P-S（Patankar-Spalding）方法對(duì)蒸發(fā)器內(nèi)部的液膜流動(dòng)進(jìn)行數(shù)值模擬，并擬合出了換熱系數(shù)計(jì)算關(guān)系式，且誤差與之前研究結(jié)果相比相差不大。

界面波的強(qiáng)化傳熱機(jī)理可以通過(guò)顯示在不同下的溫度場(chǎng)分布來(lái)說(shuō)明。Miyara[23]在流動(dòng)入口處添加小振幅正弦擾動(dòng)，數(shù)值模擬了豎直平板上帶有界面波的下降液膜的流動(dòng)特性和傳熱特性。研究結(jié)果表明：界面波通過(guò)兩種效應(yīng)來(lái)強(qiáng)化傳熱，即削薄液膜和對(duì)流效應(yīng)。對(duì)于較低的數(shù)，可以忽略對(duì)流效應(yīng)；對(duì)于較高的數(shù)，這兩種效應(yīng)都起作用。另外，時(shí)均傳熱系數(shù)沿流動(dòng)方向逐漸增大。在完全發(fā)展區(qū)域，較高數(shù)下的傳熱系數(shù)要更高一點(diǎn)。然而，在發(fā)展區(qū)域，最大的數(shù)有最低的傳熱系數(shù)。

全局傳熱系數(shù)主要由液膜的導(dǎo)熱決定，而不是可能出現(xiàn)在波內(nèi)的回流，回流對(duì)傳熱速率的影響很小。對(duì)此Jayanti等[35]給出了驗(yàn)證，他們使用CFD技術(shù)研究了波形界面上液膜的傳熱特性，得到了假想界面形狀上周期層流的溫度場(chǎng)，研究了正弦波和孤立波對(duì)傳熱的影響；同時(shí)得出界面波的存在也強(qiáng)化了傳熱系數(shù)，主要是有效地減薄了液膜厚度。

由于大部分流體是通過(guò)大波傳遞，基底液膜會(huì)更薄，有效傳熱系數(shù)會(huì)更大。薄液膜內(nèi)波對(duì)傳熱的影響是間接的，主要是由于使液膜有效地變薄，而不是強(qiáng)化液膜內(nèi)的對(duì)流。

劉晨等[36]分別從流涕黏度、表面張力、流動(dòng)阻力3個(gè)方面對(duì)流體流動(dòng)過(guò)程中的傳質(zhì)問(wèn)題進(jìn)行研究，其中流體黏度、表面張力會(huì)影響液膜厚度。液膜厚度越大，越不利于傳質(zhì)；流動(dòng)阻力也會(huì)阻礙傳質(zhì)；另外，表面張力會(huì)影響液膜流動(dòng)的波動(dòng)性，從而影響傳質(zhì)。

谷芳等[37]對(duì)豎直平板氣液兩相逆流的降膜傳質(zhì)過(guò)程進(jìn)行了數(shù)值模擬和實(shí)驗(yàn)研究，分析了在不同液相進(jìn)口濃度和不同氣相進(jìn)口速度條件下對(duì)液相出口濃度的影響。此外，劉國(guó)兵等[38]同樣對(duì)豎直平板上水和空氣的逆流模型進(jìn)行了傳質(zhì)模擬，但與谷芳不同的是其考慮了氣-液兩相間能量源項(xiàng)的條件，分析了溫度對(duì)傳質(zhì)的影響。蔣翔等[39]將立式蒸發(fā)式冷凝器中的兩相流動(dòng)簡(jiǎn)化為豎直圓管中氣液順流的降膜傳質(zhì)流動(dòng)。和氣液逆流模型相同，同樣對(duì)氣相的流速、液相的流量對(duì)傳質(zhì)的影響進(jìn)行模擬分析，另外也對(duì)進(jìn)口空氣的相對(duì)濕度對(duì)傳質(zhì)的影響進(jìn)行了分析，模擬結(jié)果與實(shí)驗(yàn)值能很好吻合。

Xu等[40]對(duì)垂直下降液膜進(jìn)行了數(shù)值模擬，研究了波的動(dòng)力學(xué)特性及相關(guān)的傳質(zhì)現(xiàn)象，并與先前實(shí)驗(yàn)得到的經(jīng)驗(yàn)關(guān)系進(jìn)行了比較。劉玉峰等[26]同樣模擬了豎壁降膜流場(chǎng)中流動(dòng)液膜的波動(dòng)特性對(duì)傳質(zhì)的影響。

模擬結(jié)果顯示，在液膜流動(dòng)過(guò)程中，液膜表面穩(wěn)定的波動(dòng)有利于氣液間的傳質(zhì)，且增強(qiáng)液膜表面的波動(dòng)程度可增強(qiáng)傳質(zhì)過(guò)程。增加氣液入口的溫差和降低進(jìn)口空氣的相對(duì)濕度也可增強(qiáng)氣-液相間熱質(zhì)交換的劇烈程度。

5 結(jié) 語(yǔ)

薄膜流動(dòng)廣泛地存在于自然界和許多重要工業(yè)領(lǐng)域中，它是一種簡(jiǎn)單的低速開(kāi)式流動(dòng)，但卻有著非常復(fù)雜豐富的動(dòng)力學(xué)行為，這就使得對(duì)薄膜流動(dòng)的研究有著重大的現(xiàn)實(shí)和理論意義。在對(duì)液膜流動(dòng)特性研究時(shí)，發(fā)現(xiàn)板面結(jié)構(gòu)、液相物性、氣液相流速等對(duì)流動(dòng)特性的影響較大。研究表面波即液膜不穩(wěn)定時(shí)，一般在入口處會(huì)有一定頻率速度擾動(dòng)，擾動(dòng)會(huì)發(fā)展成一個(gè)包含大振幅“滾波”和小振幅毛細(xì)波的孤立波。由于液膜表面存在波動(dòng)，會(huì)增大氣液接觸，從而會(huì)強(qiáng)化傳熱傳質(zhì)。

目前，對(duì)液膜流動(dòng)特性的研究已經(jīng)比較充分，但對(duì)于表面波通常都是進(jìn)行二維數(shù)值模擬，而液膜表面波波動(dòng)具有三維特征，因此模擬結(jié)果與實(shí)際存在一定差距。另外，影響液膜傳熱傳質(zhì)特性的因素種類(lèi)繁多，在某些方面，所得的認(rèn)識(shí)規(guī)律尚不統(tǒng)一，因此仍需要進(jìn)一步研究。此外，數(shù)值模擬技術(shù)作為液膜研究中的一個(gè)重要手段，屬于間接模擬研究范疇，故對(duì)其模擬結(jié)果的分析需要與理論研究有機(jī)結(jié)合，進(jìn)行互補(bǔ)分析，同時(shí)，應(yīng)進(jìn)行必要的實(shí)驗(yàn)驗(yàn)證，進(jìn)而使流動(dòng)液膜的數(shù)值模擬研究更加完善，準(zhǔn)確性更高。

[1] Albert C，Marschall H，Bothe D. Direct numerical simulation of interfacial mass transfer into falling films[J].，2014，69：343-357.

[2] Lan H，Wegener J L，Armaly B F，et al. Developing laminar gravity-driven thin liquid film flow down an inclined plane[J].，2010，132（8）：081301.

[3] 馬學(xué)虎，薄守石，蘭忠，等. 降液膜波動(dòng)的影響因素分析[J]. 高校化學(xué)工程學(xué)報(bào)，2010，24（1）：10-15.

[4] 薄守石，馬學(xué)虎，蘭忠，等. 降液膜流動(dòng)孤波中的渦流影響因素分析[J]. 工程熱物理學(xué)報(bào)，2011，32（10）：1703-1706.

[5] 薄守石. 溴化鋰溶液降膜流動(dòng)及傳熱傳質(zhì)的數(shù)值研究[D]. 大連：大連理工大學(xué)，2011.

[6] Min J K，Park I S. Numerical study for laminar wavy motions of liquid film flow on vertical wall[J].，2011，54（15）：3256-3266.

[7] Ho C D，Chang H，Chen H J，et al. CFD simulation of the two-phase flow for a falling film microreactor[J].，2011，54（15）：3740-3748.

[8] 宋建，胡珀，韋勝杰，等. 豎直平板水膜表面波演化行為實(shí)驗(yàn)研究[J]. 原子能科學(xué)技術(shù)，2012，46（7）：804-808.

[9] 谷芳. 規(guī)整填料局部流動(dòng)與傳質(zhì)的計(jì)算流體力學(xué)研究[D]. 天津：天津大學(xué)，2004.

[10] 谷芳，劉春江，袁希鋼，等. 傾斜波紋板上液膜流動(dòng)的CFD研究[J]. 化工學(xué)報(bào)，2005，56（3）：462-467.

[11] Gu F，Liu C J，Yuan X G，et al. CFD simulation of liquid film flow on inclined plates[J].，2004，27（10）：1099-1104.

[12] 李超. 規(guī)整填料塔填料內(nèi)流動(dòng)與傳質(zhì)的CFD模擬[D]. 鄭州：鄭州大學(xué)，2011.

[13] 吳治將，朱冬生，張景衛(wèi)，等. 板式蒸發(fā)式冷凝器氣-液兩相流流動(dòng)的數(shù)值模擬[J]. 工程熱物理學(xué)報(bào)，2008，29（6）：971-974.

[14] 孫荷靜. 波紋板蒸發(fā)式冷凝器流體流動(dòng)特性及傳熱實(shí)驗(yàn)研究[D]. 廣州：華南理工大學(xué)，2010.

[15] 劉鑫. 非平整基底面上薄膜涂覆流體運(yùn)動(dòng)的數(shù)值模擬研究[D]. 北京：北京工業(yè)大學(xué)，2010.

[16] 孫鳳丹. 波紋板表面流體流動(dòng)特性的數(shù)值模擬[J]. 現(xiàn)代化工，2012，32（1）：90-93.

[17] 劉梅，宋朝匣，吳正人，等. 非平整傾斜表面上流動(dòng)液膜的影響因素分析[J]. 化工進(jìn)展，2015，34（1）：75-79+112.

[18] 李相鵬，陳冰冰，高增梁. 規(guī)整填料表面液膜流動(dòng)特性的數(shù)值模擬[J]. 化工學(xué)報(bào)，2013，64（6）：1925-1933.

[19] Xu Y，Yuan J，Repke J U，et al. CFD study on liquid flow behavior on inclined flat plate focusing on effect of flow rate[J].，2012，6（2）：186-194.

[20] 樸明日，胡國(guó)輝. 壁面結(jié)構(gòu)對(duì)非定常薄膜流動(dòng)表面波特性的影響[J]. 計(jì)算物理，2011，28（6）：843-852.

[21] Kunugi T，Kino C. Numerical simulation of oscillatory falling liquid film flows[J].，2003，65（3）：387-392.

[22] Kunugi T，Kino C. DNS of falling film structure and heat transferMARS method[J].，2005，83（6）：455-462.

[23] Miyara A. Numerical simulation of wavy liquid film flowing down on a vertical wall and an inclined wall[J].，2000，39（9）：1015-1027.

[24] Valluri P，Matar O K，Hewitt G F，et al. Thin film flow over structured packings at moderate Reynolds numbers[J].，2005，60（7）：1965-1975.

[25] 許松林，趙嬋. 氣液并流垂直液膜流動(dòng)的數(shù)值模擬[J]. 天津大學(xué)學(xué)報(bào)，2014，47（12）：1039-1046.

[26] 劉玉峰，童一峻，任海剛，等. 高雷諾數(shù)豎壁降膜流動(dòng)特性的數(shù)值研究[J]. 紅外技術(shù)，2010，32（10）：567-571.

[27] Szulczewska B，Zbicinski I，Gorak A. Liquid flow on structured packing：CFD simulation and experimental study[J].，2003，26（5）：580-584.

[28] 于意奇，楊燕華，程旭，等. 降膜流動(dòng)行為的數(shù)值模擬研究[J]. 原子能科學(xué)技術(shù)，2012，46（10）：1207-1211.

[29] 孫文倩. 水平管外降膜流動(dòng)過(guò)程實(shí)驗(yàn)研究與數(shù)值模擬[D]. 北京：中國(guó)科學(xué)院研究生院工程熱物理研究所，2013.

[30] 郭常青，朱冬生，蔣翔，等. 板式蒸發(fā)式冷凝器傳熱傳質(zhì)的數(shù)值模擬[J]. 華南理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，37（3）：53-57.

[31] 張景衛(wèi)，朱冬生，吳治將，等. 板式蒸發(fā)式冷凝器 CFD 模擬與實(shí)驗(yàn)研究[C]//中國(guó)制冷學(xué)會(huì) 2007 學(xué)術(shù)年會(huì)論文集，2007.

[32] 朱冬生，張景衛(wèi)，吳治將，等. 板式蒸發(fā)式冷凝器兩相降膜流動(dòng) CFD 模擬及傳熱研究[J]. 華南理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2008，36（7）：6-10.

[33] 孫荷靜，朱冬生，吳治將，等. 波紋管管內(nèi)降膜流動(dòng)與傳熱特性的研究[J]. 石油煉制與化工，2009（10）：25-29.

[34] 周文生，韓東. 降膜蒸發(fā)器傳熱特性的數(shù)值模擬[J]. 化工進(jìn)展，2010，29（s2）：54-60.

[35] Jayanti S，Hewitt G F. Hydrodynamics and heat transfer of wavy thin film flow[J].，1996，40（1）：179-190.

[36] 劉晨，肖晶晶，王紅星，等. 環(huán)己烯水合過(guò)程汽-液兩相流動(dòng)的CFD研究[J]. 化工進(jìn)展，2013，32（11）：2555-2561.

[37] 谷芳，劉春江，余黎明，等. 氣-液兩相降膜流動(dòng)及傳質(zhì)過(guò)程的 CFD 研究[J]. 高?；瘜W(xué)工程學(xué)報(bào)，2005，19（4）：438-444.

[38] 劉國(guó)兵，惠宇，王玉璋. 沿垂直壁面氣-液降膜流動(dòng)傳質(zhì)過(guò)程的數(shù)值研究[J]. 燃?xì)廨啓C(jī)技術(shù)，2012，25（3）：34-39.

[39] 蔣翔，朱冬生，吳治將，等. 立式蒸發(fā)式冷凝器傳熱傳質(zhì)的 CFD 模擬[J]. 高?；瘜W(xué)工程學(xué)報(bào)，2009，23（4）：566-571

[40] Xu Z F，Khoo B C，Wijeysundera N E. Mass transfer across the falling film：Simulations and experiments[J].，2008，63（9）：2559-2575.

Review on the numerical simulation of the flowing liquid film

WU Zhengren1，SONG Zhaoxia1，LIU Mei2，LIU Qiusheng1

（1School of Energy Power and Mechanical Engineering，North China Electric Power University，Baoding 071003，Hebei，China；2Department of Economic Management，North China Electric Power University，Baoding 071003，Hebei，China）

As an efficient heat and mass transfer technology，liquid film flow has been extensively used in chemical engineering and many other industries. In recent years，numerical simulation methods were widely applied in the study of the properties of the liquid film flow and heat transfer characteristics. This paper analyzed the liquid film free surface tracking methods of numerical simulation. Four aspects，including different wall structure，different wall angle，liquid property，liquid flow rate and gas velocity were summarized. The features of fluctuations were investigated through changing the wall angle，the inlet Reynolds number and inlet disturbance. In addition，the research status of heat and mass transfer characteristics were summarized. The disadvantages of using numerical simulation method to study liquid film flow were also discussed. The more scientific and reasonable methods of researching the liquid film flow were prospected.

liquid film；flow dynamic；surface waves；heat and mass transfer；numerical simulation

TQ 028

A

1000–6613（2015）09–3216–06

10.16085/j.issn.1000-6613.2015.09.002

2015-01-26；修改稿日期：2015-03-12。

高等學(xué)校博士學(xué)科點(diǎn)專(zhuān)項(xiàng)科研基金（20110036110009）及國(guó)家自然科學(xué)基金（11302076）項(xiàng)目。

吳正人（1973—），男，講師，博士，研究方向?yàn)榱黧w動(dòng)力學(xué)數(shù)值模擬研究。E-mail zhengren_wu@163.com。