流動液膜的數值模擬研究綜述

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

?

流動液膜的數值模擬研究綜述

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

（1華北電力大學能源動力與機械工程學院，河北保定 071003；2華北電力大學經濟管理系，河北保定 071003）

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

液膜；流體動力學；兩相流；傳熱傳質；數值模擬

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

20世紀中期，Kapitza父子進行了液體薄膜沿傾斜壁面下落的開創性實驗。至此以后的半個多世紀，對薄膜流動的動力學性能的研究引起了許多科研工作者的關注。此外，早期對于液膜流動的研究主要集中在平整基底上，隨著研究的進展和工程生產中的需要，對液膜流動的研究就不僅僅局限于平整壁面。近年來，對非平整壁面或管內、外液膜的流動特性及傳熱傳質特性的研究成為熱點。

對液膜流動的研究主要包括理論分析、實驗研究和數值模擬等方法，本文主要對基于數值模擬技術研究液膜流動特性的現狀進行綜述分析。

1 液膜自由界面追蹤方法

對液膜流動進行數值模擬時，由于液膜自由液面隨時間變化很難準確追蹤捕捉到，從而使得數值模擬存在很大困難。有關學者經過研究提出了幾種自由邊界的處理方法，包括液膜厚度方程法（height function）、線段法（line segments）、標記粒子法（markerparticles）和VOF方法（volume of fluid）。

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

2 液膜流動特性

影響液膜流動特性的因素很多，本文主要從不同壁面結構、不同壁面傾角、液體物性、液相流量與氣相流速4個方面進行分析總結。

2.1 不同壁面結構的影響

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

由以上研究，通過對比平板和波紋板上的流動可知，液膜在長度相同的波紋板和平板上流動時，平板的液膜流動時間比波紋板上的短，這是由于波紋結構使液膜的流動路程增大，從而使液膜在板面上的停留時間變長，這樣為傳熱傳質提供了更有利的條件。

通過對不同板面結構對液膜流動的影響研究發現，波紋特征尺寸越大，液體越容易在波紋谷底形成漩渦，從而造成液體滯留，不利于板面上流體的更新。并且液膜在波紋尺寸較小的板面上容易發生斷裂，使板面上形成“干斑”，而在波紋尺寸較大的板面上容易形成連續液膜。另外，研究者還測定了不同壁面結構下剪切應力的分布，并分析得出近壁區漩渦的位置和大小可通過測定壁面剪切應力來預測。

2.2 不同壁面傾角的影響

液膜厚度、壁面剪切應力等都會隨著板面傾角變化而有一定的變化。研究結果表明，隨波紋板傾角增加，液膜厚度逐漸增加，而壁面剪切應力呈逐漸減小趨勢，但板面傾角的變化不會改變液膜的基本流動結構，即液膜自由液面與板面的相位差。同時研究發現，隨著傾角的逐漸增大，波紋板進出口壓降先增大后減小再增大。

2.3 液體物性的影響

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

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

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

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

3 表面波特性

由于液膜表面波高度影響著液膜內的傳熱傳質，因此對其的關注具有深刻意義。

3.1 入口擾動的影響

液膜表面波受入口擾動頻率及壁面結構的共同影響。當入口擾動頻率很小時，低頻率擾動發展形成孤立波；此外，表面波的振幅和波長會隨著入口擾動頻率的增大而減小。當入口擾動頻率較高時，表面波會受到很強的干擾，甚至發生合并現象。當表面波的波長和基底結構波長一致時，就會發生共振現象。共振是沿波紋壁面下落的液膜流動的獨特的物理現象。

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

研究發現，在入口邊界處的小振幅擾動會發展成包含一個大振幅“滾波”和小振幅毛細波的孤立波。沿著流動方向在滾波內產生渦流，渦流的波幅會發展成與滾波波幅大小相近，此外，渦流對于強化傳熱傳質有重要作用。另一方面，在毛細波內沒有渦流產生。高頻的擾動相互作用，滾波的振幅減小，毛細波伴隨著滾波消失，同時在已經發展成的波內會有渦流產生，在一個特定頻率下，它會由于波的相互作用重復產生和消失。

3.2 雷諾數（Re）的影響

目前，有關液膜流動的研究中涉及了有關準數的影響作用，如等。其中對于波動液膜的發展有重要影響。本文對數的影響作用進行了總結研究。

研究表明，數增大時，平均液膜厚度變大，波幅變小，孤立波逐漸消失，小振幅的波逐漸充滿整個流場；另外，數增大，沿著下游方向波幅增大，上游的波合并在一起并擴大，沒有發生衰減，最后在下游處產生了高振幅的孤立波。

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

3.3 壁面傾角的影響

當壁面傾斜角發生變化時，表面波也會發生變化，隨著傾斜角降低，波動發展緩慢，波峰處的循環流動及波谷處的渦流均逐漸消失。馬學虎[3-4]、Miyara[23]等對此進行了研究分析。

4 傳熱傳質特性

對于氣液界面傳熱傳質的研究，最早主要集中于理論和實驗方面。近年來，采用數值模擬方法研究液膜的傳熱傳質特性成為熱點。

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

界面波的強化傳熱機理可以通過顯示在不同下的溫度場分布來說明。Miyara[23]在流動入口處添加小振幅正弦擾動，數值模擬了豎直平板上帶有界面波的下降液膜的流動特性和傳熱特性。研究結果表明：界面波通過兩種效應來強化傳熱，即削薄液膜和對流效應。對于較低的數，可以忽略對流效應；對于較高的數，這兩種效應都起作用。另外，時均傳熱系數沿流動方向逐漸增大。在完全發展區域，較高數下的傳熱系數要更高一點。然而，在發展區域，最大的數有最低的傳熱系數。

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

由于大部分流體是通過大波傳遞，基底液膜會更薄，有效傳熱系數會更大。薄液膜內波對傳熱的影響是間接的，主要是由于使液膜有效地變薄，而不是強化液膜內的對流。

劉晨等[36]分別從流涕黏度、表面張力、流動阻力3個方面對流體流動過程中的傳質問題進行研究，其中流體黏度、表面張力會影響液膜厚度。液膜厚度越大，越不利于傳質；流動阻力也會阻礙傳質；另外，表面張力會影響液膜流動的波動性，從而影響傳質。

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

Xu等[40]對垂直下降液膜進行了數值模擬，研究了波的動力學特性及相關的傳質現象，并與先前實驗得到的經驗關系進行了比較。劉玉峰等[26]同樣模擬了豎壁降膜流場中流動液膜的波動特性對傳質的影響。

模擬結果顯示，在液膜流動過程中，液膜表面穩定的波動有利于氣液間的傳質，且增強液膜表面的波動程度可增強傳質過程。增加氣液入口的溫差和降低進口空氣的相對濕度也可增強氣-液相間熱質交換的劇烈程度。

5 結 語

薄膜流動廣泛地存在于自然界和許多重要工業領域中，它是一種簡單的低速開式流動，但卻有著非常復雜豐富的動力學行為，這就使得對薄膜流動的研究有著重大的現實和理論意義。在對液膜流動特性研究時，發現板面結構、液相物性、氣液相流速等對流動特性的影響較大。研究表面波即液膜不穩定時，一般在入口處會有一定頻率速度擾動，擾動會發展成一個包含大振幅“滾波”和小振幅毛細波的孤立波。由于液膜表面存在波動，會增大氣液接觸，從而會強化傳熱傳質。

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

[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] 馬學虎，薄守石，蘭忠，等. 降液膜波動的影響因素分析[J]. 高校化學工程學報，2010，24（1）：10-15.

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

[5] 薄守石. 溴化鋰溶液降膜流動及傳熱傳質的數值研究[D]. 大連：大連理工大學，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] 宋建，胡珀，韋勝杰，等. 豎直平板水膜表面波演化行為實驗研究[J]. 原子能科學技術，2012，46（7）：804-808.

[9] 谷芳. 規整填料局部流動與傳質的計算流體力學研究[D]. 天津：天津大學，2004.

[10] 谷芳，劉春江，袁希鋼，等. 傾斜波紋板上液膜流動的CFD研究[J]. 化工學報，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] 李超. 規整填料塔填料內流動與傳質的CFD模擬[D]. 鄭州：鄭州大學，2011.

[13] 吳治將，朱冬生，張景衛，等. 板式蒸發式冷凝器氣-液兩相流流動的數值模擬[J]. 工程熱物理學報，2008，29（6）：971-974.

[14] 孫荷靜. 波紋板蒸發式冷凝器流體流動特性及傳熱實驗研究[D]. 廣州：華南理工大學，2010.

[15] 劉鑫. 非平整基底面上薄膜涂覆流體運動的數值模擬研究[D]. 北京：北京工業大學，2010.

[16] 孫鳳丹. 波紋板表面流體流動特性的數值模擬[J]. 現代化工，2012，32（1）：90-93.

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

[18] 李相鵬，陳冰冰，高增梁. 規整填料表面液膜流動特性的數值模擬[J]. 化工學報，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] 樸明日，胡國輝. 壁面結構對非定常薄膜流動表面波特性的影響[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] 許松林，趙嬋. 氣液并流垂直液膜流動的數值模擬[J]. 天津大學學報，2014，47（12）：1039-1046.

[26] 劉玉峰，童一峻，任海剛，等. 高雷諾數豎壁降膜流動特性的數值研究[J]. 紅外技術，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] 于意奇，楊燕華，程旭，等. 降膜流動行為的數值模擬研究[J]. 原子能科學技術，2012，46（10）：1207-1211.

[29] 孫文倩. 水平管外降膜流動過程實驗研究與數值模擬[D]. 北京：中國科學院研究生院工程熱物理研究所，2013.

[30] 郭常青，朱冬生，蔣翔，等. 板式蒸發式冷凝器傳熱傳質的數值模擬[J]. 華南理工大學學報：自然科學版，2009，37（3）：53-57.

[31] 張景衛，朱冬生，吳治將，等. 板式蒸發式冷凝器 CFD 模擬與實驗研究[C]//中國制冷學會 2007 學術年會論文集，2007.

[32] 朱冬生，張景衛，吳治將，等. 板式蒸發式冷凝器兩相降膜流動 CFD 模擬及傳熱研究[J]. 華南理工大學學報：自然科學版，2008，36（7）：6-10.

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

[34] 周文生，韓東. 降膜蒸發器傳熱特性的數值模擬[J]. 化工進展，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] 劉晨，肖晶晶，王紅星，等. 環己烯水合過程汽-液兩相流動的CFD研究[J]. 化工進展，2013，32（11）：2555-2561.

[37] 谷芳，劉春江，余黎明，等. 氣-液兩相降膜流動及傳質過程的 CFD 研究[J]. 高校化學工程學報，2005，19（4）：438-444.

[38] 劉國兵，惠宇，王玉璋. 沿垂直壁面氣-液降膜流動傳質過程的數值研究[J]. 燃氣輪機技術，2012，25（3）：34-39.

[39] 蔣翔，朱冬生，吳治將，等. 立式蒸發式冷凝器傳熱傳質的 CFD 模擬[J]. 高校化學工程學報，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。

高等學校博士學科點專項科研基金（20110036110009）及國家自然科學基金（11302076）項目。

吳正人（1973—），男，講師，博士，研究方向為流體動力學數值模擬研究。E-mail zhengren_wu@163.com。