豎直降液膜流動在反應工程中的應用

陳世昌，馬建平，張先明，陳文興

（浙江理工大學紡織纖維材料與加工技術國家地方聯合工程實驗室，浙江 杭州 310018）

進展與述評

豎直降液膜流動在反應工程中的應用

陳世昌，馬建平，張先明，陳文興

（浙江理工大學紡織纖維材料與加工技術國家地方聯合工程實驗室，浙江 杭州 310018）

豎直元件上的降液膜流動由于具有傳熱傳質系數高、氣液接觸面積大和表面更新速度快等優點，因而在化工和聚合反應工程上有著廣泛應用。本文綜述了豎直降膜元件上液膜流動的實驗研究、理論分析和數值模擬，分別介紹了低黏流體降液膜流動在氣體吸收與化學反應等方面的應用和高黏流體降液膜流動在涂覆與聚合物脫揮等聚合物加工中的應用情況。總結表明，目前對于低黏流體在豎直元件上的降液膜流動研究已較為成熟，而高黏流體降液膜流動的大部分研究僅著重于其成膜性能上。最后分析了降液膜流動的研究趨勢，指出涉及熱質傳遞與化學反應耦合的高黏流體在豎直結構上的降液膜流動是今后的重點研究對象。

降液膜；流動不穩定性；傳熱；傳質；聚合；反應工程

豎直降膜元件上的降液膜流動是一種特殊的受迫運動，流體流動時液膜內表面黏附于固壁上，另一表面暴露于真空或特殊氣體氛圍中而被稱做自由面，因具有溫差小、傳熱傳質系數高、氣液接觸面積大、表面更新快、結構簡單以及動力消耗小等優點，作為一種高效熱質傳遞單元操作技術在工程熱物理、石油化工、核能技術等領域有著廣泛應用。降液膜流動過程具有流體動力學行為豐富多變的特點，如其自由面更新頻率快、多種波以及流場的演化等問題一直受到研究者廣泛關注。近年來，許多化工研究者采用直管等立式結構裝置，進行了氣體吸收、化學反應和降溫換熱等過程的研究，而高分子溶液涂膜、聚合物脫揮和聚酯增黏等高黏流體在豎直降膜元件上流動的同時來完成相應化工單元操作的技術也引起了越來越多研究者的興趣，且其中一些技術已在實際生產中有所應用。利用豎直降液膜的三維波動特性來強化學工程中的熱質傳遞問題，已成為傳遞科學與工程界一個十分活躍的研究領域。

1 豎直降液膜流動的不穩定性

1.1 不穩定性機理

降液膜流動的研究始于Kapitza父子[1]在豎直玻璃管上進行的實驗，他們在實驗中還觀察到了不穩定的表面波和孤波現象；Benjamin[2]對任意傾斜角（β）平板上的降液膜流動進行研究，發現豎直（β=π/2）情況下觀察不到均勻穩定的層流，液膜流動演化為復雜的運動；Yih[3]利用攝動法對沿豎直平板下落的恒溫薄膜進行了線性穩定性分析，獲得穩定流動條件下的臨界雷諾數（Rec=5cotβ/4），這說明豎直降液膜過程中，不管雷諾數多小，其流動總是長波不穩定。

因為自由表面波不穩定性是降液膜流動的主要特征，對于相應的不穩定機理及其影響因素有必要展開深入的研究。Smith[4]分析了長波不穩定性的相關機理，指出薄膜流動不穩定機理由液膜流動不穩定性的引發機理和產生界面不穩定運動的增長機理兩部分組成。對于交界面為切向應力邊界的薄膜流動來說，當界面變形時，會引發產生擾動的剪切應力，其大小與界面處速度的曲率成比例。應力產生的作用是不穩定性所需能量的主要來源，而應力擾動帶來的主要影響在于其促進擾動在薄膜縱向流動中放大，進而導致波動加劇。杜小澤等[5]根據能量變化分析流動慣性力與表面張力對薄液膜的波動影響，認為液膜蒸發時熱流率的增加導致了維持液膜穩定作用的壁面毛細吸附力減小，液膜發生斷裂。王補宣等[6]則考慮降膜流動特征與壁面熱流率密切相關，認為液膜流動狀態的變化源于薄液膜內熱狀態不平衡。馬學虎等[7]研究固壁與流體的界面作用，得出不同物理化學性質的固壁可以通過改變固液表面自由能差而影響降液膜平均厚度等流動特性。盧川等[8]則認為降液膜的不穩定流動具有區域分布的特點，其不穩定的主要原因在于降膜形成的小波、孤波與正弦波等前后波發生融合。

1.2 豎直結構上液膜流動的時空演化

Takahama和Kato[9]對豎直圓管外壁上的降液膜流動特性進行了實驗研究，隨著液膜縱向發展以及Re的增大，液膜流動的不穩定性增加，但在一定流動條件下存在著一個Rec，在此臨界值之上波動不但不會加劇，反而會降低。Karapantsios和Karabelas[10]考察了高Re下豎直管內的降液膜流動，也發現了類似的流動特征，但對膜厚的統計分析表明膜厚的變化具有隨機性，而對波峰分析發現波數幾乎不變。Liu等[11]研究了薄液膜在Re小于200時沿平板流動的表面波演化過程，認為小孤波間的相互作用及其發展是薄膜呈波形流動的根本原因，且周期性的擾動會使孤波間隔變得無序。Park和Nosoko[12]觀察了豎直平板上水膜在Re為33時從一維流動演化至二維非線性波并最終發展為馬蹄形三維波的過程（圖1）。Nosoko和Miyara[13]在液膜流動起始處施加一定頻率的擾動，發現降液膜波動方式與受迫擾動頻率相關，當擾動頻率足夠大時，初始二維波極易變為三維波，且流動方式變得極不規則，最終類似于紊流。

隨著實驗技術的發展，一些先進的流場測試技術逐漸被應用到降液膜流動的研究中來。Drosos等[14]利用高速攝像獲得了不同Re下液膜表面的波動情況，研究了不同流速下的降液膜流體動力學特性。Yu等[15]采用激光多普勒測速儀（LDA）得到不同氣液逆流動條件下降膜流動的液相速度分布，結果表明氣液逆流時的降膜流動是以近界面區域的速度波動為特征的流動，這與Nusselt模型明顯不同，其最大速度出現在液膜表面，且氣液界面區域的速度波動可以加快表面更新速率，產生的微小對流可以增強氣液界面的傳質過程。

圖1 水在垂直平板上的降液膜流動演化

1.3 加熱壁面上的降液膜流動

由于液膜傳熱在蒸發、吸收等場合中有著重要的應用，許多研究者對加熱情況下的降液膜流動進行了大量的研究。Bankoff[16]的研究結果表明，若液膜足夠薄，在對熱的壁面進行冷卻時，液膜內部出現很大的溫度梯度，其換熱及傳質強度是同等流動條件下單相流動的幾倍甚至幾十倍。Dao和Balakotaiah[17]對甘油水溶液在豎直管內的降膜吸收進行實驗研究，認為液膜的吸收過程使得降液膜自由表面波的環狀流轉變為平推流，并由此可以確定液膜的臨界流速。Samanta[18]研究了加熱豎直管外的自由降液膜流動，獲得液膜流動形態并測量了液膜厚度，給出了平均液膜厚度與Re的準則關系式。Zhang等[19-20]對均勻加熱或冷卻豎直平板上薄液膜流動的溫度分布和流動模型的研究結果表明，薄膜的流動情況與加熱條件息息相關，冷卻使其液膜表面積增加，而加熱后則會變小，認為是Marangoni效應對流動特性產生重要的影響，最終導致了液膜破裂。

1.4 降液膜流動的理論分析與數值模擬

對于復雜的降液膜流動，通過數值模擬方法建立的各種模型為液膜流動及其傳熱傳質的研究提供了極大的便利。錢煥群等[21]采用邊界層理論對自由降液膜流動的數學模型進行分析簡化，得到液膜流動的二階邊界層模型，根據線性穩定分析，將計算結果分別與Orr-Sommerfeld方程和一階邊界層模型的計算結果進行對比，發現二階邊界層模型能準確描述下降液膜流動的波的特征，且對于高Re和高頻率波動的液膜流動，二階邊界層模型近似效果更加明顯。Pavlidis等[22]研究了在豎直圓筒內表面和外表面以及平坦壁面上完全發展的一維黏彈性薄膜流動，分析并計算了速度分布和壓力對薄膜流動的幾何形狀與本構模型中的參數等產生的影響。

近年來越來越多的交界面數值模擬技術被運用到降液膜的研究中，如移動網格法、邊界積分法、粒子法、流體體積法（VOF）等。其中VOF法由于所需儲存信息少，能自動處理交叉的自由邊界，并且很容易拓展到三維計算中，用它來追蹤不可壓縮流體動力學中的自由表面特別方便，因而被廣泛應用于研究自由降液膜流動問題。Gao等[23]以VOF法模擬薄膜的自由表面，數值計算了N-S方程，通過分析入口處受迫擾動的二維豎直降液膜表面波的不穩定性，認為影響波特點的參數主要是Re和擾動頻率，并比較了高頻低Re下的近正弦形小振幅波和低頻高流速下的孤立波，出口采用自由邊界條件，使用連續表面力方法（CSF）加入表面張力進行計算，將數值結果與已有的實驗相比較，結果比較吻合。馬學虎等[24]利用CSF-VOF方法捕捉降液膜流動的自由表面，考慮表面張力對降液膜流動的作用，研究了液膜隨著入口擾動頻率、Re等因素的變化而波動的演化過程。最近，Bo等[25]在流動邊界入口處引入小振幅受迫擾動，也采用CSF-VOF法研究了表面張力和黏度對波行為與液膜流體動力學的影響。

一些特殊結構如豎直鐵絲和纖維上的降液膜流動也引起了研究者的興趣。Trifonov[26]使用積分法計算沿豎直絲線和管上的黏性降液膜變化的非線性波形結構，根據重力、黏性力和毛細力等作用項得出了有限振幅的行波解；他還對沿波紋面的黏性薄膜流動進行理論分析[27]，考慮大范圍Re和表面幾何特性，計算自由面形狀和其他流動特征，認為對于一維表面，在一定參數范圍內流動主要是受表面張力控制，這種流動只需用積分方法來處理。Grünig等[28]在沿著單根豎直鐵絲逆氣流的薄膜流動實驗中研究了液膜破裂過程，由于細絲上的流體動力學特征不同于平板表面的薄膜流動，表現為特殊的液珠流動，當薄膜遇到明顯的逆向氣流時，流動模式會發生變化，但液體滯留性和界面面積幾乎不受影響，由此液膜發生破裂。Ruyer-Quilet等[29]研究了沿單根豎直纖維上重力驅動的降液膜流動，所用模型中考慮了流動慣性和流向黏性擴散，發現黏性擴散效應對線性波具有穩定分散作用，其模擬結果與實驗相比在線性和非線性體系內都較為一致。在此基礎上，他們還研究了纖維上液膜流動的軸對稱行波的傳播[30]，根據4種不同主導作用的物理效應將薄膜流動分為4種流態：類滴流，慣性驅動流，重力驅動流以及類孤波流，分析了行波從一種流態到另一種流態的復雜轉變。Zakaria等[31]則討論了恒速繞軸轉動的豎直圓筒外表面上黏性降液膜的線性和非線性穩定性條件。

豎直降液膜流動的不穩定特性使其在需要傳熱和傳質的反應工程上能較好的發揮作用，如通過設計不同降膜元件的反應器可以使得溶液在設備內具有不同的停留時間，一般不同的化學反應則需要不同厚度的降液膜。因此，降液膜的流動特性、成膜條件、膜厚以及不同性質流體液膜的傳熱與傳質效率等對反應工程具有重要的影響作用。

2 低黏流體豎直降液膜流動的應用

2.1 氣體吸收

對于工業上氣體吸收（或濕法滌氣）過程中使用的氣液接觸器，其裝備的優化設計旨在讓液體吸收劑以薄膜的形式沿著塔形接觸器的豎直壁面流動，如一些研究者所采用內部件為豎直絲線的氣液接觸器[32-33]，溶液可以沿著豎直絲線做降膜運動，這種結構能明顯擴大了氣液交界面積，促進氣液兩相的混合，從而提高氣液兩相間的傳質與傳熱。Broniarz-Press[34]報道了水與CO2在一個同軸環形空間內氣液逆流動，考察粗糙螺旋的圓筒外表面對氣體在降液膜流動過程中吸收性能的影響。Park和Nosoko[12]研究了水膜沿著豎直壁面（圓筒內表面和外表面，平板外表面）流動時由聲振蕩器產生的表面波特點及其對CO2和O2在水膜中吸收效果的影響。Komae等[35]將柔性硅橡膠帶以多級螺旋的形式緊緊纏繞在豎直圓筒外表面，纏帶將沿著平行通道的水膜流動分為若干個流動區域，相對于在直接圓筒壁面上的降液膜流動，這種結構可以使液體和周圍的氣體相互接觸時間更長，同時離心力可以引發水的振蕩運動，提高吸收氣體性能，但是沿著這種螺旋通道形成的溪流會比在豎直平板壁面上流動的普通水膜要厚，這對氣體吸收是不利的。此外，一些研究者通過使流體在降膜設備上形成湍流、液膜表面波或進行螺旋流動，以此改善壁面液膜的氣體吸收性能[36]。

2.2 化學反應

利用溶液在重力驅動下沿著降膜設備形成薄液膜流動的同時進行化學反應的應用也越來越多。這種降膜所需的固定支撐壁面通常是豎直的薄壁或小直徑列管，該結構特別適用于溶液與氣體的化學反應。由于薄膜有利于強化傳熱與傳質的特性，它適合需要換熱的氣液化學反應，如CO2、SO2和NH3的吸收反應。Puma和Yue[37]研究了一種層流降膜的泥漿光催化反應器，用于水楊酸在氧氣下轉換為CO2和H2O的反應，其采用的反應柱是一根豎直固定表面光滑的不銹鋼管（長度為1.6m，直徑為0.108m），紫外光照射在管外流動的液膜上引發反應，效果良好。德國美茵茲微技術研究所設計的一種微型降膜反應器，其生成的液膜厚度可以小于100μm，這種反應器具有優良的排熱效果，可用于芳烴的直接氟化反應[38]；Yeong等[39]將其用于硝基苯還原為苯胺的反應，在更低濃度原液下日產量提高近一倍。Zaid等[40]將波紋壁面的管式降膜反應器用于直鏈烷基苯的磺化反應，對其波紋壁面的影響進行了研究，與光滑的降膜反應器相比，多級波紋結構明顯促進了反應器的傳熱，提高了反應物的轉化率，但波紋壁面引起的湍流和壁面的液膜停滯展開競爭，在凹進的區域產生不期望的連續反應造成產物色值較差。最近，Sun等[41]將豎直降膜柱用于LiOH和CO2反應生成Li2CO3，LiOH溶液依靠重力在豎直降膜柱內壁形成薄液膜流動，通入CO2與其反應產生沉淀物，溶液沿著豎直壁面流到出口，經過漿料槽再被輸送至降膜柱頂端，漿液在料槽和降膜之間循環流動直至碳化反應完成。隨著溶液溫度升高，Re明顯變大，膜厚和Li2CO3的暴露時間均減小，而較大流率時液膜流動的Re變大，液膜厚度增加，Li2CO3的暴露時間延長，另外，溶液結晶為碳酸根離子從界面處液膜流動到主溶液里提供了便利，使得整體吸收速率提高。一些豎直結構上的降液膜流動在反應工程中的應用如表1所示。

2.3 降溫換熱

在化工換熱器、吸收式制冷機和熱泵等強化傳熱傳質性能的設備中，豎直管式的結構設計有著廣泛的應用。如作為吸收式熱泵主要部件之一的吸收器多采用LiBr溶液在豎直管外降膜吸收來完成較為復雜的熱-質傳遞過程。Medrano等[42]進行了LiBr水溶液在光滑豎直管內壁上的波狀層流降膜過程中水蒸汽吸收的實驗研究，結果表明，在水冷卻降溫條件下降膜吸收的傳質系數比空氣冷卻時提高了3倍以上，但對于H2O-LiBr溶液在豎直管內壁的降膜吸收效果和其他研究者管外的降膜吸收相比幾乎沒有差別。尹銘等[43]考慮降膜吸收過程中工質的物性變化，通過對降膜吸收過程劃分微元進行分析，實驗研究了LiBr溶液在豎直光滑管和4種強化換熱管外降膜吸收的傳熱傳質情況，認為豎直管上低肋、翅片、斜槽等結構的存在一方面加劇了Marangoni效應和物流的擾動，另一方面可以增加換熱面積，從而使得換熱得到強化，且在一定范圍里，管內表面構型與管外相比對LiBr溶液的降膜傳熱影響更顯著。Bo等[44]利用CFD模擬豎直管外LiBr溶液的降膜吸收過程，由于界面處與吸收相關的不同加熱效果以及冷卻溶液在較小流動距離上壁面處的冷卻作用，溫度呈現指數分布，其梯度較大，隨著流動距離增加，溫度接近于線性分布，而吸收的傳熱與傳質流率由入口處的最大值逐漸變小。

表1 不同豎直結構上降膜流動在反應工程中的應用

3 高黏流體的液膜流動

降液膜流動在蒸發、冷凝、吸收、換熱、水動力等工業系統中的應用，因其主要工質是水、醇類或溴化鋰的水溶液等黏度極低的流體，對應的降液膜流動大多數發生在層流的中高Re區或湍流區。但對于一些高分子溶液的流動，如涂膜、乳膠制備以及聚合物脫揮等，流體的黏度常常很高，液膜流動時的Re并不大，尤其是對于一些伴隨化學反應的傳質過程，降膜的流體黏度可達1000 Pa?s以上，這類高黏流體的降液膜流動過程有待進一步深入研究。

3.1 涂覆中的液膜流動

旋轉涂膜中薄膜體系的宏觀不穩定會導致液膜流動呈現不均勻的擴展。Cheng等[45]采用數值逼近法研究了黏彈性聚合物流體薄膜在豎直滾筒上涂覆流動的非線性動力學的穩定性。當減小黏彈效應或降低旋轉滾筒的速度時，液膜穩定性增加；在較高Re時增加滾筒半徑，旋轉作用會使薄膜流動變得更加不穩定。Boulogne等[46]報道了高分子溶液涂覆豎直纖維的薄膜動力學實驗研究，在該實驗體系中，主要考慮的是毛細力、彈性和黏性力，忽略不計慣性力和重力。他將高分子溶液在纖維上的薄膜流動分成三段區域：在平坦薄膜區入口處，薄膜厚度均一；再進入具有規則波動的不穩定（Plateau-Rayleigh instability）區，均勻的薄膜發展為由薄膜分離的滴狀流[47]，且此區域內液膜流動形貌依賴于流體的非牛頓特性；最后進入無序流動區，液滴發生聚集擾亂了流動而產生無序的流動模式。他們分別采用剪切變稀明顯的黃原膠/水/甘油溶液和法向應力特性突出的聚丙烯酰胺/水/甘油溶液來考察剪切變稀和法向應力對液膜在纖維上流動形態的影響。實驗結果表明：剪切變稀使得牛頓流體在平膜區域的拋物線型速度分布逐漸變為塞狀分布；恒定流動速率下，由于剪切變稀效應，液膜厚度總是小于非牛頓流體的情況；這些結論對于了解復雜高彈性效應流體的液膜流動特性具有非常重要的作用。

3.2 聚合物脫揮中的降液膜流動

在一些較高溫度的降液膜流動過程中，由于傳質與化學反應同時進行，流體的黏度會隨著液膜的流動而逐漸升高，此時液膜的流動與成膜特性決定了傳質效率的高低，故而降膜設備的結構顯得尤為重要，聚合物脫揮正體現了這一過程。Stober等[48]設計了一種落條式薄膜脫揮器，聚合物熔體通過噴嘴或窄縫擠出，呈液滴、液柱或液膜狀落下，氣液接觸面積由此增大，熔體中的揮發分擴散距離縮短，從而傳熱傳質效率得到提高，加速了脫揮和縮聚過程。日本旭化成公司開發的沿絲降膜聚合裝置[49]，其熔體原料經分布器分布在集束的金屬絲或鏈等導向元件，形成落條式自由降膜流動，可以獲得極高的比表面積，在高真空環境下實現短時間內的高效脫揮。

由于高黏度的流體降液膜流動都是在一些結構特殊的設備中進行，流場結構較為復雜，液膜流動過程受設備幾何結構、流動條件和物性參數等影響較大，對其流體動力學研究較為困難。Vijayraghvan和Gupta[50]將豎直圓盤部分浸沒于真空泵油中，流體在旋轉圓盤上成膜流動，離心力、黏性力、慣性力、表面張力以及重力共同作用決定了液膜的形狀與膜厚，實驗數據表明液膜厚度與系統的轉速以及物性相關，通過積分還可以計算旋轉圓盤上液體的霧沫夾帶。在此基礎上，Cheong和Choi[51-52]研究了聚合物熔體在圓盤反應器上的流動成膜現象，探討膜厚的影響因素以及圓盤表面液膜分布的規律，對Vijayraghva所提出的膜厚表達式進行校正，其形式更能準確反映膜厚的變化規律。王良生等[53]在研究旋轉圓盤反應器的熔體成膜性能的基礎上，借助流體力學基本原理對液膜流動的成膜機理進行了分析。周賢爵等[54]則研究了籠式反應器中介質黏度、液位、轉速等對液膜流動及混合性能的影響。此外，Hasan和Naser[55]采用CFD模型對甘油溶液在轉鼓上的流動及其成膜性進行了二維模擬研究，得到了半穩態的膜厚值。

3.3 降液膜流動在聚合反應工程中應用

由于降液膜流動能促進聚合物的脫揮過程，因而在聚酯、尼龍等聚合反應工程中也有所應用。Xi等[56]研究了一種柵縫式縮聚反應器中的自由降液膜流動，在這種多層豎直的柵板上，聚合物熔體以薄膜形式靠重力驅動而自由降膜，其流動狀態近似于平推流，由于液膜厚度小，可以避免流體靜壓頭對傳質產生阻礙作用，提高了熱質傳遞效果。但是單個柵縫式降膜結構流程短，需要采用多級柵縫部件串聯來提高聚合物熔體的停留時間和表面更新頻率，而多級結構的采用會導致反應器內自清潔效果降低，容易引起滯留并造成降解。陳文興等[57-58]開發了一種應用于生產高黏聚酯的液相增黏技術，在立式反應釜內采用特殊的降膜結構，熔體從設備頂部進入并進行分配，均勻地從塔式成膜裝置的表面利用自身重力緩慢而穩定的滑落向下，形成厚度均勻的連續薄膜，可以滿足較大的成膜面積和較快的表面更新頻率，加快縮聚反應的速度，從而在較短的時間內實現了小分子的脫揮，到達聚酯增黏的目的。利用該項技術生產出的聚酯熔體特性黏度可達1.02dL/g，b值控制在3.0以下，端羧基濃度降至25mol/t，與傳統的臥式反應器相比，產品質量有所提高，生產成本明顯降低。

4 結 語

降液膜流動是一種低速的開式運動，有著極其豐富的流體動力學行為，其流動不穩定和高效熱質傳遞特性等在反應工程中有具廣泛的應用前景。對于在豎直結構上進行化學反應的降液膜流動，液膜厚度與流動特征等對反應的影響已取得了一定的進展，但降液膜的溫度分布、自由面污染、化學反應引發表面不穩定等問題都有待深入研究。由于高黏流體參與反應的降液膜流動大多數發生在較為復雜的降膜設備中，其流場結構復雜多變，目前的研究大多數僅集中于其流動成膜特性上，如膜厚及其分布規律等，而低黏流體降膜流動中的一些實驗現象，如波動性、流動時空演化等動力學特征，以及相關的數值模擬與計算等工作在高黏流體中并未得到體現；此外，關于需要考慮傳質與化學反應耦合的聚合物流體降液膜流動，相關的實驗和理論研究更是鮮見報道。隨著高效傳熱與傳質技術的發展，高黏流體的降液膜流動在反應工程上的應用必將越來越廣闊。

[1]Kapitza P L，Kapitza S P. Wave flow of thin viscous fluid layers of liquid[J].Zh Eksp Teor Fiz，1949，19：105-120.

[2]Benjamin T B. Wave formulation in laminar flow down an inclined plane[J].Journal of Fluid Mechanics，1957，2（6）：554-573.

[3]Yih C S. Stability of liquid flow down an inclined plane[J].Physics of Fluids，1963，6：321-334.

[4]Smith M K. The mechanism for the long-wave instability in thin liquid films[J].Journal of Fluid Mechanics，1990，217：469-485.

[5]杜小澤，昊少融，姜勝耀，等. 豎管內蒸發降液膜流動的不穩定性分析[J]. 原子能科學技術，2001，35（5）：406-41.

[6]王補宣，張金濤，彭曉峰. 豎壁液膜流壁面熱流率對流動的影響的實驗研究[J]. 工程熱物理學報，2001，22（1）：89-91.

[7]馬學虎，高大志，蘭忠，等. 功能表面材料與流體界面相互作用對垂直降液膜流動特性的影響[J]. 高校化學工程學報，2004，18（3）：269-274.

[8]盧川，段日強，姜勝耀.下降液膜流動不穩定性實驗研究[J]. 清華大學學報，2008，48（9）：107-109.

[9]Takahama H，Kato S. Longitudinal flow characteristics of vertically falling liquid films without concurrent gas flow[J].International Journal of Multiphase Flow，1980，6（3）：203-215.

[10]Karapantsios T D，Karabelas A J. Longitudinal characteristics of wavy falling films[J].International Journal of Multiphase Flow，1995，21（1）：119-127.

[11]Liu J，Gollub J P. Solitary wave dynamics of film flows[J].Physics of Fluids，1994，6（5）：1702-1707.

[12]Park C D，Nosoko T. Three-dimensional wave dynamics on a falling film and associated mass transfer[J].AIChE Journal，2003，49（11）：2715-2727.

[13]Nosoko T，Miyara A. The evolution and subsequent dynamics of waves on a vertically falling liquid film[J].Physics of Fluids，2004，16（4）：1118-1126.

[14]Drosos E I P，Paras S V，Karabelas A J. Characteristics of developing free falling films at intermediate Reynolds and high Kapitza numbers[J].International Journal of Multiphase Flow，2004，30（7-8）：853-876.

[15]Yu L M，Zeng A W，Yu K T. Effect of interfacial velocity fluctuations on the enhancement of the mass-transfer process in falling-film flow[J].Industrial Engineering Chemistry Research，2006，45（3）：1201-1210.

[16]Bankoff S G. Dynamics and stability of thin heated liquid films[J].Journal of Heat Transfer，1990，112（3）：538-546.

[17]Dao E K，Balakotaiah V. Experimental study of wave occlusion onfalling films in a vertical pipe[J].AIChE Journal，2000，46（7）：1300-1306.

[18]Samanta A. Stability of liquid film falling down a vertical non-uniformly heated wall[J].Physica D：Nonlinear Phenomena，2008，237（20）：2587-2598.

[19]Zhang F，Wu Y T，Geng J，et al. An investigation of falling liquid films on a vertical heated/cooled plate[J].International Journal of Multiphase Flow，2008，34（1）：13-28.

[20]Zhang F，Tang D L，Geng J，et al. Study on the temperature distribution of heated falling liquid films[J].Physica D：Nonlinear Phenomena，2008，237（7）：867-872.

[21]錢煥群，胡志華，孫賀東，等.下降液膜流動模型及穩定性分析[J].熱能動力工程，2003，18（1）：82-85.

[22]Pavlidis M，Dimakopoulos Y，Tsamopoulos J. Fully developed flow of a viscoelastic film down a vertical cylindrical or planar wall[J].Rheologica Acta，2009，48（9）：1031-1048

[23]Gao D，Morley N B，Dhir V. Numerical simulation of wavy falling film flow using VOF method[J].Journal of Computational Physics，2003，192（2）：624-642.

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

[25]Bo S S，Ma X H，Lan Z，et al. Numerical simulation on wave behavior and flow dynamics of laminar-wavy falling films：Effect of surface tension and viscosity[J].The Canadian Journal of Chemical Engineering，2012，90（1）：61-68.

[26]Trifonov Y Y. Steady-state traveling waves on the surface of a viscous liquid film falling down on vertical wires and tubes[J].AIChE Journal，1992，38（6）：821-834.

[27]Trifonov Y Y. Viscous film flow down corrugated surfaces[J].Journal of Applied Mechanics and Technical Physics，2004，45（3）：389-400.

[28]Grünig J，Skale T，Kraume M. Liquid flow on a vertical wire in a countercurrent gas flow[J].Chemical Engineering Journal，2010，164（1）：121-131.

[29]Ruyer-Quilet C，Treveleyan P，Giorgiutti-dauphiné F，et al. Modelling film flows down a fiber[J].Journal of Fluid Mechanics，2008，603：431-462.

[30]Ruyer-Quilet C，Kalliadasis S. Wavy regimes of film flow down a fiber[J].Physical Review E，2012，85（4）：6302-6324.

[31]Zakaria K，Gamiel Y. Viscous falling film instability around a vertical moving cylinder[J].Acta Mechanica Sinica，2012，28（2）：253-265.

[32]Migita H，Soga K，Mori Y H. Gas absorption in a wetted-wire column[J].AIChE Journal，2005，51（8）：2190-2198.

[33]Pakdehi S G，Taheri S. Separation of hydrazine from air by wetted wire column[J].Chemical Engineering & Technology，2010，33（10）：1687-1694.

[34]Broniarz-Press L. Enhancement of mass transfer coefficients in spiral films[J].International Journal of Heat and Mass Transfer，1997，40（17）：4197-4208.

[35]Komae N，Sato Y，Mori Y H. Helical liquid flow on a vertical cylinder and its application to gas absorption[J].AIChE Journal，2013，59（8）：3109-3118.

[36]Killion J D，Garimella S. A critical review of models of coupled heat and mass transfer in falling-film absorption[J].International Journal of Refrigeration，2001，24（8）：755-797.

[37]Puma G L，Yue P L. A laminar falling film slurry photocatalytic reactor. Part I：Model development[J].Chemical Engineering Science，1998，53（16）：2993-3006.

[38]J?hnisch K，Baerns M，Hessel V，et al. Direct fluorination of toluene using elemental fluorine in gas/liquid microreactors[J].Journal of Fluorine Chemistry，2000，15（1）：117-128.

[39]Yeong K K，Gavriilidis A，Zapf R，et al. Catalyst preparation and deactivation issues for nitrobenzene hydrogenation in a microstructured falling film reactor[J].Catalysis Today，2003，81（4）：641-651.

[40]Zaid T A，Benmaza K，Chitour C E. Sulfonation of linear alkyl benzene (LAB) in a corrugated wall falling film reactor[J].Chemical Engineering Journal，2000，76（2）：99-102.

[41]Sun Y Z，Song X F，Jin M M，et al. Gas-liquid reactive crystallization of lithium carbonate by a falling film column[J].Industrial Engineering Chemistry Research，2013，52（49）：17598-17606.

[42]Medrano M，Bourouis M，Coronas A. Absorption of water vapor in the falling film of water-lithium bromide inside a vertical tube at air-cooling thermal conditions[J].International Journal of Thermal Sciences，2002，41（9）：891-898.

[43]尹銘，陳嘉賓，陳沛，等. 豎直管外降膜吸收傳熱傳質過程強化的研究[J]. 高校化學工程學報，2002，16（6）：602-608.

[44]Bo S S，Ma X H，Lan Z，et al. Numerical simulation on the falling film absorption process in a counter-flow absorber[J].Chemical Engineering Journal，2010，156（3）：607-612.

[45]Cheng P J，Liu K C，Wang C C. Nonlinear evolution of the travelling waves in roll coating flows of thin viscoelastic polymer falling films[J].Applied Mechanics and Materials，2013，479-480：45-49.

[46]Boulogne F，Pauchard L，Giorgiutti-Dauphiné F. Instability and morphology of polymer solutions coating a fiber[J].Journal of Fluid Mechanics，1957，704：232-250.

[47]Duprat C，Ruyer-Quil C，Kalliadasis S，et al. Absolute and convective instabilities of a viscous film flowing down a vertical fiber[J].Physical Review Letter，2007，98（24）：244502-244505.

[48]Stober K E，Amos J L. Method for polymerizing styrene：US，2530409[P]. 1950-11-21.

[49]橋本博，楯進，千葉明. ポリェステルの連続重縮合方法：JP，昭64-81819[P]. 1989-3-28.

[50]Vijayraghvan K，Gupta J P. Thickness of the film on a vertically rotating disk partially immersed in a Newtonian liquid[J].Industrial Engineering Chemistry Fundamentals，1982，21（4）：333-336.

[51]Cheong S I，Choi K Y. A study on the polymer layer-forming phenomena in a rotating disk polycondensation reactor[J].Journal of Applied Polymer Science，1995，55（13）：1819-1826.

[52]Cheong S I，Choi K Y. Melt polycondensation of poly(ethylene terephthalate) in a rotating disk reactor[J].Journal of Applied Polymer Science，1995，58（9）：1473-1483.

[53]王良生，戴干策. 圓盤反應器成膜性和持液量研究[J]. 化學反應工程與工藝，2000，16（2）：127-135.

[54]周賢爵，戴干策. 籠式反應器中的流動與混合特性[J]. 高校化學工程學報，2002，16（1）：7-12.

[55]Hasan N，Naser J. Determining the thickness of liquid film in laminar condition on a rotating drum surface using CFD[J].Chemical Engineering Science，2009，64（5）：919-924.

[56]Xi Z H，Zhao L，Liu Z Y. New falling film reactor for melt polycondensation process[J].Macromolecular Symposia，2007，259（1）：10-16.

[57]陳文興，馬建平，王建輝，等. 滌綸工業絲液相增黏技術的研發[J].合成纖維工業，2013，36（3）：1-4.

[58]陳文興，馬建平，王建輝，等. 滌綸工業絲熔體直紡生產技術的研發[J]. 合成纖維工業，2013，36（4）：1-4.

Advances in application of falling liquid film for reaction engineering

CHEN Shichang，MA Jianping，ZHANG Xianming，CHEN Wenxing
（National Engineering Laboratory for Textile Fiber Materials & Processing Technology (Zhejiang)，Zhejiang Sci-Tech University，Hangzhou 310018，Zhejiang，China）

Owing to high heat and mass transfer coefficient，large gas-liquid contact area and quick surface renewal，the falling liquid film flows on vertical structure can be used in chemical and polymerization reaction engineering. Experimental investigation，theoretical analysis and numerical simulation of falling liquid film flows on vertical structure are reviewed. The application of falling liquid film of low viscosity fluid in chemical engineering and falling film of high viscosity fluid used in polymerization reaction are introduced respectively. The falling liquid film flows of low viscosity fluid are more intensely investigated，while the study on high viscosity fluid only concerns film formation. The paper proposes the possible trends of falling liquid film. Falling liquid film flows of high viscosity fluid on vertical structure which involves the transfer of heat and mass coupled with chemical reaction will be the focus of future research.

falling liquid film；flow instability；heat transfer；mass transfer；polymerization；reaction engineering

TQ 021

A

1000-6613（2014）10-2528-07

10.3969/j.issn.1000-6613.2014.10.002

2014-02-26；修改稿日期：2014-04-13。

浙江省重大科技專項計劃（2012C01007）及中國紡織工業聯合會科技指導性項目（2011139）。

陳世昌（1988—），男，博士研究生。聯系人：陳文興，教授，博士生導師，主要從事功能性纖維的研究。E-mail wxchen@zstu. edu.cn。