微通道限域空間內的氣泡破裂研究進展與展望

馮俊杰，孫冰，石寧，高正明，孫萬付

（1 中國石油化工股份有限公司青島安全工程研究院，化學品安全控制國家重點實驗室，山東 青島 266100；2 北京化工大學化學工程學院，北京 100029）

微反應器具有體積小、比表面積大、反應效率高、控制精確等優點，同時能從提高移熱速度、約束危險物質分布、減少物料囤積等方面提高化工過程的本質安全性，有望成為化工產業高效、安全、綠色的革命性技術[1-2]。微反應器的研究自20 世紀90 年代開始，杜邦公司于1993 年利用微電子加工技術制造了世界上第一個芯片化工廠，用于生產甲基異氰酸甲酯等有毒物質，麻省理工學院的研究者于1997 年提出了用于氨氧化反應的T 形薄壁微反應器，此后歐洲的美因茨（IMM）、巴斯夫（BASF）、殼牌（Shell）等公司也相繼開展了微反應器研究，日本政府在2000 年制定的“國家產業技術戰略”中把微反應器技術列為新的化工技術之一。我國“十二五”規劃也將微化工技術列為中國化學工業可持續發展的關鍵技術之一，清華大學、中國科學院、南京工業大學、四川大學等研究單位持續推進了微反應器在化學、化工、材料、能源和環境等領域的應用[3-4]。一般認為微反應器尤其適用于放熱劇烈、反應物或產物不穩定、配比要求嚴格的快速反應、高含能及精細高值化學品合成等過程，在氧化、加氫、硝化、磺化、重氮化、鹵化以及材料制備和有機合成等方面已經取得了顯著的研究進展。針對微反應器的設計、加工、化學反應、傳遞規律以及工業應用的研究，國內外正方興未艾。

微反應器中應用的絕大多數為多相體系，實現微通道尺度下多相傳遞過程的深入理解與精確控制是微反應器發揮諸多優勢的前提和廣泛應用的基礎保障[5]。微通道的特征尺寸一般在10～1000μm 之間，內部的傳遞過程一般仍滿足連續介質假設，但是諸多多相流體力學行為與宏觀尺度下存在顯著的差異，很多宏觀條件下影響較小的作用力在微通道中具有不可忽視的重要影響[6]。除微反應器外，微通道尺度的多相流動還廣泛存在于傳統反應器的局部位置（如催化劑/多孔介質內部、顆粒間隙等）以及檢測分析、生物醫藥、微換熱、微混合等重要領域[7-8]。深入開展微通道多相流研究，不僅能夠豐富并加深多相流體力學機理認識，更有助于實現微通道尺度傳遞過程的精確測定和控制，為微反應器在內的多種微流控技術提供理論支撐和科學指導。

本文介紹了針對氣泡等微尺度分散相行為的研究進展，包括主要考察體系、相關研究方法和機理分析思路等，然后重點探討了含有固相等局部限域體系的氣泡破裂研究進展，最后總結了當前研究的不足，從研究體系構建、機理分析、研究方法提升以及與工藝匹配性提升等方面探討了未來的研究與發展方向。

1 微尺度分散相行為研究進展

微通道多相流體力學主要研究微通道內多種流體同時流動、共同作用時的力學問題，主要考察流動特征、速度、相含率、壓力、停留時間、分散相特征、混合和分離效果等，普遍認為多相流體力學行為由流體特性和通道特性共同決定，而潤濕性質、流動速度、界面作用、流體物性、通道結構、材質等均具有重要影響[6]。近年來具有代表性的微通道多相流研究總結如圖1所示。其中，氣泡和液滴等分散相不但能直接影響速度場，而且能通過改變相界面性質等決定體系的傳質和傳熱效果，深入理解并調控氣泡/液滴的形變與破裂行為、優化尺寸分布與多相分散過程對微通道設備的性能具有直接而重要的影響[9]。

1.1 微尺度多相流主要研究方法

目前微通道尺度多相流的主要研究方法可以分為理論分析、實驗研究和數值模擬三類。在理論分析方面，量綱為1分析、流動剖面分析、流型分布圖［圖1h］、數學建模等方法是常用的研究手段[20]。微尺度下的多相動力學機理與宏觀尺度存在顯著差異，表面力、黏性力等作用力相對于重力、浮力往往占據更重要的作用，研究者常用毛細管數、韋伯數、雷諾數等量綱為1數描述微通道內流動情況[14]。

在實驗測量方面，研究者多年來開發出多種多樣的設備對流體的流變特性、壓力（壓強）、速度、溫度、傳質特性等相關量進行實際測定。借助顯微鏡、高速相機、激光誘導熒光（LIF）、粒子成像測速［即PIV，圖1e］、激光多普勒、核磁共振技術以及X射線斷層掃描等非接觸式測試手段具有不干擾待測流場、獲取信息量較大等優勢，已成為微通道流體力學實驗測試的主流方法[21]。其中，激光誘導熒光技術利用激光將粒子激發至高能的量子態，然后探測粒子釋放出的熒光光子，可以有效揭示微通道內的濃度場和混合效率，具有高空間分辨率、高靈敏度（探測下限可達106個粒子）、對流場干擾小等優勢，但是該方法目前常檢測的物質有限，對于其他物質的熒光分離需要較高的光譜基礎和相對昂貴的成本。micro-PIV 技術借助顯微鏡系統對示蹤粒子的瞬時捕捉，可以得到微通道內全流場的瞬時速度矢量信息，目前針對單相體系及分散相含量不太高的多相體系流場測定已較為成熟，而微觀尺度的體三維流場測試也已經實現。研究者普遍認為由于微通道的尺度特點，PIV、LIF 等方法在微尺度流體力學行為研究方面具有明顯的優勢，但是幾乎所有光學方法都需要解決的問題是如何在微小尺度下穩定提供所需要的光路條件，特別是在多相體系中光經過相界面的折射和反射會產生非常大的干擾，制約了可視化方法的應用，因此最大程度消除諸多干擾因素、提高復雜多相體系的實驗精度和捕捉關鍵信息的能力仍然是微尺度可視化實驗研究的關鍵[22-23]。圖2為微通道多相流實驗研究手段。

圖1 微通道多相流研究概覽

圖2 微通道多相流實驗研究手段[23]

在數值模擬方面，對于分散相數量較少的體系模擬研究的重點是相界面的精確計算表達，通用的CFD 方法大多是通過聯立求解Navier-Stokes（N-S）方程和界面追蹤（或界面捕捉）方法來研究氣泡的動力學行為，目前廣泛采用的包括Front-tracking方法、VOF模型、level-set 等研究方法；隨著分散相個數的增加，混合模型（Mixture Model）、雙歐拉（Euler-Euler）模型、群體平衡模型（PBM）等被廣泛應用于多相流場和分散相行為的數值模擬[24]。近年來，格子玻爾茲曼方法（LBM）作為一種新興的計算流體力學方法，得到了越來越多的關注與應用，其核心思想是將流體離散為在網格上運動的介觀粒子，通過計算粒子的碰撞和遷移規律得到粒子分布函數，進而統計計算得到宏觀變量如壓力、流速和溫度等，從而實現了模擬流體運動的連續介質模型向離散模型的轉變[25]。

1.2 氣-液及液-液兩相體系的分散過程

微通道中的氣-液及液-液兩相流動在許多方面具有相似性，均可分為連續流（continuous-phase microreactors）和分散流（dispersed-phase microreactors）兩大類，具體則分為泡狀流、彈狀流、并行流、環狀流、攪拌流及其他多種過渡流型[26]。很多研究者認為兩相分隔明顯的彈狀流（泰勒流）是較為理想的流動狀態，與其他流型相比具有操作區間寬、軸向返混小、徑向混合好等優點，一層較薄的液膜是兩相之間唯一的傳遞渠道，其厚度往往隨著毛細管數的增大而增加[10]。

分散相的生成、破裂等多相流體力學行為是氣-液及液-液兩相流研究的核心。一般認為微通道中的氣泡破裂存在擠壓（squeezing）、滴落（dripping）、射流（jetting）、穿透（threading）等機制［例見圖1e］。T（Y）形、十字形、同軸形、聚焦形等結構是微通道中研究氣泡破裂的主要單元型式，研究者已針對微通道中的氣-液兩相體系建立了多種預測氣泡破裂結果的半經驗特征模型，其中毛細管數（Ca數）等量綱為1數常被用來作為判據[27]（圖3）。Roumpea 等[28]系統考察了十字形微通道中表面力在氣泡頸拉伸變薄直至斷裂過程中的影響規律。Garstecki 等[29]針對T 形微通道中氣泡破裂過程的受力狀況進行了機理分析，進而建立了破裂結果預測模型，研究者在該模型基礎上針對具體物性參數及操作條件不斷修正，在T形結構的氣泡破裂研究方面取得了一系列成果。Fu等[12]綜述了氣泡在多種微通道結構中的完全阻塞破裂、部分阻塞破裂和無阻塞破裂等破裂形式，分析了不同破裂形式的轉變機制（圖4）。整體來看，在液相剪切或沖擊等作用顯著的無阻塞自由破裂行為方面已經有了較為深入的理解，但是針對有明顯阻塞破裂等限域條件下的破裂行為的機理仍未明確。

圖3 微通道中的流型[14]

圖4 T形微通道中的典型氣泡破裂行為[12]

微通道中氣泡/液滴之間同樣不可避免地會發生相互作用，如聚并、碰撞、反彈、滑移等。與破裂過程相似，微通道中液滴或氣泡的聚并同樣受微尺度作用力的影響，聚并與否往往取決于液滴或氣泡間的接觸時間和液膜排液時間（drainage time）的關系，由于微通道尺度液膜排液和接觸時間均遠小于傳統反應器尺度，微尺度下液滴/氣泡聚并臨界條件也與宏觀尺度有較大區別[11]。荷蘭Delft大學的Kreutzer 等[30]針對微流體液滴、微通道內界面動力學及多相微反應器的分配模式開展了許多深入的研究，特別是對氣泡在T形通道口的聚并和破裂等行為提供了很多極具意義的實驗結果和模型。研究者還發現微反應器尺度中氣泡間存在特殊的連接行為（bubble train），同時發現該流型能夠提高傳質效果[31]。

1.3 微通道并行放大

在反應器的放大過程中，與傳統反應器在尺寸上逐級放大不同，微反應器的放大可以直接通過增加反應單元數量即“數量放大”（numbering-up）來實現，數量放大的方法大大簡化了放大過程的探索及傳遞過程的復雜建模步驟，有利于低成本地快速實現產業化和規模化[32]。數量放大的關鍵是要保證每個微反應單元內的流動和混合情況一致，多個微通道的并行同時也帶來了流體分配與集成化方面的問題，相對于單相流動，多相流的并行過程更加復雜和難以控制，其中便涉及了氣泡和液滴的破裂行為，一旦相分配不均勻，在實際過程中容易造成“竄流”等微觀物料配比惡化的情況。如何實現多個微通道內多相流均勻、可控地并行分配目前尚未有統一的規律性認識，已成為阻礙微反應器大規模工業應用的瓶頸之一[33]。圖5為微通道并行放大常見結構。

圖5 微通道并行放大常見結構[1]

典型的氣液兩相體系的并行放大結構主要有梳狀和耙狀等，雖然微通道內大多為層流流動，然而其通量較小，依然存在流動的不穩定性，在多通道并行過程中通道之間的耦合作用和流體分配受到流動穩定性的影響尤為明顯[34]。針對多相流動的不穩定性與波動性，研究者通過不規則分析、混沌分析和時間頻率分析等方法進行了研究，Kakac 等[35]總結了密度波理論、壓降模型及熱振蕩性三種方法在氣液兩相流不穩定性研究的應用情況，并將實驗值與模型值進行了對比。

為使流體在各個分支管路之間分配均勻，首先需要保證上游的入口條件保持穩定，而分散相和連續相進口管徑和速度往往需要分別控制，各個分支通道的管徑往往需要有差異，毛細管數需要足夠大才能保證氣泡在分支管路口的破裂。Chen等[36]研究了多種進口流型下5 路并行T 形通道中的氣液流動，發現不同的進口流型會造成不同的分配情況，而進口流型為彈狀流時相分配相對均勻。研究表明，流體的壓力與壓降對并行處的分配過程具有重要影響，壓降主要由摩擦（主要原因）、管徑收縮、流體加速度等原因造成，Kawahara等[37]發現單相流摩擦系數符合傳統層流關系式，均相流模型預測彈狀流和環狀流壓降時會偏高，而分相流模型預測結果吻合更好。Liang 等[38]通過對雙噴嘴微通道分布器結構的研究，認為進口流型及下游通道阻力對氣液相分配起到主要作用。

研究者認為此領域需要進一步利用先進的實驗測試手段，而模擬多相分配過程需借助三維而不僅僅是傳統的二維數值模擬方法，方能得到具有普適性的機理模型[36]。Tonkovich等[39]對微反應器模塊集成放大的優勢和流體分配不均勻的問題進行了討論，并針對蒸汽甲烷重整微反應器的工業應用提出了一套可伸縮設計策略。針對多相催化過程，Jensen 等通過在氣-液兩相流體接觸前設置壓降遠大于并行通道及出口集流器的增壓通道，使得流體均勻分布程度高于90%，并將之應用于臭氧化、環己烯加氫、氫氧直接合成雙氧水、直接氟化等氣-液反應過程，Mendorf 等[40]采用類似設計，使用兩種分布器實現了8 個通道內的兩相流體均勻分布，但其設計方式操作較為復雜。

2 微通道限域結構內的氣泡破裂過程

微尺度多相流研究早期以氣-液、液-液等兩相體系為主，近年來隨著傳遞機理的理解深入以及測量與控制技術的進步，針對含有固相等復雜反應體系的研究不斷增加，大大拓寬了微化工技術的研發與應用范圍[3]。利用微反應器開展固相催化反應能夠發揮其傳遞效率高、選擇性好、安全性高等優勢，例如針對許多本征速率較快、放熱量大的加氫及氧化等反應［例見圖1中i、l、m］[15,17]，氣-液-固三相微反應器有望成為替代傳統工藝的新途徑（圖6）。同時，微反應器能夠克服傳統釜式反應器返混程度大、過程參數不易控等缺點，在許多納米材料制備、快速沉淀等產物為顆粒的工藝中具有可觀的應用前景（例見圖1中n、o）[41]。

圖6 氣-液-固三相微反應器[4]

微固定床等固定顆粒形式的多相微反應器由于能夠有效減少流動顆粒引起的堵塞、分配不均勻、催化劑和產物分離困難等問題，成為當前微化工技術研究的熱點與重點。在此基礎上，含有移動顆粒的多相體系更加復雜，也對微通道多相流的認識與控制提出了更高的要求，控制不當則會出現分配不均勻、壓降過大、死區及堵塞等問題[42]。當分散顆粒均勻有序時，往往能夠有效提高傳遞和反應速率，例如以納米顆粒作為分散相的“納米流體”作為高效傳熱介質為需要快速換熱的反應過程提供了新途徑[43]。

固體顆粒強化氣液混合與傳質一直是化工過程的重要強化手段之一，微固定床等多相反應器中氣泡的形變、破裂等行為不僅能顯著影響流場分布，還能通過改變相界面及傳遞過程直接影響反應效率。氣泡在微反應器中固定顆粒床層的形變與破裂過程是由通道結構（壁面作用）、操作條件、流體/顆粒性質，特別是顆粒層內部受限空間結構等多因素共同決定的，具有鮮明的瞬時性、非均勻性、復雜性等特征。

2.1 微尺度氣泡界面的局部不穩定性機制

由于微通道中的氣泡尺寸較小、“剛性”相對較強，往往需要通過調節通道結構及流場特征制造氣泡局部的不穩定結構（如“氣泡頸”）以提高破裂發生概率。對于微填充床等體系，氣-液界面尺寸與顆粒孔隙明顯大于孔隙尺度，固定顆粒影響下的氣泡行為實質同樣是限域結構中存在較強壁面作用的氣-液兩相過程，因此若想明確含有固體顆粒、多孔介質乃至更復雜體系中的多相分散機理，首先需要探究具有較強的局部不穩定性作用下的氣泡破裂機制[44]。

氣泡的扭曲變形和破裂等行為首先需要經過擠壓排液進而形成局部薄弱位置的過程，界面演變過程中主要作用力的時空變化規律分析是研究的重點，有研究者通過平面膜/凸面膜、軟球-硬球/壁面等模型解析動力學機制，探討并解析了壁面、顆粒及流場等多因素耦合作用下的擠壓排液機理[13]。氣泡液膜排液變薄至一定厚度后發生不穩定斷裂的過程廣泛存在于破裂、聚并等重要的流體力學行為中，而不同條件下的不穩定性機制往往存在差異，很多時候甚至呈現一定的概率分布特征[35]。綜合來看，界面不穩定性機制可以分為宏觀機制與微觀機制兩種，其中宏觀尺度包括靜態和動態不穩定性。靜態不穩定性包括流量分配不穩定性、流型轉變以及Ledinegg不穩定性等，而動態不穩定性包括密度波不穩定性、熱-聲不穩定性、壓力降脈動以及其他復合型式的不穩定性等[45]。

微觀尺度的界面不穩定性則包括Rayleigh-Taylor、Kelvin-Helmholtz、Richtmyer-Meshkov、Plateau-Rayleigh 等機制，其中Rayleigh-Taylor 不穩定性描述的是在非穩定的密度分層的狀況下流體內部密度不同區域之間產生的相互滲透[46-47]。該機制會在界面附近形成流體力學不穩定區，一定程度上能夠推動流體的混合，對于微通道中分層流等流型中的界面破裂過程也具有重要作用。梁宏[48]利用格子玻爾茲曼方法研究了二維及三維微通道內的Rayleigh-Taylor不穩定性問題，分析了雷諾數對不穩定性各發展階段的擾動增長及流體界面變化的影響，并將不穩定性的發展分為線性增長、最終速度、重加速及混沌混合階段（圖7）。

圖7 兼顧熱量傳遞的并行氣液兩相界面流體力學不穩定性機制[5]

Kelvin-Helmholtz 不穩定性描述的是界面兩側切向速度差引發的不穩定現象，Kelvin-Helmholtz不穩定性會引起界面垂直于運動方向的波動，在低于一定閾值時界面張力可以維持界面穩定，而高于此閾值時波動會逐漸發展直至界面斷裂[3]。劉佩堯等[49]采用界面跟蹤法研究分析了初始振幅、波數、表面張力、密度比、速度梯度、韋伯數等對Kelvin-Helmholtz 不穩定性發展的影響規律，鄧小龍[50]通過level-set 函數輔助Cut Cell 方法液滴表面的Kelvin-Helmholtz不穩定性機制。

Richtmyer-Meshkov 不穩定性指的是流體界面受到沖擊作用時，由于密度梯度與壓力梯度的不重合導致斜壓渦量的產生，界面上的初始擾動隨著時間不斷增長，在擾動發展的后期界面上會生成許多小尺度渦[51]。Richtmyer-Meshkov 不穩定性在渦動力學以及湍流形成機理研究方面具有重要的學術意義，相對來說在微通道層流體系中影響較小，而激波與相界面的相互作用及湍流混合過程是該不穩定性研究的經典問題。

Plateau-Rayleigh 不穩定性指的是由于局部表面張力的不平衡（從微觀的分子熱運動到宏觀的流動不對稱均能導致表面張力的不平衡）對氣泡產生微小的擾動，Plateau和Rayleigh等提出可以用一個正弦波來描述這個擾動波動，使液膜某些部分半徑增大，某些部分半徑減小，進而產生了壓強分布變化。對于氣泡來說，氣泡頸的半徑變大原本會使表面張力帶來的壓強減小，但是波動又會帶來壓強的增大，比較波峰波谷兩處壓強的變化，擾動較小時波峰處壓強較大，波谷處壓強較小，可以恢復平衡，但當擾動較大時，波谷處壓強較大，恢復平衡的難度增加；隨著時間增長，擾動波會越來越強烈（perturbations growth），當達到一個臨界半徑（critical length）時就會發生斷裂。這種現象也可用能量最小化原理給予解釋，即體系傾向于縮小表面積來使能量達到最小化[52]。Plateau-Rayleigh不穩定性機制對氣泡頸的斷裂、液柱破碎成液滴等許多過程具有重要意義，尤其在十幾到幾百個微米的微小尺度條件下慣性力往往是可忽略的，范德華力也不是很顯著，而表面力作用更加突出，Plateau-Rayleigh 不穩定性的發展就更加容易，通過對微通道限域結構中的氣泡頸及液膜Plateau-Rayleigh 不穩定性機制研究有助于深入理解并調控微尺度的氣泡/液滴破裂（圖8）。

圖8 形變氣泡的Plateau-Rayleigh不穩定性機制

2.2 含固相體系等限域結構中的氣泡形變與破裂行為

相對來說，目前針對氣-液-固三相微反應器等復雜含固體系的反應動力學研究較多，而流型及多相流體力學行為的研究則較少報道。對于尺寸小于微通道結構的氣泡破裂過程，剪切應力效果明顯，部分體系甚至出現射流（jetting）型的氣泡破裂；接近或者大于通道尺寸的氣泡破裂由于通道中的液相受到阻塞，導致連續相剪切作用影響下降，而上游的壓力累積對氣泡的形變及液膜斷裂具有重要的驅動作用[12]。

Wang 等[53]利用Plateau-Rayleigh 不穩定性機制分析了T形微通道中的氣泡破裂條件，并提出了預測氣泡破裂與否的經驗公式。針對受限狹窄空間對氣泡破裂行為的影響，Dai 等[54]實驗研究了氣泡流經縮徑等限域結構的破裂行為，發現破裂結果主要由微通道尺寸、毛細管數及連續相黏度決定；Zhang 與Wang[8,55]等分別研究了液滴/氣泡流經微縮口時的破裂過程，發現不同體系均存在發生破裂的臨界毛細管數，而子氣泡數目一般隨著流量、縮口角度、液滴尺寸的增加而增加。

在微通道中的氣泡/液滴與細小顆粒的相互作用方面，常見的作用形式有碰撞、吸附、滑移、反彈、穿透等多種，研究者已針對顆粒在氣液界面的掠過效應、對邊界層的作用以及對氣含率的影響等開展了實驗研究，并分析了氣泡對微通道中較小顆粒的捕捉和聚集效應，探討了細小顆粒對氣泡聚并的阻礙作用，見圖9[56]。近年來針對氣-液-固三相微反應器中的流型轉變規律的研究也逐漸增加[57]。Liedtke等[17]對微反應器內“漿態泰勒”流三相加氫反應進行了研究，發現該種反應器在提高了傳質效率和反應器性能的同時也存在微通道易堵塞、催化劑和產物分離困難等問題。Kobayanshi 等[15]實驗發現氣-液-固三相加氫微反應器中發生環狀流動時液膜較薄，有效降低了傳質阻力，使得反應器的產物轉化率明顯提高。謝靈丹[58]研究了質量分數為0.2%的氧化銅顆粒作為分散相對于氣液兩相流型及傳質過程的影響，發現分散顆粒的存在增加了流型分布圖上Taylor流發生的區域，減少了泡狀流和環狀流出現的區域，We數可用來預測流型的變化。

圖9 氣泡生成過程的顆粒聚集作用[56]

數值模擬方面，早先的氣液固三相數值模擬研究往往將液固兩相統一描述，即簡化為氣相-漿態兩相問題，近年來隨著計算手段的發展，三相分別討論的三相歐拉及歐拉-歐拉-拉格朗日方法越來越多應用到氣液固三相模擬中，諸如CFD-VOFDPM和CFD-DED等多相耦合模型的研究逐漸開展并優化[59]。Han 等[60]利用CFD-DEM 模型研究了三相流化床內的顆粒運動和相含率分布，Tang等[61]采用CFD-VOF-DPM 耦合模型對微通道內低固含率條件下的氣-液-固三相流動特征進行了研究，發現顆粒的存在增加了液體的有效黏度，導致氣栓脫離時間和長度都縮短。

2.3 微尺度限域體系的氣泡破裂研究展望

針對含有固相的氣-液-固、液-液-固等復雜限域體系的流體力學作用機制研究是多相微反應器研發與大規模應用的必經之路，也逐漸成為微反應器研究的重點，但是整體來看當前的研究仍不夠系統、深入：針對微通道中氣泡在縮口及狹縫等受限空間的破裂研究仍處于起步階段，對限域壁面影響下的微氣泡擠壓變形及臨界破裂等顯著有別于宏觀尺度的過程認識尚不深入，而顆粒層等復雜體系中的氣泡形變-破裂演變機制研究更是非常欠缺。

鑒于研究對象復雜、研究難度相對較大的特點，建立清晰合理、逐漸深入的研究體系是微通道多相流研究的必經之路。復雜三相輸運過程的調控需要基于對于微尺度局部限域空間內的不穩定性機理的深入理解，因此遵循微固定床-微流化床、限域結構-復雜顆粒體系的研究路徑，從微通道狹縫與“孔喉”等典型限域結構入手、逐步深入拓展到固定顆粒體系的分散相行為，有助于實現系統、深入地理解氣泡形變與破裂機理。

在實驗研究方面，侵入式測量方法由于微通道的尺寸狹小及干擾流體等問題難以普遍運用，高速攝像、PIV等可視化手段在很多微通道體系中面臨的材料不透光性、界面折射與反射干擾等難題，限制了非接觸式測量方法的測試范圍。進一步可利用顆粒追蹤、優化的PIV處理算法、精確的數字圖像分析等方法對可視化技術進行改進，例如通過進行不同相之間的折射率匹配可以大幅降低多相體系對光傳播的阻礙與干擾，目前該方法在固-液、液-液等體系中均已證明可行。此外，推進納米光鑷、原子力顯微鏡等方法在微尺度動態作用研究的應用有望對傳統的光學測試手段提供有益的補充[62]。

在數值模擬方面，分散顆粒的加入使得多相系統受力機理較為復雜，許多研究者對于顆粒引起的相界面受力變化往往追求均勻化與簡化描述，同時分散顆粒往往尺寸較小、數量較多，與其他相之間經常會出現跨尺度的情況，使得模擬設置和計算的難度大幅增加。以改進的多相流耦合模型為主要數值模擬方法對實驗進行補充參照具有重要意義，其中耦合了level-set 與VOF的CLSVOF模型在許多表面張力與自由界面形態較為重要的過程中（如氣泡的Plateau-Rayleigh 不穩定性問題）具有明顯的精度優勢，該模型通過相函數初始化、相界面重構、level-set 函數重新距離化及對流輸運方程求解等步驟進行相界面追蹤，通過定義光滑化的Heaviside函數描述相含率，能夠兼具level-set方法拓撲描述精確與VOF模型相含率高守恒性的優勢[63]；利用浸入邊界、自適應網格等方法能夠大幅提高網格質量與計算效率，結合流體力學分析和數學模型化，能夠大幅提高研究精度與計算效率。

在機理分析方面，對于各種不同的限域結構，都應當抓住氣泡擠壓排液及液膜不穩定斷裂機理這兩個關鍵科學問題。氣泡的扭曲變形和破裂等行為均經過擠壓排液進而形成局部薄弱位置的過程，而微固定床等體系中顆粒影響下的氣泡行為實質同樣可理解為限域結構中存在較強壁面作用的氣-液兩相過程；應當通過針對性的實驗與數值模擬手段，系統分析界面演變過程中主要作用力的時空變化規律，結合平面膜/凸面膜、軟球-硬球/壁面等模型分析共性動力學機制，總結凝煉壁面、顆粒及流場等多因素耦合作用下的擠壓排液機理。綜合Plateau-Rayleigh 不穩定性等微觀破裂機制，建立基于流體力學穩定性的臨界擾動波長分析等理論模型，并根據不同條件下的氣泡特性及關鍵作用因素進行拓展改進與模型化（如界面作用的曲率表達方式），有助于增進對微通道乃至宏觀尺度下的分散相界面不穩定斷裂機理認識。

此外，微通道多相流研究應注重實際需求及工藝特點的結合，將工藝匹配性與工程化應用的前景作為貫穿研究始終的考慮重點（圖10）[64]。例如根據工藝對于傳質系數的需求、通過調節氣泡破裂過程控制氣-液相界面面積；在并行放大方面開發提升多相流體分布穩定性的方法，實現氣泡及能量分散有序，同時要考慮真實操作條件如工況波動、設備加工精度、物性和溫度變化等因素；在反應器設計過程要兼顧傳遞-反應效率與安全性，同時考慮壓降等對工藝流程具有重要影響的參數。

3 結語

進入21 世紀以來，微反應器等微通道設備在許多領域得到了越來越多的關注和應用，微通道中的多相流體力學研究也取得了許多顯著的成果，但是作為一個新興的領域，由于微通道的尺度特征，同時多相流動往往存在瞬時性、非均勻性、復雜性等特征，很多機理仍未明確，實現多相微通道設備的高效設計和工業應用依舊存在困難。

本文介紹了微尺度分散相行為的研究進展，深入理解并調控氣泡的形變與破裂、優化尺寸分布與多相分散過程對于體系的熱質傳遞、反應效率乃至并行放大均具有重要意義，研究者已經開展了一系列卓有成效的研究，但是對于微通道復雜限域結構內的氣泡形變與破裂機理認識仍不夠充分。重點分析了微尺度氣泡界面的局部不穩定機制，總結了含有固相等限域體系的研究進展，然后圍繞研究體系構建、研究方法改進、相間作用解析以及工程化需求匹配等方面對下一步研究與發展方向提出了展望。通過構建從擬二維狹縫、孔喉結構逐步拓展到固定顆粒體系的研究思路，一方面可以制造典型的氣泡局部不穩定結構，以便于考察破裂過程；另一方面，能夠多維度、多層次地系統研究受限狹窄空間內的氣泡形變與破裂機理，結合多種精確高效的實驗及數值模擬手段，突破已有研究中僅單維度或單層次分析的局限，有助于豐富并加深微通道尺度的多相流體力學機理認識，推進微通道內多相流的精確控制與傳遞過程強化，為多相微反應器技術提供理論支撐和科學指導。